Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 03 November 2022
Sec. Gastrointestinal Cancers: Colorectal Cancer
This article is part of the Research Topic Precision Medicine in Gastrointestinal Cancers View all 24 articles

Establishment and validation of novel MRI radiomic feature-based prognostic models to predict progression-free survival in locally advanced rectal cancer

Fei Xie&#x;Fei Xie1†Qin Zhao&#x;Qin Zhao1†Shuqi Li&#x;Shuqi Li1†Shuangshuang WuShuangshuang Wu2Jinli LiJinli Li3Haojiang LiHaojiang Li1Shenghuan ChenShenghuan Chen4Wu JiangWu Jiang1Annan DongAnnan Dong1Liqing WuLiqing Wu5Long LiuLong Liu6Huabin HuangHuabin Huang6Shuoyu XuShuoyu Xu1Yuanzhi ShaoYuanzhi Shao2Lizhi Liu,Lizhi Liu1,7Li Li*Li Li1*Peiqiang Cai*Peiqiang Cai1*
  • 1Sun Yat-sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangdong Key Laboratory of Nasopharyngeal Carcinoma Diagnosis and Therapy, Guangzhou, Guangdong, China
  • 2School of Physics, State Key Laboratory of Optoelectronic Materials and Technologies, Sun Yat-sen University, Guangzhou, China
  • 3The Second Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou, China
  • 4Department of Radiology, The Sixth Affiliated Hospital of Guangzhou Medical University, Qingyuan People’s Hospital, Qingyuan, China
  • 5Department of Radiology, Guangzhou Concord Cancer Center, Guangzhou, China
  • 6Department of Radiology, The Seventh Affiliated Hospital of Sun Yat-sen University, Shenzhen, China
  • 7Department of Radiology, The Third People’s Hospital of Shenzhen, Shenzhen, China

In locally advanced rectal cancer (LARC), an improved ability to predict prognosis before and after treatment is needed for individualized treatment. We aimed to utilize pre- and post-treatment clinical predictors and baseline magnetic resonance imaging (MRI) radiomic features for establishing prognostic models to predict progression-free survival (PFS) in patients with LARC. Patients with LARC diagnosed between March 2014 and May 2016 were included in this retrospective study. A radiomic signature based on extracted MRI features and clinical prognostic models based on clinical features were constructed in the training cohort to predict 3-year PFS. C-indices were used to evaluate the predictive accuracies of the radiomic signature, clinical prognostic models, and integrated prognostic model (iPostM). In total, 166 consecutive patients were included (110 vs. 56 for training vs. validation). Eleven radiomic features were filtered out to construct the radiomic signature, which was significantly related to PFS. The MRI feature-derived radiomic signature exhibited better prognostic performance than the clinical prognostic models (P = 0.007 vs. 0.077). Then, we proposed an iPostM that combined the radiomic signature with tumor regression grade. The iPostM achieved the highest C-indices in the training and validation cohorts (0.942 and 0.752, respectively), outperforming other models in predicting PFS (all P < 0.05). Decision curve analysis and survival curves of the validation cohort verified that iPostM demonstrated the best performance and facilitated risk stratification. Therefore, iPostM provided the most reliable prognostic prediction for PFS in patients with LARC.

Introduction

More than 100,000 cases of rectal cancer are diagnosed worldwide annually, and approximately 70% are locally advanced rectal cancer (LARC) (1). The standard treatment for LARC is neoadjuvant chemoradiotherapy (nCRT) followed by surgery (2, 3). Adjuvant chemotherapy is recommended for all patients with stage II/III rectal cancer after neoadjuvant radiochemotherapy and surgery (4). However, the effect of postoperative chemotherapy on survival remains controversial (5, 6). Therefore, pre- and post-treatment prediction models are needed to help determine treatment strategies and identify patients who may benefit from postoperative adjuvant chemotherapy.

LARC prognosis is based on the tumor–node–metastasis (TNM) staging system. Recently, nomograms were proposed to improve LARC prognosis prediction (710), using indexes such as age, carcinoembryonic antigen (CEA), carbohydrate antigen 19–9 (CA19-9), pathological tumor stage (ypT), pathological nodal stage (ypN), and tumor regression grade (TRG). However, contradictory results have been obtained with these models. For example, in one study, TRG was more efficient than ypTN stage in predicting the outcome (11), whereas in another study, ypTN stage contributed more (10), indicating that these nomograms are not robust enough for clinical application.

