Longitudinal analyses and predictive factors of radiation-induced lung toxicity-related parameters after stereotactic radiotherapy for lung cancer

Abstract

Background and purpose

The purpose of this prospective study was to investigate changes in longitudinal parameters after stereotactic radiotherapy for lung cancer and to identify possible pretreatment factors related to radiation-induced lung toxicity and the decline in pulmonary function after radiotherapy.

Materials and methods

Protocol-specified examinations, including 4-D CT, laboratory tests, pulmonary function tests (PFTs) and body composition measurements, were performed before SRT and at 1 month, 4 months and 12 months after stereotactic radiotherapy. Longitudinal differences were tested by using repeated-measures analysis of variance. Correlations were examined by using the Pearson product-moment correlation coefficient (r).

Results

Sixteen patients were analyzed in this study. During a median follow-up period of 26.6 months, grade 1 and 2 lung toxicity occurred in 11 patients and 1 patient, respectively. The mean Hounsfield units (HU) and standard deviation (SD) of the whole lung, as well as sialylated carbohydrate antigen KL-6 (KL-6) and surfactant protein-D (SP-D), peaked at 4 months after radiotherapy (p = 0.11, p<0.01, p = 0.04 and p<0.01, respectively). At 4 months, lung V_{20 Gy} (%) and V_{40 Gy} (%) were correlated with changes in SP-D, whereas changes in the mean HU of the lung were related to body mass index and lean body mass index (r = 0.54, p = 0.02; r = 0.57, p = 0.01; r = 0.69, p<0.01; and r = 0.69, p<0.01, respectively). The parameters of PFTs gradually declined over time. When regarding the change in PFTs from pretreatment to 12 months, lung V_{5 Gy} (cc) showed significant correlations with diffusion capacity for carbon monoxide (DLCO), DLCO/alveolar volume and the relative change in DLCO (r = -0.72, p<0.01; r = -0.73, p<0.01; and r = -0.63, p = 0.01, respectively).

Conclusions

The results indicated that some parameters peaked at 4 months, but PFTs were the lowest at 12 months. Significant correlations between lung V_{5 Gy} (cc) and changes in DLCO and DLCO/alveolar volume were observed.

Introduction

Systematic therapy for non-small cell lung cancer has shown considerable progress for the past two decades, which is due to the development of new drugs, especially small molecule tyrosine kinase inhibitors and immune checkpoint inhibitors [1–3]. Therefore, systemic therapy is determined by mutations in driver oncogenes and immune checkpoint protein expression, in addition to individual factors. These targeted therapies and immunotherapies are used not only for metastatic lung cancer but also for operable locally advanced lung cancer as neoadjuvant or adjuvant therapies [4, 5]. When regarding early-stage non-small cell lung cancer, surgical resection is a standard treatment [6]. The judgment of operability is performed by thoracic surgeons based on various individual factors, including performance status, frailty, comorbidities and the results of pulmonary function tests (PFTs). In surgical trials, expected postoperative PFTs are used, such as the requirement for an expected postoperative forced expiratory volume in 1 second (FEV1) of 0.8 L or more [7]. Furthermore, stereotactic radiotherapy (SRT), which is often selected as a definitive treatment method for early-stage lung cancer in medically inoperable patients, requires only pretreatment PFTs with sometimes no cutoff value for eligibility being required, which is due to the fact that the post-SRT expectation of PFTs such as surgery is difficult in SRT trials [8]. Additionally, due to the fact that SRT has been performed in such circumstances, the toxicity rates have varied. In the RTOG 0236 trial, 16.3% of the patients had protocol-specified grade 3 or higher toxicity, and an additional 10.9% of the patients had non-protocol-specified high-grade toxicity [9]. The rates of grade 3 or higher adverse events were approximately 10% for patients with noncentral lung cancer (regardless of operability in the subsequent prospective trials); moreover, these rates were 11.9% in the RTOG 0915 trial, 9.6% in pooled analyses for randomized trials for operable patients, 10.6% in the CHISEL trial and 10.3% in the JCOG0403 trial [8–11]. The recent revised STARS trial demonstrated a lower rate of grade 3 or higher adverse events, with only 1 of 80 patients (1.2%) developing grade 3 adverse events [12]. To investigate these differences and the post-SRT effects on PFTs and other factors, more comprehensive analyses are needed. Symptoms such as dyspnea and cough are the complex result of various factors, including a decline in pulmonary function, lung consolidation and systemic or local inflammation. To determine these effects, the primary endpoint of this study was to identify possible pretreatment factors that are related to radiation-induced lung toxicity (RILT) and RILT-related markers at the radiation pneumonitis (RP) phase; in addition, the secondary endpoints were to assess changes in longitudinal parameters after SRT and to identify possible pretreatment factors that have relationships with PFTs at the 12-month follow-up after SRT.

Materials and methods

Inclusion criteria and consent from patients

