Dosiomics: Extracting 3D Spatial Features From Dose Distribution to Predict Incidence of Radiation Pneumonitis

Plan complexity and dosiomics signatures for gamma passing rate classification in volumetric modulated arc therapy: External validation across different LINACs

PURPOSE

This study aims to enhance gamma passing rate (GPR) classification by integrating plan complexity signature, dosiomics signature, and comprehensive plan parameters, and to validate this method using data from different linear accelerators (LINACs).

METHODS

This study included 235 volumetric modulated arc therapy (VMAT) treatment plans delivered using the TrueBeam LINAC as the primary dataset, along with 47 plans from the VitalBeam LINAC for external validation. The primary dataset was split into training (N = 166) and test (N = 69) subsets. Extracted features included 47 plan complexity metrics, 851 dosiomics features, and 20 plan parameters. Plan complexity score (PCscore) and dosiomics score (Doscore) were derived using the least absolute shrinkage and selection operator (LASSO). Four classification models were developed by combining PCscore, Doscore, and plan parameters according to a gamma criterion of 2 %/2 mm (γ2%/2 mm). A nomogram was constructed to combine these signatures with plan parameters. Model performance was evaluated using receiver operating characteristic (ROC) curves, calibration curves, and decision curve analysis (DCA).

RESULTS

The combined model incorporating PCscore, Doscore, and plan parameters exhibited high discriminative power, with areas under the curve (AUC) of 0.894, 0.899, and 0.904 for the training, test, and external datasets, respectively. At γ3%/2 mm, the model maintained robust performance with AUCs of 0.842 and 0.833 in the test and external datasets. Calibration curves and DCA validated the model's effectiveness.

CONCLUSIONS

Integrating plan complexity and dosiomics signatures with key plan parameters significantly improves GPR classification for VMAT treatment plans, offering a robust approach for patient-specific quality assurance (PSQA).

Radiomics and dosiomics approaches to estimate lung function after stereotactic body radiation therapy in patients with lung tumors

Lung function assessment is essential for determining the optimal treatment strategy for radiation therapy in patients with lung tumors. This study aimed to develop radiomics and dosiomics approaches to estimate pulmonary function test (PFT) results in post-stereotactic body radiation therapy (SBRT). Sixty-four patients with lung tumors who underwent SBRT were included. Models were created to estimate the PFT results at 0-6 months (Cohort 1) and 6-24 months (Cohort 2) after SBRT. Radiomics and dosiomics features were extracted from the computed tomography (CT) images and dose distributions, respectively. To estimate the PFT results, Models A (dose-volume histogram [DVH] + radiomics features) and B (DVH + radiomics + dosiomics features) were created. In the PFT results, the forced expiratory volume in one second (FEV1) and forced vital capacity (FVC) were estimated using each model, and the ratio of FEV1 to FVC (FEV1/FVC) was calculated. The Pearson's correlation coefficient (Pearson r) and area under the curve (AUC) for FEV1/FVC (< 70%) were calculated. The models were evaluated by comparing them with the conventional calculation formulae (Conventional). The Pearson r (FEV1/FVC) values were 0.30, 0.64, and 0.69 for Conventional and Models A and B (Cohort 2), respectively, and the AUC (FEV1/FVC < 70%) values were 0.63, 0.80, and 0.78, respectively. This study demonstrates the possibility of estimating lung function after SBRT using radiomics and dosiomics features based on planning CT images and dose distributions.

Developing a novel dosiomics model to predict treatment failures following lung stereotactic body radiation therapy

Purpose

The purpose of this study was to investigate the dosiomics features of the interplay between CT density and dose distribution in lung SBRT plans, and to develop a model to predict treatment failure following lung SBRT treatment.

Methods

A retrospective study was conducted involving 179 lung cancer patients treated with SBRT at the University of Nebraska Medical Center (UNMC) between October 2007 and June 2022. Features from the CT image, Biological Effective Dose (BED) and five interaction matrices between CT and BED were extracted using radiomics mathematics. Our in-house feature selection pipeline was utilized to evaluate and rank features based on their importance and redundancy, with only the selected non-redundant features being used for predictive modeling. We randomly selected 151 cases and 28 cases as training and test datasets. Four different models were trained utilizing the Balanced Random Forest framework on the same training dataset to differentiate between failure and non-failure cases. These four models utilized the same number of selected features extracted from CT-only, BED-only, a combination of CT and BED, and a composite of CT and BED including their interaction matrices, respectively.

Results

The cohort included 125 non-failure cases and 54 failure cases, with a median follow-up time of 34.4 months. We selected the top 17 important and non-redundant features (with the Pearsons's coefficient < 0.5) in each model. When evaluated on the same independent test set, the four models-CT features-only, BED features-only, a combination of CT and BED features, and a composite model including features from CT and BED that includes their interaction matrices-achieved AUC values of 0.56, 0.75, 0.73, and 0.82, respectively, with corresponding accuracies of 0.61, 0.79, 0.71, and 0.79. The composite model demonstrated the highest AUC and accuracy, indicating that incorporating interactions between CT and BED reveals more predictive capabilities in distinguishing between failure and non-failure cases.

Conclusion

The dosiomics model integrating the interaction between CT and Dose can effectively predict treatment failure following lung SBRT treatment and may serve as a useful tool to proactively evaluate and select lung SBRT treatment plans to reduce treatment failure in the future.

Combining dosiomics and machine learning methods for predicting severe cardiac diseases in childhood cancer survivors: the French Childhood Cancer Survivor Study

Background

Cardiac disease (CD) is a primary long-term diagnosed pathology among childhood cancer survivors. Dosiomics (radiomics extracted from the dose distribution) have received attention in the past few years to assess better the induced risk of radiotherapy (RT) than standard dosimetric features such as dose-volume indicators. Hence, using the spatial information contained in the dosiomics features with machine learning methods may improve the prediction of CD.

Methods

We considered the 7670 5-year survivors of the French Childhood Cancer Survivors Study (FCCSS). Dose-volume and dosiomics features are extracted from the radiation dose distribution of 3943 patients treated with RT. Survival analysis is performed considering several groups of features and several models [Cox Proportional Hazard with Lasso penalty, Cox with Bootstrap Lasso selection, Random Survival Forests (RSF)]. We establish the performance of dosiomics compared to baseline models by estimating C-index and Integrated Brier Score (IBS) metrics with 5-fold stratified cross-validation and compare their time-dependent error curves.

Results

An RSF model adjusted on the first-order dosiomics predictors extracted from the whole heart performed best regarding the C-index (0.792 ± 0.049), and an RSF model adjusted on the first-order dosiomics predictors extracted from the heart's subparts performed best regarding the IBS (0.069 ± 0.05). However, the difference is not statistically significant with the standard models (C-index of Cox PH adjusted on dose-volume indicators: 0.791 ± 0.044; IBS of Cox PH adjusted on the mean dose to the heart: 0.074 ± 0.056).

Conclusion

In this study, dosiomics models have slightly better performance metrics but they do not outperform the standard models significantly. Quantiles of the dose distribution may contain enough information to estimate the risk of late radio-induced high-grade CD in childhood cancer survivors.

Dosiomics-based detection of dose distribution variations in helical tomotherapy for prostate cancer patients: influence of treatment plan parameters

The stability of dosiomics features (DFs) and dose-volume histogram (DVH) parameters for detecting disparities in helical tomotherapy planned dose distributions was assessed. Treatment plans of 18 prostate patients were recalculated using the followings: field width (WF) (2.5 vs. 5), pitch factor (PF) (0.433 vs. 0.444), and modulation factor (MF) (2.5 vs. 3). From each of the eight plans per patient, ninety-three original and 744 wavelet-based DFs were extracted, using 3D-Slicer software, across six regions including: target volume (PTV), pelvic lymph nodes (PTV-LN), PTV + PTV-LN (PTV-All), one cm rind around PTV-All (PTV-Ring), rectum, and bladder. For the resulting DFs and DVH parameters, the coefficient of variation (CV) was calculated, and using hierarchical clustering, the features were classified into low/high variability. The significance of parameters on instability was analyzed by a three-way analysis of variance. All DF's were stable in PTV, PTV-LN, and PTV-Ring (average CV (CV ¯ )  ≤ 0.36). Only one feature in the bladder (CV ¯  = 0.9), rectum (CV ¯  = 0.4), and PTV-All (CV ¯  = 0.37) were considered unstable due to change in MF in the bladder and WF in the PTV-All. The value ofCV ¯ for the wavelet features was much higher than that for the original features. Out of 225 unstable wavelet features, 84 features hadCV ¯  ≥ 1. The CVs for all the DVHs remained very small (CV ¯ < 0.06). This study highlights that the sensitivity of DFs to changes in tomotherapy planning parameters is influenced by the region and the DFs, particularly wavelet features, surpassing the effectiveness of DVHs.

Predicting grade II-IV bone marrow suppression in patients with cervical cancer based on radiomics and dosiomics

Objective

The objective of this study is to develop a machine learning model integrating clinical characteristics with radiomics and dosiomics data, aiming to assess their predictive utility in anticipating grade 2 or higher BMS occurrences in cervical cancer patients undergoing radiotherapy.

Methods

A retrospective analysis was conducted on the clinical data, planning CT images, and radiotherapy planning documents of 106 cervical cancer patients who underwent radiotherapy at our hospital. The patients were randomly divided into training set and test set in an 8:2 ratio. The radiomic features and dosiomic features were extracted from the pelvic bone marrow (PBM) of planning CT images and radiotherapy planning documents, and the least absolute shrinkage and selection operator (LASSO) algorithm was employed to identify the best predictive characteristics. Subsequently, the dosiomic score (D-score) and the radiomic score (R-score) was calculated. Clinical predictors were identified through both univariate and multivariate logistic regression analysis. Predictive models were constructed by intergrating clinical predictors with DVH parameters, combining DVH parameters and R-score with clinical predictors, and amalgamating clinical predictors with both D-score and R-score. The predictive model's efficacy was assessed by plotting the receiver operating characteristic (ROC) curve and evaluating its performance through the area under the ROC curve (AUC), the calibration curve, and decision curve analysis (DCA).

Results

Seven radiomic features and eight dosiomic features exhibited a strong correlation with the occurrence of BMS. Through univariate and multivariate logistic regression analyses, age, planning target volume (PTV) size and chemotherapy were identified as clinical predictors. The AUC values for the training and test sets were 0.751 and 0.743, respectively, surpassing those of clinical DVH R-score model (AUC=0.707 and 0.679) and clinical DVH model (AUC=0.650 and 0.638). Furthermore, the analysis of both the calibration and the DCA suggested that the combined model provided superior calibration and demonstrated a higher net clinical benefit.

Conclusion

The combined model is of high diagnostic value in predicting the occurrence of BMS in patients with cervical cancer during radiotherapy.

Multiomics-Based Outcome Prediction in Personalized Ultra-Fractionated Stereotactic Adaptive Radiotherapy (PULSAR)

Objectives: This retrospective study aims to develop a multiomics approach that integrates radiomics, dosiomics, and delta features to predict treatment responses in brain metastasis (BM) patients undergoing PULSAR. Methods: A retrospective study encompassing 39 BM patients with 69 lesions treated with PULSAR was undertaken. Radiomics, dosiomics, and delta features were extracted from both pre-treatment and intra-treatment MRI scans alongside dose distributions. Six individual models, alongside an ensemble feature selection (EFS) model, were evaluated. The classification task focused on distinguishing between two lesion groups based on whether they exhibited a volume reduction of more than 20% at follow-up. Performance metrics, including sensitivity, specificity, accuracy, precision, F1 score, and the area under the receiver operating characteristic (ROC) curve (AUC), were assessed. Results: The EFS model integrated the features from pre-treatment radiomics, pre-treatment dosiomics, intra-treatment radiomics, and delta radiomics. It outperformed six individual models, achieving an AUC of 0.979, accuracy of 0.917, and F1 score of 0.821. Among the top nine features of the EFS model, six features came from post-wavelet transformation and three from original images. Conclusions: The study demonstrated the feasibility of employing a data-driven multiomics approach to predict treatment outcomes in BM patients receiving PULSAR treatment. Integrating multiomics with intra-treatment decision support in PULSAR shows promise for optimizing patient management and reducing the risks of under- or over-treatment.

Evaluating the Predictive Value of the Dose-volume Parameters and Vascular Endothelial Growth Factor Expression on Rectal Toxicity in Prostate Cancer Patients

Background

The response to radiation varies widely among people. Despite recent advances to conform the dose distribution to the tumor volume, it is still impossible to perform radiotherapy based on the biological characteristics of each individual. In this case, identifying a biomarker can be a step toward personalizing treatment. This research was carried out to evaluate the predictive value of dose-volume parameters and vascular endothelial growth factor (VEGF) expression on rectal proctitis toxicity in prostate cancer patients.

Materials and Methods

Eighteen patients with prostate cancer who underwent helical tomotherapy were included in the study. VEGF serum level before and after treatment was obtained. Furthermore, dosimetric parameters, including rectal volume, maximum dose, V50, V60, V65, V70, and V75 were extracted. Spearman's correlation coefficient of VEGF-related and dosimetric parameters with the grade ≥1 rectal proctitis was calculated.

Results

In the rectal toxicity group, the mean value of VEGF increased significantly after treatment compared to before (P = 0.008). Despite lower values of pre- and post-treatment VEGF in the toxicity group, this difference was not statistically significant (P > 0.05). Among the dosimetric parameters, only V65 had a significantly higher value in the toxicity group (P = 0.033). The highest correlation coefficients were obtained for pretreatment VEGF values and V65 (-0.446 and 0.450).

Conclusion

The results of the study confirm the correlation of VEGF expression with the pathobiology process of rectal radiation proctitis. However, the pathobiology process of radiation proctitis is complicated. More research is needed to prove the involvement of VEGF expression in the early detection of proctitis.

Particle Beam Radiobiology Status and Challenges: A PTCOG Radiobiology Subcommittee Report

Particle therapy (PT) represents a significant advancement in cancer treatment, precisely targeting tumor cells while sparing surrounding healthy tissues thanks to the unique depth-dose profiles of the charged particles. Furthermore, their linear energy transfer and relative biological effectiveness enhance their capability to treat radioresistant tumors, including hypoxic ones. Over the years, extensive research has paved the way for PT's clinical application, and current efforts aim to refine its efficacy and precision, minimizing the toxicities. In this regard, radiobiology research is evolving toward integrating biotechnology to advance drug discovery and radiation therapy optimization. This shift from basic radiobiology to understanding the molecular mechanisms of PT aims to expand the therapeutic window through innovative dose delivery regimens and combined therapy approaches. This review, written by over 30 contributors from various countries, provides a comprehensive look at key research areas and new developments in PT radiobiology, emphasizing the innovations and techniques transforming the field, ranging from the radiobiology of new irradiation modalities to multimodal radiation therapy and modeling efforts. We highlight both advancements and knowledge gaps, with the aim of improving the understanding and application of PT in oncology.

Multi-omics deep learning for radiation pneumonitis prediction in lung cancer patients underwent volumetric modulated arc therapy

BACKGROUND AND OBJECTIVE

To evaluate the feasibility and accuracy of radiomics, dosiomics, and deep learning (DL) in predicting Radiation Pneumonitis (RP) in lung cancer patients underwent volumetric modulated arc therapy (VMAT) to improve radiotherapy safety and management.

METHODS

Total of 318 and 31 lung cancer patients underwent VMAT from First Affiliated Hospital of Wenzhou Medical University (WMU) and Quzhou Affiliated Hospital of WMU were enrolled for training and external validation, respectively. Models based on radiomics (R), dosiomics (D), and combined radiomics and dosiomics features (R+D) were constructed and validated using three machine learning (ML) methods. DL models trained with CT (DLR), dose distribution (DLD), and combined CT and dose distribution (DL(R+D)) images were constructed. DL features were then extracted from the fully connected layers of the best-performing DL model to combine with features of the ML model with the best performance to construct models of R+DLR, D+DLD, R+D+DL(R+D)) for RP prediction.