Radiomics, a tool that reveals underlying tumor heterogeneity using medical images (12, 13), serves as a strong prognostic predictor for malignancies (14). In rectal cancer, radiomics based on magnetic resonance imaging (MRI) is highly efficient in evaluating the tumor response to nCRT (14, 15) and can help identify non-responders (16) and pathologic complete responders (1720). Prediction models based on radiomic features have added predictive ability in combined models, enhancing the accuracy by up to 74% (sensitivity 58%, specificity 77%) (18). Combined models based on radiomics and clinical data can independently predict overall survival (21), disease-free survival (2224), and progression-free survival (PFS) (15) in LARC. However, pretreatment radiomics only reflects the cancer characteristics before surgery. Considering effective treatments may revise radiomics features, a study innovatively used delta radiomics based on 4 features to predict distant metastasis (DM) in LARC, obtaining a test set balanced accuracy, sensitivity and specificity of 78.5%, 71.4% and 85.7%, respectively (25). However, other prognostic clinical and histological indicators, including TRG, associated with the efficacy of nCRT on long term survival, which is important for determining the follow-up treatment after surgery, have not been considered in the radiomic prediction model. Thus, whether the combined model with radiomics and postoperative TRG data has an improved predictive ability for risk classification remains to be determined.

Here, we investigated the abilities of prognostic models based on radiomics, pre- and post-treatment clinical factors, and combination of radiomics and pre- and post-treatment clinical factors for predicting 3-year early PFS in LARC. Furthermore, we explored the internal correlations and differences among the models to determine the effect of combining different types of markers.

Methods

Patient selection

A total of 166 patients with LARC, based on pathological examination between March 2014 to August 2016, were enrolled (Supplementary Figure S1). The inclusion criteria were: (1) age ≥ 18 years; (2) newly diagnosed with LARC (staged on MRI as cT2–4 and/or N+) without distant metastasis and other malignancies; (3) treatment with preoperative nCRT; (4) high resolution pelvic MRI examination before nCRT; and (5) availability of complete electronic medical records and imaging data. Patients who did not complete nCRT were excluded. Considering only using radiomics feautures from pretreatment, pelvic MRI after treatment was not mandatory. Basic clinical information (sex, age, weight, height, BMI, clinical T stage, and clinical N stage) and laboratory indicators (routine blood tests, liver and renal function tests, blood glucose monitoring, C-reactive protein, serum lipid level, CEA, and CA19-9 levels) were collected before nCRT. This study was approved by the institutional ethics committee of our hospital. As this study was retrospective, the requirement for informed consent was exempted.

Treatment and follow-up

The preoperative treatment regimen included intensity-modulated radiation therapy (IMRT) and concurrent chemotherapy. IMRT doses of 50 Gy for gross tumor volume (primary tumors and enlarged LNs) and 45 Gy for clinical target volume were divided into 25 fractions. Two different chemotherapy regimens were used. Of the 166 patients, 91 received oxaliplatin (OXA) and capecitabine (CAPOX) for 3 weeks; 130 mg/m² OXA was intravenously administered on the first day, and 1,000 mg/m² capecitabine (CAP) was orally administered twice a day for the first two weeks. The remaining 75 patients received oral administration of 1,000 mg/m² CAP twice a day for the first 14 days. Radical rectal resection was performed 6–8 weeks after the completion of nCRT. After surgery, ypT stage, ypN stage, and TRG were evaluated.

The PFS, calculated at the endpoint, was defined as the interval from surgery to tumor progression, including local recurrence and/or metastasis or death. Follow-up visits were performed every 3–6 months in the first 2 years, then every 6 months in the following 3 years, and once a year thereafter.

MRI scanning and segmentation

Pretreatment pelvic MRI was performed using Trio Tim 3.0T (n = 74; Siemens Healthcare GmbH Henkestr) with two body Matrix coils and two spine Matrix coils or Discovery750 3.0T (n = 92; GE Healthcare) using an 8-channel phased array body coil in the supine position. Gadolinium-diethylenetriamine pentaacetic acid was injected as the contrast agent at a dosage of 0.1 mL/kg with a flow velocity of 3.0 mL/s. The scanning protocol included the axial, coronal, and sagittal T1-weighted (T1-w) images; T2-weighted (T2-w) images, axial short-axis T2-weighted images (short-axis T2-w), and contrast-enhanced T1-weighted (T1C-w) sequences. Short-axis T2–w, a thin section (3 mm) T2-weighted fast spin echo sequence acquired in a plane perpendicular to the long axis of the tumor (26), was helpful to precisely examine the tumor and its relationship with the intestinal wall, mesorectal fascia, vessels, and adjacent organs.

MR images were retrieved from the picture archiving and communication system and loaded onto AnalyzePro1 for manual segmentation. Two radiation oncologists with more than 10 years of experience outlined the whole-tumor volumes of interest, representing the contour of the tumor on each slice of all sequences separately.

