 Research
 Open Access
 Published:
Investigating the locoregional control of simultaneous integrated boost intensitymodulated radiotherapy with different radiation fraction sizes for locally advanced nonsmallcell lung cancer: clinical outcomes and the application of an extended LQ/TCP model
Radiation Oncology volume 15, Article number: 124 (2020)
Abstract
Background
To investigate the locoregional progressionfree survival (LPFS) of intensitymodulated radiotherapy (IMRT) with different fraction sizes for locally advanced nonsmallcell lung cancer (LANSCLC), and to apply a new radiobiological model for tumor control probability (TCP).
Methods
One hundred and three LANSCLC patients treated with concurrent radiochemotherapy were retrospectively analyzed. Factors potentially predictive of LPFS were assessed in the univariate and multivariate analysis. Patients were divided into group A (2.0 ≤ fraction size<2.2Gy), B (2.2 ≤ fraction size<2.5Gy), and C (2.5 ≤ fraction size≤3.1Gy) according to the tertiles of fraction size. A novel LQRG/TCP model, incorporating four “R”s of radiobiology and Gompertzian tumor growth, was developed to predict LPFS and compared with the classical LQ/TCP model.
Results
With a median followup of 22.1 months, the median LPFS was 23.8 months. Fraction size was independently prognostic of LPFS. The median LPFS of group A, B and C was 13.8, 35.7 months and not reached, respectively. Using the new LQRG/TCP model, the average absolute and relative fitting errors for LPFS were 6.9 and 19.6% for group A, 5.5 and 8.8% for group B, 6.6 and 9.5% for group C, compared with 9.5 and 29.4% for group A, 16.6 and 36.7% for group B, 24.8 and 39.1% for group C using the conventional LQ/TCP model.
Conclusions
Hypofractionated IMRT could be an effective approach for dose intensification in LANSCLC. Compared with conventional LQ model, the LQRG model showed a better performance in predicting followup time dependent LPFS.
Background
Concurrent chemoradiotherapy (CCRT) has long been established as the standard therapy for unresectable locally advanced nonsmallcell lung cancer (LANSCLC), while the locoregional control (LRC) and overall survival (OS) were suboptimal. Patients still experienced a locoregional progression rate of approximately 30 ~ 40% at 3 years in most of the randomized trials [1]. After the impact of radiation dose for NSCLC was established [2], efforts to improve the LRC have been focused on increasing the total irradiation dose. A major finding in RTOG0617 was that total dose escalation of 74Gy with 2Gy per fraction failed to improve LRC or OS [3]. Therefore, hypofractionation radiotherapy could be another attempt to achieve better local control for LANSCLC.
Hypofractionation radiotherapy enables the delivery of an increased biologically effective dose (BED) without extending the overall treatment time (OTT). Image guided radiotherapy (IGRT) and intensitymodulated radiation therapy (IMRT) conform the radiation dose to the tumor and spare adjacent critical organs [4]. Simultaneous integrated boost (SIB)IMRT has been developed recently as a dose intensification technique, which delivers different dose prescriptions to different target volumes with promising local control. However, there are limited data on hypofractionated SIBIMRT concurrent with chemotherapy in LANSCLC, and the feasibility and efficacy needs to be investigated.
Other than that, setting up tumor control probability (TCP) model remained crucial for developing new fraction scheme for radiotherapy. The linearquadratic (LQ) model provides the theoretical background for estimating TCP via the calculation of BED. Nevertheless, the conventional LQ/TCP model was proved to be ineffective in modeling the radiobiological response of hypo and hyperfractionated IMRT for lung cancer. We recently developed a novel LQ model incorporating cell repair, redistribution, reoxygenation, regrowth and Gompertzian tumor growth (LQRG/TCP model) to improve the modeling accuracy for the local tumor control of small cell lung cancer (SCLC) [5]. Compared with LQ/TCP model and its following modifications, LQRG which incorporates all the four “R”s of radiobiology, might be more realistic assumptions for the prediction of TCP. It is interesting to explore whether the LQRG/TCP model can as well explain the clinical effect of hypofractionated IMRT with different fraction sizes in the case of NSCLC. As to our knowledge, little work has been done to modeling the time dependent locoregional progressionfree survival (LPFS) in the radiation of NSCLC and use it as a theoretical foundation to interpret the fraction size dependent clinical results of the IMRT. Some authors tried to provide more clinically relevant TCP modeling of tumor response to radiotherapy in NSCLC [6,7,8]. However, they were mainly restricted to stereotactic body radiation therapy (SBRT) and did not incorporate the redistribution and reoxygenation mechanisms of clonogenic cell and the interplay of radiation treatment between fractions into the LQ formalism.
This study retrospectively analyzed the LPFS of patients with LANSCLC, who were treated by SIBIMRT concurrent with chemotherapy at different doses per fraction, and provided evidence in support of hypofractionated SIBIMRT as an alternative of conventional fractionation. An recently developed extended LQ model, LQRG, was applied to account for the LPFS of the IMRT for NSCLC.
Methods
Acquisition of clinical data
We retrospectively reviewed the records of 103 LANSCLC patients treated from January 2012 to May 2015 in our Center. Patients included in our study had: 1.pathologically confirmed LANSCLC (stage IIIA/B); 2.treatment with SIBIMRT and concurrent chemotherapy; 3.ECOG performance status score 0–2 and life expectancy of more than 6 months. Patients who received prior thoracic surgery or thoracic radiation, or had other coexisting malignancy were excluded. The 7th edition of American Joint Committee on Cancer staging system was used to stage the diseases. This study was approved by the ethics committee of our center. Written informed consent to use clinical data was obtained.
Radiotherapy and concurrent chemotherapy
Patients were immobilized and simulated according to the standard protocol of lung cancer in our center [9]. Gross tumor volume (GTV) was defined as visible primary tumors and involved lymph nodes on CT and/or PET scans. The criteria of lymph node positivity included: short axis size ≥10 mm on pretreatment CT scan, reported positive on the pretreatment PET scan, or biopsy positive on mediastinoscopy or endobronchial ultrasoundguided biopsies. The GTVs were composite volumes from CT scans of all breathing phases. Clinical target volume (CTV) included the primary lung tumor with a 0.6 cm margin, the ipsilateral hilum, involved lymph nodes region and highrisk lymph nodal regions (for example, if 4R was involved, 4 L and 2R were included in CTV) (Supplementary Figure 1). The planning target volume for GTV (PTV1) and CTV (PTV2) covered the GTV and CTV with a 0.5cm margin, respectively.
SIBIMRT technique was used to deliver different dose prescriptions to PTV1 and PTV2 in the same treatment fraction (for example, 65Gy to PTV1 and 50Gy to PTV2 in 25 fractions). The dose and fractions were decided by planning clinic discussion and limited by normal tissue constraints. The prescription dose to PTV1 was 60 ~ 70Gy, which was the typical radical dose for LANSCLC. The prescription dose to PTV2 was 4554Gy, which was a standard prophylactic dose. Dose constraint for organs at risk (OAR) included: the maximum spinal cord dose < 46 Gy, mean lung dose (MLD) < 20Gy and the lung volumes irradiated above 20Gy (V20) < 35%, the maximum esophagus dose < 110% of prescription dose, and V40 of the heart < 50%. The actual radiation dose to PTV was defined as the minimum dose received by 95% of PTV. Radiotherapy was delivered using daily cone beam CT image guidance. The concurrent chemotherapy were platinumbased doubleagent regimen.
Followup
Patients were followed up 2 months after CCRT, then every 3 months for the first 3 years and every 6 months thereafter. Routine surveillance imaging included chest, abdomen CT and brain MRI. Patients underwent bone scan when suspected for bone metastasis, and PET scan when suspected for systemic progression. Locoregional progression, distant metastasis and survival status were recorded. LPFS was calculated from the last date of radiotherapy to the date of locoregional progression or the date of last visit before November 30th, 2018. Infield failure was defined as recurrent tumors mapped within PTV1. Common Terminology Criteria for Adverse Events Version 4.0 (CTCAE 4.0) was used to evaluate the toxicities.
Statistical analyses
Continues variables were presented as median and range. The survival analyses were performed using the KaplanMeier method. Clinical and dosimetric factors potentially predictive of LPFS were assessed using Cox proportional hazards modeling. All statistical analyses were performed with SPSS 19.0 software (IBM), and a P value< 0.05 was considered statistically significant.
Tumor control probability modeling
The conventional LQ model combined with a Binomial TCP model were described in Supplementary 1. The construction of LQRG model combined with a Gaussian TCP model are briefly presented below.
The model name “LQRG” stands for the Linear Quadratic model with 4R effects and Gompertzian tumor growth. It takes account of the four “R”s effects, and add features of delayed regrowth [10] and Gompertzian tumor growth [11].
The surviving fraction (SF) of clonogenic cells is expressed as:
where G(τ) is the generalized LeaCatcheside function [11, 12], defined by Eq. 6; τ_{R} is the average DNA repair time; τ_{S} is average resensitization time; σ is the variance of the Gaussian distribution of random variable α, immediately after an acute irradiation; T is the treatment time; T_{k} is the delayed time for regrowth; t is the elapsed time since the end of treatment. The first and second terms describe cell killing by onetrack and twotrack action (and possible repair), respectively. The third term refers to intercellular diversity of radiosensitivity and resensitization [13]. The forth term refers to the delayed tumor regrowth, while the last term involving parameter δ refers to the Gompertzian tumor growth after radiation [14]. Although quite different phenomena, redistribution and reoxygenation do share a common outcome [15], a postirradiation increase in the sensitivity of cells that survive an initial or partial exposure. Following Brenner et al. [13], we denote this common outcome resensitization and incorporates these two “R”s in the third term of the exponent of SF.
The TCP is:
with
where σ_{k} is the Gaussian width for the distribution of tumor cell numbers [7]. K_{0} is the initial average number of clonogenic cells and K_{cr} is a critical number under which the tumor would be controlled for an individual patient [7].
BED for the LQRG model is then expressed as:
To compare the α/β ratio at the same sense of the conventional LQ model, an effective β, β_{eff,} was defined for the LQRG model.
where n was the number of fractions. With this definition of β_{eff,} the quadratic term of D^{2} in SF of Eq. 5 can be cast into the same form as in Eq. 2: β_{eff}Dd. This new parameter will be utilized in later comparison of α/β between the conventional LQ model and new LQRG model in a fair way.
The relations between generalized LeaCatcheside function G(τ) and the recovery time constant τ were computed using numerical integration method. Dose and volume of GTV were extracted from the dosevolume histograms (DVHs) of radiation plans. The clinical endpoint linked to the TCP for model fitting was the time dependent LPFS. The calculated TCP(t) is fitted to the measured LPFS(t), where t is the time from the end of radiotherapy to the time of follow up. The least chisquare method was utilized in the fitting [16]. The free parameters in the function were determined by minimizing the cost function, which was the mean absolute error, using a gradient descent method. The absolute and relative fitting errors for LPFS was calculated using formulas in Additional file 1: Appendix A of our previous study [5]. The confidence intervals for estimated parameters were determined following standard procedures in nonlinear regression [17]. The maximum log likelihood and Akaike information criterion (AIC) was computed to compare the fitting goodness and synthetic quality of the models (see Additional file 1: Appendix B of our previous study [5]). In accordance with Additional file 1: Appendix A and B of [5], the minimization of the chisquare cost function, the mean absolute error (MAE), is equivalent to the maximization of log likelihood. Leaveoneout cross validation based on the timedependent LPFS samples was performed to compare the generalization capabilities of the models. All the aforementioned calculations in tumor control probability modeling were conducted on a cloudbased clinical data service platform OncoEvidance™ v1.0 (Homology Medical, Suzhou, China, 2019) and the fitted model of LQRG was embedded in the radiobiology module of Zeus Cloud TPS™ (Treatment Planning System) v1.0 (Homology Medical, Suzhou, China, 2019).
Results
Patient characteristics
The median actual radiation does to PTV1 was 63.8Gy (55.6 ~ 69.3) with a median fraction size of 2.4Gy (2.0 ~ 3.1). Patients were divided, according to tertiles of fraction size, into group A (2.0 ≤ fraction size<2.2Gy, n = 34), group B (2.2 ≤ fraction size<2.5Gy, n = 34) and group C (2.5 ≤ fraction size≤3.1Gy, n = 35). Patient demographics and disease characteristics for three groups are listed in Table 1. The radiation dosimetry data are listed in Supplementary Table 1. The relative dose intensity of concurrent chemotherapy is listed in Supplementary Table 2. Sex, age, ECOG PS, histology, clinical stage, concurrent chemotherapy, as well as GTV volumes, showed similar distribution among three groups.
Locoregional control and survival
With a median followup of 22.1 months (range, 5.4 ~ 64.3), the median OS was not reached. The 2year and 5year OS rates were 66.1 and 23.1% respectively.
The median LPFS was 23.8 months. LPFS rates dropped from 100% at time 0 to 73.9% at 1 year, 48.6% at 2 years, and 44.7% at 3 years). A total of 47 patients developed locoregional progression.
We assessed the impact of the following variables on LPFS: age, sex, ECOG PS, histology, tumor stage, GTV volume, radiation dose, conventional BED by LQ model, fraction size, OTT and concurrent chemotherapy (Table 2). In univariate analysis, radiation fraction size (p = 0.003; 2.0 ~ 2.2 vs. 2.5 ~ 3.5Gy, hazard ratio [HR] = 3.609; 2.2 ~ 2.5 vs. 2.5 ~ 3.5Gy, HR = 1.591) and OTT (p = 0.008, ≥5.3 weeks vs. < 5.3 weeks, HR = 2.246) were associated significantly with LPFS. Patients in group A and B had a median LPFS of 13.8 and 35.7 months, respectively. The median LPFS of group C was not reached (Fig. 1). LPFS rate at 2 years was 24.0, 55.8 and 67.8% for patients of group A, B and C, respectively. Patients with an OTT < 5.3 weeks had a longer LPFS than that with an OTT ≥5.3 weeks weeks (48.7 vs. 15.9 months).
Patterns of failure
Seventyfive patients (72.8%) had progression of disease: 25 (24.3%) had local locoregional progression alone, 28 (27.1%) had distant progression alone, and 22 (21.4%) had both locoregional and distant progression. The locoregional progression sites included ipsilateral lung (n = 32), mediastinum/supraclavicular lymph nodes (n = 6) and both (n = 9). There were a total of 42 (40.8%) infield locoregional failures and 5 (4.9%) outfield locoregional failures.
Toxicities
Radiation pneumonitis and esophagitis were the most common acute toxicities. Grade ≥ 3 pneumonitis and esophagitis were observed in 6 (5.7%) and 18 (17.1%) patients, respectively. There was no statistical difference in the incidence of complications among three groups. Lifethreatening hemoptysis was observed in one patient irradiated at 63Gy in 21 fractions (3.0Gy per fraction) 2 months after RT. More information on toxicities are provided in Supplementary Table 3.
Modeling
The mean total GTV doses were 62.92Gy, 64.07Gy and 62.89Gy for groups A, B and C, respectively. The mean GTV volumes were 105.83, 102.21 and 108.51cm^{3}, respectively.
Figure 1 presents the results of TCP (or LPFS) model fitting within a duration of 70 months. Compared with the excellent fit of the LQRG model, fitting of the LQ model failed immediately after 20 months. Using the new LQRG/TCP model, the average absolute and relative fitting errors for LPFS were 6.9 and 19.6% for group A, 5.5 and 8.8% for group B, 6.6 and 9.5% for group C, compared with the average absolute and relative fitting errors for LPFS as 9.5 and 29.4% for group A, 16.6 and 36.7% for group B, 24.8 and 39.1% for group C using the conventional LQ/TCP model. The values of model parameters and their confidence intervals are detailed in Table 3. Supplementary Table 4 lists the fitting errors of LQ model and LQRG model for the relation between LPFS and followup time after treatment. Figure 2 shows the BEDLPFS relation predicted at the 13th and 48th month, and suggests a better agreement between data and prediction in the LQRG model.
Supplementary Table 5 lists the intermediate results of the calculation of LPFS at the end of the treatment. These intermediate items were extracted from the exponent part of the SF from both models. Supplementary Figure 2 shows the relation between generalized LeaCatcheside function G(τ) and the recovery time constant τ for both schemes. Table 4 presents the metrics characterizing the quality of the LQ and LQRG models.
Discussion
The standard radiation scheme for LANSCLC consists of a total dose of 60 ~ 70Gy delivered in a fraction size of 2Gy per day. Various treatment schemes have been tried in an effort to improve local control by either dose escalation or shortening the total treatment time through unconventional fractionation [18]. The results from RTOG 0617 suggested that the balance between treatmentrelated toxicities and dose escalation must be taken into account seriously [3]. Asian patients physically have smaller lungs volume than western patients, which make it hard to escalate the total dose beyond 70Gy. Therefore, we tried to use SIBIMRT technique to escalate the fraction size. Our results showed that in the range of 60 ~ 70Gy, total radiation dose was not prognostic of LPFS. LPFS was improved by the escalation of radiation fraction size (HR, 0.590; CI, 0.348 ~ 1001). It implied that hypofractionated SIBIMRT could be an alternative of conventional fractionation.
Data on the efficacy of hypofractionated radiotherapy for LANSCLC remain limited. A systematic review of the literature about radicalintent hypofractionated RT for LANSCLC included 33 studies [19]. The RT dose ranged from 45.0 to 85.5 Gy, and was delivered at 15 ~ 35 fractions at a dose per fraction ranging from 2.3 to 3.5Gy. OS was found to be associated with tumor BED. A phase II study showed that escalating radiation dose to the FDGavid tumor detected by midtreatment PET provided a favorable LRC in stage II or III NSCLC [20]. RT was delivered in 30 daily fractions up to a total physical dose of 86Gy. It achieved a promising 2year infield LRC rate of 82% and overall LRC rate of 62%. Our study compared the LRC of different radiation fraction size and showed that hypofractionation (over 2.5Gy per fraction) achieved a longer LRC time than conventional fractionation. The possible advantage of hypofractionation includes the escalation of BED and shortening of treatment time, which is believed to limit tumor repopulation [21]. A work by Machtay et al. examined the relationship between BED and treatment outcome in conventionally fractionated and hyperfractionated regimens, where a 1Gy increase in BED resulted in a 3% improvement in LPFS [22]. It has been reported that there is a loss of absolute 3year survival rate of 1.6% per day OTT prolongation beyond 6 weeks [23]. Consistent with these previous findings, our data showed that fraction size and OTT were associated significantly with LPFS, and fraction size turned out to be an independent prognostic factor.
Calculated by the classical LQ model, the order of the BEDs of all three patient groups correctly predicted the height order of the LPFS curves, but the overall fitting of LPFS for NSCLC patients from followup data failed again just as in the case of SCLC, thereby reaffirming its incapability of TCP prediction. The LQRG and the new TCP model, however, succeeded in fitting and predicting the followup time dependent LPFS, which, as pointed out in the case of SCLC [5], contributed mainly to two model improvements: namely, one accounting for the 4R effect of radiobiology and the other for more reasonable assumptions for the clonogenic regrowth and cell number distribution.
As in the work on SCLC [5], an additional analysis of the intermediate results of the SF and BED further confirmed the previously revealed mechanisms behind the LQRG models, namely the positive difference between the three schemes made by the resensitization term, the slowing down of the rapid decaying rate of the LPFS as time proceeded on by the Gompertzian term and the reformulation of TCP with K_{cr}, and the characteristic distinctions between the three schemes caused by the gaps among the three curves of LeaCatchside function. Once again, these distinctions from LeaCatchside function were embodied in the second and third term of the SF in (Eq. 5) and subtly balanced by other factors to contribute to the final comparative clinical results. For example, with the optimized value of the τ_{R} = 6.32 h, G (τ_{R}) in LQRG model were 0.036, 0.040 and 0.049 for groups A, B and C, respectively, which were larger than those corresponding values of 0.033, 0.037 and 0.045 in conventional LQ model. This implied the stronger cell killing effect by twotrack actions in LQRG model which might have more accurately represented the interplay between dose fractions through its precise numerical calculation of LeaCatchside function. Moreover, the models from both conventional LQ and LQRG revealed the subtle effect of the duration of treatment through the term representing the delayed regrowth of the tumor ln2(TT_{k})/τ_{P,} which slightly favored short treatment duration with larger fraction size in hypofractionated IMRT as had been anticipated in clinical settings. Consequently, the BED of the LQRG model for groups A, B and C were 85.93Gy, 91.03Gy and 95.35Gy respectively, which was in good agreement with the clinical observation that the 2.5–3.1Gy per fraction concurrent with chemotherapy showed the best LPFS. Compared with the BED of the conventional LQ model, the relative difference between the BEDs was amplified in LQRG model, which might better explain the additional longterm benefit for schemes with larger fraction sizes.
A close examination of the optimized model parameters revealed that they fell well into the published data range. The α was 0.067, close to the value for NSCLC in previous study [7]. To comparatively evaluate the value of β with published data using conventional LQ model, we recast the Eq. 5 into the form of LQ model by introducing an effective parameter β_{eff} in Eq. 6. Rather than compute the α/β ratio as in conventional LQ model, we computed the α/β_{eff} ratio in LQRG model to compare with its conventional counterpart, α/β. The optimized value of α/β_{eff} ratio in LQRG model was 10.06, which agreed well with the published consensus of around 10 for NSCLC [24]. The cell number doubling time τ_{p} was 84.55 days, closed to the clinically agreed value of nearly 3 months for NSCLC [7, 25]. The delay regrowth time T_{k} is 533.5 h or 22.2 days, close to the clinically agreed value of nearly 1 months [7, 25]. The repair and resensitization time constants were 6.32 h and 3.04 h respectively, within the range of fitted laboratorybased experimental data from 1 h to tens of hours [13].
The preferred log likelihood of LQRG model was 117.14, which was slightly higher than that of LQG (100.52) and substantially higher than that of LQ (1.7), indicating the reasonable improvements of the fitting capability of LQRG model over those of LQG and LQ. In order to evaluate the risk of overfitting for LQRG model with six more variables than that of conventional LQ and TCP model, AIC was developed for all the three models. AIC, which acts as an operational way of trading off the complexity of an estimated model against how well the model fits the data, is a measure of the fitting quality of an estimated statistical model, with lower value indicating higher model quality. AIC not only rewards goodness of fit, but also penalizes the number of estimated parameters which is closely related to overfitting. The fact that AIC of the fitting of LQRG model to the data was − 212.27, slightly lower than that of LQG model (− 183.03) and substantially lower than that of LQ model (7.86) validated LQRG’s superior synthetic model quality with reasonably low overfitting risk. In addition, leaveoneout cross validation for the three models showed 88.2% average prediction accuracy for LQRG model, with 85.6% for LQG and 64.7% for LQ respectively. This result of leaveoneout cross validation again demonstrated LQRG’s sound predictive power and a better generalization capability with acceptable overfitting risk. Compared to conventional LQ and TCP model, although six additional variables were added to the new models, the mechanistically driven nature of these variables could potentially interpret the outstanding balance between fitting accuracy and generalization capability. As to our knowledge, this is the first theoretical work to accurately predict the followup time dependence of LPFS for NSCLC with hypofractionated SIBIMRT. Since the proposition of LQRG model in SCLC [5], this work has applied the same LQRG modelling framework in a different type of lung cancer  NSCLC, and confirmed its universality in followup time dependent LPFS modelling for lung cancer. In this study, death was not considered as an event in the LPFS estimation if the patient was not dead of locoregional progression. Distant metastasis was the main pattern of failure that affected survival in locally advanced NSCLC [26]. Our data showed 50 patients (48.5%) developed distant metastasis during followup. Therefore, considering death as an event in the LPFS might overestimate the locoregional failure and affect the following model fitting.
This study was limited by its retrospective nature. The dose and fractions were decided by planning clinic discussion and limited by normal tissue constraints. The prescription might have been influenced by doctor’s preference, tumor volume, performance status and pulmonary function. Some factors, such as tumor volume and PS, might influence the local control. Therefore, we included these potential confounding factors in univariate and multivariate analysis to minimize the confounding effect. This study was also limited by the lack of prospective collection of toxicity data. There might be an underestimation of toxicity. Therefore, dose escalation through hypofractionated SIBIMRT should always be performed with caution of toxicity. This study was also limited by a modest sample size and events which might cause overfitting of the model. A larger sample size with a longer followup time is warranted.
Conclusions
Our study shows improved locoregional control in patients treated with an escalated fraction size in the setting of concurrent chemotherapy. It supported that, for Asian patients, hypofractionated SIBIMRT could be a feasible and effective approach for dose intensification in LANSCLC. Compared with conventional LQ model, the LQRG model showed a better performance in followup time dependent LPFS modeling. A larger sample size and longer followup time was warranted before widely applying hypofractionated RT to LANSCLC in clinical practice.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 LPFS:

