Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 06 October 2022
Sec. Gastrointestinal Cancers: Colorectal Cancer

Radiomics features based on internal and marginal areas of the tumor for the preoperative prediction of microsatellite instability status in colorectal cancer

Yi Ma&#x;Yi Ma1†Changsong Lin&#x;Changsong Lin2†Song LiuSong Liu1Ying WeiYing Wei3Changfeng JiChangfeng Ji1Feng ShiFeng Shi3Fan Lin*Fan Lin4*Zhengyang Zhou*Zhengyang Zhou1*
  • 1Department of Radiology, Nanjing Drum Tower Hospital Clinical College of Nanjing Medical University, Nanjing, China
  • 2Department of Bioinformatics, Nanjing Medical University, Nanjing, China
  • 3Department of Research and Development, Shanghai United Imaging Intelligence Co., Ltd., Shanghai, China
  • 4Department of Cell Biology, Nanjing Medical University, Nanjing, China

Objectives: To explore whether the preoperative CT radiomics can predict the status of microsatellite instability (MSI) in colorectal cancer (CRC) patients and identify the region with the most stable and high-efficiency radiomics features.

Methods: This retrospective study involved 230 CRC patients with preoperative computed tomography scans and available MSI status between December 2019 and October 2021. Image segmentation and radiomic feature extraction were performed as follows. First, slices with the maximum tumor area (region of interest, ROI) were manually contoured. Subsequently, each ROI was shrunk inward by 1, 2, and 3 mm, respectively, where the remaining ROIs were considered as the internal region of the tumor (named as IROI1, IROI2, and IROI3), and the shrunk regions were considered as marginal regions of the tumor (named as MROI1, MROI2, and MROI3). Finally, radiomics features were extracted from each of the ROI. The intraclass correlation coefficient and least absolute shrinkage and selection operator method were used to choose the most reliable and relevant features of MSI status. Clinical, radiomics, and combined clinical radiomics models have been established. Calibration curve and decision curve analyses (DCA) were generated to explore the correction effect and assess the clinical applicability of the above models, respectively.

Results: In the testing cohort, the radiomics model based on IROI3 yielded the highest average area under the curve (AUC) value of 0.908, compared with the remaining radiomics models. Additionally, hypertension and N stage were considered as clinically independent factors of MSI status. The combined clinical radiomics model achieved excellent diagnostic efficacy (AUC: 0.928; sensitivity: 0.840; specificity: 0.867) in the testing cohort, as well as favorable calibration and clinical utility by calibration curve and DCA analyses.

Conclusions: The IROI3 model, which is based on a 3-mm shrink in the largest areas of the tumor, could noninvasively reflect the heterogeneity and genetic instability within the tumor. This suggests that it is an important biomarker for the preoperative prediction of MSI status. The model can extract more robust and effective radiomics features, which lays a foundation for the radiomics study of hollow organs, such as in CRC.

Introduction

Colorectal cancer (CRC) ranks the second leading cause of cancer-related mortality worldwide, with approximately 700,000 deaths each year (1, 2). Microsatellite instability (MSI), reflects the spontaneous loss or gain of nucleotides from repetitive DNA tracts, which is present in about 15% of CRC cases (3, 4). This gene replication error is usually repaired by the DNA mismatch repair (MMR) system, which helps maintain genomic stability and reduces spontaneous mutations. MSI status could be subdivided into microsatellite instability-high (MSI-H), microsatellite instability-low (MSI-L) and microsatellite stable (MSS) according to the proportion of loci with MSI (5).

MSI status is a vital predictive factor for the screening, prognosis, and therapeutic decisions of CRC patients. First, identifying MSI status is helpful for screening patients for Lynch syndrome, the most common form of hereditary CRC (6). Second, MSI is an indicator of a good prognosis for stage II CRC patients. Compared with MSS CRC, the overall and disease-free survival of stage II CRC patients with MSI is significantly prolonged (7). Third, MSI is crucial for developing treatment strategies in patients with CRC. Some randomized controlled trials have shown that adjuvant chemotherapy can improve the overall survival rate of patients with MSS (8). But MSI patients can benefit from immunotherapy (9, 10) due to their multiple tumor mutation sites (11) and extensive immunogenicity (12).

However, methods to evaluate MSI status, including polymerase chain reaction (PCR) and immunohistochemistry (IHC), are invasive and costly (13, 14). Moreover, for inoperable patients, a small biopsy sample may not be sufficient to clarify the MSI status due to tumor heterogeneity (15, 16). Therefore, the development of a noninvasive, repeatable, and effective MSI prediction method before surgery is of great significance in treatment decision-making for patients with CRC.

Radiomics is an emerging technology for acquiring high-dimensional image data about tumor phenotypes and microenvironments that cannot be detected by the naked eye (17, 18). As a non-invasive and reproducible radiological biomarker, radiomics analysis has shown great potential in tumor staging (19), prognosis evaluation (20, 21), and KRAS status prediction (22) for CRC. Recent reports have shown that peritumoral and intratumoral radiomic features provide value for evaluating tumor biological behavior (2326).