This prospective study was performed to assess RILT and its related factors after SRT for lung cancer. The main inclusion criteria were as follows: pathologically or clinically diagnosed early-stage primary lung cancer with a tumor diameter of 50 mm or less; agreement between the pulmonologists and radiation oncologists to perform SRT for lung cancer; and no past history of radiotherapy for the thorax. After approval of the Ethical Committee of Tohoku University Hospital (reference number: 2018-2-117), this study was initiated in August 2018 [13]. Written informed consent was obtained from all of the participants. All of the methods were performed in accordance with the Declaration of Helsinki.

Baseline assessment and follow-up of patients

Protocol-specified examinations, including 4-D CT, laboratory tests, PFTs and body composition measurements, were performed before SRT and at 1 month, 4 months and 12 months after SRT. Follow-up examinations were performed by radiation oncologists and pulmonologists, and adverse events were graded by using the National Cancer Institute Common Terminology Criteria for Adverse Events version 4.0. Four-dimensional CT was performed by using a SOMATOM Definition AS+ CT scanner (Siemens Medical, Iselin, NJ), and 4-DCT images were divided into 10 phases (CT00: end inspiration; CT50: end expiration). Body composition measurements were performed by using the body composition analyzer DC-217A (Tanita, Tokyo, Japan). Lean body mass (LBM) is calculated as the difference between body weight and body fat weight. In addition, body mass index (BMI) is calculated as body weight divided by the square of height in meters, and lean body mass index (LBMI) is calculated as LBM divided by the square of height in meters.

SRT procedure

Each patient was immobilized in the supine position with a VacQfix Cushion (Qfix, Avondale, PA), after which a planning CT scan at slice intervals of 2 mm was performed. Tumor and organ delineation and radiotherapy planning were performed by using Eclipse (Varian Medical Systems, Palo Alto, CA). Gross tumor volume (GTV) was determined on the basis of the visible extent of the tumor on planning CT, and internal GTV was determined by using 4-D CT images. A planning target volume (PTV) was created by adding 5 mm around the internal GTV for inter- and intrafractional uncertainty. Twenty-eight or 30 Gy in 1 fraction, 48 Gy in 4 fractions or 60 Gy in 8 fractions was prescribed to cover 95% of the PTV by using 6 MV-FFF beams. The planning algorithm included an Acuros XB; in addition, the radiotherapy technique was volumetric modulated arc therapy, and the treatment machine was a TrueBeam STx (Varian Medical Systems, Palo Alto, CA).

Dosimetric and acquisition data measurements

Lung dose-volume data were used as V_{n Gy} (%) or V_{n Gy} (cc), which were defined as the percentage or volume of the total lung volume (autosegmented lung minus GTV) exceeding n Gy, respectively. The equivalent dose in 2 Gy fractions (EQD2) was calculated by using a linear-quadratic model with the following formula: N×d×(d+α/β)/(2+α/β), wherein N is the number of fractions, d is the dose per fraction and α/β was applied to 3 Gy for normal lungs based on previous findings [14]. Subsequently, EQD2 was used to analyze V_{n Gy} of different SRT schedules, and EQD2 of 5 Gy, 20 Gy and 40 Gy were applied to lung V_{n Gy} for dose-volume analyses, due to the fact that we hypothesized that different lung V_{n Gy} would affect different parameters that may have different radiosensitivities. Mean Hounsfield units (HU) of the lung were measured via CT50 images by using the autosegmented lung of Eclipse in addition to GTV and consolidation, due to the fact that the distinction between GTV and consolidation after SRT was sometimes difficult to observe. Low attenuation area (LAA) was calculated as the magnitude of the overlap between the lung and -860 or -960 HU or lower area that was based on previous findings [15].

Outcomes for assessment and statistical analysis

Time series analyses were performed by using the series of data from each of the assessment points. Any overall longitudinal differences were tested by using repeated-measures analysis of variance. To identify possible predictive factors for RP, changes in the following parameters at the 4-month follow-up after SRT were used: mean HU of the whole lung by using CT50 images, C-reactive protein (CRP), sialylated carbohydrate antigen KL-6 (KL-6) and surfactant protein-D (SP-D). To identify possible predictive factors for the decline in pulmonary function, changes in the following parameters at 12 months after SRT were used: vital capacity (VC, L), forced vital capacity (FVC, L), FEV1 (L), FEV1/FVC (%), FVC1% of predicted, diffusion capacity for carbon monoxide (DLCO, mmol/min/kPa) and DLCO/alveolar volume (VA, L). Delta (Δ) represented the change in each parameter that was calculated as “the observed value minus the baseline value of the parameter”, and relative Δ was calculated as “Δ divided by the baseline value”. Moreover, DLCO was adjusted for hemoglobin change from baseline. Correlations between predictive pretreatment factors and parameter changes were tested. Additionally, correlations were examined by using the Pearson product-moment correlation coefficient (r). Although multiple testing was performed for some factors (such as lung dose-volume data), a p value less than 0.05 was defined as being statistically significant when considering the exploratory study design. The onset of RILT was calculated from the first day of SRT to the day that an event was confirmed. Cumulative incidences were calculated by using the Kaplan‒Meier method, and death was regarded as a competing risk. EZR version 1.54 (Saitama Medical Center, Jichi Medical University, Saitama, Japan), which is a modified version of R commander (R Foundation for Statistical Computing, Vienna, Austria), was used for the analyses [16].

