Homological radiomics analysis for prognostic prediction in lung cancer patients

Radiomics score derived from T1-w/T2-w ratio image can predict motor symptom progression in Parkinson's disease

OBJECTIVES

To clarify the association between a radiomics score (Rad-score) derived from T1-weighted signal intensity to T2-weighted signal intensity (T1-w/T2-w) ratio images and the progression of motor symptoms in Parkinson's disease (PD).

MATERIALS AND METHODS

This retrospective study included patients with PD enrolled in the Parkinson's Progression Markers Initiative. The Movement Disorders Society-Unified Parkinson's Disease Rating Scale Part III score ≥ 33 and/or Hoehn and Yahr stage ≥ 3 indicated motor function decline. The Rad-score was constructed using radiomics features extracted from T1-w/T2-w ratio images. The Kaplan-Meier analysis and Cox regression analyses were used to assess the time differences in motor function decline between the high and low Rad-score groups.

RESULTS

A total of 171 patients with PD were divided into training (n = 101, mean age at baseline, 61.6 ± 9.3 years) and testing (n = 70, mean age at baseline, 61.6 ± 10 years). The patients in the high Rad-score group had a shorter time to motor function decline than those in the low Rad-score group in the training dataset (log-rank test, p < 0.001) and testing dataset (log-rank test, p < 0.001). The multivariate Cox regression using the Rad-score and clinical factors revealed a significant association between the Rad-score and motor function decline in the training dataset (HR = 2.368, 95%CI:1.423-3.943, p < 0.001) and testing dataset (HR = 2.931, 95%CI:1.472-5.837, p = 0.002).

CONCLUSION

Rad-scores based on radiomics features derived from T1-w/T2-w ratio images were associated with the progression of motor symptoms in PD.

CLINICAL RELEVANCE STATEMENT

The radiomics score derived from the T1-weighted/T2-weighted ratio images offers a predictive tool for assessing the progression of motor symptom in patients with PD.

KEY POINTS

Radiomics score derived from T1-weighted/T2-weighted ratio images is correlated with the motor symptoms of Parkinson's disease. A high radiomics score correlated with faster motor function decline in patients with Parkinson's disease. The proposed radiomics score offers predictive insight into the progression of motor symptoms of Parkinson's disease.

Radiomics-based discriminant analysis of principal components to stratify the treatment response of lung metastases following stereotactic body radiation therapy

PURPOSE

Discriminant analysis of principal components (DAPC) was introduced to describe the clusters of genetically related individuals focusing on the variation between the groups of individuals. Borrowing this approach, we evaluated the potential of DAPC for the evaluation of clusters in terms of treatment response to SBRT of lung lesions using radiomics analysis on pre-treatment CT images.

MATERIALS AND METHODS

80 pulmonary metastases from 56 patients treated with SBRT were analyzed. Treatment response was stratified as complete, incomplete and null responses. For each lesion, 107 radiomics features were extracted using the PyRadiomics software. The concordance correlation coefficients (CCC) between the radiomics features obtained by two segmentations were calculated. DAPC analysis was performed to infer the structure of "radiomically" related lesions for treatment response assessment. The DAPC was performed using the "adegenet" package for the R software.

RESULTS

The overall mean CCC was 0.97 ± 0.14. The analysis yields 14 dimensions in order to explain 95 % of the variance. DAPC was able to group the 80 lesions into the 3 different clusters based on treatment response depending on the radiomics features characteristics. The first Linear Discriminant achieved the best discrimination of individuals into the three pre-defined groups. The greater radiomics loadings who contributed the most to the treatment response differentiation were associated with the "sphericity", "correlation" and "maximal correlation coefficient" features.

CONCLUSION

This study demonstrates that a DAPC analysis based on radiomics features obtained from pretreatment CT is able to provide a reliable stratification of complete, incomplete or null response of lung metastases following SBRT.

Topology-based radiomic features for prediction of parotid gland cancer malignancy grade in magnetic resonance images

PURPOSE

The malignancy grades of parotid gland cancer (PGC) have been assessed for a decision of treatment policies. Therefore, we have investigated the feasibility of topology-based radiomic features for the prediction of parotid gland cancer (PGC) malignancy grade in magnetic resonance (MR) images.

