Introduction
Radiation pneumonitis (RP) is one of the clinically most common toxicities of radiation therapy for lung cancer. Symptomatic RP has been reported in about 20% of patients with locally advanced non-small cell lung cancer (NSCLC). Reference Giuranno, Ient, De Ruysscher and Vooijs1,Reference Marks, Bentzen and Deasy2 For example, RP was noticed in 20/71 (28%) NSCLC patients treated with concurrent radiochemotherapy using three-dimensional conformal radiotherapy (3DCRT). Reference Tsujino, Hirota and Endo3 In another study evaluating 576 NSCLC patients treated with 3DCRT or intensity-modulated radiotherapy (IMRT) with or without chemotherapy, 117/576 (20%) patients developed symptomatic RP. Reference Jin, Tucker and Liu4 Many factors have been investigated for their association with RP risk. In the QUANTEC report, mean lung dose (MLD) and the per cent lung volume receiving ≥20 Gy (lung V20 Gy) have been established as reproducible dosimetrical predictors. Reference Marks, Bentzen and Deasy2 In a multivariate analysis evaluating clinical and treatment-related factors, V20 Gy was the strongest predictor of RP. Reference Tsujino, Hirota and Endo3 The 6-month cumulative incidence of RP ≥ G2 was 14% in patients with V20 ≤ 25 and 63% in patients with V20 ≥ 26%. Reference Tsujino, Hirota and Endo3 Other clinical and tumour characteristics, including age, gender, pulmonary function, smoking status, chemotherapy, tumour size and location, have also commonly been linked to RP risk, but findings were less consistent. Reference Giuranno, Ient, De Ruysscher and Vooijs1 The report by Jin et al. on 576 irradiated NSCLC patients identified smoking status as the only significant predictor for RP in addition to lung dose. Reference Jin, Tucker and Liu4 Never smokers had a RP risk of 37% compared to 14% for current smokers and 23% for former smokers. In this study, neither age, nor sex, performance status, the application or type of chemotherapy nor pulmonary function were related to RP development. Reference Jin, Tucker and Liu4 In contrast, the study by Robnett et al. reported on 147 patients treated with definitive radiochemotherapy. Reference Robnett, Machtay, Vines, Mckenna, Algazy and Mckenna5 Performance status was identified as the strongest predictor of RP with a 16 versus 2% risk of severe RP for performance status 1 versus 0. Reference Robnett, Machtay, Vines, Mckenna, Algazy and Mckenna5 In this study, sex (higher risk for women) and lung function (higher risk for forced expiratory volume in 1 second (FEV1) <2l) were also significantly associated with RP risk. Review reports on the RP risk have summarised the findings for different parameters. Reference Giuranno, Ient, De Ruysscher and Vooijs1,Reference Marks, Bentzen and Deasy2 While the importance of different parameters can vary between studies, dosimetrical findings have been more consistent and led to commonly accepted recommendations for radiation treatment planning. Reference Marks, Bentzen and Deasy2
Several factors have recently found increasing interest for their potential association with RP. For example, emphysema and interstitial lung disease (ILD), particularly idiopathic lung fibrosis, were identified as strong predictors for RP following stereotactic lung radiotherapy with RP ≥ G3 observed in 32% (9/28) ILD patients. Reference Bahig, Filion and Vu6 The significance of these conditions in the setting of conventionally fractionated radiotherapy is less well known. Reference Lee, Kim and Lee7
An association of RP with metabolic diseases and their treatment, such as diabetes or statin use for hyperlipidaemia, has rarely been reported. Reference Kong, Lim and Kim8,Reference Williams, Hernady and Johnston9 The pathophysiologic mechanism of diabetes mellitus causing inflammation and endothelial dysfunction is presumed to exacerbate the inflammatory effects of radiotherapy on lung tissue and thereby increase RP risk. Reference Kong, Lim and Kim8 A higher RP incidence with diabetes mellitus was only recently endorsed in clinical investigations. Reference Kong, Lim and Kim8 Statins were found to have an anti-inflammatory effect in lung tissue in animal studies. Reference Williams, Hernady and Johnston9 The effect of statins on a reduced RP risk in humans requires further evaluation.
Also, cardiac radiation dose at various levels, including mean heart dose (MHD) and heart V5 Gy–V40 Gy, is currently at the centre of investigations for its effect on overall survival. Reference Chun, Hu and Choy10,Reference Thor, Deasy and Hu11 A connection between cardiac dose and RP has been suggested but is less clear. Reference Huang, Hope and Lindsay12,Reference Shepherd, Iocolano and Leeman13 The magnitude of irradiated heart volumes may be important for the manifestation of RP.
Similarly, lung volumes vary considerably between patients. An earlier study investigated the impact of lung volumes spared from defined dose levels on the development of RP and observed a higher incidence of RP in patients with smaller spared lung volumes. Reference Briere, Krafft, Liao and Martel14 In addition, the assessment of low-dose lung volumes such as lung V5 Gy on RP risk has been controversial. Lung V5 Gy was not significantly related to RP in a large randomised study. Reference Chun, Hu and Choy10 Other observations found detrimental lung toxicity in patients with large low-dose lung volumes. Reference Shepherd, Iocolano and Leeman13,Reference McFarlane, Hochstedler and Laucis15 Further evaluation of these findings appears warranted and thus gives rise to this study.
An increasing number of investigations into RP have been including analyses of larger numbers of clinical and demographic factors to weigh the importance of factors from different categories against each other. Reference Shepherd, Iocolano and Leeman13,Reference McFarlane, Hochstedler and Laucis15–Reference Palma, Senan and Tsujino17 Given the wide spectrum of factors that might affect the RP risk, in this study we investigated a large number of parameters (clinical and radiographic characteristics, tumour features, radiation dose and technique parameters) that are either well-established or have rarely been analysed. Some factors analysed in our study were only recently found to be of interest with respect to RP development, but have not been investigated so far, such as heart volume and absolute volumes for different heart dose levels. In addition to determining risk factors for RP ≥ G3, we planned to also investigate the influence of RP ≥ G3 on overall survival in our patient cohort. A detrimental effect of RP on overall survival was observed in only few previous studies. Reference Beukema, Kawaguchi and Sijtsema18–Reference Shen, Liu and Jin20
Materials/Methods
Patient characteristics
All consecutive patients with locally advanced NSCLC who underwent conventionally fractionated radiotherapy and concomitant chemotherapy (no immunotherapy) at our institution between 2009 and beginning of 2018 were included in this analysis. Patients were staged with chest CT, FDG-18 PET-CT, and MRI or CT of the brain had pathologically confirmed diagnosis of NSCLC and underwent pulmonary function testing. Follow-up visits were performed every 3 months for the first 2 years and every 6 months for next 3 years followed by annual visits. All visits included a CT chest and additional diagnostic tests as clinically indicated. Tumour control, acute and long-term treatment-related side effects were assessed during weekly under-treatment visits and on follow-up visits.
Treatment planning
For treatment planning, patients underwent a 4D CT (Brilliance Big BoreTM, Philips Medical Systems, Andover, MA) using 3 mm slice thickness, 512 × 512 axial resolution and images sorted in 10 breathing phase bins. Depending on the available image management software at the time, tumour delineation was performed using an internal gross tumour volume (iGTV) that was created based on GTV propagation or based on a reconstructed maximum intensity projection image using commercial software (MIM MaestroTM, MIM, Cleveland, OH). Patients with >1·0 cm respiratory tumour motion underwent three repeated moderate inspiration breath hold CT scans with active breathing control (ABCTM, Elekta, Stockholm, Sweden) from which an iGTV was generated. For both scenarios, iGTV to clinical target volume (CTV) expansions were 0·6–0·8 cm. Planning target volumes (PTVs) were created through CTV expansion by 0·5 cm. Treatment planning included heterogeneity correction (PinnacleTM, Philips Radiation Oncology Systems, Fitchburg, WI or EclipseTM, Varian, Palo Alto, CA). Plans were optimised to achieve 95% PTV coverage by the prescription dose. Normal tissue dose constraints were unchanged over the years and included the following parameters: spinal cord Dmax ≤ 45 Gy, brachial plexus Dmax ≤ 66 Gy, MLD ≤ 20 Gy, lungs V20 Gy ≤ 30%, lungs V30 Gy ≤ 20%, mean oesophagus dose ≤ 34 Gy, oesophagus V60 Gy ≤ 17%, oesophagus Dmax ≤ 105% of prescription dose, heart V60 Gy ≤ 30%, heart V45 Gy ≤ 60%, MHD ≤ 35 Gy. Treatment was performed typically with 6 MV photon beams with image guidance including daily planar imaging and weekly cone beam computed tomography. Treatment-specific characteristics are listed in Table 1.
Table 1. Treatment and radiation plan characteristics
All lung doses are for Lungs-CTV.
3DCRT: three-dimensional conformal radiotherapy; cyc: cycles; IMRT: intensity-modulated radiotherapy; IQR: interquartile range; PTV: planning target volume; VMAT: volumetric modulated arc therapy.
Analysis
Based on a comprehensive database that was created from electronic medical charts and radiation therapy documentation, a retrospective IRB-approved analysis was performed to identify predictors of RP ≥ G3 including 25 clinical (gender, race, diffusing capacity for carbon monoxide—DLCO, diabetes type II, statin use, smoking history), radiographic (emphysema, ILD on pre-treatment CT) and radiotherapy dose- and technique-related factors. The latter comprised PTV, tumour location, treatment technique, lung volume (mL), MLD (Gy), lung V5 Gy (%), lung Vspared5 Gy (mL), lung V20 Gy (%), lung Vspared20 Gy (mL), heart volume (mL), MHD (Gy), heart V5 Gy/30 Gy/60 Gy (each in % and mL). Factors selected for analysis were chosen according to previously identified associations with RP risk or were newly investigated such as heart doses in ml and heart volume. For this analysis, all diagnostic pre-treatment images were analysed by a thoracic radiologist for the presence of ILD according to Fleischner Society criteria. Reference Lynch, Sverzellati and Travis21 In addition, the extent of emphysema involvement of the lung was quantified based on visual assessment of percentage lung involvement. For this analysis, emphysema affecting ≤25% of the lung volume was compared to >25% involvement.
RP was graded according to Common Terminology Criteria for Adverse Events v5. Scoring of low-grade RP can be subjective; RP ≥ G3 was therefore chosen as a more clearly defined and clinically significant endpoint for this study based on severity of symptoms and the need for steroid treatment and oxygen supplement. Survival data were obtained from the patients’ medical records and the local cancer centre tumour registry.
Statistics
Statistical analysis was performed to identify the predictive power of clinical and dosimetrical parameters towards RP. Initially, individual variables were checked using univariate logistic regression (UVA) to find marginal associations. Primary model building for association with RP ≥ G3 on this limited-size cohort set was performed using a stepwise forward logistic regression algorithm (multivariate predictive model (MVA)). Entry threshold for predictive variables was set at 10%. Many of the predictors demonstrated collinearity, especially the dosimetric factors. The multivariate forward selection method was expected to eliminate such predictors in the course of variable selection. Overall survival times were computed using Kaplan–Meier method. All analyses were performed in SAS 9.4. To demonstrate the predictive value of individual variables that showed borderline significance in MVA, receiver operating characteristic curves were created to assess accuracy of prediction. Analysis of differences between treatment techniques was performed with unpaired two-tailed t-tests.
Results
Population
One hundred and five patients with complete follow-up and clinical/radiographic data were included in this study. The patient cohort included inoperable American Joint Committee on Cancer stage IIB and III NSCLC; two-stage IV patients with solitary oligometastases outside the chest treated with curative approach were included as well. All except three had a positive smoking history or were current smokers. Patient-specific characteristics are listed in Table 2.
Table 2. Baseline patient demographics and disease characteristics (n = 105)
DLCO: diffusion capacity for carbon monoxide; FEV1: forced expiratory volume in 1 second; IQR: interquartile range.
Treatment characteristics
Patients were treated with daily conventionally fractionated radiotherapy to a median dose of 64·0 Gy with 1·8/2·0 Gy single fraction dose. Initially, 3DCRT was used in 38 (36%) patients using typically four or five static fields. Later, treatment technique was changed to IMRT delivered in 54 (51%) patients using step-and-shoot technique with five to seven fields. More recently, VMAT was employed in 14 (13%) patients using typically two to four arcs. See Table 1 for relevant radiotherapy-related parameters.
RP
RP ≥ G3 was diagnosed in 10/105 (9·5%) patients (seven RP G3, one RP G4, two RP G5). On UVA, predictors for RP ≥ G3 were diabetes, tumour location in the upper lobe, PTV, VMAT, lung V5 Gy (%), lung Vspared5 Gy (mL), lung V20 Gy (%), heart V5 Gy (% and mL). On MVA, VMAT was the only significant predictor for RP ≥ G3 (p = 0·042), see Table 3. Lung V5 Gy (%, p = 0·073) and lung V20 Gy (%, p = 0·08) were borderline significant for predicting RP ≥ G3 but showed high prediction accuracy with AUC values of 0·86 and 0·84, respectively (Figure 1).
Table 3. Results of UVA and MVA for RP ≥ G3
Cursive text indicates borderline significance.
95% CI: 95% confidence interval; 3DCRT: three-dimensional conformal radiotherapy; IMRT: intensity-modulated radiotherapy; MVA: multivariate analysis; OR: odds ratio; PTV: planning target volume; UVA: univariate analysis; VMAT: volumetric modulated arc therapy.
Figure 1. Receiver operating characteristic (ROC) analysis for RP ≥ G3. Lung V5 Gy and lung V20 Gy were borderline significantly associated with RP ≥ G3. Both parameters show high prediction accuracy.
Overall survival
Median follow-up of surviving patients was 33 months (9–132). The 3-year OS rate was 43·4%, and median survival was 28·5 months. Patients with RP ≥ 3 had a median survival of 10·0 months compared to 29·5 months with RP < 3 or no RP (p = 0·02). Corresponding 2-year OS rates were 26·7 and 56·9% (Figure 2).
Figure 2. Overall survival depending on radiation pneumonitis grade. Patients with RP ≥ G3 had significantly worse survival compared to patients with lower grade or no RP, p = 0·02.
Discussion
High-grade RP remains the most detrimental side effect of radiotherapy to the chest. In our cohort, RP ≥ G3 was observed in approximately 10%, similar to RP rates reported in other publications. Reference Huang, Hope and Lindsay12,Reference McFarlane, Hochstedler and Laucis15 Univariate logistic regression confirmed several established risk factors, such as tumour location, PTV size and lung V20 Gy (%).
Less often reported factors were also significant on UVA. Patients with diabetes mellitus were found to have a higher RP risk which agrees with few other reports. Reference Kong, Lim and Kim8,Reference Kalman, Hugo, Mahon, Deng, Mukhopadhyay and Weiss22 Kong et al. reported that the diagnosis of diabetes mellitus, HbA1c>6·15% and fasting glucose levels >121 mg/dL were all associated with increasing RP ≥ G3 by a factor>2 on MVA. Reference Kong, Lim and Kim8 Diabetes mellitus, particularly if poorly controlled, warrants further confirmation as a risk factor for RP.
In addition, UVA indicated that the absolute lung Vspared5 Gy was significantly related to RP ≥ G3. Biere et al. also observed the importance of sparing uninvolved lung volumes. Reference Briere, Krafft, Liao and Martel14 They reported that patients with smaller lung volumes spared (often female patients with smaller lungs) had a higher risk of RP. Total lung volumes in our cohort varied by a factor of five. Based on these observations, there might therefore be a minimum lung volume that should not receive any radiation dose to prevent high-grade RP. Patients with larger lung volumes might have a smaller RP risk because larger parts of their lungs can be spared from any radiation.
In our cohort, heart V5 Gy (% and mL) was associated with RP ≥ G3 on UVA. So far, little information is available on the relation between radiation dose to the heart and RP. Shepherd et al. observed that heart V16 Gy was related to RP ≥ G2 in postoperative radiotherapy for lung cancer, Reference Shepherd, Iocolano and Leeman13 whereas Huang et al. found heart V65 Gy associated with RP. Reference Huang, Hope and Lindsay12 While synergistic effects between lung and heart response to radiation have been discussed, higher lung doses are often found with lower lobe tumour locations leading also to higher heart doses because of the proximity of lung tumours to the heart. On the contrary, Wijsman et al. using IMRT and VMAT did not observe a correlation between incidental heart dose and RP risk, potentially because of the ability to selectively spare heart with modern conformal radiotherapy techniques. Reference Wijsman, Dankers and Troost23 Heart volumes in our study varied by a factor of more than four indicating that there might be large variations in the absolute heart volumes receiving certain dose levels. Further investigation into the importance of cardiac DVH parameters in addition to the relevance of functional subvolumes of the heart appears warranted.
On MVA, both lung V5 Gy and V20 Gy were borderline significant with high AUC values >0·8. While lung V20 Gy is one of the accepted standard volume thresholds for RP, lung V5Gy has been more controversial. Reference Marks, Bentzen and Deasy2 Several reports indicated that low-dose lung volumes, such as V5Gy or V10Gy, are associated with an increased RP risk. Reference Shepherd, Iocolano and Leeman13,Reference Khalil, Hoffmann, Moeller, Farr and Knap24–Reference Song, Pyo, Moon, Kim, Kim and Cho26 Contrary to a large randomised trial that did not identify lung V5Gy as a predictor of RP, a recent large prospective study with >1300 patients confirmed the importance of lung V5 Gy. Reference Chun, Hu and Choy10 Both lung V5 Gy and V20 Gy were significantly associated with RP ≥ G2, whereas for RP ≥ G3 MLD and V20 Gy were found significant. Reference McFarlane, Hochstedler and Laucis15 Similarly, in a recent analysis of IMRT-treated NSCLC patients, low and intermediate lung dose levels were associated with RP ≥ G2, but MLD, V20 Gy and particularly V30 Gy were identified as the strongest predictors of RP compared to clinical and low-dose factors. Reference Meng, Luo and Wang27 In the absence of lung V5 Gy dose constraints in current NCCN guidelines, keeping lung V5 Gy volumes lower than 60 to 75% has been shown to reduce RP risk. Reference McFarlane, Hochstedler and Laucis15,Reference Khalil, Hoffmann, Moeller, Farr and Knap24
The only significant predictor of RP ≥ G3 on MVA in our cohort was VMAT. Over the last two decades, complex treatment techniques have evolved which have enabled superior target dose conformality while reducing dose to critical normal tissues resulting in improved tumour control, survival and lower treatment-related toxicity. Reference Chun, Hu and Choy10,Reference Liao, Komaki and Thames28,Reference Jegadeesh, Liu and Gillespie29 VMAT has been used for stereotactic lung cancer treatments for several years but is also increasingly used for locally advanced lung cancer radiotherapy. Using one or more arcs, VMAT achieves highly conformal dose distributions at the target, which can be at the expense of large volumes receiving low doses if unconstrained. In a planning study, Jiang et al. compared IMRT with single and partial arc VMAT without using low lung dose constraints. Reference Jiang, Li and Liu30 Higher lung V5 Gy/10 Gy and lower lung V20 Gy/30 Gy doses were observed for both single and partial arcs compared to IMRT for total and contralateral lung. In a clinical analysis, Wijsmans et al. found lower lung V5 Gy with VMAT compared to IMRT without using low lung dose constraints. Reference Wijsman, Dankers and Troost31 Although lung V5 Gy was lower with VMAT in this study, interestingly grade 4 and 5 toxicities were only observed in the VMAT group. In addition, the number of patients with late RP ≥ G3 nearly doubled in the VMAT group. Other reports also did not find increased low doses to the lungs or increased toxicity with VMAT either with (lung V5 Gy < 65% and V10 Gy < 45%)Reference Bourbonne, Delafoy, Lucia, Quéré, Pradier and Schick 32 or without using explicit low lung dose constraints. Reference Bertelsen, Hansen and Brink33 In our cohort, RP ≥ 3 occurred only after IMRT (four G3, one G4) and VMAT treatments (three G3, two G5). No particular low-dose lung constraints were used in our study resulting in a higher likelihood for large low-dose volumes with IMRT and even more so with VMAT.
In our study, high-grade RP was significantly associated with decreased 2-year survival by a factor of 2·13. A three-fold increase in median survival was observed in patients without high-grade RP. RP as a risk factor for overall survival was observed in few studies previously. Reference Beukema, Kawaguchi and Sijtsema18–Reference Shen, Liu and Jin20 Inoue et al. observed significantly worse survival in lung cancer patients developing severe RP compared to no or mild RP following delivery of definitive radiotherapy. Reference Inoue, Kunitoh, Sekine, Sumi, Tokuuye and Saijo19 Beukema et al. observed that oesophagus cancer patients treated with 3D CRT had significantly worse survival if they developed RP. Reference Beukema, Kawaguchi and Sijtsema18 Following VMAT, Shen et al. observed that RP in NSCLC patients was significantly related to overall survival (HR 1·39). Reference Shen, Liu and Jin20 Only one patient developed RP G5 making the effect of RP-related death as a factor for worse survival less likely. It might not be the development of RP itself, but other risk factors associated with developing RP that influence survival rates. Reference Beukema, Kawaguchi and Sijtsema18 This observation might also apply to our patients, where VMAT patients who were found to have higher RP ≥ G3 risk had larger PTV, an indicator for worse outcome.
Our study has several limitations including small cohort size. We included only patients who had complete datasets and a minimum follow-up of 9 months for surviving patients to allow for observation of RP development including assessment of RP resolution after respective treatment. To avoid underreporting of toxicities and consistent grading of toxicities, patients followed in outside institutions were excluded. High-grade RP (RP ≥ G3) was selected as the endpoint of this study to allow more reliable classification compared to lower-grade toxicity. Also, low rates of high-grade RP result in low statistical power. As a consequence, the authors selected p < 0·1 as a statistically significant cut-off in this study as it is often the case in similar clinical studies. A large variety of parameters were included in this risk factor analysis. While several dose-volume parameters are likely correlated, univariate logistic regression was employed to determine the significance of each factor individually prior to multivariate analysis. An imbalance in RP events is inherent to a cohort with limited sample size. Despite low toxicity rates and small sample size, various standard RP-related risk factors (PTV size, lung V20 Gy) were identified as significant on UVA validating our approach. This study did not investigate chemotherapy, a well-known risk factor for RP, as chemotherapy agents did not change over time.
The analysis confirmed several well-established and less often reported risk factors for high-grade RP. Further investigations into low-dose lung volumes, lung volumes spared radiation dose and VMAT are warranted given this study’s findings are observed in a relatively small cohort.
Conclusions
In this study, VMAT was the only factor that was significantly correlated with RP ≥ G3. Lung V5 Gy and lung V20 Gy were reliable predictors of RP ≥ G3 with AUC values >0·8. Overall survival was significantly worse in patients who experienced RP ≥ G3. Avoiding RP ≥ G3 is therefore important not only as a toxicity per se but also as a risk factor for poor survival. When using VMAT, caution needs to be taken to reduce low-dose lung volumes in addition to other well-established dose constraints.
Acknowledgements
None.
Financial Support
Biostatistical analysis of this work was supported in part by the Massey Cancer Center institutional grant through the National Institutes of Health (grant number P30CA016059).
Conflicts of Interest
The authors declare none.
Ethical Standards
This retrospective study was approved by the local ethics committee.
Open access