Locoregional progressionfree survival
 IMRT:

Intensitymodulated radiotherapy
 LANSCLC:

Locally advanced nonsmallcell lung cancer
 TCP:

Tumor control probability
 LQRG:

Linearquadratic model incorporating cell repair, redistribution, reoxygenation, regrowth and Gompertzian tumor growth
 LQ:

Linearquadratic
 LRC:

Locoregional control
 OS:

Overall survival
 BED:

Biologically effective dose
 OTT:

Overall treatment time
 IGRT:

Image guided radiotherapy
 (SIB)IMRT:

Simultaneous integrated boost intensitymodulated radiotherapy
 SCLC:

Small cell lung cancer
 SBRT:

Stereotactic body radiation therapy
 GTV:

Gross tumor volume
 CTV:

Clinical target volume
 PTV:

Planning target volume
 CT:

Computed tomography
 PET:

Positron emission tomography
 MRI:

Magnetic resonance imaging
 SF:

Surviving fraction
 DVH:

Dosevolume histograms
 AIC:

Akaike information criterion
 CCRT:

Concurrent chemoradiotherapy
References
 1.
Auperin A, Le Pechoux C, Rolland E, Curran WJ, Furuse K, Fournel P, et al. Metaanalysis of concomitant versus sequential radiochemotherapy in locally advanced nonsmallcell lung cancer. J Clin Oncol. 2010;28:2181–90.
 2.
