- Research article
- Open Access
Prediction for early recurrence of intrahepatic mass-forming cholangiocarcinoma: quantitative magnetic resonance imaging combined with prognostic immunohistochemical markers
Cancer Imagingvolume 19, Article number: 49 (2019)
Partial hepatectomy is the first option for intrahepatic mass-forming cholangiocarcinoma (IMCC) treatment, which would prolong survival. The main reason for the poor outcome after curative resection is the high incidence of early recurrence (ER). The aim of this study was to investigate the combined predictive performance of qualitative and quantitative magnetic resonance imaging (MRI) features and prognostic immunohistochemical markers for the ER of IMCC.
Forty-seven patients with pathologically proven IMCC were enrolled in this retrospective study. Preoperative contrast-enhanced MRI and post-operative immunohistochemical staining of epidermal growth factor receptors (EGFR), vascular endothelial growth factor receptor (VEGFR), P53 and Ki67 were performed. Univariate analysis identified clinic-radiologic and pathological risk factors of ER. Radiomics analysis was performed based on four MRI sequences including fat suppression T2-weighted imaging (T2WI/FS), arterial phase (AP), portal venous phase (PVP), and delayed phase (DP) contrast enhanced imaging. A clinicoradiologic-pathological (CRP) model, radiomics model, and combined model were developed. And ROC curves were used to explore their predictive performance for ER stratification.
Enhancement patterns and VEGFR showed significant differences between the ER group and non-ER group (P = 0.001 and 0.034, respectively). The radiomics model based on AP, PVP and DP images presented superior AUC (0.889, 95% confidence interval (CI): 0.783–0.996) among seven radiomics models with a sensitivity of 0.938 and specificity of 0.839. The combined model, containing enhancement patterns, VEGFR and radiomics features, showed a preferable ER predictive performance compared to the radiomics model or CRP model alone, with AUC, sensitivity and specificity of 0.949, 0.875 and 0.774, respectively.
The combined model was the superior predictive model of ER. Combining qualitative and quantitative MRI features and VEGFR enables ER prediction, thus facilitating personalized treatment for patients with IMCC.
Intrahepatic cholangiocarcinoma (ICC) is the second most common primary liver malignancy and is arising in incidence worldwide [1, 2]. It originates from the intrahepatic biliary epithelium and can be classified into three types according to the morphologic classification system: mass forming, periductal infiltrating, and intraductal growing . Intrahepatic mass-forming cholangiocarcinoma (IMCC) accounts for a large percentage of ICC. Partial hepatectomy is the first option for IMCC curative treatment, which would prolong survival . However, even after a curative resection, the 5-year survival rate is only 20–35% [5, 6]. The main reason for the poor outcome is the incidence of recurrence, which can be as high as 54–71% [5, 7, 8].
The time interval from the resection to IMCC recurrence is an independent prognostic factor of survival . Approximately 78.8% of recurrence develops within 24 months, defined as early recurrence (ER), and the prognosis for patients with ER is worse than that for those with late recurrence. Adjuvant trans-arterial chemoembolization (TACE) or chemotherapy after surgery was associated with better survival among the IMCC patients with early recurrence [10, 11]. Therefore, patients at a high risk of ER need to be precisely determined and effective adjunctive treatment strategies and closer follow-up after operation need to be performed.
Previous studies have revealed several pathological tumor characteristics associated with postoperative ER of IMCC (e.g. tumor size, satellite lesions, lymph node metastasis, lymphatic invasion, microvascular invasion and stage) . In addition, certain immunohistochemical molecules have been reported as predictive markers of ICC prognosis. Iguchi et al. found that P53 and Ki67, markers indicating cell proliferation, were related to the overall survival of ICC patients . Epidermal growth factor receptor (EGFR) expression was demonstrated to be an independent predictor of ICC prognosis . As an anti-angiogenesis therapeutic target, vascular endothelial growth factor receptor (VEGFR) expression has been correlated with the prognosis of many cancers (e.g. breast cancer, ovarian cancer, lung cancer, lymphoma, etc.) [14,15,16,17]. However, ER prediction by these factors leads to noticeably different outcomes among studies. Therefore, it is necessary to further explore whether they are effective predictors of ER in IMCC.
Magnetic resonance imaging (MRI) is widely used in the diagnosis and treatment planning of liver tumors. Previous studies revealed some radiological features that might be predictors of IMCC prognosis, such as the degree of diffusion restriction on diffusion weighted images (DWI), enhancement pattern of contrast enhanced MRI (CE-MRI), and intensity on the hepatobiliary phase of gadoxetic acid-enhanced MRI [18,19,20]. However, these qualitative assessments were unable to quantify tumor heterogeneity. Radiomics is considered to be an emerging quantitative technique for evaluating the entire underlying intra-tumor heterogeneity by extracting numerous features from radiologic images. Studies in many cancers (e.g. colorectal cancer, breast cancer, lung cancer, esophageal cancer and hepatocellular carcinoma) have shown that radiomics has the potential for prognosis prediction [21,22,23,24,25].
We combined the above pathological characteristics and immunohistochemical molecules with both visible and invisible radiological features for better ER prediction. To our knowledge, combining immunohistochemistry and radiomics features for ER prediction in IMCC patients has not yet been investigated. Therefore, the aim of this study was to develop a nomogram based on radiological features, immunohistochemical markers, and radiomics features for predicting the ER of IMCC to create a better stratification of IMCC patients thus improving personalized treatment.
Materials and methods
This retrospective study was approved by our Institutional Review Board, and the need for informed consent was waived. We collected electronic medical records from May 2011 to August 2016 at our institution. A total of 219 patients who underwent a curative-intent resection and lymph node dissection, with histopathologically confirmed IMCC, were recruited. Among these patients, the study population was selected using the following inclusion criteria: (a) patients who underwent preoperative liver CE-MRI within 4 weeks of their surgery, (b) patients without a history of previous adjuvant treatment before the surgery, (c) patients with histopathologically proven IMCC and negative resection margin (R0), combined hepatocellular-cholangiocarcinoma were excluded, (d) patients without a history of other tumors, and (e) patients who completed at least 2 years of follow-up. Consequently, 47 patients were enrolled in our study and were divided into an ER group (n = 31) and a non-ER group (n = 16), with ER being defined as the development of intrahepatic or extrahepatic recurrence within 2 postoperative years (Fig. 1).
Clinical and pathologic characteristics
Clinical and pathological characteristics consisted of age, gender, hepatitis, carcinoembryonic antigen (CEA), carbohydrate antigen 199 (CA199), satellite lesions, maximum tumor diameter (MTD), tumor location, differential degree of tumor, stage and lymph node metastases. Lymph node metastasis was determined by the pathological results of lymph node dissection and preoperative CT/MR imaging. The threshold values chosen for CA199 and CEA levels were based on the normal ranges used at our institution (0-37 U/ml for CA199 and 0-5 ng/ml for CEA).
All pathology sections and macroscopic pictures of the resected specimen were retrospectively reviewed. Combining with immunohistochemistry, all 47 patients were confirmed IMCC by post-operative pathology. Immunohistochemical staining was detected on formalin-fixed, paraffin-embedded sections using standard immunohistochemical methods. Then five-micron-thick sections were created, and antibodies specific for EGFR, VEGFR, P53, and Ki67 (Beijing Zhongshan Golden Bridge Biotechnology Co. LTD, China) were used to perform the further immunohistochemical staining. All samples were analyzed by an anatomic pathologist with 10 years of experience, who was unaware of the patient’s outcome. Less than 10% of the positive staining was identified as negative expression, while more than 10% of the positive staining was identified as positive expression.
All patients underwent contrast enhanced CT or MRI every 3–6 months after surgery in the first 2 years. Images were analyzed to identify ER, which was determined as the presence of new intrahepatic lesions with typical imaging features of IMCC, atypical lesions with histopathological confirmation, or extrahepatic metastasis (lymph node metastases or distant metastasis) within 2 postoperative years.
Magnetic resonance imaging acquisition and analysis
All patients underwent 3.0 T MRI scans (Signa Excite HDxt, GE Healthcare, Milwaukee, USA) with an eight-element phased-array torso coil. After nonenhanced axial breath-hold T1 weighted imaging, fat suppression T2 weighted imaging (T2WI/FS) and DWI (b-values of 0 and 800 s/mm2), contrast enhanced T1-weighted three-dimensional (3D) spoiled gradient echo sequence (LAVA) was performed. Gadodiamide (Omniscan 0.5 mmol/ml; GE Healthcare, Ireland) was injected at a dose of 0.2 mL per kilogram and a rate of 2 mL per second as a bolus by an automatic pump injector and a subsequent 20 mL 0.9% sterile saline flush. Contrast enhanced imaging was performed in the arterial phase (AP) (30 s), portal venous phase (PVP) (60 s), and delayed phase (DP) (180 s). Images of AP, PVP, and DP were all used for radiomics feature extraction and analysis. Additional technical details are provided in Table 1.
MR images were reviewed independently by two radiologists who had 5 and 10 years experience in the abdominal MRI, respectively. Both radiologists were blinded to the clinical data of the patients when they evaluated the MR images. They reached a consensus by discussion when there were disagreements. The basic imaging traits potentially associated with ER included lesion shape, contour, biliary dilation, capsular retraction, DWI intensity and the enhancement pattern. The target appearance was peripheral hyperintensity compared to the center on high b-value DWI. The enhancement pattern was assessed using the following subdivisions: (a) gradual enhancement (the enhancement area gradually increased from the periphery to the center of the tumor), (b) persistent enhancement (enhancement remained through all phases), (c) wash in and wash out (hyperenhancement of the AP followed by washout), and (d) minimal or no enhancement.
Radiomics features extraction and analysis
MR images (AP, PVP, DP, T2WI/FS sequence) were loaded into ITK-SNAP software (version 2.2.0, www.itksnap.org) for 3D manual segmentation. A radiologist with 10 years of MRI experience (reader 1) performed the tumor segmentations in all 47 patients. After 2 weeks, images of all patients were segmented again by reader 1 and another radiologist (reader 2) with 5 years of experience of MRI diagnosis to assess intra−/inter-reader agreement in the feature analysis. All outcomes were based on the features extracted by the first segmentation from reader 1.
Artificial Intelligence Kit software (A.K. software; GE Healthcare, Life Sciences, Beijing, China) was used to extract 396 parameters from each sequence. Those parameters include first order histogram features (n = 42), grey-level co-occurrence matrix (GLCM) features (n = 144), grey-level run-length matrix features (n = 180), Haralick features (n = 10), morphological features (n = 9) and grey-level zone size matrix features (n = 11).
The proposed parameters were analyzed for consistency and correlation. First, the intraclass correlation coefficient was determined for each parameter for the inter-observer and intra-observer reproducibility test. Features with intraclass correlation coefficient values less than 0.8 were excluded. Second, stratified analyses were conducted using the Wilcoxon signed-rank test to discover the potential association between the remaining parameters and ER status, followed by univariate logistic regression. To keep discriminative parameters, we set a threshold of 0.1. A variance inflation factor was then used to eliminate parameters with high collinearity in a multiple mutual linear situation. Finally, multivariate logistic regression was applied to evaluate the performance of distinguishing ER status in each sequence. Different sequence combinations (AP + PVP, AP + PVP + DP, AP + PVP + DP + T2WI) were also tried to explore the best model using the same methods described above. The validation of each model was performed by using leave-one-out cross-validation. Thus, three predictive models were built: a best performance radiomics model, a clinicoradiologic-pathologic(CRP) model, and a combined model with both selected radiomics features and clinicoradiologic-pathologic features (Fig. 2).
The differences in patient characteristics between the two groups were assessed using t-test or the Mann-Whitney U test for continuous variables and the chi-square test or Fisher exact test for categorical variables. Kappa tests were used to determine inter-observer agreement for qualitative MRI features. Kappa values of 0.81–1.00 indicated excellent agreement, 0.61–0.80 signified substantial agreement, and 0.41–0.60 denoted moderate agreement. Receiver operating characteristics (ROC) curves were performed in each model. The area under the curve (AUC) of the ROC curves, accuracy, sensitivity, specificity, positive predictive values (PPVs) and negative predictive values (NPVs) were obtained and comparisons between the performance of final three models were performed using the Delong test. All statistical analyses were performed using R studio Server (Version 3.5.0; RStudio, Inc., Boston, MA, USA) and SPSS (version 20.0; IBM, Armonk, NY, USA). A two-sided P value less than 0.05 indicated a statistically significant difference.
The median age of the patients was 57 years old (range, 35 to 78 years old), and ER was found in 31 (66.0%) patients. The median MTD was 5.73 cm, ranging from 1.5 to 12.8 cm. The median follow-up was 34 months (range, 25 to 87 months). Clinical and pathological characteristics in the ER group and non-ER group are summarized in Table 2. There were no significant differences between the ER group and non-ER group in terms of these characteristics. Radiological features and immunohistochemical markers are listed in Table 3. Among these factors, enhancement patters and VEGFR expression showed significant differences between the ER group and non-ER group (P = 0.001 and 0.034, respectively). Additionally, among 33 IMCC patients with gradual enhancement pattern, 26 (78.8%) patients were in the ER group, while all 5 (100%) patients with wash in and wash out enhancement pattern were in non-ER group (P = 0.002). The inter-observer agreement for the radiological features showed excellent agreement (k = 0.811–0.849).
We began by developing radiomics models based on T2WI/FS, AP, PVP, and DP images separately. The PVP model demonstrated preferable accuracy, sensitivity, specificity, PPV and NPV (0.872, 0.75, 0.936, 0.857, and 0.879, respectively) while it presented a slightly lower AUC (0.841, 95% confidence interval (CI): 0.697–0.984) than that of the AP model (0.871, 95% CI: 0.761–0.981). Next, radiomics models based on multiple sequences were built according to the aforementioned results, including AP + PVP (two sequences with higher AUC), AP + PVP + DP (multi-phase contrast enhanced sequences), and T2WI/FS + AP + PVP + DP (all sequences) models. The AP + PVP + DP model showed superior AUC (0.889, 95% CI: 0.783–0.996) among all the radiomics models, and was used in the follow-up study. This model illustrated that the four most important parameters for predicting ER were AP_skewness and PVP_Variance, both derived from the histogram, as well as AP_ClusterShade_AllDirection_offset7_SD and AP_GLCMEntropy_angle45_offset7 derived from the GLCM. The accuracy, sensitivity, specificity, PPV, NPV and AUC of the seven radiomics models are presented in Table 4.
The clinicoradiologic-pathologic (CRP) model contained enhancement pattern and VEGFR. The combined model incorporated radiomics features, clinicoradiological features and pathological factors. The predictive performance of the CRP model, radiomics model, and combined model are listed in Table 5, and ROCs are shown in Fig. 3. Nomograms for the combined model are presented in Fig. 4. The combined model displayed the best accuracy, sensitivity, specificity, PPV, NPV and AUC (0.872, 0.938, 0.839, 0.750, 0.963 and 0.949, respectively). Also, the combined model significantly improved the predictive performance of the CRP model in predicting ER of IMCC (P = 0.009).
Combining radiomics features, enhancement patterns, and VEGFR led to significant improvements in the AUC, sensitivity, specificity and accuracy for predicting ER compared to the radiomics model or CRP model alone, which indicated that the combination of qualitative and quantitative MRI features along with immunohistochemical markers maximizes the predictive performance of ER. Radiomics models based on CE-MRI sequences (AP, PVP, or DP) showed better specificity and AUC for predicting ER than that of T2WI, though they exhibited relatively lower sensitivity.
Previous studies have investigated ER predictions of IMCC. For instance, Liang et al. developed a novel nomogram to predict the recurrence of ICC with AUC, sensitivity, and specificity values of 0.90, 0.74, 0.89, respectively . This nomogram was achieved based on radiomics and clinical stage, and radiomic signatures extracted only from AP images with a relatively low sensitivity. Jeong et al. established a predictive nomogram of IMCC recurrence based only on clinical characteristics: lymph node metastasis, tumor size, surface antigen of the hepatitis B virus, and Child–Pugh score, with a concordance C index of 0.71 (95% CI: 0.65–0.77) . Our study demonstrated preferable sensitivity and AUC compared with these studies. Furthermore, our combined model that included radiomics features, clinicoradiological and pathological factors was the superior predictive model. This suggested that combining the morphology, quantification of tumor heterogeneity and molecular pathology could better reflect aggressive malignant tumor biology.
There have been an increasing number of studies showing the potential of radiomics based on MR images for diagnosis and prognosis assessment for specific tumors. In our study, the radiomics features were extracted from T2WI/FS, AP, PVP and DP of contrast enhanced images to build the best radiomics model. Radiomics models of AP, PVP and DP all provided better AUC and specificity for predicting ER than that of T2WI/FS, which suggested that CE-MRI contain more potential tumor heterogeneous information. Catharina et al. also revealed the exceptional discriminating ability provided by CE-MRI among T1WI, T2WI, T2WI/short inversion time recovery, and contrast enhanced sequences in differentiating low-grade chondrosarcoma and enchondroma by texture analysis (TA) . Furthermore, the AP + PVP + DP model showed superior AUC among all the TA models, which was consistent with previous studies. Ueno et al.  found that a TA model based on texture parameters with several sequences led to a better predictive value than that with a single sequence, demonstrating that multi-phase CE-MRI could provide added value.
We found that VEGFR was a predictor of ER, which was consistent with previous studies. Sang et al. found that ICC with positive VEGFR expression represented aggressive malignancy owing to the mechanism that inhibition of VEGFR-2 expression increased apoptosis and decreased cell proliferation . The enhancement patterns were also related to ER due to the histopathologic basis of fibrous stroma in tumors. Gradual enhancement pattern was relevant to the large amount of fibrous stroma of the tumor and indicated a poor prognosis . Small IMCC with diameter less than 3 cm in the cirrhotic liver showed atypical wash in and wash out enhancement pattern frequently . IMCC with wash in and wash out enhancement pattern demonstrated less central fibrous stroma and more cellular areas than that with gradual enhancement pattern; also, hyperenhancement on AP was an independent factor for longer survival.
Whether MTD and CA199 levels could serve as a predictor of ER is currently controversial, perhaps due to the heterogenous population of IMCC patients. Several studies reported that no associations were found between MTD and IMCC prognosis . Nevertheless, other studies found that MTD was associated with ER [27, 34]. Similarly, a few previous studies found CA199 to be a preoperative predictor of prognosis [34, 35], while others removed CA199 as an independent prognostic factor of ER . In our study, MTD and CA199 levels displayed no significant correlation with ER.
Our study had some limitations. First, since it was a retrospective study and was performed in a single center, thus lacking the heterogeneity of MR images and the cohort population of other institutions, selection bias may exist. We will explore the prediction model using MR images from multiple centers in the future. Other limitations of this study were the relatively small sample size cohort due to the incidence of IMCC and the need to obtain pathological sections for immunohistochemistry. Finally, we only developed predictive models for ER without including long-term survival analysis. Prediction of long-term survival should be included in future studies.
Our study show that the combined model was the superior predictive model of ER compared with radiomics or CRP model alone. Combining qualitative and quantitative MRI features and VEGFR might be useful for predicting ER and guide personalized treatment in patients with IMCC.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Carbohydrate antigen 199
Contrast enhanced MRI
Epidermal growth factor receptor
Grey-level co-occurrence matrix
Intrahepatic mass-forming cholangiocarcinoma
Maximum tumor diameter
Portal venous phase
Region of interest
Vascular endothelial growth factor receptor
Razumilava N, Gores GJ. Cholangiocarcinoma. Lancet. 2014;383(9935):2168–79.
Bergquist A, von Seth E. Epidemiology of cholangiocarcinoma. Best Pract Res Clin Gastroenterol. 2015;29(2):221–32.
Chung YE, Kim MJ, Park YN, Choi JY, Pyo JY, Kim YC, et al. Varying appearances of cholangiocarcinoma: radiologic-pathologic correlation. Radiographics. 2009;29(3):683–700.
Amini N, Ejaz A, Spolverato G, Kim Y, Herman JM, Pawlik TM. Temporal trends in liver-directed therapy of patients with intrahepatic cholangiocarcinoma in the United States: a population-based analysis. J Surg Oncol. 2014;110(2):163–70.
Spolverato G, Kim Y, Alexandrescu S, Marques HP, Lamelas J, Aldrighetti L, et al. Management and outcomes of patients with recurrent intrahepatic cholangiocarcinoma following previous curative-intent surgical resection. Ann Surg Oncol. 2016;23(1):235–43.
Mavros MN, Economopoulos KP, Alexiou VG, Pawlik TM. Treatment and prognosis for patients with intrahepatic cholangiocarcinoma: systematic review and meta-analysis. Jama Surg. 2014;149(6):565–74.
Hyder O, Hatzaras I, Sotiropoulos GC, Paul A, Alexandrescu S, Marques H, et al. Recurrence after operative management of intrahepatic cholangiocarcinoma. Surgery. 2013;153(6):811–8.
Sulpice L, Rayar M, Boucher E, Pracht M, Meunier B, Boudjema K. Treatment of recurrent intrahepatic cholangiocarcinoma. Br J Surg. 2012;99(12):1711–7.
Zhang XF, Beal EW, Bagante F, Chakedis J, Weiss M, Popescu I, et al. Early versus late recurrence of intrahepatic cholangiocarcinoma after resection with curative intent. Br J Surg. 2018;105(7):848–56.
Yamashita YI, Shirabe K, Beppu T, Eguchi S, Nanashima A, Ohta M, et al. Surgical management of recurrent intrahepatic cholangiocarcinoma: predictors, adjuvant chemotherapy, and surgical therapy for recurrence: a multi-institutional study by the Kyushu study Group of Liver Surgery. Ann Gastroenterol Surg. 2017;1(2):136–42.
Li J, Wang Q, Lei Z, Wu D, Si A, Wang K, et al. Adjuvant Transarterial chemoembolization following liver resection for intrahepatic cholangiocarcinoma based on survival risk stratification. Oncologist. 2015;20(6):640–7.
Iguchi T, Yamashita N, Aishima S, Kuroda Y, Terashi T, Sugimachi K, et al. A comprehensive analysis of immunohistochemical studies in intrahepatic cholangiocarcinoma using the survival tree model. Oncology. 2009;76(4):293–300.
Yoshikawa D, Ojima H, Iwasaki M, Hiraoka N, Kosuge T, Kasai S, et al. Clinicopathological and prognostic significance of EGFR, VEGF, and HER2 expression in cholangiocarcinoma. Br J Cancer. 2008;98(2):418–25.
Santos A, Lopes C, Gartner F, Matos AJ. VEGFR-2 expression in malignant tumours of the canine mammary gland: a prospective survival study. Vet Comp Oncol. 2016;14(3):e83–92.
Sallinen H, Heikura T, Koponen J, Kosma VM, Heinonen S, Yla-Herttuala S, et al. Serum angiopoietin-2 and soluble VEGFR-2 levels predict malignancy of ovarian neoplasm and poor prognosis in epithelial ovarian cancer. BMC Cancer. 2014;14:696.
Song DH, Ko GH, Lee JH, Lee JS, Lee GW, Kim HC, et al. Myoferlin expression in non-small cell lung cancer: prognostic role and correlation with VEGFR-2 expression. Oncol Lett. 2016;11(2):998–1006.
Dimtsas GS, Georgiadi EC, Karakitsos P, Vassilakopoulos TP, Thymara I, Korkolopoulou P, et al. Prognostic significance of immunohistochemical expression of the angiogenic molecules vascular endothelial growth factor-a, vascular endothelial growth factor receptor-1 and vascular endothelial growth factor receptor-2 in patients with classical Hodgkin lymphoma. Leuk Lymphoma. 2014;55(3):558–64.
Lee J, Kim SH, Kang TW, Song KD, Choi D, Jang KT. Mass-forming intrahepatic cholangiocarcinoma: diffusion-weighted imaging as a preoperative prognostic marker. Radiology. 2016;281(1):119–28.
Kang Y, Lee JM, Kim SH, Han JK, Choi BI. Intrahepatic mass-forming cholangiocarcinoma: enhancement patterns on gadoxetic acid-enhanced MR images. Radiology. 2012;264(3):751–60.
Kim SA, Lee JM, Lee KB, Kim SH, Yoon SH, Han JK, et al. Intrahepatic mass-forming cholangiocarcinomas: enhancement patterns at multiphasic CT, with special emphasis on arterial enhancement pattern--correlation with clinicopathologic findings. Radiology. 2011;260(1):148–57.
Liang W, Yang P, Huang R, Xu L, Wang J, Liu W, et al. A combined nomogram model to preoperatively predict histologic grade in pancreatic neuroendocrine tumors. Clin Cancer Res. 2019;25(2):584–94.
Ji GW, Zhang YD, Zhang H, Zhu FP, Wang K, Xia YX, et al. Biliary tract Cancer at CT: a Radiomics-based model to predict lymph node metastasis and survival outcomes. Radiology. 2019;290(1):90–8.
Cui Y, Yang X, Shi Z, Yang Z, Du X, Zhao Z, et al. Radiomics analysis of multiparametric MRI for prediction of pathological complete response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer. Eur Radiol. 2019;29(3):1211–20.
Yuan C, Wang Z, Gu D, Tian J, Zhao P, Wei J, et al. Prediction early recurrence of hepatocellular carcinoma eligible for curative ablation using a Radiomics nomogram. Cancer Imaging. 2019;19(1):21.
Park H, Lim Y, Ko ES, Cho HH, Lee JE, Han BK, et al. Radiomics signature on magnetic resonance imaging: association with disease-free survival in patients with invasive breast Cancer. Clin Cancer Res. 2018;24(19):4705–14.
Liang W, Xu L, Yang P, Zhang L, Wan D, Huang Q, et al. Novel nomogram for preoperative prediction of early recurrence in intrahepatic cholangiocarcinoma. Front Oncol. 2018;8:360.
Jeong S, Cheng Q, Huang L, Wang J, Sha M, Tong Y, et al. Risk stratification system to predict recurrence of intrahepatic cholangiocarcinoma after hepatic resection. BMC Cancer. 2017;17(1):464.
Lisson CS, Lisson CG, Flosdorf K, Mayer-Steinacker R, Schultheiss M, von Baer A, et al. Diagnostic value of MRI-based 3D texture analysis for tissue characterisation and discrimination of low-grade chondrosarcoma from enchondroma: a pilot study. Eur Radiol. 2018;28(2):468–77.
Ueno Y, Forghani B, Forghani R, Dohan A, Zeng XZ. Chamming'S F, et al. endometrial carcinoma: MR imaging-based texture model for preoperative risk stratification-a preliminary analysis. Radiology. 2017;284(3):748–57.
Sang H, Li T, Li H, Liu J. Gab1 regulates proliferation and migration through the PI3K/Akt signaling pathway in intrahepatic cholangiocarcinoma. Tumour Biol. 2015;36(11):8367–77.
Asayama Y, Yoshimitsu K, Irie H, Tajima T, Nishie A, Hirakawa M, et al. Delayed-phase dynamic CT enhancement as a prognostic factor for mass-forming intrahepatic cholangiocarcinoma. Radiology. 2006;238(1):150–5.
Mendiratta-Lala M, Park H, Kolicaj N, Mendiratta V, Bassi D. Small intrahepatic peripheral cholangiocarcinomas as mimics of hepatocellular carcinoma in multiphasic CT. Abdom Radiol (NY). 2017;42(1):171–8.
Jiang W, Zeng ZC, Tang ZY, Fan J, Sun HC, Zhou J, et al. A prognostic scoring system based on clinical features of intrahepatic cholangiocarcinoma: the Fudan score. Ann Oncol. 2011;22(7):1644–52.
Wang Y, Li J, Xia Y, Gong R, Wang K, Yan Z, et al. Prognostic nomogram for intrahepatic cholangiocarcinoma after partial hepatectomy. J Clin Oncol. 2013;31(9):1188–95.
Shen WF, Zhong W, Liu Q, Sui CJ, Huang YQ, Yang JM. Adjuvant transcatheter arterial chemoembolization for intrahepatic cholangiocarcinoma after curative surgery: retrospective control study. World J Surg. 2011;35(9):2083–91.
This study has received funding by CAMS Innovation Fund for Medical Sciences (CIFMS) (grant number 2016-I2M-1-001), PUMC Youth Fund (grant number 2017320010), and Beijing Hope Run Special Fund of Cancer Foundation of China (grant number LC2016B15).
Ethics approval and consent to participate
This retrospective single-center study was approved by our Institutional Review Board, and the need for informed consent was waived.
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.