Stereotactic ablative radiotherapy for ultra-central lung tumors: prioritize target coverage or organs at risk?

Abstract

Background

Lung stereotactic ablative radiotherapy (SABR) is associated with low morbidity, however there is an increased risk of treatment-related toxicity in tumors directly abutting or invading the proximal bronchial tree, termed ‘ultra-central’ tumors. As there is no consensus regarding the optimal radiotherapy treatment regimen for these tumors, we performed a modeling study to evaluate the trade-offs between predicted toxicity and local control for commonly used high-precision dose-fractionation regimens.

Methods

Ten patients with ultra-central lung tumors were identified from our institutional database. New plans were generated for 3 different hypofractionated schemes: 50 Gy in 5 fractions, 60 Gy in 8 fractions and 60 Gy in 15 fractions. For each regimen, one plan was created that prioritized planning target volume (PTV) coverage, potentially at the expense of organ at risk (OAR) tolerance, and a second that compromised PTV coverage to respect OAR dose constraints. Published radiobiological models were employed to evaluate competing treatment plans based on estimates for local control and the likelihood for toxicity to OAR.

Results

The risk of esophageal or pulmonary toxicity was low (< 5%) in all scenarios. When PTV coverage was prioritized, tumor control probabilities were 92.9% for 50 Gy in 5 fractions, 92.4% for 60 Gy in 8 fractions, and 52.0% for 60 Gy in 15 fractions; however the estimated risk of grade ≥ 4 toxicity to the proximal bronchial tree was 68%, 44% and 2% respectively. When dose to OAR was prioritized, the risk of major pulmonary toxicity was reduced to < 1% in all schemes, but this compromise reduced tumor control probability to 60.3% for 50 Gy in 5 fractions, 65.7% for 60 Gy in 8 fractions and 47.8% for 60 Gy in 15 fractions.

Conclusions

The tradeoff between local control and central airway toxicity are considerable in the use of 3 commonly used hypofractionated radiotherapy regimens for ultra-central lung cancer. The results of this planning study predict that the best balance may be achieved with 60 Gy in 8 fractions compromising PTV coverage as required to maintain acceptable doses to OAR. A prospective phase I trial (SUNSET) is planned to further evaluate this challenging clinical scenario.

Background

Stereotactic ablative radiotherapy (SABR) is an established treatment for medically inoperable patients with early stage non-small cell lung cancer (NSCLC). Morbidity is low when SABR is employed for tumors that are peripherally located, and risk-adapted schemes can be employed when the lesion is near critical mediastinal structures at increased risk for treatment-related toxicity. The magnitude of such risk has been the subject of debate and the optimal SABR dose-fractionation regimen for central tumors is an active area of research [1, 2]. In the landmark Radiation Therapy and Oncology Group (RTOG) 0236 trial that demonstrated the efficacy of SABR for peripheral lung tumors; a 2 cm ‘no-fly zone’ around the proximal bronchial tree was coined as an area of trial exclusion [3]. Subsequently, the RTOG 0813 dose escalation trial was initiated to investigate the safety of 5-fraction regimens of SABR for central lung tumors, assessing the increase of total doses of 50 to 60 Gy [4]. Though the highest dose level in this trial was achieved, there have been published reports of fatal toxicity after treatment with 45–50 Gy in 5 fractions [5,6,7]. A more fractionated regimen of 60 Gy in 8 for central lung tumors was investigated by the European Organization for Research and Treatment of Cancer (EORTC) in the phase II trial, LungTech; however, efficacy and toxicity data are not yet reported [8, 9]. Notably, tumors directly abutting or invading the proximal bronchial tree – defined herein as ‘ultra-central’ lung tumors – were not specifically addressed in the RTOG 0813 or LungTech trials and represent what is likely to be a higher-risk subgroup.

