Regional differences in lung radiosensitivity after radiotherapy for non-small-cell lung cancer

Quantifying the Effect of 4-Dimensional Computed Tomography-Based Deformable Dose Accumulation on Representing Radiation Damage for Patients with Locally Advanced Non-Small Cell Lung Cancer Treated with Standard-Fractionated Intensity-Modulated Radiation Therapy

PURPOSE

The aim of this study was to investigate the dosimetric and clinical effects of 4-dimensional computed tomography (4DCT)-based longitudinal dose accumulation in patients with locally advanced non-small cell lung cancer treated with standard-fractionated intensity-modulated radiation therapy (IMRT).

METHODS AND MATERIALS

Sixty-seven patients were retrospectively selected from a randomized clinical trial. Their original IMRT plan, planning and verification 4DCTs, and ∼4-month posttreatment follow-up CTs were imported into a commercial treatment planning system. Two deformable image registration algorithms were implemented for dose accumulation, and their accuracies were assessed. The planned and accumulated doses computed using average-intensity images or phase images were compared. At the organ level, mean lung dose and normal-tissue complication probability (NTCP) for grade ≥2 radiation pneumonitis were compared. At the region level, mean dose in lung subsections and the volumetric overlap between isodose intervals were compared. At the voxel level, the accuracy in estimating the delivered dose was compared by evaluating the fit of a dose versus radiographic image density change (IDC) model. The dose-IDC model fit was also compared for subcohorts based on the magnitude of NTCP difference (|ΔNTCP|) between planned and accumulated doses.

RESULTS

Deformable image registration accuracy was quantified, and the uncertainty was considered for the voxel-level analysis. Compared with planned doses, accumulated doses on average resulted in <1-Gy lung dose increase and <2% NTCP increase (up to 8.2 Gy and 18.8% for a patient, respectively). Volumetric overlap of isodose intervals between the planned and accumulated dose distributions ranged from 0.01 to 0.93. Voxel-level dose-IDC models demonstrated a fit improvement from planned dose to accumulated dose (pseudo-R2 increased 0.0023) and a further improvement for patients with ≥2% |ΔNTCP| versus for patients with <2% |ΔNTCP|.

CONCLUSIONS

With a relatively large cohort, robust image registrations, multilevel metric comparisons, and radiographic image-based evidence, we demonstrated that dose accumulation more accurately represents the delivered dose and can be especially beneficial for patients with greater longitudinal response.

Voxel-wise analysis: A powerful tool to predict radio-induced toxicity and potentially perform personalised planning in radiotherapy

Dose - volume histograms have been historically used to study the relationship between the planned radiation dose and healthy tissue damage. However, this approach considers neither spatial information nor heterogenous radiosensitivity within organs at risk, depending on the tissue. Recently, voxel-wise analyses have emerged in the literature as powerful tools to fully exploit three-dimensional information from the planned dose distribution. They allow to identify anatomical subregions of one or several organs in which the irradiation dose is associated with a given toxicity. These methods rely on an accurate anatomical alignment, usually obtained by means of a non-rigid registration. Once the different anatomies are spatially normalised, correlations between the three-dimensional dose and a given toxicity can be explored voxel-wise. Parametric or non-parametric statistical tests can be performed on every voxel to identify the voxels in which the dose is significantly different between patients presenting or not toxicity. Several anatomical subregions associated with genitourinary, gastrointestinal, cardiac, pulmonary or haematological toxicity have already been identified in the literature for prostate, head and neck or thorax irradiation. Voxel-wise analysis appears therefore first particularly interesting to increase toxicity prediction capability by identifying specific subregions in the organs at risk whose irradiation is highly predictive of specific toxicity. The second interest is potentially to decrease the radio-induced toxicity by limiting the dose in the predictive subregions, while not decreasing the dose in the target volume. Limitations of the approach have been pointed out.

Ipsilateral lung dose as a correlative measure for radiation pneumonitis in patients treated with definitive concurrent radiochemotherapy

Objective

Mean lung dose (MLD) and percent of total lung (TL) volume that receive a dose greater than 20 Gy (V20) have been the most validated parameters in the prediction of radiation pneumonitis (RP). However, these parameters present mean values of TL parenchyma and predict the right and the left lung as a unique functional organ unit, not take into account the difference in function and dose density between the lungs. Furthermore, there have been very limited data evaluating ipsilateral lung dosimetric constraints in addition to TL parameters to predict RP in non-small cell lung cancer (NSCLC) patients treated with radiochemotherapy (RCT).

Methods

Between 2010 and 2017, clinical-radiological findings of NSCLC patients treated with RCT were evaluated in terms of RP, retrospectively. MLD, V20, and V30 values of ipsilateral lung were assessed from dose-volume histogram and registered. The primary endpoint was to assess the relation between ipsilateral lung dose constraints and RP risk.

Results

There were 75 patients. There was ≥Grade 2 RP in 33 cases (%44). In univariate analysis, ipsilateral MLD, ipsilateral V20, ipsilateral V30, and TL V30 were found to be significant. Ipsilateral MLD and PTV were found to be the independent risk factors for RP. Cutoff values for RP risk were determined as 18Gy, 35%, and 28% for ipsilateral MLD, ipsilateral V20, and ipsilateral V30, respectively. Predictive values for ipsilateral MLD and ipsilateral V20 were higher than TL.

Conclusions

In NSCLC patients treated with RCT, MLD, V20, and V30 values of ipsilateral lung parameters might increase the predictability of RP risk in addition to TL parameters.

Advances in Knowledge

Cutoff values for RP risk were determined as 18Gy, 35%, and 28% for ipsilateral MLD, ipsilateral V20, and ipsilateral V30, respectively. Predictive values for ipsilateral MLD and ipsilateral V20 were higher than TL.

Does shortening the duration of radiotherapy treatment in breast cancer increase the risk of radiation pneumonia: A retrospective study

Randomized studies evaluating hypofractionation and conventional fractionation radiotherapy treatments (RT) in patients with breast cancer have shown that hypofractionation achieves similar results to conventional fractionation in terms of survival and local control rates. It has also been shown that their long-term toxicities are similar. This study aimed to evaluate the effects of hypofractionated radiotherapy (H-RT) and conventional radiotherapy (C-RT) on lung toxicity and identify factors affecting this toxicity in patients with breast cancer. The study included 118 patients who underwent adjuvant RT following breast-conserving surgery (BCS). Out of these, 63 patients were assigned to receive C-RT, while the remaining 55 were assigned to receive H-RT. To clarify, we treated 63 patients with C-RT and 55 patients with H-RT. 60 patients were treated using 3-dimensional conformal radiotherapy (3DCRT) and 58 patients were treated using intensity modulated radiotherapy (IMRT). The patients were evaluated weekly for toxicity during radiotherapy (RT) treatment and were called every 3 months for routine controls after the treatment. The first control was performed 1 month after the treatment. Statistical analysis was performed using the SPSS20 program, and a P value of <.005 was considered statistically significant. The study found that the median age of the participants was 54.9 years and tomographic findings were observed in 70 patients. Radiological findings were detected at a median of 5 months after RT. The mean lung dose (MLD) on the treated breast side (referred to as ipsilateral lung or OAR) was 10.4 Gy for the entire group. Among patients who received 18 MV energy in RT, those with an area volume (V20) of the lung on the treated breast side >18.5%, those with a mean dose of the treated breast side lung (ipsilateral lung) >10.5 Gy, and those who received concurrent hormone therapy had significantly more tomographic findings. However, patients treated with YART had fewer tomographic findings. No symptomatic patients were observed during the follow-up period. Our findings show that the risk of lung toxicity is similar with H-RT and C-RT, and H-RT can be considered an effective and safe treatment option for breast cancer. The key factors affecting the development of lung toxicity were found to be the type of RT energy used, RT to the side breast, volume receiving 20 Gy in the side lung, side lung mean dose, and simultaneous hormonal therapy.

Is individual perfusion dose-response different than ventilation dose-response for lung cancer patients treated with radiotherapy?

OBJECTIVE

Current ventilation and perfusion dose-response studies focus on single-modalities (ventilation or perfusion) and perform pulmonary-toxicity assessment related to radiotherapy on a population-based basis. This study aims at quantitative and clinical evaluation of intrapatient differences between ventilation and perfusion dose-responses among lung cancer patients treated with radiotherapy.

METHODS

20 patients enrolled on a prospective functional avoidance protocol underwent single photon emission computed tomography-CT ventilation and perfusion scans pre- and post-radiotherapy. Relative changes in pre- to post-treatment ventilation and perfusion in lung regions receiving ≥20 Gy were calculated. In addition, the slopes of the linear fit to the relative ventilation and perfusion changes in regions receiving 0-60 Gy were calculated. A radiologist read and assigned a functional defect score to pre- and post-treatment ventilation/perfusion scans.

RESULTS

25% of patients had a difference >35% between ventilation and perfusion pre- to post-treatment changes and 20-30% of patients had opposite directions for ventilation and perfusion pre- to post-treatment changes. Using a semi-quantitative scale, radiologist assessment showed that 20% of patients had different pre- to post-treatment ventilation changes when compared to pre- to post-treatment perfusion changes.

CONCLUSION

Our data showed that ventilation dose-response can differ from perfusion dose-response for 20-30% of patients. Therefore, when performing thoracic dose-response in cancer patients, it is insufficient to look at ventilation or perfusion alone; but rather both modes of functional imaging may be needed when predicting for clinical outcomes.

ADVANCES IN KNOWLEDGE

The significance of this study can be highlighted by the differences between the intrapatient dose-response assessments of this analysis compared to existing population-based dose-response analyses. Elucidating intrapatient ventilation and perfusion dose-response differences may be valuable in predicting pulmonary toxicity in lung cancer patients post-radiotherapy.

