- 1Department of Hepatobiliary and Pancreatic Surgery, The Center of Integrated Oncology and Precision Medicine, Affiliated Hangzhou First People’s Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 2Department of Hepatobiliary and Pancreatic Surgery, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 3Institute of Translational Medicine, Zhejiang University, Hangzhou, China
- 4Department of Radiology, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 5Department of Radiology, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 6Department of General Surgery, Lishui Central Hospital, Lishui, China
- 7Department of Hepatobiliary and Pancreatic Surgery, Shulan Health Hangzhou Hospital, Hangzhou, China
- 8Nucelar & Radiological Engineering and Medical Physics Programs, Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA, United States
- 9Zhejiang University Cancer Center, Hangzhou, China
- 10NHC Key Laboratory of Combined Multi-Organ Transplantation, Hangzhou, China
- 11Institute of Organ Transplantation, Zhejiang University, Hangzhou, China
Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer and has poor prognosis. Cytokeratin (CK)19-positive (CK19+) HCC is especially aggressive; early identification of this subtype and timely intervention can potentially improve clinical outcomes. In the present study, we developed a preoperative gadoxetic acid-enhanced magnetic resonance imaging (MRI)-based radiomics model for noninvasive and accurate classification of CK19+ HCC. A multicenter and time-independent cohort of 257 patients were retrospectively enrolled (training cohort, n = 143; validation cohort A, n = 75; validation cohort B, n = 39). A total of 968 radiomics features were extracted from preoperative multisequence MR images. The maximum relevance minimum redundancy algorithm was applied for feature selection. Multiple logistic regression, support vector machine, random forest, and artificial neural network (ANN) algorithms were used to construct the radiomics model, and the area under the receiver operating characteristic (AUROC) curve was used to evaluate the diagnostic performance of corresponding classifiers. The incidence of CK19+ HCC was significantly higher in male patients. The ANN-derived combined classifier comprising 12 optimal radiomics features showed the best diagnostic performance, with AUROCs of 0.857, 0.726, and 0.790 in the training cohort and validation cohorts A and B, respectively. The combined model based on multisequence MRI radiomics features can be used for preoperative noninvasive and accurate classification of CK19+ HCC, so that personalized management strategies can be developed.
Introduction
Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer and the third leading cause of cancer-related mortality worldwide (1). Although significant progress has been made in treatment strategies, HCC has poor prognosis due to the heterogeneity of the disease and high recurrence rate after treatment (1). Clarifying the tumor biology of HCC can lead to the development of more effective and personalized management strategies.
Hepatic progenitor cells (HPCs) are bipotent liver stem cells that can differentiate into hepatocytes or cholangiocytes (2–5). HPCs express specific markers such as epithelial cell adhesion molecular (EpCAM) and cytokeratin (CK)19 and have been shown to transform into hepatic cancer stem cells that give rise to HCC (2). In vitro studies have demonstrated that CK19-positive (CK19+) HCC cell lines exhibit high invasive potential, enhanced epithelial–mesenchymal transition (EMT), and the ability to induce angiogenesis (6–8). In the clinical context, CK19+ HCC is highly invasive and resistant to chemoradiotherapy and shows high rates of lymph node metastasis and early recurrence after hepatic resection or liver transplantation (9–11). Thus, CK19+ HCC is a unique disease subtype for which more effective treatments are needed (12).
Pretreatment diagnosis of CK19+ HCC currently relies on immunohistochemical analysis of fine-needle aspiration (FNA) samples. However, as an invasive method, FNA can cause abdominal hemorrhage as well as tumor rupture and peritoneal seeding. Additionally, the small quantities of sample obtained by FNA can lead to misdiagnosis. Radiomics—which refers to the high-throughput extraction of quantitative features of medical images and their analysis and use for model construction—has been widely used in oncology research (13–17). Radiomics can fully capture tumor heterogeneity and subtle changes in texture from medical images. CK19+ HCC exhibits unique features in medical images such as arterial rim enhancement and irregular tumor margins (18–21) that are not observed in CK19-negative (CK19−) tumors. There have been attempts made to develop classifiers using radiomics approaches for the determination of CK19 status (22–24). However, these studies were conducted at a single center or were based on a small sample, thereby limiting the generalizability of the findings.
In this study, we constructed and validated a multisequence radiomics model for noninvasive and accurate identification of CK19 status in HCC patients based on multicenter and time-independent magnetic resonance imaging (MRI) data.
Materials and Methods
Study Population
The inclusion criteria were as follows: i) patients who underwent gadoxetic acid-enhanced MRI within 2 weeks before surgery and ii) patients who had a postsurgery diagnosis of HCC with definite CK19 status (+/−). The exclusion criteria were as follows: i) concurrent malignancies or distal metastasis; ii) clinical data or MR images were missing or of low quality; and iii) small tumors (maximum diameter <1 cm). The patient enrollment process is illustrated in Supplementary Figure S1.
The study included patients from the following three independent medical centers: 1) First Affiliated Hospital of Zhejiang University School of Medicine (FAH-ZJU; patients treated between January 2016 and December 2017); 2) Shulan Health (Hangzhou) Hospital (SLH; patients treated between January 2017 and December 2018); and 3) Lishui Central Hospital (LSCH; patients treated between January 2018 and December 2018). FAH-ZJU, from which the majority of patients were enrolled, was set as the training cohort, and SLH and LSCH were the independent external validation cohorts.
The following baseline (presurgery) patient data were collected: age, sex, tumor number, tumor size, hepatitis B virus infection, serum tumor markers, and serum liver and kidney function indices. Histopathologic data included liver cirrhosis, tumor differentiation grade, microvascular invasion, macrovascular invasion, and lymph node metastasis.
Study Design
A flow diagram of the study is shown in Figure 1. The study was divided into four parts: 1) patient enrollment; 2) image preprocessing and segmentation; 3) multisequence MRI radiomics feature extraction and model construction; and 4) model validation. The delineation of the region of interest (ROI) is shown in Figure 2.
Figure 1 Overview of the proposed study design. In general, this study contains four parts: (1) retrospective collection of presurgery gadoxetic acid-enhanced MRI; (2) tumor segmentation; (3) feature extraction and subsequent ranking; (4) predictive classifier construction and independent validation.
Figure 2 Delineation of the region of interest (ROI). Immunohistochemistry (IHC) staining-based grouping of cytokeratin (CK)19 expression [(A), CK19 positive; (D), CK19 negative]. Undelineated and delineated image of contrast-enhanced T2-weighted imaging [(B), CK19 positive; (E), CK19 negative]. Undelineated and delineated image of contrast-enhanced diffusion-weighted imaging [(C), CK19 positive; (F), CK19 negative].
MR Image Acquisition
Presurgery gadoxetic acid-enhanced MR images were retrieved from the Picture Archiving and Communication Systems of the participating centers. Patients fasted for 4–6 h before MR scanning. The scan was performed after injection of the contrast agent (0.025 mmol/kg gadoxetic acid, 15 ml, 7.04 g; Guangzhou Consun Pharmaceutical Co., Guangzhou, China) at a rate of 2 ml/s via an injector. The contrast agent was flushed with 20 ml saline injected at the same rate. The scanners used at the three medical centers were as follows: FAH-ZJU, 3.0-T MR scanner (Signa HDxt; GE Healthcare, Milwaukee, WI, USA); SLH, 1.5-T or 3.0-T MR scanner (Signa HDxt); and LSCH, 1.5-T MR scanner (Aera; Siemens Healthcare, Erlangen, Germany) or 3.0-T MR scanner (Ingenia; Philips Healthcare, Best, Netherlands).
Region-of-Interest Segmentation
Tumor delineation was performed using ITK-SNAP software (http://www.itksnap.orge) (25) by two senior radiologists specializing in abdominal diagnoses with 6 and 8 years of work experience. ROIs were manually contoured on each slice in both T2-weighted and diffusion-weighted imaging (DWI) phase images. The two radiologists checked each other’s delineation results for concordance, and a final check was performed by a chief radiologist.
Quantitative Feature Extraction
Radiomics features in ROIs were extracted from the MR images using MATLAB 2016a software (MathWorks, Natick, MA, USA). Image preprocessing was performed to minimize gray range variation between different MR scanners, including voxel size resampling and gray level normalization. For each patient, 968 image features were extracted including 484 T2- and 484 DWI-derived features. The features belonged to four categories: intensity, texture, wavelet, and shape. The feature extraction task was performed as described in our previous work (26, 27).
Immunohistochemistry
Histopathologic examination and diagnosis were carried out according to World Health Organization criteria. Tumor specimens were obtained along with surrounding liver tissue at a volume ratio of 1:1, fixed with formalin, and embedded in paraffin (28). Positive CK19 expression was defined as membranous or cytoplasmic immunoreactivity in ≥5% of tumor cells (29). Pathology data were obtained from the pathology department of each medical center.
Construction of the Cytokeratin 19 Status Classifier
We established three radiomics-based models, i.e., T2, DWI, and combined (both T2 and DWI radiomics features), using an artificial neural network (ANN) algorithm (hidden size = 2, initial random weight = −0.1 to 0.1, weight decay = 5e−4, and maximum number of iterations = 200). We used a two-step process to select optimal radiomics features. We first performed a univariate analysis to identify features showing a significant difference between CK19+ and CK19− groups in the training cohort. The maximum relevance minimum redundancy (MRMR) algorithm (“mRMR” package in R CRAN) was used to rank feature importance. Briefly, input features were ranked by maximizing the mutual information (MI) to class labels and minimizing the MI with other features. To obtain the best prediction model, four different machine learning algorithms, namely, multiple logistic regression (MLR), support vector machine (SVM), random forest (RF), and ANN, were applied to establish radiomics classifiers with optimal features. For each patient, a radiomics score was calculated to determine the probability of positive CK19 status. The diagnostic ability of the classifiers was evaluated based on the area under the receiver operating characteristic (ROC) curve (AUROC). Given the utility of clinical features, we also developed a clinical classifier to determine CK19 status. A model incorporating clinical factors was developed in the training cohort by multivariate logistic regression analysis. Backward stepwise selection was applied using the likelihood ratio test with Akaike’s information criterion (AIC) as the stopping rule.
Ethics Approval
This retrospective study was approved by the institutional review board of all participating centers (FAH-ZJU, SLH, and LSCH). Written informed consent was waived for the retrospective use of patients’ clinical and medical imaging data.
Statistical Analysis
Differences between groups were evaluated with the t-test or Mann–Whitney U test for continuous variables, with the chi-square test for qualitative variables, and with the Wilcoxon test for nonparametric variables. Data were analyzed using SPSS v20 software (IBM, Armonk, NY, USA) and R v3.4.1 software (R Core Team, Vienna, Austria). All statistical tests were two-sided, and a p-value <0.05 was considered statistically significant.
Results
Clinicopathologic Features of the Study Population
The study ultimately enrolled 257 HCC patients including 143 from FAH-ZJU (CK19+:CK19− = 64:79), 75 from SLH (CK19+:CK19− = 34:41), and 39 from LSCH (CK19+:CK19− = 8:31). Baseline clinical and histopathologic characteristics of the participants are shown in Table 1.
CK19 Status Radiomics Classifier Feature Selection and Construction
We identified 64 T2 and 64 DWI radiomics features that differed significantly between the CK19+ and CK19− groups in the training cohort according to the Wilcoxon test. These features were selected with the MRMR algorithm, and the top 15 were retained for further analysis (Supplementary Figure S2). In the T2 model, eight radiomics features were identified as optimal features; in the DWI model, 11 features were selected to construct the classifier; and in the combined model, the top 12 features were used to develop the predictive model. The selected features included skewness_DWI, energy_LLL_T2, uniformity_DWI, denth_T2, SZE_LHH_DWI, maxpr_LLH_T2, SZE_HLH_DWI, SZE_HLH_DWI, and idmnc_T2, Busyness_LLL_DWI, and RP_HLH_DWI.
Of the three radiomics models, the combined model showed the best predictive performance for CK19+ status with AUROCs of 0.857 [95% confidence interval (CI): 0.792–0.922] in the FAH-ZJU training cohort, 0.726 (95% CI: 0.610–0.842) in the SLH cohort, and 0.790 (95% CI: 0.639–0.941) in the LSCH cohort. We also found that the performance of the DWI model [FAH-ZJU: 0.854 (95% CI: 0.782–0.927); SLH: 0.635 (95% CI: 0.507–0.764); LSCH: 0.734 (95% CI: 0.555–0.913)] was superior to that of the T2 model [FAH-ZJU: 0.802 (95% CI: 0.729–0.875); SLH: 0.623 (95% CI: 0.492–0.754); LSCH: 0.605 (95% CI: 0.418–0.792)]. The ROC curves of the three radiomics models are shown in Figure 3.
Figure 3 The comparison of diagnostic performance of the radiomics classifiers and clinical classifier in the training cohort [(A) First Affiliated Hospital of Zhejiang University School of Medicine (FAH-ZJU)] and independent validation cohorts [(B) Shulan Health (Hangzhou) Hospital (SLH); (C) Lishui Central Hospital (LSCH)]. The performance was measured by the area under the receiver operating characteristic (AUROC) curve. The radiomics classifiers consisted of the combined classifier (purple line), the T2 classifier (green line), and the diffusion-weighted imaging (DWI) classifier (orange line). The clinical classifier was shown as the red line.
We also developed classifiers using three other machine learning methods, namely, the MLR, SVM, and RF algorithms. A comparison of the predictive performance of the corresponding models in independent cohorts showed that these three machine learning method-derived classifiers were inferior to the ANN classifier. The ROC curves of the four classifiers are shown in Figure 4.
Figure 4 The comparison of diagnostic performance of artificial neural network (ANN), multiple logistic regression (MLR), random forest (RF), and support vector machine (SVM) classifiers in the training cohort [(A) First Affiliated Hospital of Zhejiang University School of Medicine (FAH-ZJU)] and validation cohorts [(B) Shulan Health (Hangzhou) Hospital (SLH); (C) Lishui Central Hospital (LSCH)]. The performance was measured by the area under the receiver operating characteristic (AUROC) curve. ANN classifier (purple line), MLR classifier (green line), RF classifier (red line), SVM classifier (orange line).
Clinical Classifier Construction
An optimal clinical model was constructed based on the multivariate analysis that used the minimum AIC. Six clinical variables including sex, lymph node status, total bilirubin, direct bilirubin, α-fetoprotein, and log α-fetoprotein were used for model construction. The predictive ability of the clinical model was worse than that of the radiomics-based model in the training and independent validation cohorts [FAH-ZJU: 0.737 (95% CI: 0.654–0.819); SLH: 0.560 (95% CI: 0.425–0.694); LSCH: 0.585 (95% CI: 0.313–0.857)]. The ROC curves for the models are shown in Figure 3.
Discussion
In this multicenter study, we extracted radiomics features from preoperative enhanced MR images and used these to construct a classifier that can accurately distinguish between CK19+ and CK19− HCC patients.
CK19+ HCC exhibits progenitor characteristics and is associated with more malignant behavior, higher lymph node metastasis rate, higher resistance to chemoradiotherapy, and worse outcome after treatment (9–11). Transcriptomic analyses have revealed a positive correlation between HCC subtypes with CK19+ status and poor prognosis including the Hoshida S2, G1, and iClust1 subtypes (30, 31). In in vitro studies, CK19+ HCC cells showed enhanced EMT and induction of angiogenesis, which was abrogated by CK19 knockdown (6, 7, 32). These results indicate that the features of CK19+ HCC fundamentally differ from those of CK19− tumors. Irregular tumor margin, arterial rim enhancement, lower tumor-to-liver signal intensity ratio on hepatobiliary phase MRI, and lower tumor-to-liver apparent diffusion coefficient were found to be significant independent predictors of CK19+ HCC (20). In a multiparametric MRI heterogeneity analysis of HCC, correlations were observed between dynamic contrast-enhanced MRI findings and CK19 status (33). These studies indicate that the unique tumor biology of CK19+ HCC is reflected in MRI-based radiomics features.
Classifiers in previous studies were constructed based on subjective characteristics, which could limit their reliability and generalizability. In contrast, we used time-independent data from different medical centers for the radiomics analysis. We compared the prediction strength of two MRI sequences, namely, DWI and T2, and found that a model combining features from both sequences had the best predictive performance. This is likely because the two models revealed different structural information about the tumors. We identified 128 features differing significantly between patients with opposite CK19 status (i.e., positive and negative). We used four machine learning algorithms to determine the best classifier and found that the one based on the ANN algorithm showed the best prediction performance for CK19 status. The top 12 features were identified as the optimal combination for construction of the ANN classifier, which achieved AUROCs of 0.857 (95% CI: 0.792–0.922) in the FAH-ZJU training cohort and 0.726 (95% CI: 0.610–0.842) and 0.790 (95% CI: 0.639–0.941) in the SLH and LSCH validation cohorts, respectively, outperforming classifiers derived using the other algorithms. The SVM, RF, and MLR algorithms have been used in other previous radiomics analyses. A classifier developed using SVM was able to identify glypican 3-positive HCC (34); one that was based on RF effectively predicted the malignant potential of gastrointestinal stromal tumors (35); and an MLR-based classifier predicted coronavirus disease (COVID)-infected pulmonary lesions (36). We applied the ANN method because it can automatically approximate nonlinear functions with arbitrary accuracy from the training data, outperforming other algorithms by capturing subtle functional relationships in the data (37). Thus, the ANN algorithm enhanced the robustness of our radiomics model.
There were some limitations to our study. Firstly, the sample size was relatively small—especially the LSCH cohort—that could have caused sample bias. Secondly, we used MRI data from three medical centers equipped with different MRI scanners; although the images were processed before they were combined and the features extracted, there may have been bias related to the instruments. Lastly, CK19 is expressed in liver tumors other than HCC (e.g., cholangiocarcinoma, combined cholangiocarcinoma–HCC, and hepatoid adenocarcinoma); whether our classifier can accurately differentiate CK19+ HCC from these tumors remains to be determined.
Conclusion
In conclusion, based on multicenter and time-independent preoperative gadoxetic acid-enhanced MRI data and radiomics analysis, we established a noninvasive and accurate radiomics classifier to determine CK19 status in HCC patients. Our findings can guide the development of personalized management strategies for CK19+ HCC patients that can improve their prognosis.
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 studies involving human participants were reviewed and approved by First Affiliated Hospital Zhejiang University School of Medicine Lishui Central Hospital Shulan Health (Hangzhou) Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
XX and TN designed the study. FY and YCW collected the data. FY, YDW, and LX analyzed the data. FY and YDW wrote the manuscript. XS, JW, and DL did the ROI delineation and recheck. CS and SZ revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Science and Technology Major Project (No. 2017ZX10203205), the National Natural Science Funds for Distinguished Young Scholar of China (No. 81625003), Key Program, National Natural Science Foundation of China (No. 81930016), Key Research & Development Plan of Zhejiang Province (No. 2019C03050), and Natural Science Foundation of China (No. 81801824, No. 81802889, No. 82003248).
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.
The reviewer JQ declared a shared affiliation with several of the authors, FY, LX, YDW, YCW, JW, XS, SZ, and XX, to the handling editor at time of review.
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.
Acknowledgments
We thank Dr. Ruihui Wang and for his help on image check. We also want to express our thanks to Dr. Xue Wen, who worked as our consultant, for her professional pathology advice.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.672126/full#supplementary-material
Supplementary Figure 1 | The patient enrollment of the study.
Supplementary Figure 2 | The mRMR algorithm importance ranking of top 15 radiomics features. The length of the histogram represents the score value of each feature.
References
1. Forner A, Reig M, Bruix J. Hepatocellular Carcinoma. Lancet (2018) 391(10127):1301–14. doi: 10.1016/S0140-6736(18)30010-2
2. Lee JS, Heo J, Libbrecht L, Chu IS, Kaposi-Novak P, Calvisi DF, et al. A Novel Prognostic Subtype of Human Hepatocellular Carcinoma Derived From Hepatic Progenitor Cells. Nat Med (2006) 12(4):410–6. doi: 10.1038/nm1377
3. Tummala KS, Brandt M, Teijeiro A, Grana O, Schwabe RF, Perna C, et al. Hepatocellular Carcinomas Originate Predominantly From Hepatocytes and Benign Lesions From Hepatic Progenitor Cells. Cell Rep (2017) 19(3):584–600. doi: 10.1016/j.celrep.2017.03.059
4. Sia D, Villanueva A, Friedman SL, Llovet JM. Liver Cancer Cell of Origin, Molecular Class, and Effects on Patient Prognosis. Gastroenterology (2017) 152(4):745–61. doi: 10.1053/j.gastro.2016.11.048
5. Matsumoto T, Takai A, Eso Y, Kinoshita K, Manabe T, Seno H, et al. Proliferating EpCAM-Positive Ductal Cells in the Inflamed Liver Give Rise to Hepatocellular Carcinoma. Cancer Res (2017) 77(22):6131–43. doi: 10.1158/0008-5472.CAN-17-1800
6. Takano M, Shimada K, Fujii T, Morita K, Takeda M, Nakajima Y, et al. Keratin 19 as a Key Molecule in Progression of Human Hepatocellular Carcinomas Through Invasion and Angiogenesis. BMC Cancer (2016) 16(1):903. doi: 10.1186/s12885-016-2949-y
7. Govaere O, Komuta M, Berkers J, Spee B, Janssen C, de Luca F, et al. Keratin 19: A Key Role Player in the Invasion of Human Hepatocellular Carcinomas. Gut (2014) 63(4):674–85. doi: 10.1136/gutjnl-2012-304351
8. Kawai T, Yasuchika K, Ishii T, Katayama H, Yoshitoshi EY, Ogiso S, et al. Keratin 19, A Cancer Stem Cell Marker in Human Hepatocellular Carcinoma. Clin Cancer Res (2015) 21(13):3081–91. doi: 10.1158/1078-0432.CCR-14-1936
9. Kim H, Choi GH, Na DC, Ahn EY, Kim GI, Lee JE, et al. Human Hepatocellular Carcinomas With “Stemness”-Related Marker Expression: Keratin 19 Expression and a Poor Prognosis. Hepatology (2011) 54(5):1707–17. doi: 10.1002/hep.24559
10. Lee SH, Lee JS, Na GH, You YK, Kim DG. Immunohistochemical Markers for Hepatocellular Carcinoma Prognosis After Liver Resection and Liver Transplantation. Clin Transplant (2017) 31(1). doi: 10.1111/ctr.12852
11. Fatourou E, Koskinas J, Karandrea D, Palaiologou M, Syminelaki T, Karanikolas M, et al. Keratin 19 Protein Expression Is an Independent Predictor of Survival in Human Hepatocellular Carcinoma. Eur J Gastroenterol Hepatol (2015) 27(9):1094–102. doi: 10.1097/MEG.0000000000000398
12. Guan DX, Shi J, Zhang Y, Zhao JS, Long LY, Chen TW, et al. Sorafenib Enriches Epithelial Cell Adhesion Molecule-Positive Tumor Initiating Cells and Exacerbates a Subtype of Hepatocellular Carcinoma Through TSC2-AKT Cascade. Hepatology (2015) 62(6):1791–803. doi: 10.1002/hep.28117
13. Yuan M, Pu XH, Xu XQ, Zhang YD, Zhong Y, Li H, et al. Lung Adenocarcinoma: Assessment of Epidermal Growth Factor Receptor Mutation Status Based on Extended Models of Diffusion-Weighted Image. J Magn Reson Imaging (2017) 46(1):281–9. doi: 10.1002/jmri.25572
14. Hu XX, Yang ZX, Liang HY, Ding Y, Grimm R, Fu CX, et al. Whole-Tumor MRI Histogram Analyses of Hepatocellular Carcinoma: Correlations With Ki-67 Labeling Index. J Magn Reson Imaging (2017) 46(2):383–92. doi: 10.1002/jmri.25555
15. Moriya T, Saito K, Tajima Y, Harada TL, Araki Y, Sugimoto K, et al. 3D Analysis of Apparent Diffusion Coefficient Histograms in Hepatocellular Carcinoma: Correlation With Histological Grade. Cancer Imaging (2017) 17(1):1. doi: 10.1186/s40644-016-0103-3
16. Choi Y, Kim SH, Youn IK, Kang BJ, Park WC, Lee A. Rim Sign and Histogram Analysis of Apparent Diffusion Coefficient Values on Diffusion-Weighted MRI in Triple-Negative Breast Cancer: Comparison With ER-Positive Subtype. PloS One (2017) 12(5):e0177903. doi: 10.1371/journal.pone.0177903
17. Shen X, Yang F, Yang P, Yang M, Xu L, Zhuo J, et al. A Contrast-Enhanced Computed Tomography Based Radiomics Approach for Preoperative Differentiation of Pancreatic Cystic Neoplasm Subtypes: A Feasibility Study. Front Oncol (2020) 10:248. doi: 10.3389/fonc.2020.00248
18. Joo I, Kim H, Lee JM. Cancer Stem Cells in Primary Liver Cancers: Pathological Concepts and Imaging Findings. Korean J Radiol (2015) 16(1):50–68. doi: 10.3348/kjr.2015.16.1.50
19. Jeong HT, Kim MJ, Kim YE, Park YN, Choi GH, Choi JS, et al. MRI Features of Hepatocellular Carcinoma Expressing Progenitor Cell Markers. Liver Int (2012) 32(3):430–40. doi: 10.1111/j.1478-3231.2011.02744.x
20. Choi SY, Kim SH, Park CK, Min JH, Lee JE, Choi YH, et al. Imaging Features of Gadoxetic Acid-Enhanced and Diffusion-Weighted MR Imaging for Identifying Cytokeratin 19-Positive Hepatocellular Carcinoma: A Retrospective Observational Study. Radiology (2018) 286(3):897–908. doi: 10.1148/radiol.2017162846
21. Hu XX, Wang WT, Yang L, Yang ZX, Liang HY, Ding Y, et al. MR Features Based on LI-RADS Identify Cytokeratin 19 Status of Hepatocellular Carcinomas. Eur J Radiol (2019) 113:7–14. doi: 10.1016/j.ejrad.2019.01.036
22. Wang W, Gu D, Wei J, Ding Y, Yang L, Zhu K, et al. A Radiomics-Based Biomarker for Cytokeratin 19 Status of Hepatocellular Carcinoma With Gadoxetic Acid-Enhanced MRI. Eur Radiol (2020) 30(5):3004–14. doi: 10.1007/s00330-019-06585-y
23. Wang HQ, Yang C, Zeng MS, Rao SX, Ji Y, Weng X, et al. Magnetic Resonance Texture Analysis for the Identification of Cytokeratin 19-Positive Hepatocellular Carcinoma. Eur J Radiol (2019) 117:164–70. doi: 10.1016/j.ejrad.2019.06.016
24. Huang X, Long L, Wei J, Li Y, Xia Y, Zuo P, et al. Radiomics for Diagnosis of Dual-Phenotype Hepatocellular Carcinoma Using Gd-EOB-DTPA-Enhanced MRI and Patient Prognosis. J Cancer Res Clin Oncol (2019) 145(12):2995–3003. doi: 10.1007/s00432-019-03062-3
25. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-Guided 3D Active Contour Segmentation of Anatomical Structures: Significantly Improved Efficiency and Reliability. Neuroimage (2006) 31(3):1116–28. doi: 10.1016/j.neuroimage.2006.01.015
26. 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. doi: 10.3389/fonc.2018.00360
27. Xu L, Yang P, Liang W, Liu W, Wang W, Luo C, et al. A Radiomics Approach Based on Support Vector Machine Using MR Images for Preoperative Lymph Node Status Evaluation in Intrahepatic Cholangiocarcinoma. Theranostics (2019) 9(18):5374–85. doi: 10.7150/thno.34149
28. Cong WM, Bu H, Chen J, Dong H, Zhu YY, Feng LH, et al. Practice Guidelines for the Pathological Diagnosis of Primary Liver Cancer: 2015 Update. World J Gastroenterol (2016) 22(42):9279–87. doi: 10.3748/wjg.v22.i42.9279
29. Durnez A, Verslype C, Nevens F, Fevery J, Aerts R, Pirenne J, et al. The Clinicopathological and Prognostic Relevance of Cytokeratin 7 and 19 Expression in Hepatocellular Carcinoma. A Possible Progenitor Cell Origin. Histopathology (2006) 49(2):138–51. doi: 10.1111/j.1365-2559.2006.02468.x
30. Miltiadous O, Sia D, Hoshida Y, Fiel MI, Harrington AN, Thung SN, et al. Progenitor Cell Markers Predict Outcome of Patients With Hepatocellular Carcinoma Beyond Milan Criteria Undergoing Liver Transplantation. J Hepatol (2015) 63(6):1368–77. doi: 10.1016/j.jhep.2015.07.025
31. Llovet JM, Montal R, Sia D, Finn RS. Molecular Therapies and Precision Medicine for Hepatocellular Carcinoma. Nat Rev Clin Oncol (2018) 15(10):599–616. doi: 10.1038/s41571-018-0073-4
32. Yoneda N, Sato Y, Kitao A, Ikeda H, Sawada-Kitamura S, Miyakoshi M, et al. Epidermal Growth Factor Induces Cytokeratin 19 Expression Accompanied by Increased Growth Abilities in Human Hepatocellular Carcinoma. Lab Invest (2011) 91(2):262–72. doi: 10.1038/labinvest.2010.161
33. Hectors SJ, Wagner M, Bane O, Besa C, Lewis S, Remark R, et al. Quantification of Hepatocellular Carcinoma Heterogeneity With Multiparametric Magnetic Resonance Imaging. Sci Rep (2017) 7(1):2452. doi: 10.1038/s41598-017-02706-z
34. Gu D, Xie Y, Wei J, Li W, Ye Z, Zhu Z, et al. MRI-Based Radiomics Signature: A Potential Biomarker for Identifying Glypican 3-Positive Hepatocellular Carcinoma. J Magn Reson Imaging (2020) 52(6):1679–87. doi: 10.1002/jmri.27199
35. Wang M, Feng Z, Zhou L, Zhang L, Hao X, Zhai J. Computed-Tomography-Based Radiomics Model for Predicting the Malignant Potential of Gastrointestinal Stromal Tumors Preoperatively: A Multi-Classifier and Multicenter Study. Front Oncol (2021) 11:582847. doi: 10.3389/fonc.2021.582847
36. Qiu J, Peng S, Yin J, Wang J, Jiang J, Li Z, et al. A Radiomics Signature to Quantitatively Analyze COVID-19-Infected Pulmonary Lesions. Interdiscip Sci (2021) 13(1):61–72. doi: 10.1007/s12539-020-00410-7
Keywords: cytokeratin 19 (CK19), progenitor attributes, enhanced magnetic resonance imaging, radiomics, hepatocellular carcinoma
Citation: Yang F, Wan Y, Xu L, Wu Y, Shen X, Wang J, Lu D, Shao C, Zheng S, Niu T and Xu X (2021) MRI-Radiomics Prediction for Cytokeratin 19-Positive Hepatocellular Carcinoma: A Multicenter Study. Front. Oncol. 11:672126. doi: 10.3389/fonc.2021.672126
Received: 25 February 2021; Accepted: 20 July 2021;
Published: 12 August 2021.
Edited by:
Xin-Wu Cui, Huazhong University of Science and Technology, ChinaReviewed by:
Jiale Qin, Zhejiang University, ChinaSwapnil Ulhas Rane, Advanced Centre for Treatment, Research and Education in Cancer, India
Copyright © 2021 Yang, Wan, Xu, Wu, Shen, Wang, Lu, Shao, Zheng, Niu and Xu. 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: Xiao Xu, zjxu@zju.edu.cn; Tianye Niu, tyniu@gatech.edu
†These authors have contributed equally to this work