Results

Seventeen patients were enrolled between August 2018 and February 2021, and one patient later withdrew from the data analysis of the study. As a result, 16 patients were analyzed, and the date of data cutoff was March 31, 2022. Patient characteristics and baseline data are shown in Tables 1 and 2. The pretreatment mean and median (range) values of CRP (mg/dL) were 0.20 and 0.06 (0.01–0.87), respectively. All of the patients completed 12-month examinations at the date of data cutoff, but 1 patient could not undergo spirometry after surgery for oral cancer. Additionally, four sets of 4-D CT images could not be divided into 10 phases. When regarding SRT doses, 3 patients, 2 patients, 6 patients and 5 patients received 28 Gy in 1 fraction, 30 Gy in 1 fraction, 48 Gy in 4 fractions and 60 Gy in 8 fractions, respectively. Details of the PTV and dosimetric parameters are shown in Table 1.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

During a median follow-up period of 26.6 months (range: 10.8–37.6 months), one patient died from lung cancer at 16 months after SRT. Thoracic radiotherapy, thoracic surgery or systemic chemotherapy was not performed within 12 months after SRT. Grade 1 or higher RP occurred in 12 patients, including 1 patient with grade 2 RP. Moreover, the one-year cumulative incidences of grade 1 or higher and grade 2 RP were 75.0% (95% confidence interval [CI]: 52.8–92.2%) and 5.2% (95% CI: 0.9–36.8%), respectively. Due to the fact that grade 2 or higher RP occurred in only 1 patient, a predictive factor analysis for RP was not performed.

Longitudinal analyses are shown in Table 2. The mean HU and standard deviation (SD) of the whole lung, KL-6 and SP-D peaked at 4 months after SRT, and these changes were significant, excluding the mean HU of the lung (p = 0.11, p<0.01, p = 0.04 and p<0.01, respectively). All parameters of PFTs excluding FVC gradually declined over time; therefore, they showed the lowest values at the 12-month follow-up after SRT, but only changes in FEV1 and DLCO were significant (p = 0.02 for both variables). Longitudinal changes in Δmean HU of the whole lung, ΔCRP, ΔKL-6 and ΔSP-D from baseline to each follow-up period according to RP are shown in Fig 1.

[Figure omitted. See PDF.]

Although Δmean HU of the whole lung, ΔKL-6 and ΔSP-D peaked at 4 months, ΔKL-6 in the patients with grade 2 radiation pneumonitis was decreased at 4 months.

Pearson’s r for predictive factors of Δmean HU of the whole lung by using CT50 images, ΔCRP, ΔKL-6 and ΔSP-D at the 4-month follow-up after SRT are shown in Table 3. Among these factors, only ΔSP-D had a significant correlation with dosimetric factors: lung V_{20 Gy} (%) and V_{40 Gy} (%) (r = 0.54, 95% CI: 0.06‐0.81, p = 0.02 and r = 0.57, 95% CI: 0.11‐0.83, p = 0.01, respectively). When regarding Δmean HU of the whole lung using CT50 images, ΔCRP and ΔKL-6 showed significant correlations with pretreatment BMI, LBMI or both (Table 3); specifically, BMI was related to Δmean HU of the lung and ΔKL-6 (r = 0.69, 95% CI: 0.26‐0.89, p<0.01 and r = 0.54, 95% CI: 0.07‐0.82, p = 0.02, respectively), and LBMI was related to Δmean HU of the lung and ΔCRP (r = 0.69, 95% CI: 0.25‐0.89, p<0.01 and r = 0.55, 95% CI: 0.07‐0.82, p = 0.02, respectively).

[Figure omitted. See PDF.]

Pearson’s r of absolute and relative changes in parameters of pulmonary function tests from baseline to 12 months after SRT are shown in Tables 4 and 5. Among dosimetric factors, lung V_{5 Gy} (cc) showed the strongest correlations with ΔDLCO and ΔDLCO/VA (r = -0.72, 95% CI: -0.90‐-0.34, p<0.01 and r = -0.73, 95% CI: -0.90‐-0.35, p<0.01, respectively) and the second strongest correlation with relative ΔDLCO (r = -0.63, 95% CI: -0.86‐-0.18, p = 0.01). When regarding ΔVC, ΔFVC, ΔFEV1, ΔFEV1/FVC and ΔFEV1% of predicted, there were no correlations of dosimetric and other pretreatment factors with them, except for correlations between tumor diameter and ΔVC (r = -0.55, 95% CI: -0.83‐-0.06, p = 0.03) and between LBM and ΔFEV1% of predicted (r = 0.54, 95% CI: 0.04‐0.82, p = 0.03). The correlation of ΔFEV1% predicted indicated that patients with a higher pretreatment LBM had a smaller decline in FEV1% predicted at 12 months after SRT. Scatter plots of the highest correlations in Table 4 of ΔVC, ΔFEV1% of predicted, ΔDLCO and ΔDLCO/VA are shown in Fig 2.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Discussion