Technical note: rtdsm-An open-source software for radiotherapy dose-surface map generation and analysis

BACKGROUND

Dose-outcome studies in radiation oncology have historically excluded spatial information due to dose-volume histograms being the most dominant source of dosimetric information. In recent years, dose-surface maps (DSMs) have become increasingly popular for characterization of spatial dose distributions and identification of radiosensitive subregions for hollow organs. However, methodological variations and lack of open-source, publicly offered code-sharing between research groups have limited reproducibility and wider adoption.

PURPOSE

This paper presents rtdsm, an open-source software for DSM calculation with the intent to improve the reproducibility of and the access to DSM-based research in medical physics and radiation oncology.

METHODS

A literature review was conducted to identify essential functionalities and prevailing calculation approaches to guide development. The described software has been designed to calculate DSMs from DICOM data with a high degree of user customizability and to facilitate DSM feature analysis. Core functionalities include DSM calculation, equivalent dose conversions, common DSM feature extraction, and simple DSM accumulation.

RESULTS

A number of use cases were used to qualitatively and quantitatively demonstrate the use and usefulness of rtdsm. Specifically, two DSM slicing methods, planar and noncoplanar, were implemented and tested, and the effects of method choice on output DSMs were demonstrated. An example comparison of DSMs from two different treatments was used to highlight the use cases of various built-in analysis functions for equivalent dose conversion and DSM feature extraction.

CONCLUSIONS

We developed and implemented rtdsm as a standalone software that provides all essential functionalities required to perform a DSM-based study. It has been made freely accessible under an open-source license on Github to encourage collaboration and community use.

Lung Subregion Partitioning by Incremental Dose Intervals Improves Omics-Based Prediction for Acute Radiation Pneumonitis in Non-Small-Cell Lung Cancer Patients

Purpose: To evaluate the effectiveness of features obtained from our proposed incremental-dose-interval-based lung subregion segmentation (IDLSS) for predicting grade ≥ 2 acute radiation pneumonitis (ARP) in lung cancer patients upon intensity-modulated radiotherapy (IMRT). (1) Materials and Methods: A total of 126 non-small-cell lung cancer patients treated with IMRT were retrospectively analyzed. Five lung subregions (SRs) were generated by the intersection of the whole lung (WL) and five sub-regions receiving incremental dose intervals. A total of 4610 radiomics features (RF) from pre-treatment planning computed tomographic (CT) and 213 dosiomics features (DF) were extracted. Six feature groups, including WL-RF, WL-DF, SR-RF, SR-DF, and the combined feature sets of WL-RDF and SR-RDF, were generated. Features were selected by using a variance threshold, followed by a Student t-test. Pearson’s correlation test was applied to remove redundant features. Subsequently, Ridge regression was adopted to develop six models for ARP using the six feature groups. Thirty iterations of resampling were implemented to assess overall model performance by using the area under the Receiver-Operating-Characteristic curve (AUC), accuracy, precision, recall, and F1-score. (2) Results: The SR-RDF model achieved the best classification performance and provided significantly better predictability than the WL-RDF model in training cohort (Average AUC: 0.98 ± 0.01 vs. 0.90 ± 0.02, p < 0.001) and testing cohort (Average AUC: 0.88 ± 0.05 vs. 0.80 ± 0.04, p < 0.001). Similarly, predictability of the SR-DF model was significantly stronger than that of the WL-DF model in training cohort (Average AUC: 0.88 ± 0.03 vs. 0.70 ± 0.030, p < 0.001) and in testing cohort (Average AUC: 0.74 ± 0.08 vs. 0.65 ± 0.06, p < 0.001). By contrast, the SR-RF model significantly outperformed the WL-RF model only in the training set (Average AUC: 0.93 ± 0.02 vs. 0.85 ± 0.03, p < 0.001), but not in the testing set (Average AUC: 0.79 ± 0.05 vs. 0.77 ± 0.07, p = 0.13). (3) Conclusions: Our results demonstrated that the IDLSS method improved model performance for classifying ARP with grade ≥ 2 when using dosiomics or combined radiomics-dosiomics features.

Radiation and immune checkpoint inhibitor-mediated pneumonitis risk stratification in patients with locally advanced non-small cell lung cancer: role of functional lung radiomics?

BACKGROUND

Patients undergoing chemoradiation and immune checkpoint inhibitor (ICI) therapy for locally advanced non-small cell lung cancer (NSCLC) experience pulmonary toxicity at higher rates than historical reports. Identifying biomarkers beyond conventional clinical factors and radiation dosimetry is especially relevant in the modern cancer immunotherapy era. We investigated the role of novel functional lung radiomics, relative to functional lung dosimetry and clinical characteristics, for pneumonitis risk stratification in locally advanced NSCLC.

METHODS

Patients with locally advanced NSCLC were prospectively enrolled on the FLARE-RT trial (NCT02773238). All received concurrent chemoradiation using functional lung avoidance planning, while approximately half received consolidation durvalumab ICI. Within tumour-subtracted lung regions, 110 radiomics features (size, shape, intensity, texture) were extracted on pre-treatment [^99mTc]MAA SPECT/CT perfusion images using fixed-bin-width discretization. The performance of functional lung radiomics for pneumonitis (CTCAE v4 grade 2 or higher) risk stratification was benchmarked against previously reported lung dosimetric parameters and clinical risk factors. Multivariate least absolute shrinkage and selection operator Cox models of time-varying pneumonitis risk were constructed, and prediction performance was evaluated using optimism-adjusted concordance index (c-index) with 95% confidence interval reporting throughout.

RESULTS

Thirty-nine patients were included in the study and pneumonitis occurred in 16/39 (41%) patients. Among clinical characteristics and anatomic/functional lung dosimetry variables, only the presence of baseline chronic obstructive pulmonary disease (COPD) was significantly associated with the development of pneumonitis (HR 4.59 [1.69-12.49]) and served as the primary prediction benchmark model (c-index 0.69 [0.59-0.80]). Discrimination of time-varying pneumonitis risk was numerically higher when combining COPD with perfused lung radiomics size (c-index 0.77 [0.65-0.88]) or shape feature classes (c-index 0.79 [0.66-0.91]) but did not reach statistical significance compared to benchmark models (p > 0.26). COPD was associated with perfused lung radiomics size features, including patients with larger lung volumes (AUC 0.75 [0.59-0.91]). Perfused lung radiomic texture features were correlated with lung volume (adj R² = 0.84-1.00), representing surrogates rather than independent predictors of pneumonitis risk.

CONCLUSIONS

In patients undergoing chemoradiation with functional lung avoidance therapy and optional consolidative immune checkpoint inhibitor therapy for locally advanced NSCLC, the strongest predictor of pneumonitis was the presence of baseline chronic obstructive pulmonary disease. Results from this novel functional lung radiomics exploratory study can inform future validation studies to refine pneumonitis risk models following combinations of radiation and immunotherapy. Our results support functional lung radiomics as surrogates of COPD for non-invasive monitoring during and after treatment. Further study of clinical, dosimetric, and radiomic feature combinations for radiation and immune-mediated pneumonitis risk stratification in a larger patient population is warranted.

Radiomics-guided radiation therapy: opportunities and challenges

Radiomics is an advanced image-processing framework, which extracts image features and considers them as biomarkers towards personalized medicine. Applications include disease detection, diagnosis, prognosis, and therapy response assessment/prediction. As radiation therapy aims for further individualized treatments, radiomics could play a critical role in various steps before, during and after treatment. Elucidation of the concept of radiomics-guided radiation therapy (RGRT) is the aim of this review, attempting to highlight opportunities and challenges underlying the use of radiomics to guide clinicians and physicists towards more effective radiation treatments. This work identifies the value of RGRT in various steps of radiotherapy from patient selection to follow-up, and subsequently provides recommendations to improve future radiotherapy using quantitative imaging features.

Symptomatic Radiation Pneumonitis in NSCLC Patients Receiving EGFR-TKIs and Concurrent Once-daily Thoracic Radiotherapy: Predicting the Value of Clinical and Dose-volume Histogram Parameters

BACKGROUND

The incidence of symptomatic radiation pneumonitis (RP) and its relationship with dose-volume histogram (DVH) parameters in non-small cell lung cancer (NSCLC) patients receiving epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs) and concurrent once-daily thoracic radiotherapy (TRT) remain unclear. We aim to analyze the values of clinical factors and dose-volume histogram (DVH) parameters to predict the risk for symptomatic RP in these patients.

METHODS

Between 2011 and 2019, we retrospectively analyzed and identified 85 patients who had received EGFR-TKIs and once-daily TRT simultaneously (EGFR-TKIs group) and 129 patients who had received concurrent chemoradiotherapy (CCRT group). The symptomatic RP was recorded according to the Common Terminology Criteria for Adverse Event (CTCAE) criteria (grade 2 or above). Statistical analyses were performed using SPSS 26.0.

RESULTS