Extraction of radiomic features and radiomic signature construction

Feature extraction in this study was performed using PyRadiomics2 (27). In total, 14,089 radiomic features were extracted from the axial T1C-w, T1-w, T2-w, and short-axis T2-w scans. The parameter settings for image preprocessing and radiomic feature extraction are presented in the Supplementary Information. To evaluate the effect of semi-automatic segmentations on the values of radiomic features, the inter-class correlation coefficient (ICC) was utilized. For this, 30 patients in the training cohort were randomly selected and segmented by two other radiation oncologists with more than 10 years of experience (Supplementary Figure S2A). The stability of each extracted feature was assessed by different expert radiologists. Stable radiomic features were defined as ICCs > 0.8 (Supplementary Figure S2B).

The RAD score was computed for each patient using a linear combination of selected features, weighted by their respective coefficients, and used to construct a radiomic signature. The potential association of the radiomic signature with 3-year early PFS was evaluated in the training cohort and then validated in the validation cohort.

Statistical analysis

The LASSO Cox regression method (Supplementary Information) was applied to select the most effective combination of prognostic features. The models were determined using the backward stepwise Akaike information criterion method, in which the least significant variables were removed one by one after fitting a full model with the candidate variables. The Mann-Whitney U test for continuous variables and the chi-square test for categorical variables were used to compare clinical characteristics between the training and validation cohorts. Considering multiple factors (including clinical T stage, clinical N stage, TRG, CEA, CA19-9, GLO, ypT stage, ypN stage), univariate and multivariate analyses were performed with the Cox proportional hazards model, and the hazard ratios (HRs) and 95% confidence intervals (CIs) were calculated. Harrell’s concordance indices were used to assess the predictive power of each model. Statistical analyses were performed using R version 4.0.23. R packages, including glmnet, caret, survival, rms, Hmisc, corrplot, pheatmap, and rmda, were used. All statistical tests were two-sided, and P < 0.05 was considered statistically significant.

Results

A total of 166 consecutive patients (median age, 58 [interquartile range, 49‒64] years; 117 [70.48%] men) were included in the analysis. The patients were randomly divided in a 2:1 ratio, with 110 in the training cohort and 56 in the validation cohort. Table 1 summarizes the clinical characteristics of the patients with LRAC in the training and validation cohorts. There were no differences in those characteristics between the training and validation cohorts (P = 0.202–0.930).

TABLE 1
www.frontiersin.org

Table 1 Clinical and pathological characteristics of patients in the training and validation cohorts.

The median follow-up duration was 32.3 months (range, 2.6–58.8 months). During the last follow-up, disease progression was confirmed in 24 patients (21.8%) in the training cohort and 14 (25.0%) in the validation cohort (P = 0.708).

After univariate analysis in the training cohort, candidate variables with a P value < 0.1 were included in the multivariate model analysis. Pretreatment variables, including clinical T stage, CEA, CA 19-9, and globulin, and post-treatment predictors, including pathologic T stage and TRG, were selected as summarized in Supplementary Table S1. The pretreatment clinical prognostic model (PreM) for 3-year PFS prediction was established based on the four pretreatment variables. Multivariate analysis identified CEA and globulin as independent predictors (Supplementary Table S2). The clinical stage prognostic model was built based on the clinical T and N stages.

The C-indices of PreM and the clinical stage prognostic model were 0.627 (95% confidence interval [CI]: 0.572–0.682) and 0.578 (95% CI: 0.522–0.634) in the training cohort, and 0.552 (95% CI: 0.482–0.622) and 0.611 (95% CI: 0.531–0.691) in the validation cohort, respectively (Table 2). The nomograms and corresponding calibration curves for the probability of 3-year PFS generated using the pretreatment clinical prognostic models are shown in Supplementary Figure S3.

TABLE 2
www.frontiersin.org

Table 2 C-index for each prognostic model for survival prediction.

Additionally, by integrating the clinical variables for pre- and post-treatment, we established post-treatment prognostic models. The pathologic stage prognostic model was built based on the pathologic T and N stages. Stepwise multivariable analyses identified pathologic T stage, CEA, and globulin as independent predictors for PFS in the first post-treatment prognostic model (PostM1). In the validation cohort, PostM1 showed a higher predictive power for PFS than PreM (P = 0.307) and the pathologic stage prognostic model (P = 0.156), but it was slightly lower than that of the clinical stage prognostic model (P = 0.429). When the influences of clinical T stage and pathologic T stage were excluded, TRG, CEA, and globulin remained significant for PFS after performing multivariate Cox regression. Consequently, PostM2 was built with these three factors. The C-indices of the three post-treatment prognostic models are summarized in Table 2. The nomograms and corresponding calibration curves for the 3-year PFS probability generated using the post-treatment clinical prognostic models are shown in Figure 1 and Supplementary Figure S4. The multivariable Cox regression results of the pre- and post-treatment clinical prognostic models are summarized in Supplementary Table S2. Notably, compared with PreM without TRG, PostM2 achieved a better predictive performance (P = 0.356), which indicates that TRG is an important prognostic factor for predicting PFS.