In this study, possible RILT-related parameters were prospectively analyzed. Analyses were performed for longitudinal changes in parameters, correlations between possible RILT-related parameters at 4 months after SRT (RP phase) and pretreatment factors, as well as correlations between PFTs at 12 months after SRT and pretreatment factors. RILT consists of these factors and various other factors. Therefore, analyses for each parameter would be beneficial, especially in the current era of radiotherapy combined with immune checkpoint inhibitors to reduce grade 2 or more RILT [17]. In the analyses of longitudinal changes after SRT, the mean HU and SD of the whole lung, KL-6 and SP-D peaked at the 4-month follow-up after SRT, which would reflect the emergence of RP. Prior to the changes in the parameters at the 4-month follow-up after SRT, a decline in lymphocytes and an increase in SP-D were observed at the 1-month follow-up. Although these early changes are interesting, it is not clear as to whether early changes in the parameters reflect the severity of RILT. Although the mean HU and SD of the whole lung were not changed at the 1-month follow-up in this analysis, an early graphical change after SRT was reported to correlate with severe RP [18]. Another study showed that RP on CT occurred after a median period of 2.9 months after SRT for the symptomatic RP group in comparison with 5.1 months for the asymptomatic RP group [19]. An early emergence of an RP shadow would be needed for a careful follow-up. When regarding CRP, no statistical change was observed, although there was no abnormal distribution of CPR; specifically, the range of pretreatment CRP was 0.01–0.87. It was reported that pretreatment CRP level and maximum CRP after radiotherapy may predict or reflect symptomatic RILT [20, 21]. Actually, the CRP of grade 2 RILT case peaked at 4 months in Fig 1; therefore, the result of this study did not contradict the importance of CRP.

When regarding longitudinal changes in parameters at the 4-month follow-up, only SP-D had relationships with dosimetric parameters, and the relationship between SP-D and lung V_{40 Gy} (%) was the strongest (r = 0.57). Although SP-D peaked at 4 months after SRT in patients with grade 1 and grade 2 RILT, KL-6 in the patient with grade 2 RILT was decreased at the 4-month follow-up, and the results of the studies are therefore different (Fig 1). This discrepancy in the results may be because KL-6 and SP-D have been reported to only exhibit an elevation in 70% of patients with grade 4–5 RILT [22]. At the least, the change in the mean HU of the whole lung peaked at 4 months after SRT. Although the change in mean HU had no correlation with dosimetric factors, this parameter had strong correlations with BMI and LBMI (r = 0.69 and r = 0.69, respectively). The interpretation of these correlations was difficult because BMI was affected by many factors. For example, smoking was associated with lower BMI, but smoking cessation was associated with higher BMI [23]. When patients develop lung emphysema or chronic obstructive pulmonary disease, cachexia can develop [24]. Patients with LAA, which is often caused by chronic obstructive pulmonary disease, tend to show less HU change [15]. There were no significant lung dose-volume correlations with changes in lung HU on CT images, whereas various relationships with quantified CT images have been reported. De Ruysscher et al. reported that they quantified RILT by using the average HU change in each of the irradiated lung dose bins, and they demonstrated that the mean change in HU per Gy was 1.7 ± 2.0 [25]. Individual radiosensitivity differences were also suggested, with HU/Gy ranging from 0 to 10. This difference may also reflect different RILT shadow patterns after radiotherapy [26]. In another study, the normal tissue complication probability of the lung after SRT was investigated by using HU changes in each lung pixel on CT [27]. They reported that D₅₀, which was defined as the dose with half of all subjects presenting with a symptom, was approximately 35 Gy using 4 NTCP models. Due to the fact that SRT is highly focused radiation, the lung dose-volume effect on the mean HU of the whole lung may be relatively weak, and analyses of local HU changes in the lung will show stronger relationships with lung dose-volume parameters.