MATERIALS AND METHODS

Two-dimensional T1- and T2-weighted MR images of 39 patients with PGC were selected for this study. Imaging properties of PGC can be quantified using the topology, which could be useful for assessing the number of the k-dimensional holes or heterogeneity in PGC regions using invariants of the Betti numbers. Radiomic signatures were constructed from 41,472 features obtained after a harmonization using an elastic net model. PGC patients were stratified using a logistic classification into low/intermediate- and high-grade malignancy groups. The training data were increased by four times to avoid the overfitting problem using a synthetic minority oversampling technique. The proposed approach was assessed using a 4-fold cross-validation test.

RESULTS

The highest accuracy of the proposed approach was 0.975 for the validation cases, whereas that of the conventional approach was 0.694.

CONCLUSION

This study indicated that topology-based radiomic features could be feasible for the noninvasive prediction of the malignancy grade of PGCs.

Three-dimensional topological radiogenomics of epidermal growth factor receptor Del19 and L858R mutation subtypes on computed tomography images of lung cancer patients

OBJECTIVES

To elucidate a novel radiogenomics approach using three-dimensional (3D) topologically invariant Betti numbers (BNs) for topological characterization of epidermal growth factor receptor (EGFR) Del19 and L858R mutation subtypes.

METHODS

In total, 154 patients (wild-type EGFR, 72 patients; Del19 mutation, 45 patients; and L858R mutation, 37 patients) were retrospectively enrolled and randomly divided into 92 training and 62 test cases. Two support vector machine (SVM) models to distinguish between wild-type and mutant EGFR (mutation [M] classification) as well as between the Del19 and L858R subtypes (subtype [S] classification) were trained using 3DBN features. These features were computed from 3DBN maps by using histogram and texture analyses. The 3DBN maps were generated using computed tomography (CT) images based on the Čech complex constructed on sets of points in the images. These points were defined by coordinates of voxels with CT values higher than several threshold values. The M classification model was built using image features and demographic parameters of sex and smoking status. The SVM models were evaluated by determining their classification accuracies. The feasibility of the 3DBN model was compared with those of conventional radiomic models based on pseudo-3D BN (p3DBN), two-dimensional BN (2DBN), and CT and wavelet-decomposition (WD) images. The validation of the model was repeated with 100 times random sampling.

RESULTS

The mean test accuracies for M classification with 3DBN, p3DBN, 2DBN, CT, and WD images were 0.810, 0.733, 0.838, 0.782, and 0.799, respectively. The mean test accuracies for S classification with 3DBN, p3DBN, 2DBN, CT, and WD images were 0.773, 0.694, 0.657, 0.581, and 0.696, respectively.

CONCLUSION

3DBN features, which showed a radiogenomic association with the characteristics of the EGFR Del19/L858R mutation subtypes, yielded higher accuracy for subtype classifications in comparison with conventional features.

CT-derived radiomic analysis for predicting the survival rate of patients with non-small cell lung cancer receiving radiotherapy

BACKGROUND

Radiomics provides an opportunity to minimize adverse effects and optimize the efficacy of treatments noninvasively. This study aims to develop a computed tomography (CT) derived radiomic signature to predict radiological response for the patients with non-small cell lung cancer (NSCLC) receiving radiotherapy.

METHODS

Total 815 NSCLC patients receiving radiotherapy were sourced from public datasets. Using CT images of 281 NSCLC patients, we adopted genetic algorithm to establish a predictive radiomic signature for radiotherapy that had optimal C-index value by Cox model. Survival analysis and receiver operating characteristic curve were performed to estimate the predictive performance of the radiomic signature. Furthermore, radiogenomics analysis was performed in a dataset with matched images and transcriptome data.

RESULTS

Radiomic signature consisting of three features was established and then validated in the validation dataset (log-rank P = 0.0047) including 140 patient, and showed a significant predictive power in two independent datasets totaling 395 NSCLC patients with binary 2-year survival endpoint. Furthermore, the novel proposed radiomic nomogram significantly improved the prognostic performance (concordance index) of clinicopathological factors. Radiogenomics analysis linked our signature with important tumor biological processes (e.g. Mismatch repair, Cell adhesion molecules and DNA replication) associated with clinical outcomes.