Machtay M, Paulus R, Moughan J, Komaki R, Bradley JE, Choy H, et al. Defining localregional control and its importance in locally advanced nonsmall cell lung carcinoma. J Thorac Oncol. 2012;7:716–22.
 3.
Bradley JD, Paulus R, Komaki R, Masters G, Blumenschein G, Schild S, et al. Standarddose versus highdose conformal radiotherapy with concurrent and consolidation carboplatin plus paclitaxel with or without cetuximab for patients with stage IIIA or IIIB nonsmallcell lung cancer (RTOG 0617): a randomised, twobytwo factorial phase 3 study. Lancet Oncol. 2015;16:187–99.
 4.
Liao ZX, Komaki RR, Thames HD Jr, Liu HH, Tucker SL, Mohan R, et al. Influence of technologic advances on outcomes in patients with unresectable, locally advanced nonsmallcell lung cancer receiving concomitant chemoradiotherapy. Int J Radiat Oncol Biol Phys. 2010;76:775–81.
 5.
Li QW, Qiu B, Wang B, Zhang J, Chen L, Zhou Y, et al. Comparison of hyper and hypofractionated radiation schemes with IMRT technique in small cell lung cancer: clinical outcomes and the introduction of extended LQ and TCP models. Radiother Oncol. 2019;136:98–105.
 6.
