- Open Access
Dose distribution and tumor control probability in out-of-field lymph node stations in intensity modulated radiotherapy (IMRT) vs 3D-conformal radiotherapy (3D-CRT) of non-small-cell lung cancer: an in silico analysis
Radiation Oncologyvolume 10, Article number: 178 (2015)
The advent of IMRT and image-guided radiotherapy (IGRT) in combination with involved-field radiotherapy (IF-RT) in inoperable non-small-cell lung cancer results in a decreased incidental dose deposition in elective nodal stations. While incidental nodal irradiation is considered a relevant by-product of 3D-CRT to control microscopic disease this planning study analyzed the impact of IMRT on dosimetric parameters and tumor control probabilities (TCP) in elective nodal stations in direct comparison with 3D-CRT.
Methods and materials
The retrospective planning study was performed on 41 patients with NSCLC (stages II-III). The CTV was defined as the primary tumor (GTV + 3 mm) and all FDG-PET-positive lymph node stations. As to the PTV (CTV + 7 mm), both an IMRT plan and a 3D-CRT plan were established. Plans were escalated until the pre-defined dose-constraints of normal tissues (spinal cord, lung, esophagus and heart) were reached. Additionally, IMRT plans were normalized to the total dose of the corresponding 3D-CRT. For two groups of out-of-field mediastinal node stations (all lymph node stations not included in the CTV (LNall_el) and those directly adjacent to the CTV (LNadj_el)) the equivalent uniform dose (EUD) and the TCP (for microscopic disease a D50 of 36.5 Gy was assumed) for the treatment with IMRT vs 3D-CRT were calculated.
In comparison, a significantly higher total dose for the PTV could be achieved with the IMRT planning as opposed to conventional 3D-CRT planning (74.3 Gy vs 70.1 Gy; p = 0.03). In identical total reference doses, the EUD of LNadj_el is significantly lower with IMRT than with 3D-CRT (40.4 Gy vs. 44.2 Gy. P = 0.05) and a significant reduction of TCP with IMRT vs 3D-CRT was demonstrated for LNall_el and LNadj_el (12.6 % vs. 14.8 %; and 23.6 % vs 27.3 %, respectively).
In comparison with 3D-CRT, IMRT comes along with a decreased EUD in out-of-field lymph node stations. This translates into a statistically significant decrease in TCP-values. Yet, the combination of IF-RT and IMRT leads to a significantly better sparing of normal tissues and higher total doses whereas the potential therapeutic drawback of decreased incidental irradiation of elective lymph nodes is moderate.
Radiochemotherapy is the standard treatment for inoperable, non-metastasized non-small-cell lung cancer (NSCLC). Radiotherapy planning of NSCLC underwent major changes during the past decade. Both the integration of FDG-positron emission tomography (FDG-PET) and the shift to involved-field radiotherapy (IF-RT) became a new therapeutic standard . It was shown in a series of trials that IF-RT (in the context of FDG-PET based treatment planning) does not go along with an unacceptable risk of isolated failure in ‘out-of-field’ lymph node stations but allows for further dose escalation [2–6]. The favorable low rates of ‘out-of-field’-failure with IF-RT were also attributed to a significant dose-deposition in elective nodal areas with the use of 3D-conformal radiotherapy (3D-CRT) [7, 8]. Meanwhile, a transition from 3D-CRT to intensity modulated radiotherapy (IMRT) in NSCLC takes place [9, 10]. Given the fact that IMRT is of additional benefit with respect to dose escalation and the sparing of relevant organs at risk (OAR) [11–13], it is yet unclear if the high dose conformality of IMRT implies an adverse effect on tumor control probability (TCP) for microscopic disease in the “rind”-region outside the clinical target volume. This therapeutic backdrop was the incentive to perform the presented in silico-analysis, which is a head-to-head comparison of IMRT and 3D-CRT in terms of dose distribution and hypothetical TCP in out-of-field mediastinal and hilar lymph node stations. It was carried out on a group of 41 patients who had participated in the ‘PET-PLAN pilot trial’, which examined the rate of isolated ‘out-of-field’ nodal failures with FDG-PET based IF-RT .
Patient selection and planning conditions
Forty-one patients who had been enrolled in the previously published ‘PET-PLAN pilot trial’ were selected for the presented in silico analysis . It was approved by the local ethics committee (Ärztekammer des Saarlandes), and all participating patients had given written informed consent. All selected patients had pathologically confirmed, medically inoperable stage II-III NSCLC and were candidates for definitive radio(chemo)therapy. Further inclusion criteria and the detailed modality of both the computed tomography and 18F-FDG-PET acquisition as prerequisites for the treatment planning are described in the preceding publication .
Target volumes and treatment planning (IMRT and 3D-CRT)
FDG-PET-based target volumes were defined by a radiation oncologist with the support of a radiation physicist and a nuclear medicine physician. The registration of CT and FDG-PET-image sets was scrutinized for anatomical plausibility and – if necessary - manually adjusted. The gross tumor volume (GTV) of the primary tumor was initially autocontoured with a contrast-oriented source-to- background (‘S/B’) algorithm for FDG-PET–based delineation of tumor volumes as described by Nestle et al.  and subsequently adjusted to the tumor borders as shown in the CT scan displayed in soft tissue and lung window. FDG-PET–negative atelectasis was excluded from the GTV. The GTV was expanded to the clinical target volume (CTV) by 3 mm in all dimensions; then, the CTV was expanded by another 7 mm to create the planning target volume (PTV) for the primary tumor.
As to the involved lymph nodes, only American Joint Committee of Cancer (AJCC)-lymph node stations containing FDG-PET-positive lymph nodes were contoured as CTV, referring to an anatomical contouring guide , whereas CT-positive but FDG-PET-negative lymph nodes were disregarded. The nodal CTV was then expanded by 7 mm to obtain the PTV for the involved lymph node stations. Finally, the definitive PTV was obtained from the unification of the PTVs of the primary tumor and the involved lymph node stations.
Dose constraints for normal tissues were defined as follows: for the whole lung (the GTV of primary tumor had been subtracted from the contoured lung volume) V20 ≤ 35 %, and mean lung dose (MLD) ≤ 20 Gy; the spinal cord 45 Gy as maximum; the esophagus V50 ≤ 50 % or mean dose ≤34 Gy and maximum dose ≤70 Gy; and the heart V45 ≤ 67 % and V60 ≤ 33 %.
For each patient, both an IMRT- and a 3D-CRT-treatment plan were generated and optimized by an experienced radiation physicist, using the PINNACLE3 treatment planning system (Version 8.2, Philips Medical Systems, Best, The Netherlands). IMRT “step-and shoot” plans (photon energy: 6 MV) were calculated using ‘direct machine parameter optimization’. The final dose distribution was calculated with a collapsed cone algorithm. Typically, six or seven coplanar beams with 70 segments were used; the beam angles were individually adapted. Inverse planning was started with a default prescription dose of 70 Gy and standardized objectives for the coverage of the planning target volume (PTV) as provided by the ICRU Report 83 were applied . To allow for potential dose escalation, pre-defined lower OAR constraints were used as default values for the initial 70 Gy plan: whole lung – V20 < 28 %; heart – V40 < 20 %; spinal cord – maximum dose of 40 Gy; esophagus – V50 < 50 % and maximum dose of 70 Gy. The total prescription dose was escalated or deescalated by 2 Gy starting with 70 Gy until all dosimetric premises with respect to PTV coverage and OARs were met. To compare IMRT with 3D-CRT plans at equal dose levels with respect to dose distribution in out-of-field lymph node stations and sparing of OARs, IMRT plans were additionaly normalized to the same maximum dose level (defined as total reference dose, Dref) as it was reached with 3D-CRT plans (adapted doses are labeled as IMRTnorm). The normalization was achieved by matching the number of fractions of the IMRT plans with those of the 3D-CRT plans (no de novo IMRT-planning was performed for dose normalization).
For 3D-CRT planning at least three beams had to be used and the number and arrangement of beams and their weights, the use of wedges as well as the photon energy (6 and 18 MV) were optimized by an experienced medical physicist.
A minimal total prescription dose of 60 Gy was aspired and depended on the dose constraints for normal tissues. The prescribed single dose was 2 Gy and the total prescription dose was incrementally escalated by 2 Gy as long as the dose constraints outlined above were not surpassed. The maximum doses for both the IMRT and the 3D-CRT plan were limited to 110 Gy.
Contouring and grouping of out-of-field lymph node stations
All hilar and mediastinal AJCC-lymph node stations were contoured based on the individual CT-anatomy according to the atlas provided by Chapet et al. . Specifically, two nodal ‘out-of-field’-volumes were defined for each patient. The first volume, LNall_el, comprised the entire elective hilar and mediastinal lymph node stations, i.e. all uninvolved levels outside the CTV. The second volume, LNadj_el, included only the elective ‘out-of-field’-lymph node levels 7 (infracarinal, if uninvolved), the ipsilateral hilum (levels 10/11, if uninvolved) and those uninvolved and directly adjacent to (involved) lymph node levels included in the CTV. For example, if only level 4R (deep paratracheal lymph nodes, right side) was involved, the volume Lnadj_el contained the bordering lymph node levels 2R, 4 L, 7 and 10_11R.
Dose distribution in out-of-field nodal regions
Dose volume histograms (DVH) were created for both delineated out-of-field lymph node volumes (LNall_el, LNadj_el) for three plans in each patient (3D-CRT, IMRT and IMRTnorm). For LNall_el and LNadj_el mean doses were calculated and additionally (in order to correct for heterogeneity in dose distribution), the equivalent uniform dose (EUD) was derived from DVH-based values by using the formula as given by Wu et al. :
where N is the number of voxels in the anatomic region of interest, D i the dose in the i th voxel, and a is a tumor-dependent value, which describes the dose-volume effect. The EUD-concept was designed to translate a heterogeneous dose distribution in a given volume into an isoeffective (hypothetical) homogeneous dose distribution with the compared biologic effect mostly being tumor control. One advantage of the EUD-concept is that it simplifies the comparison of volumes with various heterogeneous dose-distribution while it remains at the same time an objective measure.
Tumor control probability
TCP values were computed for LNall_el and LNadj_el in all three plan variants (3D-CRT, IMRT and IMRTnorm). For TCP calculation, the formula given by Martel et al. was applied :
where D50 is the dose needed to obtain a 50 % TCP, D is the actual total dose deposited in the irradiated volume, and γ is the normalized slope of the sigmoidal response curve at D50. Parameters to control microscopic disease in lung cancer were used as provided by Okunieff et al. . Precisely, D was derived from the computed EUD, D50 was 36.5 Gy and γ was 0.72.
The paired sample t-test was used to test for statistical differences between data sets in case of normal distribution, otherwise the Wilcoxon signed rank test was applied. Differences were considered significant when the p-value was <0.05. Origin Pro 9.0 (OriginLab Corporation, Northampton, MA) was used for the statistical analysis.
The patient characteristics are presented in Table 1. Primary tumors were located as follows: centrally (left: n = 10; right: n = 16), left upper lobe (n = 6), left lower lobe (n = 3), right upper lobe (n = 3) and right lower lobe (n = 3). Thirty-nine patients had lymph node involvment as diagnosed with FDG-PET. Six patients had a ‘single level’ involvement and 33 patients had a ‘multilevel’ spread. The mean number of involved lymph node stations was 2.5 ± 0.9 sd (standard deviation).
Comparison of IMRT and 3D-CRT planning
Relevant PTV-based planning parameters are shown in Table 2. With IMRT-planning a significangly higher mean total dose could be achieved as compared with 3D-CRT (74.3 vs. 70.1 Gy, p < 0.00001). Also, 7 of 41 patients could be administered a total dose of 66 Gy or more with IMRT. These results can be ascribed to the superior dose conformality of IMRT vs. 3D-CRT (CI: 0.79 vs. 0.50, p < 0.00001) resulting in an improved sparing of OAR, which is shown in Table 3. At normalized dose levels, IMRT was superior to 3D-CRT in sparing the spinal cord (Dmax) and limiting the exposure of the esophagus with higher doses (V60). Differences between IMRT and 3D-CRT in regard to PTV coverage and OAR-exposure are exemplified in one patient in Fig. 1a.
Dose distribution and TCP in out-of-field lymph node stations
Values referring to mean dose, EUD, and TCP in LNall_el and LNadj_el are presented in Table 4. If regarded at normalized dose levels (IMRTnorm vs. 3D-CRT), mean dose in LNall_el and LNadj_el was significantly decreased with IMRTnorm. EUD was significantly decreased with IMRT and IMRTnorm in comparison with 3D-CRT. The significantly lower doses in LNall_el and LNadj_el with IMRTnorm vs. 3D-CRT translated into significantly, but overall moderately reduced values in TCP. TCP-values did not differ significantly between IMRT and 3D-CRT if IMRT is used for additional dose-escalation.
In Fig. 1b the delineated out-of-field lymph node stations and the resulting DVH for both volumes LNall_el and LNadj_el are depicted for one patient.
The presented planning study revealed a potential hazard in the use of IMRT in IF-RT of locally advanced NSCLC: IMRT goes along with a decreased dose deposition in out-of-field lymph node stations and may thus involve the risk of an increased rate of out-of-field nodal failure. However, this holds true only under the condition that 3D-CRT and IMRT are compared at equalized dose levels. Furthermore, the calculated differences in TCP were slim between IMRTnorm and 3D-CRT (12.6 % vs. 14.8 %). As our results show: if the full potential of IMRT for dose escalation were exploited, no signficant difference would be observed in TCP-values compared to 3D-CRT, while EUD in out-of-field lymph node stations still remained lower. These findings indicate that one cannot count as much on IMRT as on 3D-CRT to compensate for inaccurate nodal staging by providing effective incidental dose-coverage in elective nodal regions. It is of note that only an increased rate of isolated nodal failure (as opposed to a combined in-field and out-of-field recurrence) should be considered as detrimental for IMRT. However, some meaningful strengths of IMRT in combination with IF-RT have to be weighed against its potential adverse effects on tumor control in elective nodal regions. As we could demonstrate, IMRT opens the door to a significant dose-escalation or alternatively embodies the benefit of substantial sparing of normal tissues at an equal total dose to that of 3D-CRT.
To our knowledge this is the first study examining the dose distribution in elective lymph node stations adjacent to the CTV in the context of IF-IMRT in NSCLC. Several planning studies examined the potential benefits of incidental nodal irradiation of uninvolved lymph node regions when using 3D-CRT [7, 8, 20, 21]. Even with PET-based IF-RT of stage III NSCLC, rates of out-of-field nodal recurrences can be as high as 10 % and therefore cannot be considered as insignificant . It is controversial if these recurrences can be attributed either to the small a priori risk of microscopic involvement in elective nodal regions or the benefits of incidental nodal irradiation. Yet, thorough dosimetric analyses of isolated out-of-field nodal failures indicate a dose–response relationship, and a critical cutoff may be in the range of 40–50 Gy [7, 21]. Depending on the total burden of subclinical disease a wide range of therapeutic total doses can be assumed . For a collective of patients with node-negative head and neck cancer a reduction of neck relapses by more than 90 % with total doses of 50 Gy and by less than 50 % with less than 30 Gy (2 Gy single fraction) was reported . In that context it has to be stressed that the calculation of the presented TCP values is inaccurate inasmuch as neither the ‘true’ load of subclinical disease nor the ‘true’ tumoricidal doses remain known for the average patient. Nevertheless, the use of somewhat arbitrary radiobiologic parameters for TCP-estimation still provides some guidance when comparing individual dose disbributions. Also, the TCP-values should be regarded as complementary to the by far more objective EUD-concept.
As the use of IMRT has been widely adopted for inoperable stage III NSCLC our results underline the importance of keeping track of out-of-field recurrences at times when highly conformal dose distributions and tight margins around the CTV are applied in the era of image-guided radiotherapy.
Incidental irradiation of out-of-field lymph node stations in NSCLC may be beneficial with respect to control of microscopic lymph node involvement. With IMRT, smaller doses are deposited in lymph node stations adjacent to the CTV due to its higher conformality as compared to 3D-conformal radiotherapy. The planning study presented here shows a moderate yet statistically significant reduction in hypothetical tumor control probability for microscopic involvement in out-of-field lymph node stations for IMRT vs. 3D-CRT at equal dose levels. This potential drawback may be outweighed by the superiority of IMRT in sparing organs at risk and by its potential to achieve higher total doses in the GTV. Yet, accurate – i.e. FDG-PET based – lymph node staging seems to be an important prerequisite for IMRT.
De Ruysscher D, Faivre-Finn C, Nestle U, Hurkmans CW, Le Péchoux C, Price A, et al. European Organisation for Research and Treatment of Cancer Recommendations for Planning and Delivery of High-Dose, High-Precision Radiotherapy for Lung Cancer. J Clin Oncol. 2010;28:5301–10.
Bradley J, Bae K, Choi N, Forster K, Siegel BA, Brunetti J, et al. A phase II comparative study of gross tumor volume definition with or without PET/CT fusion in dosimetric planning for non-small-cell lung cancer (NSCLC): primary analysis of Radiation Therapy Oncology Group (RTOG) 0515. Int J Radiat Oncol Biol Phys. 2012;82:435–41.
Rosenzweig KE, Sura S, Jackson A, Yorke E. Involved-field radiation therapy for inoperable non small-cell lung cancer. J Clin Oncol. 2007;35:5557–61.
Sulman EP, Komaki R, Klopp AH, Cox JD, Chang JY. Exclusion of selective nodal irradiation is associated with minimal elective nodal failure in non-small cell lung cancer. Radiat Oncol. 2009;4:5–11.
Belderbos JS, Kepka L, Kong FM, Martel MK, Videtic GM, Jeremic B. Report from the International Atomic Energy Agency (IAEA) consultants’ meeting on elective nodal irradiation in lung cancer: non-small cell lung cancer (NSCLC). Int J Radiat Oncol Biol Phys. 2008;72:335–42.
Fleckenstein J, Hellwig D, Kremp S, Grgic A, Gröschel A, Kirsch CM, et al. F-18-FDG-PET confined radiotherapy of locally advanced NSCLC with concomitant chemotherapy: results of the PET-PLAN pilot trial. Int J Radiat Oncol Biol Phys. 2011;81:e283–9.
Kimura T, Togami T, Nishiyama Y, Okawa M, Takashima H. Impact of incidental irradiation on clinically uninvolved nodal regions in patients with advanced non-small-cell lung cancer treated with involved-field radiation therapy: does incidental irradiation contribute to the low incidence of elective nodal failure? Int J Radiat Oncol Biol Phys. 2010;77:337–43.
Zhao L, Chen M, Ten Haken R, Chetty I, Chapet O, Hayman JA, et al. Three-dimensional conformal radiation may deliver considerable dose of incidental nodal irradiation in patients with early stage node-negative non-small cell lung cancer when the tumor is large and centrally located. Radiother Oncol. 2007;82:153–9.
Harris JP, Murphy JD, Hanlon Al, Le QT, Loo BW Jr, Diehn M. A population based comparative effectiveness study of radiotherapy techniques in stage III non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2014;88:872–84.
Shirvani SM, Jiang J, Gomez DR, Chang JY, Buchholz TA, Smith BD. Intensity modulated radiotherapy for stage III non-small cell lung cancer in the United States: predictors of use and association with toxicities. Lung Cancer. 2013;82:252–9.
Grills IS, Yan D, Martinez AA, Vicini FA, Wong JW, Kestin LL. Potential for reduced toxicity and dose escalation in the treatment of inoperable non-small-cell lung cancer: a comparison of intensity-modulated radiation therapy (IMRT), 3D conformal radiation, and elective nodal irradiation. Int J Radiat Oncol Biol Phys. 2003;57:875–90.
Murshed H, Liu HH, Liao Z, Barker JL, Wang X, Tucker SL, et al. Dose and volume reduction for normal lung using intensity-modulated radiotherapy for advanced-stage non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2004;58:1258–67.
Lievens Y, Nulens A, Gaber MA, Defraene G, De Wever W, Stroobants S, et al. Intensity-modulated radiotherapy for locally advanced non-small-cell lung cancer: a dose-escalation planning study. Int J Radiat Oncol Biol Phys. 2011;80:306–13.
Nestle U, Kremp S, Schaefer-Schuler A, Sebastian-Welsch C, Hellwig D, Rübe C, et al. Comparison of different methods for delineation of 18F-FDG PET-positive tissue for target volume definition in radiotherapy of patients with non-small cell lung cancer. J Nucl Med. 2005;46:1342–8.
Chapet O, Kong FM, Quint LE, Chang AC, Ten Haken RK, Eisbruch A, et al. CT-based definition of thoracic lymph node stations: an atlas from the University of Michigan. Int J Radiat Oncol Biol Phys. 2005;63:170–8.
ICRU. International Commission on Radiation Units and Measurements. Prescribing, recording, and reporting photon-beam intensity-modulated radiation therapy (IMRT). ICRU Report 83. J ICRU. 2010;10:1–106.
Wu Q, Mohan R, Niemierko A, Schmidt-Ullrich R. Optimization of intensity-modulated radiotherapy plans based on the equivalent uniform dose. Int J Radiat Oncol Biol Phys. 2002;52:224–35.
Martel M, Ten Haken R, Hazuka M, Kessler M, Strawderman M, Turrisi AT, et al. Estimation of tumor control probability model parameters from 3-D dose distributions of non-small cell lung cancer patients. Lung Cancer. 1999;24:31–7.
Okunieff P, Morgan D, Niemierko A, Suit HD. Radiation dose response of human tumors. Int J Radiat Oncol Biol Phys. 1995;32:1227–37.
Kepka L, Bujko K, Zolciak-Siwinska A, Garmol D. Incidental irradiation of mediastinal and hilar lymph node stations during 3D-conformal radiotherapy for non-small cell lung cancer. Acta Oncol. 2008;47:954–61.
Kepka L, Maciejewski B, Withers R. Does incidental irradiation with doses below 50 Gy effectively reduce isolated nodal failures in non-small-cell lung cancer: dose–response relationship. Int J Radiat Oncol Biol Phys. 2009;73:1391–6.
Rajpara RS, Schreibmann E, Fox T, Stapleford LJ, Beitler JJ, Curran WJ, et al. Locoregional tumor failure after definitive radiation for patients with stage III non-small cell lung cancer. Radiat Oncol. 2014;9:187.
Withers HR, Suwinski R. Radiation dose response for subclinical metastases. Semin Radiat Oncol. 1998;8:224–8.
Suwinski R, Maciejewski B, Withers HR. Dose–response relationship for elective neck irradiation of head and neck cancer – facts and controversies. Neoplasma. 1998;45:107–12.
We thank all the patients who participated in the study and Dr. Scott M Campbell (Nazareth College, Rochester, NY, U.S.A.) for his valuable editorial support.
The authors declare that they have no competing interests.
JF, SK, and CR participated in the design of the analysis and drafted the manuscript. SK, AE, and KK performed the target volume delineation, the plan-calculation and data acquisition; JF performed the statistical analysis. All authors read and approved the final manuscript.