CONCLUSIONS

The radiomic signature, reflecting tumor biological processes, could noninvasively predict therapeutic efficacy of NSCLC patients receiving radiotherapy and demonstrate unique advantage for clinical application.

Machine Learning and Deep Learning in Cardiothoracic Imaging: A Scoping Review

Machine-learning (ML) and deep-learning (DL) algorithms are part of a group of modeling algorithms that grasp the hidden patterns in data based on a training process, enabling them to extract complex information from the input data. In the past decade, these algorithms have been increasingly used for image processing, specifically in the medical domain. Cardiothoracic imaging is one of the early adopters of ML/DL research, and the COVID-19 pandemic resulted in more research focus on the feasibility and applications of ML/DL in cardiothoracic imaging. In this scoping review, we systematically searched available peer-reviewed medical literature on cardiothoracic imaging and quantitatively extracted key data elements in order to get a big picture of how ML/DL have been used in the rapidly evolving cardiothoracic imaging field. During this report, we provide insights on different applications of ML/DL and some nuances pertaining to this specific field of research. Finally, we provide general suggestions on how researchers can make their research more than just a proof-of-concept and move toward clinical adoption.

Can Persistent Homology Features Capture More Intrinsic Information about Tumors from 18F-Fluorodeoxyglucose Positron Emission Tomography/Computed Tomography Images of Head and Neck Cancer Patients?

This study hypothesized that persistent homology (PH) features could capture more intrinsic information about the metabolism and morphology of tumors from ¹⁸F-fluorodeoxyglucose positron emission tomography (PET)/computed tomography (CT) images of patients with head and neck (HN) cancer than other conventional features. PET/CT images and clinical variables of 207 patients were selected from the publicly available dataset of the Cancer Imaging Archive. PH images were generated from persistent diagrams obtained from PET/CT images. The PH features were derived from the PH PET/CT images. The signatures were constructed in a training cohort from features from CT, PET, PH-CT, and PH-PET images; clinical variables; and the combination of features and clinical variables. Signatures were evaluated using statistically significant differences (p-value, log-rank test) between survival curves for low- and high-risk groups and the C-index. In an independent test cohort, the signature consisting of PH-PET features and clinical variables exhibited the lowest log-rank p-value of 3.30 × 10⁻⁵ and C-index of 0.80, compared with log-rank p-values from 3.52 × 10⁻² to 1.15 × 10⁻⁴ and C-indices from 0.34 to 0.79 for other signatures. This result suggests that PH features can capture the intrinsic information of tumors and predict prognosis in patients with HN cancer.

RadiomicsJ: a library to compute radiomic features

Despite the widely recognized need for radiomics research, the development and use of full-scale radiomics-based predictive models in clinical practice remains scarce. This is because of the lack of well-established methodologies for radiomic research and the need to develop systems to support radiomic feature calculations and predictive model use. Several excellent programs for calculating radiomic features have been developed. However, there are still issues such as the types of image features, variations in the calculated results, and the limited system environment in which to run the program. Against this background, we developed RadiomicsJ, an open-source radiomic feature computation library. RadiomicsJ will not only be a new research tool to enhance the efficiency of radiomics research but will also become a knowledge resource for medical imaging feature studies through its release as an open-source program.

Relapse predictability of topological signature on pretreatment planning CT images of stage I non-small cell lung cancer patients before treatment with stereotactic ablative radiotherapy

Background

This study aimed to explore the predictability of topological signatures linked to the locoregional relapse (LRR) and distant metastasis (DM) on pretreatment planning computed tomography images of stage I non‐small cell lung cancer (NSCLC) patients before treatment with stereotactic ablative radiotherapy (SABR).

Methods