In total, the incidences of symptomatic (grade≥2) and severe RP (grade≥3) were 43.5% (37/85) and 16.5% (14/85) in EGFR-TKIs group vs 27.1% (35/129) and 10.1% (13/129) in CCRT group respectively. After 1:1 ratio between EGFR-TKIs group and CCRT group was matched by propensity score matching, chi-square test suggested that the incidence of symptomatic RP in the MATCHED EGFR-TKIs group was higher than that in the matched CCRT group (χ2=4.469, P=0.035). In EGFR-TKIs group, univariate and multivariate analyses indicated that the percentage of ipsilateral lung volume receiving ≥30 Gy (ilV30) [odds ratio (OR): 1.163, 95%CI: 1.036-1.306, P=0.011] and the percentage of total lung volume receiving ≥20 Gy (tlV20) (OR: 1.171, 95%CI: 1.031-1.330, P=0.015), with chronic obstructive pulmonary disease (COPD) or not (OR: 0.158, 95%CI: 0.041-0.600, P=0.007), were independent predictors of symptomatic RP. Compared to patients with lower ilV30/tlV20 values (ilV30 and tlV20<cut-off point values) and without COPD, patients with higher ilV30/tlV20 values (ilV30 and tlV20>cut-off point values) and COPD had a significantly higher risk for developing symptomatic RP, with a hazard ratio (HR) of 1.350 (95%CI: 1.190-1.531, P<0.001).

CONCLUSIONS

Patients receiving both EGFR-TKIs and once-daily TRT were more likely to develop symptomatic RP than patients receiving concurrent chemoradiotherapy. The ilV30, tlV20, and comorbidity of COPD may predict the risk of symptomatic RP among NSCLC patients receiving EGFR-TKIs and conventionally fractionated TRT concurrently.

Combining Serial and Parallel Functionality in Functional Lung Avoidance Radiation Therapy

PURPOSE

Functional lung avoidance (FLA) radiation therapy (RT) aims to minimize post-RT pulmonary toxicity by preferentially avoiding dose to high-functioning lung (HFL) regions. A common limitation is that FLA approaches do not consider the conducting architecture for gas exchange. We previously proposed the functionally weighted airway sparing (FWAS) method to spare airways connected to HFL regions, showing that it is possible to substantially reduce risk of radiation-induced airway injury. Here, we compare the performance of FLA and FWAS and propose a novel method combining both approaches.

METHODS

We used breath-hold computed tomography (BHCT) and simulation 4-dimensional computed tomography (4DCT) from 12 lung stereotactic ablative radiation therapy patients. Four planning strategies were examined: (1) Conventional: no sparing other than clinical dose-volume constraints; (2) FLA: using a 4DCT-based ventilation map to delineate the HFL, plans were optimized to reduce mean dose and V_13.50 in HFL; (3) FWAS: we autosegemented 11 to 13 generations of individual airways from each patient's BHCT and assigned priorities based on the relative contribution of each airway to total ventilation. We used these priorities in the optimization along with airway dose constraints, estimated as a function of airway diameter and 5% probability of collapse; and (4) FLA + FWAS: we combined information from the 2 strategies. We prioritized clinical dose constraints for organs at risk and planning target volume in all plans. We performed the evaluation in terms of ventilation preservation accounting for radiation-induced damage to both lung parenchyma and airways.

RESULTS

We observed average ventilation preservation for FLA, FWAS, and FLA + FWAS as 3%, 8.5%, and 14.5% higher, respectively, than for Conventional plans for patients with ventilation preservation in Conventional plans <90%. Generalized estimated equations showed that all improvements were statistically significant (P ≤ .036). We observed no clinically relevant improvements in outcomes of the sparing techniques in patients with ventilation preservation in Conventional plans ≥90%.

CONCLUSIONS

These initial results suggest that it is crucial to consider the parallel and the serial nature of the lung to improve post-radiation therapy lung function and, consequently, quality of life for patients.

Multi-institutional phase II study on the safety and efficacy of dynamic tumor tracking-stereotactic body radiotherapy for lung tumors

BACKGROUND AND PURPOSE

This study aimed to evaluate the safety and efficacy of dynamic tumor tracking-stereotactic body radiotherapy (DTT-SBRT) for lung tumors.

MATERIALS AND METHODS

Patients with cStage I primary lung cancer or metastatic lung cancer with an expected range of respiratory motion of ≥10 mm were eligible for the study. The prescribed dose was 50 Gy in four fractions. A gimbal-mounted linac was used for DTT-SBRT delivery. The primary endpoint was local control at 2 years.

RESULTS

Forty-eight patients from four institutions were enrolled in this study. Forty-two patients had primary non-small-cell lung cancer, and six had metastatic lung tumors. DTT-SBRT was delivered for 47 lesions in 47 patients with a median treatment time of 28 min per fraction. The median respiratory motion during the treatment was 13.7 mm (range: 4.5-28.1 mm). The motion-encompassing method was applied for the one remaining patient due to the poor correlation between the abdominal wall and tumor movement. The median follow-up period was 32.3 months, and the local control at 2 years was 95.2% (lower limit of the one-sided 85% confidence interval [CI]: 90.3%). The overall survival and progression-free survival at 2 years were 79.2% (95% CI: 64.7%-88.2%) and 75.0% (95% CI: 60.2%-85.0%), respectively. Grade 3 toxicity was observed in one patient (2.1%) with radiation pneumonitis. Grade 4 or 5 toxicity was not observed.

CONCLUSION

DTT-SBRT achieved excellent local control with low incidences of severe toxicities in lung tumors with respiratory motion.

Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis

Simple Summary

The pathophysiology of radiation pneumonitis (RP) after thoracic cancer radiation treatments is still not completely understood although the identification of underlying RP mechanisms may improve the therapeutic window of thoracic cancer patients. The aim of our retrospective study was to explore the dose–response patterns associated with RP by a multi-center voxel-based analysis. In a heterogeneously treated population of 382 thoracic cancer patients, we confirmed the previously described heart–lung interaction in the development of RP. The empowerment of VBA with a novel description of dose map spatial properties based on probabilistic independent component analysis (PICA) and connectograms provided valuable additional and independent information on the radiobiology of RP.

Abstract

This study investigates the dose–response patterns associated with radiation pneumonitis (RP) in patients treated for thoracic malignancies with different radiation modalities. To this end, voxel-based analysis (VBA) empowered by a novel strategy for the characterization of spatial properties of dose maps was applied. Data from 382 lung cancer and mediastinal lymphoma patients from three institutions treated with different radiation therapy (RT) techniques were analyzed. Each planning CT and biologically effective dose map (α/β = 3 Gy) was spatially normalized on a common anatomical reference. The VBA of local dose differences between patients with and without RP was performed and the clusters of voxels with dose differences that significantly correlated with RP at a p-level of 0.05 were generated accordingly. The robustness of VBA inference was evaluated by a novel characterization for spatial properties of dose maps based on probabilistic independent component analysis (PICA) and connectograms. This lays robust foundations to the obtained findings that the lower parts of the lungs and the heart play a prominent role in the development of RP. Connectograms showed that the dataset can support a radiobiological differentiation between the main heart and lung substructures.

Integrating Multiomics Information in Deep Learning Architectures for Joint Actuarial Outcome Prediction in Non-Small Cell Lung Cancer Patients After Radiation Therapy

PURPOSE

Novel actuarial deep learning neural network (ADNN) architectures are proposed for joint prediction of radiation therapy outcomes-radiation pneumonitis (RP) and local control (LC)-in stage III non-small cell lung cancer (NSCLC) patients. Unlike normal tissue complication probability/tumor control probability models that use dosimetric information solely, our proposed models consider complex interactions among multiomics information including positron emission tomography (PET) radiomics, cytokines, and miRNAs. Additional time-to-event information is also used in the actuarial prediction.

METHODS AND MATERIALS

Three architectures were investigated: ADNN-DVH considered dosimetric information only; ADNN-com integrated multiomics information; and ADNN-com-joint combined RP2 (RP grade ≥2) and LC prediction. In these architectures, differential dose-volume histograms (DVHs) were fed into 1D convolutional neural networks (CNN) for extracting reduced representations. Variational encoders were used to learn representations of imaging and biological data. Reduced representations were fed into Surv-Nets to predict time-to-event probabilities for RP2 and LC independently and jointly by incorporating time information into designated loss functions.

RESULTS

Models were evaluated on 117 retrospective patients and were independently tested on 25 newly accrued patients prospectively. A multi-institutional RTOG0617 data set of 327 patients was used for external validation. ADNN-DVH yielded cross-validated c-indexes (95% confidence intervals) of 0.660 (0.630-0.690) for RP2 prediction and 0.727 (0.700-0.753) for LC prediction, outperforming a generalized Lyman model for RP2 (0.613 [0.583-0.643]) and a generalized log-logistic model for LC (0.569 [0.545-0.594]). The independent internal test and external validation yielded similar results. ADNN-com achieved an even better performance than ADNN-DVH on both cross-validation and independent internal test. Furthermore, ADNN-com-joint, which yielded performance similar to ADNN-com, realized joint prediction with c-indexes of 0.705 (0.676-0.734) for RP2 and 0.740 (0.714-0.765) for LC and achieved an area under a free-response receiving operator characteristic curve (AU-FROC) of 0.729 (0.697-0.773) for the joint prediction of RP2 and LC.

CONCLUSION

Novel deep learning architectures that integrate multiomics information outperformed traditional normal tissue complication probability/tumor control probability models in actuarial prediction of RP2 and LC.

Spatial descriptions of radiotherapy dose: normal tissue complication models and statistical associations

For decades, dose-volume information for segmented anatomy has provided the essential data for correlating radiotherapy dosimetry with treatment-induced complications. Dose-volume information has formed the basis for modelling those associations via normal tissue complication probability (NTCP) models and for driving treatment planning. Limitations to this approach have been identified. Many studies have emerged demonstrating that the incorporation of information describing the spatial nature of the dose distribution, and potentially its correlation with anatomy, can provide more robust associations with toxicity and seed more general NTCP models. Such approaches are culminating in the application of computationally intensive processes such as machine learning and the application of neural networks. The opportunities these approaches have for individualising treatment, predicting toxicity and expanding the solution space for radiation therapy are substantial and have clearly widespread and disruptive potential. Impediments to reaching that potential include issues associated with data collection, model generalisation and validation. This review examines the role of spatial models of complication and summarises relevant published studies. Sources of data for these studies, appropriate statistical methodology frameworks for processing spatial dose information and extracting relevant features are described. Spatial complication modelling is consolidated as a pathway to guiding future developments towards effective, complication-free radiotherapy treatment.

