- Research article
- Open Access
- Open Peer Review
Competing risk analyses of overall survival and cancer-specific survival in patients with combined hepatocellular cholangiocarcinoma after surgery
BMC Cancer volume 19, Article number: 178 (2019)
Our objective was to identify risk factors affecting overall survival (OS) and cancer-specific survival (CSS) and build nomograms to predict survival based on a large population-based cohort.
Two hundred and thirty patients diagnosed with CHCC between 2004 and 2015 were retrospectively extracted from the Surveillance, Epidemiology, and End Results (SEER) database as a training cohort. In addition, Ninety-nine patients diagnosed with CHCC between 2000 and 2017 were retrospectively extracted from Sun Yat-Sen University Cancer Center (SYSUCC) as an external validation. Nomograms for predicting probability of OS and CSS were established. Performance of the nomograms was measured by concordance index (C-index) and the area under receiver operating characteristic (ROC) curve (AUC).
In training cohort, the 1-, 2 and 3-year OS were 67.7, 46.8 and 37.9%, and the 1-, 2 and 3-year CSS were 73.1, 52.0 and 43.0%, respectively. The established nomograms were well calibrated in both training and validation cohort, with concordance indexes (C-index) of 0.652 and 0.659, respectively for OS prediction; 0.706 and 0.763, respectively for CSS prediction. Nomograms also displayed better discriminatory compared with 8th edition tumor-node-metastasis (TNM) stage system for predicting OS and CSS.
We constructed nomograms to predict OS and CSS based on a relatively large cohort. The established nomograms were well validated and could serve to improve predictions of survival risks and guide management of patients with CHCC after surgery.
Combined hepatocellular cholangiocarcinoma (CHCC) is a rare primary liver cancer, which is composed of mixed elements of both hepatocellular carcinoma (HCC) and cholangiocarcinoma (ICC)  and accounts for only 0.4–14.5% of the primary liver cancer . Regarding the treatment of CHCC, patients can obtain the best chance to the greatest survival benefit from surgery . However, CHCC has worse prognosis compared with HCC or ICC . CHCC was firstly described in 1949 by Allen and Lisa . However, due to the low morbidity of CHCC and the absence of unified diagnostic criterion, the clinical and pathological features of CHCC remain unclear. Moreover, different from HCC for which many preoperative prognostic prediction systems have been established [6,7,8], the prognostic stage system of patients with CHCC remained unclear, varying considerably from different reports [9,10,11]. The 8th tumor-node-metastasis (TNM) stage system, although it was the most frequently used stage system, it contained some common prognostic factors, such as tumor size, lymph node (LN) metastasis and distant metastasis. In addition, there was no a TNM stage system which is specially designed for CHCC. It was reported that the differences of clinical features among CHCC, HCC and ICC could lead to the variations of prognostic factors [12, 13]. There were also many factors, such as age, gender and tumor grade, which were shown to have great impact on survival. However, they were not included in the TNM stage system. Therefore, it is necessary to establish a stage system which is specially designed for prognostic prediction in patients with CHCC.
In addition, with the improvement of survival of cancer patients, most of patients with CHCC are faced with advanced ages, which are associated with an increasing high rate of comorbidities. Moreover, ta high risk of competing events, which might contribute to more competing deaths, was observed in patients with CHCC as the age increases [14, 15]. Thus, when prognosis is evaluated, competing risks are worthy of being considered. However, most prognostic analyses only focused on overall survival (OS) and ignored the impact of survival from competing events [8, 9, 16]. Competing risk analysis evaluates the informative nature of censoring and the occurrence rates of a particular event, which is more suitable for prognostic analysis. Misleading conclusions might be obtained due to the failure to recognize the presence of competing risks in survival analysis .
The present study was to build nomograms to predict 1-, 2-, and 3-year OS and cancer-specific survival (CSS) of these patients based on the Surveillance, Epidemiology, and End Results (SEER) database. Also, another large cohort of patients with CHCC from China was used to externally validate the established nomograms.
The study population was identified from SEER database from 2004 to 2015. We focused on cases pathologically confirmed CHCC after surgery [International Classification of Diseases for Oncology, Third Edition (ICD-O-3) site code C22.0 and C22.1; histology code: 8180/3]. In addition, consecutive patients with pathological diagnosis of CHCC after surgery between 2000 and 2017 at the department of Hepatobiliary and Pancreatic Surgery of Sun Yat-Sen University Cancer Center (SYSUCC) were also enrolled in the present study. The exclusion criteria are the same as those described in our previous study .
Records for the age at diagnosis, gender, tumor site, tumor grade, tumor size, TNM stage, follow-up information and cause of death were retrospectively retrieved from SEER database and the medical management system of SYSUCC. Survival time was defined as the duration from the date of diagnosis to last follow-up or death due to all causes (OS) or CHCC (CSS).
Nomogram construction and validation
Nomograms were constructed based on cohort from SEER database and externally validated based on cohort from SYSUCC database. Student’s t test and chi-square test or Fisher’s exact test were used to compare continuous variables and categorical variables, respectively. The Kaplan-Meier curves were analyzed by log-rank tests. Univariate analysis and multivariate analysis were constructed using the Cox regression model and hazard ratio (HR) and the associated 95% confidence interval (CI) for each variable were determined. Clinical and pathological factors were analyzed by the Fine and Grey’s model for their cumulative incidence function (CIF) on cancer-specific mortality and non-cancer-specific mortality. Independent prognostic factors identified in the multivariate analysis were used to build nomograms to predict the 1-, 2- and 3-year OS and CSS rates.
As two important aspects of the performance of the established nomograms, the discrimination and calibration power were evaluated by concordance index (C-index) and calibration curves, respectively . Bootstraps with 1000 resamples were used in the validation of the nomogram. In addition, the area under receiver operating characteristic (ROC) curve (AUC) was used to evaluate the precision of the survival predictions.
R version 3.4.2 software (The R Foundation for Statistical Computing, Vienna, Austria. http://www.r-project.org), along with SPSS version 22 (SPSS Inc., Chicago, IL, USA), was used to conduct statistical analyses. A two tailed P-value < 0.05 was considered statistically significant.
Two hundred and thirty patients with CHCC and another ninety-nine patients with CHCC were retrospectively identified from SEER database as training cohort and SYSUCC database as external validation cohort, respectively in the present study. The baseline characteristics of the training cohort and the validation cohort were shown in Table 1. Among these patients, the mean age was 59.8 years and 49.7 years for the patients in the training cohort and validation cohort, respectively. Most patients were male and had tumor origin from liver in both train and validation cohort. Poor differentiation (130, 56.7%) was the most common tumor grade, while most patients had tumors which were moderately differentiated in the validation cohort. The proportions of patients were comparable between two cohorts in terms of T stage (8th), N stage (8th) and 8th edition TNM stage system.
OS and CSS of patients
During the follow-up period, deaths were observed in 142 out of 230 (61.7%) patients in the training cohort and 43 out of 99 (43.4%) patients in the validation cohort. In the training cohort, CHCC contributed to deaths of 118 (51.3%) patients and competing risk events contributed to deaths of 24 (10.4%) patients. In the validation cohort, there were 31 cancer-specific death and 12 non-cancer-specific death during the follow-up period. Table 2 outlined the comparisons of 1-, 2- and 3-year OS rates, cancer-specific mortalities and non-cancer-specific mortalities of patients. It was shown that older age, larger tumor and advanced T stage (8th) were responsible for higher cumulative rates of cancer-specific-mortality. Earlier N stage (8th), well differentiation and origin from intrahepatic bile duct seemed to be related to the decreased cancer-special-mortalities while the differences were not significant (Fig. 1).
The median OS and CSS for patients were 22.0 (95% CI: 18.0–29.0) months and 27.0 (95%CI: 20.0–37.0) months, respectively. The 1-, 2 and 3-year OS were 67.7, 46.8 and 37.9%, and the 1-, 2 and 3-year CSS were 73.1, 52.0 and 43.0%, respectively. The Kaplan-Meier curves of OS analyses were shown in Fig. 2. Patients who were younger than 60 years old or had smaller tumor (≤ 5 cm) had significant longer OS. Tumor originated from intrahepatic bile ducts, well differentiated tumor, or earlier T stage (8th) also indicated better OS.
Construction and validation of nomograms
Univariate analyses were performed to filter prognostic factors. It was revealed that age, tumor site, tumor grade, tumor size and T stage (8th) were significantly associated with OS. After a stepwise removal of variables, age (HR = 1.031, 95% CI = 1.006–1.056, P = 0.015), tumor grade (HR = 2.049, 95% CI = 1.277–3.288, P = 0.003) and T stage (8th) (HR = 1.410, 95% CI = 1.071–1.855, P = 0.014) remained significant predictors for OS (Table 3). Proportional subdistribution hazard assumption for CSS analysis also showed that Age (HR = 1.038, 95% CI = 1.010–1.067, P = 0.008), tumor grade (HR = 2.027, 95% CI = 1.195–3.439, P = 0.009), tumor size (HR = 1.849, 95% CI = 1.001–3.427, P = 0.049) and T stage (8th) (HR = 1.429, 95% CI = 1.038–1.969, P = 0.029) were all independently associated with CSS.
Nomograms for predicting OS and CSS were constructed with all of the independent predictors of patients in the training cohort (Fig. 3). The C-indexes for OS and CSS prediction were 0.652 (95% CI = 0.579–0.725) and 0.706 (95% CI = 0.630–0.782), respectively, showing good accuracy of the established nomograms for survival prediction. In addition, the comparison of C-indexes of the established nomograms and the 8th edition TNM stage system showed that the established nomograms had enhanced discriminatory ability in predicting OS and CSS (OS, C-index = 0.652, 95%CI = 0.579–0.725 vs C-index = 0.567, 95%CI = 0.492–0.642, P = 0.015; CSS, C-index = 0.706, 95%CI = 0.630–0.782 vs C-index = 0.553, 95%CI = 0.469–0.637, P < 0.001). The accuracy of nomogram was verified by bootstrapped resamples via the validation cohort. Fair agreement between the nomogram-predicted survival and the actual survival was observed (Fig. 4) and it was indicated that discrimination of nomogram with regard to the SYSUCC validation cohort was also higher than that of 8th edition TNM stage system even though it did not exhibit independent significance (Table 4).
Furthermore, two ROC models of OS and CSS regarding the prediction ability were compared (Table 5). In the training cohort, the values of AUC of the nomogram for predicting 1-, 2 and 3-year OS and CSS were 0.703, 0.675 and 0.753; 0.752, 0.702 and 0.791, respectively, which were all higher than those of 8th edition TNM stage system (Fig. 5). Regarding to the validation cohort, the values of AUC of the nomogram for predicting 1-, 2 and 3-year OS and CSS were 0.638, 0.647 and 0.600; 0.775, 0.800 and 0.785, respectively, whereas the AUC values of the 8th edition TNM stage system for predicting 1-, 2 and 3-year OS and CSS were 0.630, 0.638 and 0.575; 0.722, 0.720 and 0.689, respectively (Fig. 6). The established nomograms showed superior discriminatory capacity than 8th TNM stage system for predicting OS and CSS in both training and validation cohort.
CHCC is a primary malignant tumor and represents a small proportion of all liver cancers. Due to the rarity of CHCC, most previous studies of CHCC were only limited to single-center cohorts with small sample sizes. The clinicopathological predictors of CHCC remained unclear and the special predictive system was unavailable for the personal treatment. Moreover, most previous studies mainly focused on OS, other than CSS, which reflected the nature of causes of deaths in cancer patients, especially those with increasing ages . Thus, we tried to evaluate the mortality of patients and built nomograms to predict OS and CSS for patients with CHCC after surgery in this study.
It was observed that the increasing ages had a negative effect of survival in patients with CHCC after surgery, which was more obvious on CSS than OS. Moreover, similar with other studies [20, 21], it was indicated that the increasing ages were shown to be independent prognostic factors of survival in this study. Thus, maybe considering age was more appropriate when prognosis of patients with CHCC after surgery was evaluated.
In the presence of competing risk model, other independent prognostic factors included tumor grade, tumor size and T stage (8th). Tumor size is the predominant feature of T stage (8th) and an important component of the 8th edition TNM stage system. It was shown that advanced T stage (8th) represented greater risks of lower OS and CSS in this study. In addition, heavier weight from T stage (8th) in predicting CSS than OS was observed, showing cancer-specific mortalities were more largely depended on inherent feature of tumor. Another factor reflected the intrinsic nature of tumor, tumor grade, was also associated with changes of prognoses of patients with CHCC, which was in accordance with many previous studies [12, 22, 23]. The addition of tumor grade, which was independent of other prognostic factors, such as tumor size and LN metastasis, might contribute to more accurate estimation of tumor behavior and survival outcomes of patients .
The differences of origin and the complex nature may lead to the unique features of CHCC compared with HCC and ICC. The predictive significance was not observed for LN metastasis in patients with CHCC in this study. This result was similar with that from a large-scale study . The proportion of patients who were accompanied with LN metastasis was extremely low. In this study, LN metastasis was depended on surgical resection and pathologic confirmation, other than imaging scan. This criterion could contribute to the lower rates of LN metastasis. In addition, similar with other similar studies [23, 25], as an important indicator of advanced TNM stages, LN metastasis was failed to indicate inferior survival in our study, which could partly explain why the loss of monotonous gradient for survival prediction of TNM stage and the superior predictive power of the established nomograms in our study.
With the increasing occurrence and concern of competing risk events, more and more focuses have been paid on competing analyses, such as lung cancer , breast cancer  and gastric cancer . Considering the non-cancer events contributed to 16.9% of deaths, competing interests were taken into account in survival analyses in this study. As far as we know, it was the first time to build prognostic nomograms to specially predict OS and CSS for patients with CHCC after surgery based on competing risk analysis. Significantly elevated predictive power was observed for the established nomograms in this study. The inclusion of additional variables guaranteed that nomograms were better in predicting OS and CSS, compared with the 8th edition TNM stage system. In addition, the nomograms were established based on a population-based dataset and cross-validated from an external dataset, making our results more generable than those from studies of small cohort or single center. Thus, a diverse range of parameters of CHCC patients are assessed by doctors more objectively and precisely based on the established nomograms. In addition, this newly established system can be used to identify subgroups of patients with a more homogeneous prognosis, estimate individual survival, and then to specialize personal treatment.
There were several limitations for this study. The major limitation of the present study is that not all risk factors were included to construct the nomograms. Some important tumor biomarker, such as carbohydrate antigen 19–9 (CA19–9), and some positive prognostic variables, such as surgical margin status and vascular invasion, were unavailable in SEER dataset. Maybe the additional inclusion of these variables might elevate the predictive power. This is also the major part of our future research. Another limitation is that although the established nomograms showed good discrimination and validation, the values of C-index and AUC are not relatively high. Further validation based on large-scale cohort is needed for these nomograms.
In conclusion, competing risk analyses were conducted and nomograms specially to predict OS and CSS for these patients were established for the first time in this study. The established nomograms can be used to accurately provide valuable prognostic information, allowing tailed treatments for patients with CHCC after surgery.
Area under ROC curve
Carbohydrate antigen 19–9
Combined hepatocellular cholangiocarcinoma
Cumulative incidence function
Receiver operating characteristic
Surveillance, Epidemiology, and End Results
Sun Yat-Sen University Cancer Center
Bosman FT. HH: WHO classifcation of tumours of the digestive system, 4th edn. Lyon: International Agency for Research on Cancer; 2010.
Kassahun WT, Hauss J. Management of combined hepatocellular and cholangiocarcinoma. Int J Clin Pract. 2008;62(8):1271–8.
Yin X, Zhang BH, Qiu SJ, Ren ZG, Zhou J, Chen XH, Zhou Y, Fan J. Combined hepatocellular carcinoma and cholangiocarcinoma: clinical features, treatment modalities, and prognosis. Ann Surg Oncol. 2012;19(9):2869–76.
Yap AQ, Chen CL, Yong CC, Kuo FY, Wang SH, Lin CC, Liu YW, Lin TL, Li WF, Millan CA, et al. Clinicopathological factors impact the survival outcome following the resection of combined hepatocellular carcinoma and cholangiocarcinoma. Surg Oncol. 2013;22(1):55–60.
Allen RA, Lisa JR. Combined liver cell and bile duct carcinoma. Am J Pathol. 1949;25(4):647–55.
Xu L, Peng ZW, Chen MS, Shi M, Zhang YJ, Guo RP, Lin XJ, Lau WY. Prognostic nomogram for patients with unresectable hepatocellular carcinoma after transcatheter arterial chemoembolization. J Hepatol. 2015;63(1):122–30.
He CB, Lin XJ. Inflammation scores predict the survival of patients with hepatocellular carcinoma who were treated with transarterial chemoembolization and recombinant human type-5 adenovirus H101. PLoS One. 2017;12(3):e0174769.
He CB, Lao XM, Lin XJ. Transarterial chemoembolization combined with recombinant human adenovirus type 5 H101 prolongs overall survival of patients with intermediate to advanced hepatocellular carcinoma: a prognostic nomogram study. Chinese J Cancer. 2017;36(1):59.
He C, Mao Y, Lao X, Li S, Lin X. Neutrophil-to-lymphocyte ratio predicts overall survival of patients with combined hepatocellular cholangiocarcinoma. Oncol Lett. 2018;15(4):4262–8.
He C, Mao Y, Wang J, Song Y, Huang X, Lin X, Li S. The predictive value of staging systems and inflammation scores for patients with combined hepatocellular cholangiocarcinoma after surgical resection: a retrospective study. J Gastrointest Surg. 2018;22(7):1239–50.
Tian MX, He WJ, Liu WR, Yin JC, Jin L, Tang Z, Jiang XF, Wang H, Zhou PY, Tao CY, et al. A novel risk prediction model for patients with combined hepatocellular-cholangiocarcinoma. J Cancer. 2018;9(6):1025–32.
Kim SH, Park YN, Lim JH, Choi GH, Choi JS, Kim KS. Characteristics of combined hepatocelluar-cholangiocarcinoma and comparison with intrahepatic cholangiocarcinoma. Eur J Surg Oncol. 2014;40(8):976–81.
Ma MC, Chen YJ, Chiu TJ, Lan J, Liu CT, Chen YC, Tien HH, Chen YY. Positive expression of Midkine predicts early recurrence and poor prognosis of initially resectable combined hepatocellular cholangiocarcinoma. BMC Cancer. 2018;18(1):227.
Schiergens TS, Lindenthaler A, Thomas MN, Rentsch M, Mittermeier L, Brand K, Kuchenhoff H, Lee S, Guba M, Werner J, et al. Time-dependent impact of age and comorbidities on long-term overall survival after liver resection. Liver Int. 2016;36(9):1340–50.
Szpakowski JL, Tucker LY. Causes of death in patients with hepatitis B: a natural history cohort study in the United States. Hepatology (Baltimore, Md). 2013;58(1):21–30.
Shindoh J, Kawamura Y, Kobayashi Y, Kiya Y, Sugawara T, Akuta N, Kobayashi M, Suzuki Y, Ikeda K, Hashimoto M. Platelet-albumin score as a sensitive measure for surgical risk prediction and survival outcomes of patients with hepatocellular carcinoma. J Gastrointest Surg. 2018;23(1):76–83.
Noordzij M, Leffondre K, van Stralen KJ, Zoccali C, Dekker FW, Jager KJ. When do we need competing risks methods for survival analysis in nephrology? Nephrol Dial Transplant. 2013;28(11):2670–7.
Harrell FE Jr, Lee KL, Mark DB. Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996;15(4):361–87.
Norris CM, Ghali WA, Saunders LD, Brant R, Galbraith D, Faris P, Knudtson ML. Ordinal regression model and the linear regression model were superior to the logistic regression models. J Clin Epidemiol. 2006;59(5):448–56.
Li D, Zhong C, Tang X, Yu L, Ding K, Yuan Y. Competing nomograms help in the selection of elderly patients with colon cancer for adjuvant chemotherapy. J Cancer Res Clin Oncol. 2018;144(5):909–23.
Zhou H, Zhang Y, Qiu Z, Chen G, Hong S, Chen X, Zhang Z, Huang Y, Zhang L. Nomogram to predict cause-specific mortality in patients with surgically resected stage I non-small-cell lung cancer: a competing risk analysis. Clin Lung Cancer. 2018;19(2):e195–203.
Chantajitr S, Wilasrusmee C, Lertsitichai P, Phromsopha N. Combined hepatocellular and cholangiocarcinoma: clinical features and prognostic study in a Thai population. J Hepato-Biliary-Pancreat Surg. 2006;13(6):537–42.
Bergquist JR, Groeschl RT, Ivanics T, Shubert CR, Habermann EB, Kendrick ML, Farnell MB, Nagorney DM, Truty MJ, Smoot RL. Mixed hepatocellular and cholangiocarcinoma: a rare tumor with a mix of parent phenotypic characteristics. HPB (Oxford). 2016;18(11):886–92.
Moeini A, Sia D, Zhang Z, Camprecios G, Stueck A, Dong H, Montal R, Torrens L, Martinez-Quetglas I, Fiel MI, et al. Mixed hepatocellular cholangiocarcinoma tumors: Cholangiolocellular carcinoma is a distinct molecular entity. J Hepatol. 2017;66(5):952–61.
Chu KJ, Lu CD, Dong H, Fu XH, Zhang HW, Yao XP. Hepatitis B virus-related combined hepatocellular-cholangiocarcinoma: clinicopathological and prognostic analysis of 390 cases. Eur J Gastroenterol Hepatol. 2014;26(2):192–9.
Lian M, Perez M, Liu Y, Schootman M, Frisse A, Foldes E, Jeffe DB. Neighborhood socioeconomic deprivation, tumor subtypes, and causes of death after non-metastatic invasive breast cancer diagnosis: a multilevel competing-risk analysis. Breast Cancer Res Treat. 2014;147(3):661–70.
Shin DW, Suh B, Park Y, Lim H, Suh YS, Yun JM, Cho B, Yang HK. Risk of coronary heart disease and ischemic stroke incidence in gastric cancer survivors: a Nationwide study in Korea. Ann Surg Oncol. 2018;25(11):3248–56.
The authors would like to thank SEER for open access to the database.
This work was supported by the Sun Yat-sen University Grant for Medical Humanities Practice and Teaching (No. 23000–18008023). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in the writing of the manuscript.
Availability of data and materials
The dataset from SEER database generated and/or analyzed during the current study are available in the SEER dataset repository (https://seer.cancer.gov/). Another dataset from SYSUCC dataset has been validated by uploading the key raw data onto the Research Data Deposit public platform (http://www.researchdata.org.cn), with the Approval Number as RDDA2018000870.
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of Sun Yat-sen University Cancer Center. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Written informed consent was obtained from individual participants included in the study. All authors signed authorization forms and received permission from SEER to access and use the dataset.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
He, C., Zhang, Y., Cai, Z. et al. Competing risk analyses of overall survival and cancer-specific survival in patients with combined hepatocellular cholangiocarcinoma after surgery. BMC Cancer 19, 178 (2019) doi:10.1186/s12885-019-5398-6
- Combined hepatocellular cholangiocarcinoma
- Overall survival
- Cancer-specific survival