Mehta N, King CR, Agazaryan N, Steinberg M, Hua A, Lee P. Stereotactic body radiation therapy and 3dimensional conformal radiotherapy for stage I nonsmall cell lung cancer: a pooled analysis of biological equivalent dose and local control. Pract Radiat Oncol. 2012;2:288–95.
 7.
Tai A, Liu F, Gore E, Li XA. An analysis of tumor control probability of stereotactic body radiation therapy for lung cancer with a regrowth model. Phys Med Biol. 2016;61:3903–13.
 8.
Liu F, Tai A, Lee P, Biswas T, Ding GX, El Naqa I, et al. Tumor control probability modeling for stereotactic body radiation therapy of earlystage lung cancer using multiple biophysical models. Radiother Oncol. 2017;122:286–94.
 9.
Ma L, Qiu B, Li Q, Chen L, Wang B, Hu Y, et al. An esophagussparing technique to limit radiation esophagitis in locally advanced nonsmall cell lung cancer treated by simultaneous integrated boost intensitymodulated radiotherapy and concurrent chemotherapy. Radiat Oncol. 2018;13:130.
 10.
Fowler JF. The linearquadratic formula and progress in fractionated radiotherapy. Br J Radiol. 1989;62:679–94.
 11.
Lea DE. Actions of radiations on living cells. London: Cambridge University Press; 1946.
 12.
