- 1Department of Radiology, Key Laboratory of Clinical Cancer Pharmacology and Toxicology Research of Zhejiang Province, Affiliated Hangzhou First People’s Hospital, Cancer Center, Zhejiang University School of Medicine, Hangzhou, China
- 2Department of Radiology, 903rd Hospital of PLA, Hangzhou, China
- 3Department of Radiology, Liangzhu Hospital, Hangzhou, China
- 4GE Healthcare, Precision Health Institution, Shanghai, China
- 5Zhejiang Chinese Medical University, Hangzhou, China
- 6Department of Radiology, Sir Run Run Shaw Hospital Affiliated to Medical College Zhejiang University, Hangzhou, China
- 7Department of Head and Neck Radiotherapy, Zhejiang Cancer Hospital/Zhejiang Province Key Laboratory of Radiation Oncology, Hangzhou, China
Objective: We aimed to establish an MRI radiomics model and a Delta radiomics model to predict tumor retraction after induction chemotherapy (IC) combined with concurrent chemoradiotherapy (CCRT) for primary nasopharyngeal carcinoma (NPC) in non-endemic areas and to validate its efficacy.
Methods: A total of 272 patients (155 in the training set, 66 in the internal validation set, and 51 in the external validation set) with biopsy pathologically confirmed primary NPC who were screened for pretreatment MRI were retrospectively collected. The NPC tumor was delineated as a region of interest in the two sequenced images of MRI before treatment and after IC, followed by radiomics feature extraction. With the use of maximum relevance minimum redundancy (mRMR) and least absolute shrinkage and selection operator (LASSO) algorithms, logistic regression was performed to establish pretreatment MRI radiomics and pre- and post-IC Delta radiomics models. The optimal Youden’s index was taken; the receiver operating characteristic (ROC) curve, calibration curve, and decision curve were drawn to evaluate the predictive efficacy of different models.
Results: Seven optimal feature subsets were selected from the pretreatment MRI radiomics model, and twelve optimal subsets were selected from the Delta radiomics model. The area under the ROC curve, accuracy, sensitivity, specificity, negative predictive value (NPV), and positive predictive value (PPV) of the MRI radiomics model were 0.865, 0.827, 0.837, 0.813, 0.776, and 0.865, respectively; the corresponding indicators of the Delta radiomics model were 0.941, 0.883, 0.793, 0.968, 0.833, and 0.958, respectively.
Conclusion: The pretreatment MRI radiomics model and pre- and post-IC Delta radiomics models could predict the IC-CCRT response of NPC in non-epidemic areas.
Introduction
Nasopharyngeal carcinoma (NPC) is a malignant tumor originating from the nasopharyngeal mucosal epithelium, which is sensitive to radiotherapy (1). Global Cancer Statistics of 2020 indicated 133,354 new cases of NPC and 80,008 deaths worldwide (2, 3). Meanwhile, the incidence and mortality of NPC in China are higher than the global average estimate. Radiotherapy is the primary treatment for NPC, and concurrent chemoradiotherapy (CCRT) can improve the radiotherapy effect by shrinking the tumor, increasing radiosensitivity, and reducing the radiation dose (4). The 2019 National Comprehensive Cancer Network Clinical Practice Guidelines have recommended induction chemotherapy (IC) combined with CCRT (IC-CCRT) as a class 2A modality for the treatment of advanced NPC (5). In recent years, there has been increasing evidence that IC-CCRT or radiotherapy has clinical value in improving progression-free survival (PFS) and relapse-free survival (RFS) of NPC patients (6). IC-CCRT is also effective in the treatment of non-endemic areas of NPC (7). Although NPC treatment has improved with the advancement of chemoradiotherapy strategies, the 5-year survival rate of some patients with advanced disease is about 60%–85%, and the therapeutic efficacy remains unsatisfactory (8, 9).
The development of treatment options and the evaluation of prognosis for NPC mainly depend on the tumor node metastasis (TNM) stage. However, anatomy-based TNM staging only reflects the tumor shape and invasion into surrounding structures and ignores the internal characteristics of the tumor with the same stage. Hence, despite receiving similar treatment regimens, about 20% of patients show unsatisfactory results due to individual differences and tumor heterogeneity (10). Radiochemotherapy resistance remains one of the main causes of poor prognosis and treatment failure in NPC (11), while residual mass is an independent factor for poor prognosis (12, 13). Treatment of residual disease is associated with better survival outcomes compared to the treatment of recurrent tumors (14). As the tumor shrinks during treatment, adjacent normal brain tissue, skull base bone, and other tissues will fall into the high-dose irradiated tumor areas, increasing the risk of radiation-related injury (15). Therefore, it is necessary to reveal the heterogeneity of tumors as early as possible, facilitating to predict tumor shrinkage in individualized and precise treatment and prognosis of NPC patients.
Radiomics has become a popular method to study tumor heterogeneity in recent years. It can describe tumor heterogeneity and other features by mining the high-dimensional quantitative characteristics of standard medical images (CT, MRI, PET, etc.), providing clinical and high-throughput quantitative information and more personalized treatment options (16). Radiomics characteristics are usually defined in two ways (17), including single-time point radiology and Delta radiomics. Single-time point radiology is mostly used before or during treatment to establish a genomics characteristics model for diagnosis (18, 19), tumor risk stratification (20–23), and prognosis prediction (24, 25), associated with higher powers compared with the TNM staging system. Delta radiomics uses radiological features during or after treatment to provide a wealth of information to identify and quantify treatment-induced changes to guide clinical decisions. It may be more suitable for assessing tumor treatment efficacy (26). Some studies have shown that Delta radiomics-based models yield higher powers than single-time-point-based models (27–29).
Existing radiomics guidelines have recommended the use of multicenter data to ensure the generalizability of the findings (30, 31). However, for many studies, there are very little external validation data. In addition, Delta radiomics studies in predicting adverse events in head and neck squamous cell carcinoma are mostly based on CT imaging (32, 33). The aim of this study was to construct a pretreatment MRI radiomics model and Delta radiomics models before and after IC and to explore their application values in dynamically predicting chemoradiotherapy efficacy for the treatment of NPC in non-epidemic regions.
Materials and Methods
Patients
The study was approved by the local scientific research ethics committee, and informed consent was waived due to its retrospective nature. The study process was in accordance with the Declaration of Helsinki. The information of 668 included patients with pathologically confirmed NPCat the Cancer Hospital of the University of Chinese Academy of Sciences (Zhejiang Cancer Hospital) was collected between January 2007 and June 2012. Then screening was performed according to the following conditions: 1) NPC diagnosis was pathologically confirmed; 2) patients were treated with IC-CCRT; 3) MRI examination was performed within 2 weeks before and after CCRT treatment; 3) fat-suppressed (FS) T2-weighted imaging (T2WI) and FS contrast-enhanced T2-weighted imaging (CE-T1WI) images were available. The patient selection process is shown in Figure 1. Then the data of patients with primary NPC in Sir Run Run Shaw Hospital, Zhejiang University, during the same period were re-screened as external validation data. A total of 51 qualified cases were identified before treatment, including 33 qualified cases after IC.
Figure 1 The case screening process in the training set and internal validation set. A total of 221 qualified cases were screened before the final treatment, including 96 qualified cases after IC. IC, induction chemotherapy.
According to the Response Evaluation Criteria in Solid Tumors (RECIST Version 1.1), NPC patients with response to IC-CCRT were assigned to the complete tumor retraction group. The total retraction group was defined as no evidence of residual disease on reexamination of MRI, that is, complete response (CR). The residual group was defined as residual disease on reexamination of MRI images after IC-CCRT, including partial response (PR), stable disease (SD), or progressive disease (PD).
MRI Acquisition
MRI images of all patients were obtained using two different MRI scanners (Siemens Magnetom Symphony 1.5T and Siemens SKyra 3.0T, Munich, Germany). Axial FS T2WI images were initially obtained, followed by axial FS CE-T1WI imaging after gadolinium-based contrast agent administration at 0.01 mmol/kg. The acquisition protocol for neck MRI was slightly different but mainly consisted of the following parameters: 1) axial FS T2WI, repetition time/echo time (TR/TE) 6,360 ms/95 ms, 90° flip angle, 256 × 168 matrix, slice thickness 4.68 mm, slice spacing 4.68 mm; and 2) axial FS CE-T1WI TR/TE 450 ms/8.8 ms, 90° flip angle, 256 × 168 matrix, slice thickness 4.68 mm, slice spacing 4.68 mm.
MRI Preprocessing and Image Segmentation
The nasopharyngeal MRI examinations before treatment, after IC, and after CCRT of each case were searched from the hospital PACS, and the DICOM format images of FS T2WI and FS CE-T1WI sequences were exported. DICOM format images were imported using ITK-SNAP (www.itksnap.org, version 3.8.0) software. We first selected axial FS CE-T1WI images, in order to improve the accuracy of lesion delineation, and axial FS T2WI images were imported in the “Add Another DICOM Series” option. Two attending radiologists with 7 and 8 years of experience in head and neck radiological diagnosis manually segmented the primary NPC tumor bodies layer by layer. The same approach was used when tumor bodies were delineated on the axial FS T2WI images. Finally, the original images and segmented images were stored according to the format requirements. The same method was applied for tumor segmentation in 96 patients after IC.
Radiomics Feature Extraction and Delta Radiomics Feature Calculation
The study used the software of Artificial Intelligence Kit (AK) V3.4.0.R issued by GE Company (Chicago, IL, USA). This software has been registered, approved, and applied to medical radiomics research (34–36). First, the unsegmented original data were imported into AK software, the original data were resampled, the resolution was adjusted to 1 mm × 1 mm × 1 mm, the layer thickness was 1 mm, the image gray level was uniformly adjusted to 0–255, and then the region of interest (ROI) image was imported. First Order, Shape, GLCM, GLDM, GLSZM, GLRLM, and NGTDM were used for feature selection, and LoG, Wavelet, and LBP were used for filter selection. Radiomics features were extracted from MRI images of 221 patients before treatment and 96 patients after IC.
For Delta radiology profile estimation, the change in each radiology profile was calculated by the following equation:
Here, Feature value2 represents the post-IC MRI value, and Feature value1 represents the pretreatment MRI value.
Feature Selection, Model Establishment, and Statistical Analysis
All statistical analyses were performed using R statistical software (version 4.0.3). The intraclass correlation coefficient (ICC) was used to evaluate inter- and intra-observer agreement. Two attending radiologists randomly selected 40 patients and segmented the tumor once and twice again to calculate observer stability for each feature. Radiological features with ICC greater than 0.75 were defined as stable features. In order to ensure that the images from different devices were comparable, the data were normalized before feature selection, that is, Z-score transformation, and the calculation formula was as follows:
Here, μ is the mean, σ is the standard deviation, and n is the sample size.
The subjects were divided into the training group and the validation cohort in a ratio of 7:3. Then, maximum relevance minimum redundancy (mRMR) and least absolute shrinkage and selection operator (LASSO) were used to select features. We first used mRMR to retain 20 features that were maximally correlated with tags and least redundant with each other and then used LASSO to select the optimal subset of features for model construction. In order to avoid overfitting, we used 10 times of cross-validation to select the adjustment parameter λ. According to the screened characteristics and corresponding coefficients, a logistic regression model of FS T2WI combined with FS CE-T1WI was established, and the radiomics signature (Radscore) was obtained. The area under the receiver operating characteristic (ROC) curve (AUC), diagnostic accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) of the training set and the internal validation set were calculated. ROC analysis was used to evaluate the diagnostic efficacy of the model. Calibration analysis was used to evaluate the goodness-of-fit of the model, and the decision curve was used to evaluate the clinical value of the model. The radiomics flowchart is shown in Figure 2.
Results
Analysis of Radiomics Results
The ICC results were all greater than 0.75, and finally, the segmentation results of the senior doctor were selected. First, we selected 20 features that were most relevant to the tags and least redundant with each other among more than 2,000 features of the two sequences of FS T2WI and FS CE-T1WI by the mRMR method (Figures 3A, B; 4A, B). Then LASSO was performed to select the optimal feature subset for constructing the model. After 10 cross-validations, six and twelve optimal subsets were finally retained (Figures 3C, 4C). The Radscore was then calculated for each patient, and the results were detailed in the Appendix.
Figure 3 Feature extraction before treatment. (A) Uses ten-fold cross-validation and minimum criteria to select the adjustment parameter (λ). (B) LASSO coefficients of 20 texture features, drawing vertical lines at the selected optimal value in λ sequence. (C) Screens seven features. LASSO, least absolute shrinkage and selection operator.
Figure 4 Feature extraction of Delta radiomics. (A) Uses ten-fold cross-validation and minimum criteria to select the adjustment parameter (λ). (B) LASSO coefficients of 20 texture features, drawing vertical lines at the selected optimal value in λ sequence. (C) Screens twelve features. LASSO, least absolute shrinkage and selection operator.
Power Analysis of the Model to Predict the Efficacy of Induction Chemotherapy Combined With Concurrent Chemoradiotherapy
The diagnostic effects of the pretreatment MRI radiomics model and Delta radiomics model in predicting the efficacy of IC-CCRT for NPC are summarized in Table 1. The AUC values (Figure 5), calibration curves (Figure 6), and decision curves (Figure 7) of both models showed good performance in the training set, internal validation set, and external validation set.
Table 1 The diagnostic effects of the pretreatment MRI radiomics model and Delta radiomics model in predicting the efficacy of IC-CCRT for NPC.
Figure 5 ROC characteristic curves and corresponding AUC values for the pretreatment MRI radiomics model (A) and the Delta radiomics model (B). ROC, receiver operating characteristic; AUC, area under the receiver operating characteristic curve.
Figure 6 Calibration curves for MRI radiomics model (A) and Delta radiomics model (B). The thin gray line represents the ideal reference line, the blue line represents the training set, the yellow line represents the internal validation set, and the thick gray line represents the external validation set. In this reference line, regardless of the training set, internal validation set, or external validation set, the predicted probability matches the observed proportion, indicating that both models have good performance in judging the tumor retraction of NPC primary tumors after IC-CCRT treatment. NPC, nasopharyngeal carcinoma; IC-CCRT, induction chemotherapy combined with concurrent chemoradiotherapy.
Figure 7 Decision curves of MRI radiomics model (A) and Delta radiomics model (B). The x-axis represents the probability of response to IC-CCRT, ranging from 0% to 100%; the y-axis measures the net benefit; the blue, yellow, and coarse gray lines represent nomograms of different groups. IC-CCRT, induction chemotherapy combined with concurrent chemoradiotherapy.
Discussion
We selected the seven most relevant radiological features (six from FS T2WI images and one from FS CE-T1WI images) from the pretreatment MRI radiomics model of NPC and twelve most relevant radiological features (seven from FS CE-T1WI images and five from FS T2WI images) from the Delta radiomics model. The clinical predictive values of the two models after chemoradiotherapy for NPC were also analyzed. The results showed that the two constructed models had a high diagnostic performance for NPC, and there were some differences between the two models.
In this study, we selected two MRI sequences, FS T2WI and FS CE-T1WI, for the following three reasons. First, MRI was the most effective and commonly used method for the diagnosis and staging of NPC (37). MRI-based radiomics could better define tumor biology and had immense potential to support oncological decisions (38, 39). Second, these two sequences were commonly obtained during conventional MRI scans and might be more universal. Furthermore, FS sequence and enhanced examination could enhance the image contrast between the tumor and surrounding tissues, and lesion delineation was more accurate. The radiomics characteristics extracted from the combined images were reported to have better predictive performance compared with a single sequence.
In recent years, radiomics has been increasingly used in the diagnosis, tumor stratification, and prognosis prediction of NPC. For example, our previous study showed that the NPC radiomics model based on 18F-FDG PET/MRI images had immense value in the staging evaluation of primary NPC (36). The NPC could be divided into subtypes with different survival patterns based on the radiological characteristics of MRI (40). The early metastasis of NPC could be predicted based on Epstein–Barr virus (EBV), clinical data, and MRI nomogram. Radiomics could predict the PFS of NPC (41), and some radiomics characteristics could identify the patients who needed adaptive radiotherapy (42) and those who would benefit the most from adjuvant therapy or IC (43, 44). A combined model constructed based on EBV, clinical data, and radiological characteristics can predict NPC progression (45) and distant metastasis. However, few of these studies have external validation data, and there have been no reports on Delta radiomics of tumor retraction after IC-CCRT for NPC tumors. Our constructed MRI radiomics model at a single time point before treatment and Delta radiomics before and after IC could predict the efficacy of IC-CCRT for NPC for early prediction during treatment.
The optimal features of the pretreatment MRI radiomics model and Delta radiomics model focused on first-order characteristics, log and texture features after wavelet filtering, and morphological features of the original data. First-order statistics describe the distribution of voxel intensity in the image area defined by the mask through commonly used and basic metrics. About 90th percentile of first-order features and the changes in median and skewness are highly correlated with tumor retraction. The median value represents the gray intensity in the ROI, the asymmetry of the distribution of bias measurement, and the correlation value of the mean value. Sphericity characteristics of morphology in the MRI radiomics model, the roundness of the tumor area relative to the sphere morphology, were significantly correlated with tumor retraction. Our study showed that GLSM, GLCM, and NGTDM were also highly correlated with tumor retraction. GLSZM-ZE represented the uncertainty/randomness of the size of the measurement area and the gray level; the greater the value, the higher the heterogeneity. GLSZM-GLZE represented the distribution of the lower gray level area, and the higher value, the lower the gray value and the proportion of the size area in the image. GLCM describes the texture by studying the spatial correlation of the gray level. Energy represents the uncertainty and randomness of the image. The value of GLCM Joint Entropy represents the complexity of the co-occurrence matrix. In addition to first-order features in the Delta radiomics model, and SALGLE and LALGLE features in GLCM-IMC1, GLSZM is highly correlated with tumor retraction, which can quantify the complexity of texture. SALGLE measures the proportion of smaller sizes with lower gray values in regional joint distribution images. LALGLE measures the proportion of larger sizes with lower gray values in regional joint distribution images. These characteristics were most pronounced after chemoradiation for NPC in this study. Therefore, both the pretreatment MRI radiomics model and Delta radiomics model could predict tumor retraction after IC-CCRT.
However, this study had several limitations. First, this study was a retrospective analysis, and there was inevitable selectivity bias. Second, after IC, the cycle of CCRT was not fixed. We plan to collect more relevant data at a later stage and do more in-depth studies in the future. Third, because the data after treatment of IC were not uniform and the performance of the two models could not be compared, there was a general problem of poor refolding in radiomics. We are going to collect more cases to train the stability of the Delta radiomics model and increase its general applicability.
Conclusion
Pretreatment MRI radiomics at a single time point and Delta radiomics before and after IC could predict tumor retraction after IC-CCRT for the treatment of NPC in non-endemic areas. This study provided a quantitative basis for early intervention and timely optimization of treatment options.
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 the Cancer Hospital of the University of Chinese Academy of Sciences (Zhejiang Cancer Hospital), Affiliated Hangzhou First People’s Hospital and Sir Run Run Shaw Hospital affifiliated to Medical College Zhejiang University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
ZD designed and planned the study. XG, MW, HC, YX, HH, and FJ collected the clinical and image data. YX and HJ performed ROI segmentation. YX, SD, and LW analyzed the data. YX participated in drafting and writing the manuscript. All authors contributed to this study and approved the submitted version.
Funding
This study was supported by grants from the National Natural Science Foundation of China (No. 81871337), Natural Science Foundation of Zhejiang Province (Y22H185692, LGF18H160018, and LY16H180007), Zhejiang Provincial Health Science and Technology Project (No. 2018KY304, No. 2020RC092, and No. 2021RC108), and Hangzhou Medical and Health Science and Technology Project (No. A20200507).
Conflict of Interest
Author SD was employed by GE Healthcare.
The remaining 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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.824509/full#supplementary-material
References
1. De Martel C, Georges D, Bray F, Ferlay J, Clifford GM. Global Burden of Cancer Attributable to Infections in 2018: A Worldwide Incidence Analysis. Lancet Glob Health (2020) 8(2):e180–90. doi: 10.1016/S2214-109X(19)30488-7
2. 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
3. International Agency for Research on Cancer. World Health Organization (2020). Available at: https://gco.iarc.fr/today/home.
4. Kai D, Renjie C. Progress in the Treatment of Nasopharyngeal Carcinoma. J Nanjing Med Univ (Nat Sci Ed) (2021) 41(06):921–926 + 936.
5. Pfister DG, Spencer S, Adelstein D, Adkins D, Anzai Y, Brizel DM, et al. Head and Neck Cancers, Version 2.2020, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw (2020) 18(7):873–98. doi: 10.6004/jnccn.2020.0031
6. Yang Q, Cao SM, Guo L, Hua YJ, Huang PY, Zhang XL, et al. Induction Chemotherapy Followed by Concurrent Chemoradiotherapy Versus Concurrent Chemoradiotherapy Alone in Locoregionally Advanced Nasopharyngeal Carcinoma: Long-Term Results of a Phase III Multicentre Randomised Controlled Trial. Eur J Canc (2019) 119:87–96. doi: 10.1016/j.ejca.2019.07.007
7. Zhao L, Gong J, Xi Y, Xu M, Li C, Kang X, et al. MRI-Based Radiomics Nomogram may Predict the Response to Induction Chemotherapy and Survival in Locally Advanced Nasopharyngeal Carcinoma. Eur Radiol (2020) 30(1):537–46. doi: 10.1007/s00330-019-06211-x
8. Li WF, Chen NY, Zhang N, Hu GQ, Xie FY, Sun Y, et al. Concurrent Chemoradiotherapy With/Without Induction Chemotherapy in Locoregionally Advanced Nasopharyngeal Carcinoma: Long-Term Results of Phase 3 Randomized Controlled Trial. Int J Canc (2019) 145(1):295–305. doi: 10.1002/ijc.32099
9. Sun X, Su S, Chen C, Han F, Zhao C, Xiao W, et al. Long-Term Outcomes of Intensity-Modulated Radiotherapy for 868 Patients With Nasopharyngeal Carcinoma: An Analysis of Survival and Treatment Toxicities. Radiother Oncol (2014) 110(3):398–403. doi: 10.1016/j.radonc.2013.10.020
10. Zhang L, Dong D, Li H, Tian J, Ouyang F, Mo X, et al. Development and Validation of a Magnetic Resonance Imaging-Based Model for the Prediction of Distant Metastasis Before Initial Treatment of Nasopharyngeal Carcinoma: A Retrospective Cohort Study. EBioMedicine (2019) 40:327–35. doi: 10.1016/j.ebiom.2019.01.013
11. Sala E, Mema E, Himoto Y, Veeraraghavan H, Brenton JD, Snyder A, et al. Unravelling Tumour Heterogeneity Using Next-Generation Imaging: Radiomics, Radiogenomics, and Habitat Imaging. Clin Radiol (2017) 72(1):3–10. doi: 10.1016/j.crad.2016.09.013
12. He Y, Zhou Q, Shen L, Zhao Y, Lei M, Wei R, et al. A Retrospective Study of the Prognostic Value of MRI-Derived Residual Tumors at the End of Intensity-Modulated Radiotherapy in 358 Patients With Locally-Advanced Nasopharyngeal Carcinoma. Radiat Oncol (2015) 10:89. doi: 10.1186/s13014-015-0401-0
13. Lv JW, Zhou GQ, Li JX, Tang LL, Mao YP, Lin AH, et al. Magnetic Resonance Imaging-Detected Tumor Residue After Intensity-Modulated Radiation Therapy and its Association With Post-Radiation Plasma Epstein-Barr Virus Deoxyribonucleic Acid in Nasopharyngeal Carcinoma. J Canc (2017) 8(5):861–9. doi: 10.7150/jca.17957
14. Stoker SD, van Diessen JN, de Boer JP, Karakullukcu B, Leemans CR, Tan IB. Current Treatment Options for Local Residual Nasopharyngeal Carcinoma. Curr Treat Opt Oncol (2013) 14(4):475–91. doi: 10.1007/s11864-013-0261-5
15. Yang Y, Yang W, Li FL, Liu MD, LiY. The Effect of Tumor Regression During Radiotherapy on Gross Tumor Volume in Patients With Nasopharyngeal Carcinoma. . J Chin J Otorhinolaryngol Skull Base Surg (2021) 27(03):289–93. doi: 10.11798/j.issn.1007-1520.202103204
16. O’Connor JPB. Cancer Heterogeneity and Imaging. Semin Cell Dev Biol (2017) 64:48–57. doi: 10.1016/j.semcdb.2016.10.001
17. Chang Y, Lafata K, Sun W, Wang C, Chang Z, Kirkpatrick JP, et al. An Investigation of Machine Learning Methods in Delta-Radiomics Feature Analysis. PloS One (2019) 13:14(12):e0226348. doi: 10.1371/journal.pone.0226348
18. Du D, Feng H, Lv W, Ashrafinia S, Yuan Q, Wang Q, et al. Machine Learning Methods for Optimal Radiomics-Based Differentiation Between Recurrence and Inflammation: Application to Nasopharyngeal Carcinoma Post-Therapy PET/CT Images. Mol Imaging Biol (2020) 22(3):730–8. doi: 10.1007/s11307-019-01411-9
19. Zhong X, Li L, Jiang H, Yin J, Lu B, Han W, et al. Cervical Spine Osteoradionecrosis or Bone Metastasis After Radiotherapy for Nasopharyngeal Carcinoma? The MRI-Based Radiomics for Characterization. BMC Med Imag (2020) 20(1):104. doi: 10.1186/s12880-020-00502-2
20. Peng H, Dong D, Fang MJ, Li L, Tang LL, Chen L, et al. Prognostic Value of Deep Learning PET/CT-Based Radiomics: Potential Role for Future Individual Induction Chemotherapy in Advanced Nasopharyngeal Carcinoma. Clin Cancer Res (2019) 25(14):4271–9. doi: 10.1158/1078-0432.CCR-18-3065
21. Kim MJ, Choi Y, Sung YE, Lee YS, Kim YS, Ahn KJ, et al. Early Risk-Assessment of Patients With Nasopharyngeal Carcinoma: The Added Prognostic Value of MR-Based Radiomics. Transl Oncol (2021) 14(10):101180. doi: 10.1016/j.tranon.2021.101180
22. Xu H, Lv W, Feng H, Du D, Yuan Q, Wang Q, et al. Subregional Radiomics Analysis of PET/CT Imaging With Intratumor Partitioning: Application to Prognosis for Nasopharyngeal Carcinoma. Mol Imaging Biol (2020) 22(5):1414–26. doi: 10.1007/s11307-019-01439-x
23. Zhuo EH, Zhang WJ, Li HJ, Zhang GY, Jing BZ, Zhou J, et al. Radiomics on Multi-Modalities MR Sequences can Subtype Patients With non-Metastatic Nasopharyngeal Carcinoma (NPC) Into Distinct Survival Subgroups. Eur Radiol (2019) 29(10):5590–9. doi: 10.1007/s00330-019-06075-1
24. Qiang M, Li C, Sun Y, Sun Y, Ke L, Xie C, et al. A Prognostic Predictive System Based on Deep Learning for Locoregionally Advanced Nasopharyngeal Carcinoma. J Natl Cancer Inst (2021) 113(5):606–15. doi: 10.1093/jnci/djaa149
25. Zhang B, Lian Z, Zhong L, Zhang X, Dong Y, Chen Q, et al. Machine-Learning Based MRI Radiomics Models for Early Detection of Radiation-Induced Brain Injury in Nasopharyngeal Carcinoma. BMC Canc (2020) 20(1):502. doi: 10.1186/s12885-020-06957-4
26. Holscher HC, Bloem JL, van der Woude HJ, Hermans J, Nooy MA, Taminiau AH, et al. Can MRI Predict the Histopathological Response in Patients With Osteosarcoma After the First Cycle of Chemotherapy? Clin Radiol (1995) 50(6):384–90. doi: 10.1016/s0009-9260(05)83135-6
27. Lin P, Yang PF, Chen S, Shao YY, Xu L, Wu Y, et al. A Delta-Radiomics Model for Preoperative Evaluation of Neoadjuvant Chemotherapy Response in High-Grade Osteosarcoma. Cancer Imag (2020) 20(1):7. doi: 10.1186/s40644-019-0283-8
28. Ma Y, Ma W, Xu X, Cao F. How Does the Delta-Radiomics Better Differentiate Pre-Invasive GGNs From Invasive GGNs? Front Oncol (2020) 10:1017. doi: 10.3389/fonc.2020.01017
29. Shayesteh S, Nazari M, Salahshour A, Sandoughdaran S, Hajianfar G, Khateri M, et al. Treatment Response Prediction Using MRI-Based Pre-, Post-, and Delta-Radiomic Features and Machine Learning Algorithms in Colorectal Cancer. Med Phys (2021) 48(7):3691–701. doi: 10.1002/mp.14896
30. Vallières M, Zwanenburg A, Badic B, Cheze Le Rest C, Visvikis D, Hatt M. Responsible Radiomics Research for Faster Clinical Translation. J Nucl Med (2018) 59(2):189–93. doi: 10.2967/jnumed.117.200501
31. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: The Bridge Between Medical Imaging and Personalized Medicine. Nat Rev Clin Oncol (2017) 14(12):749–62. doi: 10.1038/nrclinonc.2017.141
32. Liu Y, Shi H, Huang S, Chen X, Zhou H, Chang H, et al. Early Prediction of Acute Xerostomia During Radiation Therapy for Nasopharyngeal Cancer Based on Delta Radiomics From CT Images. Quant Imaging Med Surg (2019) 9(7):1288–302. doi: 10.21037/qims.2019.07.08
33. Barua S, Elhalawani H, Volpe S, Al Feghali KA, Yang P, Ng SP, et al. Computed Tomography Radiomics Kinetics as Early Imaging Correlates of Osteoradionecrosis in Oropharyngeal Cancer Patients. Front Artif Intell (2021) 4:618469. doi: 10.3389/frai.2021.618469
34. Shao Y, Chen Z, Ming S, Ye Q, Shu Z, Gong C, et al. Predicting the Development of Normal-Appearing White Matter With Radiomics in the Aging Brain: A Longitudinal Clinical Study. Front Aging Neurosci (2018) 10:393. doi: 10.3389/fnagi.2018.00393
35. Feng Q, Song Q, Wang M, Pang P, Liao Z, Jiang H, et al. Hippocampus Radiomic Biomarkers for the Diagnosis of Amnestic Mild Cognitive Impairment: A Machine Learning Method. Front Aging Neurosci (2019) 11:323. doi: 10.3389/fnagi.2019.00323
36. Feng Q, Liang J, Wang L, Niu J, Ge X, Pang P, et al. Radiomics Analysis and Correlation With Metabolic Parameters in Nasopharyngeal Carcinoma Based on PET/MR Imaging. Front Oncol (2020) 10:1619. doi: 10.3389/fonc.2020.01619
37. Dai YL, King AD. State of the Art MRI in Head and Neck Cancer. Clin Radiol (2018) 73(1):45–59. doi: 10.1016/j.crad.2017.05.020
38. Choi ER, Lee HY, Jeong JY, Choi YL, Kim J, Bae J, et al. Quantitative Image Variables Reflect the Intratumoral Pathologic Heterogeneity of Lung Adenocarcinoma. Oncotarget (2016) 7(41):67302–13. doi: 10.18632/oncotarget.11693
39. Cuocolo R, Caruso M, Perillo T, Ugga L, Petretta M. Machine Learning in Oncology: A Clinical Appraisal. Cancer Lett (2020) 481:55–62. doi: 10.1016/j.canlet.2020.03.032
40. Zhang B, Tian J, Dong D, Gu D, Dong Y, Zhang L, et al. Radiomics Features of Multiparametric MRI as Novel Prognostic Factors in Advanced Nasopharyngeal Carcinoma. Clin Cancer Res (2017) 23(15):4259–69. doi: 10.1158/1078-0432.CCR-16-2910
41. Shen H, Wang Y, Liu D, Lv R, Huang Y, Peng C, et al. Predicting Progression-Free Survival Using MRI-Based Radiomics for Patients With Nonmetastatic Nasopharyngeal Carcinoma. Front Oncol (2020) 10:618. doi: 10.3389/fonc.2020.00618
42. Yu TT, Lam SK, To LH, Tse KY, Cheng NY, Fan YN, et al. Pretreatment Prediction of Adaptive Radiation Therapy Eligibility Using MRI-Based Radiomics for Advanced Nasopharyngeal Carcinoma Patients. Front Oncol (2019) 9:1050. doi: 10.3389/fonc.2019.01050
43. Zheng DW, Xu SG, Lai GJ, Chen JY, Ren W, Chen YB. Preliminary Study on Predicting the Efficacy of Neoadjuvant Chemotherapy for Nasopharyngeal Carcinoma Based on Pretreatment Magnetic Resonance Imaging Nomogram. MRI (2021) 12(04):23–9.
44. Zhong L, Dong D, Fang X, Zhang F, Zhang N, Zhang L, et al. A Deep Learning-Based Radiomic Nomogram for Prognosis and Treatment Decision in Advanced Nasopharyngeal Carcinoma: A Multicentre Study. EBioMedicine (2021) 70:103522. doi: 10.1016/j.ebiom.2021.103522
45. Li Q, Wang T, Huang Y, Li Q, Liu P, Grimm R, et al. Whole-Tumor Histogram and Texture Imaging Features on Magnetic Resonance Imaging Combined With Epstein-Barr Virus Status to Predict Disease Progression in Patients With Nasopharyngeal Carcinoma. Front Oncol (2021) 11:610804. doi: 10.3389/fonc.2021.610804
Keywords: nasopharyngeal carcinoma, magnetic resonance imaging, radiomics, induction chemotherapy, concurrent chemoradiotherapy
Citation: Xi Y, Ge X, Ji H, Wang L, Duan S, Chen H, Wang M, Hu H, Jiang F and Ding Z (2022) Prediction of Response to Induction Chemotherapy Plus Concurrent Chemoradiotherapy for Nasopharyngeal Carcinoma Based on MRI Radiomics and Delta Radiomics: A Two-Center Retrospective Study. Front. Oncol. 12:824509. doi: 10.3389/fonc.2022.824509
Received: 30 November 2021; Accepted: 23 March 2022;
Published: 22 April 2022.
Edited by:
Oliver Diaz, University of Barcelona, SpainReviewed by:
Caterina Giannitto, Humanitas Research Hospital, ItalyEros Montin, NYU Grossman School of Medicine, United States
Copyright © 2022 Xi, Ge, Ji, Wang, Duan, Chen, Wang, Hu, Jiang and Ding. 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: Feng Jiang, amlhbmdmZW5nQHpqY2Mub3JnLmNu; Zhongxiang Ding, aGFuZ3pob3Vkeng3M0AxMjYuY29t
†These authors share first authorship