When regarding changes in PFT parameters, all of the PFT values gradually declined with time, as was previously reported [28]. Although only FEV1 and DLCO showed significant changes in this study, significant reductions in FVC, FVC % of predicted and FEV1% of predicted after SRT have also been reported [19, 28]. In this study, the median and mean ± SD of relative ΔFEV1 and those of relative ΔDLCO from pretreatment to 12-month follow-up after SRT were -4.2% and -6.8% ± 10.5 and -10.2% and -8.6% ± 11.1, respectively; additionally, those findings were almost the same as those from a previous study [28, 29]. Absolute values are also important because baseline values in SRT patients may be lower than those in surgical candidates, and the median and mean ± SD of ΔFEV1 and those of ΔDLCO were -0.06 L and -0.12 L ± 0.16 and -1.33 and -1.02 ± 1.47, respectively. In surgical patients, it was reported that the mean ± SD values of FEV1 and the transfer factor of the lungs for carbon monoxide (mL/min/mmHg) declined from 2.38 ± 0.79 and 22.3 ± 6.7 to 2.17 ± 0.73 and 21.4 ± 5.9, respectively, at 6 months after lobectomy; additionally, these parameters declined from 2.50 ± 0.47 and 22.8 ± 6.7 to 1.65 ± 0.29 and 16.4 ± 3.7, respectively, at 6 months after pneumonectomy [30]. When considering these results, the magnitude of the reductions in PFT parameters after SRT was relatively small, but caution is needed when treating patients with extremely low pulmonary function.

Correlations in the parameters of PFTs at the 12-month follow-up are shown in Tables 4 and 5. Lung dose-volume parameters, especially lung V_{5 Gy} (cc), showed strong correlations with ΔDLCO and ΔDLCO/VA. These results suggested that DLCO may be strongly affected by a low dose distribution in contrast to other PFTs. Therefore, DLCO was thought to have the highest radiosensitivity; additionally, to the best of our knowledge, significant lung dose-volume relationships after SRT with ΔDLCO and ΔDLCO/VA have not been previously reported. Although there were conflicting results for DLCO after SRT, it may be difficult to measure ΔDLCO in a retrospective design because DLCO is a vulnerable parameter and can be affected by chemotherapy agents [31, 32]. Stephans et al. reported correlations of lung ΔFEV1% predicted with lung V5 _Gy (cc) and V10 _Gy (cc), but there was no correlation between ΔDLCO and lung dose-volume parameters in SRT [33]. Some dose-volume relationships with DLCO have been reported in conventional fraction data series. In Hodgkin’s lymphoma patients, a lower mean lung dose in patients treated with bleomycin-based chemotherapy alone or bleomycin-based chemotherapy and mediastinal radiotherapy was significantly related to a lower decline in the percentage of predicted DLCO at 1 year [34]. In patients with non-small cell lung cancer treated with 3-dimensional conformal radiotherapy or intensity-modulated radiotherapy or proton beam therapy with or without chemotherapy, several lung dosimetric parameters were significantly correlated with DLCO [35]. In this study, LBM was correlated with ΔFEV1% of predicted significance and ΔFEV1 with marginal significance, but lung dose-volume had no significant correlation. Patients with a high LBM may compensate for the decline in FEV1 and FEV1% predicted by using respiratory muscles, when assuming that patients with a high LBM have increased respiratory muscles. In another study, it was shown that there were significant correlations between lung dose-volume parameters and the relative ΔFEV1 in subgroups of patients divided by the severity of RP [19]. Dose-volume relationships appeared as early as 1 month after SRT in the symptomatic RP subgroup but did not appear until 6 months in the asymptomatic RP subgroup; in addition, there was no significant correlation in the no RP subgroup.

In recent radiotherapy techniques, some lung areas can be intentionally avoided, and there is a new attempt at CT ventilation functional image-guided radiotherapy [36]. A CT ventilation functional image was obtained by 4-D CT images, and this ventilation image had a correlation with FEV1% of predicted and FEV1/FVC [37]. Due to the fact that we can see where it is functional on CT images, functional lung areas can be avoided in radiotherapy planning, and it may be beneficial to avoid functional lung areas leading to a lower incidence of severe RILT. When regarding DLCO, some CT findings have been reported to be associated with DLCO [38]. Although we do not know if there are higher and lower functional areas of DLCO in the lung, further investigations, including radiomics and dosiomics, may aid in visualizing functional DLCO areas in CT images [39, 40]. More precisely, by considering overall pulmonary function, radiotherapy planning can be offered to reduce symptomatic or severe RILT.

There were several limitations in this study. For example, the number of patients and the number of events were relatively small. Only one patient developed grade 2 or more RILT. Due to the fact that this study was an exploratory study, many comparisons were performed; therefore, there may be some false-positives, and further studies are needed to confirm these findings. In addition, the SRT treatment schedule varied; therefore, there was a limitation of dose correction to compare the different radiation schedules. Specifically, there was a limitation of the linear-quadratic model estimation.

In conclusion, different lung dose-volume parameters affected RP-related parameters and parameters of PFTs. Some possible RILT markers peaked at 4 months, but PFTs declined over time and were the lowest at 12 months. The analyses of PFT parameter changes from pretreatment to 12 months after SRT found that DLCO and DLCO/VA were significantly correlated with lung V_{5 Gy} (cc), which indicated that DLCO and DLCO/VA were possibly affected by a relatively lower radiation dose distribution.