Probing thoracic dose patterns associated to pericardial effusion and mortality in patients treated with photons and protons for locally advanced non-small-cell lung cancer

PURPOSE

To investigate thoracic dose-response patterns for pericardial effusion (PCE) and mortality in patients treated for locally advanced Non-Small-Cell Lung Cancer (NSCLC) by Intensity Modulated RT (IMRT) or Passive-Scattering Proton Therapy (PSPT).

METHODS

Among 178 patients, 43.5% developed grade ≥ 2 PCE. Clinical and dosimetric factors associated with PCE or overall survival (OS) were identified via multi-variable Cox proportional hazards modeling. The Voxel-Based Analyses (VBAs) of local dose differences between patients with and without PCE and mortality was performed. The robustness of VBA results was assessed by a novel characterization of spatial properties of dose distributions based on probabilistic independent component analysis (PICA) and connectograms.

RESULTS

Several non-dosimetric variables were selected by the multivariable analysis for the considered outcomes, while the time-dependent PCE onset was uncorrelated with the OS (p = 0.34) at a multi-variable Cox analysis. Despite the significant PSPT dosimetric advantage, the RT technique did not affect the occurrence of PCE or OS. VBAs highlighted largely overlapping clusters significantly associated with PCE endpoints in heart and lungs. No significant dosimetric patterns related to mortality endpoints were found. PICA identified 43 components homogeneously scattered within thorax, while connectograms showed modest correlations between doses in main cardio-pulmonary substructures.

CONCLUSIONS

Spatially resolved analysis highlighted dose patterns related to radiation-induced cardiac toxiciy and the observed organ-based dose-response mismatch in PSPT and IMRT. Indeed, the thoracic regions spared by PSPT poorly overlapped with the areas involved in PCE development, as highlited by VBA. PICA and connectograms proved valuable tools for assessing the robusteness of obtained VBA inferences.

Computed Tomography-Based Delta-Radiomics Analysis for Discriminating Radiation Pneumonitis in Patients With Esophageal Cancer After Radiation Therapy

PURPOSE

Our purpose was to construct a computed tomography (CT)-based delta-radiomics nomogram and corresponding risk classification system for individualized and accurate estimation of severe acute radiation pneumonitis (SARP) in patients with esophageal cancer (EC) after radiation therapy.

METHODS AND MATERIALS

Four hundred patients with EC were enrolled from 2 independent institutions and were divided into the training (n = 200) and validation (n = 200) cohorts. Eight hundred fifty radiomics features of lung were extracted from treatment planning images, including the positioning CT before radiation therapy (CT₁) and the resetting CT after receiving 40 to 45 Gy (CT₂). The longitudinal net changes in radiomics features from CT₁ to CT₂ were calculated and defined as delta-radiomics features. Least absolute shrinkage and selection operator algorithm was performed to features selection and delta-radiomics signature building. Integrating the signature with multidimensional clinicopathologic, dosimetric, and hematological predictors of SARP, a novel CT-based delta-radiomics nomogram was established according to multivariate analysis. The clinical application values of nomogram were both evaluated in the training and validation cohorts by concordance index, calibration curves, and decision curve analysis. Recursive partitioning analysis was used to generate a risk classification system.

RESULTS

The delta-radiomics signature consisting of 24 features was significantly associated with SARP status (P < .001). Incorporating it with other high-risk factors, Subjective Global Assessment score, pulmonary fibrosis score, mean lung dose, and systemic immune inflammation index, the developed delta-radiomics nomogram showed increased improvement in SARP discrimination accuracy with concordance index of 0.975 and 0.921 in the training and validation cohorts, respectively. Calibration curves and decision curve analysis confirmed the satisfactory clinical feasibility and utility of nomogram. The risk classification system displayed excellent performance on identifying SARP occurrence (P < .001).

CONCLUSIONS

The delta-radiomics nomogram and risk classification system as low-cost and noninvasive means exhibited superior predictive accuracy and provided individualized probability of SARP stratification for patients with EC.

Investigating the SPECT Dose-Function Metrics Associated With Radiation-Induced Lung Toxicity Risk in Patients With Non-small Cell Lung Cancer Undergoing Radiation Therapy

Purpose

Dose to normal lung has commonly been linked with radiation-induced lung toxicity (RILT) risk, but incorporating functional lung metrics in treatment planning may help further optimize dose delivery and reduce RILT incidence. The purpose of this study was to investigate the impact of the dose delivered to functional lung regions by analyzing perfusion (Q), ventilation (V), and combined V/Q single-photon-emission computed tomography (SPECT) dose-function metrics with regard to RILT risk in patients with non-small cell lung cancer (NSCLC) patients who received radiation therapy (RT).

Methods and Materials

SPECT images acquired from 88 patients with locally advanced NSCLC before undergoing conventionally fractionated RT were retrospectively analyzed. Dose was converted to the nominal dose equivalent per 2 Gy fraction, and SPECT intensities were normalized. Regional lung segments were defined, and the average dose delivered to each lung region was quantified. Three functional categorizations were defined to represent low-, normal-, and high-functioning lungs. The percent of functional lung category receiving ≥20 Gy and mean functional intensity receiving ≥20 Gy (iV₂₀) were calculated. RILT was defined as grade 2+ radiation pneumonitis and/or clinical radiation fibrosis. A logistic regression was used to evaluate the association between dose-function metrics and risk of RILT.

Results

By analyzing V/Q normalized intensities and functional distributions across the population, a wide range in functional capability (especially in the ipsilateral lung) was observed in patients with NSCLC before RT. Through multivariable regression models, global lung average dose to the lower lung was found to be significantly associated with RILT, and Q and V iV₂₀ were correlated with RILT when using ipsilateral lung metrics. Through a receiver operating characteristic analysis, combined V/Q low-function receiving ≥20 Gy (low-functioning V/Q₂₀) in the ipsilateral lung was found to be the best predictor (area under the curce: 0.79) of RILT risk.

Conclusions

Irradiation of the inferior lung appears to be a locational sensitivity for RILT risk. The multivariable correlation between ipsilateral lung iV₂₀ and RILT, as well as the association of low-functioning V/Q₂₀ and RILT, suggest that irradiating low-functioning regions in the lung may lead to higher toxicity rates.

Scarred Lung. An Update on Radiation-Induced Pulmonary Fibrosis

Radiation-induced pulmonary fibrosis is a common severe long-time complication of radiation therapy for tumors of the thorax. Current therapeutic options used in the clinic include only supportive managements strategies, such as anti-inflammatory treatment using steroids, their efficacy, however, is far from being satisfactory. Recent studies have demonstrated that the development of lung fibrosis is a dynamic and complex process, involving the release of reactive oxygen species, activation of Toll-like receptors, recruitment of inflammatory cells, excessive production of nitric oxide and production of collagen by activated myofibroblasts. In this review we summarized the current state of knowledge on the pathophysiological processes leading to the development of lung fibrosis and we also discussed the possible treatment options.

Dyspnea in Patients Receiving Radical Radiotherapy for Non-Small Cell Lung Cancer: A Prospective Study

Background and Purpose

Dyspnea is an important symptomatic endpoint for assessment of radiation-induced lung injury (RILI) following radical radiotherapy in locally advanced disease, which remains the mainstay of treatment at the time of significant advances in therapy including combination treatments with immunotherapy and chemotherapy and the use of local ablative radiotherapy techniques. We investigated the relationship between dose-volume parameters and subjective changes in dyspnea as a measure of RILI and the relationship to spirometry.

Material and Methods

Eighty patients receiving radical radiotherapy for non-small cell lung cancer were prospectively assessed for dyspnea using two patient-completed tools: EORTC QLQ-LC13 dyspnea quality of life assessment and dyspnea visual analogue scale (VAS). Global quality of life, spirometry and radiation pneumonitis grade were also assessed. Comparisons were made with lung dose-volume parameters.

Results

The median survival of the cohort was 26 months. In the evaluable group of 59 patients there were positive correlations between lung dose-volume parameters and a change in dyspnea quality of life scale at 3 months (V₃₀ p=0.017; V₄₀ p=0.026; V₅₀ p=0.049; mean lung dose p=0.05), and a change in dyspnea VAS at 6 months (V₃₀ p=0.05; V₄₀ p=0.026; V₅₀ p=0.028) after radiotherapy. Lung dose-volume parameters predicted a 10% increase in dyspnea quality of life score at 3 months (V₄₀; p=0.041, V₅₀; p=0.037) and dyspnea VAS score at 6 months (V₄₀; p=0.027) post-treatment.

Conclusions

Worsening of dyspnea is an important symptom of RILI. We demonstrate a relationship between lung dose-volume parameters and a 10% worsening of subjective dyspnea scores. Our findings support the use of subjective dyspnea tools in future studies on radiation-induced lung toxicity, particularly at doses below conventional lung radiation tolerance limits.

Risk factors for symptomatic radiation pneumonitis after stereotactic body radiation therapy (SBRT) in patients with non-small cell lung cancer

BACKGROUND AND PURPOSE

Radiation pneumonitis (RP) can be a potential fatal toxicity of stereotactic body radiation therapy (SBRT) for medically inoperable non-small cell lung cancer (NSCLC). This study aimed to examine the risk factors that predict RP and explore dosimetric tolerance for safe practice in a large institutional series of NSCLC patients.

MATERIALS AND METHODS