Considering the particularity of intestinal tumors, their edges are easily affected by air and feces in the intestinal cavity, resulting in artifacts (27, 28). We considered whether this would further affect the stability of the extracted radiomics features. Therefore, we retrospectively collected the clinical and MSI status information of patients with CRC. Through image segmentation and feature extraction from the internal and marginal regions of the tumors, we constructed a clinical model, seven radiomics models, and a visual nomogram to predict the preoperative MSI status. This study aimed to evaluate whether analyzing the imaging characteristics from different regions within the tumor is more helpful in predicting MSI status and identifying the region with the most stable and high-efficiency radiomics features.

Materials and methods

Patients

This retrospective study was approved by the institutional ethical board at our hospital and exempted from informed consent. Data of 323 patients with CRC proved by surgical pathology were continuously collected from December 2019 to October 2021. Patient inclusion details were as follows: (a) patient underwent the abdominal enhanced computed tomography (CT) examination before surgery, (b) CRC confirmed by postoperative pathology, and (c) available MSI status tested by IHC. Patient exclusion details were as follows: (a) insufficient image quality to identify the tumor delineation with motion or metal artifacts (n = 16); (b) the interval between CT examination and operation exceeded 2 weeks (n = 9); (c) the maximum cross-sectional short diameter of the lesion was< 1 cm on CT images (n = 47); and (d) any anticancer therapy was performed before CT imaging (n = 21). The subject inclusion and exclusion criteria are presented in Figure 1. Finally, 50 patients with MSI and 180 with MSS CRC were enrolled in our study.

FIGURE 1
www.frontiersin.org

Figure 1 Flowchart of patient selection and grouping process. MSS Microsatellite stability, MSI Microsatellite instability.

Clinical data, including sex, age, comorbidities, tumor location, TNM stage, histologic grade, carcinoembryonic antigen (CEA) level, carbohydrate antigen (CA) 125 level, and CA199 level, were abstracted from the medical records. These findings were unanimously confirmed by both clinicians.

MSI status assessment

IHC staining of postoperative pathological tissues of CRC was first performed by a standard streptavidin-biotin-peroxidase procedure (29), then MSI status was determined by evaluating the results of the four MMR proteins (MLH1, MSH2, MSH6, and PMS2). Any of the four MMR protein negative expressions were identified as MSI, while all four MMR proteins positive expressions were identified as MSS (30).

CT scan protocol

Patients underwent fasting for more than 4 h before the examination; during the scan, the patients were instructed to hold their breath, and each one received a flat scan first, followed by a three-phase enhanced scan. The scanning range was from the diaphragmatic apex to the bilateral suprapubic level. All studies were completed on the same CT scanner and received the same examination protocol. Parameters of the CT scan protocol are listed in Table S1 of the Supplementary Material.

Image segmentation and feature extraction

The images of patients, which were selected from the picture archiving and communication system with a 5-mm venous layer thickness, were downloaded in DICOM format and uploaded to a research platform, the uAI Research Portal (Shanghai United Imaging Intelligence, Co., Ltd.). The flowchart contained image segmentation, feature extraction and selection, model building, and evaluation (Figure 2). All lesions were manually segmented by a senior abdominal radiologist (reader 1 [M.Y.] with 9 years of experience) blinded to MSI status. The slice with the maximum tumor area was selected, including bleeding and necrosis within the tumor, avoiding perienteric fat, vascular, air, and feces, and labeled as regions of interest (ROI). For the patients with multiple lesions, the largest lesion was selected for ROI delineation according to the endoscopic findings. We applied a set of morphological operations to the tumor region to analyze the information of the different areas inside the tumor. Since the outside of the tumor is mostly gas which may not provide useful information, we opted to use the morphological shrinkage operation (31) to automatically shrink the tumor boundary inward by 1, 2, and 3 mm. In this way, we define the remaining tumor regions as interior areas of the tumor (regarded as IROI1, IROI2, and IROI3), and assign the shrunken ring regions as the margin areas of the tumor (called MROI1, MROI2, and MROI3). Thus, the original tumor ROI is the sum of the marginal and corresponding internal regions (Figure 3).

FIGURE 2
www.frontiersin.org

Figure 2 Workflow of MSI status prediction of colorectal cancer patients including image segmentation, feature extraction and selection, model building and evaluation.

FIGURE 3
www.frontiersin.org

Figure 3 The tumors were segmented on venous phase, three groups of different marginal and internal regions were formed by retracting the tumor margin by 1, 2, and 3 mm on the representative section. The red areas represent the maximum slice of the tumor (ROI). The blue areas represent the rings with the thickness of 1-3mm in the marginal regions of the tumor (MROI), respectively. The hot-pink areas represent the remaining internal regions of the tumor (IROI) after removing the rings, respectively. ROI, regions of interest. MROI, the marginal areas of tumor. IROI, the interior areas of tumor.