FIGURE 1
www.frontiersin.org

Figure 1 Nomogram for 3-year progression-free survival (PFS) in (A) iPostM and (C) PostM2. The nomogram allows the user to determine the probability of 3-year PFS corresponding to a patient’s combination of covariates. Calibration curves for predicting 3-year PFS in (B) iPostM and (D) PostM2 in the training and validation cohorts. The y-axis shows observed survival estimated using the Kaplan-Meier method, and the x-axis shows predicted survival calculated using the prognostic model. The closer fit to the diagonal dotted line indicates a better assessment. iPostM, integrated prognostic model combining tumor regression grade (TRG) and radiomic signature; PostM2, post-treatment clinical prognostic model without pathologic stage.

In the training cohort, we selected 11 radiomic features based on MRI that were significantly associated with PFS (Supplementary Table S3). The detailed selection process and LASSO results are shown in Supplementary Figures S5 and S6, respectively. The formula for the Rad score calculation of the radiomic prognostic model is shown in the Supplementary Information. In the training cohort, the radiomic signature yielded a C-index of 0.937 (95% CI: 0.917–0.957). The good prognostic performance of this radiomic signature was validated with a corresponding C-index of 0.730 (95% CI: 0.651‒0.809) in the validation cohort. The radiomic nomogram showed significant improvement compared to the clinical prognostic models (P = 0.007‒0.039), except for PostM2 (P = 0.077) (Table 2). Thus, the developed radiomic signature was more accurate than the clinical prognostic models for evaluating 3-year PFS.

Next, we built two integrated prognostic models that combined the radiomic signature based on MRI features with important pre- and post-treatment clinical factors. Using the multivariate Cox proportional hazard model based on pretreatment clinical factors and radiomic signature, we found that only the radiomic signature remained significant for PFS after adjusting for various cofactors. The integrated PostM (iPostM) was constructed using the radiomic signature and TRG (Table 3). The iPostM showed significant improvement compared to the radiomic signature in terms of evaluating 3-year PFS (C-index: 0.752; 95% CI: 0.684‒0.820), with a P value < 0.05 (Table 2).

TABLE 3
www.frontiersin.org

Table 3 Multivariate Cox regression analysis of the final integrated model.

The result of iPostM is visually represented by a nomogram, as shown in Figure 1A. The calibration curve for the 3-year PFS probability showed good agreement between the evaluation based on nomogram and actual survival (Figure 1B).

Then, we calculated the risk scores of all the prognostic models for each patient in both the training and validation cohorts and then classified the patients into categories: low-risk (patients with a score < 0) and high-risk (patients with a score ≥ 0), with zero as the risk score cutoff. The survival curves between patients in the low- and high-risk categories, generated using clinical and radiomic prognostic models in the training and validation cohorts, are shown in Supplementary Figures S7 and S8, and Figure 2. The patients with disease progression after treatment were concentrated in the high score area, and the survival curve showed good prognostic stratification of patients in the low-and high-risk groups in the validation cohort of iPostM. However, such trends were not observed in the validation cohort of other clinical prognostic models and radiomic signature.

FIGURE 2
www.frontiersin.org

Figure 2 Stratified Kaplan-Meier analyses of the prognostic models to estimate 3-year progression-free survival (PFS) in various risk stratification subgroups in the validation cohort. Patients with high and low risks of PFS were stratified by the prognostic models (A) Pathologic TN, (B) PostM2, (C) Radiomics signature, (D) iPostM models. Only the iPostM could stratify patients into high- and low-risk groups based on significantly different 3-year PFS rates (P < 0.05). The log-rank test was used to calculate P values. Pathologic TN, pathologic stage prognostic model; PostM2, post-treatment clinical prognostic model without pathologic stage; iPostM, integrated prognostic model combining tumor regression grade (TRG) and radiomic signature.

Decision curve analysis of the validation cohort in all the prognostic models showed that iPostM was the most efficient (Figure 3). The iPostM exhibited the highest efficacy especially in high risk areas. Using a heatmap to determine the association between the radiomic signature and clinical data, we found no significant correlation between the radiomic signature and clinical factors (Supplementary Figure 9). In contrast, TRG and pathologic T stage showed a strong correlation (r = 0.69).

FIGURE 3
www.frontiersin.org

