- 1Department of Radiology, Affiliated Dongyang Hospital of Wenzhou Medical University, Dongyang, China
- 2Department of Radiology, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 3Department of Radiology, Jiaxing Hospital of Traditional Chinese Medicine, Jiaxing, China
- 4Key Laboratory of Diagnosis and Treatment of Digestive System Tumors of Zhejiang Province, Ningbo, China
- 5Cancer Center, Zhejiang University, Hangzhou, China
Background: This study aimed to establish an effective model for preoperative prediction of tumor deposits (TDs) in patients with rectal cancer (RC).
Methods: In 500 patients, radiomic features were extracted from magnetic resonance imaging (MRI) using modalities such as high-resolution T2-weighted (HRT2) imaging and diffusion-weighted imaging (DWI). Machine learning (ML)-based and deep learning (DL)-based radiomic models were developed and integrated with clinical characteristics for TD prediction. The performance of the models was assessed using the area under the curve (AUC) over five-fold cross-validation.
Results: A total of 564 radiomic features that quantified the intensity, shape, orientation, and texture of the tumor were extracted for each patient. The HRT2-ML, DWI-ML, Merged-ML, HRT2-DL, DWI-DL, and Merged-DL models demonstrated AUCs of 0.62 ± 0.02, 0.64 ± 0.08, 0.69 ± 0.04, 0.57 ± 0.06, 0.68 ± 0.03, and 0.59 ± 0.04, respectively. The clinical-ML, clinical-HRT2-ML, clinical-DWI-ML, clinical-Merged-ML, clinical-DL, clinical-HRT2-DL, clinical-DWI-DL, and clinical-Merged-DL models demonstrated AUCs of 0.81 ± 0.06, 0.79 ± 0.02, 0.81 ± 0.02, 0.83 ± 0.01, 0.81 ± 0.04, 0.83 ± 0.04, 0.90 ± 0.04, and 0.83 ± 0.05, respectively. The clinical-DWI-DL model achieved the best predictive performance (accuracy 0.84 ± 0.05, sensitivity 0.94 ± 0. 13, specificity 0.79 ± 0.04).
Conclusions: A comprehensive model combining MRI radiomic features and clinical characteristics achieved promising performance in TD prediction for RC patients. This approach has the potential to assist clinicians in preoperative stage evaluation and personalized treatment of RC patients.
1 Introduction
Colorectal cancer (CRC) is the third most common malignancy and second leading cause of death worldwide. In particular, rectal cancer (RC) accounts for approximately one-third of CRC (1). A tumor deposit (TD) is defined as a discontinuous cancerous nodule located in the mesenteric fascia without obvious nodal or vascular features. The median incidence of TDs in CRC patients is approximately 21.3% (2). Positive TDs can elevate clinical stages of RC patients. RC patients with positive TDs are classified as N1c and treated as clinical stage III, in the absence of nodal metastases. The efficacy of adjuvant chemotherapy in stage III colon cancer had been widely recognized in previous studies, For the TD patients, with the receipt of chemotherapy had decreased risk of cancer-specific mortality compared with those not (3). However, only 52% of TD-positive and lymph node (LN)-negative patients receive preoperative adjuvant chemotherapy (4). Therefore, early identification of TDs is important and valuable for stage evaluation and treatment planning.
Magnetic resonance imaging (MRI) is considered the most reliable imaging modality for the initial pretreatment evaluation of patients with RC, including the assessment of TN staging, circumferential resection margin (CRM), and extramural vascular invasion (EMVI). Moreover, it can assist in the clinical preoperative management of RC patients, determination of surgical scope, and assessment of treatment response to neoadjuvant therapy (5, 6). However, MRI is limited in determining tumor spread in the mesorectum. Gröne et al. used a small diameter of 5 mm as the statistical threshold to determine N staging. With this value, the sensitivity, specificity, and accuracy of MRI staging were 72%, 45.7%, and 56.7%, respectively (7). Langman et al. showed that mesenteric nodules <3 mm had a 28% probability of malignancy (8). These studies focused on the presence of tumor spread in mesorectal nodules, either LN metastasis or TDs. However, the evaluation of malignant LNs alone is insufficient to reflect the actual spread of RC in the mesentery. According to the current European Society for Medical Oncology (ESMO) preoperative risk assessment criteria, patients with TDs are classified into a high-risk group with a worse prognosis (9). A previous study confirmed that the presence of TDs is an independent risk factor for the prognosis of patients with RC (10). An analysis of two prognostic studies in N0 and N1c stages showed a significant difference in the five-year survival rates (N0, 91.5%; N1c, 37%) (11). Therefore, preoperative determination of tumor deposition status in patients with RC is essential for optimal treatment.
Currently, the presence of TDs is determined by pathological analysis after radical tumor resection. However, this method is invasive and can be performed only postoperatively. According to a previous study, MRI can help preoperatively and identify TDs and LN metastases, as these lesions appear to have distinguishable imaging characteristics on MRI (12). However, TDs are usually less than 5 mm in diameter, and identifying such small nodules and accurately assess the characteristics of the nodule can be challenging for radiologists who are already overburdened in reading MRI in daily practice.
By extracting vast amounts of quantitative features from imaging and providing non-visual information that indicates the biological behavior of tumors, radiomics has gained popularity for the non-invasive prediction of clinical or prognostic features of tumors, such as T staging of RCs, LN status, vascular and nerve invasion, distant metastasis, and pathological complete response to neoadjuvant chemotherapy (13–20). Meanwhile, artificial intelligence including deep neural networks has demonstrated high performance in the analysis of medical images (21–23), providing cancer risk assessment, recurrence, and survival predictions with higher accuracy than human experts. Recently, several radiomic models have been developed based on ultrasound (US), computed tomography (CT), and MRI to preoperatively predict TDs in patients with RC (10, 24, 25). However, the sample sizes in these studies were relatively small. Furthermore, the MRI study carried out by Yang et al. (25) only extracted the radiomic features from high-resolution T2 weighted (HRT2) MRI, whereas functional MRI, such as diffusion-weighted imaging (DWI), carries more information on the heterogeneity of tumors. Currently, there is a lack of functional MRI-based deep-learning (DL) radiomics research in this field.
This study aimed to develop an MR-based DL radiomic model for preoperative TD prediction in a larger cohort with higher prediction accuracy. This model extracts radiomic features from both HRT2 and DWI images and integrates clinical factors into TD prediction.
2 Methods
2.1 Patient characteristics
The records of 784 consecutive RC patients who underwent preoperative MRI and radical surgery between 2013 and 2020 at Sir Run Run Shaw Hospital affiliated with Zhejiang University School of Medicine were reviewed retrospectively. The local institutional review board approved this study and provided a waiver of consent. The inclusion criteria were: (a) pathologically confirmed primary RC; (b) no neoadjuvant chemotherapy or radiotherapy before surgery; and (c) tumor visible in at least three sequential slices of HRT2 MRI. The exclusion criteria were as follows: (a) inadequate MRI quality due to intractable artifacts, including HRT2 (n = 113) and DWI (n = 78); (b) tumors not visible in HRT2 images (n = 14); (c) carcinoembryonic antigen (CEA) and carbohydrate antigen 19-9 (CA19-9) levels not obtained (n = 44); (d) lack of tissue differentiation grading in pathology reports (n = 5); and (e) co-occurrence of other digestive system malignancies (n = 30). Ultimately, 500 patients were enrolled in this study (Figure 1).
Figure 1 Flowchart of patient selection and TD distribution in the study. TD, tumor deposition; CEA, carcinoembryonic antigen; CA19-9, carbohydrate antigen 19-9.
2.2 Clinical characteristics and pathological criteria
Clinical characteristics, including sex, age, body mass index (BMI), and carbohydrate antigen 19-9 (CA19-9) and carcinoembryonic antigen (CEA) levels, were collected from electronic medical records. Histological grading, pathological tumor node metastasis (pTNM) staging, LNs, TDs, vascular and nerve invasion, and other clinical data were obtained from pathological reports. The eighth edition of the AJCC staging system was used as a reference for the pTNM staging. TDs are defined as discrete tumor foci in the pericolic or perirectal fat, without histological evidence of residual lymph node or identifiable vascular or neural structures.The distance between the tumor and anus, TN staging, CRM, and EMVI based on the MRI were obtained from standardized reports of the picture archiving and communication system (PACS). The criteria for determining LN positivity on MRI were based on the latest recommendations of the 2016 European Society of Gastrointestinal and Abdominal Radiology consensus meeting (26). MRI tumor length was defined as the T2 sagittal tumor length. The distance from the anus was defined as the distance from the most inferior boundary of the tumor to the subcutaneous edge of the anus. The criterion for EMVI positivity was tumor invasion of the extramural vessels, with or without vessel dilatation (27). CRM positivity was defined as a tumor location within 1 mm of the mesorectal fascia, including suspicious LNs, TDs, tumor expansion, and EMVI (28). Unclear or missing information in the MRI reports were labeled and finally confirmed by a radiologist with nine years of working experience.
2.3 MRI scanning
MRI acquisitions were performed using the following 3.0-T MRI scanners: Signa HDxt (GE Healthcare, Chicago, IL), Discovery MR750w (GE Healthcare), and MAGNETOM Skyra (Siemens Healthineers, Erlangen, Germany). The MRI protocol consisted of one axial HRT2 MRI sequence and one DWI acquisition obtained using b-values of 0 and 1,000 (or 800) s/mm². No intravenous contrast agents were administered. Details of the MRI acquisition parameters are listed in Table 1.
2.4 Tumor segmentation and processing
Before image segmentation, patient-sensitive information was anonymized. The primary tumor region (3D volume) was semi-manually segmented on axial HRT2 and DWI images by a junior radiologist (with more than three years of experience in radiology) using an open-source software tool (ITK-SNAP 3.8; www.itksnap.org) (29). Automatic tumor segmentation using a CE-net-based DL segmentation model (30) was performed on the axial HRT2 images to assist radiologists. All segmentation masks were reviewed by a senior radiologist (with more than five years of experience in radiology) and finally confirmed by another senior radiologist (with more than 10 years of experience in radiology). Disagreements were resolved through discussion.
2.5 Extraction of features
International Biomarker Standardization Initiative (IBSI)-compliant radiomic features were extracted separately for the HRT2 images and DWI images using PyRadiomics, an open-source Python package (version 2.1.2, https://pyradiomics.readthedocs.io) (31). Before feature extraction, z-score normalization of the MRI signal intensities for both the HRT2 and DWI images was performed using PyRadiomics. Consequently, 564 features were obtained for each HRT2 and DWI image, including 13 first-order statistics, 35 shape features, 9 orientation features, and 507 texture features, such as the gray-level co-occurrence matrix, gray-level size zone matrix, gray-level run-length matrix, gray-level dependence matrix, neighborhood gray-tone difference matrix, Gabor filter, Laplacian of Gaussian filter, local binary patterns, and local phase and vascularity filters. A variance test was performed on the extracted features to remove features with a low variance (<0.01). A t-test was used to estimate the radiomic features that were significantly correlated with TDs. Features with p < 0.05 were considered significant features for model development. A detailed flowchart of this process is shown in Figure 2.
Figure 2 Flowchart describing the methods used in developing the clinical-DL radiomics model for TD prediction in patients with RC. DL, deep learning; TD, tumor deposition; RC, rectal cancer.
2.6 Pre-processing of features
Each radiomic feature was standardized using z-score normalization to improve the robustness of the model. Missing information on the clinical characteristics was replaced with the mean value of the corresponding feature. The number of positive and negative samples was balanced using an up-sampling method within the open-source Python package Imbalanced-learn (version 0.9.0) (32).
2.7 Development of radiomic models
Radiomic models based on common machine learning (ML) techniques and DL methods were developed and compared to predict TDs. Three ML models and three DL models were constructed using features from HRT2 images (HRT2-ML and HRT2-DL models), DWI images (DWI-ML and DWI-DL models), and joint HRT2-DWI images (Merged-ML and Merged-DL models).
Integrated models combining radiomic information and clinical characteristics were developed to further improve the predictive performance. Three ML models and three DL models were constructed using clinical characteristics and features from HRT2 images (Clinical-HRT2-ML and Clinical-HRT2-DL models), DWI images (Clinical-DWI-ML and Clinical-DWI-DL models), and joint HRT2-DWI images (Clinical-Merged-ML and Clinical-Merged-DL models). For comparison, a clinical model that analyzed only clinical characteristics was also developed.
ML models used the least absolute shrinkage and selection operator (LASSO) technique (33) to select the optimized subset of features from 221 preprocessed features, followed by a support vector machine (SVM) to construct a prediction model.
DL models used a four-layer multi-layer perceptron (MLP) model, in which 221 preprocessed features were directly input. The feature numbers at each layer were 256, 128, 64, and 2, respectively, and the softmax activation function was used for the final output. The network model was actualized using the open-source deep learning framework PyTorch (34), where the batch size was set to 16, the learning rate was 0.001, and the Adam algorithm was used as the optimizer.
The predictive performance of each model was evaluated using five-fold repeated cross-validation. Each cross-validated split of the data was used to perform feature selection techniques to avoid bias in the estimation of the predictive performance.
2.8 Statistical analyses
Statistical analysis was performed with SPSS software (version 26.0; IBM, Armonk, NY) and R (version 3.5.1; R Foundation, Vienna, Austria). Differences in categorical characteristics between RC patients with and without TDs were compared using Pearson’s chi-squared test and Fisher’s exact test. Continuous variables are expressed as means ± standard deviations. Differences in continuous characteristics between the two groups were compared using the Mann-Whitney U test. For all statistical analyses, P < 0.05 (two-sided test) was considered statistically significant. The predictive performance of the models was evaluated using the area under the receiver operating characteristic (ROC) curve (AUC), over five-fold cross-validation.
3 Results
3.1 Clinical characteristics of patients
The final 500 patients included 315 males and 185 females (mean age, 64.59 ± 10.7 years). According to the pathologically confirmed TD results, the patients were divided into TD+ (n = 133) and TD− (n = 367) subgroups. The baseline characteristics of the patients are summarized in Table 2.
3.2 Performance of the radiomic models
In total, six radiomics models were developed, and the Merge-ML and DWI-DL models both demonstrated comparable performance, with AUCs of 0.69 ± 0.04 and 0.68 ± 0.03, respectively (P<0.05). The other models’ AUCs were lower than those of the aforementioned two models. The AUCs of the HRT2-ML, HRT2-DL, DWI-ML, and Merged-DL models were 0.62 ± 0.02, 0.57 ± 0.06, 0.64 ± 0.08 and 0.59 ± 0.04, respectively.
3.3 Performance of the integrated models combining radiomic information and clinical characteristics
TN staging, tumor length (measured in the sagittal view), tumor index CEA, CRM, and EMVI as assessed in the MRI report were significantly different between the TD+ and TD− groups. These clinical markers were used to establish a clinical model and were introduced into the integrated models. The performance of all models are listed in Table 3. Both the Clinical-ML and Clinical-DL models performed similarly, with AUC values of 0.81 ± 0.04 and 0.81 ± 0.06, respectively. Among the integrated models, the Clinical-DWI-DL model achieved the highest performance, with a diagnostic accuracy of 0.84 ± 0.05, an AUC score of 0.90 ± 0.04, sensitivity of 0.94 ± 0. 03, and specificity of 0.79 ± 0.08 (Figure 3). The Clinical-Merged-ML and Clinical-Merged-DL models achieved similar performances, with AUC scores of 0.83 ± 0.01 and 0.83 ± 0.05, which were both lower than the Clinical-DWI-DL model.
Figure 3 ROC curve of the combined model of DWI and clinical factors. The blue dotted lines represent the ROC curve of five-fold cross-validation (CV), and the red line represents the mean ROC curve of the five-fold CV. ROC, receiver operating characteristic; DWI, diffusion-weighted imaging.
3.4 Robustness of the model on different scanners
We compared the performance of our DWI-DL model on one scanner (GE-Signa HDxt) and all three scanners to assess the robustness of the DL-model on different scanners. We obtained AUC scores of 0.69 ± 0.08 when used with a single scanner and 0.68 ± 0.03 when used with three scanners. We did not detect a statistically significant difference (p = 0.69), even though the model performed better on a single scanner than on three scanners.
4 Discussion
In this study, we developed and validated a variety of models for non-invasive preoperative prediction of TDs in patients with RC, based on radiomic features, clinical factors, and a combination of both. Among all the models, the integrated DL-based model using a combination of DWI radiomic features and clinical characteristics was the most effective and achieved promising predictive performance. This approach can serve as a potential preoperative assessment tool to assist clinicians in preoperative stage evaluation and personalized treatment of patients with RC.
Of the 500 included RC patients, 26% presented with TDs, which is slightly higher than the median incidence (21.3%) of TDs in patients with CRC, as previously reported (2). TDs are an important prognostic factor in CRC, as a significantly worse prognosis has been found in patients with TDs, regardless of the sub-staging of the LNs (35). TDs are also an independent risk factor for liver, lung, and peritoneal metastases (36). Moreover, patients with TDs have a higher risk of LN metastasis and lymphovascular and perineural invasion (37). In this study, the proportion of patients with LN metastases and vascular and nerve invasion in the TD+ group (73.7%, 21.8%, and 30.1%, respectively) was also significantly higher than that in the TD− group (25.9%, 5.2%, and 7.6%, respectively), indicating a possible correlation between TDs and LN metastases, neurovascular invasion, and multi-channel tumor metastases, which are also associated with worse prognosis in patients with CRC.
While TDs cannot be reliably assessed preoperatively using traditional imaging techniques that depend on the naked eye, previous studies have shown that they may be predicted using radiomics, which provides implicit information on tumor heterogeneity far beyond the capability of visual inspection. Radiomic models based on US, CT, and MRI have been established for TD prediction. However, functional MRI (e.g., DWI) provides more information on tumor heterogeneity. Therefore, we established radiomic models based on DWI, which demonstrated higher predictive performance than HRT2-only radiomic models.
The Clinical-ML and Clinical-DL models perform similarly in Table 3. This is because there are only 7 clinical features, which is a relatively small number, and both DL and ML work well with such low-dimensional data. We can also observe that the Clinical-DWI-DL model improves by roughly 9% over the Clinical-DL model, while the Clinical-DWI-ML model barely improves. This could be because deep learning models outperform ML-based models in high-dimensional data situations.
A study on MRI evaluation demonstrated that a joint-modality (HRT2 and DWI) radiomic model achieved higher diagnostic performance than HRT2-only and DWI-only models (38). However, in our study, the joint-modality model did not outperform single-modality models. This result is similar to the findings of Shin et al. who predicted the complete pathological response in RC, and their joint-modality model using features from T2-weighted and DWI images had a classification performance similar to that of the T2-only model (39). We also developed an integrated model that combined radiomic features and clinical characteristics to improve the predictive performance. Among the integrated DL-based models combined with clinical factors, the model utilizing DWI-only radiomic features achieved the highest performance. This may be due to inconsistent baselines and different spatial resolutions between HRT2 and DWI scans of RC, which cannot be reduced by spatial resampling prior to feature extraction.
In addition to investigating the model’s performance, we investigated its robustness due to the complexity of clinical data collection. We chose the model with the best performance (DWI-DL) for this investigation because clinical information is independent of the scanner, allowing us to test the model’s robustness across a range of scanners. We found no significant differences in the radiomics model across machines (AUCs of the single- and multi-scanner models: 0.69 ± 0.08 and 0.68 ± 0.03, respectively, with P = 0.69), indicating good robustness of our radiomics model.
This study has some limitations. First, to our knowledge, this is the largest study to date on TD radiomics research, but it is still not large enough to avoid selection bias that compromises the generalization ability of our models. Second, this retrospective study excluded patients who had received neoadjuvant chemotherapy or radiotherapy before surgery, which introduced a further selection bias. Third, this was a single-center study, and the difference in sample sizes between the TD+ and TD− groups was large. Therefore, further prospective, multicenter studies with larger cohorts are warranted to improve prediction outcomes and define the potential standardization of our models.
5 Conclusions
Our integrated model combining clinical variables (tumor markers and MRI reporting status) and MRI radiomic features in a DL model can non-invasively and preoperatively predict TDs in patients with RC. In particular, the model that used DWI and clinical features showed the highest predictive performance. This model can serve as a potential preoperative assessment tool in clinical practice for more effective tumor staging and risk stratification to provide optimal treatment for patients with RC.
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.
Author contributions
CF, LZ and FZ conceived of the presented idea. MH, JQ, PL, ZY, KS, MW, XW and FL collected the data. LZ, CF, JZ, WL and TS analyzed the data. FL drafted the manuscript. TS and JS improved the quality of English. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by Health Commission of Zhejiang Province (WKJ-ZJ-2210), Science Technology Department of Zhejiang Province (2019C03014), and the National Natural Science Foundation of China (81871403).
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. Delattre JF, Cohen R, Henriques J, Falcoz A, Emile JF, Fratte S, et al. Prognostic value of tumor deposits for disease-free survival in patients with stage III colon cancer: A Post hoc analysis of the IDEA France phase III trial (PRODIGE-GERCOR). J Clin Oncol (2020) 38(15):1702–10. doi: 10.1200/JCO.19.01960
2. Lord AC, D'Souza N, Pucher PH, Moran BJ, Abulafi AM, Wotherspoon A, et al. Significance of extranodal tumour deposits in colorectal cancer: A systematic review and meta-analysis. Eur J Cancer (2017) 82:92–102. doi: 10.1016/j.ejca.2017.05.027
3. Shi M, Zhang H, Yao G, Wu J, Zhu C, Zhang X, et al. The role of tumor deposits in predicting the efficacy of chemotherapy in stage III colon cancer. Front Oncol (2020) 10:586603. doi: 10.3389/fonc.2020.586603
4. Wong-Chong N, Motl J, Hwang G, Nassif GJ Jr, Albert MR, Monson JRT, et al. Impact of tumor deposits on oncologic outcomes in stage III colon cancer. Dis Colon Rectum (2018) 61(9):1043–52. doi: 10.1097/DCR.0000000000001152
5. Kaur H, Choi H, You YN, Rauch GM, Jensen CT, Hou P, et al. MR imaging for preoperative evaluation of primary rectal cancer: Practical considerations. Radiographics (2012) 32(2):389–409. doi: 10.1148/rg.322115122
6. Horvat N, Carlos Tavares Rocha C, Clemente Oliveira B, Petkovska I, Gollub MJ. MRI Of rectal cancer: Tumor staging, imaging techniques, and management. Radiographics (2019) 39(2):367–87. doi: 10.1148/rg.2019180114
7. Gröne J, Loch FN, Taupitz M, Schmidt C, Kreis ME. Accuracy of various lymph node staging criteria in rectal cancer with magnetic resonance imaging. J Gastrointest Surg (2018) 22(1):146–53. doi: 10.1007/s11605-017-3568-x
8. Langman G, Patel A, Bowley DM. Size and distribution of lymph nodes in rectal cancer resection specimens. Dis Colon Rectum (2015) 58(4):406–14. doi: 10.1097/DCR.0000000000000321
9. Glynne-Jones R, Wyrwicz L, Tiret E, Brown G, Rödel C, Cervantes A, et al. Rectal cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol (2017) 28(suppl_4):iv22–40. doi: 10.1093/annonc/mdx224
10. Jin Y, Li M, Zhao Y, Huang C, Liu S, Liu S, et al. Computed tomography-based radiomics for preoperative prediction of tumor deposits in rectal cancer. Front Oncol (2021) 11:710248. doi: 10.3389/fonc.2021.710248
11. 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(7):806–15. doi: 10.6004/jnccn.2020.0032
12. Lord AC, D'Souza N, Shaw A, Rokan Z, Moran B, Abulafi M, et al. MRI-Diagnosed tumor deposits and EMVI status have superior prognostic accuracy to current clinical TNM staging in rectal cancer. Ann Surg (2022) 276(2):334–44. doi: 10.1097/SLA.0000000000004499
13. You J, Yin J. Performances of whole tumor texture analysis based on MRI: Predicting preoperative T stage of rectal carcinomas. Front Oncol (2021) 11:678441. doi: 10.3389/fonc.2021.678441
14. Lu H, Yuan Y, Zhou Z, Ma X, Shen F, Xia Y, et al. Assessment of MRI-based radiomics in preoperative T staging of rectal cancer: Comparison between minimum and maximum delineation methods. BioMed Res Int (2021) 2021:5566885. doi: 10.1155/2021/5566885
15. Huang YQ, Liang CH, He L, Tian J, Liang CS, Chen X, et al. Development and validation of a radiomics nomogram for preoperative prediction of lymph node metastasis in colorectal cancer. J Clin Oncol (2016) 34(18):2157–64. doi: 10.1200/JCO.2015.65.9128
16. Li C, Yin J. Radiomics based on T2-weighted imaging and apparent diffusion coefficient images for preoperative evaluation of lymph node metastasis in rectal cancer patients. Front Oncol (2021) 11:671354. doi: 10.3389/fonc.2021.671354
17. Zhang K, Ren Y, Xu S, Lu W, Xie S, Qu J, et al. A clinical-radiomics model incorporating T2-weighted and diffusion-weighted magnetic resonance images predicts the existence of lymphovascular invasion / perineural invasion in patients with colorectal cancer. Med Phys (2021) 48(9):4872–82. doi: 10.1002/mp.15001
18. Albaradei S, Thafar M, Alsaedi A, Van Neste C, Gojobori T, Essack M, et al. Machine learning and deep learning methods that use omics data for metastasis prediction. Comput Struct Biotechnol J (2021) 19:5008–18. doi: 10.1016/j.csbj.2021.09.001
19. 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(23):7253–62. doi: 10.1158/1078-0432.CCR-17-1038
20. 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(3):833–43. doi: 10.1148/radiol.2018172300
21. Esteva A, Kuprel B, Novoa RA, Ko J, Swetter SM, Blau HM, et al. Dermatologist-level classification of skin cancer with deep neural networks. Nature (2017) 542(7639):115–8. doi: 10.1038/nature21056
22. Ragab DA, Sharkas M, Marshall S, Ren J. Breast cancer detection using deep convolutional neural networks and support vector machines. PeerJ (2019) 7:e6201. doi: 10.7717/peerj.6201
23. Kumar Y, Gupta S, Singla R, Hu YC. A systematic review of artificial intelligence techniques in cancer prediction and diagnosis. Arch Comput Methods Eng (2022) 29(4):2043–70. doi: 10.1007/s11831-021-09648-w
24. Chen LD, Li W, Xian MF, Zheng X, Lin Y, Liu BX, et al. Preoperative prediction of tumour deposits in rectal cancer by an artificial neural network-based US radiomics model. Eur Radiol (2020) 30(4):1969–79. doi: 10.1007/s00330-019-06558-1
25. Yang YS, Feng F, Qiu YJ, Zheng GH, Ge YQ, Wang YT. High-resolution MRI-based radiomics analysis to predict lymph node metastasis and tumor deposits respectively in rectal cancer. Abdom Radiol (NY) (2021) 46(3):873–84. doi: 10.1007/s00261-020-02733-x
26. Beets-Tan RGH, Lambregts DMJ, Maas M, Bipat S, Barbaro B, Curvo-Semedo L, et al. Magnetic resonance imaging for clinical management of rectal cancer: Updated recommendations from the 2016 European society of gastrointestinal and abdominal radiology (ESGAR) consensus meeting. Eur Radiol (2018) 28(4):1465–75. doi: 10.1007/s00330-017-5026-2
27. Kalisz KR, Enzerra MD, Paspulati RM. MRI Evaluation of the response of rectal cancer to neoadjuvant chemoradiation therapy. Radiographics (2019) 39(2):538–56. doi: 10.1148/rg.2019180075
28. 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(2):330–44. doi: 10.1148/radiol.13121361
29. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage (2006) 31(3):1116–28. doi: 10.1016/j.neuroimage.2006.01.015
30. Gu Z, Cheng J, Fu H, Zhou K, Hao H, Zhao Y, et al. CE-net: Context encoder network for 2D medical image segmentation. IEEE Trans Med Imaging (2019) 38(10):2281–92. doi: 10.1109/TMI.2019.2903562
31. 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(21):e104–e7. doi: 10.1158/0008-5472.CAN-17-0339
32. Lemaître G, Nogueira F, Aridas CK. Imbalanced-learn: A Python toolbox to tackle the curse of imbalanced datasets in machine learning. ArXiv (2017). doi: 10.48550/arXiv.1609.06570
33. Sauerbrei W, Royston P, Binder H. Selection of important variables and determination of functional form for continuous predictors in multivariable model building. Stat Med (2007) 26(30):5512–28. doi: 10.1002/sim.3148
34. Paszke A, Gross S, Massa F, Lerer A, Bradbury J, Chanan G, et al. PyTorch: An imperative style, high-performance deep learning library. NeurIPS. Vancouver CANADA: Curran Associates, Inc. (2019).
35. Bai R, Tan Y, Li D, Yang M, Yu L, Yuan Y, et al. Development and validation of a novel prognostic nomogram including tumor deposits could better predict survival for colorectal cancer: A population-based study. Ann Transl Med (2021) 9(8):620. doi: 10.21037/atm-20-4728
36. Cohen R, Shi Q, Meyers J, Jin Z, Svrcek M, Fuchs C, et al. Combining tumor deposits with the number of lymph node metastases to improve the prognostic accuracy in stage III colon cancer: A post hoc analysis of the CALGB/SWOG 80702 phase III study (Alliance)(☆). Ann Oncol (2021) 32(10):1267–75. doi: 10.1016/j.annonc.2021.07.009
37. Nagtegaal ID, Knijn N, Hugen N, Marshall HC, Sugihara K, Tot T, et al. Tumor deposits in colorectal cancer: Improving the value of modern staging-a systematic review and meta-analysis. J Clin Oncol (2017) 35(10):1119–27. doi: 10.1200/JCO.2016.68.9091
38. Park SH, Cho SH, Choi SH, Jang JK, Kim MJ, Kim SH, et al. MRI Assessment of complete response to preoperative chemoradiation therapy for rectal cancer: 2020 guide for practice from the Korean society of abdominal radiology. Korean J Radiol (2020) 21(7):812–28. doi: 10.3348/kjr.2020.0483
Keywords: deep learning, rectal cancer, tumor deposit, magnetic resonance imaging, diffusion-weighted imaging
Citation: Fu C, Shao T, Hou M, Qu J, Li P, Yang Z, Shan K, Wu M, Li W, Wang X, Zhang J, Luo F, Zhou L, Sun J and Zhao F (2023) Preoperative prediction of tumor deposits in rectal cancer with clinical-magnetic resonance deep learning-based radiomic models. Front. Oncol. 13:1078863. doi: 10.3389/fonc.2023.1078863
Received: 24 October 2022; Accepted: 06 February 2023;
Published: 20 February 2023.
Edited by:
Alla Reznik, Lakehead University, CanadaReviewed by:
Manuel Conson, University of Naples Federico II, ItalyJoost Nederend, Catharina Hospital, Netherlands
Copyright © 2023 Fu, Shao, Hou, Qu, Li, Yang, Shan, Wu, Li, Wang, Zhang, Luo, Zhou, Sun and Zhao. 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: Fenhua Zhao, zhfenhua@163.com; Jihong Sun, sunjihong@zju.edu.cn; Long Zhou, zhoulong21@zju.edu.cn
†These authors have contributed equally to this work and share first authorship