Before Radiomics feature extraction, we first resampled all images into 1.5×1.5×1.5 mm3 to reduce any heterogeneity in acquisition voxel size. Second, we used an abdominal imaging window (window width (WW): 400, window level (WL): 40) for CT images normalization by min-max normalization method and adjusted the intensity range of each image to 0-255. Radiomics features were subsequently extracted from the widely used PyRadiomics (version 3.0.1) package (32), which contains seven feature categories: first-order statistics, shape, gray-level co-occurrence matrix (GLCM), gray-level dependence matrix (GLDM), gray-level size zone matrix (GLSZM), gray-level run length matrix (GLRLM), and neighboring gray-level tone difference matrix (NGTDM). Fourteen image filters were implemented with SimpleITK (version 2.1.0) package in Python were adopted (33), and Radiomics features were obtained from the filtered images as well as the original images. The image filters include additive Gaussian noise, binomial blur image, box mean, box sigma image, curvature flow, Laplacian sharpening, discrete Gaussian, mean, normalize, recursive Gaussian, shot noise, speckle noise, Laplacian of Gaussian and wavelet. Ultimately, 2,600 radiomics features were extracted from each annotation of the ROI, MROI1, MROI2, MROI3, IROI1, IROI2, and IROI3.

Feature selection

Two months after image annotation, the above procedures were repeated by reader 1 and another abdominal radiologist (reader 2 [C.J.F.] with six years of experience) using 30 (34) randomly selected images. First, the ROIs of 30 patients were manually contoured, and the radiomic features of MROI1, MROI2, MROI3, IROI1, IROI2, and IROI3 were automatically obtained for each patient. The inter-/intra- delineator reproducibility was evaluated by the inter-/intra-class correlation coefficients (ICCs) with a pairwise correlation method, and features with both inter-class correlation coefficient and intra-class correlation coefficient greater than 0.75 are considered to have good or excellent reliability and used for subsequent feature selections. Z-score standardized normalization method was utilized to guarantee the comparability of different features by rescale the features of different levels into the same level. After the normalization, we used the least absolute shrinkage and selection operator (LASSO) algorithm (35), an approach that calculates the regression coefficients and successively shrinks them to avoid overfitting, to select the most helpful features in distinguishing between MSI and MSS. The corresponding coefficients were evaluated to calculate the Rad-score (RS) for each patient. RS was calculated using the following equation:

Radscore =i=1nCi × Xi+ b

where Xi is the ith feature, Ci is the coefficient of the ith feature from the LASSO, b is the intercept of the LASSO regression algorithm, and n is the number of selected features.

Model and nomogram construction

In this study, a stratified 5-fold cross-validation strategy was used to randomly but equally divided the patients into five partitions to ensure that the same percentage of each class (i.e., MSI/MSS) is preserved in each partition. One partition was taken each time without repetition as the testing cohort, and the remaining four partitions as the training cohort. The above steps were repeated five times to obtain five different sets of training-testing cohorts. The average value was obtained for a more stable and accurate sample evaluation (Figure 1).

The clinical parameters with P< 0.1 from the univariate analysis were enrolled into the multivariate regression analysis with P< 0.05 for identifying independent predictors of MSI status and constructing the clinical model (26, 36). Logistic regression analysis was used to develop seven radiomics models as well as a combined clinical radiomics model. The combined model was built with incorporating the RS of the optimal radiomics model (26) and the clinical independent predictors selected by clinical model, and it was presented in the form of an individualized nomogram.

Statistical analysis