We divided 125 primary stage I NSCLC patients (LRR: 34, DM: 22) into training (n = 60) and test datasets (n = 65), and the training dataset was augmented to 260 cases using a synthetic minority oversampling technique. The relapse predictabilities of the conventional wavelet‐based features (WF), topology‐based features [BF, Betti number (BN) map features; iBF, inverted BN map features], and their combined features (BWF, iBWF) were compared. The patients were stratified into high‐risk and low‐risk groups using the medians of the radiomics scores in the training dataset.

Results

For the LRR in the test, the iBF, iBWF, and WF showed statistically significant differences (p < 0.05), and the highest nLPC was obtained for the iBF. For the DM in the test, the iBWF showed a significant difference and the highest nLPC.

Conclusion

The iBF indicated the potential of improving the LRR and DM prediction of stage I NSCLC patients prior to undergoing SABR.

Radiomics in Oncological PET Imaging: A Systematic Review-Part 1, Supradiaphragmatic Cancers

Radiomics is an upcoming field in nuclear oncology, both promising and technically challenging. To summarize the already undertaken work on supradiaphragmatic neoplasia and assess its quality, we performed a literature search in the PubMed database up to 18 February 2022. Inclusion criteria were: studies based on human data; at least one specified tumor type; supradiaphragmatic malignancy; performing radiomics on PET imaging. Exclusion criteria were: studies only based on phantom or animal data; technical articles without a clinically oriented question; fewer than 30 patients in the training cohort. A review database containing PMID, year of publication, cancer type, and quality criteria (number of patients, retrospective or prospective nature, independent validation cohort) was constructed. A total of 220 studies met the inclusion criteria. Among them, 119 (54.1%) studies included more than 100 patients, 21 studies (9.5%) were based on prospectively acquired data, and 91 (41.4%) used an independent validation set. Most studies focused on prognostic and treatment response objectives. Because the textural parameters and methods employed are very different from one article to another, it is complicated to aggregate and compare articles. New contributions and radiomics guidelines tend to help improving quality of the reported studies over the years.

A Nomogram Based on CT Deep Learning Signature: A Potential Tool for the Prediction of Overall Survival in Resected Non-Small Cell Lung Cancer Patients

PURPOSE

To develop and further validate a deep learning signature-based nomogram from computed tomography (CT) images for prediction of the overall survival (OS) in resected non-small cell lung cancer (NSCLC) patients.

PATIENTS AND METHODS

A total of 1792 deep learning features were extracted from non-enhanced and venous-phase CT images for each NSCLC patient in training cohort (n=231). Then, a deep learning signature was built with the least absolute shrinkage and selection operator (LASSO) Cox regression model for OS estimation. At last, a nomogram was constructed with the signature and other independent clinical risk factors. The performance of nomogram was assessed by discrimination, calibration and clinical usefulness. In addition, in order to quantify the improvement in performance added by deep learning signature, the net reclassification improvement (NRI) was calculated. The results were validated in external validation cohort (n=77).

RESULTS

A deep learning signature with 9 selected features was significantly associated with OS in both training cohort (hazard ratio [HR]=5.455, 95% CI: 3.393-8.769, P<0.001) and external validation cohort (HR=3.029, 95% CI: 1.673-5.485, P=0.004). The nomogram combining deep learning signature with clinical risk factors of TNM stage, lymphatic vessel invasion and differentiation grade showed favorable discriminative ability with C-index of 0.800 as well as a good calibration, which was validated in external validation cohort (C-index=0.723). Additional value of deep learning signature to the nomogram was statistically significant (NRI=0.093, P=0.027 for training cohort; NRI=0.106, P=0.040 for validation cohort). Decision curve analysis confirmed the clinical usefulness of this nomogram in predicting OS.

CONCLUSION

The deep learning signature-based nomogram is a robust tool for prognostic prediction in resected NSCLC patients.

Exploring the variability of radiomic features of lung cancer lesions on unenhanced and contrast-enhanced chest CT imaging

PURPOSE

The aim of this methods work is to explore the different behavior of radiomic features resulting by using or not the contrast medium in chest CT imaging of non-small cell lung cancer.

METHODS