Ultra-central lung tumors are a relatively uncommon clinical scenario and reports evaluating SABR outcomes in this context are limited. One study found that 30–40 Gy over 5 fractions resulted in 1-year local control rate of 63% of hilar tumors; even with this moderate dose, a SABR-related death due to bronchial fistula was noted [10]. A similar investigation of 35–40 Gy for lung tumors abutting or invading the mainstem bronchus revealed 70% local control at 1 year, but with no reports of treatment-related death and 1 report of grade 4 atelectasis [11]. Likewise, a retrospective study identified 7 NSCLC patients treated with 50 Gy in 4–5 fractions who all achieved local tumor control at 2 years, without any grade 2 or higher toxicities [12]. In contrast to these results, Haseltine et al. estimated a 22% rate of fatal complications attributable to a 45–50 Gy in 5 fraction SABR regimen when tumors abut the proximal bronchial tree [6]. Finally, a large cohort (n = 47), ranging from stage IA-IIIA and with or without nodal disease, was reported by Tekatli et al. where the actuarial incidence of treatment-related death was 16.3% at 12 months for patients who received 60 Gy in 12 fractions [13]. Considering these findings, concern remains over the balance between local control and toxicity in treating ultra-central lung tumors and there is no consensus regarding the most appropriate SABR or hypofractionated RT scheme.

In the absence of robust evidence for the ideal dose-fractionation regimen, radiobiological modeling can be a helpful tool to evaluate the risks and benefits of different treatment plans. Herein, we employ published radiobiological models to evaluate the likelihood of toxicity and estimate local control for SABR dose-fractionation regimens that are commonly used in the treatment of ultra-central lung tumors.

Methods

Cohort characteristics

Ten patients with ultra-central lung tumors, defined as tumors with a gross tumor volume (GTV) directly abutting or invading the proximal bronchial tree, treated with RT between November 2013 and November 2016 were identified from our institutional database and included in this research ethics board approved planning study. Clinically delivered doses varied substantially. For two patients, ≥95% of the PTV was covered by 60 Gy in 8 fractions; this dose was prescribed in three additional patients, with compromised PTV coverage. For three patients, ≥98% of the PTV was covered by 95% of 60 Gy in 15 fractions. One patient received 54 Gy in 8 fractions with compromised PTV coverage, and one patient received 40 Gy in 15 fractions. The location of an ultra-central tumor relative to the bronchial tree, and other important organs at risk (OAR), is illustrated in Fig. 1 with axial, sagittal, and coronal computed tomography (CT) slices.

Fig. 1

An illustrative example of a contoured ultra-central lung tumor case. The relative locations of tumor and OAR are indicated; IGTV (red), PTV (blue), proximal bronchial tree (yellow), esophagus (green), heart (maroon), and healthy lung (cyan). Here, the PTV volume is 83.1 cm³ and 2.2 cm³ overlap with the proximal bronchial tree

CT simulation

Patients were immobilized in the supine position, with their arms above their head. Fast helical CT Simulation (CT-SIM) scans were acquired, as well as a 4DCT scan to assess respiratory motion. The untagged average, maximum inspiration, and maximum expiration CT datasets were transferred to the Pinnacle treatment planning system (software version 9.10, Philips Healthcare, Fitchburg, WI, USA). Delineation of target volumes and OARs were completed by a thoracic radiation oncologist (JL). The GTV was contoured on the maximum inspiration and maximum expiration 4DCT images and combined to form the internal gross tumor volume (IGTV) to account for tumor motion. A 5 mm expansion on the IGTV was applied to create the planning target volume (PTV) for treatment planning.

Treatment planning

Experienced dosimetrists (AE, BM) generated new radiation treatment plans for each of the 10 cases using three commonly used hypofractionated schemes: 50 Gy in 5 fractions, 60 Gy in 8 fractions, and 60 Gy in 15 fractions. An inhomogeneous dose prescription was used to create SABR plans for 50 Gy in 5 fraction and 60 Gy in 8 fraction regimens; whereas a traditional homogeneous prescription was used for 60 Gy in 15 fractions. Two plans were generated for each of the three dose regimens: one prioritizing tumor coverage (with dose to OARs as low as reasonably achievable) and the other plan compromising PTV coverage in order to respect the dose constraints for OARs. OARs included the proximal bronchial tree, esophagus, great vessels, heart, spinal cord and lungs. All plans were created using volumetric modulated arc therapy (VMAT) with two arcs covering 180–360° using 6 or 10 MV photon beams. The dose constraints for OAR were based on Timmerman’s summary for 5-fraction SBRT and were corrected for dose-fractionation using biologically effective dose calculations (Table 1) [14]. Descriptive statistics were generated to illustrate the dosimetric trade-offs of each scenario.

