- 1Department of Radiation Oncology, The Affiliated Cancer Hospital of Zhengzhou University and Henan Cancer Hospital, Zhengzhou, China
- 2Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong SAR, China
Purpose: This study investigates the impact of lung function on radiation pneumonitis prediction using a dual-omics analysis method.
Methods: We retrospectively collected data of 126 stage III lung cancer patients treated with chemo-radiotherapy using intensity-modulated radiotherapy, including pre-treatment planning CT images, radiotherapy dose distribution, and contours of organs and structures. Lung perfusion functional images were generated using a previously developed deep learning method. The whole lung (WL) volume was divided into function-wise lung (FWL) regions based on the lung perfusion functional images. A total of 5,474 radiomics features and 213 dose features (including dosiomics features and dose-volume histogram factors) were extracted from the FWL and WL regions, respectively. The radiomics features (R), dose features (D), and combined dual-omics features (RD) were used for the analysis in each lung region of WL and FWL, labeled as WL-R, WL-D, WL-RD, FWL-R, FWL-D, and FWL-RD. The feature selection was carried out using ANOVA, followed by a statistical F-test and Pearson correlation test. Thirty times train-test splits were used to evaluate the predictability of each group. The overall average area under the receiver operating characteristic curve (AUC), accuracy, precision, recall, and f1-score were calculated to assess the performance of each group.
Results: The FWL-RD achieved a significantly higher average AUC than the WL-RD group in the training (FWL-RD: 0.927 ± 0.031, WL-RD: 0.849 ± 0.064) and testing cohorts (FWL-RD: 0.885 ± 0.028, WL-RD: 0.762 ± 0.053, p < 0.001). When using radiomics features only, the FWL-R group yielded a better classification result than the model trained with WL-R features in the training (FWL-R: 0.919 ± 0.036, WL-R: 0.820 ± 0.052) and testing cohorts (FWL-R: 0.862 ± 0.028, WL-R: 0.750 ± 0.057, p < 0.001). The FWL-D group obtained an average AUC of 0.782 ± 0.032, obtaining a better classification performance than the WL-D feature-based model of 0.740 ± 0.028 in the training cohort, while no significant difference was observed in the testing cohort (FWL-D: 0.725 ± 0.064, WL-D: 0.710 ± 0.068, p = 0.54).
Conclusion: The dual-omics features from different lung functional regions can improve the prediction of radiation pneumonitis for lung cancer patients under IMRT treatment. This function-wise dual-omics analysis method holds great promise to improve the prediction of radiation pneumonitis for lung cancer patients.
Introduction
Lung cancer is the leading cause of cancer-related death worldwide (Sung et al., 2021). Radiation therapy or radiotherapy (RT) is one of the golden-standard treatment techniques for patients with locally advanced non-small-cell lung cancer (NSCLC) (Kong et al., 2005; Chang et al., 2016). Study shows a higher radiation dose can achieve better tumor control and improve the treatment outcome (Kong et al., 2005). However, dose escalation of lung cancer is greatly limited by radiation-induced side effects, such as radiation pneumonitis (RP). RP may occur in up to 30% of lung RT patients and is lethal in 2% of them (Zhang et al., 2012; Kipritidis et al., 2015). Hence, predicting RP is highly desirable for better dose optimization and personalization to maximize the treatment outcome in lung cancer RT.
At present treatment planning of lung cancer RT, several dosimetric factors from the dose-volume histogram (DVH) were found to be associated with RP, such as V5, V20, and Dmean (Baisden et al., 2007; Barriger et al., 2012; Bongers et al., 2013; Palma et al., 2013; Cai et al., 2014; Pinnix et al., 2015). These parameters are commonly used as dose constraints in clinical plan evaluation (Ganti et al., 2021). Meanwhile, several prevalent models using DVH parameters, such as normal tissue complication probability (NTCP), were proposed to predict high risk RP patients (Begosh-Mayne et al., 2020; Wang et al., 2020). However, DVH parameters can only distinguish statistical one-dimensional dose information rather than characterizing the dose distribution heterogeneity. With the aid of the radiomics definition (Lambin et al., 2017), dosiomics features were calculated based on the three-dimensional dose distribution to describe the dose spatial information (Liang et al., 2019). Several studies also have demonstrated significantly superior models with the dosiomics feature compared to the DVH-based model or the NTCP model for predicting RP (Liang et al., 2019; Palma et al., 2019; Adachi et al., 2021). Meanwhile, CT-based radiomics features describe the statistical information, shaped, and textual characteristics in a certain volume. The dual-omics combines the radiomics and dose features and is able to further improves the prediction for RP (Adachi et al., 2021; Jiang et al., 2021; Puttanawarut et al., 2022). However, those radiomics or dose features utilized in current studies were calculated from the whole lung region, rather than considering the heterogeneity inside the lung, for example, the difference in high- and low- functional lung regions.
Lung function information has been proven to be associated with RP, which promises to improve the RP prediction accuracy (Bucknell et al., 2018; Lee et al., 2018; Weller et al., 2019; Bourbonne et al., 2020; O’Reilly et al., 2020). O’Reilly et al. demonstrated the RP prediction improvement using the DVH factor (V20) from three high functional lung regions and compared these biomarkers to the entire lung region (O’Reilly et al., 2020). Lee et al. evaluated the correlation between several DVH factors (V5, V20, and Dmean) between the high functional lung region and the whole lung region, showing the potential of stratifying patients for pneumonitis prediction (Lee et al., 2018). Owen et al. demonstrated that irradiating to low functional lung region may increase radiation toxicity (Owen et al., 2021). Several studies also showed the potential of using dosimetry parameters based on functional lung images in predicting RP (Wang et al., 2012; Farr et al., 2015; Kimura et al., 2015; Xiao et al., 2018; Owen et al., 2021). However, these studies only explored the association between the dose factors and the RP without investigating the correlation between anatomical CT images and the RP. Besides, most studies focused on the dose features in the high functional lung region rather than the low functional lung region.
In this study, we developed a function-wise lung (FWL) analysis approach by integrating radiomics and dose features from both whole lung (WL) and FWL (including separated high- and low- functional lung) to predict RP for NSCLC patients. The radiomics features, dose features, and combined dual-omics features of each group were utilized for analysis. The feature selection metrics are the ANOVA followed by the statistical F-test and the Pearson correlation test. Thirty times train-test splits were used to evaluate the predictability of each group. The overall average area under the receiver operating characteristic curve (AUC), accuracy, precision, recall, and f1-score were calculated to assess the performance of each group.
Materials and methods
Data characteristics
The inclusion criteria are as follows: 1) diagnosed as primary locally-advanced lung cancer (stages IIIA/IIIC (AJCC 8th)); 2) having no distant metastasis; 3) treated with curative intensity-modulated radiotherapy (IMRT); 4) receiving contrast-enhanced CT for RT; 5) 18–70 years old. And the exclusion criteria are as follows:1) received chest RT or surgery or chemotherapy previously; 2) having previous chest malignancies; 3) received RT < 2 weeks; 4) incomplete RT treatment due to factors other than acute RP; 5) incomplete RT data.
The study was approved by the Institutional Review Board of the Affiliated Cancer Hospital of Zhengzhou University. Initially, a total of 162 pathological confirmed NSCLC patients staging IIIA/IIIC between 2015 and 2019 were retrospectively collected from the hospital. Considering the excluded criteria (shown in below), 126 cases were final enrolled in the study (Figure 1). All patients were treated by the 6 MV IMRT with a 50–70 Gy total prescription dose and 1.8–2.2 Gy fractional dose for 5 days per week. The radiation pneumonitis (RP) case was consecutively followed up at least 6 months after the first radiotherapy, and then graded with the Common Terminology Criteria for Adverse Events (CTCAE) V4.0. by one qualified imaging physician based on the CT scans. In this study, RP patients with grading
Image acquisition
Three types of image data were involved in this study, including planning CT images, three-dimension dose distribution images, and organs-at-risk (OAR) structures. All planning CT images were acquired from a 16-slice Brilliance Big Bore CT (Philips Medical System, Cleveland, OH, U.S.). The scanning parameters were as follows: scanning X-ray tube voltage = 120 kV, current = 321 mA, thickness = 3 mm, slice pixels = 512 × 512 and spacing = 1.152 mm × 1.152 mm. The scanning range was from the level of the cricoid cartilage to the lower border of the 12th thoracic vertebra covering the WL volume (Bradley et al., 2020). The 3D dose was calculated with a grid of 3 mm in the treatment planning system (TPS). The gross tumor volume was excluded from the lung volume with manually contouring by a qualified physician.
Function-wise lung region
In this study, the functional images were generated using a previously developed deep learning neural network, which can translate the pulmonary anatomy information into functional information (Ren et al., 2021a; Ren et al., 2021b). In general, a 3D attention residual neural network was utilized to extract high level features from CT images and synthesize the perfusion images. This model was trained with CT and single-photon emission computerized tomography (SPECT) perfusion images of lung disease patients. This model used a 3D encoding-decoding structure to capture the hierarchical texture features of the input CT images with two attention modules to help focus on the defect regions, which is able to achieve a medium-to-high approximation with the ground truth SPECT perfusion images.
After image synthesis, the functional image was normalized to the range of 0–1 by subtracting the minimum value and then divided by the maximum pixel value of function images. A threshold of 0.3 was used to divide the lung into high- and low- functional regions. Then the high- and low- functional lung regions were segmented on the CT image. This procedure is illustrated in Figure 2. After thresholding, the FWL was defined as the combination of the high- and low- functional lung regions. Besides, the WL region was also utilized as the basic comparison model.
Feature extraction
In the study, radiomics features (Lambin et al., 2017) and dose features of the previous regions were extracted from the CT images and 3D dose distribution.
For radiomics features, the first-order and high-order radiomics features were extracted based on the original image and 11 filter-based images. The details of these radiomics features were described in the study (Lam et al., 2021). The only difference was the settings of bin counts, in the range of [20, 50, 100, 150, 200] and a total of 5,474 radiomics features were extracted from a region of interest (ROI).
The dose features can be categorized into three types: 1) scale-invariant 3D dose statistical moments (Pham et al., 2011), 2) DVH parameters (Marks et al., 2010; Faught et al., 2017), and 3) dosiomics features (Liang et al., 2019). The scale-invariant 3D dose statistical moments described the dose spatial distribution along three directions of anterior-posterior, medial-lateral, and craniocaudal (Pham et al., 2011). Except for the constant value of the order of [0, 0, 0], a total of 63 dose statistical moments were employed in the dose features. The DVH parameters consisted of Dx and Vx, where Dx is the dose larger than x% volume, and Vx is the volume larger than the x Gy or x% of the prescription dose. A total of 59 DVH parameters were included. The dosiomics features were radiomics features based on the image of 3D dose distribution. In the study, only the original image type was adopted in extracting dosiomics features. A number of 91 dosiomics features were extracted from the original 3D dose distribution in an ROI. Eventually, a total of 213 dose features were included.
In the study, two kinds of regions (WL and FWL) were used to extract features with a total of 5,687 features and 11,374 features, respectively.
Feature selection
The feature dimension reduction is a crucial step to avoid model overfitting or underfitting. A combination of the F-test and the Pearson correlation test was utilized for the feature selection on the scikit-learn package in Python (version 1.0.1) (Pedregosa et al., 2011; Buitinck et al., 2013). Besides, the randomly under-sampling method was performed for comprehensively screening out the optimal feature group, which was described in the study (Yu et al., 2019; Lam et al., 2021).
The detail of feature selection is illustrated in Figure 3. 70% patients were randomly under-sampled from the whole patient cohorts by 100 times. At each sampling, all features with a variance of zero were filtered out to reduce the feature dimensions and the subsequent computational complexity. After that, the F-score of all features was calculated by combing the label data based on the F-test, and an F-score with a p-value smaller than 0.01 was marked as 1, otherwise as 0. Through 100 times sampling, a matric with 100×N (N: feature quantity with variance >0) was obtained. It is followed by the frequency filtering process to acquire more stable and robust features. Then, 10% of the quantity of all features or at least 40 features were reserved. Finally, the primary feature group was chosen by the Pearson correlation test with the threshold of coefficient of 0.5 as keeping the higher frequency one for two correlated features.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971849/fphar-13-971849-HTML/image_m/fphar-13-971849-g003.jpg)
FIGURE 3. The scheme of randomly under-sampled feature selection method using unsupervised ans supervised feature selection algorithms.
Model construction and evaluation
In the study, two single-omics models, radiomics model (R) and dosiomics model (D), and the combined model (RD) were developed for WL and FWL regions separately, producing six models in total, labeled as WL-R, WL-D, WL-RD, FWL-R, FWL-D, and FWL-RD.
The schematic diagram of the model development and evaluation is shown in Figure 4. All patient cohort was randomly divided into training and testing cohorts with a ratio of 3:1 across a repeat stratified splitting process of 30 times with different randomization, which simulated various patients’ data distributions to assess the model performance. At each split, training cohorts were sent to the procedure of feature selection, and the relevant primary feature group was obtained. Then, different feature combinations owning from one to all primary features were explored. The finally optimal feature group was determined by the maximum of the following overall average area under the receiver operating characteristic curve (ROC) curve (AUC) in the testing cohort. With the optimal feature combination, a classification regression algorithm of Ridge was utilized to develop a classification model using 10-fold cross-validation and hyper-parameters optimization search in the training cohort. The loss function for the Ridge classifier is
Model comparison and statistic analysis
For each omics feature, the model performance using the corresponding omics features extracted from the function-wise lung regions was compared against the WL region in the training and testing cohorts using the five evaluation metrics. For each involved lung region, the model using dual-omics features was compared against the single-omics features in both training and testing cohorts using five evaluation metrics. Besides, net clinical benefits for all models were investigated using decision curve analyses (DCA) (Vickers and Elkin, 2006; Vickers et al., 2019). The DCA is a method to evaluate the clinical valuation of models overcoming the limitations of both traditional statistical metrics, such as discrimination and calibration.
The two-sided paired student t-test was utilized to compare the above-mentioned models with a group of features. On the other hand, the two-sided paired student t-test was also performed for the continuous clinical characteristics, while the Chi-square test was applied for the categorical variables. A p-value smaller than 0.05 was considered statistically significant. Except for the previously mentioned five evaluation metrics, the 95% confidence interval (CI) by the Delong method with (DeLong et al., 1988) 2000 times for all metrics was provided to access the ability to discriminate between severe RP cases and non-RP cases. Statistical analysis was performed with Python 3.7 and Pingouin 0.5.0 (Vallat, 2018).
Results
Patients characteristics
A total of 126 NSCLC patients were retrospectively involved in the study. The main characteristics of the patients are listed in Table 1. As shown in the table, 50.8% of patients (64 cases) developed the radiation pneumonitis with a grade ≥2. Except for the gender with a p-value of 0.04, the other clinical factors had no statistically significant difference between severe RP cases and non-RP cases.
Optimal feature group
The final optimal features for six sets of WL-R, WL-D, WL-RD, FWL-R, FWL-D, and FWL-RD were listed in Supplementary Table S1. The model performance with different feature numbers was plotted in Supplementary Figure S1. A total of 39, 24, and 34 features were kept in the final optimal feature groups for FWL-R, FWL-D, and FWL-RD, respectively. The FWL-RD features consisted of 6 dosiomics and 28 radiomics features. For the region of WL, a total of 31, 4, and 29 features were utilized in the final optimal feature group for the R, D, and RD sets, respectively. The RD features consist of 6 dosiomics and 23 radiomics features. The feature number in the model of WL-R set was 35 with the maximum testing AUC. However, only 31 features were utilized in the final optimal feature group.
Model performance
Table 2 shows the average model performance for six feature sets of WL-D, WL-R, WL-RD, FWL-D, FWL-R, and FWL-RD in training and testing cohorts. Figure 5 shows model performance comparison between the WL and FWL models using each feature modality by considering 30 times data separations. For using dual-omics, the model using FWL-RD achieved significantly higher performance than the model using WL-RD in both training and testing cohorts, with an average AUC
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971849/fphar-13-971849-HTML/image_m/fphar-13-971849-t002.jpg)
TABLE 2. The average model performance in the training and testing cohorts using six feature sets of WL-D, WL-R, WL-RD, FWL-D, FWL-R, FWL-RD. The dark red color represents higher values.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971849/fphar-13-971849-HTML/image_m/fphar-13-971849-g005.jpg)
FIGURE 5. The comparison of model performance in the training and testing cohorts by using two region features of the whole lung and function-wise lung regions. The star means the
The decision curve analysis for all models is shown in Figure 6. The models using the FWL region’s feature performed a better clinical value than the models using WL region’s feature for both single and dual-omics. And, the model using FWL-RD achieved the highest overall net benefit across the majority of the range of reasonable threshold probabilities in both training and testing cohorts compared with the other feature group. ROC in the training and testing cohorts for all six models and their comparison in each feature modality and lung region were shown in Supplementary Figures S2, S3. For the best model with FWL-RD feature set, the weights of each final optimal features are displayed in Supplementary Table S3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971849/fphar-13-971849-HTML/image_m/fphar-13-971849-g006.jpg)
FIGURE 6. The comparison of clinical application values using decision curve analysis. The left and right plots showed the results of training and testing cohorts, respectively. The solid cyan line, solid brown line, solid light green line, solid pink line, solid green line, and solid red line represent the results of WL-D, FWL-D, WL-R, FWL-R, WL-RD, FWL-RD sets, respectively. The horizontal solid black line denotes that all patients didn’t suffer from RP. On the contrary, the dashed black line represents a condition that all patients occurred RP.
Discussion
In the study, we proposed an FWL sub-region generation method to benefit the prediction of acute radiation pneumonitis using pre-treatment imaging data. The predictability of each single omics and dual-omics (radiomics, dosiomics, and their combination) from the FWL were investigated and compared with the features from the WL region. As shown in Table 2 and Figure 6, the evaluation metrics and the decision curve analysis revealed that the FWL subregion generation method presented a significant prediction improvement in terms of radiomics and dual-omics features than using the WL region (p < 0.001), but not for dosiomics features.
For the models using FWL feature sets, the prediction accuracy has significant improvement as compared with the models using WL feature sets. This may suggest the features from both high and low functional lung regions have better prognostic power than the WL region. In FWL-R, FWL-D, and FWL-RD final features sets, there are 16, 7, and 15 features from the high functional regions, while they are 23, 2, and 19 for low functional regions. It should be noted that the high FWL dosiomics features played a more critical role in the FWL-D signatures. Several studies have showed the same conclusion for predicting RP when using the dose features from the high functional region (Yorke et al., 2002; Hunt et al., 2006; Wang et al., 2012; Hoover et al., 2014; Faught et al., 2017; Bucknell et al., 2018; Lee et al., 2018; O’Reilly et al., 2020). The low functional lung radiomics occupied a slightly more quantity than high FWL regions. The lower FWL’s radiomics signatures further lead to the more low FWL’s dual-omics signatures. It may imply that the heterogeneity of lung tissue, characterized by radiomics feature, presented both in high and low functional regions. The improvement for FWL-R to WL-R and FWL-RD to WL-RD may come from the smaller volume of region (high or low functional regions) than the WL region can benefit from extracting and distinguishing more heterogeneous radiomics features. Palm et al. (Palma et al., 2019) also found that the lower right lung has a significant correlation with radiation pneumonitis. In general, the lower functional region covers part of the lower right lung region. This may be one reason for the improvement by integrating the features from the low functional lung region.
For the FWL-RD signatures, except for one radiomics feature, the other radiomics signatures come from filtered CT images. The other radiomics signatures of FWL-RD are high-order features to describe gray level textural information of the lung region. For dosiomics signatures of FWL-RD, most of the signatures come from high-order features describing the dose distribution in the lung region or subregion. The selected final features are dominated by high-order omics features, which are also similar to previous studies (Hirose et al., 2020; Bourbonne et al., 2021; Jiang et al., 2021; Puttanawarut et al., 2022). None of DVH parameters (such as V5, V20, Dmean) were included in our data study, which is inconsistent with the previous studies (Palma et al., 2013; Glick et al., 2018; Onishi et al., 2018).
In our dataset, the threshold of 0.3 only was adopted in dividing the lung into high and low functional lung regions. Previous studies report that the threshold can be different, ranging from 20% to the value of the maximum functional lung image pixel (Seppenwoolde et al., 2000; Kawakami et al., 2007; Lavrenkov et al., 2007; Ohno et al., 2011; Ding et al., 2018). Following their method, we have investigated the model performance using three omics features from the FWL regions using a list threshold from 0.2 to 0.8 with a step of 0.1, as shown in Supplementary Figure S4. Besides, we statistically analyzed the difference for the testing AUC between the threshold of 0.3 and the others by using the t-test, as shown in Supplementary Table S4. As shown in the figure and table, except for dosiomics, the threshold of 0.3 achieved a statistical higher classification result in testing cohorts for the majority of feature groups of radiomics and dual-omics (except for the threshold of 0.2, 0.6 and 0.8 in RD feature groups with p = 0.107, 0.054 and p = 0.343 respectively), which is consistent with a previous study (Seppenwoolde et al., 2000). Based on the previous observations, we determine the threshold of 0.3 as an optimal threshold by considering three kinds of features. For dosiomics, the threshold of 0.2 obtained a maximum AUC value, which agreed with the study (Ding et al., 2018). In addition, an optimal threshold only using the high functional lung region’s omics features was also assessed with the threshold list, as shown in Supplementary Figure S5. And the corresponding statistical analysis was shown in Supplementary Table S5. However, non-ignificant improvement (p > 0.05) was observed by comparing the high-functional lung regions to the whole lung region.
The current study still faces several limitations. First, the functional lung images generated by DL model may have uncertainties. Even though the DL-based approach can make the image acquisition convenient and less costly for patients, the uncertainty caused by the DL model can cause a discrepancy in high and low functional lung regions. This variance may finally affect the correlation relation between some omics features and the RP, thus changing final signature features. This proposed FWL approach should be verified using the real perfusion lung image. Second, all involved patients were treated by the IMRT technique. The other radiotherapy, such as volumetric modulated arc therapy and proton radiotherapy, should be investigated for our proposed FWL region method to further explore its feasibility and capability. Third, the unbalanced between the small sample cohort and a large number of features can induce overfitting both in training and testing cohorts (Hawkins, 2004). To minimize this effect, we adopted randomly under-sampling method in the feature selection to enhance the stability of final feature signatures. However, a large prospective cohort should be carried out to access the validation of our proposed FWL method. Finally, the reproducibility and stability of omics features were not validated against disturbance. Some studies have demonstrated that the reproducibility and stability of features can be affected by the dose calculation grid size and algorithm (Placidi et al., 2020), CT image acquisition, ROI segmentation (Zwanenburg et al., 2019), and time or volume change in 4D-CT (Larue et al., 2017; Lafata et al., 2018), etc. Therefore, it is important to validate the feature robustness before clinical application.
Conclusion
In the study, we proposed an FWL approach to deeply explore the heterogeneous lung tissue and omics features and evaluated the approach in improvement of the prediction of the RP for lung cancer IMRT patients. The dual-omics features from different functional regions can improve the prediction of radiation pneumonitis for lung cancer patients under IMRT treatment. This function-wise dual-omics analysis method holds great promise to improve the prediction of radiation pneumonitis for lung cancer patients.
Data availability statement
The datasets generated for this study are available on request to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Institutional Review Board of the Affiliated Cancer Hospital of Zhengzhou University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
BL and GR contributed to the study for data analyzing and writing. HG and JC designed and organized the study. WG, CC, YW and, YY collected the image data. XZ, QD, and LM collected clinical data. JZ, XT, and ZM contributed to the study for model construction and evaluation. S-KL, HT, and HL reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was partly supported by the Provincial and Ministry Co-constructed Project of Henan Province Medical Science and Technology Research (SBGJ202103038 and SBGJ202102056), Henan Province Key R&D and Promotion Project (Science and Technology Research) (222102310015), research grants of General Research Fund (GRF 15103520), the Research Grants Council, and Health and Medical Research Fund (HMRF 07183266), the Food and Health Bureau, the Government of the Hong Kong Special Administrative Region.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.971849/full#supplementary-material
References
Adachi, T., Nakamura, M., Shintani, T., Mitsuyoshi, T., Kakino, R., Ogata, T., et al. (2021). Multi-institutional dose-segmented dosiomic analysis for predicting radiation pneumonitis after lung stereotactic body radiation therapy. Med. Phys. 48, 1781–1791. doi:10.1002/mp.14769
Baisden, J. M., Romney, D. A., Reish, A. G., Cai, J., Sheng, K., Jones, D. R., et al. (2007). Dose as a function of lung volume and planned treatment volume in helical tomotherapy intensity-modulated radiation therapy-based stereotactic body radiation therapy for small lung tumors. Int. J. Radiat. Oncol. Biol. Phys. 68, 1229–1237. doi:10.1016/J.IJROBP.2007.03.024
Barriger, R. B., Forquer, J. A., Brabham, J. G., Andolino, D. L., Shapiro, R. H., Henderson, M. A., et al. (2012). A dose-volume analysis of radiation pneumonitis in non-small cell lung cancer patients treated with stereotactic body radiation therapy. Int. J. Radiat. Oncology*Biology*Physics 82, 457–462. doi:10.1016/j.ijrobp.2010.08.056
Begosh-Mayne, D., Kumar, S. S., Toffel, S., Okunieff, P., and O’Dell, W. (2020). The dose-response characteristics of four NTCP models: Using a novel CT-based radiomic method to quantify radiation-induced lung density changes. Sci. Rep. 10, 1–9. doi:10.1038/s41598-020-67499-0
Bongers, E. M., Botticella, A., Palma, D. A., Haasbeek, C. J. A., Warner, A., Verbakel, W. F. A. R., et al. (2013). Predictive parameters of symptomatic radiation pneumonitis following stereotactic or hypofractionated radiotherapy delivered using volumetric modulated arcs. Radiotherapy Oncol. 109, 95–99. doi:10.1016/j.radonc.2013.10.011
Bourbonne, V., Da-Ano, R., Jaouen, V., Lucia, F., Dissaux, G., Bert, J., et al. (2021). Radiomics analysis of 3D dose distributions to predict toxicity of radiotherapy for lung cancer. Radiother. Oncol. 155, 144–150. doi:10.1016/j.radonc.2020.10.040
Bourbonne, V., Lucia, F., Dissaux, G., Bert, J., Visvikis, D., Pradier, O., et al. (2020). Pulmonary and esophageal toxicity in lung cancer treated by (Chemo)-radiotherapy: A radiomics-based prediction model. Int. J. Radiat. Oncology*Biology*Physics 108, S31. doi:10.1016/j.ijrobp.2020.07.2127
Bradley, J. D., Hu, C., Komaki, R. R., Masters, G. A., Blumenschein, G. R., Schild, S. E., et al. (2020). Long-Term results of NRG oncology RTOG 0617: Standard- versus high-dose chemoradiotherapy with or without cetuximab for unresectable stage III non-small-cell lung cancer. Jco 38, 706–714. doi:10.1200/JCO.19.01162
Bucknell, N. W., Hardcastle, N., Bressel, M., Hofman, M. S., Kron, T., Ball, D., et al. (2018). Functional lung imaging in radiation therapy for lung cancer: A systematic review and meta-analysis. Radiother. Oncol. 129, 196–208. doi:10.1016/j.radonc.2018.07.014
Buitinck, L., Louppe, G., Blondel, M., Pedregosa, F., Mueller, A., Grisel, O., et al. (2013). API design for machine learning software: Experiences from the scikit-learn project, 1–15. arXiv.
Cai, J., Malhotra, H. K., and Orton, C. G. (2014). A 3D-conformal technique is better than IMRT or VMAT for lung SBRT. Med. Phys. 41, 040601–040604. doi:10.1118/1.4856175
Chang, C.-K., Lin, H.-W., Ip, W.-H., Prince, T. A., Kulkarni, S. R., Levitan, D., et al. (2016). Large super-fast rotator hunting using the intermediate palomar transient factory. ApJS 227, 20. doi:10.3847/0067-0049/227/2/20
DeLong, E. R., DeLong, D. M., and Clarke-Pearson, D. L. (1988). Comparing the areas under two or more correlated receiver operating characteristic curves: A nonparametric approach. Biometrics 44, 837. doi:10.2307/2531595
Ding, X., Man, X., Sun, M., Xing, J., Fan, T., Gao, M., et al. (2018). Which is the optimal threshold for defining functional lung in single-photon emission computed tomography lung perfusion imaging of lung cancer patients? Nucl. Med. Commun. 39, 103–109. doi:10.1097/MNM.0000000000000774
Farr, K. P., Kallehauge, J. F., Møller, D. S., Khalil, A. A., Kramer, S., Bluhme, H., et al. (2015). Inclusion of functional information from perfusion SPECT improves predictive value of dose-volume parameters in lung toxicity outcome after radiotherapy for non-small cell lung cancer: A prospective study. Radiotherapy Oncol. 117, 9–16. doi:10.1016/j.radonc.2015.08.005
Faught, A. M., Yamamoto, T., Castillo, R., Castillo, E., Zhang, J., Miften, M., et al. (2017). Evaluating which dose-function metrics are most critical for functional-guided radiation therapy. Int. J. Radiat. Oncol. Biol. Phys. 99, 202–209. doi:10.1016/j.ijrobp.2017.03.051
Ganti, A. K. P., Loo, B. W., Bassetti, M., Blakely, C., Chiang, A., D'Amico, T. A., et al. (2021). Small cell lung cancer, version 2.2022, NCCN clinical practice guidelines in oncology. JNCCN J. Natl. Compr. Cancer Netw. 19, 1441–1464. doi:10.6004/JNCCN.2021.0058
Glick, D., Lyen, S., Kandel, S., Shapera, S., Le, L., Lindsay, P., et al. (2018). Impact of pretreatment interstitial lung disease on radiation pneumonitis and survival in patients treated with lung stereotactic body radiation therapy (SBRT). Clin. Lung Cancer 19, e219–e226. doi:10.1016/j.cllc.2017.06.021
Hawkins, D. M. (2004). The problem of overfitting. J. Chem. Inf. Comput. Sci. 44, 1–12. doi:10.1021/ci0342472
Hirose, T., Arimura, H., Ninomiya, K., Yoshitake, T., Fukunaga, J., and Shioyama, Y. (2020). Radiomic prediction of radiation pneumonitis on pretreatment planning computed tomography images prior to lung cancer stereotactic body radiation therapy. Sci. Rep. 10, 1. doi:10.1038/s41598-020-77552-7
Hoover, D. A., Reid, R. H., Wong, E., Stitt, L., Sabondjian, E., Rodrigues, G. B., et al. (2014). SPECT-based functional lung imaging for the prediction of radiation pneumonitis: A clinical and dosimetric correlation. J. Med. Imaging Radiat. Oncol. 58, 214–222. doi:10.1111/1754-9485.12145
Hunt, M. A., Jackson, A., Narayana, A., and Lee, N. (2006). Geometric factors influencing dosimetric sparing of the parotid glands using IMRT. Int. J. Radiat. Oncology*Biology*Physics 66, 296–304. doi:10.1016/j.ijrobp.2006.05.028
Jiang, W., Song, Y., Sun, Z., Qiu, J., and Shi, L. (2021). Dosimetric factors and radiomics features within different regions of interest in planning CT images for improving the prediction of radiation pneumonitis. Int. J. Radiat. Oncol. Biol. Phys. 110, 1161–1170. doi:10.1016/j.ijrobp.2021.01.049
Kawakami, Y., Suga, K., Zaki, M., Iwanaga, H., Yamashita, T., Hayashi, N., et al. (2007). Semi-quantitation of pulmonary perfusion heterogeneity on respiratory-gated inspiratory and expiratory perfusion SPECT in patients with pulmonary emphysema. Ann. Nucl. Med. 21, 577–584. doi:10.1007/s12149-007-0071-6
Kimura, T., Doi, Y., Nakashima, T., Imano, N., Kawabata, H., Nishibuchi, I., et al. (2015). Combined ventilation and perfusion imaging correlates with the dosimetric parameters of radiation pneumonitis in radiation therapy planning for lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 93, 778–787. doi:10.1016/j.ijrobp.2015.08.024
Kipritidis, J., Hofman, M. S., Siva, S., Callahan, J., Le Roux, P. Y., Woodruff, H. C., et al. (2015). Estimating lung ventilation directly from 4D CT Hounsfield unit values. Med. Phys. 43, 33–43. doi:10.1118/1.4937599
Kong, F. M., Ten Haken, R. K., Schipper, M. J., Sullivan, M. A., Chen, M., Lopez, C., et al. (2005). High-dose radiation improved local tumor control and overall survival in patients with inoperable/unresectable non-small-cell lung cancer: Long-term results of a radiation dose escalation study. Int. J. Radiat. Oncol. Biol. Phys. 63, 324–333. doi:10.1016/j.ijrobp.2005.02.010
Lafata, K., Cai, J., Wang, C., Hong, J., Kelsey, C. R., and Yin, F. F. (2018). Spatial-temporal variability of radiomic features and its effect on the classification of lung cancer histology. Phys. Med. Biol. 63, 225003. doi:10.1088/1361-6560/AAE56A
Lam, S.-K., Zhang, Y.-P., Zhang, J., Li, B., Sun, J.-C., Liu, Y.-T. C., et al. (2021). Multi-organ omics-based prediction for adaptive radiation therapy eligibility in nasopharyngeal carcinoma patients undergoing concurrent chemoradiotherapy. Front. Oncol. 11, 5406. doi:10.3389/FONC.2021.792024
Lambin, P., Leijenaar, R. T. H., Deist, T. M., Peerlings, J., de Jong, E. E. C., van Timmeren, J., et al. (2017). Radiomics: The bridge between medical imaging and personalized medicine. Nat. Rev. Clin. Oncol. 14, 749–762. doi:10.1038/nrclinonc.2017.141
Larue, R. T. H. M., Van De Voorde, L., van Timmeren, J. E., Leijenaar, R. T. H., Berbée, M., Sosef, M. N., et al. (2017). 4DCT imaging to assess radiomics feature stability: An investigation for thoracic cancers. Radiother. Oncol. 125, 147–153. doi:10.1016/j.radonc.2017.07.023
Lavrenkov, K., Christian, J. A., Partridge, M., Niotsikou, E., Cook, G., Parker, M., et al. (2007). A potential to reduce pulmonary toxicity: The use of perfusion SPECT with IMRT for functional lung avoidance in radiotherapy of non-small cell lung cancer. Radiother. Oncol. 83, 156–162. doi:10.1016/j.radonc.2007.04.005
Lee, H. J., Zeng, J., Vesselle, H. J., Patel, S. A., Rengan, R., and Bowen, S. R. (2018). Correlation of functional lung heterogeneity and dosimetry to radiation pneumonitis using perfusion SPECT/CT and FDG PET/CT imaging. Int. J. Radiat. Oncol. Biol. Phys. 102, 1255–1264. doi:10.1016/j.ijrobp.2018.05.051
Liang, B., Yan, H., Tian, Y., Chen, X., Yan, L., Zhang, T., et al. (2019). Dosiomics: Extracting 3D spatial features from dose distribution to predict incidence of radiation pneumonitis. Front. Oncol. 9, 269. doi:10.3389/fonc.2019.00269
Marks, L. B., Bentzen, S. M., Deasy, J. O., Kong, F. M., Bradley, J. D., Vogelius, I. S., et al. (2010). Radiation dose-volume effects in the lung. Int. J. Radiat. Oncology*Biology*Physics 76, S70–S76. doi:10.1016/j.ijrobp.2009.06.091
Ohno, Y., Koyama, H., Nogami, M., Takenaka, D., Onishi, Y., Matsumoto, K., et al. (2011). State-of-the-art radiological techniques improve the assessment of postoperative lung function in patients with non-small cell lung cancer. Eur. J. Radiology 77, 97–104. doi:10.1016/j.ejrad.2009.07.024
Onishi, H., Yamashita, H., Shioyama, Y., Matsumoto, Y., Takayama, K., Matsuo, Y., et al. (2018). Stereotactic body radiation therapy for patients with pulmonary interstitial change: High incidence of fatal radiation pneumonitis in a retrospective multi-institutional study. Cancers 10, 257. doi:10.3390/cancers10080257
O’Reilly, S., Jain, V., Huang, Q., Cheng, C., Teo, B. K. K., Yin, L., et al. (2020). Dose to highly functional ventilation zones improves prediction of radiation pneumonitis for proton and photon lung cancer radiation therapy. Int. J. Radiat. Oncology*Biology*Physics 107, 79–87. doi:10.1016/j.ijrobp.2020.01.014
Owen, D. R., Sun, Y., Boonstra, P. S., McFarlane, M., Viglianti, B. L., Balter, J. M., et al. (2021). Investigating the SPECT dose-function metrics associated with radiation-induced lung toxicity risk in patients with non-small cell lung cancer undergoing radiation therapy. Adv. Radiat. Oncol. 6, 100666. doi:10.1016/J.ADRO.2021.100666
Palma, D. A., Senan, S., Tsujino, K., Barriger, R. B., Rengan, R., Moreno, M., et al. (2013). Predicting radiation pneumonitis after chemoradiation therapy for lung cancer: An international individual patient data meta-analysis. Int. J. Radiat. Oncology*Biology*Physics 85, 444–450. doi:10.1016/j.ijrobp.2012.04.043
Palma, G., Monti, S., Xu, T., Scifoni, E., Yang, P., Hahn, S. M., et al. (2019). Spatial dose patterns associated with radiation pneumonitis in a randomized trial comparing intensity-modulated photon therapy with passive scattering proton therapy for locally advanced non-small cell lung cancer. Int. J. Radiat. Oncology*Biology*Physics 104, 1124–1132. doi:10.1016/j.ijrobp.2019.02.039
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., et al. (2011) Scikit-learn: Machine learning in Python, J. Mach. Learn. Res. 12 2825–2830.
Pham, M. T., Woodford, O. J., Perbet, F., Maki, A., Stenger, B., and Cipolla, R. (2011). A new distance for scale-invariant 3D shape recognition and registration. Proc. IEEE Int. Conf. Comput. Vis. 2011. 145–152. doi:10.1109/ICCV.2011.6126236
Pinnix, C. C., Smith, G. L., Milgrom, S., Osborne, E. M., Reddy, J. P., Akhtari, M., et al. (2015). Predictors of radiation pneumonitis in patients receiving intensity modulated radiation therapy for Hodgkin and non-hodgkin lymphoma. Int. J. Radiat. Oncol. Biol. Phys. 92, 175–182. doi:10.1016/j.ijrobp.2015.02.010
Placidi, L., Lenkowicz, J., Cusumano, D., Boldrini, L., Dinapoli, N., and Valentini, V. (2020). Stability of dosomics features extraction on grid resolution and algorithm for radiotherapy dose calculation. Phys. Medica 77, 30–35. doi:10.1016/j.ejmp.2020.07.022
Puttanawarut, C., Sirirutbunkajorn, N., Tawong, N., Jiarpinitnun, C., Khachonkham, S., Pattaranutaporn, P., et al. (2022). Radiomic and dosiomic features for the prediction of radiation pneumonitis across esophageal cancer and lung cancer. Front. Oncol. 12, 1–9. doi:10.3389/fonc.2022.768152
Ren, G., Lam, S. K., Zhang, J., Xiao, H., Cheung, A. L., Ho, W. Y., et al. (2021). Investigation of a novel deep learning-based computed tomography perfusion mapping framework for functional lung avoidance radiotherapy. Front. Oncol. 11, 1–12. doi:10.3389/fonc.2021.644703
Ren, G., Zhang, J., Li, T., Xiao, H., Cheung, L. Y., Ho, W. Y., et al. (2021). Deep learning-based computed tomography perfusion mapping (DL-CTPM) for pulmonary CT-to-perfusion translation. Int. J. Radiat. Oncol. Biol. Phys. 110, 1508–1518. doi:10.1016/j.ijrobp.2021.02.032
Seppenwoolde, Y., Muller, S. H., Theuws, J. C. M., Baas, P., Belderbos, J. S. A., Boersma, L. J., et al. (2000). Radiation dose-effect relations and local recovery in perfusion for patients with non-small-cell lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 47, 681–690. doi:10.1016/S0360-3016(00)00454-5
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J. Clin. 71, 209–249. doi:10.3322/caac.21660
Vickers, A. J., and Elkin, E. B. (2006). Decision curve analysis: A novel method for evaluating prediction models. Med. Decis. Mak. 26, 565–574. doi:10.1177/0272989X06295361
Vickers, A. J., van Calster, B., and Steyerberg, E. W. (2019). A simple, step-by-step guide to interpreting decision curve analysis. Diagn Progn. Res. 3, 1–8. doi:10.1186/s41512-019-0064-7
Wang, D., Sun, J., Zhu, J., Li, X., Zhen, Y., and Sui, S. (2012). Functional dosimetric metrics for predicting radiation-induced lung injury in non-small cell lung cancer patients treated with chemoradiotherapy. Radiat. Oncol. 7, 1. doi:10.1186/1748-717X-7-69
Wang, Z., Chen, M., Sun, J., Jiang, S., Wang, L., Wang, X., et al. (2020). Lyman-Kutcher-Burman normal tissue complication probability modeling for radiation-induced esophagitis in non-small cell lung cancer patients receiving proton radiotherapy. Radiother. Oncol. 146, 200–204. doi:10.1016/j.radonc.2020.03.003
Weller, A., Dunlop, A., Oxer, A., Gunapala, R., Murray, I., Gray, M. J., et al. (2019). Spect perfusion imaging versus CT for predicting radiation injury to normal lung in lung cancer patients. Br. J. Radiol. 92, 20190184–20190189. doi:10.1259/bjr.20190184
Xiao, L., Yang, G., Chen, J., Yang, Y., Meng, X., Wang, X., et al. (2018). Comparison of predictive powers of functional and anatomic dosimetric parameters for radiation-induced lung toxicity in locally advanced non-small cell lung cancer. Radiother. Oncol. 129, 242–248. doi:10.1016/j.radonc.2018.09.005
Yorke, E. D., Jackson, A., Rosenzweig, K. E., Merrick, S. A., Gabrys, D., Venkatraman, E. S., et al. (2002). Dose-volume factors contributing to the incidence of radiation pneumonitis in non-small-cell lung cancer patients treated with three-dimensional conformal radiation therapy. Int. J. Radiat. Oncol. Biol. Phys. 54, 329–339. doi:10.1016/S0360-3016(02)02929-2
Yu, T.-T., Lam, S.-K., Kar, F., Lee, -H., Au, K. H., Wai-Yi Yip, C., et al. (2019). Constructing novel prognostic biomarkers of advanced nasopharyngeal carcinoma from multiparametric MRI radiomics using ensemble-model based iterative feature selection. arXiv.
Zhang, X. J., Sun, J. G., Sun, J., Ming, H., Wang, X. X., Wu, L., et al. (2012). Prediction of radiation pneumonitis in lung cancer patients: A systematic review. J. Cancer Res. Clin. Oncol. 138, 2103–2116. doi:10.1007/S00432-012-1284-1
Keywords: lung functional imaging, radiation pneumonitis, radiomics, dosiomics, radiotherapy
Citation: Li B, Ren G, Guo W, Zhang J, Lam S-K, Zheng X, Teng X, Wang Y, Yang Y, Dan Q, Meng L, Ma Z, Cheng C, Tao H, Lei H, Cai J and Ge H (2022) Function-Wise Dual-Omics analysis for radiation pneumonitis prediction in lung cancer patients. Front. Pharmacol. 13:971849. doi: 10.3389/fphar.2022.971849
Received: 17 June 2022; Accepted: 04 August 2022;
Published: 19 September 2022.
Edited by:
Khairunnisa Hasikin, University of Malaya, MalaysiaReviewed by:
Jun Dang, Chinese Academy of Medical Sciences and Peking Union Medical College, ChinaXiangkun Xu, University of Texas Southwestern Medical Center, United States
Qingyi Liu, Shandong University of Science and Technology, China
Copyright © 2022 Li, Ren, Guo, Zhang, Lam, Zheng, Teng, Wang, Yang, Dan, Meng, Ma, Cheng, Tao, Lei, Cai and Ge. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hong Ge, gehong666@126.com; Jing Cai, jing.cai@polyu.edu.hk