Sachs RK, Hlatky LR, Hahnfeldt P. Simple ODE models of tumour growth and antiangiogenic or radiation treatment. Mathl Comput Model. 2001;33:1297–305.
 13.
Brenner DJ, Hlatky LR, Hahnfeldt PJ, Hall EJ, Sachs RK. A convenient extension of the linearquadratic model to include redistribution and reoxygenation. Int J Radiat Oncol Biol Phys. 1995;32:379–90.
 14.
Speer JF, Petrosky VE, Retsky MW, Wardwell RH. A stochastic numerical model of breast cancer growth that simulates clinical data. Cancer Res. 1984;44:4124–30.
 15.
Thames HD Jr, Peters LJ, Withers HR, Fletcher GH. Accelerated fractionation vs hyperfractionation: rationales for several treatments per day. Int J Radiat Oncol Biol Phys. 1983;9:127–38.
 16.
Moraru IC, Tai A, Erickson B, Li XA. Radiation dose responses for chemoradiation therapy of pancreatic cancer: an analysis of compiled clinical data using biophysical models. Pract Radiat Oncol. 2014;4:13–9.
 17.
Bates DM, Watts DG. Nonlinear regression analysis and its applications; 1988.
 18.
Mauguen A, Le Pechoux C, Saunders MI, Schild SE, Turrisi AT, Baumann M, et al. Hyperfractionated or accelerated radiotherapy in lung cancer: an individual patient data metaanalysis. J Clin Oncol. 2012;30:2788–97.
 19.