Figure 3 The DCA curves of the nomograms compared for 3-year progression-free survival (PFS) in the validation cohort. The x-axis represents the threshold probabilities, and the y-axis indicates the net benefit. The net benefit is calculated by adding the benefits (true-positive results) and subtracting the risks (false-positive results), with the latter weighted by a factor related to the harm of an undetected cancer relative to the harm of unnecessary treatment. Clinical TN, clinical stage prognostic model; PreM, pretreatment clinical prognostic model; ypTN, pathologic stage prognostic model; PostM1, post-treatment clinical prognostic model; PostM2, post-treatment clinical prognostic model without pathologic stage; iPostM, integrated prognostic model combining tumor regression grade (TRG) and radiomic signature.

Discussion

In the present study, we established effective prognostic models for predicting 3-year PFS in patients with LARC. Further, we proposed an integrated prognostic model that combined the radiomic signature and TRG and yielded the highest C-index with a value of 0.752 in the validation cohort, outperforming the radiomic signature and all other clinical prognostic models. There was a strong correlation between TRG and pathologic T stage (r = 0.69). Notably, only the integrated post-treatment prognostic model could stratify patients into high- and low-risk groups based on significantly different 3-year PFS rates.

The prognostic models, based on only clinical and pathological factors, showed relatively weak predictive performance, indicating the need to find more useful markers. Additionally, the C-indices of PostMs were higher than those of PreM, possibly because pathological information can better reflect the preoperative state of the tumor. The models with TRG, used as a qualitative evaluation of tumor cells replaced by fibrosis that can reflect the sensitivity of tumors to nCRT, showed higher C-indexes, confirming the findings of previous studies (11, 2830). TRG was significantly correlated with ypT, which may lead to the absence of ypT in multi-parameter models. Different from some previous studies (10, 24, 31), the discriminatory power of ypN was weakened and not included in the final model, which may attributed to the heterogenous distribution of ypN stage, the adding of TRG and radiomics in our model. For example, Cui et al. (24) included ypN in the prediction model instead of TRG. In their study, the percentages of ypN 0 in the training and validation groups were 60.3% and 50.9%, respectively, which were lower than those in our study (81.33% and 80.36%) and previous findings (96.9% and 94.6%) (9). Thus, to make radiomics-based model more robust, larger sample size and external validation are needed to confirm the results.

Our study indicates that radiomics is an independent prognostic factor for predicting 3-year early PFS in LARC, which is in accordance with a previous study (23). MRI is the standard imaging method for post-nCRT evaluation in LARC, and MRI-based radiomics, can provide minable information from conventional medical images and dig out quantitative features which can reflect tumor heterogeneity, other intrinsic characteristics and microenvironment related with individualized biological behavior of tumor (1214, 25, 32). It is reported that higher levels of radiomics heterogeneity (ie, higher entropy) was associated with worse response to treatment and/or survival (33, 34), which may caused by the occurrence of constant complex mutations within a tumor to remain resistant to treatment (34). The correlation coefficient between the RAD score and clinical data was low, indicating that radiomic features include details derived from images rather than the clinical TNM stage derived from the macro level (35). Further, compared clinical and pathological models with underfitting, the iPostM obtained enough fitting in training cohort and achieved the highest C-index in the validation cohort, in which the RAD score contributed significantly to the nomogram. The reason may be that the RAD score comprises multiple underlying tumor characteristics associated with disease risks (32, 36, 37), whereas TRG is merely a pathological signature. Additionally, radiomics may decrease the discriminatory power of previously proven independent prognostic indicators (8, 38, 39) as shown in Supplementary Table S4. Thus, radiomics was the most influential indicator in the nomogram, followed by TRG and others. Based on the relatively comprehensive information, iPostM may identify patients in high-risk groups and suggest the administration of adjuvant chemotherapy after TME (4) to reduce the risk of occurrence of distant metastasis and recurrence.

Based on the strengths of similar prognostic prediction models for LARC (23, 24), we included multiparametric MR images for extracting radiomic features, with the addition of T1-w, T2-w, and short-axis T2-w, while other studies included relatively fewer sequences. Multiple MRI sequences can detect anatomical details and provide more specific histological information, such as necrosis, cystic degeneration, hemorrhage, and tumor angiogenesis (40). Notably, previous studies have proven that the predictive performance of a radiomic model derived from multi-modal MRI is superior to that based on mono-modal MRI (8, 38). Hence, our model could improve the prognosis ability for patients with LARC.