Patients with early-stage and locally recurrent NSCLC who received lung SBRT between 2002 and 2015 formed the study population. The primary endpoint was grade 2 or above radiation pneumonitis (RP2). Lungs were re-contoured consistently by one radiation oncologist according to the RTOG atlas for organs at risk. Dosimetric factors were computed consistently with exclusion of gross tumor volume of either ipsilateral, contralateral, or total lungs.

RESULTS

A total of 339 patients were eligible. With a median follow-up of 47 months, RP2 was recorded in 10% patients. History of respiratory comorbidity, previous thoracic radiation, right lung location, mean lung doses of total or ipsilateral lung, and total lung volume receiving 20 Gy were all significantly associated with the risk of RP2. The dosimetric parameters of contralateral lung, including mean dose and volume receiving more than 5, 10, and 20 Gy, were not significantly associated with RP2 (ps > 0.05). A model of combining significant clinical and dosimetric factors had a predictive accuracy AUC of 0.76. According to this model, RP2 can be limited to <10% should the patient have no previous lung radiation and the mean dose of total and ipsilateral lungs be kept less than 6 Gy and 20 Gy, respectively.

CONCLUSION

Dosimetric factors of total or ipsilateral lung together with important clinical factors were significant risk factors for symptomatic radiation pneumonitis after SBRT. Constraining mean lung dose can limit clinically significant lung toxicity.

Dose-sparing effect of deep inspiration breath hold technique on coronary artery and left ventricle segments in treatment of breast cancer

BACKGROUND AND PURPOSE

The risk of radiation-induced cardiac injury remains a challenging problem in the treatment of breast cancer. Certain cardiac structures receive higher doses than others, which results in variable frequencies of radiation-induced injuries across these structures. Radiation dose can be reduced using the deep inspiration breath hold (DIBH) technique. We aimed to investigate the dose reductions from DIBH in individual cardiac segments.

MATERIALS AND METHODS

A dosimetric analysis was performed on left-sided breast cancer patients who underwent breast-conserving surgery and whole breast irradiation. Radiation doses to the cardiac structures were compared between the DIBH and free-breathing (FB) techniques and the dose reductions with DIBH were correlated to the lung expansion.

RESULTS

