- Research article
- Open access
- Published:
Habitat-based radiomics analysis for evaluating immediate response in colorectal cancer lung metastases treated by radiofrequency ablation
Cancer Imaging volume 24, Article number: 44 (2024)
Abstract
Purpose
To create radiomics signatures based on habitat to assess the instant response in lung metastases of colorectal cancer (CRC) after radiofrequency ablation (RFA).
Methods
Between August 2016 and June 2019, we retrospectively included 515 lung metastases in 233 CRC patients who received RFA (412 in the training group and 103 in the test group). Multivariable analysis was performed to identify independent risk factors for developing the clinical model. Tumor and ablation regions of interest (ROI) were split into three spatial habitats through K-means clustering and dilated with 5 mm and 10 mm thicknesses. Radiomics signatures of intratumor, peritumor, and habitat were developed using the features extracted from intraoperative CT data. The performance of these signatures was primarily evaluated using the area under the receiver operating characteristics curve (AUC) via the DeLong test, calibration curves through the Hosmer-Lemeshow test, and decision curve analysis.
Results
A total of 412 out of 515 metastases (80%) achieved complete response. Four clinical variables (cancer antigen 19–9, simultaneous systemic treatment, site of lung metastases, and electrode type) were utilized to construct the clinical model. The Habitat signature was combined with the Peri-5 signature, which achieved a higher AUC than the Peri-10 signature in the test set (0.825 vs. 0.816). The Habitat+Peri-5 signature notably surpassed the clinical and intratumor radiomics signatures (AUC: 0.870 in the test set; both, p < 0.05), displaying improved calibration and clinical practicality.
Conclusions
The habitat-based radiomics signature can offer precise predictions and valuable assistance to physicians in developing personalized treatment strategies.
Introduction
The lung ranks as the second most frequent location for metastases from colorectal cancer (CRC) [1]. Despite the importance of surgery, published reports indicate varying 5-year survival rates of 24 to 56% following surgical removal of CRC lung metastases [2]. European Society for Medical Oncology (ESMO) guidelines for metastatic CRC suggest considering local ablation alongside resection, based on factors like tumor size, number, location, lung tissue loss, comorbidity, and other relevant considerations [3].
Computed tomography (CT)-guided percutaneous radiofrequency ablation (RFA) stands as a viable choice for treating small lung metastases, especially those under 3 cm [4, 5]. However, due to the lack of histological evidence confirming complete ablation and a reported recurrence rate of 32.6% [6], assessing ablation outcomes after RFA remains challenging. The usual assessment of ground-glass opacity (GGO) [7], linked to residual tumor and recurrence [8], can be impacted by intraoperative intra-alveolar hemorrhage (IAH) or atelectasis [4, 9, 10]. The initial 3-month inflammatory response around the ablated lesion complicates early efficacy evaluation [11]. Hence, there’s a need to objectively define and evaluate the immediate ablation outcomes for CRC lung metastases.
In the context of personalized and precise treatments, radiomics can extract high-dimensional quantitative features from medical images. These features encompass data about tumor heterogeneity and the microenvironment [12, 13], enabling a more accurate assessment of traits and treatment response [14, 15]. Peritumoral features have been proposed to enhance radiomic models’ predictive abilities [16,17,18,19]. Unlike previous methods, a novel approach divides tumors into subregions known as habitats, containing voxels with similar attributes and consistent tumor biology [20]. This approach improved intratumoral heterogeneity quantification [21, 22].
To the best of our knowledge, no studies have aimed to develop habitat-based radiomics analysis to predict early RFA efficacy in CRC lung metastases. In this study, we employed a new radiomics approach to identify imaging biomarkers within intratumoral, peritumoral, and sub-regional zones. This enables assessing the immediate response to RFA in CRC lung metastases.
Methods
Patient selection and clinicopathological information
Due to the retrospective nature of this study, patient informed consent was waived. We included 233 consecutive CRC patients with lung metastases who underwent initial RFA between August 2016 and June 2019. Inclusion criteria were: (1) confirmed CRC through histology; (2) lung metastases treated with RFA, ≤ 3 cm in maximum diameter; (3) comprehensive medical records with clinical variables and CT data from procedure and follow-up; (4) technically successful ablation; (5) adequate normal organ function. Exclusion criteria were: (1) receiving other local treatments like radiotherapy or re-ablation; (2) inability to tolerate RFA; (3) concurrent malignant tumors or extrapulmonary metastasis. Patients with multiple nodules were included, analyzing each nodule individually [10, 23, 24]. A cohort of 515 lung metastases in 233 CRC patients who underwent RFA was enrolled and randomly divided into training and test cohorts at a 4:1 ratio (Fig. 1). Supplementary data 1 provided detailed RFA equipment and procedure information.
The clinical variables including age, gender, serum tumor markers (carcinoembryonic antigen (CEA) and cancer antigen 19–9 (CA19–9)), lymphadenopathy at diagnosis, concomitant systemic treatment, and primary tumor location were collected within 1 week before RFA. Radiological data from intraoperative CT scans included pulmonary metastases size, location, proximity to vital structures like the heart or major blood vessels (> 3 mm in diameter), distance to the pleura or diaphragm (within 1 cm), electrode type, and complications such as IAH or pneumothorax.
CT examination protocol and local efficacy assessment
Pre- and immediate post-ablation CT scans were conducted using the United Imaging uCT 760 (United Imaging Medical Technology Inc.) and Philips Brilliance 64 slice (Philips Medical Systems Inc.) machines. Settings were: 200 mA, 120 kVp, 0.5 s/round, with 1 mm or 1.5 mm section thickness. The images were reconstructed using iterative reconstruction, and the resulting CT data was stored in the. DICOM format.
A contrast-enhanced chest CT was conducted as the baseline 1 month after ablation [25], followed by additional scans every 2 to 3 months. The treatment’s local efficacy was assessed by two experienced radiologists who were unaware of clinical data. Evaluation was grounded in chest-enhanced CT scans performed at least 6 months after RFA, adhering to the modified response evaluation criteria in solid tumors (mRECIST) criteria [11, 26]. Should there be differences in interpretation between the radiologists, consultation with a senior expert boasting over 20 years of experience was pursued. Complete response (CR) was determined by the presence of cavity, fibrosis, or nodule without enhancement. In contrast, the presence of irregular nodular, scattered, or eccentric patterns of peripheral enhancement within 1 cm of the ablation area in two consecutive CT scans denoted a non-complete response (non-CR).
Workflow of radiomics analysis
The radiomics analysis was executed through a series of steps: image segmentation, feature extraction, feature selection, signature construction, and evaluation (Fig. 2).
Advanced image processing and mask segmentation
To enhance the robustness of medical image analysis, preprocessing techniques were applied. The CT images were uniformed to a common resolution of 1 mm × 1 mm × 1 mm by the B-spline interpolation algorithm, and then the window width was adjusted within the range of - 1200 Hu to 600 Hu and the intensity was scaled within the range of 0 ~ 255.
Recently, numerous studies have demonstrated that the improved methods based on U-Net performed well in the segmentation of pulmonary nodules, which were trained on the Image Database Resource Initiative (IDRI) created by the US Institutes of Health based on the Lung Image Database Consortium (LIDC) [27,28,29,30,31]. We also have trained a 3D UNet model based on the open-source dataset to segment the target lesions and immediate ablation regions from pre- and postoperative CT images, with a Dice coefficient of 83.04% [32]. These masks were subsequently verified by two junior radiologists (HZH and HC, 8 years of specialized chest imaging) and the necessary adjustments have been made to guarantee accuracy and repeatability using the ITK-SNAP (version 3.8.0, http://www.itksnap.org). If they had disagreements, it would be determined in consultation with the senior expert (WTL, 30 years of specialized chest imaging).
Peritumoral region dilation
The regions of interest (ROI) were expanded using the mask padding toolkit provided by the Onekey AI platform. We evaluated the impact of different peritumoral sizes on model predictability by applying dilation intervals of 5 mm. Any ROIs extending beyond the lungs or overlapping with the heart, major blood vessels, or diaphragm were manually adjusted.
Habitat generation
Local features, including local entropy and energy values, were extracted from each voxel within VOI. These feature vectors represented diverse aspects of voxel properties. A 77-dimensional feature vector (Supplementary data 2) was generated for each block using a 3 × 3 × 3 non-overlapping moving window. The Calinski-Harabasz (CH) value selection method [33] was used to determine the optimal number of clusters. Subsequently, the K-means method was employed to cluster sub-regions for each sample.
Feature extraction and selection
Handcrafted features extracted using the Pyradiomics tool (version 3.0.1) were categorized into three groups: geometry, intensity, and texture (Supplementary data 2), following the guidelines of the imaging biomarker standardization initiative (IBSI). Unsupervised clustering yielded varying physical meanings for the habitat extracted from the subregions with identical centers. To mitigate this, mean feature values were computed.
Robustness was assessed through test-retest and inter-rater analyses, with the intraclass correlation coefficient (ICC) set at a threshold of ≥0.85. Nonetheless, ICC was not suitable for assessing the unsupervised habitat signature. All features were standardized using Z-scores to maintain a normal distribution, followed by a t-test to retain radiomic features with p-value < 0.05. Pearson’s correlation coefficient was computed to pinpoint highly consistent features with a coefficient > 0.9. The greedy recursive deletion strategy was then employed to filter out the highly redundant features. To curb overfitting, the minimum redundancy maximum relevance (mRMR) algorithm selected the top 8 features for each modality. The final features set was determined using the least absolute shrinkage and selection operator (LASSO) regression. LASSO adjusted parameter λ to assign zero regression coefficients to irrelevant features. Optimal λ value selection involved 10-fold cross-validation with minimum criteria, aiming for the lowest mean square error (MSE).
Signature construction
Several radiomics signatures were formulated based on distinct regions: intratumor and ablated area (Intra), intratumor with expanded tumor and ablated regions (5 mm and 10 mm, Peri-X), intratumor and ablated area subregions (Habitat), and intratumor combined with peritumoral regions. Additionally, the optimal peritumoral region was integrated with the tumor microenvironment habitat, termed Habitat + Peri-X. The Clinical signature was created from independent risk factors identified via multivariate logistic analysis. Commonly used machine learning models, including logistic regression (LR), support vector machine (SVM), K-nearest neighbor (KNN), random forest (RF), extremely randomized trees (ExtRa Trees), eXtreme gradient boosting (XGBoost), light gradient boosting machine (LightGBM), and multi-layer perceptron (MLP), were employed for model construction. Optimal hyperparameters for each model were determined using five-fold cross-validation and the Grid-search algorithm.
Performance evaluation
The performance of various signatures was verified using an independent test dataset, generating receiver operating characteristic (ROC) curves to calculate the corresponding area under the curve (AUC). The Delong test was used to compare predictive performance differences among the models [34]. Additionally, accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were computed. The Youden index determined the optimal cut-off value maximizing the sum of sensitivity and specificity [35]. Calibration curves were plotted to assess calibration accuracy, alongside the Hosmer-Lemeshow (HL) test [36] (A significant test statistic implies that the model does not calibrate perfectly.). Moreover, decision curve analysis (DCA) gauged the clinical utility of predictive signatures [37].
Statistical analysis
Statistical analyses were performed using IBM SPSS (version 26.0). Continuous variables were presented as mean ± standard deviation (SD) and compared using the Man-Whitney U test. Categorical variables were expressed as counts with percentages and compared using the Chi-square or Fisher test. Variables with a P-value < 0.05 in univariate regression analysis were included in multivariable analysis. Variables with a P-value < 0.05 in multivariable analysis were considered independent predictors linked to immediate response. All statistical tests were two-sided with a significance level set at P < 0.05.
Results
Baseline characteristics
The training dataset consisted of 412 lesions (330 CR and 82 non-CR) selected via random division, while 103 lesions (82 CR and 21 non-CR) contributed the independent test dataset (Table 1).
Establishment of clinical models
Univariate logistic regression analysis (Table 2) revealed significant factors for identifying ablated lesions: CEA, CA19–9, concomitant systemic treatment, lung metastases location, and electrode type (P < 0.05). Multivariate regression analysis indicated that CA19–9 (odds ratio [OR] = 1.002, P < 0.001), concomitant systemic treatment (OR = 0.916, P = 0.042), lung metastases location (OR = 1.032, P = 0.019), and electrode type (OR = 0.778, P = 0.004) were independent factors influencing ablation effect, forming the basis for clinical models construction.
Sub-region cluster and feature selection
Optimal CH value emerged when tumors were clustered into three sub-regions in the entire cohort (Fig. 3). To enhance the feature set, features from both pre- and post-ablation images were fused, resulting in 3668 features. Features were extracted via Pyradiomics (http://pyradiomics.readthedocs.io).
Radiomics features with non-zero coefficients were selected using the LASSO method with the best lambda (Supplementary data 3). For habitat-based radiomics signature, lambda of 0.0095 yielded the best, selecting 71 features. These comprised 20 habitat features (8 pre- and 12 post-RFA) and 51 peritumoral features (29 pre- and 22 post-RFA).
Performance and comparison of signatures
Based on the analysis of prediction performance (Supplementary data 3), the Intra radiomics signature employed RF, the Peri-X radiomics signature utilized LightGBM, the Habitat radiomics signature was developed using ExtRa Trees, while the Habitat+Peri-5 radiomics and clinical signatures were constructed with XGBoost.
Summary of predictive performance for clinical and radiomics signatures was presented in Table 3 and Fig. 4. The clinical signature achieved AUC values of 0.827 (training) and 0.697 (test), hinting at potentially limited generalizability to unseen data. The Peri-5 signature outperformed Peri-10 in test AUC, prompting its combination with the Habitat signature. The Habitat+Peri-5 signature excelled in both training (AUC 0.972) and test (AUC 0.846) sets (Fig. 4 a and d). Moreover, calibration curves of the Habitat + Peri-5 signature exhibited robust concurrence between observed and projected probabilities (Fig. 4 b and e). The HL test yielded insignificant deviations (training cohort: P = 0.212; test cohort: P = 0.283), signifying conformity. DCA portrayed the enhanced clinical benefit of the fusion signature over other signatures (Fig. 4 c and f).
Discussion
This study systematically explored the tumor microenvironment through a comprehensive evaluation of intra- and peritumoral regions in pre- and post-ablation CT scans. As a result, we devised predictive signatures for early ablation efficacy in CRC lung metastases treated with RFA, leveraging sub-regional radiomics features. These signatures provided effective tools for tailoring treatment strategies in CRC patients with lung metastases.
Radiomics analysis enables the extraction and characterization of a broad range of information reflecting underlying biological diversity in a cost-effective manner [12]. A radiomics signature, composed of multiple features, serves as a robust prognostic biomarker that could complement clinical data [38]. The overestimation of the completely necrotic region as GGO on CT [39] indicates the need to ablate peritumoral lung parenchyma within 5 to 10 mm for complete ablation [9, 40]. Previous studies have highlighted the utility of combining intra- and peritumoral radiomics features for improved treatment response prediction [16, 19]. These findings underscore the importance of peritumoral features [16,17,18, 41]. Therefore, our study aimed to explore the predictive performance of peritumoral radiomics features by assessing the impact of peritumoral region sizes. The results demonstrated that the radiomics signature derived from a 5 mm dilated distance outside both the tumor and ablated area yielded the highest prediction performance compared to the intra- and peritumoral 10 mm regions. This signature exhibited high AUC values and low over-fitting, emphasizing the significance of peritumoral features in radiomics and the influence of peritumoral size on prediction performance.
The sub-region cluster analysis conducted in this study highlighted the importance of sub-region analysis in capturing CRC lung metastases’ tumor heterogeneity. By optimizing the number of clusters, we identified three spatially distinct habitats through K-means clustering. These habitats within the ablation area represented various characteristics: the inner region possibly indicated necrotic areas, the middle region potentially signified effusion areas, and the edge region likely represented congested areas. These findings were consistent with the histopathological results, demonstrating an inner-to-outer transition [39]. As anticipated, the habitat signature, derived from radiomics features within these unique subregions, contributed valuable information in assessing early effects after RFA. Furthermore, by combining habitat and peritumoral signatures, we achieved higher AUC values (training cohort: 0.972, test cohort: 0.870). This outcome validated the superiority of the Habitat+Peri-5 signature over the Intra radiomics signature (P = 0.039, Fig. 5).
Furthermore, multivariate regression analyses identified several clinical variables, including CA19–9 levels, lung metastases location, concomitant systemic treatment, and electrode type. Notably, lung metastases location, particularly in the lower lobes, emerged as an independent risk factor, possibly due to the influence of respiratory movements on accurate probe positioning. Conversely, concomitant systemic treatment and the use of expandable electrodes were protective factors against treatment failure, evident from OR values below 1. Importantly, the Habitat + Peri-5 signature demonstrated a significantly improved AUC value (P = 0.020, Fig. 5) compared to the clinical model (training cohort: 0.827, test cohort: 0.697), underscoring enhanced precision and clinical utility, as supported by the HL test and DCA results.
A study by Markich et al. [10] examined clinical, radiological, technical, and radiomics features examined before and after RFA to evaluate local control in 48 CRC patients with 119 lung metastases. However, their reliance on CT scans taken 48 hours post-ablation hindered real-time procedure assessment, delaying timely interventions. To address this limitation, our recent study [42] integrated relative radiomics features from pre- and immediate post-RFA CT scans with clinical and radiological variables from 479 lung metastases in 198 CRC patients. This aimed to establish a novel multimodal data fusion model for evaluating immediate RFA efficacy. Liu et al. [43] explored intratumor density heterogeneity changes following microwave ablation (MWA) of pulmonary tumors, utilizing radiomics-based CT features for prognostic value in predicting treatment response. Additionally, Zhu et al. [44] developed intra- and peritumoral radiomics models based on post-operative CT images to predict early MWA efficacy in malignant lung tumors, validating the predictive ability. However, neither study explored the optimal peritumoral region nor assessed the impact of peritumoral region sizes. Also, the inclusion of primary and metastatic lung tumors in both studies raised concerns about disease heterogeneity.
This study bears several limitations that warrant acknowledgment. Firstly, it was a single-center retrospective study with a limited sample size. Therefore, a more extensive, multi-center prospective study is essential to validate the generalizability of our signature and ascertain its utility in clinical settings. Secondly, the spatially distinct habitats identified through K-means clustering could not be pathologically confirmed due to technical challenges and ethical considerations. However, pursuing such a correlation might entail unnecessary surgeries.
Conclusion
In summary, this study introduces a pioneering signature that combines habitat and peritumoral radiomics to access immediate response and predict outcomes of RFA in CRC patients with lung metastases. The habitat-based radiomics signature holds the potential for advancing precision medicine and shaping treatment strategies.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- AUC:
-
Area under ROC curve
- ACC:
-
Accuracy
- CA199:
-
Cancer antigen 199
- CEA:
-
Carcinoembryonic antigen
- CH:
-
Calinski–Harabasz
- CI:
-
Confidence interval
- CR:
-
Complete response
- CRC:
-
Colorectal cancer
- CT:
-
Computed tomography
- DCA:
-
Decision curve analysis
- ESMO:
-
European Society for Medical Oncology
- ExtRa:
-
Trees Extremely randomized trees
- GGO:
-
Ground-glass opacity
- HL:
-
Hosmer-Lemeshow
- IAH:
-
Intra-alveolar hemorrhage
- IBSI:
-
Imaging biomarker standardization initiative
- ICC:
-
Intraclass correlation coefficient
- IDRI:
-
Image Database Resource Initiative
- KNN:
-
K-nearest neighbor
- LASSO:
-
Least absolute shrinkage and selection operator
- LIDC:
-
Lung Image Database Consortium
- LightGBM:
-
Light gradient boosting machine
- LR:
-
Logistic regression
- MLP:
-
Multi-layer perceptron
- mRECIST:
-
Modified response evaluation criteria in solid tumors
- mRMR:
-
Minimum redundancy maximum relevance
- MSE:
-
Mean square error
- MWA:
-
Microwave ablation
- NPV:
-
Negative predictive value
- OR:
-
Odds ratio
- PPV:
-
Positive predictive value
- RF:
-
Random forest
- RFA:
-
Radiofrequency ablation
- ROC:
-
Receiver operating characteristic
- ROI:
-
Region of interest
- SD:
-
Standard deviation
- SVM:
-
Support vector machine
- VOI:
-
Volumes of interest
- XGBoost:
-
eXtreme gradient boosting
References
Mitry E, Guiu B, Cosconea S, Jooste V, Faivre J, Bouvier AM. Epidemiology, management and prognosis of colorectal cancer with lung metastases: a 30-year population-based study. Gut. 2010;59:1383–8. https://doi.org/10.1136/gut.2010.211557.
Pfannschmidt J, Dienemann H, Hoffmann H. Surgical resection of pulmonary metastases from colorectal cancer: a systematic review of published series. Ann Thorac Surg. 2007;84:324–38. https://doi.org/10.1016/j.athoracsur.2007.02.093.
Cervantes A, Adam R, Roselló S, Arnold D, Normanno N, Taïeb J, et al. Metastatic colorectal cancer: ESMO clinical practice guideline for diagnosis, treatment and follow-up. Ann oncol: offic j Europ Soc Med Oncol ESMO. 2023;34:10–32. https://doi.org/10.1016/j.annonc.2022.10.003.
de Baère T, Aupérin A, Deschamps F, Chevallier P, Gaubert Y, Boige V, et al. Radiofrequency ablation is a valid treatment option for lung metastases: experience in 566 patients with 1037 metastases. Ann Oncol. 2015;26:987–91. https://doi.org/10.1093/annonc/mdv037.
Ibrahim T, Tselikas L, Yazbeck C, Kattan J. Systemic versus local therapies for colorectal Cancer pulmonary metastasis: what to choose and when? J Gastrointest Cancer. 2016;47:223–31. https://doi.org/10.1007/s12029-016-9818-4.
Hiyoshi Y, Miyamoto Y, Kiyozumi Y, Sawayama H, Eto K, Nagai Y, et al. CT-guided percutaneous radiofrequency ablation for lung metastases from colorectal cancer. Int J Clin Oncol. 2019;24:288–95. https://doi.org/10.1007/s10147-018-1357-5.
Abtin FG, Eradat J, Gutierrez AJ, Lee C, Fishbein MC, Suh RD. Radiofrequency ablation of lung tumors: imaging features of the postablation zone. Radiographics. 2012;32:947–69. https://doi.org/10.1148/rg.324105181.
de Baère T, Palussière J, Aupérin A, Hakime A, Abdel-Rehim M, Kind M, et al. Midterm local efficacy and survival after radiofrequency ablation of lung tumors with minimum follow-up of 1 year: prospective evaluation. Radiol. 2006;240:587–96. https://doi.org/10.1148/radiol.2402050807.
Yang Q, Qi H, Zhang R, Wan C, Song Z, Zhang L, et al. Risk factors for local progression after percutaneous radiofrequency ablation of lung tumors: evaluation based on a review of 147 tumors. J Vasc Interv Radiol. 2017;28:481–9. https://doi.org/10.1016/j.jvir.2016.11.042.
Markich R, Palussière J, Catena V, Cazayus M, Fonck M, Bechade D, et al. Radiomics complements clinical, radiological, and technical features to assess local control of colorectal cancer lung metastases treated with radiofrequency ablation. Eur Radiol. 2021; https://doi.org/10.1007/s00330-021-07998-4.
Abtin FG, Eradat J, Gutierrez AJ, Lee C, Fishbein MC, Suh RD. Radiofrequency ablation of lung tumors: imaging features of the postablation zone. Radiograph. 2012;32:947–69. https://doi.org/10.1148/rg.324105181.
Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006. https://doi.org/10.1038/ncomms5006.
Rios Velazquez E, Parmar C, Liu Y, Coroller TP, Cruz G, Stringfield O, et al. Somatic mutations drive distinct imaging phenotypes in lung Cancer. Cancer Res. 2017;77:3922–30. https://doi.org/10.1158/0008-5472.Can-17-0122.
Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures. They Are Data Radiol. 2016;278:563–77. https://doi.org/10.1148/radiol.2015151169.
Limkin EJ, Sun R, Dercle L, Zacharaki EI, Robert C, Reuzé S, et al. Promises and challenges for the implementation of computational medical imaging (radiomics) in oncology. Ann Oncol. 2017;28:1191–206. https://doi.org/10.1093/annonc/mdx034.
Braman N, Prasanna P, Whitney J, Singh S, Beig N, Etesami M, et al. Association of Peritumoral Radiomics with Tumor Biology and Pathologic Response to preoperative targeted therapy for HER2 (ERBB2)-positive breast Cancer. JAMA Netw Open. 2019;2:e192561. https://doi.org/10.1001/jamanetworkopen.2019.2561.
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. Radiol. 2019;290:783–92. https://doi.org/10.1148/radiol.2018180910.
Shi J, Dong Y, Jiang W, Qin F, Wang X, Cui L, et al. MRI-based peritumoral radiomics analysis for preoperative prediction of lymph node metastasis in early-stage cervical cancer: a multi-center study. Magn Reson Imaging. 2022;88:1–8. https://doi.org/10.1016/j.mri.2021.12.008.
Braman NM, Etesami M, Prasanna P, Dubchuk C, Gilmore H, Tiwari P, et al. Intratumoral and peritumoral radiomics for the pretreatment prediction of pathological complete response to neoadjuvant chemotherapy based on breast DCE-MRI. Breast Cancer Res. 2017;19:57. https://doi.org/10.1186/s13058-017-0846-1.
O'Connor JP, Rose CJ, Waterton JC, Carano RA, Parker GJ, Jackson A. Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome. Clin cancer res: an official j Am Assoc Cancer Res. 2015;21:249–57. https://doi.org/10.1158/1078-0432.Ccr-14-0990.
Park JE, Kim HS, Kim N, Park SY, Kim YH, Kim JH. Spatiotemporal heterogeneity in multiparametric physiologic MRI is associated with patient outcomes in IDH-wildtype glioblastoma. Clin cancer res: an official j Am Assoc Cancer Res. 2021;27:237–45. https://doi.org/10.1158/1078-0432.Ccr-20-2156.
Kim J, Ryu SY, Lee SH, Lee HY, Park H. Clustering approach to identify intratumour heterogeneity combining FDG PET and diffusion-weighted MRI in lung adenocarcinoma. Eur Radiol. 2019;29:468–75. https://doi.org/10.1007/s00330-018-5590-0.
Staal FCR, Taghavi M, van der Reijd DJ, Gomez FM, Imani F, Klompenhouwer EG, et al. Predicting local tumour progression after ablation for colorectal liver metastases: CT-based radiomics of the ablation zone. Eur J Radiol. 2021;141:109773. https://doi.org/10.1016/j.ejrad.2021.109773.
Taghavi M, Staal F, Gomez Munoz F, Imani F, Meek DB, Simões R, et al. CT-based Radiomics analysis before thermal ablation to predict local tumor progression for colorectal liver metastases. Cardiovasc Intervent Radiol. 2021;44:913–20. https://doi.org/10.1007/s00270-020-02735-8.
Lencioni R, Crocetti L, Cioni R, Suh R, Glenn D, Regge D, et al. Response to radiofrequency ablation of pulmonary tumours: a prospective, intention-to-treat, multicentre clinical trial (the RAPTURE study). Lancet Oncol. 2008;9:621–8. https://doi.org/10.1016/s1470-2045(08)70155-4.
Wei Z, Ye X, Yang X, Zheng A, Huang G, Li W, et al. Microwave ablation in combination with chemotherapy for the treatment of advanced non-small cell lung cancer. Cardiovasc Intervent Radiol. 2015;38:135–42. https://doi.org/10.1007/s00270-014-0895-0.
Roy R, Mazumdar S, Chowdhury AS. ADGAN: Attribute-driven generative adversarial network for synthesis and multiclass classification of pulmonary nodules. IEEE Trans Neural Netw Learn Syst. 2024;35:2484–95. https://doi.org/10.1109/tnnls.2022.3190331.
Yu H, Li J, Zhang L, Cao Y, Yu X, Sun J. Design of lung nodules segmentation and recognition algorithm based on deep learning. BMC Bioinformat. 2021;22:314. https://doi.org/10.1186/s12859-021-04234-0.
Jalali Y, Fateh M, Rezvani M, Abolghasemi V, Anisi MH. ResBCDU-net: a deep learning framework for lung CT image segmentation. Sensors (Basel). 2021:21. https://doi.org/10.3390/s21010268.
Dong X, Xu S, Liu Y, Wang A, Saripan MI, Li L, et al. Multi-view secondary input collaborative deep learning for lung nodule 3D segmentation. Cancer Imag. 2020;20:53. https://doi.org/10.1186/s40644-020-00331-0.
Usman M, Lee BD, Byon SS, Kim SH, Lee BI, Shin YG. Volumetric lung nodule segmentation using adaptive ROI with multi-view residual learning. Sci Rep. 2020;10:12839. https://doi.org/10.1038/s41598-020-69817-y.
Ni Y, Xie Z, Zheng D, Yang Y, Wang W. Two-stage multitask U-net construction for pulmonary nodule segmentation and malignancy risk prediction. Quantit imag med surg. 2022; 12: 292–309 https://doi.org/10.21037/qims-21-19.
Liu Y, Li Z, Xiong H, Gao X, Wu J, Wu S. Understanding and enhancement of internal clustering validation measures. IEEE transact cybernet. 2013;43:982–94. https://doi.org/10.1109/tsmcb.2012.2220543.
DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44:837–45.
Hu T, Wang S, Huang L, Wang J, Shi D, Li Y, et al. A clinical-radiomics nomogram for the preoperative prediction of lung metastasis in colorectal cancer patients with indeterminate pulmonary nodules. Eur Radiol. 2019;29:439–49. https://doi.org/10.1007/s00330-018-5539-3.
Kramer AA, Zimmerman JE. Assessing the calibration of mortality benchmarks in critical care: the Hosmer-Lemeshow test revisited. Crit Care Med. 2007;35:2052–6. https://doi.org/10.1097/01.Ccm.0000275267.64078.B0.
Vickers AJ, Cronin AM, Elkin EB, Gonen M. Extensions to decision curve analysis, a novel method for evaluating diagnostic tests, prediction models and molecular markers. BMC Med Inform Decis Mak. 2008;8:53. https://doi.org/10.1186/1472-6947-8-53.
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:2157–64. https://doi.org/10.1200/jco.2015.65.9128.
Yamamoto A, Nakamura K, Matsuoka T, Toyoshima M, Okuma T, Oyama Y, et al. Radiofrequency ablation in a porcine lung model: correlation between CT and histopathologic findings. AJR Am J Roentgenol. 2005;185:1299–306. https://doi.org/10.2214/ajr.04.0968.
Yan P, Tong AN, Nie XL, Ma MG. Assessment of safety margin after microwave ablation of stage I NSCLC with three-dimensional reconstruction technique using CT imaging. BMC Med Imaging. 2021;21:96. https://doi.org/10.1186/s12880-021-00626-z.
Ding J, Chen S, Serrano Sosa M, Cattell R, Lei L, Sun J, et al. Optimizing the Peritumoral region size in Radiomics analysis for sentinel lymph node status prediction in breast Cancer. Acad Radiol. 2022;29(Suppl 1):S223–s28. https://doi.org/10.1016/j.acra.2020.10.015.
Huang H, Zheng D, Chen H, Chen C, Wang Y, Xu L, et al. A CT-based radiomics approach to predict immediate response of radiofrequency ablation in colorectal cancer lung metastases. Front. Oncol. 2023:13:1107026. https://doi.org/10.3389/fonc.2023.1107026.
Liu B, Li C, Sun X, Zhou W, Sun J, Liu H, et al. Assessment and prognostic value of immediate changes in post-ablation Intratumor density heterogeneity of pulmonary tumors via Radiomics-based computed tomography features. Front Oncol. 2021;11:615174. https://doi.org/10.3389/fonc.2021.615174.
Zhu F, Yang C, Xia Y, Wang J, Zou J, Zhao L, et al. CT-based radiomics models may predict the early efficacy of microwave ablation in malignant lung tumors. Cancer Imaging. 2023;23:60. https://doi.org/10.1186/s40644-023-00571-w.
Acknowledgements
We thank the Department of Medical Imaging, Mental Health Center, Shanghai Jiao Tong University School of Medicine, Laboratory for Medical Imaging Informatics, Shanghai Institute of Technical Physics, Chinese Academy of Science, and OnekeyAI platform, for their assistance with this work.
Authorship principles
Haozhe Huang, Hong Chen, Dezhong Zheng, and Chao Chen contributed equally to this work. Wentao Li and Yuanyuan Yang were both corresponding authors.
Funding
This work was supported by the National Natural Science Foundation of China (No. 82272095) and Shanghai Anticancer Association EYAS PROJECT (SACA-CY23B03).
Author information
Authors and Affiliations
Contributions
Haozhe Huang, Wentao Li and Yuanyuan Yang conceived of the study. Haozhe Huang, Hong Chen, Dezhong Zheng, and Chao Chen participated in the design of the study and the initial synthesis of results. Yaohui Wang, Xinhong He, and Lichao Xu performed the ablation operations. Haozhe Huang, Hong Chen, Chao Chen, and Ying Wang collected the clinical and radiological data. Haozhe Huang and Hong Chen were major contributors in writing the manuscript. Dezhong Zheng and Chao Chen built models and made statistical analyses. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of Fudan University Shanghai Cancer Center (Approval Number: 1612167–18).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Huang, H., Chen, H., Zheng, D. et al. Habitat-based radiomics analysis for evaluating immediate response in colorectal cancer lung metastases treated by radiofrequency ablation. Cancer Imaging 24, 44 (2024). https://doi.org/10.1186/s40644-024-00692-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40644-024-00692-w