- 1Radiation Oncology Department Sichuan Cancer Hospital and Institute, Affiliated Cancer Hospital of University of Electronic Science and Technology of China, Sichuan Clinical Research Center for Cancer, Chengdu, China
- 2School of Medicine, University of Electronic Science and Technology of China, Chengdu, China
- 3Clinical Medical College, Chengdu University of Traditional Chinese Medicine, Chengdu, China
Magnetic resonance–guided adaptive radiotherapy (MRgART) represents the latest frontier in precision radiotherapy. It is distinguished from other modalities by the possibility of acquiring high-contrast soft tissue images, combined with the ability to recalculate and re-optimize the plan on the daily anatomy. The extensive database of available images offers ample scope for using disciplines such as radiomics to try to correlate features and outcomes. This study aimed to correlate the change of radiomics feature along the treatment to pathological complete response (pCR) for locally advanced rectal cancer (LARC) patients. Twenty-eight LARC patients undergoing neoadjuvant chemoradiotherapy (nCRT) with a short course (25 Gy, 5 Gy × 5f) MRgART at 1.5 Tesla MR-Linac were enrolled. The T2-weighted images acquired at each fraction, corresponding target delineation, pCR result of the surgical specimen, and clinical variables were collected. Seven families of features [First Order, Shape, Gray-level Co-occurrence Matrix (GLCM), Gray-level Dependence Matrix (GLDM), Gray-level Run Length Matrix (GLRLM), Gray-level Size Zone Matrix (GLSZM), and Neighborhood Gray Tone Difference Matrix (NGTDM)] were extracted, and delta features were calculated from the ratio of features of each successive fraction to those of the first fraction. Mann-Whitney U test and LASSO were utilized to reduce the dimension of features and select those features that are most significant to pCR. At last, the radiomics signatures were established by linear regression with the final set of features and their coefficients. A total of 581 radiomics features were extracted, and 2,324 delta features were calculated for each patient. Nineteen features and delta features, and one clinical variable (cN) were significant (p< 0.05) to pCR; seven predictive features were further selected and included in the linear regression to construct the radiomics signature significantly discriminating pCR and non-pCR groups (p< 0.05). Delta features based on MRI images acquired during a short course MRgART could potentially be used to predict treatment response in LARC patients undergoing nCRT.
Introduction
Colon-rectum cancer ranks as the second leading cause of cancer-related mortality worldwide (1, 2), representing approximately 9.4% (935,173) of all new cancer cases. Currently, the standard treatment approach for locally advanced rectal cancer (LARC) involves neoadjuvant radiochemotherapy (nCRT), followed by total mesorectal excision (TME) (3–7). nCRT induces a pathological complete response (pCR) in approximately 11%–42% of LARC patients (8). Extensive research has demonstrated that patients achieving pCR generally experience a more favorable prognosis, characterized by a reduced risk of local failure, lower propensity for local or distant recurrence, improved overall survival, and metastasis-free survival (9, 10). Recent studies have proposed conservative treatment strategies for patients who achieve clinical complete response (cCR) after nCRT, such as local excision or Watch and Wait (W&W), aiming to minimize morbidity and long-term effects associated with unnecessary TME surgery on anorectal and sexual function (11, 12). However, the concept of cCR remains controversial, as it relies on the absence of detectable tumor during clinical and radiological re-evaluation following nCRT. Notably, some research found all patients with pCR also exhibit cCR, not all patients with cCR achieve pCR (13); nevertheless, other studies found different results showing that 61.3% of patients with a pCR had evidence of non-cCR (14). Therefore, early identification of patients who will achieve pCR before TME surgery has become increasingly vital for their accurate multidisciplinary management and remains a topic of ongoing research. Radiomics is an emerging field that includes the extraction of concealed data from medical images and the development of the classification model by machine learning or deep learning techniques. This field aims to provide diagnostic and prognostic insights by detecting features of the tumor or healthy tissue that may not be easily discernible through visual examination (15, 16). Radiomics encompasses the extraction and analysis of spatial layout and geometric shape across various imaging techniques, including computed tomography (CT), magnetic resonance (MR), and positron emission tomography (PET) (17–20). Leveraging these characteristics, numerous studies (21–29) have employed radiomics based on clinical medical images to predict pCR or cCR in LARC, facilitating accurate diagnosis and appropriate treatment selection for patients. Magnetic resonance imaging (MRI) serves as the gold standard in rectal cancer diagnosis, delivering excellent soft tissue contrast and high spatial resolution for assessing LARC response (30). Consequently, several radiomics approaches have been applied to MRI, enabling the prediction of pCR or cCR response to nCRT (22, 31), long-term survival in LARC patients (32), and recognition of different stages of rectal cancer (33). In addition to classical radiomics, delta radiomics is also gaining ground (34); compared to radiomics, which is based on clinical images acquired at a single time point, delta radiomics studies the temporal variation of radiomic features extracted from a dataset of images acquired throughout the treatment. Multiple studies have proposed delta radiomics as an alternative approach for predicting treatment response, tumor prognosis, toxicity, and treatment guidance across various cancer types (35–41). In the specific field of LARC patients, delta radiomics has been proven effective in predicting the response of metastatic renal cell cancer and colorectal cancer liver metastases to chemotherapies (36, 42), cCR or pCR after nCRT (8, 43), disease-free overall survival (44). To the best of our knowledge, the application of delta radiomics based on MR-Linac imaging to predict pCR in LARC patients undergoing a short course of radiotherapy (SCRT) has not been explored previously.
Our study is focused on the treatment of LARC patients treated at 1.5 Tesla Magnetic Resonance Linac (MR-Linac). The use of such a hybrid machine enables a precise magnetic resonance-guided adaptive radiotherapy (MRgART) (45) by the acquisition of diagnostic-quality MRI scans before each daily fraction, facilitating accurate dose delivery, hypo-fractionated dose escalation to the tumor, and improving local control rates with reduced toxicity rates. Importantly, via the use of an MR-Linac, we can collect multiple MRI images throughout the treatment process and, therefore, there is ample room for delta radiomics analysis.
This study aims to predict the response of LARC patients undergoing a SCRT with a 1.5 Tesla MR-Linac. Specifically, we aim to assess the feasibility of using delta radiomics as a predictive tool to determine pCR of LARC patients after nCRT. By investigating the potential of delta radiomics, this research strives to enhance patient management and treatment outcomes.
Materials and methods
Patients
This retrospective study includes a cohort of 28 patients diagnosed with pathologically confirmed LARC. Specifically, patients with clinical stage T3-4N0 or T3-4N1-2 who received MRgART using a 1.5 Tesla MR-Linac were included in this retrospective study. To maintain the study’s focus, exclusion criteria comprised patients with distant metastases, a history of prior chemotherapy or radiotherapy for rectal cancer, concurrent or previous malignancies, and individuals with known allergies to intravenous contrast agents or contraindications to MRI. This study was approved by the Institutional Ethics Committee of our hospital.
Patient clinical workflow
Patients underwent a SCRT followed after one week by chemotherapy, then evaluated as cCR or non-cCR according to the criteria described in the next section. In the absence of cCR, surgery was scheduled, otherwise, patients followed the W&W strategy with close follow-up consisting of regular restaging by MRI and digital rectal examination.
SCRT was administered over 1 week delivering 25 Gy in 5 Gy per fraction. Each treatment fraction starts with the acquisition of an online MRI used for the plan adaptation. The pre-treatment CT reference image, contours, and plan, together with the daily online MRI are used as input to adapt the plan for that specific session. The reference image is matched with the online MRI through rigid registration, then the workflow proceeds with an adapt to shape (ATS) approach based on the new patient anatomy; target and organs at risk were adapted to the daily anatomy using deformable registration with normalized mutual information of Monaco (Elekta AB, Stockholm, Sweden) v.5.4, and then manually adjusted (if needed) by the radiation oncologist. The CT reference treatment plan includes all the information needed to generate the synthetic CT (sCT) including average electron density and a given layering priority for each contour; the adapted plan is then recalculated on the daily MR-based sCT.
Chemotherapy was administered for 12 weeks, comprising six cycles of mFOLFOX6 regimen. This regimen involved intravenous administration of Oxaliplatin (85 mg/m2) and folinic acid (400 mg/m2) on day 1, intravenous push of 5-fluorouracil (400 mg/m2) on day 1, and continuous intravenous pumping of 5-fluorouracil (2400 mg/m2) over 46h–48h, repeated every 14 days.
TME surgery was performed through either anterior resection or abdominoperineal resection, depending on the individual patient’s requirements and circumstances.
cCR and pCR evaluation criteria
In this study, a patient was considered as cCR when the following criteria were fulfilled: absence of lymph nodes on restaging MRI imaging, no residual primary tumor evident on morphological and diffusion-weighted imaging with intact rectal wall layers, absence of detectable mass during manual rectal examination, and no residual lesions or presence of a flat scar during endoscopic examination.
The pathologic staging was determined following the tumor node metastasis (TNM) classification system, recommended by the American Joint Committee on Cancer (AJCC), 8th edition (46). Histopathological examination of surgically resected specimens was performed by two pathologists with 10 years of experience in rectal tumor pathology who were blinded to the MRI data. Tumor response to nCRT was classified using the tumor regression grade (TRG) classification system proposed by Mandard et al. (47). Specifically, patients were categorized as pCR if TRG = 1, or non-pCR if TRG > 1.
Imaging dataset and feature extraction
Patient’s MR images acquired with the 1.5 Tesla MR scanner at Unity MR-Linac during the online adaptive clinical workflow were used for features extraction. The routine scanning protocol consists of 3D-T2-weighted transversal sequences, with image resolution of 0.83 mm × 0.83 mm, slice thickness of 1.5 mm, and acquisition time of 117 s. For each patient, a total of five MR scans, corresponding to fractions 1 to 5, were included in the analysis. Each MRI image was normalized using Z-score normalization method. In addition, to address the issue of high-frequency MRI signal noise and mitigate the impact of significant signal variation, the Laplacian of Gaussian (LoG) filters were applied to raw MRI images. The LoG filters utilized a range of standard deviations (σ) from 0.2 to 1.0 mm, with a step size of 0.2 mm. The calculation is as follows:
The segmentation of the gross tumor volume (GTV) used for features extraction is obtained at the daily plan adaptation using first the deformable registration module of Monaco treatment planning system and then a slice-by-slice manual adjustment by the radiation oncologist. The radiomics analysis was focused on the whole GTV volume for each of the five fractions. GTV delineation for a representative patient along the treatment course is shown in Figure 1.
Figure 1 Gross tumor volume delineation (in red) on sagittal planes for the five ((A–E), respectively) MRgART fractions of a representative SCRT.
Seven families of radiomics features, including First Order, Shape, GLCM, GLDM, GLRLM, GLSZM, and NGTDM, were extracted from both the raw images and the LoG-filtered images. In order to capture changes in the features throughout the treatment, the delta features (ΔF) were calculated as the ratio of the features acquired at fractions 2–5 (F2–F5) to the corresponding features at fraction 1 (F1) as follows, resulting in four datasets (ΔF2, ΔF3, ΔF4, and ΔF5). Additionally, three potential clinical variables, namely clinical T-stage (cT), clinical N-stage (cN), and age, were also collected for subsequent analysis.
Where F1 represents the feature extracted from MRI of the first fraction, F1 is the features from MRI of ith fraction.
Feature selection and radiomics signature
Firstly, the Mann-Whitney U test was employed to identify potential predictive features (p< 0.05) capable of discriminating between patients showing pCR and non-pCR. The test was applied to features, delta features extracted and on the three clinical variables (cN, cT, and age of the patient).
To further reduce the dimensionality of the feature set, the least absolute shrinkage and selection operator (LASSO) was utilized to select the most significant features; the regularization process shrinks the coefficients associated with each variable toward zero, thus retaining only the most relevant features while discarding the less important ones. This was achieved by increasing the regularization factor λ and applying the leave-one-out cross-validation (LOOCV) to determine the optimal set of features. Ultimately, the features and their corresponding coefficients that yielded the lowest “Binomial Deviance” in LASSO were used to construct the radiomics signature (Rad-score) as follows:
Where n is the number of features, Xi and Ci are the ith feature and their corresponding coefficients, and b is the intercept.
Statistical analysis
The MRI images were preprocessed in a Python environment, and features extraction was performed using the PyRadiomics V3.0.1 Python package. Subsequently, the Mann-Whitney U test of the image features, LASSO regression, and evaluation were conducted using R statistical software (version 3.6.1, R Foundation for Statistical Computing, Vienna, Austria).
Results
Patients and treatment characteristics
The median age of the 28 patients included in this study was 52.5 years with a higher proportion of males than females (78.6% vs. 21.4%, respectively). Among the patients enrolled, 26 (92.9%) were identified as non-cCR patients and underwent TME surgery at the end of chemotherapy. Of the two cCR patients (7.1%), one remained cCR for 22 weeks after the end of nCRT and then tumor recurrence was detected and underwent surgery; the second one, 45 weeks after the end of nCRT was still cCR at the time of the study, so no pCR information was available. By the end, 27 patients had already undergone surgery and histopathological examination revealed that five of them (18.5%) were pCR. For this radiomics study, 135 MRI datasets and corresponding GTV segmentations of 27 of the 28 LARC patients treated with MRgART at Unity MR-Linac for a SCRT from August 2021 to date were analyzed. The clinical characteristics and response to treatment of the patients are summarized in Table 1.
Features selection
A total of 581 features were extracted from each MRI image, including 111 features from the raw image and 470 features from the LoG-filtered images. As for delta features, a total of 2,324 delta features were calculated from the 581 features of each of the four fractions following the first.
The Mann-Whitney U test revealed that 18 variables had a significant association with pCR (p< 0.05). These variables included six features, one clinical variable (cN), and 11 delta features (six delta features at second fraction, four delta features at third fraction, one delta feature at fourth fraction). No significant delta features were found at fifth fraction. A standard deviation of 0.2 mm or 0.4 mm of the LoG-filtered images contributed to the extraction of significant features; the results obtained are shown in Table 2.
Using LASSO regression, the dimension of the 18 variables was further reduced, thus retaining only the most relevant features and discarding the less important ones. Figure 2 illustrates the reduction in feature dimensionality for increasing values of Lambda. The binomial deviance of the LASSO regression varies with the regularization parameter lambda. The optimal lambda value (Λ.min) of 0.0245 was obtained when the LASSO regression achieved the lowest binomial deviance using LOOCV validation as shown in Figure 3.
Figure 3 Optimal lambda value determination by minimization of the binomial deviance. The optimal lambda value and corresponding dimension of predictive features with non-zero coefficients are highlighted by the vertical dotted line.
With this optimal lambda value, 11 features had a zero coefficient in the LASSO regression and were eliminated; of the seven remaining ones, we obtained one clinical variable (cN), one feature from first fraction, four delta features at second fraction, and one delta feature at third fraction, respectively. The final set of significant features and the corresponding coefficients to construct the rad score for patients are listed in Table 3. The distributions of all final predictive features characterising the pCR and non-pCR patient groups are shown in Figure 4.
Table 3 The optimal set of features and corresponding coefficients used to construct the rad score for patients.
Figure 4 The distributions of the final predictive features and delta features characterizing the pCR and non-pCR patient group. The lower and upper borders of the box indicate the 25th and 75th percentile of the data, respectively, while the whiskers indicate the fifth and 95th percentile of the data.
Radiomics signature
With the seven variables and their corresponding coefficients, the Rad-Score for each patient was constructed by LASSO regression. Figure 5 shows that the Rad-Score can significantly discriminate pCR and non-pCR patients (p< 0.05).
Figure 5 Capability of radiomics signatures (rad score) to significantly distinguish pCR and non-pCR patients.
Discussion
Artificial Intelligence has gained significant attention in recent years for its potential contribution to achieving full-personalized medicine (48); in the treatment of cancer patients, its role is fundamental throughout the whole process, not only to automatize and speed up the segmentation process, to set up automatic prediction models for the gamma passing rate of the patient’s treatment plan (49, 50), but also in the massive use and analysis of diagnostic imaging to elaborate predictive models for patient’s outcomes. In the management of LARC, personalized medicine is essential; the W&W approach is considered a favorable option for patients assessed as cCR after nCRT to avoid morbidity and toxicity associated with unnecessary surgery. However, the definitive determination of residual tumors relies on the evaluation of resected specimens. Therefore, the ability to predict pCR without the need for an actual total TME procedure would greatly facilitate the personalized treatment approach. Patients predicted to achieve pCR based on reliable models would have a stronger rationale to choose the W&W strategy, whereas those anticipated as non-pCR cases may contemplate additional medical interventions, such as consolidation chemotherapy and/or radiotherapy boost.
Delta radiomics is developed upon conventional radiomics analysis while taking the variation of features throughout the treatment process into consideration. Since radiomics features before or after treatment are not sufficient to describe all characteristics of tumor, the delta radiomics approach might be potentially more powerful for prognostic purposes. MRgART workflow at MR-Linac provides favorable resources for delta radiomics study because multiple MRI images are available for each radiotherapy fraction.
In this study, the delta radiomics approach using the Mann-Whitney U test, Pearson coefficient analysis and LASSO regression was applied to construct the radiomics signature, to identify pCR in a cohort of LARC patients treated in a SCRT at 1.5 Tesla MR. Finally, one original feature and five delta features are selected as the final predictors. Our finding suggests that delta features have the potential to be powerful predictors of pCR for LARC. Other studies also found the potential of delta radiomics when applied to various clinical end points for LARC patients. Boldrini et al. also found that delta features have better discrimination ability than standard radiomics features, even if their study is conducted on the prediction of cCR for LARC. Chiloiro et al. (44) found that the change of area/surface ratio is most significant in identifying 2-year-DFS for LARC. Shayesteh et al. (43) also proved that their delta radiomics model outperformed both pre- and post-treatment features in pCR prediction.
Repeatability and generalization are major concerns in radiomics studies, particularly for delta radiomics studies; heterogeneity in the acquisition of images and timing are important factors to consider (8, 39, 41, 44, 51–54). Those facts might limit the transferability of the results. In our study, delta radiomics analysis was based on MRI from an MR-Linac, allowing reproducibility of the MRI sequences used during the treatment itself and a timing that is guided by the treatment delivery days, and therefore easily reproducible within the same center, and hopefully for multicenter studies (8, 37).
The study is subject to limitations. The number of patients enrolled is too small to set up a fine-tuned prediction model; there is difficulty in performing internal and external validations, reflecting the limited number of LARC patients treated with the relatively new MRgART technique. Moreover, considering that, for each patient, the MRI images were acquired during the SCRT over 1 week, the change in tumor heterogeneity and texture might not be fully exhibited; thus, the change of features might be still conservative. Despite these limitations, the study is the first to investigate and highlight the feasibility of using delta radiomics with promising prediction results on MR-Linac in SCRT schedule for LARC. For a more comprehensive delta radiomics study, future works are under design to include in the analysis one or more MRIs at the end of the SCRT, the continuous updating of the case history, and including internal and external validations as soon as numbers allow.
Conclusion
Delta features approach based on MRI images acquired during MRgART at MR-Linac could potentially be used for treatment response prediction in LARC patients undergoing nCRT with SCRT.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
This research was approved by Ethics Committee of Sichuan Cancer Hospital with the approval number SCCHEC-02-2022-164.
Author contributions
BT, LO, and QP contributed to conception and design of the study. JX and LO organized the database. YL, FW, and CL performed the statistical analysis. JW and JX wrote the first draft of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by Deng Feng Precision Radiotherapy Research Program (2021DF016) and Sichuan Cancer Hospital Youth Project (YB2021032, YB2023025).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
2. Siegel RL, Wagle NS, Cercek A, Smith RA, Jemal A. Colorectal cancer statistics, 2023. CA Cancer J Clin (2023) 73(3):233–54. doi: 10.3322/caac.21772
3. Bosset JF, Collette L, Calais G, Mineur L, Maingon P, Radosevic-Jelic L, et al. Chemotherapy with preoperative radiotherapy in rectal cancer. N Engl J Med (2006) 355(11):1114–23. doi: 10.1056/NEJMoa060829
4. Kapiteijn E, Marijnen C, Nagtegaal ID, Putter H, Velde CJVD. Preoperative radiotherapy combined with total mesorectal excision for resectable rectal cancer. New Engl J Med (2001) 345(9):638–46. doi: 10.1056/NEJMoa010580
5. Sauer R, Liersch T, Merkel S, Fietkau R, Hohenberger W, Hess C, et al. Preoperative versus postoperative chemoradiotherapy for rectal cancer. J Clin Oncol (2004) 351(16):1926–33. doi: 10.1056/NEJMoa040694
6. Kang MK. Implications of recent neoadjuvant clinical trials on the future practice of radiotherapy in locally advanced rectal cancer. World J Gastroenterol (2023) 29(6):1011. doi: 10.3748/wjg.v29.i6.1011
7. Yuan Z, Weng S, Ye C, Hu H, Zhang S, Yuan Y. CSCO guidelines for colorectal cancer version 2022: updates and discussions. Chin J Cancer Res (2022) 34(2):67. doi: 10.21147/j.issn.1000-9604.2022.02.01
8. Boldrini L, Cusumano D, Chiloiro G, Casa C, Masciocchi C, Lenkowicz J, et al. Delta radiomics for rectal cancer response prediction with hybrid 0.35 T magnetic resonance-guided radiotherapy (MRgRT): a hypothesis-generating study for an innovative personalized medicine approach. Radiol Med (2019) 124(2):145–53. doi: 10.1007/s11547-018-0951-y
9. Macchia G, Gambacorta MA, Masciocchi C, Chiloiro G, Mantello G, di Benedetto M, et al. Time to surgery and pathologic complete response after neoadjuvant chemoradiation in rectal cancer: a population study on 2094 patients. Clin Trans Radiat Oncol (2017) 4:8–14. doi: 10.1016/j.ctro.2017.04.004
10. Vecchio FM, Valentini V, Minsky BD, Padula GDA, Venkatraman ES, Balducci M, et al. The relationship of pathologic tumor regression grade (TRG) and outcomes after preoperative therapy in rectal cancer. Int J Radiat Oncol Biol Phys (2005) 62(3):752–60. doi: 10.1016/j.ijrobp.2004.11.017
11. Renehan AG, Malcomson L, Emsley R, Gollins S, O’Dwyer ST. Watch-and-wait approach versus surgical resection after chemoradiotherapy for patients with rectal cancer (the OnCoRe project): a propensity-score matched cohort analysis. Lancet Oncol (2016) 17(2):174–83. doi: 10.1016/S1470-2045(15)00467-2
12. Beets-Tan R, Leijtens J, Beets GL. Wait-and-see policy for clinical complete responders after chemoradiation for rectal cancer. J Clin Oncol (2011) 29:4633–40. doi: 10.1200/JCO.2011.37.7176
13. Hiotis SP, Weber SM, Cohen AM, Minsky BD, Paty PB, Guillem JG, et al. Assessing the predictive value of clinical complete response to neoadjuvant therapy for rectal cancer: an analysis of 488 patients. J Am Coll Surgeons (2002) 194(2):131–5. doi: 10.1016/S1072-7515(01)01159-0
14. Smith F, Chang K, Sheahan K, Hyland J, O’Connell P, Winter D. The surgical significance of residual mucosal abnormalities in rectal cancer following neoadjuvant chemoradiotherapy. J Br Surg (2012) 99(7):993–1001. doi: 10.1002/bjs.8700
15. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology (2016) 278(2):563–77. doi: 10.1148/radiol.2015151169
16. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, Aerts HJWL. Radiomics: Extracting more information from medical images using advanced feature analysis. Eur J Cancer (2012) 43(4):441–6. doi: 10.1016/j.ejca.2011.11.036
17. Giannini V, Mazzetti S, Bertotto I, Chiarenza C, Cauda S, Delmastro E, et al. Predicting locally advanced rectal cancer response to neoadjuvant therapy with18F-FDG PET and MRI radiomics features. Eur J Nucl Med Mol Imaging (2019) 46(4):878–88. doi: 10.1007/s00259-018-4250-6
18. Cook GJR, Siddique M, Taylor BP, Yip C, Chicklore S, Goh V. Radiomics in PET: principles and applications. Clin Trans Imaging (2014) 2(3):269–76. doi: 10.1007/s40336-014-0064-0
19. Mcgarry SD, Hurrell SL, Kaczmarowski AL, Cochran EJ, Laviolette PS. Magnetic resonance imaging-based radiomic profiles predict patient prognosis in newly diagnosed glioblastoma before therapy. Tomography (2016) 2(3):223–8. doi: 10.18383/j.tom.2016.00250
20. Coroller TP, Grossmann P, Hou Y, Velazquez ER, Leijenaar RT, Hermann G, et al. CT-based radiomic signature predicts distant metastasis in lung adenocarcinoma. Radiotherapy Oncol (2015) 114(3):345–50. doi: 10.1016/j.radonc.2015.02.015
21. Yi X, Pei Q, Zhang Y, Zhu H, Wang Z, Chen C, et al. MRI-based radiomics predicts tumor response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer. Front Oncol (2019) 9:552. doi: 10.3389/fonc.2019.00552
22. Zhou X, Yi Y, Liu Z, Cao W, Lai B, Sun K, et al. Radiomics-based pretherapeutic prediction of non-response to neoadjuvant therapy in locally advanced rectal cancer. Ann Surg Oncol (2019) 26(6):1676–84. doi: 10.1245/s10434-019-07300-3
23. Bulens P, Couwenberg A, Intven M, Debucquoy A, Vandecaveye V, Van Cutsem E, et al. Predicting the tumor response to chemoradiotherapy for rectal cancer: Model development and external validation using MRI radiomics. Radiother Oncol (2020) 142:246–52. doi: 10.1016/j.radonc.2019.07.033
24. Cui Y, Yang W, Ren J, Li D, Du X, Zhang J, et al. Prognostic value of multiparametric MRI-based radiomics model: Potential role for chemotherapeutic benefits in locally advanced rectal cancer. Radiother Oncol (2020) 154:161–9. doi: 10.1016/j.radonc.2020.09.039
25. Cusumano D, Meijer G, Lenkowicz J, Chiloiro G, Boldrini L, Masciocchi C, et al. A field strength independent MR radiomics model to predict pathological complete response in locally advanced rectal cancer. Radiol Med (2020) 126(3):421–9. doi: 10.1007/s11547-020-01266-z
26. Petresc B, Lebovici A, Caraiani C, Feier DS, Graur F, Buruian MM. Pre-treatment T2-WI based radiomics features for prediction of locally advanced rectal cancer non-response to neoadjuvant chemoradiotherapy: A preliminary study. Cancers (Basel) (2020) 12(7):1894. doi: 10.3390/cancers12071894
27. Shayesteh SP, Alikhassi A, Farhan F, Ghalehtaki R, Soltanabadi M, Haddad P, et al. Prediction of response to neoadjuvant chemoradiotherapy by MRI-based machine learning texture analysis in rectal cancer patients. J Gastrointest Cancer (2020) 51(2):601–9. doi: 10.1007/s12029-019-00291-0
28. Delli Pizzi A, Chiarelli AM, Chiacchiaretta P, d’Annibale M, Croce P, Rosa C, et al. MRI-based clinical-radiomics model predicts tumor response before treatment in locally advanced rectal cancer. Sci Rep (2021) 11(1):5379. doi: 10.1038/s41598-021-84816-3
29. Boldrini L, Lenkowicz J, Orlandini LC, Yin G, Cusumano D, Chiloiro G, et al. Applicability of a pathological complete response magnetic resonance-based radiomics model for locally advanced rectal cancer in intercontinental cohort. Radiat Oncol (2022) 17(1):78. doi: 10.1186/s13014-022-02048-9
30. Barbaro B, Fiorucci C, Tebala C, Valentini V, Gambacorta MA, Vecchio FM, et al. Locally advanced rectal cancer: MR imaging in prediction of response after preoperative chemotherapy and radiation therapy. Radiology (2009) 250(3):730–9. doi: 10.1148/radiol.2503080310
31. Cusumano D, Dinapoli N, Boldrini L, Chiloiro G, Gatta R, Masciocchi C, et al. Fractal-based radiomic approach to predict complete pathological response after chemo-radiotherapy in rectal cancer. Radiol Med (2018) 123(4):286–95. doi: 10.1007/s11547-017-0838-3
32. Jalil O, Afaq A, Ganeshan B, Patel UB, Boone D, Endozo R, et al. Magnetic resonance based texture parameters as potential imaging biomarkers for predicting long-term survival in locally advanced rectal cancer treated by chemoradiotherapy. Colorectal Dis (2017) 19(4):349–62. doi: 10.1111/codi.13496
33. Liu L, Liu Y, Xu L, Li Z, Lv H, Dong N, et al. Application of texture analysis based on apparent diffusion coefficient maps in discriminating different stages of rectal cancer. J Magnetic Resonance Imaging (2017) 45(6):1798–808. doi: 10.1002/jmri.25460
34. Nardone V, Reginelli A, Grassi R, Boldrini L, Vacca G, D’Ippolito E, et al. Delta radiomics: a systematic review. Radiol Med (2021) 126(12):1571–83. doi: 10.1007/s11547-021-01436-7
35. Cunliffe A, Armato SG, Castillo R, Pham N, Guerrero T, Al-Hallaq HA. Lung texture in serial thoracic computed tomography scans: correlation of radiomics-based features with radiation therapy dose and radiation pneumonitis development. Int J Radiat Oncology Biology Phys (2015) 91(5):1048–56. doi: 10.1016/j.ijrobp.2014.11.030
36. Rao S-X, Lambregts DM, Schnerr RS, Beckers RC, Maas M, Albarello F, et al. CT texture analysis in colorectal liver metastases: a better way than size and volume measurements to assess response to chemotherapy? United Eur Gastroenterol J (2016) 4(2):257–63. doi: 10.1177/2050640615601603
37. Cusumano D, Boldrini L, Yadav P, Casà C, Valentini V. Delta radiomics analysis for local control prediction in pancreatic cancer patients treated using magnetic resonance guided radiotherapy. Diagnostics (Basel Switzerland) (2021) 11(1):72. doi: 10.3390/diagnostics11010072
38. Jeon SH, Song C, Chie EK, Kim B, Kim YH, Chang W, et al. Delta-radiomics signature predicts treatment outcomes after preoperative chemoradiotherapy and surgery in rectal cancer. Radiat Oncol (2019) 14(1):43. doi: 10.1186/s13014-019-1246-8
39. Fave X, Zhang L, Yang J, Mackin D, Balter P, Gomez D, et al. Delta-radiomics features for the prediction of patient outcomes in non-small cell lung cancer. Sci Rep (2017) 7(1):588. doi: 10.1038/s41598-017-00665-z
40. 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. Quantitative Imaging Med Surg (2019) 9(7):1288–302. doi: 10.21037/qims.2019.07.08
41. Basler L, Gabrys HS, Hogan SA, Pavic M, Bogowicz M, Vuong D, et al. Radiomics, tumor volume, and blood biomarkers for early prediction of pseudoprogression in patients with metastatic melanoma treated with immune checkpoint inhibition. Clin Cancer Res (2020) 26(16):4414–25. doi: 10.1158/1078-0432.CCR-20-0020
42. Vicky G, Balaji G, Paul N, Juttla JK, Anup V, Miles KA. Assessment of response to tyrosine kinase inhibitors in metastatic renal cell cancer: CT texture as a predictive biomarker. Int J Med Radiol (2011) 261(1):165–71.
43. 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
44. Chiloiro G, Boldrini L, Preziosi F, Cusumano D, Yadav P, Romano A, et al. A predictive model of 2yDFS during MR-guided RT neoadjuvant chemoradiotherapy in locally advanced rectal cancer patients. Front Oncol (2022) 12:831712. doi: 10.3389/fonc.2022.831712
45. Kerkmeijer LGW, Fuller CD, Verkooijen HM, Marcel V, Ananya C, Harrington KJ, et al. The MRI-linear accelerator consortium: evidence-based clinical introduction of an innovation in radiation oncology connecting researchers, methodology, data collection, quality assurance, and technical development. Front Oncol (2016) 6:215. doi: 10.3389/fonc.2016.00215
46. Weiser MR. AJCC 8th edition: colorectal cancer. Ann Surg Oncol (2018) 25:1454–5. doi: 10.1245/s10434-018-6462-1
47. Mandard AM, Dalibard F, Mandard JC, Marnay J, Henry-Amar M, Petiot JF, et al. Pathologic assessment of tumor regression after preoperative chemoradiotherapy of esophageal carcinoma. Clinicopathologic correlations. Cancer (1994) 73(11):2680–6. doi: 10.1002/1097-0142(19940601)73:11<2680::AID-CNCR2820731105>3.0.CO;2-C
48. Cusumano D, Boldrini L, Dhont J, Fiorino C, Indovina L. Artificial Intelligence in magnetic Resonance guided Radiotherapy: Medical and physical considerations on state of art and future perspectives. Physica Med (2021) 85:175–91. doi: 10.1016/j.ejmp.2021.05.010
49. Shiba E, Saito A, Furumi M, Kawahara D, Miki K, Murakami Y, et al. Predictive gamma passing rate for three-dimensional dose verification with finite detector elements via improved dose uncertainty potential accumulation model. Med Phys (2020) 47(3):1349–56. doi: 10.1002/mp.13985
50. Lam D, Zhang X, Li H, Deshan Y, Schott B, Zhao T, et al. Predicting gamma passing rates for portal dosimetry-based IMRT QA using machine learning. Med Phys (2019) 46(10):4666–75. doi: 10.1002/mp.13752
51. Fan M, Chen H, You C, Liu L, Gu Y, Peng W, et al. Radiomics of tumor heterogeneity in longitudinal dynamic contrast-enhanced magnetic resonance imaging for predicting response to neoadjuvant chemotherapy in breast cancer. Front Mol Biosci (2021) 8:622219. doi: 10.3389/fmolb.2021.622219
52. Wang L, Gao Z, Li C, Sun L, Li J, Yu J, et al. Computed tomography-based delta-radiomics analysis for discriminating radiation pneumonitis in patients with esophageal cancer after radiation therapy. Int J Radiat Oncology Biology Phys (2021) 111(2):443–55. doi: 10.1016/j.ijrobp.2021.04.047
53. Crimi F, Capelli G, Spolverato G, Bao QR, Florio A, Milite Rossi S, et al. MRI T2-weighted sequences-based texture analysis (TA) as a predictor of response to neoadjuvant chemo-radiotherapy (nCRT) in patients with locally advanced rectal cancer (LARC). Radiol Med (2020) 125(12):1216–24. doi: 10.1007/s11547-020-01215-w
Keywords: delta radiomics, pathological complete response, rectal cancer, MRgART, neoadjuvant chemoradiotherapy, MR-Linac
Citation: Wu J, Xiao J, Li Y, Wu F, Peng Q, Li C, Tang B and Orlandini LC (2023) Feasibility of delta radiomics–based pCR prediction for rectal cancer patients treated with magnetic resonance–guided adaptive radiotherapy. Front. Oncol. 13:1230519. doi: 10.3389/fonc.2023.1230519
Received: 29 May 2023; Accepted: 03 November 2023;
Published: 24 November 2023.
Edited by:
Lorenzo Placidi, Agostino Gemelli University Polyclinic (IRCCS), ItalyReviewed by:
Martina Murr, University of Tübingen, GermanyHuong Elena Tran, Agostino Gemelli University Polyclinic (IRCCS), Italy
Marica Vagni, Agostino Gemelli University Policlinic, Italy
Copyright © 2023 Wu, Xiao, Li, Wu, Peng, Li, Tang and Orlandini. 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: Bin Tang, jackytang86@163.com
†These authors have contributed equally to this work and share first authorship