RESULTS

The R+D model achieved a best area under curve (AUC) of 0.84, 0.73, and 0.73 in the internal validation cohorts with Support Vector Machine (SVM), XGBoost, and Logistic Regression (LR), respectively. The DL(R+D) model achieved a best AUC of 0.89 and 0.86 using ResNet-34 in training and internal validation cohorts, respectively. The R+D+DL(R+D) model achieved a best performance in the external validation cohorts with an AUC, accuracy, sensitivity, and specificity of 0.81(0.62-0.99), 0.81, 0.84, and 0.67, respectively.

CONCLUSIONS

The integration of radiomics, dosiomics, and DL features is feasible and accurate for the RP prediction to improve the management of lung cancer patients underwent VMAT.

Comparative Analysis of Repeatability in CT Radiomics and Dosiomics Features under Image Perturbation: A Study in Cervical Cancer Patients

This study aims to evaluate the repeatability of radiomics and dosiomics features via image perturbation of patients with cervical cancer. A total of 304 cervical cancer patients with planning CT images and dose maps were retrospectively included. Random translation, rotation, and contour randomization were applied to CT images and dose maps before radiomics feature extraction. The repeatability of radiomics and dosiomics features was assessed using intra-class correlation of coefficient (ICC). Pearson correlation coefficient (r) was adopted to quantify the correlation between the image characteristics and feature repeatability. In general, the repeatability of dosiomics features was lower compared with CT radiomics features, especially after small-sigma Laplacian-of-Gaussian (LoG) and wavelet filtering. More repeatable features (ICC > 0.9) were observed when extracted from the original, Large-sigma LoG filtered, and LLL-/LLH-wavelet filtered images. Positive correlations were found between image entropy and high-repeatable feature number in both CT and dose (r = 0.56, 0.68). Radiomics features showed higher repeatability compared to dosiomics features. These findings highlight the potential of radiomics features for robust quantitative imaging analysis in cervical cancer patients, while suggesting the need for further refinement of dosiomics approaches to enhance their repeatability.

Utilizing radiomics and dosiomics with AI for precision prediction of radiation dermatitis in breast cancer patients

PURPOSE

This study explores integrating clinical features with radiomic and dosiomic characteristics into AI models to enhance the prediction accuracy of radiation dermatitis (RD) in breast cancer patients undergoing volumetric modulated arc therapy (VMAT).

MATERIALS AND METHODS

This study involved a retrospective analysis of 120 breast cancer patients treated with VMAT at Kaohsiung Veterans General Hospital from 2018 to 2023. Patient data included CT images, radiation doses, Dose-Volume Histogram (DVH) data, and clinical information. Using a Treatment Planning System (TPS), we segmented CT images into Regions of Interest (ROIs) to extract radiomic and dosiomic features, focusing on intensity, shape, texture, and dose distribution characteristics. Features significantly associated with the development of RD were identified using ANOVA and LASSO regression (p-value < 0.05). These features were then employed to train and evaluate Logistic Regression (LR) and Random Forest (RF) models, using tenfold cross-validation to ensure robust assessment of model efficacy.

RESULTS

In this study, 102 out of 120 VMAT-treated breast cancer patients were included in the detailed analysis. Thirty-two percent of these patients developed Grade 2+ RD. Age and BMI were identified as significant clinical predictors. Through feature selection, we narrowed down the vast pool of radiomic and dosiomic data to 689 features, distributed across 10 feature subsets for model construction. In the LR model, the J subset, comprising DVH, Radiomics, and Dosiomics features, demonstrated the highest predictive performance with an AUC of 0.82. The RF model showed that subset I, which includes clinical, radiomic, and dosiomic features, achieved the best predictive accuracy with an AUC of 0.83. These results emphasize that integrating radiomic and dosiomic features significantly enhances the prediction of Grade 2+ RD.

CONCLUSION

Integrating clinical, radiomic, and dosiomic characteristics into AI models significantly improves the prediction of Grade 2+ RD risk in breast cancer patients post-VMAT. The RF model analysis demonstrates that a comprehensive feature set maximizes predictive efficacy, marking a promising step towards utilizing AI in radiation therapy risk assessment and enhancing patient care outcomes.

A dosiomics approach to treatment outcome modeling in carbon ion radiotherapy for skull base chordomas

PURPOSE

To investigate the role of dosiomics features extracted from physical dose (DPHYS), RBE-weighted dose (DRBE) and dose-averaged Linear Energy Transfer (LETd), to predict the risk of local recurrence (LR) in skull base chordoma (SBC) treated with Carbon Ion Radiotherapy (CIRT). Thus, define and evaluate dosiomics-driven tumor control probability (TCP) models.

MATERIALS AND METHODS

54 SBC patients were retrospectively selected for this study. A regularized Cox proportional hazard model (r-Cox) and Survival Support Vector Machine (s-SVM) were tuned within a repeated Cross Validation (CV) and patients were stratified in low/high risk of LR. Models' performance was evaluated through Harrell's concordance statistic (C-index), and survival was represented through Kaplan-Meier (KM) curves. A multivariable logistic regression was fit to the selected feature sets to generate a dosiomics-driven TCP model for each map. These were compared to a reference model built with clinical parameters in terms of f-score and accuracy.

RESULTS

The LETd maps reached a test C-index of 0.750 and 0.786 with r-Cox and s-SVM, and significantly separated KM curves. DPHYS maps and clinical parameters showed promising CV outcomes with C-index above 0.8, despite a poorer performance on the test set and patients stratification. The LETd-based TCP showed a significatively higher f-score (0.67[0.52-0.70], median[IQR]) compared to the clinical model (0.4[0.32-0.63], p < 0.025), while DPHYS achieved a significatively higher accuracy (DPHYS: 0.73[0.65-0.79], Clinical: 0.6 [0.52-0.72]).