Supporting information

S1 File.

https://doi.org/10.1371/journal.pone.0278707.s001

(XLSX)

Acknowledgments

We are grateful to the radiation oncologists, pulmonologists, medical physicists and radiation technologists at Tohoku University Hospital who contributed to the recruitment, treatment and follow-up of the patients. When creating the protocol, advice was provided to us by Dr. Soichiro Toda from the Clinical Research, Innovation and Education Center, Tohoku University Hospital.

Citation: Yamamoto T, Katsuta Y, Sato K, Tsukita Y, Umezawa R, Takahashi N, et al. (2022) Longitudinal analyses and predictive factors of radiation-induced lung toxicity-related parameters after stereotactic radiotherapy for lung cancer. PLoS ONE 17(12): e0278707. https://doi.org/10.1371/journal.pone.0278707

About the Authors:

Takaya Yamamoto

Roles: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Validation, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

https://orcid.org/0000-0003-3562-1037

Yoshiyuki Katsuta

Roles: Data curation, Investigation, Supervision, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

Kiyokazu Sato

Roles: Investigation, Resources, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

Yoko Tsukita

Roles: Investigation, Resources, Writing – review & editing

Affiliation: Department of Respiratory Medicine, Tohoku University Graduate School of Medicine, Sendai, Japan

Rei Umezawa

Roles: Investigation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

Noriyoshi Takahashi

Roles: Investigation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

https://orcid.org/0000-0002-7653-0046

Yu Suzuki

Roles: Investigation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

Kazuya Takeda

Roles: Investigation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

https://orcid.org/0000-0002-7354-4063

Keita Kishida

Roles: Investigation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

So Omata

Roles: Investigation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

Eisaku Miyauchi

Roles: Investigation, Resources, Writing – review & editing

Affiliation: Department of Respiratory Medicine, Tohoku University Graduate School of Medicine, Sendai, Japan

Ryota Saito

Roles: Investigation, Resources, Writing – review & editing

Affiliation: Department of Respiratory Medicine, Tohoku University Graduate School of Medicine, Sendai, Japan

Noriyuki Kadoya

Roles: Investigation, Resources, Validation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

Keiichi Jingu

Roles: Investigation, Project administration, Supervision, Validation, Writing – review & editing

Affiliation: Department of Radiation Oncology, Tohoku University Graduate School of Medicine, Sendai, Japan

References

1. Rizzo A, Cusmai A, Giovannelli F, Acquafredda S, Rinaldi L, Misino A, et al. Impact of Proton Pump Inhibitors and Histamine-2-Receptor Antagonists on Non-Small Cell Lung Cancer Immunotherapy: A Systematic Review and Meta-Analysis. Cancers (Basel). 2022;14:1404. pmid:35326555

2. Lamberti G, Andrini E, Sisi M, Rizzo A, Parisi C, Di Federico A, et al. Beyond EGFR, ALK and ROS1: Current evidence and future perspectives on newly targetable oncogenic drivers in lung adenocarcinoma. Crit Rev Oncol Hematol. 2020;156:103119. pmid:33053439

3. Rizzo A. Chemoimmunotherapy versus immune checkpoint inhibitors monotherapy as first-line treatment for advanced non-small cell lung cancer. Thorac Cancer. 2022;13:656. pmid:35029065

4. Wu YL, John T, Grohe C, Majem M, Goldman JW, Kim SW, et al. Postoperative Chemotherapy Use and Outcomes From ADAURA: Osimertinib as Adjuvant Therapy for Resected EGFR-Mutated NSCLC. J Thorac Oncol. 2022;17:423–33. pmid:34740861

5. Felip E, Altorki N, Zhou C, Csőszi T, Vynnychenko I, Goloborodko O, et al. Adjuvant atezolizumab after adjuvant chemotherapy in resected stage IB-IIIA non-small-cell lung cancer (IMpower010): a randomised, multicentre, open-label, phase 3 trial. Lancet. 2021;398:1344–57. pmid:34555333

6. NCCN Guidelines Version 3.2022- Non-Small Cell Lung Cancer, accessed April 30 2022.

7. Saji H, Okada M, Tsuboi M, Nakajima R, Suzuki K, Aokage K, et al. Segmentectomy versus lobectomy in small-sized peripheral non-small-cell lung cancer (JCOG0802/WJOG4607L): a multicentre, open-label, phase 3, randomised, controlled, non-inferiority trial. Lancet. 2022;399:1607–17. pmid:35461558.

8. Ball D, Mai GT, Vinod S, Babington S, Ruben J, Kron T, et al. Stereotactic ablative radiotherapy versus standard radiotherapy in stage 1 non-small-cell lung cancer (TROG 09.02 CHISEL): a phase 3, open-label, randomised controlled trial. Lancet Oncol. 2019;20:494–503. pmid:30770291.