Nevertheless, there are some limitations to our study. First, it was a relatively small retrospective study performed and validated at a single center with short follow-up duration; however, we enrolled the patients consecutively to reduce underlying selective bias. Second, we only extracted signatures from preoperative primary tumors, lacking lymph nodes, functional MRI images, such as diffusion weighted imaging and apparent diffusion coefficient, which would provide more signatures and inner information. In addition, the value of post-treatment radiomics features need to be further confirmed. Third, we focused on random combinations of imaging features with clinical data rather than genetic heterogeneity. Furthermore, TRG, which is included in our best performing model, was not available to predict the prognosis for the “wait and see” patients. Future studies need to be carried out to validate the prognostic value of our iPostM model in multiple centers with longer follow-up duration.

In summary, we developed pre- and post-treatment prediction models based on clinical and radiomic features. Post-treatment prognostic models with postoperative pathological factors showed better predictive performance than pretreatment prognostic models, and TRG was important for predicting the 3-year PFS of LARC. The multi-modal MRI radiomic model with improved predictive ability could act as a pretreatment-independent prognostic factor for LARC and assist clinicians in determining appropriate neoadjuvant chemoradiotherapy regimens. In addition, the integrated post-treatment prognostic model has potential for recognizing high-risk patients who may benefit from postoperative adjuvant therapy. In future studies, we will validate the performance of our model and explore its clinical applications.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the institutional ethics committee of Sun Yat-sen University Cancer Center. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

Conception and design of this study: FX, QZ and SL. Patients providing: FX, QZ, SC and WJ. Data collection: SL, AD, JL, LW, LLiu and HH. Data analysis and interpretation: SW, HL. Study guidance: SX, YS, LZL, LL and PC. Manuscript writing and revision: All authors. Final approval of the manuscript: All authors.

Funding

This study has received funding by the Key-Area Research and Development of Guangdong Province (Grant No. 2020B1111190001), Science and Technology Planning Project of Guangzhou City, China (Grant No. 201907010043) and the National Natural Science Foundation of China (Grant No. 61975244, 82171906, and 81902638).

Acknowledgments

We also acknowledge the help of English language editing from Editage (www.editage.com).

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.901287/full#supplementary-material

Footnotes

  1. ^ https://analyzedirect.com/analyzepro/
  2. ^ https://github.com/Radiomics/pyradiomics
  3. ^ http://www.r-project.org/

References

1. Liu Z, Zhang XY, Shi YJ, Wang L, Zhu HT, Tang Z, et al. Radiomics analysis for evaluation of pathological complete response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer. Clin Cancer Res (2017) 23:7253–62. doi: 10.1158/1078-0432.CCR-17-1038

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Zhu S, Brodin NP, English K, Ohri N, Chuy JW, Rajdev LN, et al. Comparing outcomes following total neoadjuvant therapy and following neoadjuvant chemoradiation therapy in patients with locally advanced rectal cancer. EClinicalmedicine (2019) 16:23–9. doi: 10.1016/j.eclinm.2019.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Benson AB, Venook AP, Al-Hawary MM, Arain MA, Chen YJ, Ciombor KK, et al. NCCN Guidelines Insights: Rectal Cancer, Version 6.2020. J Natl Compr Canc Netw (2020) 18:806–15. doi: 10.6004/jnccn.2020.0032

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Benson AB, Venook AP, Al-Hawary MM, Cederquist L, Chen YJ, Ciombor KK, et al. Rectal Cancer, Version 2.2018, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw (2018) 16:874–901. doi: 10.6004/jnccn.2018.0061

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Biagi JJ, Raphael MJ, Mackillop WJ, Kong W, King WD, Booth CM. Association between time to initiation of adjuvant chemotherapy and survival in colorectal cancer: a systematic review and meta-analysis. JAMA (2011) 305:2335–42. doi: 10.1001/jama.2011.749

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Maas M, Nelemans PJ, Valentini V, Crane CH, Capirci C, Rödel C, et al. Adjuvant chemotherapy in rectal cancer: Defining subgroups who may benefit after neoadjuvant chemoradiation and resection: a pooled analysis of 3,313 patients. Int J Cancer (2015) 137:212–20. doi: 10.1002/ijc.29355

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Liu J, Huang X, Yang W, Li C, Li Z, Zhang C, et al. Nomogram for predicting overall survival in stage II-III colorectal cancer. Cancer Med (2020) 9:2363–71. doi: 10.1002/cam4.2896

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Song J, Chen Z, Huang D, Wu Y, Lin Z, Chi P, et al. Nomogram predicting overall survival of resected locally advanced rectal cancer patients with neoadjuvant chemoradiotherapy. Cancer Manag Res (2020) 12:7375–82. doi: 10.2147/CMAR.S255981

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sun Y, Zhang Y, Wu X, Lin H, Lu X, Huang Y, et al. Prognostic significance of neoadjuvant rectal score in locally advanced rectal cancer after neoadjuvant chemoradiotherapy and construction of a prediction model. J Surg Oncol (2018) 117:737–44. doi: 10.1002/jso.24907

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Li ZH, Xie PY, Zhang DF, Li YJ, Wu L, Dong J, et al. Nomogram for predicting disease-free survival among a multicenter cohort of Chinese patients with locally advanced rectal cancer. Cancer Manag Res (2019) 11:2471–83. doi: 10.2147/CMAR.S196614

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Xu L, Cai S, Xiao T, Chen Y, Qiu H, Wu B, et al. Prognostic significance of tumour regression grade after neoadjuvant chemoradiotherapy for a cohort of patients with locally advanced rectal cancer: an 8-year retrospective single-institutional study. Colorectal Dis (2017) 19:O263–71. doi: 10.1111/codi.13757

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Capobianco E, Dominietto M. From medical imaging to radiomics: role of data science for advancing precision health. J Pers Med (2020) 10:15. doi: 10.3390/jpm10010015