Table 1 Dose objectives for each fractionation regimen, based on standard linear-quadratic model biologically effective doses calculated with α/β=3

Radiobiological modeling

Treatment plans were evaluated using published radiobiological models for normal tissue complication probability (NTCP) and tumor control probability (TCP). Differential dose volume data was exported from Pinnacle in 1 cGy bins for modeling. To correct for differences in dose fractionation, physical doses were converted to isoeffective doses in 2 Gy per fraction (also termed the equivalent dose in 2 Gy per fraction; EQD₂) using α/β=3 for late responding normal tissue and α/β=10 for tumor (Eq. 1). For each fractional volume, D_i is the total physical dose and d_i is the physical dose per fraction.

$$ {EQD}_2={D}_i\frac{1+\frac{d_i}{\raisebox{1ex}{$\alpha $}\!\left/ \!\raisebox{-1ex}{$\beta $}\right.}}{1+\frac{2}{\raisebox{1ex}{$\alpha $}\!\left/ \!\raisebox{-1ex}{$\beta $}\right.}} $$

(1)

Lyman-Kutcher-Burman models were used to estimate the risk for acute esophageal toxicity greater than grade 2 by RTOG scoring, symptomatic pneumonitis of all grades, and pericarditis or pericardial effusion, using the Mohan formulation (Eqs. 2, 3 and 4) [15].

$$ NTCP=\frac{1}{\sqrt{2\pi }}\kern0.5em {\int}_{-\infty}^t{e}^{\frac{-{x}^2}{2}}\ dx $$

(2)

$$ t=\frac{D_{eff}-{TD}_{50}}{m{TD}_{50}} $$

(3)

$$ {D}_{eff}={\left(\sum \limits_i{v}_i{D_i}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$n$}\right.}\right)}^n $$

(4)

In this model, TD₅₀ is the uniform dose to an organ associated with 50% complication risk, m describes slope of the dose-response, n describes volume effect, and D_i is the dose (EQD₂) to fractional volume, v_i. NTCP were estimated using published parameters by Belderbos et al. for esophagitis (TD₅₀ = 47 Gy, n = 0.69, m = 0.36), Semenko et al. for pneumonitis (TD₅₀ = 29.9 Gy, n = 1, m = 0.41) and Martel et al. for pericarditis (TD₅₀ = 50.6 Gy, n = 0.64, m = 0.13) [16,17,18]. Additional model parameters published by Chapet, Zehentmayr, Seppenwoolde, and Burman were also evaluated for comparison (Additional file 1: Table S1) [19,20,21,22].

To our knowledge, there are no published model parameters for determining risk for proximal bronchial tree toxicity. We present risk estimates based on average minimum dose (EQD₂) to 1cm³, 2cm³, and 3cm³ of the proximal bronchial tree for each treatment regimen compared against published probability curves for grade 4 or 5 proximal bronchial tree toxicity at two years, which included outcomes such as fatal hemoptysis [23].

TCP was calculated using the Martel model (Eq. 5) which describes the dose-response curve for NSCLC and estimates progression-free survival at 30 months [24].

$$ TCP\ (D)=\frac{1}{1+{\left(\raisebox{1ex}{${D}_{50}$}\!\left/ \!\raisebox{-1ex}{${D}_i$}\right.\right)}^{4\gamma }} $$

(5)

In this model, D_i is the uniform dose irradiated to a fractional volume, D₅₀ is the dose needed to achieve a 50% probability of tumor control, and γ is the normalized slope of the sigmoid-shaped dose response curve at D₅₀. Input parameters published by Fowler et al. (D₅₀ = 70 Gy, γ=1.94) were used to improve applicability of the Martel model to SABR regimens [25].

The probability for uncomplicated tumor control (UTCP) was calculated, as previously described, to quantify the most desirable treatment regimen (Eq. 6) [26].

$$ UTCP= TCP\ast \left(1- NTCP\right) $$

(6)

Results