Clinical data were analyzed using the Mann–Whitney U test for continuous variables and χ2 test for categorical variables. A receiver operating characteristic (ROC) curve was used to evaluate the predictive effectiveness of each model. The areas under the ROC curve (AUC) of the radiomics models were quantified and compared. Calibration plots and decision curve analysis (DCA) were created to explore the correction effect and assess the clinical applicability of clinical, radiomics, and combined clinical radiomics models, respectively. The net reclassification index (NRI) was used to assess and compare the predictive power of the three models: NRI > 0 is a positive improvement, indicating that the predictive efficiency of the model has improved; NRI< 0 is a negative improvement, implying that the predictive efficiency of the model has decreased. All statistical analyses were conducted with the R software (version 3.5.2; http://www.Rproject.org) and IBM SPSS Statistics for Windows, version 26 (IBM Corp., Armonk, N.Y., USA). All statistical tests were two-sided, and statistical significance was set at P< 0.05.

Results

Patient characteristic

There were 230 CRC patients enrolled for analysis, consisting of 139 males (60.4%) and 91 females (39.6%) with an average age of 62.7 ± 11.8 years (range, 27–93 years); There were 203 patients with BMI records. The average BMI of the patients was 23.45 kg/m2 (range, 14.06 – 41.29 kg/m2). 180 patients (78.3%) were confirmed as having MSS and 50 (21.7%) were confirmed as having MSI. A comparison of the clinical characteristics between the two groups is summarized in Table 1. Univariate analysis revealed significant differences in patient gender, hypertension, tumor location, and N stage (all P< 0.05).

TABLE 1
www.frontiersin.org

Table 1 Characteristics of patients [median (Q1, Q3) or no. (%)].

Clinical model building

Logistic regression analysis identified hypertension [β= -0.971, OR = 0.379 (95% confidence intervals (CI), 0.175–0.822), P = 0.014] and N stage [β= -1.338, OR = 0.262 (95% CI, 0.123–0.561), P = 0.001] as independent factors of MSI status. A clinical model incorporating the above predictors had been developed and it produced moderate performance with an AUC of 0.695 (95% CI, 0.523–0.867) in the testing cohort. (Supplementary Table S2).

Radiomics feature selection and model analysis

The features with low reproducibility (ICC values< 0.75) were removed as described earlier, so the number of features extracted from ROI MROI1, MROI2, MROI3, IROI1, IROI2, and IROI3 were reduced to 2069, 1260, 1485, 1600, 2029, 2003, and 1828, respectively. After Z-score normalization, the most representative radiomics features were screened by LASSO to build the logistic regression model. The number of selected features and the efficiency of the seven different radiomics models are listed in Table 2. This indicated that the established 7 radiomics models could perform well (all AUC ≥ 0.820) to predict the MSI status of CRC preoperatively. In the testing cohort, the internal regions of the tumor labeled as IROI1 [AUC: 0.858 (95% CI, 0.728–0.989)], IROI2 [AUC: 0.861 (95% CI, 0.751–0.969)], and IROI3 [AUC: 0.908 (95% CI, 0.821–0.991)] produced higher AUC values than the three corresponding marginal areas of tumor regarded as MROI1 [AUC: 0.820 (95% CI, 0.674–0.965)], MROI2 [AUC: 0.840 (95% CI, 0.714–0.961)], and MROI3 [AUC: 0.846 (95% CI, 0.722–0.969)]. In addition, by analyzing the internal regions of the tumor labeled as IROI1, IROI2, and IROI3, an improvement in predictive performance was observed. Moreover, the radiomics model based on IROI3 yielded the best predictive performance (accuracy: 0.813; sensitivity: 0.800; specificity: 0.817). The ROC analysis to assess the performance of different radiomics models is shown in Figure 4. Ultimately, IROI3 was chosen as the final radiomics model because of its best prediction ability compared to the remaining models. (Figure 4, Table 2). The features used to build the IROI3 model are shown in Figure S1 in the Supplementary Material.

TABLE 2
www.frontiersin.org

Table 2 Predictive performance of seven radiomics models in training and testing cohorts.

FIGURE 4
www.frontiersin.org

Figure 4 The receiver operating characteristic (ROC) curves of seven radiomics models in training cohort (A) and testing cohort (B). IROI3 model [area under the curve (AUC) = 0.960 and 0.908 in the training and testing cohort, respectively] achieved better performance than the other radiomics models.

Clinical application

The rad-score of IROI3 and selected clinically independent factors were used to develop a combined model based on logistic regression, which was presented as a quantitative nomogram (Figure 5). In the testing cohort, the combined model (AUC: 0.928 [95% CI, 0.860–0.991)] had a better performance than the clinical or radiomics models. NRIs further indicated significant improvements in the combined model compared to the clinical model (NRI: 0.490, P< 0.001), with no significant improvements in the combined model compared to the radiomics model (NRI: 0.090, P = 0.097) (Table 3). The calibration curve showed that the radiomics and combined models had a better agreement between observation and prediction to evaluate MSI status than the clinical model. The decision curves of the combined model gained the highest net benefit compared to the other two models at ranges of 12–70% for the radiomics model and 0–100% for the clinical model (Figure 6).

FIGURE 5
www.frontiersin.org

Figure 5 An individualized nomogram for preoperative predicting the status of MSI in colorectal cancer patients. The Nomogram was built based on the N stage, hypertension and the Rad-score of the optimal radiomics model.

TABLE 3
www.frontiersin.org

Table 3 Pairwise comparisons of AUCs of the clinical model, radiomics model, and combined model.

FIGURE 6
www.frontiersin.org

Figure 6 Calibration curves of the clinical model, radiomics model and combined model (A). The diagonal dotted line indicated perfect prediction, and the solid lines indicated the prediction performance of the three models respectively. The closer to the dotted line, the better prediction of the model. The calibration curve showed the radiomics and combined models with favorable performance for predicting MSI status, which was better than clinical model. Decision curve analysis (DCA) of the three model (B). The x-axis represented the threshold probability. The y-axis was the net benefit. The higher curve at any range of threshold probability was the best prediction to maximize the net benefit. The DCA indicated that the combined model provides a better clinical utility than the other two models.

Discussion

In this retrospective study, radiomics features were extracted from different regions of the tumor for the prediction of MSI status to explore the effects of air and feces around the tumors on the model. We constructed three groups of different marginal and internal models by retracting the tumor margin by 1, 2, and 3 mm on the representative section. Seven radiomics models were developed to evaluate the MSI status. The clinical features were added to build a visual nomogram. Our results showed that the model based on the 3-mm adduction had the highest AUC value to predict the MSI state. The nomogram showed outstanding prediction with an AUC value of 0.928 in the testing cohort.

This study included 12 clinical indicators. The incidence of MSI was 21.74% (50/230), and it mostly occurred in the right colon at a rate of 64% (32/50), which is consistent with previous literature (37, 38). We also found that N stage and hypertension were clinically independent predictors and closely related to MSI status. Some studies (7, 39) have shown that patients with MSI CRC have a better prognosis, which may be associated with the lower rate of lymph node metastasis in the current study. The incidence rate of hypertension was lower in patients with MSI; these intriguing findings may help to clarify the status of MSI before operation and save costs. Although a recent study (40) based on 100 patients showed that there was no significant relationship between hypertension and MSI status, this conclusion needs to be further validated on a larger data.

Consistent with most relevant studies in recent years (4144), we adopted the IHC method to determine MSI status. Although PCR is the gold standard for diagnosing MSI status, the operation process is complex and the cost is high. IHC has a high correlation with the detection results of PCR. It provides a cost-effective, sensitive (92.3%), and extremely specific (100%) method for screening for DNA mismatch repair defects (29).

CRC is surrounded by irrelevant information, such as air and feces in the intestinal cavity, and the radiomics features extracted from adjacent areas may be affected. Accordingly, in our study, the maximum cross-section of the lesion was delineated, and the adduction technique was used to analyze radiomics features from different intratumoral regions. The short diameters of the included tumors were > 1 cm, which not only optimize the clinical efficacy (45), but also avoid image annihilation during processing. The tumor margin was retracted by 1–3 mm respectively, and seven radiomics models were generated. By comparing the AUC values of the models, 23 quantitative radiomic features were selected to calculate the RS of the optimal radiomic model. Among them, 19 texture features were obtained from five categories of texture features (i.e., GLCM, GLDM, GLSZM, GLRLM, and NGTDM). These features were used as a measure of the grayscale non-uniformity of the images to reflect the inherent heterogeneity of the tumor. This finding was consistent with a previous study in which the commonest radiomics features to predict MSI status were the texture features (46).

Our research found that the model based on a 3 mm adduction was determined as the optimal radiomics model. The combined clinical radiomics model showed excellent prediction with an AUC of 0.928 in the testing cohort, which was higher than that reported in previous studies (4144). Furthermore, the internal regions of the tumor produced a higher AUC value than the corresponding marginal areas of the tumor. The closer to the tumor center, the greater the improvement in the predictive performance.

The excellent ability of our model to predict MSI status may be explained as follows. Pathologically, MSI CRC tends to manifest as mixed morphological characteristics, such as mucinous, glandular, and solid components (47), which causes tumor heterogeneity and made it possible to enable the analysis of preoperative MSI status with radiomics. Moreover, the invasive ability of cancer cells differs between the center and edge of the tumor. A study by Zhao et al. (48) found that compared with cancer cells at the edge of the tumor, cancer cells in the center of the tumor were easier to metastasize and spread. Furthermore, the marginal areas of the tumor may be affected by adjacent air or feces. The generated artifacts may make the extracted features unstable, thus affecting the model performance. Compared to the radiomics model, we found that the performance of the combined model was not significantly improved after adding clinical features in the testing cohort (P = 0.097). This result also implied the independent value of the radiomics model we created in the preoperative prediction of MSI status in CRC patients.

In order to promote clinical practice, we built a radiomics nomogram based on the characteristics of radiology and clinical features to realize the non-invasive and individualized prediction of the preoperative MSI status of CRC patients by clinicians. Moreover, considering the influence of intestinal contents and feces around the colorectal cancer, we measured adduction of the primary lesions in 2D plane to determine the ROI and build models to predict the preoperative microsatellite status. In this study, the ROI region was an objective and quantitative analysis and processing based on the location of the primary lesion. The simplicity and repeatability of the prediction models were further verified in the testing cohort. Therefore, the results of this study provide a new method to determine ROI for future imaging radiomics research on CRC.

Nevertheless, this study still had several limitations. Firstly, it was a single-center study with limited samples. Hence, further verification is necessary by conducting an external and multicenter study. Second, there might have been a selection bias owing to the retrospective study. Third, all images were obtained using the same CT scanner, which might have affected the generalizability of our results. Fourth, for small tumors with a 3 mm adduction, the remaining internal area was too small, affecting the modeling performance, although we chose tumors with a short diameter greater than 1 cm to avoid image annihilation. Further exploration is needed to develop an equal-proportion adduction processing software according to the tumor size.

Conclusion

Three different groups of marginal and internal models to predict MSI status, were constructed by adducting the tumor edge by 1–3 mm. Our study confirmed that a model based on the 3-mm adduction can noninvasively reflect tumor heterogeneity and genetic instability. It is an important biomarker for the preoperative prediction of MSI status in CRC patients. The model can extract more stable and effective radiomics features, which lays a foundation for the radiomics study of hollow organs, such as in CRC.

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 Medical Ethics Committee of the Nanjing Drum Tower Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

Data curation, study design, manuscript writing, and manuscript approval were performed by YM and CL, they were accountable for all aspects of the work. CT data analysis and manuscript approval were performed by SL and CJ. Statistical analysis and manuscript approval were performed by YW and FS. FL and ZZ conceptualized and designed the study. All authors contributed to the article and approved the submitted version.

Conflict of interest

YW and FS are employees of Shanghai United Imaging Intelligence Co., Ltd. The company has no role in performing the surveillances and interpreting the data.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

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

References

1. Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Piñeros M, et al. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer (2019) 144(8):1941–53. doi: 10.1002/ijc.31937

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Siegel RL, Miller KD, Fedewa SA, Ahnen DJ, Meester RGS, Barzi A, et al. Colorectal cancer statistics, 2017. CA Cancer J Clin (2017) 67(3):177–93. doi: 10.3322/caac.21395

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Dekker E, Tanis PJ, Vleugels JLA, Kasi PM, Wallace MB. Colorectal cancer. Lancet (2019) 394(10207):1467–80. doi: 10.1016/S0140-6736(19)32319-0

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Chang L, Chang M, Chang HM, Chang F. Microsatellite instability: A predictive biomarker for cancer immunotherapy. Appl Immunohistochem Mol Morphol (2018) 26(2):e15–21. doi: 10.1097/PAI.0000000000000575

PubMed Abstract | CrossRef Full Text | Google Scholar

5. De' Angelis GL, Bottarelli L, Azzoni C, De' Angelis N, Leandro G, Di Mario F, et al. Microsatellite instability in colorectal cancer. Acta BioMed (2018) 89(9-S):97–101. doi: 10.23750/abm.v89i9-S.7960

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Vasen HF. Clinical description of the lynch syndrome [hereditary nonpolyposis colorectal cancer (HNPCC)]. Fam Cancer (2005) 4(3):219–25. doi: 10.1007/s10689-004-3906-5

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Merok MA, Ahlquist T, Røyrvik EC, Tufteland KF, Hektoen M, Sjo OH, et al. Microsatellite instability has a positive prognostic impact on stage II colorectal cancer after complete resection: results from a large, consecutive Norwegian series. Ann Oncol (2013) 24(5):1274–82. doi: 10.1093/annonc/mds614

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ribic CM, Sargent DJ, Moore MJ, Thibodeau SN, French AJ, Goldberg RM, et al. Tumor microsatellite-instability status as a predictor of benefit from fluorouracil-based adjuvant chemotherapy for colon cancer. N Engl J Med (2003) 349(3):247–57. doi: 10.1056/NEJMoa022289

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD, et al. PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med (2015) 372(26):2509–20. doi: 10.1056/NEJMoa1500596

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chalabi M, Fanchi LF, Dijkstra KK, Van den Berg JG, Aalbers AG, Sikorska K, et al. Neoadjuvant immunotherapy leads to pathological responses in MMR-proficient and MMR-deficient early-stage colon cancers. Nat Med (2020) 26(4):566–76. doi: 10.1038/s41591-020-0805-8

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Hegde M, Ferber M, Mao R, Samowitz W, Ganguly A, Working Group of the American College of Medical Genetics and Genomics (ACMG) Laboratory Quality Assurance Committee. ACMG technical standards and guidelines for genetic testing for inherited colorectal cancer (Lynch syndrome, familial adenomatous polyposis, and MYH-associated polyposis). Genet Med (2014) 16(1):101–16. doi: 10.1038/gim.2013.166

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Fischer F, Baerenfaller K, Jiricny J. 5-fluorouracil is efficiently removed from DNA by the base excision and mismatch repair systems. Gastroenterology (2007) 133(6):1858–68. doi: 10.1053/j.gastro.2007.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sepulveda AR, Hamilton SR, Allegra CJ, Grody W, Cushman-Vokoun AM, Funkhouser WK, et al. Molecular biomarkers for the evaluation of colorectal cancer guideline from the American society for clinical pathology, college of American pathologists, association for molecular pathology, and American society of clinical oncology. J Clin Oncol (2017) 35(13):1453–86. doi: 10.1200/JCO.2016.71.9807

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Umar A, Boland CR, Terdiman JP, Syngal S, Srivastava S. Revised Bethesda guidelines for hereditary nonpolyposis colorectal cancer (Lynch syndrome) and microsatellite instability. J Natl Cancer Inst (2004) 96(4):261–8. doi: 10.1093/jnci/djh034

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Sacher AG, Dahlberg SE, Heng J, Mach S, Jn PA, Oxnard GR. Association between younger age and targetable genomic alterations and prognosis in non–Small-Cell lung cancer. JAMA Oncol (2016) 2(3):313–20. doi: 10.1001/jamaoncol.2015.4482

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Itakura H, Achrol AS, Mitchell LA, Loya JJ, Liu T, Westbroek EM, et al. Magnetic resonance image features identify glioblastoma phenotypic subtypes with distinct molecular pathway activities. Sci Transl Med (2015) 7(303):303ra138–303ra138. doi: 10.1126/scitranslmed.aaa7582

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Gardin I, Grégoire V, Gibon D, Kirisli H, Pasquier D, Thariat J, et al. Radiomics: Principles and radiotherapy applications. Crit Rev Oncol Hematol (2019) 138:44–50. doi: 10.1016/j.critrevonc.2019.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

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(3):365–75. doi: 10.1007/s00259-017-3855-5

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Cui Y, Liu H, Ren J, Du X, Xin L, Li D, et al. Development and validation of a MRI-based radiomics signature for prediction of KRAS mutation in rectal cancer. Eur Radiol (2020) 30(4):1948–58. doi: 10.1007/s00330-019-06572-3

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Prasanna P, Patel J, Partovi S, Madabhushi A, Tiwari P. Erratum to: Radiomic features from the peritumoral brain parenchyma on treatment-naïve multi-parametric MR imaging predict long versus short-term survival in glioblastoma multiforme: Preliminary findings. Eur Radiol (2017) 27(10):4198–9. doi: 10.1007/s00330-017-4815-y

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Beig N, Khorrami M, Alilou M, Prasanna P, Braman N, Orooji M, et al. Perinodular and intranodular radiomic features on lung CT images distinguish adenocarcinomas from granulomas. Radiology (2019) 290(3):783–92. doi: 10.1148/radiol.2018180910

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wang X, Zhao X, Li Q, Xia W, Peng Z, Zhang R, et al. Can peritumoral radiomics increase the efficiency of the prediction for lymph node metastasis in clinical stage T1 lung adenocarcinoma on CT? Eur Radiol (2019) 29(11):6049–58. doi: 10.1007/s00330-019-06084-0

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chong H, Gong Y, Pan X, Liu A, Chen L, Yang C, et al. Peritumoral dilation radiomics of gadoxetate disodium-enhanced MRI excellently predicts early recurrence of hepatocellular carcinoma without macrovascular invasion after hepatectomy. J Hepatocell Carcinoma (2021) 8:545–63. doi: 10.2147/JHC.S309570

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lee JS, Kim SH, Kim JJ, Kim BS, Choi GM, Kim DR. Gastrointestinal air motion artifact which can be mistaken for active gastrointestinal bleeding in multidetector computed tomography: Phantom and clinical study. J Comput Assist Tomogr (2020) 44(1):145–52. doi: 10.1097/RCT.0000000000000972

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Winklhofer S, Lambert JW, Wang ZJ, Sun Y, Gould RG, Zagoria RJ, et al. Reduction of peristalsis-related gastrointestinal streak artifacts with dual-energy CT: a patient and phantom study. Abdom Radiol (NY). (2016) 41(8):1456–65. doi: 10.1007/s00261-016-0702-2

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Lindor NM, Burgart LJ, Leontovich O, Goldberg RM, Cunningham JM, Sargent DJ, et al. Immunohistochemistry versus microsatellite instability testing in phenotyping colorectal tumors. J Clin Oncol (2002) 20(4):1043–8. doi: 10.1200/JCO.2002.20.4.1043

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Gelsomino F, Barbolini M, Spallanzani A, Pugliese G, Cascinu S. The evolving role of microsatellite instability in colorectal cancer: A review. Cancer Treat Rev (2016) 51:19–26. doi: 10.1016/j.ctrv.2016.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Li J, Zhang C, Wei J, Zheng P, Zhang H, Xie Y, et al. Intratumoral and peritumoral radiomics of contrast-enhanced CT for prediction of disease-free survival and chemotherapy response in stage II/III gastric cancer. Front Oncol (2020) 10:552270. doi: 10.3389/fonc.2020.552270

PubMed Abstract | CrossRef Full Text | Google Scholar

32. 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–7. doi: 10.1158/0008-5472.CAN-17-0339

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Yaniv Z, Lowekamp BC, Johnson HJ, Beare R. SimpleITK image-analysis notebooks: a collaborative environment for education and reproducible research. J Digit Imaging (2018) 31(3):290–303. doi: 10.1007/s10278-017-0037-8

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Liu D, Zhang W, Hu F, Yu P, Zhang X, Yin H, et al. A bounding box-based radiomics model for detecting occult peritoneal metastasis in advanced gastric cancer: A multicenter study. Front Oncol (2021) 11:777760. doi: 10.3389/fonc.2021.777760

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Tibshirani R. Regression shrinkage and selection via the lasso. J R Stat Soc Ser B Methodol (1996) 58(1):267–88.

Google Scholar

36. Chong HH, Yang L, Sheng RF, Yu YL, Wu DJ, Rao SX, et al. Multi-scale and multi-parametric radiomics of gadoxetate disodium-enhanced MRI predicts microvascular invasion and outcome in patients with solitary hepatocellular carcinoma ≤ 5 cm. Eur Radiol (2021) 31(7):4824–38. doi: 10.1007/s00330-020-07601-2

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Boland CR, Goel A. Microsatellite instability in colorectal cancer. Gastroenterology (2010) 138(6):2073–2087.e3. doi: 10.1053/j.gastro.2009.12.064

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Shin US, Cho SS, Moon SM, Park SH, Jee SH, Jung EJ, et al. Is microsatellite instability really a good prognostic factor of colorectal cancer? Ann Coloproctol (2014) 30(1):28–34. doi: 10.3393/ac.2014.30.1.28

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Taieb J, Shi Q, Pederson L, Alberts S, Wolmark N, Van Cutsem E, et al. Prognosis of microsatellite instability and/or mismatch repair deficiency stage III colon cancer patients after disease recurrence following adjuvant treatment: results of an ACCENT pooled analysis of seven studies. Ann Oncol (2019) 30(9):1466–71. doi: 10.1093/annonc/mdz208

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Atef N, Alieldin N, Sherif G, Loay I, Mostafa Mahmoud A, Mohamed G. Microsatellite instability and life style factors in sporadic colorectal cancer. Asian Pac J Cancer Prev (2020) 21(5):1471–80. doi: 10.31557/APJCP.2020.21.5.1471

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Cao Y, Zhang G, Zhang J, Yang Y, Ren J, Yan X, et al. Predicting microsatellite instability status in colorectal cancer based on triphasic enhanced computed tomography radiomics signatures: A multicenter study. Front Oncol (2021) 11:687771. doi: 10.3389/fonc.2021.687771

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Li Z, Zhong Q, Zhang L, Wang M, Xiao W, Cui F, et al. Computed tomography-based radiomics model to preoperatively predict microsatellite instability status in colorectal cancer: A multicenter study. Front Oncol (2021) 11:666786. doi: 10.3389/fonc.2021.666786

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ying M, Pan J, Lu G, Zhou S, Fu J, Wang Q, et al. Development and validation of a radiomics-based nomogram for the preoperative prediction of microsatellite instability in colorectal cancer. BMC Cancer (2022) 22(1):524. doi: 10.1186/s12885-022-09584-3

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wu J, Zhang Q, Zhao Y, Liu Y, Chen A, Li X, et al. Radiomics analysis of iodine-based material decomposition images with dual-energy computed tomography imaging for preoperatively predicting microsatellite instability status in colorectal cancer. Front Oncol (2019) 9:1250. doi: 10.3389/fonc.2019.01250

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Pickhardt PJ, Graffy PM, Weigman B, Deiss-Yehiely N, Hassan C, Weiss JM. Diagnostic performance of multitarget stool DNA and CT colonography for noninvasive colorectal cancer screening. Radiology (2020) 297(1):120–9. doi: 10.1148/radiol.2020201018

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Fan S, Li X, Cui X, Zheng L, Ren X, Ma W, et al. Computed tomography-based radiomic features could potentially predict microsatellite instability status in stage II colorectal cancer: A preliminary study. Acad Radiol (2019) 26(12):1633–40. doi: 10.1016/j.acra.2019.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

47. De Smedt L, Lemahieu J, Palmans S, Govaere O, Tousseyn T, Van Cutsem E, et al. Microsatellite instable vs stable colon carcinomas: analysis of tumour heterogeneity, inflammation and angiogenesis. Br J Cancer (2015) 113(3):500–9. doi: 10.1038/bjc.2015.213

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Zhao Y, Fu X, Lopez JI, Rowan A, Au L, Fendler A, et al. Selection of metastasis competent subclones in the tumour interior. Nat Ecol Evol (2021) 5(7):1033–45. doi: 10.1038/s41559-021-01456-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: microsatellite instability, radiomics, colorectal cancer, internal and marginal, computed tomography

Citation: Ma Y, Lin C, Liu S, Wei Y, Ji C, Shi F, Lin F and Zhou Z (2022) Radiomics features based on internal and marginal areas of the tumor for the preoperative prediction of microsatellite instability status in colorectal cancer. Front. Oncol. 12:1020349. doi: 10.3389/fonc.2022.1020349

Received: 16 August 2022; Accepted: 20 September 2022;
Published: 06 October 2022.

Edited by:

Chunxiao Guo, University of Minnesota, United States

Reviewed by:

Jinjin Cao, Massachusetts General Hospital and Harvard Medical School, United States
Yuhao Zeng, Cleveland Clinic, United States
Cassiano Crespo Santiago, University of Minnesota Twin Cities, United States

Copyright © 2022 Ma, Lin, Liu, Wei, Ji, Shi, Lin and Zhou. 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: Fan Lin, linfan@njmu.edu.cn; Zhengyang Zhou, zyzhou@nju.edu.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.