9. Timmerman R, Paulus R, Galvin J, Michalski J, Straube W, Bradley J, et al. Stereotactic body radiation therapy for inoperable early stage lung cancer. JAMA. 2010;303:1070–6. pmid:20233825.

10. Videtic GM, Hu C, Singh AK, Chang JY, Parker W, Olivier KR, et al. A Randomized Phase 2 Study Comparing 2 Stereotactic Body Radiation Therapy Schedules for Medically Inoperable Patients With Stage I Peripheral Non-Small Cell Lung Cancer: NRG Oncology RTOG 0915 (NCCTG N0927). Int J Radiat Oncol Biol Phys. 2015;93:757–64. pmid:26530743.

11. Nagata Y, Hiraoka M, Shibata T, Onishi H, Kokubo M, Karasawa K, et al. Prospective Trial of Stereotactic Body Radiation Therapy for Both Operable and Inoperable T1N0M0 Non-Small Cell Lung Cancer: Japan Clinical Oncology Group Study JCOG0403. Int J Radiat Oncol Biol Phys. 2015;93:989–96. pmid:26581137.

12. Chang JY, Mehran RJ, Feng L, Verma V, Liao Z, Welsh JW, et al. Stereotactic ablative radiotherapy for operable stage I non-small-cell lung cancer (revised STARS): long-term results of a single-arm, prospective trial with prespecified comparison to surgery. Lancet Oncol. 2021;22:1448–57. pmid:34529930.

13. https://upload.umin.ac.jp/cgi-open-bin/icdr_e/ctr_view.cgi?recptno=R000036774, accessed April 30 2022.

14. Fowler JF, Tomé WA, Fenwick JD, Mehta MP. A challenge to traditional radiation oncology. Int J Radiat Oncol Biol Phys. 2004;60:1241–56. pmid:15519797.

15. Yamamoto T, Kadoya N, Sato Y, Matsushita H, Umezawa R, Kubozono M, et al. Prognostic Value of Radiation Pneumonitis After Stereotactic Body Radiotherapy: Effect of Pulmonary Emphysema Quantitated Using CT Images. Clin Lung Cancer. 2018;19:e85–e90. pmid:28655592.

16. Kanda Y. Investigation of the freely available easy-to-use software ’EZR’ for medical statistics. Bone Marrow Transplant. 2013;48:452–8. pmid:23208313.

17. Yamamoto T, Tsukita Y, Katagiri Y, Matsushita H, Umezawa R, Ishikawa Y, et al. Durvalumab after chemoradiotherapy for locally advanced non-small cell lung cancer prolonged distant metastasis-free survival, progression-free survival and overall survival in clinical practice. BMC Cancer. 2022;22:364. pmid:35379201.

18. Takeda A, Ohashi T, Kunieda E, Enomoto T, Sanuki N, Takeda T, et al. Early graphical appearance of radiation pneumonitis correlates with the severity of radiation pneumonitis after stereotactic body radiotherapy (SBRT) in patients with lung tumors. Int J Radiat Oncol Biol Phys. 2010;77:685–90. pmid:20510193.

19. Berg J, Ramberg C, Haugstvedt JOS, Bengtson MB, Gabrielsen AM, Brustugun OT, et al. Lung Function After Stereotactic Body Radiation Therapy for Early-Stage Non-Small Cell Lung Cancer, Changes and Predictive Markers. Front Oncol. 2021;11:674731. pmid:34109123.

20. Inoue A, Kunitoh H, Sekine I, Sumi M, Tokuuye K, Saijo N. Radiation pneumonitis in lung cancer patients: a retrospective study of risk factors and the long-term prognosis. Int J Radiat Oncol Biol Phys. 2001;49:649–55. pmid:11172945.

21. Liu F, Qiu B, Xi Y, Luo Y, Luo Q, Wu Y, et al. Efficacy of Thymosin α1 in Management of Radiation Pneumonitis in Patients With Locally Advanced Non-Small Cell Lung Cancer Treated With Concurrent Chemoradiotherapy: A Phase 2 Clinical Trial (GASTO-1043). Int J Radiat Oncol Biol Phys. 2022;114:433–43. pmid:35870709.

22. Yamashita H, Kobayashi-Shibata S, Terahara A, Okuma K, Haga A, Wakui R, et al. Prescreening based on the presence of CT-scan abnormalities and biomarkers (KL-6 and SP-D) may reduce severe radiation pneumonitis after stereotactic radiotherapy. Radiat Oncol. 2010;5:32. pmid:20459699.

23. Piirtola M, Jelenkovic A, Latvala A, Sund R, Honda C, Inui F, et al. Association of current and former smoking with body mass index: A study of smoking discordant twin pairs from 21 twin cohorts. PLoS One. 2018;13:e0200140. pmid:30001359