Kaster TS, Yaremko B, Palma DA, Rodrigues GB. Radicalintent hypofractionated radiotherapy for locally advanced nonsmallcell lung cancer: a systematic review of the literature. Clin Lung Cancer. 2015;16:71–9.
 20.
Kong FM, Ten Haken RK, Schipper M, Frey KA, Hayman J, Gross M, et al. Effect of Midtreatment PET/CTadapted radiation therapy with concurrent chemotherapy in patients with locally advanced nonsmallcell lung cancer: a phase 2 clinical trial. JAMA Oncol. 2017;3:1358–65.
 21.
Withers HR, Taylor JM, Maciejewski B. The hazard of accelerated tumor clonogen repopulation during radiotherapy. Acta Oncol. 1988;27:131–46.
 22.
Machtay M, Bae K, Movsas B, Paulus R, Gore EM, Komaki R, et al. Higher biologically effective dose of radiotherapy is associated with improved outcomes for locally advanced nonsmall cell lung carcinoma treated with chemoradiation: an analysis of the radiation therapy oncology group. Int J Radiat Oncol Biol Phys. 2012;82:425–34.
 23.
Fowler JF, Chappell R. Nonsmall cell lung tumors repopulate rapidly during radiation therapy. Int J Radiat Oncol Biol Phys. 2000;46:516–7.
 24.