CrossRef Full Text | Google Scholar

13. Lubner MG, Smith AD, Sandrasegaran K, Sahani DV, Pickhardt PJ. CT texture analysis: definitions, applications, biologic correlates, and challenges. Radiographics (2017) 37:1483–503. doi: 10.1148/rg.2017170056

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong E, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol (2017) 14:749–62. doi: 10.1038/nrclinonc.2017.141

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wang J, Liu X, Hu B, Gao Y, Chen J, Li J. Development and validation of an MRI-based radiomic nomogram to distinguish between good and poor responders in patients with locally advanced rectal cancer undergoing neoadjuvant chemoradiotherapy. Abdom Radiol (NY) (2020) 46:1805–15. doi: 10.1007/s00261-020-02846-3

PubMed Abstract | CrossRef Full Text | Google Scholar

16. 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 (2020) 12:1894. doi: 10.3390/cancers12071894

CrossRef Full Text | Google Scholar

17. Shaish H, Aukerman A, Vanguri R, Spinelli A, Armenta P, Jambawalikar S, et al. Radiomics of MRI for pretreatment prediction of pathologic complete response, tumor regression grade, and neoadjuvant rectal score in patients with locally advanced rectal cancer undergoing neoadjuvant chemoradiation: an international multicenter study. Eur Radiol (2020) 30:6263–73. doi: 10.1007/s00330-020-06968-6

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Petkovska I, Tixier F, Ortiz EJ, Golia Pernicka JS, Paroder V, Bates DD, et al. Clinical utility of radiomics at baseline rectal MRI to predict complete response of rectal cancer after chemoradiation therapy. Abdom Radiol (NY) (2020) 45:3608–17. doi: 10.1007/s00261-020-02502-w

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Horvat N, Veeraraghavan H, Khan M, Blazic I, Zheng J, Capanu M, et al. MR imaging of rectal cancer: radiomics analysis to assess treatment response after neoadjuvant therapy. Radiology (2018) 287:833–43. doi: 10.1148/radiol.2018172300

PubMed Abstract | CrossRef Full Text | Google Scholar

20. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lovinfosse P, Polus M, Van Daele D, Martinive P, Daenen F, Hatt M, et al. FDG PET/CT radiomics for predicting the outcome of locally advanced rectal cancer. Eur J Nucl Med Mol Imaging (2018) 45:365–75. doi: 10.1007/s00259-017-3855-5

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Wang J, Shen L, Zhong H, Zhou Z, Hu P, Gan J, et al. Radiomics features on radiotherapy treatment planning CT can predict patient survival in locally advanced rectal cancer patients. Sci Rep (2019) 9:15346. doi: 10.1038/s41598-019-51629-4

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Meng Y, Zhang Y, Dong D, Li C, Liang X, Zhang C, et al. Novel radiomic signature as a prognostic biomarker for locally advanced rectal cancer. J Magn Reson Imaging (2018). doi: 10.1002/jmri.25968

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (2021) 154:161–9. doi: 10.1016/j.radonc.2020.09.039

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Chiloiro G, Rodriguez-Carnero P, Lenkowicz J, Casà C, Masciocchi C, Boldrini L, et al. Delta Radiomics Can Predict Distant Metastasis in Locally Advanced Rectal Cancer: The Challenge to Personalize the Cure. Front Oncol (2020) 10:595012. doi: 10.3389/fonc.2020.595012

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Nougaret S, Reinhold C, Mikhael HW, Rouanet P, Bibeau F, Brown G. The use of MR imaging in treatment planning for patients with rectal carcinoma: have you checked the “DISTANCE”? Radiology (2013) 268:330–44. doi: 10.1148/radiol.13121361

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. Cancer Res (2017) 77:e104–7. doi: 10.1158/0008-5472.CAN-17-0339

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Park IJ, You YN, Agarwal A, Skibber JM, Rodriguez-Bigas MA, Eng C, et al. Neoadjuvant treatment response as an early response indicator for patients with rectal cancer. J Clin Oncol (2012) 30:1770–6. doi: 10.1200/JCO.2011.39.7901