For the 75 patients included in our study, DIBH effectively reduced doses to the heart, left lung, left anterior descending coronary artery (LAD) and left ventricle (LV), but the degree of dose reductions was variable across different structures. The absolute dose reductions were greatest in the distal LAD (14.4 Gy) and apical LV (12.1 Gy) segments, compared with the other LAD (middle 9.7 Gy, proximal 1.6 Gy) and LV (anterior 5.3 Gy, lateral 2.9 Gy, septal 2.0 Gy, inferior 0.2 Gy) segments. Left lung expansion was significantly correlated with the dose reductions in the LAD (Spearman's rank correlation coefficient, ρ, 0.304) and LV (ρ, 0.420) segments.

CONCLUSIONS

Our study demonstrates the dose-sparing effects of DIBH in various cardiac structures, especially the distal LAD and apical LV segments. The large dose reductions seen in the distal LAD and apical LV segments could potentially translate into clinical benefit of reduced cardiac toxicity, as these structures have been previously shown to receive the highest doses and are associated with radiation-induced injury.

Variations Between Dose-Ventilation and Dose-Perfusion Metrics in Radiation Therapy Planning for Lung Cancer

Purpose

Currently, several active clinical trials of functional lung avoidance radiation therapy using different imaging modalities for ventilation or perfusion are underway. Patients with lung cancer often show ventilation-perfusion mismatch, whereas the significance of dose-function metric remains unclear. The aim of the present study was to compare dose-ventilation metrics with dose-perfusion metrics for radiation therapy plan evaluation.

Methods and Materials

Pretreatment 4-dimensional computed tomography and ^99mTc-macroaggregated albumin single-photon emission computed tomography perfusion images of 60 patients with lung cancer treated with radiation therapy were analyzed. Ventilation images were created using the deformable image registration of 4-dimensional computed tomography image sets and image analysis for regional volume changes as a surrogate for ventilation. Ventilation and perfusion images were converted into percentile distribution images. Analyses included Pearson’s correlation coefficient and comparison of agreements between the following dose-ventilation and dose-perfusion metrics: functional mean lung dose and functional percent lung function receiving 5, 10, 20, 30, and 40 Gy (fV₅, fV₁₀, fV₂₀, fV₃₀, and fV₄₀, respectively).

Results

Overall, the dose-ventilation metrics were greater than the dose-perfusion metrics (ie, fV₂₀, 26.3% ± 9.9% vs 23.9% ± 9.8%). Correlations between the dose-ventilation and dose-perfusion metrics were strong (range, r = 0.94-0.97), whereas the agreements widely varied among patients, with differences as large as 6.6 Gy for functional mean lung dose and 11.1% for fV₂₀. Paired t test indicated that the dose-ventilation and dose-perfusion metrics were significantly different.

Conclusions

Strong correlations were present between the dose-ventilation and dose-perfusion metrics. However, the agreement between the dose-ventilation and dose-perfusion metrics widely varied among patients, suggesting that ventilation-based radiation therapy plan evaluation may not be comparable to that based on perfusion. Future studies should elucidate the correlation of dose-function metrics with clinical pulmonary toxicity metrics.

Hypofractionated Radiotherapy Dose Scheme and Application of New Techniques Are Associated to a Lower Incidence of Radiation Pneumonitis in Breast Cancer Patients

Purpose: Radiation pneumonitis (RP) is one of the most severe toxicities experienced by patients with breast cancer after radiotherapy (RT). RT fractionation schemes and techniques for breast cancer have undergone numerous changes over the past decades. This study aimed to investigate the incidence of RP as a function of such changes in patients with breast cancer undergoing RT and to identify dosimetric markers that predict the risk of this adverse event.

Methods and Materials: We identified 1,847 women with breast cancer who received adjuvant RT at our institution between 2015 and 2017. The RT technique was individually tailored based on each patient's clinicopathological features. Deep inspiration breath hold technique or prone positioning were used for patients who underwent left whole-breast irradiation for cardiac sparing, while those requiring regional lymph node irradiation underwent volumetric-modulated arc therapy (VMAT).

Results: Of 1,847 patients who received RT, 21.2% received the conventional dose scheme, while 78.8% received the hypofractionated dose scheme (mostly 40 Gy in 15 fractions). The median follow-up period was 14.5 months, and the overall RP rate was 2.1%. The irradiated organ at risk was corrected concerning biologically equivalent dose. The ipsilateral lung V₃₀ in equivalent dose in 2 Gy (EQD2) was the most significant dosimetric factor associated with RP development. Administering RT using VMAT, and hypofractionated dose scheme significantly reduced ipsilateral lung V₃₀.

Conclusions: Application of new RT techniques and hypofractionated scheme significantly reduce the ipsilateral lung dose. Our data demonstrated that ipsilateral lung V₃₀ in EQD2 is the most relevant dosimetric predictor of RP in patients with breast cancer.

Voxel-based analysis in radiation oncology: A methodological cookbook

Recently, 2D or 3D methods for dose distribution analysis have been proposed as evolutions of the Dose Volume Histogram (DVH) approaches. Those methods, collectively referred to as pixel- or voxel-based (VB) methods, evaluate local dose response patterns and go beyond the organ-based philosophy of Normal Tissue Complication Probability (NTCP) modelling. VB methods have been introduced in the context of radiation oncology in the very last years following the virtuous example of neuroimaging experience. In radiation oncology setting, dose mapping is a suitable scheme to compare spatial patterns of local dose distributions between patients who develop toxicity and who do not. In this critical review, we present the methods that include spatial dose distribution information for evaluating different toxicity endpoints after radiation therapy. The review addresses two main topics. First, the critical aspects in dose map building, namely the spatial normalization of the dose distributions from different patients. Then, the issues related to the actual dose map comparison, i.e. the viable options for a robust VB statistical analysis and the potential pitfalls related to the adopted solutions. To elucidate the different theoretical and technical issues, the covered topics are illustrated in relation to practical applications found in the existing literature. We conclude the overview on the VB philosophy in radiation oncology by introducing new phenomenological approaches to NTCP modelling that accounts for inhomogeneous organ radiosensitivity.

Acute severe radiation pneumonitis among non-small cell lung cancer (NSCLC) patients with moderate pulmonary dysfunction receiving definitive concurrent chemoradiotherapy: Impact of pre-treatment pulmonary function parameters

PURPOSE

Severe acute radiation pneumonitis (SARP) is a life-threatening complication of thoracic radiotherapy. Pre-treatment pulmonary function (PF) may influence its incidence. We have previously reported on the incidence of SARP among patients with moderate pulmonary dysfunction who received definitive concurrent chemoradiotherapy (dCCRT) for non-small cell lung cancer (NSCLC).

METHODS

The clinical outcomes, dose-volume histograms (DVH), and PF parameters of 122 patients (forced expiratory volume in 1 s [FEV1%]: 60-69%) receiving dCCRT between 2013 and 2019 were recorded. SARP was defined as grade ≥3 RP occurring during or within 3 months after CCRT. Logistic regression, receiver operating characteristics curves (ROC), and hazard ratio (HR) analyses were performed to evaluate the predictive value of each factor for SARP.

RESULTS

Univariate and multivariate analysis indicated that the ratio of carbon monoxide diffusing capacity (DLCO%; odds ratio [OR]: 0.934, 95% confidence interval [CI] 0.896-0.974, p = 0.001) and mean lung dose (MLD; OR: 1.002, 95% CI 1.001-1.003, p = 0.002) were independent predictors of SARP. The ROC AUC of combined DLCO%/MLD was 0.775 (95% confidence interval [CI]: 0.688-0.861, p = 0.001), with a sensitivity and specificity of 0.871 and 0.637, respectively; this was superior to DLCO% (0.656) or MLD (0.667) alone. Compared to the MLD-low/DLCO%-high group, the MLD-high/DLCO%-low group had the highest risk for SARP, with an HR of 9.346 (95% CI: 2.133-40.941, p = 0.003).

CONCLUSION

The DLCO% and MLD may predict the risk for SARP among patients with pre-treatment moderate pulmonary dysfunction who receive dCCRT for NSCLC. Prospective studies are needed to validate our findings.

Comparison of regional lung perfusion response on longitudinal MAA SPECT/CT in lung cancer patients treated with and without functional tissue-avoidance radiation therapy

OBJECTIVE

The effect of functional lung avoidance planning on radiation dose-dependent changes in regional lung perfusion is unknown. We characterized dose-perfusion response on longitudinal perfusion single photon emission computed tomography (SPECT)/CT in two cohorts of lung cancer patients treated with and without functional lung avoidance techniques.

METHODS

The study included 28 primary lung cancer patients: 20 from interventional (NCT02773238) (FLARE-RT) and eight from observational (NCT01982123) (LUNG-RT) clinical trials. FLARE-RT treatment plans included perfused lung dose constraints while LUNG-RT plans adhered to clinical standards. Pre- and 3 month post-treatment macro-aggregated albumin (MAA) SPECT/CT scans were rigidly co-registered to planning four-dimensional CT scans. Tumour-subtracted lung dose was converted to EQD2 and sorted into 5 Gy bins. Mean dose and percent change between pre/post-RT MAA-SPECT uptake (%ΔPERF), normalized to total tumour-subtracted lung uptake, were calculated in each binned dose region. Perfusion frequency histograms of pre/post-RT MAA-SPECT were analyzed. Dose-response data were parameterized by sigmoid logistic functions to estimate maximum perfusion increase (%ΔPERFmaxincrease), maximum perfusion decrease (%ΔPERFmaxdecrease), dose midpoint (Dmid), and dose-response slope (k).

RESULTS

Differences in MAA perfusion frequency distribution shape between time points were observed in 11/20 (55%) FLARE-RT and 2/8 (25%) LUNG-RT patients (p < 0.05). FLARE-RT dose response was characterized by >10% perfusion increase in the 0-5 Gy dose bin for 8/20 patients (%ΔPERFmaxincrease = 10-40%), which was not observed in any LUNG-RT patients (p = 0.03). The dose midpoint Dmid at which relative perfusion declined by 50% trended higher in FLARE-RT compared to LUNG-RT cohorts (35 GyEQD2 vs 21 GyEQD2, p = 0.09), while the dose-response slope k was similar between FLARE-RT and LUNG-RT cohorts (3.1-3.2, p = 0.86).

CONCLUSION

Functional lung avoidance planning may promote increased post-treatment perfusion in low dose regions for select patients, though inter-patient variability remains high in unbalanced cohorts. These preliminary findings form testable hypotheses that warrant subsequent validation in larger cohorts within randomized or case-matched control investigations.

ADVANCES IN KNOWLEDGE

This novel preliminary study reports differences in dose-response relationships between patients receiving functional lung avoidance radiation therapy (FLARE-RT) and those receiving conventionally planned radiation therapy (LUNG-RT). Following further validation and testing of these effects in larger patient populations, individualized estimation of regional lung perfusion dose-response may help refine future risk-adaptive strategies to minimize lung function deficits and toxicity incidence.

Radiation-Induced Lung Injury (RILI)

Radiation pneumonitis (RP) and radiation fibrosis (RF) are two dose-limiting toxicities of radiotherapy (RT), especially for lung, and esophageal cancer. It occurs in 5-20% of patients and limits the maximum dose that can be delivered, reducing tumor control probability (TCP) and may lead to dyspnea, lung fibrosis, and impaired quality of life. Both physical and biological factors determine the normal tissue complication probability (NTCP) by Radiotherapy. A better understanding of the pathophysiological sequence of radiation-induced lung injury (RILI) and the intrinsic, environmental and treatment-related factors may aid in the prevention, and better management of radiation-induced lung damage. In this review, we summarize our current understanding of the pathological and molecular consequences of lung exposure to ionizing radiation, and pharmaceutical interventions that may be beneficial in the prevention or curtailment of RILI, and therefore enable a more durable therapeutic tumor response.

Normal tissue complication probability (NTCP) models for modern radiation therapy

Mathematical models of normal tissue complication probability (NTCP) able to robustly predict radiation-induced morbidities (RIM) play an essential role in the identification of a personalized optimal plan, and represent the key to maximizing the benefits of technological advances in radiation therapy (RT). Most modern RT techniques pose, however, new challenges in estimating the risk of RIM. The aim of this report is to schematically review NTCP models in the framework of advanced radiation therapy techniques. Issues relevant to hypofractionated stereotactic body RT and ion beam therapy are critically reviewed. Reirradiation scenarios for new or recurrent malignances and NTCP are also illustrated. A new phenomenological approach to predict RIM is suggested.

Development and internal validation of a multinomial NTCP model for the severity of acute dyspnea after radiotherapy for lung cancer

BACKGROUND AND PURPOSE

Dyspnea evolution after radiotherapy for lung cancer is complex with potential symptom deterioration and improvement from baseline. We developed and internally validated a multinomial normal tissue complication probability (NTCP) model predicting dyspnea grade.

MATERIALS AND METHODS

Patient-reported dyspnea was collected pre-treatment and during 6 months follow-up for 182 stage I-IV lung cancer patients treated with radical (chemo)radiotherapy. Dyspnea changes (ΔDys) from the baseline grade (Dys₀) to the follow-up grade (Dys) were evaluated. A multinomial logistic regression model simultaneously predicting 3 grades of Dys (Dys ≥ 3, Dys = 2 and Dys ≤ 1 (reference level)) was optimized. Reference NTCP models predicting Dys ≥ 2 and Dys ≥ 3 risks irrespective of Dys₀ were generated for comparison. Models were shrunken and performance was assessed using optimism-corrected AUC (bootstrapping).

RESULTS

Rates of ΔDys ≥ 1 (deterioration) and ΔDys ≤ -1 (improvement) at 6 months were 31.9% and 12.6%. Dys ≥ 3, Dys = 2 and Dys ≤ 1 rates were 13.7%, 20.9% and 65.4%, respectively. The multinomial model (combining the risk factors Dys₀ and MLD and the protective factor chemotherapy treatment) predicted Dys ≥ 3, Dys = 2 and Dys ≤ 1 with AUC (95% CI) of 0.72 (0.65-0.75) 0.76 (0.72-0.79) and 0.78 (0.74-0.80), respectively. Reference Dys ≥ 2 and Dys ≥ 3 models showed worse AUC: 0.64 (0.59-0.67) and 0.66 (0.50-0.70), respectively.

CONCLUSIONS

Dyspnea grade could be predicted with high accuracy using a multinomial NTCP model, yielding personalized dyspnea symptom improvement and deterioration risks.

Machine Learning to Build and Validate a Model for Radiation Pneumonitis Prediction in Patients with Non-Small Cell Lung Cancer

PURPOSE

Radiation pneumonitis is an important adverse event in patients with non-small cell lung cancer (NSCLC) receiving thoracic radiotherapy. However, the risk of radiation pneumonitis grade ≥ 2 (RP2) has not been well predicted. This study hypothesized that inflammatory cytokines or the dynamic changes during radiotherapy can improve predictive accuracy for RP2.

EXPERIMENTAL DESIGN

Levels of 30 inflammatory cytokines and clinical information in patients with stages I-III NSCLC treated with radiotherapy were from our prospective studies. Statistical analysis was used to select predictive cytokine candidates and clinical covariates for adjustment. Machine learning algorithm was used to develop the generalized linear model for predicting risk RP2.

RESULTS

A total of 131 patients were eligible and 17 (13.0%) developed RP2. IL8 and CCL2 had significantly (Bonferroni) lower expression levels in patients with RP2 than without RP2. But none of the changes in cytokine levels during radiotherapy was significantly associated with RP2. The final predictive GLM model for RP2 was established, including IL8 and CCL2 at baseline level and two clinical variables. Nomogram was constructed based on the GLM model. The model's predicting ability was validated in the completely independent test set (AUC = 0.863, accuracy = 80.0%, sensitivity = 100%, specificity = 76.5%).

CONCLUSIONS

By machine learning, this study has developed and validated a comprehensive model integrating inflammatory cytokines with clinical variables to predict RP2 before radiotherapy that provides an opportunity to guide clinicians.

CT-based ventilation imaging in radiation oncology

A form of lung function imaging is emerging that uses phase-resolved four-dimensional CT (4DCT or breath-hold CT) images along with image processing techniques to generate lung function maps that provide a surrogate of lung ventilation. CT-based ventilation (referred to as CT-ventilation) research has gained momentum in Radiation Oncology because many lung cancer patients undergo four-dimensional CT simulation as part of the standard treatment planning process. Therefore, generating CT-ventilation images provides functional information without burdening the patient with an extra imaging procedure. CT-ventilation has progressed from an image processing calculation methodology, to validation efforts, to retrospective demonstration of clinical utility in Radiation Oncology. In particular, CT-ventilation has been proposed for two main clinical applications: functional avoidance radiation therapy and thoracic dose-response assessment. The idea of functional avoidance radiation therapy is to preferentially spare functional portions of the lung (as measured by CT-ventilation) during radiation therapy with the hypothesis that reducing dose to functional portions of the lung will lead to reduced rates of radiation-related thoracic toxicity. The idea of imaging-based dose-response assessment is to evaluate pre- to post-treatment CT-ventilation-based imaging changes. The hypothesis is that early, imaging-change-based response can be an early predictor of subsequent thoracic toxicity. Based on the retrospective evidence, the clinical applications of CT-ventilation have progressed from the retrospective setting to on-going prospective clinical trials. This review will cover basic CT-ventilation calculation methodologies, validation efforts, presentation of clinical applications, summarize on-going clinical trials, review potential uncertainties and shortcomings of CT-ventilation, and discuss future directions of CT-ventilation research.

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

PURPOSE

Radiation pneumonitis (RP) is commonly associated with thoracic radiation therapy, and its incidence is related to dose and volume of the normal lung in the path of radiation. Our aim was to investigate dose patterns associated with RP in patients enrolled in a randomized trial of intensity modulated radiation therapy (IMRT) versus passive scattering proton therapy (PSPT) for locally advanced non-small cell lung cancer.

METHODS

We analyzed 178 patients prospectively treated with PSPT or IMRT for non-small cell lung cancer to a prescribed dose of 66 or 74 Gy in conventional daily fractionation with concurrent chemotherapy. Forty patients (22%) developed clinically symptomatic RP. Voxel-based analysis of local dose differences was done with a nonparametric permutation test accounting for multiple comparisons. From the obtained 3-dimensional significance maps, we derived clusters of voxels that exhibited dose differences between groups at a statistical significance level of 0.05.

RESULTS

The voxel-based analysis highlighted that (1) significant dose differences between patients with and without RP were found in the lower part of the lungs and in the heart and (2) the anatomic regions significantly spared by PSPT and the clusters in which doses were significantly correlated with RP development were disjoint.

CONCLUSIONS

The analyzed trial data provide an unprecedented opportunity to substantiate previous hypotheses regarding the role of the heart and the lower lungs in the development of RP. Knowledge of this relationship between RP and thoracic regional radiosensitivity should be considered in clinical practice and in the design of future trials.

SPECT V/Q in Lung Cancer Radiotherapy Planning

Curative-intent lung cancer radiation therapy either alone (RT) or combined with immuno-chemotherapy is associated with potential risk of serious radiation-induced lung injury. This review provides a summary of the role of SPECT ventilation perfusion (V/Q) imaging as an emerging adjunct to lung cancer RT planning and treatment dosimetry. Denoted "functional lung avoidance RT" it is hypothesized that preferential dosimetric avoidance of physiologically functional lung may reduce the frequency of radiation-induced lung injury. SPECT V/Q imaging datasets available during the planning process allows the prioritization (or "personalization') of RT dose to minimize the volume of functional lung probabilistically exposed to injurious radiation dose. Selective escalation of target dose and adaptive planning and replanning is also enabled. The emergent importance of the tumor-lung microenvironment and its biologic relationship to local immune effectors in lung cancer provides further incentive to individualize RT planning and delivery. This review examines important normal tissue dosimetric constraints that are part of current standards-of-care and the new dosimetric parameters associated with functional lung avoidance RT. SPECT V/Q has been a valuable tool in investigating the feasibility and efficacy of functional lung avoidance RT but is yet to become main stream due to the lack of large clinical trials. It is encouraging however that functional lung avoidance is feasible in RT dose-target delineation and some of the more promising studies are discussed.

[Radiation-induced lung toxicity predictors after stereotactic radiation therapy for non-small cell lung carcinoma stage I]

In case of refusal or contraindication for surgical management of a stage I non-small cell lung carcinoma, the validated alternative therapy is stereotactic irradiation. This technique reaches an equivalent tumour control rate than surgery and significantly higher than conventional radiotherapy. One of the dreaded complications is radiation induced lung toxicity (radiation pneumonitis and lung fibrosis), especially when it is symptomatic, occurring in about 10 % of cases. This article is a literature review of this complication's predictive factors.

Functional lung imaging in radiation therapy for lung cancer: A systematic review and meta-analysis

RATIONALE

Advanced imaging techniques allow functional information to be derived and integrated into treatment planning.

METHODS

A systematic review was conducted with the primary objective to evaluate the ability of functional lung imaging to predict risk of radiation pneumonitis. Secondary objectives were to evaluate dose-response relationships on post treatment functional imaging and assess the utility in including functional lung information into treatment planning. A structured search for publications was performed following PRISMA guidelines and registered on PROSPERO.

RESULTS

814 articles were screened against review criteria and 114 publications met criteria. Methods of identifying functional lung included using CT, MRI, SPECT and PET to image ventilation or perfusion. Six studies compared differences between functional and anatomical lung imaging at predicting radiation pneumonitis. These found higher predictive values using functional lung imaging. Twenty-one studies identified a dose-response relationship on post-treatment functional lung imaging. Nineteen planning studies demonstrated the ability of functional lung optimised planning techniques to spare regions of functional lung. Meta-analysis of these studies found that mean (95% CI) functional volume receiving 20 Gy was reduced by 4.2% [95% CI: 2.3: 6.0] and mean lung dose by 2.2 Gy [95% CI: 1.2: 3.3] when plans were optimised to spare functional lung. There was significant variation between publications in the definition of functional lung.

CONCLUSION

Functional lung imaging may have potential utility in radiation therapy planning and delivery, although significant heterogeneity was identified in approaches and reporting. Recommendations have been made based on the available evidence for future functional lung trials.

Correlation of Regional Lung Ventilation and Gas Transfer to Red Blood Cells: Implications for Functional-Avoidance Radiation Therapy Planning

PURPOSE

To investigate the degree to which lung ventilation and gas exchange are regionally correlated, using the emerging technology of hyperpolarized (HP)-129Xe magnetic resonance imaging (MRI).

METHODS AND MATERIALS

Hyperpolarized-129Xe MRI studies were performed on 17 institutional review board-approved human subjects, including 13 healthy volunteers, 1 emphysema patient, and 3 non-small cell lung cancer patients imaged before and approximately 11 weeks after radiation therapy (RT). Subjects inhaled 1 L of HP-129Xe mixture, followed by the acquisition of interleaved ventilation and gas exchange images, from which maps were obtained of the relative HP-129Xe distribution in three states: (1) gaseous, in lung airspaces; (2) dissolved interstitially, in alveolar barrier tissue; and (3) transferred to red blood cells (RBCs), in the capillary vasculature. The relative spatial distributions of HP-129Xe in airspaces (regional ventilation) and RBCs (regional gas transfer) were compared. Further, we investigated the degree to which ventilation and RBC transfer images identified similar functional regions of interest (ROIs) suitable for functionally guided RT. For the RT patients, both ventilation and RBC functional images were used to calculate differences in the lung dose-function histogram and functional effective uniform dose.

RESULTS

The correlation of ventilation and RBC transfer was ρ = 0.39 ± 0.15 in healthy volunteers. For the RT patients, this correlation was ρ = 0.53 ± 0.02 before treatment and ρ = 0.39 ± 0.07 after treatment; for the emphysema patient it was ρ = 0.24. Comparing functional ROIs, ventilation and RBC transfer demonstrated poor spatial agreement: Dice similarity coefficient = 0.50 ± 0.07 and 0.26 ± 0.12 for the highest-33%- and highest-10%-function ROIs in healthy volunteers, and in RT patients (before treatment) these were 0.58 ± 0.04 and 0.40 ± 0.04. The average magnitude of the differences between RBC- and ventilation-derived functional effective uniform dose, fV20Gy, fV10Gy, and fV5Gy were 1.5 ± 1.4 Gy, 4.1% ± 3.8%, 5.0% ± 3.8%, and 5.3% ± 3.9%, respectively.

CONCLUSION

Ventilation may not be an effective surrogate for true regional lung function for all patients.

Design and development of a nonrigid phantom for the quantitative evaluation of DIR-based mapping of simulated pulmonary ventilation

PURPOSE

The validation of deformable image registration (DIR)-based pulmonary ventilation mapping is time consuming and prone to inaccuracies and is also affected by deformation parameters. In this study, we developed a nonrigid phantom as a quality assurance (QA) tool that simulates ventilation to evaluate DIR-based images quantitatively.

METHODS

The phantom consists of an acrylic cylinder filled with polyurethane foam designed to simulate pulmonic alveoli. A polyurethane membrane is attached to the inferior end of the phantom to simulate the diaphragm. In addition, tracheobronchial-tree-shaped polyurethane tubes are inserted through the foam and converge outside the phantom to simulate the trachea. Solid polyurethane is also used to model arteries, which closely follow the model airways. Two three-dimensional (3D) CT scans were performed during exhalation and inhalation phases using xenon (Xe) gas as the inhaled contrast agent. The exhalation 3D-CT image is deformed to an inhalation 3D-CT image using our in-house program based on the NiftyReg open-source package. The target registration error (TRE) between the two images was calculated for 16 landmarks located in the simulated lung volume. The DIR-based ventilation image was generated using Jacobian determinant (JD) metrics. Subsequently, differences in the Hounsfield unit (HU) values between the two images were measured. The correlation coefficient between the JD and HU differences was calculated. In addition, three 4D-CT scans are performed to evaluate the reproducibility of the phantom motion and Xe gas distribution.

RESULTS

The phantom exhibited a variety of displacements for each landmark (range: 1-20 mm). The reproducibility analysis indicated that the location differences were <1 mm for all landmarks, and the HU variation in the Xe gas distribution was close to zero. The mean TRE in the evaluation of spatial accuracy according to the DIR software was 1.47 ± 0.71 mm (maximum: 2.6 mm). The relationship between the JD and HU differences had a large correlation (R = -0.71) for the DIR software.

CONCLUSION

The phantom implemented new features, namely, deformation and simulated ventilation. To assess the accuracy of the DIR-based mapping of the simulated pulmonary ventilation, the phantom allows for simulation of Xe gas wash-in and wash-out. The phantom may be an effective QA tool, because the DIR algorithm can be quickly changed and its accuracy evaluated with a high degree of precision.

Functional-guided radiotherapy using knowledge-based planning

BACKGROUND AND PURPOSE

There are two significant challenges when implementing functional-guided radiotherapy using 4DCT-ventilation imaging: (1) lack of knowledge of realistic patient specific dosimetric goals for functional lung and (2) ensuring consistent plan quality across multiple planners. Knowledge-based planning (KBP) is positioned to address both concerns.

MATERIAL AND METHODS

A KBP model was created from 30 previously planned functional-guided lung patients. Standard organs at risk (OAR) in lung radiotherapy and a ventilation contour delineating areas of high ventilation were included. Model validation compared dose-metrics to standard OARs and functional dose-metrics from 20 independent cases that were planned with and without KBP.

RESULTS

A significant improvement was observed for KBP optimized plans in V20Gy and mean dose to functional lung (p = 0.005 and 0.001, respectively), V20Gy and mean dose to total lung minus GTV (p = 0.002 and 0.01, respectively), and mean doses to esophagus (p = 0.005).

CONCLUSION

The current work developed a KBP model for functional-guided radiotherapy. Modest, but statistically significant, improvements were observed in functional lung and total lung doses.

Anatomic, functional and molecular imaging in lung cancer precision radiation therapy: treatment response assessment and radiation therapy personalization

This article reviews key imaging modalities for lung cancer patients treated with radiation therapy (RT) and considers their actual or potential contributions to critical decision-making. An international group of researchers with expertise in imaging in lung cancer patients treated with RT considered the relevant literature on modalities, including computed tomography (CT), magnetic resonance imaging (MRI) and positron emission tomography (PET). These perspectives were coordinated to summarize the current status of imaging in lung cancer and flag developments with future implications. Although there are no useful randomized trials of different imaging modalities in lung cancer, multiple prospective studies indicate that management decisions are frequently impacted by the use of complementary imaging modalities, leading both to more appropriate treatments and better outcomes. This is especially true of ¹⁸F-fluoro-deoxyglucose (FDG)-PET/CT which is widely accepted to be the standard imaging modality for staging of lung cancer patients, for selection for potentially curative RT and for treatment planning. PET is also more accurate than CT for predicting survival after RT. PET imaging during RT is also correlated with survival and makes response-adapted therapies possible. PET tracers other than FDG have potential for imaging important biological process in tumors, including hypoxia and proliferation. MRI has superior accuracy in soft tissue imaging and the MRI Linac is a rapidly developing technology with great potential for online monitoring and modification of treatment. The role of imaging in RT-treated lung cancer patients is evolving rapidly and will allow increasing personalization of therapy according to the biology of both the tumor and dose limiting normal tissues.

Functional lung avoidance and response-adaptive escalation (FLARE) RT: Multimodality plan dosimetry of a precision radiation oncology strategy

PURPOSE

Nonsmall cell lung cancer (NSCLC) patient radiation therapy (RT) is planned without consideration of spatial heterogeneity in lung function or tumor response. We assessed the dosimetric and clinical feasibility of functional lung avoidance and response-adaptive escalation (FLARE) RT to reduce dose to [^99m Tc]MAA-SPECT/CT perfused lung while redistributing an escalated boost dose within [¹⁸ F]FDG-PET/CT-defined biological target volumes (BTV).

METHODS

Eight stage IIB-IIIB NSCLC patients underwent FDG-PET/CT and MAA-SPECT/CT treatment planning scans. Perfused lung objectives were derived from scatter/collimator/attenuation-corrected MAA-SPECT uptake relative to ITV-subtracted lung to maintain < 20 Gy mean lung dose (MLD). Prescriptions included 60 Gy to the planning target volume (PTV) and concomitant boost of 74 Gy mean to biological target volumes (BTV = GTV + PET gradient segmentation) scaled to each BTV voxel by relative FDG-PET SUV. Dose-painting-by-numbers prescriptions were integrated into commercial treatment planning systems via uptake threshold discretization. Dose constraints for lung, heart, cord, and esophagus were defined. FLARE RT plans were optimized with volumetric modulated arc therapy (VMAT), proton pencil beam scanning (PBS) with 3%-3 mm robust optimization, and combination of PBS (avoidance) plus VMAT (escalation). The high boost dose region was evaluated within a standardized SUV_peak structure. FLARE RT plans were compared to reference VMAT plans. Linear regression between radiation dose to BTV and normalized FDG PET SUV at every voxel was conducted, from which Pearson linear correlation coefficients and regression slopes were extracted. Spearman rank correlation coefficients were estimated between radiation dose to lung and normalized SPECT uptake. Dosimetric differences between treatment modalities were evaluated by Friedman nonparametric paired test with multiple sampling correction.

RESULTS

No unacceptable violations of PTV and normal tissue objectives were observed in 24 FLARE RT plans. Compared to reference VMAT plans, FLARE VMAT plans achieved a higher mean dose to BTV (73.7 Gy 98195. 61.3 Gy), higher mean dose to SUV_peak (89.7 Gy vs. 60.8 Gy), and lower mean dose to highly perfused lung (7.3 Gy vs. 14.9 Gy). These dosimetric gains came at the expense of higher mean heart dose (9.4 Gy vs. 5.8 Gy) and higher maximum cord dose (50.1 Gy vs. 44.6 Gy) relative to the reference VMAT plans. Between FLARE plans, FLARE VMAT achieved higher dose to the SUV_peak ROI than FLARE PBS (89.7 Gy vs. 79.2 Gy, P = 0.01), while FLARE PBS delivered lower dose to lung than FLARE VMAT (11.9 Gy vs. 15.6 Gy, P < 0.001). Voxelwise linear dose redistribution slope between BTV dose and FDG PET uptake was higher in magnitude for FLARE PBS + VMAT (0.36 Gy per %SUV_max ) compared to FLARE VMAT (0.27 Gy per %SUV_max ) or FLARE PBS alone (0.17 Gy per %SUV_max ).

CONCLUSIONS

FLARE RT is clinically feasible with VMAT and PBS. A combination of PBS for functional lung avoidance and VMAT for FDG PET dose escalation balanced target and normal tissue objective tradeoffs. These results provide a technical platform for testing of FLARE RT safety and efficacy within a precision radiation oncology trial.

A complete 4DCT-ventilation functional avoidance virtual trial: Developing strategies for prospective clinical trials

Introduction

4DCT‐ventilation is an exciting new imaging modality that uses 4DCT data to calculate lung‐function maps. Because 4DCTs are acquired as standard of care for lung cancer patients undergoing radiotherapy, 4DCT‐ventiltation provides functional information at no extra dosimetric or monetary cost to the patient. The development of clinical trials is underway to use 4DCT‐ventilation imaging to spare functional lung in patients undergoing radiotherapy. The purpose of this work was to perform a virtual trial using retrospective data to develop the practical aspects of a 4DCT‐ventilation functional avoidance clinical trial.

Methods

The study included 96 stage III lung cancer patients. A 4DCT‐ventilation map was calculated using the patient's 4DCT‐imaging, deformable registration, and a density‐change‐based algorithm. Clinical trial inclusion assessment used quantitative and qualitative metrics based on the patient's spatial ventilation profile. Clinical and functional plans were generated for 25 patients. The functional plan aimed to reduce dose to functional lung while meeting standard target and critical structure constraints. Standard and dose‐function metrics were compared between the clinical and functional plans.

Results

Our data showed that 69% and 59% of stage III patients have regional variability in function based on qualitative and quantitative metrics, respectively. Functional planning demonstrated an average reduction of 2.8 Gy (maximum 8.2 Gy) in the mean dose to functional lung.

Conclusions

Our work demonstrated that 60–70% of stage III patients would be eligible for functional planning and that a typical functional lung mean dose reduction of 2.8 Gy can be expected relative to standard clinical plans. These findings provide salient data for the development of functional clinical trials.

Reliability of dose volume constraint inference from clinical data

Dose volume histogram points (DVHPs) frequently serve as dose constraints in radiotherapy treatment planning. An experiment was designed to investigate the reliability of DVHP inference from clinical data for multiple cohort sizes and complication incidence rates. The experimental background was radiation pneumonitis in non-small cell lung cancer and the DVHP inference method was based on logistic regression. From 102 NSCLC real-life dose distributions and a postulated DVHP model, an 'ideal' cohort was generated where the most predictive model was equal to the postulated model. A bootstrap and a Cohort Replication Monte Carlo (CoRepMC) approach were applied to create 1000 equally sized populations each. The cohorts were then analyzed to establish inference frequency distributions. This was applied to nine scenarios for cohort sizes of 102 (1), 500 (2) to 2000 (3) patients (by sampling with replacement) and three postulated DVHP models. The Bootstrap was repeated for a 'non-ideal' cohort, where the most predictive model did not coincide with the postulated model. The Bootstrap produced chaotic results for all models of cohort size 1 for both the ideal and non-ideal cohorts. For cohort size 2 and 3, the distributions for all populations were more concentrated around the postulated DVHP. For the CoRepMC, the inference frequency increased with cohort size and incidence rate. Correct inference rates  >[Formula: see text] were only achieved by cohorts with more than 500 patients. Both Bootstrap and CoRepMC indicate that inference of the correct or approximate DVHP for typical cohort sizes is highly uncertain. CoRepMC results were less spurious than Bootstrap results, demonstrating the large influence that randomness in dose-response has on the statistical analysis.