CONCLUSION

This analysis supports the role of LETd as relevant source of prognostic factors for LR in SBC treated with CIRT. This is reflected in the TCP modeling, where LETd and DPHYS showed an improved performance with respect to clinical models.

Applications of artificial intelligence for machine- and patient-specific quality assurance in radiation therapy: current status and future directions

Machine- and patient-specific quality assurance (QA) is essential to ensure the safety and accuracy of radiotherapy. QA methods have become complex, especially in high-precision radiotherapy such as intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT), and various recommendations have been reported by AAPM Task Groups. With the widespread use of IMRT and VMAT, there is an emerging demand for increased operational efficiency. Artificial intelligence (AI) technology is quickly growing in various fields owing to advancements in computers and technology. In the radiotherapy treatment process, AI has led to the development of various techniques for automated segmentation and planning, thereby significantly enhancing treatment efficiency. Many new applications using AI have been reported for machine- and patient-specific QA, such as predicting machine beam data or gamma passing rates for IMRT or VMAT plans. Additionally, these applied technologies are being developed for multicenter studies. In the current review article, AI application techniques in machine- and patient-specific QA have been organized and future directions are discussed. This review presents the learning process and the latest knowledge on machine- and patient-specific QA. Moreover, it contributes to the understanding of the current status and discusses the future directions of machine- and patient-specific QA.

Synergizing the interaction of single nucleotide polymorphisms with dosiomics features to build a dual-omics model for the prediction of radiation pneumonitis

OBJECTIVE

Radiation pneumonitis (RP) is the major dose-limiting toxicity of thoracic radiotherapy. This study aimed to developed a dual-omics (single nucleotide polymorphisms, SNP and dosiomics) prediction model for symptomatic RP.

MATERIALS AND METHODS

The potential SNPs, which are of significant difference between the RP grade ≥ 3 group and the RP grade ≤ 1 group, were selected from the whole exome sequencing SNPs using the Fisher's exact test. Patients with lung cancer who received thoracic radiotherapy at our institution from 2009 to 2016 were enrolled for SNP selection and model construction. The factorization machine (FM) method was used to model the SNP epistasis effect, and to construct the RP prediction model (SNP-FM). The dosiomics features were extracted, and further selected using the minimum redundancy maximum relevance (mRMR) method. The selected dosiomics features were added to the SNP-FM model to construct the dual-omics model.

RESULTS

For SNP screening, peripheral blood samples of 28 patients with RP grade ≥ 3 and the matched 28 patients with RP grade ≤ 1 were sequenced. 81 SNPs were of significant difference (P < 0.015) and considered as potential SNPs. In addition, 21 radiation toxicity related SNPs were also included. For model construction, 400 eligible patients (including 108 RP grade ≥ 2) were enrolled. Single SNP showed no strong correlation with RP. On the other hand, the SNP-SNP interaction (epistasis effect) of 19 SNPs were modeled by the FM method, and achieved an area under the curve (AUC) of 0.76 in the testing group. In addition, 4 dosiomics features were selected and added to the model, and increased the AUC to 0.81.

CONCLUSIONS

A novel dual-omics model by synergizing the SNP epistasis effect with dosiomics features was developed. The enhanced the RP prediction suggested its promising clinical utility in identifying the patients with severe RP during thoracic radiotherapy.

Radiation pneumonitis prediction with dual-radiomics for esophageal cancer underwent radiotherapy

BACKGROUND

To integrate radiomics and dosiomics features from multiple regions in the radiation pneumonia (RP grade ≥ 2) prediction for esophageal cancer (EC) patients underwent radiotherapy (RT).

METHODS

Total of 143 EC patients in the authors' hospital (training and internal validation: 70%:30%) and 32 EC patients from another hospital (external validation) underwent RT from 2015 to 2022 were retrospectively reviewed and analyzed. Patients were dichotomized as positive (RP+) or negative (RP-) according to CTCAE V5.0. Models with radiomics and dosiomics features extracted from single region of interest (ROI), multiple ROIs and combined models were constructed and evaluated. A nomogram integrating radiomics score (Rad_score), dosiomics score (Dos_score), clinical factors, dose-volume histogram (DVH) factors, and mean lung dose (MLD) was also constructed and validated.

RESULTS

Models with Rad_score_Lung&Overlap and Dos_score_Lung&Overlap achieved a better area under curve (AUC) of 0.818 and 0.844 in the external validation in comparison with radiomics and dosiomics models with features extracted from single ROI. Combining four radiomics and dosiomics models using support vector machine (SVM) improved the AUC to 0.854 in the external validation. Nomogram integrating Rad_score, and Dos_score with clinical factors, DVH factors, and MLD further improved the RP prediction AUC to 0.937 and 0.912 in the internal and external validation, respectively.

CONCLUSION

CT-based RP prediction model integrating radiomics and dosiomics features from multiple ROIs outperformed those with features from a single ROI with increased reliability for EC patients who underwent RT.

Improvement of Prediction Performance for Radiation Pneumonitis by Using 3-Dimensional Dosiomic Features

INTRODUCTION

Patients with early non-small-cell lung cancer (NSCLC) have a relatively long survival time after stereotactic body radiation therapy (SBRT). Predicting radiation-induced pneumonia (RP) has important clinical and social implications for improving the quality of life of such patients. This study developed an RP prediction model by using 3-dimensional (3D) dosiomic features. The model can be used to guide radiation therapy to reduce toxicity.

METHODS

Radiomic features were extracted from pre-treatment CT, dose-volume histogram (DVH) parameters and dosiomic features were extracted from the 3D dose distribution of 140 lung cancer patients. Four predictive models: (1) CT; (2) CT + DVH; (3) CT + Rtdose; and (4) Hybrid, CT + DVH + Rtdose, were trained to predict symptomatic RP by extremely randomized trees. Accuracy, sensitivity, specificity, and area under the receiver operator characteristic curve were evaluated.

RESULT

Results showed that the fraction regimen was correlated with symptomatic RP (P < .001). The proposed model achieved promising prediction results. The performance metrics for CT, CT + DVH, CT + Rtdose, and Hybrid were as follows: accuracy: 0.786, 0.821, 0.821, and 0.857; sensitivity: 0.625, 1, 0.875, and 1; specificity: 0.8, 0.565, 0.5, and 0.875; and area under the receiver operator characteristic curve: 0.791, 0.809, 0.907, and 0.920, respectively.

CONCLUSION

Dosiomic features can improve the performance of the predictive model for symptomatic RP compared with that obtained with the CT + DVH model. The model proposed in this study can help radiation oncologists individually predict the incidence rate of RP.

A dosiomics model for prediction of radiation-induced acute skin toxicity in breast cancer patients: machine learning-based study for a closed bore linac

BACKGROUND

Radiation induced acute skin toxicity (AST) is considered as a common side effect of breast radiation therapy. The goal of this study was to design dosiomics-based machine learning (ML) models for prediction of AST, to enable creating optimized treatment plans for high-risk individuals.

METHODS

Dosiomics features extracted using Pyradiomics tool (v3.0.1), along with treatment plan-derived dose volume histograms (DVHs), and patient-specific treatment-related (PTR) data of breast cancer patients were used for modeling. Clinical scoring was done using the Common Terminology Criteria for Adverse Events (CTCAE) V4.0 criteria for skin-specific symptoms. The 52 breast cancer patients were grouped into AST 2 + (CTCAE ≥ 2) and AST 2 - (CTCAE < 2) toxicity grades to facilitate AST modeling. They were randomly divided into training (70%) and testing (30%) cohorts. Multiple prediction models were assessed through multivariate analysis, incorporating different combinations of feature groups (dosiomics, DVH, and PTR) individually and collectively. In total, seven unique combinations, along with seven classification algorithms, were considered after feature selection. The performance of each model was evaluated on the test group using the area under the receiver operating characteristic curve (AUC) and f1-score. Accuracy, precision, and recall of each model were also studied. Statistical analysis involved features differences between AST 2 - and AST 2 + groups and cutoff value calculations.

RESULTS

Results showed that 44% of the patients developed AST 2 + after Tomotherapy. The dosiomics (DOS) model, developed using dosiomics features, exhibited a noteworthy improvement in AUC (up to 0.78), when spatial information is preserved in the dose distribution, compared to DVH features (up to 0.71). Furthermore, a baseline ML model created using only PTR features for comparison with DOS models showed the significance of dosiomics in early AST prediction. By employing the Extra Tree (ET) classifiers, the DOS + DVH + PTR model achieved a statistically significant improved performance in terms of AUC (0.83; 95% CI 0.71-0.90), accuracy (0.70), precision (0.74) and sensitivity (0.72) compared to other models.

CONCLUSIONS

This study confirmed the benefit of dosiomics-based ML in the prediction of AST. However, the combination of dosiomics, DVH, and PTR yields significant improvement in AST prediction. The results of this study provide the opportunity for timely interventions to prevent the occurrence of radiation induced AST.

A deep learning-based method for the prediction of temporal lobe injury in patients with nasopharyngeal carcinoma

PURPOSE

To establish a deep learning-based model to predict radiotherapy-induced temporal lobe injury (TLI).

MATERIALS AND METHODS

Spatial features of dose distribution within the temporal lobe were extracted using both the three-dimensional convolution (C3D) network and the dosiomics method. The Minimal Redundancy-Maximal-Relevance (mRMR) method was employed to rank the extracted features and select the most relevant ones. Four machine learning (ML) classifiers, including logistic regression (LR), k-nearest neighbors (kNN), support vector machines (SVM) and random forest (RF), were used to establish prediction models. Nested sampling and hyperparameter tuning methods were applied to train and validate the prediction models. For comparison, a prediction model base on the conventional D0.5cc of the temporal lobe obtained from dose volume (DV) histogram was established. The area under the receiver operating characteristic (ROC) curve (AUC) was utilized to compare the predictive performance of the different models.

RESULTS

A total of 127 nasopharyngeal carcinoma (NPC) patients were included in the study. In the model based on C3D deep learning features, the highest AUC value of 0.843 was achieved with 5 features. For the dosiomics features model, the highest AUC value of 0.715 was attained with 1 feature. Both of these models demonstrated superior performance compared to the prediction model based on DV parameters, which yielded an AUC of 0.695.

CONCLUSION

The prediction model utilizing C3D deep learning features outperformed models based on dosiomics features or traditional parameters in predicting the onset of TLI. This approach holds promise for predicting radiation-induced toxicities and guide individualized radiotherapy.

Contralateral Hypertrophy Post Yttrium-90 Transarterial Radioembolization in Patients With Hepatocellular Carcinoma and Portal Vein Tumor Thrombus

