- 1Department of Hepatopancreatobiliary Surgery, The Affiliated Hospital of Qinghai University, Xining, Qinghai, China
- 2Department of Hepatopancreatobiliary Surgery, The Chongqing University Fuling Hospital, Fuling, Chongqing, China
Objective: This study aimed to investigate the clinical characteristics and risk factors of patients with hepatocellular carcinoma (HCC) with extrahepatic metastases (EHM) and to establish an effective predictive nomogram.
Methods: Clinical and pathological data from 607 patients with hepatocellular carcinoma admitted to the Affiliated Hospital of Qinghai University between 1 January 2015 and 31 May 2018 were documented, as well as demographics, clinical pathological characteristics, and tumor-related parameters to clarify clinical risk factors for HCC EHM. These risks were selected to build an R-based clinical prediction model. The predictive accuracy and discriminating ability of the model were determined by the concordance index (C-index) and the calibration curve. The results were validated with a bootstrap resample and 151 patients from 1 June 2018 to 31 December 2019 at the same facility.
Results: In multivariate analysis, independent factors for EHM were neutrophils, prothrombin time, tumor number, and size, all of which were selected in the model. The C-index in the EHM prediction model was 0.672 and in the validation cohort was 0.694. In the training cohort and the validation cohort, the calibration curve for the probability of EHM showed good agreement between the nomogram prediction and the actual observation.
Conclusion: The extrahepatic metastasis prediction model of hepatocellular carcinoma constructed in this study has some evaluation capability.
Introduction
Hepatocellular carcinoma (HCC) is the leading malignant tumor from the liver, which is the seventh most common, and has the second highest death rate (in all 36 tumors) (1). It has already been reported in the article that the risk factor for poor prognosis is the lack of diagnosis of extrahepatic metastasis (EHM) (2). Therefore, an accurate evaluation of HCC metastases is critical for improved prognosis.
HCC is a kind of cancer that develops as a result of a secondary liver illness [such as viral hepatitis (HBV or HCV), alcoholic or fatty liver disease]. Liver function indices are intimately linked to the occurrence and development of HCC (3, 4). Besides, studies have reported that primary tumor progression characteristics (such as vascular invasion, tumor size and number, etc.) are independent risk factors for EHM (5, 6). The above parameters can be risk factors for metastatic HCC and will be included in this study as observational information. Since patients with HCC usually receive antitumor treatment during the clinical course, it may give a confounding effect on the analysis of metastatic factors. As a result, we evaluated the clinical features and risk factors of patients with HCC and EHM who did not receive anti-tumor therapy, and we created an effective EHM diagnostic nomogram.
Patients and methods
Patients and study design
A retrospective study was conducted on patients who were diagnosed with HCC from 1 January 2015 to 31 December 2019 at the Affiliated Hospital of Qinghai University (Xining, China). Inclusion criteria were as follows: 1) according to the Guidelines for diagnosis and treatment of primary liver cancer, the patient was diagnosed with HCC (2021 Edition) (7) and 2) with a complete medical record. Exclusion criteria were as follows: 1) no prior history of anticancer treatment; 2) no priors for other cancers; and 3) without other confirmed or suspected cancers. The training cohort consisted of patients between 1 January 2015, and 31 May 2018, and the validation cohort consisted of patients between 1 June 2018 and 31 December 2019. Depending on whether EHM was present at the time of the first diagnosis, the training cohort was further split into extrahepatic metastatic (observation group) and non-extrahepatic metastatic (control group) groups. Age and sex-related demographic data as well as clinicopathological characteristics such as body mass index, smoking and drinking history, blood tests, assessments of HBV and HCV infections, results of liver function tests, and tumor-related parameters were prospectively gathered.
Diagnosis and definitions
The appearance of a newly detected tumor confirmed on two radiologic images, with or without an elevation of serum tumor markers, was defined as metastasis. A patient with a smoking history was defined as having smoked continuously or cumulatively for 6 months or more in the past. Drinking more than three standard glasses of alcohol per day or more than seven standard glasses of alcohol per week for 1 month or more, either continuously or cumulatively, is considered a drinking history.
Follow-up
During the 2 years following diagnosis, all patients were seen once every 3 months. An abdominal ultrasound, blood test, and liver function test were performed at each of the follow-up visits. When a tumor recurrence or metastasis was suspected, a contrast-enhanced CT or MRI was performed, and the results were reviewed individually by two experienced doctors.
Statistical analyses
Statistical analyses to identify risk factors were performed using IBM SPSS Statistics 23.0 for Windows (SPSS Inc., Chicago, IL, USA). Continuous variables were compared using the Mann–Whitney U test for variables with an abnormal distribution. Categorical variables were compared using the chi-squared or Fisher exact test. In the univariate analyses, p < 0.05 was considered statistically significant. Multivariate logistic regression analysis was used to evaluate the independent risk factors of extrahepatic metastases. In the multivariate analyses, p < 0.05 was considered statistically significant.
A nomogram was formulated based on the results of multivariate logistic regression analysis and by using the ‘rms’ package in R version 4.12 (http://www.r-project.org/). A final model selection was performed by a forward conditional selection process. The predictive performance of the nomogram was measured by concordance index (C-index) and the calibration curve. Bootstraps with 1000 resample were used for these activities (Figure 1).
Figure 1 After the screening of medical records for the study, 607 cases were identified for final enrolment based on exclusion criteria. The training cohort was used to identify factors that were able to predict EHM, thereby establishing a nomogram for this study. This nomogram was then validated regarding its accuracy in the evaluation of EHM risk by using both the training cohort and validation cohort.
Results
Presentation of patients
In this study, during the defined study period (January 1st,2015 to December 31st,2019), 1673 cases were identified as HCC in the Affiliated Hospital of Qinghai University. 1066 cases were excluded according to the exclusion criteria, and 607 cases were finally enrolled, including 456 patients in the training cohort from 1 January 2015 to 31 May 2018, and 151 patients in the validation cohort from 1 June 2018 to 31 December 2019 (Figure 1). Among all the 170 EMH patients enrolled at the time of diagnosis, 122 patients (71.8%) had lymph node metastasis, 114 patients (67.1%) had lung metastasis, 12 patients (7.1%) had bone metastasis, 7 patients (4.1%) had diaphragm metastasis, 4 patients (2.9%) had peritoneal or pelvic metastasis, 3 patients (2.2%) had splenic metastasis, and pancreatic metastases in 2 cases (1.5%), adrenal metastases in 2 cases (1.5%) and gastric metastases in 1 case (0.8%).
Baseline of characteristics and multivariate analysis
The training cohort was concentrated on 136 cases in the EHM group, of which 112 (81.4%) were male and 24 (17.6%) were female, aged 53.42 (48.15, 62.69) years; 320 cases in the control group, of which 255 (79.7%) were male and 65 (20.3%) were female, aged 53.83 (46.82, 63.32) years. Table 1 shows further characteristics of training and validation patients. According to univariate analysis and multivariate analyses, neutrophils, prothrombin time, tumor count, and size have been shown to be independent risk factors for EHM in initial patients (Table 2 and Figure 2).
Figure 2 Multivariate analysis of the clinical characteristics for predicting EHM in the training cohort. NE, neutrophil; PT, prothrombin time.
Construction and validation of the initial patients EHM nomogram
Based on the results of the multivariate logistic regression analysis in the training cohort, neutrophil, prothrombin time, tumor number, and tumor size were used as variables to construct the nomogram (Figure 3). When the ROC curve was plotted using the training cohort, the area under the ROC curve was calculated as 0.672 (Figure 4A). In the validation cohort, the nomogram displayed a C index of 0.694 (Figure 5A). This result indicates that there is some discrimination in the mode and can be relied on to accurately predict the extrahepatic metastases of hepatocellular carcinoma. Model consistency was assessed by drawing calibration curves using data from the Training and Validation cohorts (Figures 4B, 5B). The diagonal line represents the precise match between the expected and real circumstances, the dashed line represents the model’s theoretical forecast, and the solid line represents the actual prediction obtained by repeated sampling. The two curves are less discontinuous from the diagonal line, suggesting that the model-anticipated results are more consistent with what actually occurred and the model-predicted results are more credible.
Figure 3 The nomogram predicting initial patient EHM probability. NE, neutrophil; PT, prothrombin time.
Figure 4 (A): Receiver operating characteristic curves for EHM of the patients in the training cohort. The receiver operating characteristic curves at the initial diagnosis are shown, and its area ROC curves are provided. (B): Calibration plots of EHM in the training cohorts. The calibration curves derived from the training cohorts are almost a diagonal line that would represent perfectly reliable prediction.
Figure 5 (A): Receiver operating characteristic curves for EHM of the patients in the validation cohort. The receiver operating characteristic curves at the initial diagnosis are shown, and its area ROC curves are provided. (B): Calibration plots of EHM in the validation cohort. The calibration curves derived from the validation cohort are almost a diagonal line that would represent perfectly reliable prediction.
Follow-up
The 415 patients were fully monitored for 23.694 person-months (median, 15.118 months; range, 3.7 to 62.6 months), in which 308 (74.2%) died. In all 415 patients enrolled, the overall survival (OS) of patients with EHM was significantly worse than non-EHM patients (Figure 6, 6.7 months vs. 23.1 months, p = 0.00).
A total of 372 patients were enrolled in this study without metastasis at the time of initial diagnosis and all were followed up. A total of 372 patients were enrolled in this study without metastasis at the time of initial diagnosis, all of whom underwent follow-up. Based on the EHM occurring during the follow-up period, the patients were divided into observational metastases and non-methane groups. The characteristics of the patients are listed in Table 3. According to univariate analysis and multivariate analyses, neutrophils, Total Protein, Carcinoembryonic Antigen and tumor sizes have been shown to be independent risk factors for EHM in followed up patients (Table 4).
Discussion
The incidence of HCC has increased in many countries in recent years. The primary risk factors for HCC worldwide include chronic hepatitis B virus (HBV) or hepatitis C virus (HCV) and the consumption of aflatoxin-contaminated food. The prevalence of HCC caused by metabolic syndrome, obesity, diabetes, excessive alcohol consumption, and non-alcoholic fatty liver disease (NAFLD) is gradually increasing (8). Thus, as the etiology of the disease has changed, the risk factors for EHM of HCC may also be changed, and further research is required.
For patients with HCC with EHM, most previous studies examined only the relationship between clinicopathological characteristics and the prognosis of patients with HCC (5, 9–11). However, the common shortcoming was that the patients enrolled previously received anti-tumor treatment, and the lab results were incidentally altered (to exhibit low white blood cell counts, low platelet counts, poor liver or kidney function, and so on) (9), which may bring inevitable interference to the analysis of metastases. Therefore, this study investigated the clinical characteristics and risk factors of patients with HCC with EHM who did not receive anti-tumor treatment, and established an effective diagnostic nomogram for EHM. In our study, 28% (170) of all 607 HCC patients included had extrahepatic metastases, consistent with previous study results (9, 12). In the initially diagnosed patients, the metastatic sites included the lymph node (122, 71.8%), lung (104, 61.2%), bone (12, 7.1%), diaphragm (7, 4.1%), peritoneum (4, 2.3%), spleen (3, 1.8%), pancreas (2, 1.2%), adrenal gland (2, 1.2%), and stomach (1, 0.6%). The lung may be the most common site of EHM from HCC speciously, but in our study, the proportion of lymph node metastases was the highest (122, 71.8%). Additionally, a similar finding was obtained by another Chinese study (13). In that research study, among the 132 patients with extrahepatic metastases from hepatocellular carcinoma diagnosed by whole-body PET/CT, 72 (54.5%) had metastases in the lymph node, 32 (24.2%) had metastases in the bones, and 28 (21.2%) had metastases in the lungs. This may due to the main symptom of patients with simple lymph node metastases rarely present with clinical symptoms, and only a small proportion of patients are noted when the enlarged lymph nodes caused a compression effect, for example, the jaundice caused by bile duct compression. As a result, the rate of missed lymph node metastases was high and easily ignored in a relatively earlier stage of metastasis.
According to the univariate and multivariate analysis of EHM in our research (Table 2 and Figure 2), there may be a potential relation between HCC patients with EHM and the tumor number. This relationship was also mentioned in another report which showed that the number of tumors >2 can be easier found in patients with EHM (14). They therefore concluded that the number of tumors could be associated with aggressive biological features. Coincidentally, in our research, the tumor counts are independent predictors for EHM (Figure 2). Similar findings for HCC metastasis were noted previously, which revealed that a multiple tumor number was a risk factor for EHM in HCC (15, 16).
A close relationship between HCC size and EHM seems to exist according to our results (Figure 2). One possible explanation is that the lager lesion contains more tumor stem cells, which are frequently identified as the source of malignant phenotypes such as aggressive growth, portal vein thrombosis, or metastasis. Another explanation is that the tumor biology changes beyond a certain mass, just as tumors in general cannot grow beyond a critical small size (9, 17, 18).
In this study, we also demonstrated that an elevated neutrophil was a significant independent predictor for EHM of HCC. Neutrophil is one of the most simple and effective markers of inflammation and is associated with poor prognosis in various cancers (11, 19, 20). Therefore, the immanent reason may be that the high neutrophil count is a symbol of an adequate environment for tumor progression, which has been shown to promote tumor growth and metastasis by secreting chemokines, vascular endothelial growth factor, and matrix metalloproteinase-9, which are involved in the development of local inflammation and angiogenesis (21–23).
Currently, angiopoietin-2 (Ang-2), microRNAs (miRNAs), and lncRNA BACE1-AS are available as a test for the evaluation of EHM in HCC (24–27). However, the above parameters are mainly laboratory-based and are not available in most hospitals. Therefore, our results identified the risk factors for EHM of HCC which are based on noninvasively clinically readily available data and developed a nomogram with some predictive ability. We constructed a predictive model of EHM of HCC by taking the above independent risk factors as variables. It is used to predict the initial probability of patient EHM, and validated internally and externally, confirming its certain predictive capacity for EHM. The calibration curves were drawn and showed that the nomogram predictions overlapped well with the actual clinical situation, with good agreement and credible prediction results.
In previous studies (9, 12, 14, 28), most researchers focused on the observation of the patient’s prognosis. Therefore, only a preliminary analysis of the reference indicators associated with EHM was performed. This study further explored the independent risk factors for EHM and developed a nomogram with some predictive ability, building on the previous work. It has a role to play in reducing missed EHM and designing optimal therapies for those patients. This study is a single-center study only and the model could be further improved with a larger sample size and multi-center data.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The study was censored on October 19,2020 and was approved by the Ethics Committee of The Affiliated Hospital of Qinghai University. All subjects signed an informed consent form. The patients/participants provided their written informed consent to participate in this study.
Author contributions
LZ is responsible for writing, LR, MQ, WW and XS are responsible for document of patients data. FY, MD and HW are responsible for document. following up, ZW, HF are responsible for guiding research and revising papers. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (No. 82160466) and Research team for minimally invasive diagnosis and treatment of biliary and pancreatic diseases (The Affiliated Hospital of Qinghai University).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
2. Le Linn Y, Guo Y, Koh YX, Chow PKH, Chan CY, Chung AYF, et al. Preoperative predictors of futile resection of intraabdominal extrahepatic metastases from hepatocellular carcinoma. World J Surg (2021) 45(4):1144–51. doi: 10.1007/s00268-020-05907-2
3. Peng W, Shen J, Dai J, Leng S, Xie F, Zhang Y, et al. Preoperative aspartate aminotransferase to albumin ratio correlates with tumor characteristics and predicts outcome of hepatocellular carcinoma patients after curative hepatectomy: a multicenter study. BMC Surg (2022) 22(1):307. doi: 10.1186/s12893-022-01751-4
4. Kwo PY, Cohen SM, Lim JK. ACG clinical guideline: Evaluation of abnormal liver chemistries. Am J Gastroenterol (2017) 112(1):18–35. doi: 10.1038/ajg.2016.517
5. Jun L, Zhenlin Y, Renyan G, Yizhou W, Xuying W, Feng X, et al. Independent factors and predictive score for EHM of hepatocellular carcinoma following curative hepatectomy. Oncologist (2012) 17(7):963–9. doi: 10.1634/theoncologist.2011-0447
6. Ochiai T, Ikoma H, Okamoto K, Kokuba Y, Sonoyama T, Otsuji E. Clinicopathologic features and risk factors for extrahepatic recurrences of hepatocellular carcinoma after curative resection. World J Surg (2012) 36(1):136–43. doi: 10.1007/s00268-011-1317-y
7. Department of Medical Administration, National Health and Health Commission of the People's Republic of China. Guidelines for diagnosis and treatment of primary liver cancer in China (2019 edition). Chin J Hepatology (2020) 28(2):112–28. doi: 10.3760/cma.j.issn.1007-3418.2020.02.004
8. McGlynn KA, Petrick JL, El-Serag HB. Epidemiology of hepatocellular carcinoma. Hepatol (2021) 73 Suppl 1:4–13. doi: 10.1002/hep.31288
9. Hsu CY, Liu PH, Ho SY, Huang YH, Lee YH, Lee RC, et al. Metastasis in patients with hepatocellular carcinoma: Prevalence, determinants, prognostic impact and ability to improve the Barcelona clinic liver cancer system. Liver Int (2018) 38(10):1803–11. doi: 10.1111/liv.13748
10. Carr BI, Guerra V. Hepatocellular carcinoma EHM in relation to tumor size and alkaline phosphatase Levels[J]. Oncology (2016) 90(3):136–42. doi: 10.1159/000443480
11. Ding PR, An X, Zhang RX, Fang YJ, Li LR, Chen G, et al. Elevated preoperative neutrophil to lymphocyte ratio predicts risk of recurrence following curative resection for stage IIA colon cancer. Int J Colorectal Dis (2010) 25(12):1427–33. doi: 10.1007/s00384-010-1052-0
12. Katyal S, Oliver JH 3rd, Peterson MS, Ferris JV, Carr BS, Baron RL. Extrahepatic metastases of hepatocellular carcinoma. Radiology (2000) 216(3):698–703. doi: 10.1148/radiology.216.3.r00se24698
13. Li G, Wu P, Huan H, Ma K, Li X, Bie P, et al. Extrahepatic metastasis markedly impact on the survival in hepatocellular carcinoma patients. Chin J Anat Clin Sci (2014) 19(04):307–9. doi: 10.3760/cma.j.issn.2095-7041.2014.04.011
14. Kanda M, Tateishi R, Yoshida H, Sato T, Masuzaki R, Ohki T, et al. EHM of hepatocellular carcinoma: incidence and risk factors. Liver Int (2008) 28(9):1256–63. doi: 10.1111/j.1478-3231.2008.01864.x
15. Llovet JM, Brú C, Bruix J. Prognosis of hepatocellular carcinoma: the BCLC staging classification. Semin Liver Dis (1999) 19(3):329–38. doi: 10.1055/s-2007-1007122
16. The Cancer of the Liver Italian Program (CLIP) Investigators. A new prognostic system for hepatocellular carcinoma: a retrospective study of 435 patients: the cancer of the liver Italian program (CLIP) investigators. Hepatology (1998) 28(3):751–5. doi: 10.1002/hep.510280322
17. Myszczyszyn A, Czarnecka AM, Matak D, Szymanski L, Lian F, Kornakiewicz A, et al. The role of hypoxia and cancer stem cells in renal cell carcinoma pathogenesis. Stem Cell Rev Rep (2015) 11(6):919–43. doi: 10.1007/s12015-015-9611-y
18. Tian Q, Xue Y, Zheng W, Sun R, Ji W, Wang X, et al. Overexpression of hypoxia-inducible factor 1α induces migration and invasion through notch signaling. Int J Oncol (2015) 47(2):728–38. doi: 10.3892/ijo.2015.3056
19. Jung MR, Park YK, Jeong O, Seon JW, Ryu SY, Kim DY, et al. Elevated preoperative neutrophil to lymphocyte ratio predicts poor survival following resection in late stage gastric cancer. J Surg Oncol (2011) 104(5):504–10. doi: 10.1002/jso.21986
20. Stotz M, Gerger A, Eisner F, Szkandera J, Loibner H, Ress AL, et al. Increased neutrophil-lymphocyte ratio is a poor prognostic factor in patients with primary operable and inoperable pancreatic cancer. Br J Cancer (2013) 109(2):416–21. doi: 10.1038/bjc.2013.332
21. Strieter RM, Burdick MD, Mestas J, Gomperts B, Keane MP, Belperio JA. Cancer CXC chemokine networks and tumor angiogenesis. Eur J Cancer (2006) 42(6):768–78. doi: 10.1016/j.ejca.2006.01.006
22. Gong Y, Koh DR. Neutrophils promote inflammatory angiogenesis via release of preformed VEGF in an in vivo corneal model. Cell Tissue Res (2010) 339(2):437–48. doi: 10.1007/s00441-009-0908-5
23. Aino H, Sumie S, Niizeki T, Kuromatsu R, Tajiri N, Nakano M, et al. Clinical characteristics and prognostic factors for advanced hepatocellular carcinoma with EHM. Mol Clin Oncol (2014) 2(3):393–8. doi: 10.3892/mco.2014.259
24. Liu SP, Hu YH, Zhang H, Xia GB, Chen Q. Association of quantitative peripheral blood angiopoietin-2 mRNA assay with invasive metastasis and prognosis of hepatocellular carcinoma. J Clin Dig Dis (2019) 31(05):324–8. doi: CNKI:SUN:LCXH.0.2019-05-018
25. Yang YL, Tsai MC, Chang YH, Wang CC, Chu PY, Lin HY, et al. MIR29A impedes metastatic behaviors in hepatocellular carcinoma via targeting LOX, LOXL2, and VEGFA. Int J Mol Sci (2021) 22(11):6001. doi: 10.3390/ijms22116001
26. Liu C, Wang H, Tang L, Huang H, Xu M, Lin Y, et al. LncRNA BACE1-AS enhances the invasive and metastatic capacity of hepatocellular carcinoma cells through mediating miR-377-3p/CELF1 axis. Life Sci (2021) 275:119288. doi: 10.1016/j.lfs.2021.119288
27. Lv E, Sheng J, Yu C, Rao D, Huang W. LncRNA influence sequential steps of hepatocellular carcinoma metastasis. Biomed Pharmacother 136:111224. doi: 10.1016/j.biopha.2021.111224
Keywords: primary hepatic carcinoma, extrahepatic metastases, risk factors, clinical features, the prediction model
Citation: Zhou L, Ren L, Yu W, Qi M, Yuan J, Wang W, Su X, Yin F, Deng M, Wang H, Long H, Zeng J, Yu J, Fan H and Wang Z (2022) Construction and validation of a prediction model of extrahepatic metastasis for hepatocellular carcinoma based on common clinically available data. Front. Oncol. 12:961194. doi: 10.3389/fonc.2022.961194
Received: 04 June 2022; Accepted: 25 October 2022;
Published: 16 November 2022.
Edited by:
Boris Gala-Lopez, Dalhousie University, CanadaReviewed by:
Debanjali Dasgupta, Mayo Clinic, United StatesYan Zhou, Tianjin Third Central Hospital, China
Takashi Karako, National Center For Global Health and Medicine, Japan
Copyright © 2022 Zhou, Ren, Yu, Qi, Yuan, Wang, Su, Yin, Deng, Wang, Long, Zeng, Yu, Fan and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhixin Wang, emhpeGluX3dhbmcwMDFAc2luYS5jb20=; Haining Fan, ZmFuaGFpbmluZ0BtZWRtYWlsLmNvbS5jbg==
†These authors have contributed equally to this work