Chest CT scans, unenhanced and contrast-enhanced, of 17 patients were selected from images collected as part of the staging process. The major T1-T3 lesion was contoured through a semi-automatic approach. These lesions formed the lesion phantoms to study features behavior. The stability of 94 features of the 3D-Slicer package Radiomics was analyzed. Feature discrimination power was quantified by means of Gini's coefficient. Correlation between distance matrices was evaluated through Mantel statistic. Heatmap, cluster and silhouette plots were applied to find well-structured partitions of lesions.

RESULTS

The Gini's coefficient evidenced a low discrimination power, <0.05, for four features and a large discrimination power, around 0.8, for five features. About 90% of features was affected by the contrast medium, masking tumor lesions variability; thirteen features only were found stable. On 8178 combinations of stable features, only one group of four features produced the same partition of lesions with the silhouette width greater than 0.51, both on unenhanced and contrast-enhanced images.

CONCLUSIONS

Gini's coefficient highlighted the features discrimination power in both CT series. Many features were sensitive to the use of the contrast medium, masking the lesions intrinsic variability. Four stable features produced, on both series, the same partition of cancer lesions with reasonable structure; this may merit being objects of further validation studies and interpretative investigations.

Robust radiogenomics approach to the identification of EGFR mutations among patients with NSCLC from three different countries using topologically invariant Betti numbers

Objectives

To propose a novel robust radiogenomics approach to the identification of epidermal growth factor receptor (EGFR) mutations among patients with non-small cell lung cancer (NSCLC) using Betti numbers (BNs).

Materials and methods

Contrast enhanced computed tomography (CT) images of 194 multi-racial NSCLC patients (79 EGFR mutants and 115 wildtypes) were collected from three different countries using 5 manufacturers’ scanners with a variety of scanning parameters. Ninety-nine cases obtained from the University of Malaya Medical Centre (UMMC) in Malaysia were used for training and validation procedures. Forty-one cases collected from the Kyushu University Hospital (KUH) in Japan and fifty-four cases obtained from The Cancer Imaging Archive (TCIA) in America were used for a test procedure. Radiomic features were obtained from BN maps, which represent topologically invariant heterogeneous characteristics of lung cancer on CT images, by applying histogram- and texture-based feature computations. A BN-based signature was determined using support vector machine (SVM) models with the best combination of features that maximized a robustness index (RI) which defined a higher total area under receiver operating characteristics curves (AUCs) and lower difference of AUCs between the training and the validation. The SVM model was built using the signature and optimized in a five-fold cross validation. The BN-based model was compared to conventional original image (OI)- and wavelet-decomposition (WD)-based models with respect to the RI between the validation and the test.

Results

The BN-based model showed a higher RI of 1.51 compared with the models based on the OI (RI: 1.33) and the WD (RI: 1.29).

Conclusion

The proposed model showed higher robustness than the conventional models in the identification of EGFR mutations among NSCLC patients. The results suggested the robustness of the BN-based approach against variations in image scanner/scanning parameters.

Radiomic features based on Hessian index for prediction of prognosis in head-and-neck cancer patients

This study demonstrated the usefulness of radiomic features based on the Hessian index of differential topology for the prediction of prognosis prior to treatment in head-and-neck (HN) cancer patients. The Hessian index, which can indicate tumor heterogeneity with convex, concave, and other points (saddle points), was calculated as the number of negative eigenvalues of the Hessian matrix at each voxel on computed tomography (CT) images. Three types of signatures were constructed in a training cohort (n = 126), one type each from CT conventional features, Hessian index features, and combined features from the conventional and index feature sets. The prognostic value of the signatures were evaluated using statistically significant difference (p value, log-rank test) to compare the survival curves of low- and high-risk groups. In a test cohort (n = 68), the p values of the models built with conventional, index, combined features, and clinical variables were 2.95 [Formula: see text] 10-2, 1.85 [Formula: see text] 10-2, 3.17 [Formula: see text] 10-2, and 1.87 [Formula: see text] 10-3, respectively. When the features were integrated with clinical variables, the p values of conventional, index, and combined features were 3.53 [Formula: see text] 10-3, 1.28 [Formula: see text] 10-3, and 1.45 [Formula: see text] 10-3, respectively. This result indicates that index features could provide more prognostic information than conventional features and further increase the prognostic value of clinical variables in HN cancer patients.