Objectives Contralateral hypertrophy of non-irradiated liver following Yttrium-90 (90Y) transarterial radioembolization (TARE) is increasingly recognized as an option to facilitate curative surgical resection in patients that would otherwise not be surgical candidates due to a small future liver remnant (FLR). This study aimed to investigate the correlation between patient features and liver hypertrophy and identify potential predictors for liver growth in patients with hepatocellular carcinoma (HCC) and portal vein tumor thrombus (PVTT) undergoing TARE. Methodology Twenty-three patients with HCC and PVTT were included. Contralateral liver hypertrophy was assessed at six months posttreatment based on CT or MRI imaging. Thirteen patient features were selected for statistical and prediction analysis. Univariate Spearman correlation and analysis of variance (ANOVA) tests were performed. Subsequently, four feature-selection methods based on multivariate analysis were used to improve model generalization performance. The selected features were applied to train linear regression models, with fivefold cross-validation to assess the performance of the predicted models. Results The ratio of disease-free target liver volume to spared liver volume and total liver volume showed the highest correlations with contralateral hypertrophy (P-values = 0.03 and 0.05, respectively). In three out of four feature-selection methods, the feature of disease-free target liver volume to total liver volume ratio was selected, having positive correlations with the outcome and suggesting that more hypertrophy may be expected when more volume of disease-free liver is irradiated. Conclusions Contralateral hypertrophy post-90Y TARE can be an option for facilitating surgical resection in patients with otherwise small FLR.

From pixels to prognosis: unveiling radiomics models with SHAP and LIME for enhanced interpretability

Radiomics-based prediction models have shown promise in predicting Radiation Pneumonitis (RP), a common adverse outcome of chest irradiation. Τhis study looks into more than just RP: it also investigates a bigger shift in the way radiomics-based models work. By integrating multi-modal radiomic data, which includes a wide range of variables collected from medical images including cutting-edge PET/CT imaging, we have developed predictive models that capture the intricate nature of illness progression. Radiomic features were extracted using PyRadiomics, encompassing intensity, texture, and shape measures. The high-dimensional dataset formed the basis for our predictive models, primarily Gradient Boosting Machines (GBM)-XGBoost, LightGBM, and CatBoost. Performance evaluation metrics, including Multi-Modal AUC-ROC, Sensitivity, Specificity, and F1-Score, underscore the superiority of the Deep Neural Network (DNN) model. The DNN achieved a remarkable Multi-Modal AUC-ROC of 0.90, indicating superior discriminatory power. Sensitivity and specificity values of 0.85 and 0.91, respectively, highlight its effectiveness in detecting positive occurrences while accurately identifying negatives. External validation datasets, comprising retrospective patient data and a heterogeneous patient population, validate the robustness and generalizability of our models. The focus of our study is the application of sophisticated model interpretability methods, namely SHAP (SHapley Additive exPlanations) and LIME (Local Interpretable Model-Agnostic Explanations), to improve the clarity and understanding of predictions. These methods allow clinicians to visualize the effects of features and provide localized explanations for every prediction, enhancing the comprehensibility of the model. This strengthens trust and collaboration between computational technologies and medical competence. The integration of data-driven analytics and medical domain expertise represents a significant shift in the profession, advancing us from analyzing pixel-level information to gaining valuable prognostic insights.

Dosiomics and radiomics-based prediction of pneumonitis after radiotherapy and immune checkpoint inhibition: The relevance of fractionation

OBJECTIVES

Post-therapy pneumonitis (PTP) is a relevant side effect of thoracic radiotherapy and immunotherapy with checkpoint inhibitors (ICI). The influence of the combination of both, including dose fractionation schemes on PTP development is still unclear. This study aims to improve the PTP risk estimation after radio(chemo)therapy (R(C)T) for lung cancer with and without ICI by investigation of the impact of dose fractionation on machine learning (ML)-based prediction.

MATERIALS AND METHODS

Data from 100 patients who received fractionated R(C)T were collected. 39 patients received additional ICI therapy. Computed Tomography (CT), RT segmentation and dose data were extracted and physical doses were converted to 2-Gy equivalent doses (EQD2) to account for different fractionation schemes. Features were reduced using Pearson intercorrelation and the Boruta algorithm within 1000-fold bootstrapping. Six single (clinics, Dose Volume Histogram (DVH), ICI, chemotherapy, radiomics, dosiomics) and four combined models (radiomics + dosiomics, radiomics + DVH + Clinics, dosiomics + DVH + Clinics, radiomics + dosiomics + DVH + Clinics) were trained to predict PTP. Dose-based models were tested using physical dose and EQD2. Four ML-algorithms (random forest (rf), logistic elastic net regression, support vector machine, logitBoost) were trained and tested using 5-fold nested cross validation and Synthetic Minority Oversampling Technique (SMOTE) for resampling in R. Prediction was evaluated using the area under the receiver operating characteristic curve (AUC) on the test sets of the outer folds.

RESULTS

The combined model of all features using EQD2 surpassed all other models (AUC = 0.77, Confidence Interval CI 0.76-0.78). DVH, clinical data and ICI therapy had minor impact on PTP prediction with AUC values between 0.42 and 0.57. All EQD2-based models outperformed models based on physical dose.

CONCLUSIONS

Radiomics + dosiomics based ML models combined with clinical and dosimetric models were found to be suited best for PTP prediction after R(C)T and could improve pre-treatment decision making. Different RT dose fractionation schemes should be considered for dose-based ML approaches.

Advancing equitable and personalized cancer care: Novel applications and priorities of artificial intelligence for fairness and inclusivity in the patient care workflow

Patient care workflows are highly multimodal and intertwined: the intersection of data outputs provided from different disciplines and in different formats remains one of the main challenges of modern oncology. Artificial Intelligence (AI) has the potential to revolutionize the current clinical practice of oncology owing to advancements in digitalization, database expansion, computational technologies, and algorithmic innovations that facilitate discernment of complex relationships in multimodal data. Within oncology, radiation therapy (RT) represents an increasingly complex working procedure, involving many labor-intensive and operator-dependent tasks. In this context, AI has gained momentum as a powerful tool to standardize treatment performance and reduce inter-observer variability in a time-efficient manner. This review explores the hurdles associated with the development, implementation, and maintenance of AI platforms and highlights current measures in place to address them. In examining AI's role in oncology workflows, we underscore that a thorough and critical consideration of these challenges is the only way to ensure equitable and unbiased care delivery, ultimately serving patients' survival and quality of life.

Predicting radiation pneumonitis in lung cancer: a EUD-based machine learning approach for volumetric modulated arc therapy patients

Purpose