Relevant doses to the proximal bronchial tree for each treatment regimen are reported in Table 2. When dose to the organs at risk was prioritized, PTV coverage by the 100% isodose line for the 50 Gy in 5 fraction regimen was reduced from 96.0% (95.0–100) to 67.8% (19.6–89.2); for 60 Gy in 8, from 96.0% (95.0–99.9) to 68.5% (18.8–89.5); and for 60 Gy in 15 fractions, coverage by the 95% isodose line was reduced from 97.3% (95.5–100) to 91.6% (82.0–97.8). One patient with a large (85.8 cm³ PTV) ultra-central tumor located close to the spinal cord required significant compromise of tumor coverage in order to respect the dose constraints for organs at risk. The range of PTV coverage for the other nine patients ranged from 56.8–89.2% for 50 Gy in 5, 59.4–89.5% for 60 Gy in 8 and 82.0–97.8% for 60 Gy in 15 fractions.

Table 2 Doses received by the proximal bronchial tree across competing treatment plans

Cumulative dose-volume relationships for PTV, esophagus, healthy lung, and proximal bronchial tree are shown in Fig. 2 for the six competing treatment plans. SABR, when delivered as 50 Gy in 5 fractions or 60 Gy in 8 fractions, required substantial compromise of PTV coverage to satisfy tolerance doses to OAR. A comparison of the OAR-prioritized with PTV-prioritized plans illustrates an appreciable difference in dose to the proximal bronchial tree and esophagus for all fractionation regimens, whereas dose to the lung was similar regardless of plan or prioritization.

Fig. 2

Dose-volume relationships in 2 Gy per fraction for competing dose-fractionation regimens, based on α/β =3 for OAR and α/β = 10 for tumor. Data are graphed at 50 cGy resolution and are expressed as the median (PTV) or average (OAR) fractional volumes to illustrate the relative differences in PTV coverage and OAR sparing between the treatment plans

A simple evaluation of EQD₂ (based on α/β = 10) for each fractionation schedule, if homogeneous dose is assumed, results in 87.5 Gy, 83.3 Gy, and 70 Gy to tumor for 60 Gy in 8 fractions, 50 Gy in 5 fractions, and 60 Gy in 15 fractions, respectively. The derived TCP values are 85, 80, and 50, suggesting TCP is most favorable for 60 Gy in 8 fractions.

The NTCP and TCP results are presented in Table 3. The mean probability of acute grade ≥ 2 esophageal toxicity, pneumonitis, and pericarditis or pericardial effusion was less than 5% in all dose-fractionation regimens, regardless of PTV or OAR prioritization. If PTV coverage was not compromised with SABR, the mean estimated risk of grade 4 or 5 toxicity to the proximal bronchial tree based on D_1cc was high at 49.7% and 42.7% for 50 Gy in 5 fractions and 60 Gy in 8 fractions, respectively, and exceeded 70% for 4 of 10 patients. The PTV-prioritized 60 Gy in 15 fraction regimen had a risk of 4.0%. When dose to OAR was prioritized, the risk of toxicity to the proximal bronchial tree was reduced to 1% or less for all 3 fractionation schemes. This compromise resulted in a reduction in mean TCP compared to PTV-prioritized plans from 92.9% to 60.3% for 50 Gy in 5 fractions, 92.4% to 65.7% for 60 Gy in 8 fractions, and 52.0% to 47.8% for 60 Gy in 15 fractions.

Table 3 Radiobiological modeling results for tumor control probability and normal tissue complication probability across competing dose-fractionation regimens

The most favourable UTCP regarding proximal bronchial tree injury was 65.7% with 60 Gy in 8 fractions prioritizing OAR. Notably, the previously mentioned case with substantial compromise of PTV coverage had especially low UTCP for this regimen; when excluded, the favourability of this treatment plan increases with minimum UTCP greater than 44% for the remaining cases and 69% on average. The mean UTCP for the other treatment plans were 46.7%, 53.5%, and 49.9% for 50 Gy in 5 fractions, 60 Gy in 8 fractions, and 60 Gy in 15 fractions, respectively, when PTV coverage was prioritized. When OARs were prioritized, the UTCP increased to 60.3% for 50 Gy in 5, but decreased to 47.2% for 60 Gy in 15 fractions.