PubMed Abstract | CrossRef Full Text | Google Scholar

29. McCoy MJ, Hemmings C, Hillery S, Penter C, Bulsara MK, Zeps N, et al. Neoadjuvant chemoradiotherapy for rectal cancer: how important is tumour regression? ANZ J Surg (2017) 87:E233–9. doi: 10.1111/ans.13394

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mace AG, Pai RK, Stocchi L, Kalady MF. American Joint Committee on Cancer and College of American Pathologists regression grade: a new prognostic factor in rectal cancer. Dis Colon Rectum (2015) 58:32–44. doi: 10.1097/DCR.0000000000000266

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Duchalais E, Glyn Mullaney T, Spears GM, Kelley SR, Mathis K, Harmsen WS, et al. Prognostic value of pathological node status after neoadjuvant radiotherapy for rectal cancer. Br J Surg (2018) 105(11):1501–9. doi: 10.1002/bjs.10867

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology (2016) 278:563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Dagogo-Jack I, Shaw AT. Tumour heterogeneity and resistance to cancer therapies. Nat Rev Clin Oncol (2017) 15:81–94. doi: 10.1038/nrclinonc.2017.166

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Marusyk A, Almendro V, Polyak K. Intra-tumour heterogeneity: a looking glass for cancer? Nat Rev Cancer (2012) 12:323. doi: 10.1038/nrc3261

PubMed Abstract | CrossRef Full Text | Google Scholar

35. 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:43. doi: 10.1186/s13014-019-1246-8

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Rödel C, Sauer R. Integration of novel agents into combined-modality treatment for rectal cancer patients. Strahlenther Onkol (2007) 183:227–35. doi: 10.1007/s00066-007-9000-9

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Rödel C, Graeven U, Fietkau R, Hohenberger W, Hothorn T, Arnold D, et al. Oxaliplatin added to fluorouracil-based preoperative chemoradiotherapy and postoperative chemotherapy of locally advanced rectal cancer (the German CAO/ARO/AIO-04 study): final results of the multicentre, open-label, randomised, phase 3 trial. Lancet Oncol (2015) 16:979–89. doi: 10.1016/S1470-2045(15)00159-X

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Zheng Z, Wang X, Huang Y, Lu X, Chi P. Predictive value of carbohydrate antigen 19-9 level changes in patients with locally advanced rectal cancer treated with neoadjuvant chemoradiotherapy. Colorectal Dis (2020) 22:2068–77. doi: 10.1111/codi.15355

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Moon SH, Kim DY, Park JW, Oh JH, Chang HJ, Kim SY, et al. Can the new American Joint Committee on Cancer staging system predict survival in rectal cancer patients treated with curative surgery following preoperative chemoradiotherapy? Cancer (2012) 118:4961–8. doi: 10.1002/cncr.27507

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Horvat N, Rocha CCT, Oliveria BC, Petkovska I, Gollub MJ. MRI of rectal cancer: tumor staging, imaging techniques, and management. Radiographics (2019) 39:367–87. doi: 10.1148/rg.2019180114

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: magnetic resonance imaging, radiomics, prognostic models, locally advanced rectal cancer, clinical predictors

Citation: Xie F, Zhao Q, Li S, Wu S, Li J, Li H, Chen S, Jiang W, Dong A, Wu L, Liu L, Huang H, Xu S, Shao Y, Liu L, Li L and Cai P (2022) Establishment and validation of novel MRI radiomic feature-based prognostic models to predict progression-free survival in locally advanced rectal cancer. Front. Oncol. 12:901287. doi: 10.3389/fonc.2022.901287

Received: 21 March 2022; Accepted: 20 October 2022;
Published: 03 November 2022.

Edited by:

Afsaneh Barzi, City of Hope National Medical Center, United States

Reviewed by:

Seung-Bum Ryoo, Seoul National University Hospital, South Korea
Giuditta Chiloiro, Agostino Gemelli University Polyclinic (IRCCS), Italy

Copyright © 2022 Xie, Zhao, Li, Wu, Li, Li, Chen, Jiang, Dong, Wu, Liu, Huang, Xu, Shao, Liu, Li and Cai. 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: Peiqiang Cai, caipq@sysucc.org.cn; Li Li, lil@sysucc.org.cn

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.