This study aims to develop an optimal machine learning model that uses lung equivalent uniform dose (lung EUD to predict radiation pneumonitis (RP) occurrence in lung cancer patients treated with volumetric modulated arc therapy (VMAT).

Methods

We analyzed a cohort of 77 patients diagnosed with locally advanced squamous cell lung cancer (LASCLC) receiving concurrent chemoradiotherapy with VMAT. Patients were categorized based on the onset of grade II or higher radiation pneumonitis (RP 2+). Dose volume histogram data, extracted from the treatment planning system, were used to compute the lung EUD values for both groups using a specialized numerical analysis code. We identified the parameter α, representing the most significant relative difference in lung EUD between the two groups. The predictive potential of variables for RP2+, including physical dose metrics, lung EUD, normal tissue complication probability (NTCP) from the Lyman-Kutcher-Burman (LKB) model, and lung EUD-calibrated NTCP for affected and whole lung, underwent both univariate and multivariate analyses. Relevant variables were then employed as inputs for machine learning models: multiple logistic regression (MLR), support vector machine (SVM), decision tree (DT), and K-nearest neighbor (KNN). Each model's performance was gauged using the area under the curve (AUC), determining the best-performing model.

Results

The optimal α-value for lung EUD was 0.3, maximizing the relative lung EUD difference between the RP 2+ and non-RP 2+ groups. A strong correlation coefficient of 0.929 (P< 0.01) was observed between lung EUD (α = 0.3) and physical dose metrics. When examining predictive capabilities, lung EUD-based NTCP for the affected lung (AUC: 0.862) and whole lung (AUC: 0.815) surpassed LKB-based NTCP for the respective lungs. The decision tree (DT) model using lung EUD-based predictors emerged as the superior model, achieving an AUC of 0.98 in both training and validation datasets.

Discussions

The likelihood of developing RP 2+ has shown a significant correlation with the advancements in RT technology. From traditional 3-D conformal RT, lung cancer treatment methodologies have transitioned to sophisticated techniques like static IMRT. Accurately deriving such a dose-effect relationship through NTCP modeling of RP incidence is statistically challenging due to the increased number of degrees-of-freedom. To the best of our knowledge, many studies have not clarified the rationale behind setting the α-value to 0.99 or 1, despite the closely aligned calculated lung EUD and lung mean dose MLD. Perfect independence among variables is rarely achievable in real-world scenarios. Four prominent machine learning algorithms were used to devise our prediction models. The inclusion of lung EUD-based factors substantially enhanced their predictive performance for RP 2+. Our results advocate for the decision tree model with lung EUD-based predictors as the optimal prediction tool for VMAT-treated lung cancer patients. Which could replace conventional dosimetric parameters, potentially simplifying complex neural network structures in prediction models.

The relationship of body mass index to setup errors, dosimetric parameters and incidence of radiation pneumonitis in non-small cell lung cancer patients undergoing intensity-modulated radiation therapy: a single-center observational study

BACKGROUND

The relationship among body mass index (BMI), setup error and radiation pneumonitis is not clearly illustrated.

OBJECTIVE

The present study aimed to investigate the role of BMI in non-small cell lung cancer (NSCLC) patients' radiation treatment, focusing on its relationship with setup error of patient positioning, the dosimetric parameters of intensity-modulated radiation therapy (IMRT) and the incidence of radiation pneumonitis.

METHODS

This prospective observational study included 523 cases of NSCLC patients during 2020-2022. Patients were divided into different groups by different BMI. The setup error was obtained by cone beam CT (CBCT) at three positions, lateral (LAT), longitudinal (LNG) and vertical (VRT). IMRT dosimetric parameters of V5, V20, and mean dose were collected.

RESULTS

Patients with BMI ≥28 kg/m2 showed significantly higher absolute values of LAT, LNG and VRT, higher V5, V20, mean dose, as well as higher total incidence of radiation pneumonitis and grade III radiation pneumonitis compared with patients with BMI <24 kg/m2 or 24-28 kg/m2. Spearman's analysis demonstrated that the absolute values of LAT, LNG and VRT were positively correlated with BMI, and positive correlation existed among BMI, dosimetric parameters and setup errors. ROC curves showed that LAT in setup errors and V5 in dosimetric parameters had the best diagnostic value for prediction of radiation pneumonitis. Only BMI, LAT, V5 and V20 were the independent risk factors for radiation pneumonitis.

CONCLUSIONS

Setup error caused by higher BMI might be associated with the dosimetric parameters, as well as the incidence of radiation pneumonitis in NSCLC patients.

Multi-omics and Multi-VOIs to predict esophageal fistula in esophageal cancer patients treated with radiotherapy

OBJECTIVE

This study aimed to develop a prediction model for esophageal fistula (EF) in esophageal cancer (EC) patients treated with intensity-modulated radiation therapy (IMRT), by integrating multi-omics features from multiple volumes of interest (VOIs).

METHODS

We retrospectively analyzed pretreatment planning computed tomographic (CT) images, three-dimensional dose distributions, and clinical factors of 287 EC patients. Nine groups of features from different combination of omics [Radiomics (R), Dosiomics (D), and RD (the combination of R and D)], and VOIs [esophagus (ESO), gross tumor volume (GTV), and EG (the combination of ESO and GTV)] were extracted and separately selected by unsupervised (analysis of variance (ANOVA) and Pearson correlation test) and supervised (Student T test) approaches. The final model performance was evaluated using five metrics: average area under the receiver-operator-characteristics curve (AUC), accuracy, precision, recall, and F1 score.

RESULTS

For multi-omics using RD features, the model performance in EG model shows: AUC, 0.817 ± 0.031; 95% CI 0.805, 0.825; p < 0.001, which is better than single VOI (ESO or GTV).

CONCLUSION

Integrating multi-omics features from multi-VOIs enables better prediction of EF in EC patients treated with IMRT. The incorporation of dosiomics features can enhance the model performance of the prediction.

Dosiomics for intensity-modulated radiotherapy in patients with prostate cancer: survival analysis stratified by baseline prostate-specific antigen and Gleason grade group in a 2-institutional retrospective study

OBJECTIVE

This study evaluated the prognostic impact of the quality of dose distribution using dosiomics in patients with prostate cancer, stratified by pretreatment prostate-specific antigen (PSA) levels and Gleason grade (GG) group.

METHODS

A total of 721 patients (Japanese Foundation for Cancer Research [JFCR] cohort: N = 489 and Tokyo Radiation Oncology Clinic [TROC] cohort: N = 232) with localized prostate cancer treated by intensity-modulated radiation therapy were enrolled. Two predictive dosiomic features for biochemical recurrence (BCR) were selected and patients were divided into certain groups stratified by pretreatment PSA levels and GG. Freedom from biochemical failure (FFBF) was estimated using the Kaplan-Meier method based on each dosiomic feature and univariate discrimination was evaluated using the log-rank test. As an exploratory analysis, a dosiomics hazard (DH) score was developed and its prognostic power for BCR was examined.

RESULTS

The dosiomic feature extracted from planning target volume (PTV) significantly distinguished the high- and low-risk groups in patients with PSA levels >10 ng/mL (7-year FFBF: 86.7% vs 76.1%, P < .01), GG 4 (92.2% vs 76.9%, P < .01), and GG 5 (83.1% vs 77.8%, P = .04). The DH score showed significant association with BCR (hazard score: 2.04; 95% confidence interval: 1.38-3.01; P < .001).

CONCLUSION

The quality of planned dose distribution on PTV may affect the prognosis of patients with poor prognostic factors, such as PSA levels >10 ng/mL and higher GGs.

ADVANCES IN KNOWLEDGE

The effects of planned dose distribution on prognosis differ depending on the patient's clinical background.

Development and validation of survival prognostic models for head and neck cancer patients using machine learning and dosiomics and CT radiomics features: a multicentric study

BACKGROUND

This study aimed to investigate the value of clinical, radiomic features extracted from gross tumor volumes (GTVs) delineated on CT images, dose distributions (Dosiomics), and fusion of CT and dose distributions to predict outcomes in head and neck cancer (HNC) patients.

METHODS

A cohort of 240 HNC patients from five different centers was obtained from The Cancer Imaging Archive. Seven strategies, including four non-fusion (Clinical, CT, Dose, DualCT-Dose), and three fusion algorithms (latent low-rank representation referred (LLRR),Wavelet, weighted least square (WLS)) were applied. The fusion algorithms were used to fuse the pre-treatment CT images and 3-dimensional dose maps. Overall, 215 radiomics and Dosiomics features were extracted from the GTVs, alongside with seven clinical features incorporated. Five feature selection (FS) methods in combination with six machine learning (ML) models were implemented. The performance of the models was quantified using the concordance index (CI) in one-center-leave-out 5-fold cross-validation for overall survival (OS) prediction considering the time-to-event.

RESULTS

The mean CI and Kaplan-Meier curves were used for further comparisons. The CoxBoost ML model using the Minimal Depth (MD) FS method and the glmnet model using the Variable hunting (VH) FS method showed the best performance with CI = 0.73 ± 0.15 for features extracted from LLRR fused images. In addition, both glmnet-Cindex and Coxph-Cindex classifiers achieved a CI of 0.72 ± 0.14 by employing the dose images (+ incorporated clinical features) only.

CONCLUSION

Our results demonstrated that clinical features, Dosiomics and fusion of dose and CT images by specific ML-FS models could predict the overall survival of HNC patients with acceptable accuracy. Besides, the performance of ML methods among the three different strategies was almost comparable.

Deep-Learning Model Prediction of Radiation Pneumonitis Using Pretreatment Chest Computed Tomography and Clinical Factors

Objectives: This study aimed to build a comprehensive deep-learning model for the prediction of radiation pneumonitis using chest computed tomography (CT), clinical, dosimetric, and laboratory data. Introduction: Radiation therapy is an effective tool for treating patients with lung cancer. Despite its effectiveness, the risk of radiation pneumonitis limits its application. Although several studies have demonstrated models to predict radiation pneumonitis, no reliable model has been developed yet. Herein, we developed prediction models using pretreatment chest CT and various clinical data to assess the likelihood of radiation pneumonitis in lung cancer patients. Methods: This retrospective study analyzed 3-dimensional (3D) lung volume data from chest CT scans and 27 features including dosimetric, clinical, and laboratory data from 548 patients who were treated at our institution between 2010 and 2021. We developed a neural network, named MergeNet, which processes lung 3D CT, clinical, dosimetric, and laboratory data. The MergeNet integrates a convolutional neural network with subsequent fully connected layers. A support vector machine (SVM) and light gradient boosting machine (LGBM) model were also implemented for comparison. For comparison, the convolution-only neural network was implemented as well. Three-dimensional Resnet-10 network and 4-fold cross-validation were used. Results: Classification performance was quantified by using the area under the receiver operative characteristic curve (AUC) metrics. MergeNet showed the AUC of 0.689. SVM, LGBM, and convolution-only networks showed AUCs of 0.525, 0.541, and 0.550, respectively. Application of DeLong test to pairs of receiver operating characteristic curves respectively yielded P values of .001 for the MergeNet-SVM pair and 0.001 for the MergeNet-LGBM pair. Conclusion: The MergeNet model, which incorporates chest CT, clinical, dosimetric, and laboratory data, demonstrated superior performance compared to other models. However, since its prediction performance has not yet reached an efficient level for clinical application, further research is required. Contribution: This study showed that MergeNet may be an effective means to predict radiation pneumonitis. Various predictive factors can be used together for the radiation pneumonitis prediction task via the MergeNet.

Integration of dosimetric parameters, clinical factors, and radiomics to predict symptomatic radiation pneumonitis in lung cancer patients undergoing combined immunotherapy and radiotherapy

PURPOSE

This study aimed to combine clinical/dosimetric factors and handcrafted/deep learning radiomic features to establish a predictive model for symptomatic (grade ≥ 2) radiation pneumonitis (RP) in lung cancer patients who received immunotherapy followed by radiotherapy.

MATERIALS AND METHODS

This study retrospectively collected data of 73 lung cancer patients with prior receipt of ICIs who underwent thoracic radiotherapy (TRT). Of these 73 patients, 41 (56.2 %) developed symptomatic grade ≥ 2 RP. RP was defined per multidisciplinary clinician consensus using CTCAE v5.0. Regions of interest (ROIs) (from radiotherapy planning CT images) utilized herein were gross tumor volume (GTV), planning tumor volume (PTV), and PTV-GTV. Clinical/dosimetric (mean lung dose and V5-V30) parameters were collected, and 107 handcrafted radiomic (HCR) features were extracted from each ROI. Deep learning-based radiomic (DLR) features were also extracted based on pre-trained 3D residual network models. HCR models, Fusion HCR model, Fusion HCR + ResNet models, and Fusion HCR + ResNet + Clinical models were built and compared using the receiver operating characteristic (ROC) curve with measurement of the area under the curve (AUC). Five-fold cross-validation was performed to avoid model overfitting.

RESULTS

HCR models across various ROIs and the Fusion HCR model showed good predictive ability with AUCs from 0.740 to 0.808 and 0.740-0.802 in the training and testing cohorts, respectively. The addition of DLR features improved the effectiveness of HCR models (AUCs from 0.826 to 0.898 and 0.821-0.898 in both respective cohorts). The best performing prediction model (HCR + ResNet + Clinical) combined HCR & DLR features with 7 clinical/dosimetric characteristics and achieved an average AUC of 0.936 and 0.946 in both respective cohorts.

CONCLUSIONS

In patients undergoing combined immunotherapy/RT for lung cancer, integrating clinical/dosimetric factors and handcrafted/deep learning radiomic features can offer a high predictive capacity for RP, and merits further prospective validation.

Dosiomics and radiomics improve the prediction of post-radiotherapy neutrophil-lymphocyte ratio in locally advanced non-small cell lung cancer

PURPOSE

To develop and validate a dosiomics and radiomics model based on three-dimensional (3D) dose distribution map and computed tomography (CT) images for the prediction of the post-radiotherapy (post-RT) neutrophil-to-lymphocyte ratio (NLR).

METHODS

This work retrospectively collected 242 locally advanced non-small cell lung cancer (LA-NSCLC) patients who were treated with definitive radiotherapy from 2012 to 2016. The NLR collected one month after the completion of RT was defined as the primary outcome. Clinical characteristics and two-dimensional dosimetric factors calculated from the dose-volume histogram (DVH) were included. A total of 4165 dosiomics and radiomics features were extracted from the 3D dose maps and CT images within five different anatomical regions of interest (ROIs), respectively. Then, a three-step feature selection method was proposed to progressively filter features from coarse to fine: (i) model-based ranking according to individual feature's performance, (ii) maximum relevance and minimum redundancy (mRMR), (iii) select from model based on feature importance calculated with an ensemble of several decision trees. The selected feature subsets were utilized to develop the prediction model with GBDT. All patients were divided into a development set and an independent testing set (2:1). Five-fold cross-validation was applied to the development set for both feature selection and model training procedure. Finally, a fusion model combining dosiomics, radiomics and clinical features was constructed to further improve the prediction results. The area under receiver operating characteristic curve (ROC) were used to evaluate the model performance.

RESULTS

The clinical-based and DVH-based models showed limited predictive power with AUCs of 0.632 (95% CI: 0.490-0.773) and 0.634 (95% CI: 0.497-0.771), respectively, in the independent testing set. The 9 feature-based dosiomics and 3 feature-based radiomics models showed improved AUCs of 0.738 (95% CI: 0.628-0.849) and 0.689 (95% CI: 0.566-0.813), respectively. The dosiomics & radiomics & clinical fusion model further improved the model's generalization ability with an AUC of 0.765 (95% CI: 0.656-0.874).

CONCLUSIONS

Dosiomics and radiomics can benefit the prediction of post-RT NLR of LA-NSCLC patients. This can provide a reference for evaluating radiotherapy-related inflammation.

Radiomics and dosiomics for predicting radiation-induced hypothyroidism and guiding intensity-modulated radiotherapy

To guide individualized intensity-modulated radiotherapy (IMRT), we developed and prospectively validated a multiview radiomics risk model for predicting radiation-induced hypothyroidism in patients with nasopharyngeal carcinoma. And simulated radiotherapy plans with same dose-volume-histogram (DVH) but different dose distributions were redesigned to explore the clinical application of the multiview radiomics risk model. The radiomics and dosiomics were built based on selected radiomics and dosiomics features from planning computed tomography and dose distribution, respectively. The multiview radiomics risk model that integrated radiomics, dosiomics, DVH parameters, and clinical factors had better performance than traditional normal tissue complication probability models. And multiview radiomics risk model could identify differences of patient hypothyroidism-free survival that cannot be stratified by traditional models. Besides, two redesigned simulated plans further verified the clinical application and advantage of the multiview radiomics risk model. The multiview radiomics risk model was a promising method to predict radiation-induced hypothyroidism and guide individualized IMRT.

Prediction of radiation pneumonia after radiotherapy for esophageal cancer using a unified fractional dosiomics combined model

OBJECTIVE

This study aimed to construct an optimal model to predict radiation pneumonia (RP) after radiotherapy for esophageal cancer using unified fractional dosiomics and to investigate the improvements in the prediction efficiency of each model for RP.

METHODS

The clinical data, DVH, pre-treatment CT, and dose distribution of 182 patients were retrospectively analyzed.The independent risk factors were screened using univariate and multivariate logistic regression. The mutual information (MI),least absolute shrinkage and selection operator (LASSO), and recursive feature elimination (RFE) methods were used to screen the omics features. The AUC values of ROC, calibration curves, and clinical decision curves were calculated to evaluate the efficacy and trends of each model.

RESULTS

The AUC of dosiomics model were 0.783 and 0.760 in the training and test cohorts, higher than 0.585 and 0.579 in the training and test cohorts of the DVH model. The AUC value of the R + D combination was the highest, reaching 0.833. The combined R + D model had a better calibration degree than the other models (mean absolute error = 0.018) and better net benefit in clinical decision-making.

CONCLUSIONS

The radiomics combined dosiomics model was the best combined model to predict RP after radiotherapy for esophageal cancer. The dosiomics model could cover the efficiency of the DVH model and significantly improve the efficiency of the combined model.In the future, we will include other centers for further verification.

ADVANCES IN KNOWLEDGE

For the first time, this study used CT images combined dose distribution to predict the occurrence of radiation pneumonitis after radiotherapy for esophageal cancer.

Dosiomic-based prediction of dysgeusia in head & neck cancer patients treated with radiotherapy

PURPOSE

To investigate the potential of dosiomics in predicting radiotherapy-induced taste distortion (dysgeusia) in head & neck (H&N) cancer.

METHODS

A cohort of 80 H&N cancer patients treated with radical or adjuvant radiotherapy and with a follow-up of at least 24 months was enrolled. Treatment information, as well as tobacco and alcohol consumption were also collected. The whole tongue was manually delineated on the planning CT and mapped to the dose map retrieved from the treatment planning system. For every patient, 6 regions of the tongue were examined; for each of them, 145 dosiomic features were extracted from the dose map and fed to a logistic regression model to predict the grade of dysgeusia at follow-up, with and without including clinical features. A mean dose-based model was considered for reference.

RESULTS

Both dosiomics and mean dose models achieved good prediction performance for acute dysgeusia with AUC up to 0.88. For the dosiomic model, the central and anterior ⅔ regions of the tongue were the most predictive. For all models, a gradual reduction in the performance was observed at later times for chronic dysgeusia prediction, with higher values for dosiomics. The inclusion of smoke and alcohol habits did not improve model performances.

CONCLUSION

The dosiomic analysis of the dose to the tongue identified features able to predict acute dysgeusia. Dosiomics resulted superior to the conventional mean dose-based model for chronic dysgeusia prediction. Larger, prospective studies are needed to support these results before integrating dosiomics in radiotherapy planning.

A radiomic- and dosiomic-based machine learning regression model for pretreatment planning in 177 Lu-DOTATATE therapy

BACKGROUND

Standardized patient-specific pretreatment dosimetry planning is mandatory in the modern era of nuclear molecular radiotherapy, which may eventually lead to improvements in the final therapeutic outcome. Only a comprehensive definition of a dosage therapeutic window encompassing the range of absorbed doses, that is, helpful without being detrimental can lead to therapy individualization and improved outcomes. As a result, setting absorbed dose safety limits for organs at risk (OARs) requires knowledge of the absorbed dose-effect relationship. Data sets of consistent and reliable inter-center dosimetry findings are required to characterize this relationship.

PURPOSE

We developed and standardized a new pretreatment planning model consisting of a predictive dosimetry procedure for OARs in patients with neuroendocrine tumors (NETs) treated with 177 Lu-DOTATATE (Lutathera). In the retrospective study described herein, we used machine learning (ML) regression algorithms to predict absorbed doses in OARs by exploiting a combination of radiomic and dosiomic features extracted from patients' imaging data.

METHODS

Pretreatment and posttreatment data for 20 patients with NETs treated with 177 Lu-DOTATATE were collected from two clinical centers. A total of 3412 radiomic and dosiomic features were extracted from the patients' computed tomography (CT) scans and dose maps, respectively. All dose maps were generated using Monte Carlo simulations. An ML regression model was designed based on ML algorithms for predicting the absorbed dose in every OAR (liver, left kidney, right kidney, and spleen) before and after the therapy and between each therapy session, thus predicting any possible radiotoxic effects.

RESULTS

We evaluated nine ML regression algorithms. Our predictive model achieved a mean absolute dose error (MAE, in Gy) of 0.61 for the liver, 1.58 for the spleen, 1.30 for the left kidney, and 1.35 for the right kidney between pretherapy 68 Ga-DOTATOC positron emission tomography (PET)/CT and posttherapy 177 Lu-DOTATATE single photon emission (SPECT)/CT scans. Τhe best predictive performance observed was based on the gradient boost for the liver, the left kidney and the right kidney, and on the extra tree regressor for the spleen. Evaluation of the model's performance according to its ability to predict the absorbed dose in each OAR in every possible combination of pretherapy 68 Ga-DOTATOC PET/CT and any posttherapy 177 Lu-DOTATATE treatment cycle SPECT/CT scans as well as any 177 Lu-DOTATATE SPECT/CT treatment cycle and the consequent 177 Lu-DOTATATE SPECT/CT treatment cycle revealed mean absorbed dose differences ranges from -0.55 to 0.68 Gy. Incorporating radiodosiomics features from the 68 Ga-DOTATOC PET/CT and first 177 Lu-DOTATATE SPECT/CT treatment cycle scans further improved the precision and minimized the standard deviation of the predictions in nine out of 12 instances. An average improvement of 57.34% was observed (range: 17.53%-96.12%). However, it's important to note that in three instances (i.e., Ga,C.1 → C3 in spleen and left kidney, and Ga,C.1 → C2 in right kidney) we did not observe an improvement (absolute differences of 0.17, 0.08, and 0.05 Gy, respectively). Wavelet-based features proved to have high correlated predictive value, whereas non-linear-based ML regression algorithms proved to be more capable than the linear-based of producing precise prediction in our case.

CONCLUSIONS

The combination of radiomics and dosiomics has potential utility for personalized molecular radiotherapy (PMR) response evaluation and OAR dose prediction. These radiodosiomic features can potentially provide information on any possible disease recurrence and may be highly useful in clinical decision-making, especially regarding dose escalation issues.

CT-based dosiomics and radiomics model predicts radiation-induced lymphopenia in nasopharyngeal carcinoma patients

Purpose

This study aims to develop and validate a model predictive for the incidence of grade 4 radiation-induced lymphopenia (G4RIL), based on dosiomics features and radiomics features from the planning CT of nasopharyngeal carcinoma (NPC) treated by radiation therapy.

Methods

The dataset of 125 NPC patients treated with radiotherapy from August 2018 to March 2019 was randomly divided into two sets-an 85-sample training set and a 40-sample test set. Dosiomics features and radiomics features of the CT image within the skull bone and cervical vertebrae were extracted. A feature selection process of multiple steps was employed to identify the features that most accurately forecast the data and eliminate superfluous or insignificant ones. A support vector machine learning classifier with correction for imbalanced data was trained on the patient dataset for prediction of RIL (positive classifier for G4RIL, negative otherwise). The model's predictive capability was gauged by gauging its sensitivity (the likelihood of a positive test being administered to patients with G4RIL) and specificity in the test set. The area beneath the ROC curve (AUC) was utilized to explore the association of characteristics with the occurrence of G4RIL.

Results

Three clinical features, three dosiomics features, and three radiomics features exhibited significant correlations with G4RIL. Those features were then used for model construction. The combination model, based on nine robust features, yielded the most impressive results with an ACC value of 0.88 in the test set, while the dosiomics model, with three dosiomics features, had an ACC value of 0.82, the radiomics model, with three radiomics features, had an ACC value of 0.82, and the clinical model, with its initial features, had an ACC value of 0.6 for prediction performance.

Conclusion

The findings show that radiomics and dosiomics features are correlated with the G4RIL of NPC patients. The model incorporating radiomics features and dosiomics features from planning CT can predict the incidence of G4RIL in NPC patients.

Radiation pneumonia predictive model for radiotherapy in esophageal carcinoma patients

BACKGROUND

The machine learning models with dose factors and the deep learning models with dose distribution matrix have been used to building lung toxics models for radiotherapy and achieve promising results. However, few studies have integrated clinical features into deep learning models. This study aimed to explore the role of three-dimension dose distribution and clinical features in predicting radiation pneumonitis (RP) in esophageal cancer patients after radiotherapy and designed a new hybrid deep learning network to predict the incidence of RP.

METHODS

A total of 105 esophageal cancer patients previously treated with radiotherapy were enrolled in this study. The three-dimension (3D) dose distributions within the lung were extracted from the treatment planning system, converted into 3D matrixes and used as inputs to predict RP with ResNet. In total, 15 clinical factors were normalized and converted into one-dimension (1D) matrixes. A new prediction model (HybridNet) was then built based on a hybrid deep learning network, which combined 3D ResNet18 and 1D convolution layers. Machine learning-based prediction models, which use the traditional dosiomic factors with and without the clinical factors as inputs, were also constructed and their predictive performance compared with that of HybridNet using tenfold cross validation. Accuracy and area under the receiver operator characteristic curve (AUC) were used to evaluate the model effect. DeLong test was used to compare the prediction results of the models.

RESULTS

The deep learning-based model achieved superior prediction results compared with machine learning-based models. ResNet performed best in the group that only considered dose factors (accuracy, 0.78 ± 0.05; AUC, 0.82 ± 0.25), whereas HybridNet performed best in the group that considered both dose factors and clinical factors (accuracy, 0.85 ± 0.13; AUC, 0.91 ± 0.09). HybridNet had higher accuracy than that of Resnet (p = 0.009).

CONCLUSION

Based on prediction results, the proposed HybridNet model could predict RP in esophageal cancer patients after radiotherapy with significantly higher accuracy, suggesting its potential as a useful tool for clinical decision-making. This study demonstrated that the information in dose distribution is worth further exploration, and combining multiple types of features contributes to predict radiotherapy response.

Enhanced prediction of postoperative radiotherapy-induced esophagitis in non-small cell lung cancer: Dosiomic model development in a real-world cohort and validation in the PORT-C randomized controlled trial

BACKGROUND

Radiotherapy-induced esophagitis (RE) diminishes the quality of life and interrupts treatment in patients with non-small cell lung cancer (NSCLC) undergoing postoperative radiotherapy. Dosimetric models showed limited capability in predicting RE. We aimed to develop dosiomic models to predict RE.

METHODS

Models were trained with a real-world cohort and validated with PORT-C randomized controlled trial cohort. Patients with NSCLC undergoing resection followed by postoperative radiotherapy between 2004 and 2015 were enrolled. The endpoint was grade ≥2 RE. Esophageal three-dimensional dose distribution features were extracted using handcrafted and convolutional neural network (CNN) methods, screened using an entropy-based method, and selected using minimum redundancy and maximum relevance. Prediction models were built using logistic regression. The areas under the receiver operating characteristic curve (AUC) and precision-recall curve were used to evaluate prediction model performance. A dosimetric model was built for comparison.

RESULTS

A total of 190 and 103 patients were enrolled in the training and validation sets, respectively. Using handcrafted and CNN methods, 107 and 4096 features were derived, respectively. Three handcrafted, four CNN-extracted and three dosimetric features were selected. AUCs of training and validation sets were 0.737 and 0.655 for the dosimetric features, 0.730 and 0.724 for handcrafted features, and 0.812 and 0.785 for CNN-extracted features, respectively. Precision-recall curves revealed that CNN-extracted features outperformed dosimetric and handcrafted features.

CONCLUSIONS

Prediction models may identify patients at high risk of developing RE. Dosiomic models outperformed the dosimetric-feature model in predicting RE. CNN-extracted features were more predictive but less interpretable than handcrafted features.

Analysis of the risk factors of radiation pneumonitis and the predictive ability of dosiomics in non-small-cell lung cancer

Purpose: This prospective study investigated the incidence of radiation pneumonitis (RP) after immunotherapy followed by radiotherapy in non-small-cell lung cancer, analyzed the risk factors for RP, and explored the predictive performance of dosimetry and dosiomics. Methods & materials: Risk factors for grade ≥2 RP were calculated by using a logistic regression model. Predictive performance was compared on the basis of area under the curve values. Results: Grade ≥2 RP occurred in 16 cases (26.7%). The AUC values of V5 Gy, gray-level dependence matrix-small dependence high gray-level emphasis (GLDM-SDHGLE) and combined features were 0.685, 0.724 and 0.734, respectively. Conclusion: Smoking history, bilateral lung V5 Gy and GLDM-SDHGLE were independent risk factors for RP. Dosiomics can effectively predict RP.

A multiomics approach-based prediction of radiation pneumonia in lung cancer patients: impact on survival outcome

PURPOSE

To predict the risk of radiation pneumonitis (RP), a multiomics model was built to stratify lung cancer patients. Our study also investigated the impact of RP on survival.

METHODS

This study retrospectively collected 100 RP and 99 matched non-RP lung cancer patients treated with radiotherapy from two independent centres. They were divided into training (n = 175) and validation cohorts (n = 24). The radiomics, dosiomics and clinical features were extracted from planning CT and electronic medical records and were analysed by LASSO Cox regression. A multiomics prediction model was developed by the optimal algorithm. Overall survival (OS) between the RP, non-RP, mild RP, and severe RP groups was analysed by the Kaplan‒Meier method.

RESULTS

Sixteen radiomics features, two dosiomics features, and one clinical feature were selected to build the best multiomics model. The optimal performance for predicting RP was the area under the receiver operating characteristic curve (AUC) of the testing set (0.94) and validation set (0.92). The RP patients were divided into mild (≤ 2 grade) and severe (> 2 grade) RP groups. The median OS was 31 months for the non-RP group compared with 49 months for the RP group (HR = 0.53, p = 0.0022). Among the RP subgroup, the median OS was 57 months for the mild RP group and 25 months for the severe RP group (HR = 3.72, p < 0.0001).

CONCLUSIONS

The multiomics model contributed to improving the accuracy of RP prediction. Compared with the non-RP patients, the RP patients displayed longer OS, especially the mild RP patients.

Prediction of toxicity outcomes following radiotherapy using deep learning-based models: A systematic review

PURPOSE

This study aims to perform a comprehensive systematic review of deep learning (DL) models in predicting RT-induced toxicity.

MATERIALS AND METHODS

A literature review was performed following the PRISMA guidelines. Systematic searches were performed in PubMed, Scopus, Cochrane and Embase databases from the earliest record up to September 2022. Related studies on deep learning models for radiotherapy toxicity prediction were selected based on predefined PICOS criteria.

RESULTS

Fourteen studies of radiotherapy-treated patients on different types of cancer [prostate (n=2), HNC (n=4), liver (n=2), lung (n=4), cervical (n=1), and oesophagus (n=1)] were eligible for inclusion in the systematic review. Information regarding patient characteristics and model development was summarized. Several approaches, such as ensemble learning, data augmentation, and transfer learning, that were utilized by selected studies were discussed.

CONCLUSION

Deep learning techniques are able to produce a consistent performance for toxicity prediction. Future research using large and diverse datasets and standardization of the study methodologies are required to improve the consistency of the research output.

Application of error classification model using indices based on dose distribution for characteristics evaluation of multileaf collimator position errors

This study aims to evaluate the specific characteristics of various multileaf collimator (MLC) position errors that are correlated with the indices using dose distribution. The dose distribution was investigated using the gamma, structural similarity, and dosiomics indices. Cases from the American Association of Physicists in Medicine Task Group 119 were planned, and systematic and random MLC position errors were simulated. The indices were obtained from distribution maps and statistically significant indices were selected. The final model was determined when all values of the area under the curve, accuracy, precision, sensitivity, and specificity were higher than 0.8 (p < 0.05). The dose-volume histogram (DVH) relative percentage difference between the error-free and error datasets was examined to investigate clinical relations. Seven multivariate predictive models were finalized. The common significant dosiomics indices (GLCM Energy and GLRLM_LRHGE) can characterize the MLC position error. In addition, the finalized logistic regression model for MLC position error prediction showed excellent performance with AUC > 0.9. Furthermore, the results of the DVH were related to dosiomics analysis in that it reflects the characteristics of the MLC position error. It was also shown that dosiomics analysis could provide important information on localized dose-distribution differences in addition to DVH information.

Dosiomics-Based Prediction of Radiation-Induced Valvulopathy after Childhood Cancer

Valvular Heart Disease (VHD) is a known late complication of radiotherapy for childhood cancer (CC), and identifying high-risk survivors correctly remains a challenge. This paper focuses on the distribution of the radiation dose absorbed by heart tissues. We propose that a dosiomics signature could provide insight into the spatial characteristics of the heart dose associated with a VHD, beyond the already-established risk induced by high doses. We analyzed data from the 7670 survivors of the French Childhood Cancer Survivors' Study (FCCSS), 3902 of whom were treated with radiotherapy. In all, 63 (1.6%) survivors that had been treated with radiotherapy experienced a VHD, and 57 of them had heterogeneous heart doses. From the heart-dose distribution of each survivor, we extracted 93 first-order and spatial dosiomics features. We trained random forest algorithms adapted for imbalanced classification and evaluated their predictive performance compared to the performance of standard mean heart dose (MHD)-based models. Sensitivity analyses were also conducted for sub-populations of survivors with spatially heterogeneous heart doses. Our results suggest that MHD and dosiomics-based models performed equally well globally in our cohort and that, when considering the sub-population having received a spatially heterogeneous dose distribution, the predictive capability of the models is significantly improved by the use of the dosiomics features. If these findings are further validated, the dosiomics signature may be incorporated into machine learning algorithms for radiation-induced VHD risk assessment and, in turn, into the personalized refinement of follow-up guidelines.

Radiomics and dosiomics-based prediction of radiotherapy-induced xerostomia in head and neck cancer patients

BACKGROUND AND AIM

Dose-response modeling for radiotherapy-induced xerostomia in head and neck cancer (HN) patients is a promising frontier for personalized therapy. Feature extraction from diagnostic and therapeutic images (radiomics and dosiomics features) can be used for data-driven response modeling. The aim of this study is to develop xerostomia predictive models based on radiomics-dosiomics features.

METHODS

Data from the cancer imaging archive (TCIA) for 31 HN cancer patients were employed. For all patients, parotid CT radiomics features were extracted, utilizing Lasso regression for feature selection and multivariate modeling. The models were developed by selected features from pretreatment (CT1), mid-treatment (CT2), post-treatment (CT3), and delta features (ΔCT2-1, ΔCT3-1, ΔCT3-2). We also considered dosiomics features extracted from the parotid dose distribution images (Dose model). Thus, combination models of radio-dosiomics (CT + dose & ΔCT + dose) were developed. Moreover, clinical, and dose-volume histogram (DVH) models were built. Nested 10-fold cross-validation was used to assess the predictive classification of patients into those with and without xerostomia, and the area under the receiver operative characteristic curve (AUC) was used to compare the predictive power of the models. The sensitivity and accuracy of models also were obtained.

RESULTS

In total, 59 parotids were assessed, and 13 models were developed. Our results showed three models with AUC of 0.89 as most predictive, namely ΔCT2-1 + Dose (Sensitivity 0.99, Accuracy 0.94 & Specificity 0.86), CT3 model (Sensitivity 0.96, Accuracy 0.94 & Specificity 0.86) and DVH (Sensitivity 0.93, Accuracy 0.89 & Specificity 0.84). These models were followed by Clinical (AUC 0.89, Sensitivity 0.81, Accuracy 0.97 & Specificity 0.89) and CT2 & Dose (AUC 0.86, Sensitivity 0.97, Accuracy 0.87 & Specificity 0.82). The Dose model (developed by dosiomics features only) had AUC, Sensitivity, Specificity, and Accuracy of 0.72, 0.98, 0.33, and 0.79 respectively.

CONCLUSION

Quantitative features extracted from diagnostic imaging during and after radiotherapy alone or in combination with dosiomics markers obtained from dose distribution images can be used for radiotherapy response modeling, opening up prospects for personalization of therapies toward improved therapeutic outcomes.

Dosiomics Risk Model for Predicting Radiation Induced Temporal Lobe Injury and Guiding Individual Intensity-Modulated Radiation Therapy

PURPOSE

We aimed to assess the value of dose distribution-based dosiomics and planning computed tomography-based radiomics to predict radiation-induced temporal lobe injury (TLI) and guide individualized intensity modulated radiation therapy.

METHODS AND MATERIALS

A total of 5599 nasopharyngeal carcinoma patients were enrolled, including 2503, 1072, 988, and 1036 patients in the training, validation, prospective test, and external test cohorts, respectively. The concordance index (C-index) was used to compare the performance of the radiomics and dosiomics models with that of the quantitative analyses of normal tissue effects in the clinic and Wen's models. The predicted TLI-free survival rates of redesigned simulated plans with the same dose-volume histogram but different dose distributions for same patient in a cohort of 30 randomly selected patients were compared by the Wilcoxon matched-pairs signed-rank test.

RESULTS

The radiomics and dosiomics signatures were constructed based on 30 selected computed tomography features and 10 selected dose distribution features, respectively, which were important predictors of TLI-free survival (all P <.001). However, the radiomics signature had a low C-index. The dosiomics risk model combining the dosiomics signature, D_1cc, and age had favorable performance, with C-index values of 0.776, 0.811, 0.805, and 0.794 in the training, validation, prospective test, and external test cohorts, respectively, which were better than those of the quantitative analyses of normal tissue effects in the clinic model and Wen's model (all P <.001). The dosiomics risk model can further distinguish patients in a same risk category divided by other models (all P <.05). Conversely, the other models were unable to separate populations classified by the dosiomics risk model (all P > .05). Two simulated plans with the same dose-volume histogram but different dose distributions had different TLI-free survival rates predicted by dosiomics risk model (all P ≤ .002).

CONCLUSIONS

The dosiomics risk model was superior to traditional models in predicting the risk of TLI. This is a promising approach to precisely predict radiation-induced toxicities and guide individualized intensity modulated radiation therapy.

Multi-omics to predict acute radiation esophagitis in patients with lung cancer treated with intensity-modulated radiation therapy

PURPOSE

The study aimed to predict acute radiation esophagitis (ARE) with grade  ≥ 2 for patients with locally advanced lung cancer (LALC) treated with intensity-modulated radiation therapy (IMRT) using multi-omics features, including radiomics and dosiomics.

METHODS

161 patients with stage IIIA-IIIB LALC who received chemoradiotherapy (CRT) or radiotherapy by IMRT with a prescribed dose from 45 to 70 Gy from 2015 to 2019 were enrolled retrospectively. All the toxicity gradings were given following the Common Terminology Criteria for Adverse Events V4.0. Multi-omics features, including radiomics, dosiomics (including dose-volume histogram dosimetric parameters), were extracted based on the planning CT image and three-dimensional dose distribution. All data were randomly divided into training cohorts (N = 107) and testing cohorts (N = 54). In the training cohorts, features with reliably high outcome relevance and low redundancy were selected under random patient subsampling. Four classification models (using clinical factors (CF) only, using radiomics features (RFs) only, dosiomics features (DFs) only, and the hybrid features (HFs) containing clinical factors, radiomics and dosiomics) were constructed employing the Ridge classifier using two-thirds of randomly selected patients as the training cohort. The remaining patient was treated as the testing cohort. A series of models were built with 30 times training-testing splits. Their performances were assessed using the area under the ROC curve (AUC) and accuracy.

RESULTS

Among all patients, 51 developed ARE grade  ≥ 2, with an incidence of 31.7%. Next, 8990 radiomics and 213 dosiomics features were extracted, and 3, 6, 12, and 13 features remained after feature selection in the CF, DF, RF and DF models, respectively. The RF and HF models achieved similar classification performance, with the training and testing AUCs of 0.796 ± 0.023 (95% confidence interval (CI [0.79, 0.80])/0.744 ± 0.044 (95% CI [0.73, 0.76]) and 0.801 ± 0.022 (95% CI [0.79, 0.81]) (p = 0.74), respectively. The model performances using CF and DF features were poorer, with training and testing AUCs of 0.573 ± 0.026 (95% CI [0.56, 0.58])/ 0.509 ± 0.072 (95% CI [0.48, 0.53]) and 0.679 ± 0.027 (95% CI [0.67, 0.69])/0.604 ± 0.041 (95% CI [0.53, 0.63]) compared with the above two models (p < 0.001), respectively.

CONCLUSIONS

In LALC patients treated with CRT IMRT, the ARE grade  ≥ 2 can be predicted using the pretreatment radiotherapy image features. To predict ARE, the multi-omics features had similar predictability with radiomics features; however, the dosiomics features and clinical factors had a limited classification performance.

Machine learning in non-small cell lung cancer radiotherapy: A bibliometric analysis

Background

Machine learning is now well-developed in non-small cell lung cancer (NSCLC) radiotherapy. But the research trend and hotspots are still unclear. To investigate the progress in machine learning in radiotherapy NSCLC, we performed a bibliometric analysis of associated research and discuss the current research hotspots and potential hot areas in the future.

Methods

The involved researches were obtained from the Web of Science Core Collection database (WoSCC). We used R-studio software, the Bibliometrix package and VOSviewer (Version 1.6.18) software to perform bibliometric analysis.

Results

We found 197 publications about machine learning in radiotherapy for NSCLC in the WoSCC, and the journal Medical Physics contributed the most articles. The University of Texas MD Anderson Cancer Center was the most frequent publishing institution, and the United States contributed most of the publications. In our bibliometric analysis, “radiomics” was the most frequent keyword, and we found that machine learning is mainly applied to analyze medical images in the radiotherapy of NSCLC.

Results

The research we identified about machine learning in NSCLC radiotherapy was mainly related to the radiotherapy planning of NSCLC and the prediction of treatment effects and adverse events in NSCLC patients who were under radiotherapy. Our research has added new insights into machine learning in NSCLC radiotherapy and could help researchers better identify hot research areas in the future.

Dosiomics and radiomics to predict pneumonitis after thoracic stereotactic body radiotherapy and immune checkpoint inhibition

Introduction

Pneumonitis is a relevant side effect after radiotherapy (RT) and immunotherapy with checkpoint inhibitors (ICIs). Since the effect is radiation dose dependent, the risk increases for high fractional doses as applied for stereotactic body radiation therapy (SBRT) and might even be enhanced for the combination of SBRT with ICI therapy. Hence, patient individual pre-treatment prediction of post-treatment pneumonitis (PTP) might be able to support clinical decision making. Dosimetric factors, however, use limited information and, thus, cannot exploit the full potential of pneumonitis prediction.

Methods

We investigated dosiomics and radiomics model based approaches for PTP prediction after thoracic SBRT with and without ICI therapy. To overcome potential influences of different fractionation schemes, we converted physical doses to 2 Gy equivalent doses (EQD2) and compared both results. In total, four single feature models (dosiomics, radiomics, dosimetric, clinical factors) were tested and five combinations of those (dosimetric+clinical factors, dosiomics+radiomics, dosiomics+dosimetric+clinical factors, radiomics+dosimetric+clinical factors, radiomics+dosiomics+dosimetric+clinical factors). After feature extraction, a feature reduction was performed using pearson intercorrelation coefficient and the Boruta algorithm within 1000-fold bootstrapping runs. Four different machine learning models and the combination of those were trained and tested within 100 iterations of 5-fold nested cross validation.

Results

Results were analysed using the area under the receiver operating characteristic curve (AUC). We found the combination of dosiomics and radiomics features to outperform all other models with AUCradiomics+dosiomics, D = 0.79 (95% confidence interval 0.78-0.80) and AUCradiomics+dosiomics, EQD2 = 0.77 (0.76-0.78) for physical dose and EQD2, respectively. ICI therapy did not impact the prediction result (AUC ≤ 0.5). Clinical and dosimetric features for the total lung did not improve the prediction outcome.

Conclusion

Our results suggest that combined dosiomics and radiomics analysis can improve PTP prediction in patients treated with lung SBRT. We conclude that pre-treatment prediction could support clinical decision making on an individual patient basis with or without ICI therapy.

3D deep convolution neural network for radiation pneumonitis prediction following stereotactic body radiotherapy

In this study, we investigated 3D convolutional neural networks (CNNs) with input from radiographic and dosimetric datasets of primary lung tumors and surrounding lung volumes to predict the likelihood of radiation pneumonitis (RP). Pre-treatment, 3- and 6-month follow-up computed tomography (CT) and 3D dose datasets from one hundred and ninety-three NSCLC patients treated with stereotactic body radiotherapy (SBRT) were retrospectively collected and analyzed for this study. DenseNet-121 and ResNet-50 models were selected for this study as they are deep neural networks and have been proven to have high accuracy for complex image classification tasks. Both were modified with 3D convolution and max pooling layers to accept 3D datasets. We used a minority class oversampling approach and data augmentation to address the challenges of data imbalance and data scarcity. We built two sets of models for classification of three (No RP, Grade 1 RP, Grade 2 RP) and two (No RP, Yes RP) classes as outputs. The 3D DenseNet-121 models performed better (F1 score [0.81], AUC [0.91] [three class]; F1 score [0.77], AUC [0.84] [two class]) than the 3D ResNet-50 models (F1 score [0.54], AUC [0.72] [three-class]; F1 score [0.68], AUC [0.71] [two-class]) (p = 0.017 for three class predictions). We also attempted to identify salient regions within the input 3D image dataset via integrated gradient (IG) techniques to assess the relevance of the tumor surrounding volume for RP stratification. These techniques appeared to indicate the significance of the tumor and surrounding regions in the prediction of RP. Overall, 3D CNNs performed well to predict clinical RP in our cohort based on the provided image sets and radiotherapy dose information.