The mean IGTV and PTV volumes were 30.1 cm³ (range 7.1–83.1 cm³) and 78.8cm³ (25.0–171.5 cm³), respectively. The mean volume of IGTV and PTV overlap with the proximal bronchial tree were 0.2 cm³ (0.01–0.5 cm³) and 1.1 cm³ (0.4–2.4 cm³), respectively. A comparison of PTV overlap volume with likelihood of proximal bronchial tree toxicity suggests a trend for low risk (< 10%) across all overlap volumes studied (up to 2.5 cm³) for 60 Gy in 15 fractions. In contrast, substantially greater variability existed for the other dose-fractionation regimens and risk may exceed 70% when overlap volume is greater than 1.5 cm³ (Fig. 3).

Fig. 3

Relationship between likelihood of proximal bronchial tree injury and the tumor overlap volume. The risk of grade 4 or 5 toxicity to the proximal bronchial tree based on D_1cc vs. PTV overlap volume with the proximal bronchial tree for plans that prioritize PTV coverage

Discussion

Consideration of the competing risks for local tumor failure and acute treatment toxicity is required when using hypofractionated radiotherapy or SABR to treat ultra-central lung tumors. This planning study suggests that SABR for ultra-central tumors is feasible, and a risk-adapted approach of 60 Gy in 8 fractions prioritizing OAR tolerance over PTV coverage may provide an acceptable balance between tumor control and toxicity. Although PTV coverage was compromised in this regimen, TCP was reasonable and risk to the proximal bronchial tree was negligible; however, given that proximal bronchial tree injury is correlated with high doses to small volumes (1-3 cm³), it should be emphasized that favorable toxicity heavily relies on good image-guidance to ensure the high-dose gradient adequately spares this structure. Alternatively, if treatment set up is difficult, or tumors are large and multiple OARs are a concern, a conservative approach of 60 Gy in 15 fractions prioritizing PTV coverage is reasonable, acknowledging this provides safe doses to OAR at the cost of low TCP compared to other regimens.,

There are a few important limitations to this work. First, only three dose-fractionation regimens were compared and there are other potential dose fractionation schemes which may be appropriate in this setting; however, the three regimens selected for this study are commonly used in practice. Secondly, the small sample size in this exploratory study precludes rigorous statistical analysis and a larger confirmatory study is needed; still, given the paucity of data for this clinical situation, we believe this study provides useful information for physicians treating ultra-central lung tumors. In addition, while careful consideration was given to select models that were relevant to our cohort and the practice at our institution, interpreting results from model-based predictions should be done cautiously as these estimates are derived from outcome data at other centers and can be sensitive to host factors. The model-based predictions presented here are meant to give an estimate for tumor control and toxicity; further prospective clinical research is needed to better assess the true therapeutic ratio of different dose fractionation schemes for ultra-central lung tumors and inform clinical decision-making.

To prospectively address this clinical issue, the Canadian Pulmonary Radiotherapy Group (CAPRI) has launched the Stereotactic body radiotherapy for ultra-central non-small cell lung cancer: A safety and efficacy trial (SUNSET trial, NCT03306680). SUNSET is a multi-centre phase 1 dose escalation study using the continuous re-assessment method, to determine the maximally tolerated SABR dose associated with 30% or lower rate of grade ≥ 3 toxicity within 2 years. This trial will investigate a dose of 60 Gy starting in 8 fractions and escalate to 5 fractions (with de-escalation down to 15 fractions, if necessary). Our modeling results suggest only 4 of the 10 patients included in our study would meet the SUNSET outcome objectives at level 1 (60 Gy in 8 fractions). Additional data from the EORTC LungTech phase II trial of SABR (60 Gy in 8 fractions) for central early stage lung cancer are also awaited (NCT01795521). Until the efficacy and toxicity analysis from these trials are complete, the best dose-fractionation regimen for ultra-central tumors remains unclear.

Conclusions

In balancing the competing concerns over acute toxicity and local control when treating an ultra-central tumor with SABR, we propose a regimen of 60 Gy in 8 fractions, compromising PTV coverage to respect OAR dose. Alternatively, 60 Gy in 15 fractions represents a conservative and safe approach with the caveat of lower tumor control. Further prospective research is needed to determine the true dose-response relationships associated with ultra-central tumors and elucidate the optimal dose-fractionation regimen.