24. Wagner PD. Possible mechanisms underlying the development of cachexia in COPD. Eur Respir J. 2008;31:492–501. pmid:18310396.

25. De Ruysscher D, Sharifi H, Defraene G, Kerns SL, Christiaens M, De Ruyck K, et al. Quantification of radiation-induced lung damage with CT scans: the possible benefit for radiogenomics. Acta Oncol. 2013;52:1405–10. pmid:23957564.

26. Yamamoto T, Kadoya N, Morishita Y, Sato Y, Matsushita H, Umezawa R, et al. Assessment and agreement of the CT appearance pattern and its severity grading of radiation-induced lung injury after stereotactic body radiotherapy for lung cancer. PLoS One. 2018;13:e0204734. pmid:30286105.

27. Begosh-Mayne D, Kumar SS, Toffel S, Okunieff P, O’Dell W. The dose-response characteristics of four NTCP models: using a novel CT-based radiomic method to quantify radiation-induced lung density changes. Sci Rep. 2020;10:10559. pmid:32601297.

28. Stone B, Mangona VS, Johnson MD, Ye H, Grills IS. Changes in Pulmonary Function Following Image-Guided Stereotactic Lung Radiotherapy: Neither Lower Baseline Nor Post-SBRT Pulmonary Function Are Associated with Worse Overall Survival. J Thorac Oncol. 2015;10:1762–9. pmid:26334751.

29. Niezink AGH, de Jong RA, Muijs CT, Langendijk JA, Widder J. Pulmonary Function Changes After Radiotherapy for Lung or Esophageal Cancer: A Systematic Review Focusing on Dose-Volume Parameters. Oncologist. 2017;22:1257–64. pmid:28550029.

30. Bolliger CT, Jordan P, Solèr M, Stulz P, Tamm M, Wyser C, et al. Pulmonary function and exercise capacity after lung resection. Eur Respir J. 1996;9:415–21. pmid:8729998.

31. Guckenberger M, Klement RJ, Kestin LL, Hope AJ, Belderbos J, Werner-Wasik M, et al. Lack of a dose-effect relationship for pulmonary function changes after stereotactic body radiation therapy for early-stage non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2013;85:1074–81. pmid:23154077.

32. Landman Y, Stemmer SM, Sulkes A, Neiman V, Granot T, Hendler D, et al. Prospective Long-Term Follow-Up of Pulmonary Diffusion Capacity Reduction Caused by Dose-Dense Chemotherapy in Patients with Breast Cancer. J Oncol. 2019;2019:2584859. pmid:31772578.

33. Stephans KL, Djemil T, Reddy CA, Gajdos SM, Kolar M, Machuzak M, et al. Comprehensive analysis of pulmonary function Test (PFT) changes after stereotactic body radiotherapy (SBRT) for stage I lung cancer in medically inoperable patients. J Thorac Oncol. 2009;4:838–44. pmid:19487961.

34. Ng AK, Li S, Neuberg D, Chi R, Fisher DC, Silver B, et al. A prospective study of pulmonary function in Hodgkin’s lymphoma patients. Ann Oncol. 2008;19:1754–8. pmid:18467315.

35. Lopez Guerra JL, Gomez DR, Zhuang Y, Levy LB, Eapen G, Liu H, et al. Changes in pulmonary function after three-dimensional conformal radiotherapy, intensity-modulated radiotherapy, or proton beam therapy for non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2012;83:e537–43. pmid:22420964.

36. Yamamoto T, Kabus S, Bal M, Keall P, Benedict S, Daly M. The first patient treatment of computed tomography ventilation functional image-guided radiotherapy for lung cancer. Radiother Oncol. 2016;118:227–31. pmid:26687903.

37. Yamamoto T, Kabus S, Lorenz C, Mittra E, Hong JC, Chung M, et al. Pulmonary ventilation imaging based on 4-dimensional computed tomography: comparison with pulmonary function tests and SPECT ventilation images. Int J Radiat Oncol Biol Phys. 2014;90:414–22. pmid:25104070.

38. Grydeland TB, Thorsen E, Dirksen A, Jensen R, Coxson HO, Pillai SG, et al. Quantitative CT measures of emphysema and airway wall thickness are related to D(L)CO. Respir Med. 2011;105:343–51. pmid:21074394.

39. Jiang W, Song Y, Sun Z, Qiu J, Shi L. Dosimetric Factors and Radiomics Features Within Different Regions of Interest in Planning CT Images for Improving the Prediction of Radiation Pneumonitis. Int J Radiat Oncol Biol Phys. 2021;110:1161–70. pmid:33548340.

40. Adachi T, Nakamura M, Shintani T, Mitsuyoshi T, Kakino R, Ogata T, et al. Multi-institutional dose-segmented dosiomic analysis for predicting radiation pneumonitis after lung stereotactic body radiation therapy. Med Phys. 2021;48:1781–91. pmid:33576510.