Steel GG. Basic clinical radiobiology, 3rd edn. London: Arnold; 2002.
 25.
Fowler JF. Biological factors influencing optimum fractionation in radiation therapy. Acta Oncol. 2001;40:712–7.
 26.
Chun SG, Hu C, Choy H, Komaki RU, Timmerman RD, Schild SE, et al. Impact of intensitymodulated radiation therapy technique for locally advanced nonsmallcell lung cancer: a secondary analysis of the NRG oncology RTOG 0617 randomized clinical trial. J Clin Oncol. 2017;35:56–62.
Acknowledgements
Not applicable.
Funding
This work was supported by Science and Technology Planning Project of Guangdong Province, China (Grant Number 2016A020215190, 2016ZC0030),the National Key R&D Program of China [Grant Number 2018YFC0116800], Suzhou Municipal Science and Technology Program [Grant Number SYS2018009], Suzhou New and HighTech Zone Healthcare Program [Grant Number 2017Z001] and Huzhou Healthcare Program [Grant Number 2018GY09].
Author information
Affiliations
Contributions
All authors should have made substantial contributions to all of the following: (1) HL, YZ contributed to the conception and design of the study. (2) BQ, QWL, XLA, JZ, LC and JYG contributed to acquisition of data. (3) BQ, QWL, BW, JH, ZFZ, GHY, MJ, HHJ and CL Analyzed and interpreted data. (4) BQ, QWL, XLA drafted the article. (5) HL and YZ made the final approval of the version to be submitted. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The treatment of all patients was discussed by the thoracic multidisciplinary team, and participant information collection was approved by the Ethics Committee of Sun Yatsen University Cancer Center. Written informed consent was obtained from the patient for the publication of this report.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing financial interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Qiu, B., Li, Q.W., Ai, X.L. et al. Investigating the locoregional control of simultaneous integrated boost intensitymodulated radiotherapy with different radiation fraction sizes for locally advanced nonsmallcell lung cancer: clinical outcomes and the application of an extended LQ/TCP model. Radiat Oncol 15, 124 (2020). https://doi.org/10.1186/s1301402001555x
Received:
Accepted:
Published:
Keywords
 Locally advanced nonsmallcell lung cancer
 Locoregional progressionfree survival
 Fraction size
 Hypofractionation
 Tumor control probability model