Comparing different NTCP models that predict the incidence of radiation pneumonitis. Normal tissue complication probability

An evaluation of planning techniques for stereotactic body radiation therapy in lung tumors

PURPOSE

To evaluate four planning techniques for stereotactic body radiation therapy (SBRT) in lung tumors.

METHODS AND MATERIALS

Four SBRT plans were performed for 12 patients with stage I/II non-small-cell lung cancer under the following conditions: (1) conventional margins on free-breathing CT (plan 1), (2) generation of an internal target volume (ITV) using 4DCT with beam delivery under free-breathing conditions (plan 2), (3) gating at end-exhale (plan 3), and (4) gating at end-inhale (plan 4). Planning was performed following the RTOG 0236 protocol with a prescription dose of 54 Gy (3 fractions). For each plan 4D dose was calculated using deformable-image registration.

RESULTS

There was no significant difference in tumor dose delivered by the 4 plans. However, compared with plan 1, plans 2-4 reduced total lung BED by 1.9+/-1.2, 3.1+/-1.6 and 3.5+/-2.1 Gy, reduced mean lung dose by 0.8+/-0.5, 1.5+/-0.8, and 1.6+/-1.0 Gy, reduced V20 by 1.5+/-1.0%, 2.7+/-1.4%, and 2.8+/-1.8%, respectively, with p<0.01. Compared with plan 2, plans 3-4 reduced lung BED by 1.2+/-1.0 and 1.6+/-1.5 Gy, reduced mean lung dose by 0.6+/-0.5 and 0.8+/-0.7 Gy, reduced V20 by 1.2+/-1.1% and 1.3+/-1.5%, respectively, with p<0.01. The differences in lung BED, mean dose and V20 of plan 4 compared with plan 3 were insignificant.

CONCLUSIONS

Tumor dose coverage was statistically insignificant between all plans. However, compared with plan 1, plans 2-4 significantly reduced lung doses. Compared with plan 2, plan 3-4 also reduced lung toxicity. The difference in lung doses between plan 3 and plan 4 was not significant.

Intensity-modulated radiation therapy (IMRT) for inoperable non-small cell lung cancer: the Memorial Sloan-Kettering Cancer Center (MSKCC) experience

INTRODUCTION

Intensity-modulated radiation therapy (IMRT) is an advanced treatment delivery technique that can improve the therapeutic dose ratio. Its use in the treatment of inoperable non-small cell lung cancer (NSCLC) has not been well studied. This report reviews our experience with IMRT for patients with inoperable NSCLC.

METHODS AND MATERIALS

We performed a retrospective review of 55 patients with stage I-IIIB inoperable NSCLC treated with IMRT at our institution between 2001 and 2005. The study endpoints were toxicity, local control, and overall survival.

RESULTS

With a median follow-up of 26 months, the 2-year local control and overall survival rates for stage I/II patients were 50% and 55%, respectively. For the stage III patients, 2-year local control and overall survival rates were 58% and 58%, respectively, with a median survival time of 25 months. Six patients (11%) experienced grade 3 acute pulmonary toxicity. There were no acute treatment-related deaths. Two patients (4%) had grade 3 or worse late treatment-related pulmonary toxicity.

CONCLUSIONS

IMRT treatment resulted in promising outcomes for inoperable NSCLC patients.

Hybrid IMRT for treatment of cancers of the lung and esophagus

PURPOSE

To report on a hybrid intensity-modulated radiation therapy (IMRT; static plus IMRT beams treated concurrently) technique for lung and esophageal patients to reduce the volume of lung treated to low doses while delivering a conformal dose distribution.

METHODS

Treatment plans were analyzed for 18 patients (12 lung and 6 esophageal). Patients were treated with a hybrid technique that concurrently combines static (approximately two-thirds dose) and IMRT (approximately one-third dose) beams. These plans were compared with conventional three-dimensional (3D; non-IMRT) plans and all IMRT plans using custom four- and five-field arrangements and nine equally spaced coplanar beams. Plans were optimized to reduce V13 and V5 values. Dose-volume histograms were calculated for the planning target volume, heart, and the ipsilateral, contralateral, and total lung. Lung volumes V5, V13, V20, V30; mean lung dose (MLD); and the generalized equivalent uniform dose (gEUD) were calculated for each plan.

RESULTS

Hybrid plans treated significantly smaller total and contralateral lung volumes with low doses than nine-field IMRT plans. Largest reductions were for contralateral lung V5, V13, and V20 values for lung (-11%, -15%, -7%) and esophageal (-16%, -20%, -7%) patients. Smaller reductions were found also for 3D and four- and five-field IMRT plans. MLD and gEUDs were similar for all plan types. The 3D plans treated much larger extra planning target volumes to prescribed dose levels.

CONCLUSIONS

Hybrid IMRT demonstrated advantages for reduction of low-dose lung volumes in the thorax for reducing low dose to lung while also reducing the potential magnitude of dose deviations due to intrafraction motion and small field calculation accuracy.

Comparison of intensity-modulated radiotherapy planning based on manual and automatically generated contours using deformable image registration in four-dimensional computed tomography of lung cancer patients

PURPOSE

To evaluate the implications of differences between contours drawn manually and contours generated automatically by deformable image registration for four-dimensional (4D) treatment planning.

METHODS AND MATERIALS

In 12 lung cancer patients intensity-modulated radiotherapy (IMRT) planning was performed for both manual contours and automatically generated ("auto") contours in mid and peak expiration of 4D computed tomography scans, with the manual contours in peak inspiration serving as the reference for the displacement vector fields. Manual and auto plans were analyzed with respect to their coverage of the manual contours, which were assumed to represent the anatomically correct volumes.

RESULTS

Auto contours were on average larger than manual contours by up to 9%. Objective scores, D(2%) and D(98%) of the planning target volume, homogeneity and conformity indices, and coverage of normal tissue structures (lungs, heart, esophagus, spinal cord) at defined dose levels were not significantly different between plans (p = 0.22-0.94). Differences were statistically insignificant for the generalized equivalent uniform dose of the planning target volume (p = 0.19-0.94) and normal tissue complication probabilities for lung and esophagus (p = 0.13-0.47). Dosimetric differences >2% or >1 Gy were more frequent in patients with auto/manual volume differences > or =10% (p = 0.04).

CONCLUSIONS

The applied deformable image registration algorithm produces clinically plausible auto contours in the majority of structures. At this stage clinical supervision of the auto contouring process is required, and manual interventions may become necessary. Before routine use, further investigations are required, particularly to reduce imaging artifacts.

Four-dimensional Stereotactic Radiotherapy for Early Stage Non-Small Cell Lung Cancer: A Comparative Planning Study

In this study we sought to assess the potential of the respiratory tumor tracking system of the CyberKnife to administer 3 fractions of 15 Gy in the treatment of early stage non-small cell lung cancer (NSCLC). The CyberKnife plans were compared to those developed for 3-D conformal radiotherapy (3-D CRT) administering 20 fractions of 3 Gy based on a slow CT. Ten patients with stage I NSCLC, who were previously treated with 3-D CRT, were re-planned with the CyberKnife treatment planning system. In the 3-D CRT plan, the planning target volume (PTV) included the gross tumor volume (GTV)slow and a 15-mm margin, whereas in the CyberKnife plan the margin was 8 mm. The physical doses from both treatment plans were converted to normalized total doses (NTD) using the linear quadratic model with an α/βtumor of 10 Gy and α/βorgansatrisk(OAR) of 3 Gy. The average minimal and mean doses administered to the PTV with the CyberKnife and 3-D CRT were 93 and 115.8 Gy and 61 and 66 Gy, respectively (p<0.0001). The mean V20 of the CyberKnife and 3-D CRT plans were 8.2% and 6.8%, respectively (p=0.124). Both plans complied with the OAR constraints. In conclusion, 4-dimensional stereotactic radiotherapy can increase the minimal and mean biological dose with 51% and 75%, in comparison with 3-D CRT without significantly increasing the V20, respectively.

Predictive factors for radiation-induced pulmonary toxicity after three-dimensional conformal chemoradiation in locally advanced non-small-cell lung cancer

BACKGROUND AND PURPOSE

Radiation pneumonitis (RP) is a restricting complication of non-small-cell lung cancer irradiation. Three-dimensional conformal radiotherapy (3D-CRT) represents an advance because exposure of normal tissues is minimised. This study tries to identify prognostic factors associated with severe RP.

MATERIALS AND METHODS

Eighty patients with stage IIIA (20%) and IIIB (80%) NSCLC treated with cisplatin- based induction chemotherapy followed by concurrent chemotherapy and hyperfractionated 3D-CRT (median dose: 72.4 Gy, range: 54.1-85.9) were retrospectively evaluated. Acute and late RP were scored using RTOG glossary. Potential predictive factors evaluated included clinical, therapeutic and dosimetric factors. The lungs were defined as a whole organ. Univariate and multivariate analyses were performed.

RESULTS

Early and late RP grade>or=3 were observed in two patients (2%) and 10 patients (12%), respectively. Five patients (6%) died of pulmonary toxicity, 3 of whom had pre-existing chronic obstructive pulmonary disease (COPD). Median time to occurrence of late RP was 4.5 months (range: 3-8). Multivariate analysis showed that COPD (OR=10.1, p=0.01) and NTCPkwa>30% (OR=10.5, p=0.007) were independently associated with late grade>or=3 RP. Incidence of RP>or=3 grade for patients with COPD and/or NTCPkwa>30% was 25% vs. 4% for patients without COPD and NTCPkwa<30% (p=0.01). Risk of severe RP was higher for patients with COPD and/or NTCPkwa>30% (OR=7.3; CI 95%=1.4-37.3, p=0.016).

CONCLUSIONS

COPD and NTCP are predictive of severe RP. Careful medical evaluation and meticulous treatment planning are of paramount importance to decrease the incidence of severe RP.

Using patient data similarities to predict radiation pneumonitis via a self-organizing map

This work investigates the use of the self-organizing map (SOM) technique for predicting lung radiation pneumonitis (RP) risk. SOM is an effective method for projecting and visualizing high-dimensional data in a low-dimensional space (map). By projecting patients with similar data (dose and non-dose factors) onto the same region of the map, commonalities in their outcomes can be visualized and categorized. Once built, the SOM may be used to predict pneumonitis risk by identifying the region of the map that is most similar to a patient's characteristics. Two SOM models were developed from a database of 219 lung cancer patients treated with radiation therapy (34 clinically diagnosed with Grade 2+ pneumonitis). The models were: SOM(all) built from all dose and non-dose factors and, for comparison, SOM(dose) built from dose factors alone. Both models were tested using ten-fold cross validation and Receiver Operating Characteristics (ROC) analysis. Models SOM(all) and SOM(dose) yielded ten-fold cross-validated ROC areas of 0.73 (sensitivity/specificity = 71%/68%) and 0.67 (sensitivity/specificity = 63%/66%), respectively. The significant difference between the cross-validated ROC areas of these two models (p < 0.05) implies that non-dose features add important information toward predicting RP risk. Among the input features selected by model SOM(all), the two with highest impact for increasing RP risk were: (a) higher mean lung dose and (b) chemotherapy prior to radiation therapy. The SOM model developed here may not be extrapolated to treatment techniques outside that used in our database, such as several-field lung intensity modulated radiation therapy or gated radiation therapy.

A new formula for normal tissue complication probability (NTCP) as a function of equivalent uniform dose (EUD)

To facilitate the use of biological outcome modeling for treatment planning, an exponential function is introduced as a simpler equivalent to the Lyman formula for calculating normal tissue complication probability (NTCP). The single parameter of the exponential function is chosen to reproduce the Lyman calculation to within approximately 0.3%, and thus enable easy conversion of data contained in empirical fits of Lyman parameters for organs at risk (OARs). Organ parameters for the new formula are given in terms of Lyman model m and TD(50), and conversely m and TD(50) are expressed in terms of the parameters of the new equation. The role of the Lyman volume-effect parameter n is unchanged from its role in the Lyman model. For a non-homogeneously irradiated OAR, an equation relates d(ref), n, v(eff) and the Niemierko equivalent uniform dose (EUD), where d(ref) and v(eff) are the reference dose and effective fractional volume of the Kutcher-Burman reduction algorithm (i.e. the LKB model). It follows in the LKB model that uniform EUD irradiation of an OAR results in the same NTCP as the original non-homogeneous distribution. The NTCP equation is therefore represented as a function of EUD. The inverse equation expresses EUD as a function of NTCP and is used to generate a table of EUD versus normal tissue complication probability for the Emami-Burman parameter fits as well as for OAR parameter sets from more recent data.

Restricted Field IMRT Dramatically Enhances IMRT Planning for Mesothelioma

PURPOSE

To improve the target coverage and normal tissue sparing of intensity-modulated radiotherapy (IMRT) for mesothelioma after extrapleural pneumonectomy.

METHODS AND MATERIALS

Thirteen plans from patients previously treated with IMRT for mesothelioma were replanned using a restricted field technique. This technique was novel in two ways. It limited the entrance beams to 200 degrees around the target and three to four beams per case had their field apertures restricted down to the level of the heart or liver to further limit the contralateral lung dose. New constraints were added that included a mean lung dose of <9.5 Gy and volume receiving >or=5 Gy of <55%.

RESULTS

In all cases, the planning target volume coverage was excellent, with an average of 97% coverage of the planning target volume by the target dose. No change was seen in the target coverage with the new technique. The heart, kidneys, and esophagus were all kept under tolerance in all cases. The average mean lung dose, volume receiving >or=20 Gy, and volume receiving >or=5 Gy with the new technique was 6.6 Gy, 3.0%, and 50.8%, respectively, compared with 13.8 Gy, 15%, and 90% with the previous technique (p < 0.0001 for all three comparisons). The maximal value for any case in the cohort was 8.0 Gy, 7.3%, and 57.5% for the mean lung dose, volume receiving >or=20 Gy, and volume receiving >or=5 Gy, respectively.

CONCLUSION

Restricted field IMRT provides an improved method to deliver IMRT to a complex target after extrapleural pneumonectomy. An upcoming Phase I trial will provide validation of these results.

A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy

PURPOSE

Provide a simple research tool that may be used to calculate the NCTP or TCP of a particular treatment plan. Illustrate the implementation of the EUD-based NTCP and TCP models as a research tool.

METHODS AND MATERIALS

A high-level computing language was chosen to implement Niemierko's EUD-based NTCP and TCP mathematical models. The necessary treatment planning software requirements were clearly defined.

RESULTS

The computer code is presented and explained. Six simple examples were created to quickly troubleshoot the reader's code implementation. A table of model parameters based on the Emami data was generated.

A neural network model to predict lung radiation-induced pneumonitis

A feed-forward neural network was investigated to predict the occurrence of lung radiation-induced Grade 2+ pneumonitis. The database consisted of 235 patients with lung cancer treated using radiotherapy, of whom 34 were diagnosed with Grade 2+ pneumonitis at follow-up. The network was constructed using an algorithm that alternately grew and pruned it, starting from the smallest possible network, until a satisfactory solution was found. The weights and biases of the network were computed using the error back-propagation approach. Momentum and variable leaning techniques were used to speed convergence. Using the growing/pruning approach, the network selected features from 66 dose and 27 non-dose variables. During network training, the 235 patients were randomly split into ten groups of approximately equal size. Eight groups were used to train the network, one group was used for early stopping training to prevent overfitting, and the remaining group was used as a test to measure the generalization capability of the network (cross-validation). Using this methodology, each of the ten groups was considered, in turn, as the test group (ten-fold cross-validation). For the optimized network constructed with input features selected from dose and non-dose variables, the area under the receiver operating characteristics (ROC) curve for cross-validated testing was 0.76 (sensitivity: 0.68, specificity: 0.69). For the optimized network constructed with input features selected only from dose variables, the area under the ROC curve for cross-validation was 0.67 (sensitivity: 0.53, specificity: 0.69). The difference between these two areas was statistically significant (p = 0.020), indicating that the addition of non-dose features can significantly improve the generalization capability of the network. A network for prospective testing was constructed with input features selected from dose and non-dose variables (all data were used for training). The optimized network architecture consisted of six input nodes (features), four hidden nodes, and one output node. The six input features were: lung volume receiving > 16 Gy (V16), generalized equivalent uniform dose (gEUD) for the exponent a = 1 (mean lung dose), gEUD for the exponent a = 3.5, free expiratory volume in 1 s (FEV1), diffusion capacity of carbon monoxide (DLCO%), and whether or not the patient underwent chemotherapy prior to radiotherapy. The significance of each input feature was individually evaluated by omitting it during network training and gauging its impact by the consequent deterioration in cross-validated ROC area. With the exception of FEV1 and whether or not the patient underwent chemotherapy prior to radiotherapy, all input features were found to be individually significant (p < 0.05). The network for prospective testing is publicly available via internet access.

Sensitivity of NTCP parameter values against a change of dose calculation algorithm

Optimization of radiation treatment planning requires estimations of the normal tissue complication probability (NTCP). A number of models exist that estimate NTCP from a calculated dose distribution. Since different dose calculation algorithms use different approximations the dose distributions predicted for a given treatment will in general depend on the algorithm. The purpose of this work is to test whether the optimal NTCP parameter values change significantly when the dose calculation algorithm is changed. The treatment plans for 17 breast cancer patients have retrospectively been recalculated with a collapsed cone algorithm (CC) to compare the NTCP estimates for radiation pneumonitis with those obtained from the clinically used pencil beam algorithm (PB). For the PB calculations the NTCP parameters were taken from previously published values for three different models. For the CC calculations the parameters were fitted to give the same NTCP as for the PB calculations. This paper demonstrates that significant shifts of the NTCP parameter values are observed for three models, comparable in magnitude to the uncertainties of the published parameter values. Thus, it is important to quote the applied dose calculation algorithm when reporting estimates of NTCP parameters in order to ensure correct use of the models.

Toxicity associated to radiotherapy treatment in lung cancer patients

Radiation therapy in combination with other treatments, such as chemotherapy, increases loco-regional control and survival in patients with lung cancer. Nevertheless, the subsequent toxicity of this treatment occurs in up to 37% of the irradiated patients. Some factors related to the patient, including performance status, pulmonary function tests (FEV1, DCLO), tumour site, as well as treatment-related factors such as radiation dose, fractionation and addition of chemotherapy, can be related to the risk of pulmonary toxicity. With the advent of tridimensional conformal radiotherapy (3DCRT), dose-volume histograms can be generated to assess the dose received by the organs at risk. Volume dose (Vdose), mean lung dose (MLD) and normal tissue complication probability (NTCP) are the dosimetric parameters most frequently used. The possible relationship between these parameters and clinical and anatomical factors has to be considered. Steroid treatment should be started soon in case of pneumonitis to avoid the development of late pulmonary fibrosis. Finally, some pharmacological agents to prevent radiation-related pneumonitis are under investigation.

Predicting Risk of Radiation-Induced Lung Injury

Radiation-induced lung injury (RILI) is the most common, dose-limiting complication of thoracic radio- and radiochemotherapy. Unfortunately, predicting which patients will suffer from this complication is extremely difficult. Ideally, individual phenotype- and genotype-based risk profiles should be able to identify patients who are resistant to RILI and who could benefit from dose escalation in chemoradiotherapy. This could result in better local control and overall survival. We review the risk predictors that are currently in clinical use--dosimetric parameters of radiotherapy such as normal tissue complication probability, mean lung dose, V20 and V30--as well as biomarkers that might individualize risk profiles. These biomarkers comprise a variety of proinflammatory and profibrotic cytokines and molecules including transforming growth factor beta1 that are implicated in development and persistence of RILI. Dosimetric parameters of radiotherapy show a low negative predictive value of 60% to 80%. Depending on the studied molecule, negative predictive value of biomarkers is approximately 50%. The predictive power of biomarkers might be increased if they are coupled with radiogenomics, e.g., genotyping analysis of single nucleotide polymorphisms in transforming growth factor beta1, transforming growth factor beta1 pathway genes, and other cytokines. Genetic variability and the complexity of RILI and its underlying molecular mechanisms make identification of biological risk predictors challenging. Further investigations are needed to develop more effective risk predictors of RILI.

A nomogram to predict radiation pneumonitis, derived from a combined analysis of RTOG 9311 and institutional data

PURPOSE

To test the Washington University (WU) patient dataset, analysis of which suggested that superior-to-inferior tumor position, maximum dose, and D35 (minimum dose to the hottest 35% of the lung volume) were valuable to predict radiation pneumonitis (RP), against the patient database from Radiation Therapy Oncology Group (RTOG) trial 9311.

METHODS AND MATERIALS

The entire dataset consisted of 324 patients receiving definitive conformal radiotherapy for non-small-cell lung cancer (WU = 219, RTOG 9311 = 129). Clinical, dosimetric, and tumor location parameters were modeled to predict RP in the individual datasets and in a combined dataset. Association quality with RP was assessed using Spearman's rank correlation (r) for univariate analysis and multivariate analysis; comparison between subgroups was tested using the Wilcoxon rank sum test.

RESULTS

The WU model to predict RP performed poorly for the RTOG 9311 data. The most predictive model in the RTOG 9311 dataset was a single-parameter model, D15 (r = 0.28). Combining the datasets, the best derived model was a two-parameter model consisting of mean lung dose and superior-to-inferior gross tumor volume position (r = 0.303). An equation and nomogram to predict the probability of RP was derived using the combined patient population.

CONCLUSIONS

Statistical models derived from a large pool of candidate models resulted in well-tuned models for each subset (WU or RTOG 9311), which did not perform well when applied to the other dataset. However, when the data were combined, a model was generated that performed well on each data subset. The final model incorporates two effects: greater risk due to inferior lung irradiation, and greater risk for increasing normal lung mean dose. This formula and nomogram may aid clinicians during radiation treatment planning for lung cancer.

Self-consistent tumor control probability and normal tissue complication probability models based on generalized EUDa)

Traditional methods to compute the tumor control probability (TCP) or normal tissue complication probability (NTCP) typically require a heterogeneous radiation dose distribution to be converted into a simple uniform dose distribution with an equivalent biological effect. Several power-law type dose-volume-histogram reduction schemes, particularly Niemierko's generalized equivalent uniform dose model [Med. Phys. 26, 1000 (1999)], have been proposed to achieve this goal. In this study, we carefully examine the mathematical outcome of these schemes. We demonstrate that (1) for tumors, with each tumor cell independently responding to local radiation dose, a closed-form analytical solution for tumor survival fraction and TCP can be obtained; (2) for serial structured normal tissues, an exponential power-law form relating survival to functional sub-unit (FSU) radiation is required, and a closed-form analytical solution for the related NTCP is provided; (3) in the case of a parallel structured normal tissue, when NTCP is determined solely by the number of the surviving FSUs, a mathematical solution is available only when there is a non-zero threshold dose and/or a finite critical dose defining the radiotherapy response. Some discussion is offered for the partial irradiation effect on normal tissues in this category; (4) for normal tissues with alternative architectures, where the radiation response of FSU is inhomogeneous, there is no exact global mathematical solution for SF or NTCP within the available schemes. Finally, numerical fits of our models to some experimental data are also presented.

The impact of breathing motion versus heterogeneity effects in lung cancer treatment planning

The purpose of this study is to investigate the effects of tissue heterogeneity and breathing-induced motion/deformation on conformal treatment planning for pulmonary tumors and to compare the magnitude and the clinical importance of changes induced by these effects. Treatment planning scans were acquired at normal exhale/inhale breathing states for fifteen patients. The internal target volume (ITV) was defined as the union of exhale and inhale gross tumor volumes uniformly expanded by 5 mm. Anterior/posterior opposed beams (AP/PA) and three-dimensional (3D)-conformal plans were designed using the unit-density exhale ("static") dataset. These plans were further used to calculate (a) density-corrected ("heterogeneous") static dose and (b) heterogeneous cumulative dose, including breathing deformations. The DPM Monte Carlo code was used for dose computations. For larger than coin-sized tumors, relative to unit-density plans, tumor and lung doses increased in the heterogeneity-corrected plans. In comparing cumulative and static plans, larger normal tissue complication probability changes were observed for tumors with larger motion amplitudes and uncompensated breathing-induced hot/cold spots in lung. Accounting for tissue heterogeneity resulted in average increases of 9% and 7% in mean lung dose (MLD) for the 6 MV and 15 MV photon beams, respectively. Breathing-induced effects resulted in approximately 1% and 2% average decreases in MLD from the static value, for the 6 and 15 MV photon beams, respectively. The magnitude of these effects was not found to correlate with the treatment plan technique, i.e., AP/PA versus 3D-CRT. Given a properly designed ITV, tissue heterogeneity effects are likely to have a larger clinical significance on tumor and normal lung treatment evaluation metrics than four-dimensional respiratory-induced changes.

Physical models and simpler dosimetric descriptors of radiation late toxicity

Predicting radiation damage to specific organs is becoming ever more challenging with the use of intensity-modulated beams, nonuniform dose distributions, partial organ irradiation, and interpatient and even intraorgan variations in radiation sensitivity. Data-based physical models can be of use in summarizing complicated dose-volume data to help describe clinical outcomes and ultimately aid in the prediction of clinical toxicity. This article attempts to provide a brief overview of the use of normal tissue complication probability (NTCP) models and other simple dose/volume metrics to describe a few clinically significant complications (either frequent or serious) associated with radiation therapy of the head and neck, thorax, and abdominal-pelvic regions. Specifically, it reviews the application of these methods for late toxicities of the parotid, lung, heart, spinal cord, liver, and rectum. It focuses on organ-specific NTCP parameters as well as simple dosimetric descriptors that might be used to help treatment plan evaluation in clinical practice.

Is there a selection bias in radiotherapy dose-escalation protocols?

BACKGROUND

To investigate the existence of a selection bias using a virtual radiotherapy dose-escalation trial. In dose-escalation trials, normal tissue constraints generally remain constant while the tumor dose is increased. Since tumor dose and normal tissue constraints are competing demands, a point will be reached at which the tumor dose cannot be increased without exceeding normal tissue constraints.

METHODS AND MATERIALS

In 9 patients with non-small-cell lung cancer, the tumor dose was escalated from 66 Gy to 78 Gy in 4-Gy dose levels using intensity-modulated radiotherapy planning, while the limiting normal tissue dose constraints remained constant. Dosimetric, radiobiologic, and other planning parameters were compared at the 66-Gy dose level for patients eligible for all dose levels and for those eligible only for lower dose levels.

RESULTS

Seven of 9 patients were eligible for all dose levels (Group E). Two of 9 patients ("ineligible" or Group I) qualified only for lower total doses (95% confidence interval, 0.075-0.6, significant). In Group E, mean planning target volumes were smaller (132 vs. 404 cm(3), nonsignificant), monitor units per fraction were significantly lower (448 vs. 802, p = 0.0008), and the average composite score for plan quality was better than in Group I (0.012 vs. 0.068, nonsignificant). Average tumor-control probabilities were higher (0.33 vs. 0.23, nonsignificant), and normal tissue-complication probabilities were lower for Group E than for Group I.

CONCLUSIONS

Patients eligible for higher dose levels had significantly superior estimated outcome parameters. A method to eliminate this source of error in the interpretation of dose-escalation trials is suggested.

[Approaches to individual radiotherapy in breast cancer: individual risk estimation and individualized techniques]

INTRODUCTION

Radiotherapy comprises an integral part of the curative therapy of breast cancer by improving the locoregional control and survival when given on an individualized basis. Conformal radiotherapy and three-dimensional radiation treatment planning enhance the safety of radiotherapy by adjusting the irradiated volume to the shape of the target volume, and providing control of the radiation dose to the organs at risk (OARs).

PATIENTS AND METHODS

The methods introduced at the authors' institute in 2002 are demonstrated. The breast/chest wall and lymph node areas were irradiated provided that there was a minimum risk of local or locoregional relapse of 10%. CT-based 3D radiotherapy treatment planning and individual patient-positioning were applied, with thermoplastic mask-fixation in the second part of the study. The dose constraints of the OARs were given in accordance with the literature recommendations. In the first group of patients, individually shaped blocks, in the second group, multileaf collimator, and in the third group, with the aim of a more homogenous dose-distribution in the target volume, intensity-modulated beams were applied.

RESULTS

During the study, 737 breast cancer patients received conformal radiotherapy based on individual risk estimation. In 372 cases only local, while in 365 cases locoregional radiotherapy was delivered. The dose-homogeneity in the target volume was significantly improved in the second period of the study, when segments were superposed on the radiotherapy fields. The proportions of the target volumes irradiated with +/-10% of the planned dose in the breast/chest wall, axillary and supraclavicular lymph nodes and internal mammary lymph nodes varied between 90.5-94.2%, 84.1-93.8% and 86.7-91.6%, respectively, depending on the radiation technique used. The parameters indicating the dose to the ipsilateral lung or to the heart were significantly higher when locoregional radiotherapy was applied compared to that in case of local radiotherapy. Radiation dose to the ipsilateral lung and the heart was significantly reduced in the second part of the study when locoregional, but not when local radiotherapy was delivered. The introduction of individual immobilization by means of thermoplastic mask-fixation resulted in a relevant decrease in the uncertainty due to breathing motion and daily positioning errors, and also in a significant reduction of the dose to the contralateral breast.

CONCLUSIONS

Adjuvant radiotherapy should be based on individual risk-benefit features. The need of the introduction of special techniques may be decided after the dose-volume analysis of the conformal radiotherapy plan based on 3D radiation treatment planning.

A potential to reduce pulmonary toxicity: The use of perfusion SPECT with IMRT for functional lung avoidance in radiotherapy of non-small cell lung cancer

BACKGROUND AND PURPOSE

The study aimed to examine specific avoidance of functional lung (FL) defined by a single photon emission computerized tomography (SPECT) lung perfusion scan, using intensity modulated radiotherapy (IMRT) and three-dimensional conformal radiotherapy (3-DCRT) in patients with non-small cell lung cancer (NSCLC).

MATERIALS AND METHODS

Patients with NSCLC underwent planning computerized tomography (CT) and lung perfusion SPECT scan in the treatment position using fiducial markers to allow co-registration in the treatment planning system. Radiotherapy (RT) volumes were delineated on the CT scan. FL was defined using co-registered SPECT images. Two inverse coplanar RT plans were generated for each patient: 4-field 3-DCRT and 5-field step-and-shoot IMRT. 3-DCRT plans were created using automated AutoPlan optimisation software, and IMRT plans were generated employing Pinnacle(3) treatment planning system (Philips Radiation Oncology Systems). All plans were prescribed to 64 Gy in 32 fractions using data for the 6 MV beam from an Elekta linear accelerator. The objectives for both plans were to minimize the volume of FL irradiated to 20 Gy (fV(20)) and dose variation within the planning target volume (PTV). A spinal cord dose was constrained to 46 Gy. Volume of PTV receiving 90% of the prescribed dose (PTV(90)), fV(20), and functional mean lung dose (fMLD) were recorded. The PTV(90)/fV(20) ratio was used to account for variations in both measures, where a higher value represented a better plan.

RESULTS

Thirty-four RT plans of 17 patients with stage I-IIIB NSCLC suitable for radical RT were analysed. In 6 patients with stage I-II disease there was no improvement in PTV(90), fV(20), PTV/fV(20) ratio and fMLD using IMRT compared to 3-DCRT. In 11 patients with stage IIIA-B disease, the PTV was equally well covered with IMRT and 3-DCRT plans, with IMRT producing better PTV(90)/fV(20) ratio (mean ratio - 7.2 vs. 5.3, respectively, p=0.001) and reduced fMLD figures compared to 3-DCRT (mean value - 11.5 vs. 14.3 Gy, p=0.001). This was due to reduction in fV(20) while maintaining PTV coverage.

CONCLUSION

The use of IMRT compared to 3-DCRT improves the avoidance of FL defined by perfusion SPECT scan in selected patients with locally advanced NSCLC. If the dose to FL is shown to be the primary determinant of lung toxicity, IMRT would allow for effective dose escalation by specific avoidance of FL.

Intensity modulation with respiratory gating for radiotherapy of the pleural space

We wanted to describe a technique for the implementation of intensity-modulated radiotherapy (IMRT) with a real-time position monitor (RPM) respiratory gating system for the treatment of pleural space with intact lung. The technique is illustrated by a case of pediatric osteosarcoma, metastatic to the pleura of the right lung. The patient was simulated in the supine position where a breathing tracer and computed tomography (CT) scans synchronized at end expiration were acquired using the RPM system. The gated CT images were used to define target volumes and critical structures. Right pleural gated IMRT delivered at end expiration was prescribed to a dose of 44 Gy, with 55 Gy delivered to areas of higher risk via simultaneous integrated boost (SIB) technique. IMRT was necessary to avoid exceeding the tolerance of intact lung. Although very good coverage of the target volume was achieved with a shell-shaped dose distribution, dose over the targets was relatively inhomogeneous. Portions of target volumes necessarily intruded into the right lung, the liver, and right kidney, limiting the degree of normal tissue sparing that could be achieved. The radiation doses to critical structures were acceptable and well tolerated. With intact lung, delivering a relatively high dose to the pleura with acceptable doses to surrounding normal tissues using respiratory gated pleural IMRT is feasible. Treatment delivery during a limited part of the respiratory cycle allows for reduced CT target volume motion errors, with reduction in the portion of the planning margin that accounts for respiratory motion, and subsequent increase in the therapeutic ratio.

Normal Tissue Tolerance Dose Metrics for Radiation Therapy of Major Organs

Late organ toxicity from therapeutic radiation is a function of many confounding variables. The total dose delivered to the organ and the volumes of organ exposed to a given dose of radiation are 2 important variables that can be used to predict the risk of late toxicity. Three-dimensional radiation planning enables accurate calculation of the volume of tissue exposed to a given dose of radiation, graphically depicted as a dose-volume histogram. Dose metrics obtained from this 3-dimensional dataset can be used as a quantitative measure to predict late toxicity. This review summarizes the published clinical data on the risk of late toxicity as a function of quantitative dose metrics and attempts to offer suggested dose constraints for radiation treatment planning.

Reduction of normal lung irradiation in locally advanced non-small-cell lung cancer patients, using ventilation images for functional avoidance

PURPOSE

To investigate the ability of four-dimensional computed tomography (4D-CT)-derived ventilation images to identify regions of highly functional lung for avoidance in intensity-modulated radiotherapy (IMRT) planning in locally advanced non-small-cell lung cancer (NSCLC).

METHODS AND MATERIALS

The treatment-planning records from 21 patients with Stage III NSCLC were selected. Ventilation images were generated from the 4D-CT sets, and each was imported into the treatment-planning system. Ninetieth percentile functional volumes (PFV90), constituting the 10% of the lung volume where the highest ventilation occurs, were generated. Baseline IMRT plans were generated using the lung volume constraint on V20 (<35%), and two additional plans were generated using constraints on the PFV90 without a volume constraint. Dose-volume and dose-function histograms (DVH, DFH) were generated and used to evaluate the planning target volume coverage, lung volume, and functional parameters for comparison of the plans.

RESULTS

The mean dose to the PFV90 was reduced by 2.9 Gy, and the DFH at 5 Gy (F5) was reduced by 9.6% (SE = 2.03%). The F5, F10, V5, and V10 were all significantly reduced from the baseline values. We identified a favorable subset of patients for whom there was a further significant improvement in the mean lung dose.

CONCLUSIONS

Four-dimensional computed tomography-derived ventilation regions were successfully used as avoidance structures to reduce the DVH and DFH at 5 Gy in all cases. In a subset, there was also a reduction in the F10 and V10 without a change in the V20, suggesting that this technique could be safely used.

The risk of early and late lung sequelae after conformal radiotherapy in breast cancer patients

PURPOSE

To study the risks of early and late radiogenic lung damage in breast cancer patients after conformal radiotherapy.

METHODS AND MATERIALS

Radiogenic lung sequelae were assessed prospectively in 119 patients by means of clinical signs, radiologic abnormalities, and the mean density change (MDC) of the irradiated lung on CT.

RESULTS

Significant positive associations were detected between the development of lung abnormalities 3 months or 1 year after the radiotherapy and the age of the patient, the ipsilateral mean lung dose (MLD), the radiation dose to 25% of the ipsilateral lung (D(25%)) and the volume of the ipsilateral lung receiving 20 Gy (V(20 Gy)). The irradiation of the axillary and supraclavicular lymph nodes favored the development of pneumonitis but not that of fibrosis. No relation was found between the preradiotherapy plasma TGF-beta level and the presence of radiogenic lung damage. At both time points, MDC was strongly related to age. Significant positive associations were demonstrated between the risks of pneumonitis or fibrosis and the age of the patient, MLD, D(25%), and V(20 Gy). A synergistic effect of MLD, D(25%), and V(20 Gy) with age in patients older than 59 years is suggested.

CONCLUSION

Our analyses indicate that the risks of early and late radiogenic lung sequelae are strongly related to the age of the patient, the volume of the irradiated lung, and the dose to it.

Comparison of inverse-planned three-dimensional conformal radiotherapy and intensity-modulated radiotherapy for non–small-cell lung cancer

PURPOSE

Lungs are the major dose-limiting organ during radiotherapy (RT) for non-small-cell lung cancer owing to the development of pneumonitis. This study compared intensity-modulated RT (IMRT) with three-dimensional conformal RT (3D-CRT) in reducing the dose to the lungs.

METHODS

Ten patients with localized non-small-cell lung cancer underwent computed tomography (CT). The planning target volume (PTV) was defined and the organs at risk were outlined. An inverse-planning program, AutoPlan, was used to design the beam angle-optimized six-field noncoplanar 3D-CRT plans. Each 3D-CRT plan was compared with a series of five IMRT plans per patient. The IMRT plans were created using a commercial algorithm and consisted of a series of three, five, seven, and nine equidistant coplanar field arrangements and one six-field noncoplanar plan. The planning objectives were to minimize the lung dose while maintaining the dose to the PTV. The percentage of lung volume receiving >20 Gy (V20) and the percentage of the PTV covered by the 90% isodose (PTV90) were the primary endpoints. The PTV90/V20 ratio was used as the parameter accounting for both the reduction in lung volume treated and the PTV coverage.

RESULTS

All IMRT plans, except for the three-field coplanar plans, improved the PTV90/V20 ratio significantly compared with the optimized 3D-CRT plan. Nine coplanar IMRT beams were significantly better than five or seven coplanar IMRT beams, with an improved PTV90/V20 ratio.

CONCLUSION

The results of our study have shown that IMRT can reduce the dose to the lungs compared with 3D-CRT by improving the conformity of the plan.

Prediction of radiation-induced changes in the lung after stereotactic body radiation therapy of non–small-cell lung cancer

PURPOSE

To estimate the risk of radiation-induced changes in the lung before single-dose treatment (stereotactic body radiation therapy [SBRT]) of lung cancer, the quantitative dose-response and volume-response relations must be known.

METHODS AND MATERIALS

A total of 64 patients treated for non-small-cell lung cancer with single doses of 20-30 Gy were classified according to the occurrence or nonoccurrence of perifocal changes in the lung detected by CT. Patients without toxic events in the lung were required to have >or=6 months of follow-up. The mean dose (D(mean)) in the ipsilateral lung and the volume receiving >7 or 10 Gy (V7 and V10, respectively) were used to calculate the dose-response and volume-response curves. The predictive value of additional variables was also investigated.

RESULTS

Of the 64 patients, 83% exhibited the selected endpoint. The tolerance values at a 50% probability of toxic events were 1.2 +/- 0.7 Gy for the D(mean) and 5.8 +/- 3.0% and 3.1 +/- 2.0% for V7 and V10, respectively. A nonsignificant shift to higher doses was seen for the dose-response curve for the upper compared with the lower part of the lung.

CONCLUSION

The D(mean), V7, and V10 can be used to predict the risk of lung toxicity after SBRT treatment of non-small-cell lung cancer. Because of the lack of patients with low prescribed doses, however, the related uncertainty of this prediction is still relatively large. The D(mean), V7, and V10 are equally well suited. The additional investigated variables did not provide significant advantages. The lower part of the lung appears to be more radiosensitive than the upper.

Equivalent uniform dose concept evaluated by theoretical dose volume histograms for thoracic irradiation

BACKGROUND AND PURPOSE

The goal of our study was to quantify the limits of the EUD models for use in score functions in inverse planning software, and for clinical application.

MATERIALS AND METHODS

We focused on oesophagus cancer irradiation. Our evaluation was based on theoretical dose volume histograms (DVH), and we analyzed them using volumetric and linear quadratic EUD models, average and maximum dose concepts, the linear quadratic model and the differential area between each DVH.

RESULTS

We evaluated our models using theoretical and more complex DVHs for the above regions of interest. We studied three types of DVH for the target volume: the first followed the ICRU dose homogeneity recommendations; the second was built out of the first requirements and the same average dose was built in for all cases; the third was truncated by a small dose hole. We also built theoretical DVHs for the organs at risk, in order to evaluate the limits of, and the ways to use both EUD(1) and EUD/LQ models, comparing them to the traditional ways of scoring a treatment plan. For each volume of interest we built theoretical treatment plans with differences in the fractionation.

CONCLUSION

We concluded that both volumetric and linear quadratic EUDs should be used. Volumetric EUD(1) takes into account neither hot-cold spot compensation nor the differences in fractionation, but it is more sensitive to the increase of the irradiated volume. With linear quadratic EUD/LQ, a volumetric analysis of fractionation variation effort can be performed.

How extensive of a 4D dataset is needed to estimate cumulative dose distribution plan evaluation metrics in conformal lung therapy?

The purpose of this study was to investigate the number of intermediate states required to adequately approximate the clinically relevant cumulative dose to deforming/moving thoracic anatomy in four-dimensional (4D) conformal radiotherapy that uses 6 MV photons to target tumors. Four patients were involved in this study. For the first three patients, computed tomography images acquired at exhale and inhale were available; they were registered using B-spline deformation model and the computed transformation was further used to simulate intermediate states between exhale and inhale. For the fourth patient, 4D-acquired, phase-sorted datasets were available and each dataset was registered with the exhale dataset. The exhale-inhale transformation was also used to simulate intermediate states in order to compare the cumulative doses computed using the actual and the simulated datasets. Doses to each state were calculated using the Dose Planning Method (DPM) Monte Carlo code and dose was accumulated for scoring on the exhale anatomy via the transformation matrices for each state and time weighting factors. Cumulative doses were estimated using increasing numbers of intermediate states and compared to simpler scenarios such as a "2-state" model which used only the exhale and inhale datasets or the dose received during the average phase of the breathing cycle. Dose distributions for each modeled state as well as the cumulative doses were assessed using dose volume histograms and several treatment evaluation metrics such as mean lung dose, normal tissue complication probability, and generalized uniform dose. Although significant "point dose" differences can exist between each breathing state, the differences decrease when cumulative doses are considered, and can become less significant yet in terms of evaluation metrics depending upon the clinical end point. This study suggests that for certain "clinical" end points of importance for lung cancer, satisfactory predictions of accumulated total dose to be received by the distorting anatomy can be achieved by calculating the dose to but a few (or even simply the average) phases of the breathing cycle.

Incidence of radiation pneumonitis after thoracic irradiation: Dose-volume correlates

PURPOSE

To define clinical and dosimetric parameters correlated with the risk of clinically relevant radiation pneumonitis (RP) after thoracic radiotherapy.

METHODS AND MATERIALS

Records of consecutive patients treated with definitive thoracic radiotherapy were retrospectively reviewed for the incidence of RP of Grade 2 or greater by the Common Toxicity Criteria. Dose-volume histograms using total lung volume (TL) and TL minus gross tumor volume (TL-G) were created with and without heterogeneity corrections. Mean lung dose (MLD), effective lung volume (V(eff)), and percentage of TL or TL-G receiving greater than or equal to 10, 13, 15, 20, and 30 Gy (V10-V30, respectively) were analyzed by logistic regression. Receiver operating characteristic (ROC) curves were generated to estimate RP predictive values.

RESULTS

Twelve cases of RP were identified in 92 eligible patients. Mean lung dose, V10, V13, V15, V20, and V(eff) were significantly correlated to RP. Combinations of MLD, V(eff), V20, and V30 lost significance using TL-G and heterogeneity corrections. Receiver operating characteristic analysis determined V10 and V13 as the best predictors of RP risk, with a decrease in predictive value above those volumes.

CONCLUSIONS

Intrathoracic radiotherapy should be planned with caution when using radiotherapy techniques delivering doses of 10 to 15 Gy to large lung volumes.

Time trends in the maximal uptake of FDG on PET scan during thoracic radiotherapy. A prospective study in locally advanced non-small cell lung cancer (NSCLC) patients

BACKGROUND AND PURPOSE

18F-fluoro-2-deoxy-glucose (FDG) uptake on PET scan is a prognostic factor for outcome in NSCLC. We investigated changes in FDG uptake during fractionated radiotherapy in relation to metabolic response with the ultimate aim to adapt treatment according to early response.

METHODS AND MATERIALS

Twenty-three patients, medically inoperable or with advanced NSCLC, underwent four repeated PET-CT scans before, during and after radiotherapy. Changes in maximal standardized uptake value (SUVmax) were described. Patients were treated with accelerated radiotherapy with a total tumour-dose depending on normal tissue dose constraints.

RESULTS

The most striking result was the large intra-individual heterogeneity in the evolution of SUVmax. For the total group a non-significant increase in the first week (p=0.05), and a decrease in the second week (p=0.02) and after radiotherapy (p<0.01) was observed. Different time trends were shown for responders (no change during radiotherapy) and non-responders (48% increase during first week, p=0.02 and 15% decrease in the second week, p=0.04). Non-responders had a higher SUVmax on all time points investigated.

CONCLUSIONS

Time trends in SUVmax showed a large intra-individual heterogeneity and different patterns for metabolic responders and non-responders. These new findings may reflect intrinsic tumour characteristics and might finally be useful to adapt treatment.

Prospective assessment of dosimetric/physiologic-based models for predicting radiation pneumonitis

PURPOSE

Clinical and 3D dosimetric parameters are associated with symptomatic radiation pneumonitis rates in retrospective studies. Such parameters include: mean lung dose (MLD), radiation (RT) dose to perfused lung (via SPECT), and pre-RT lung function. Based on prior publications, we defined pre-RT criteria hypothesized to be predictive for later development of pneumonitis. We herein prospectively test the predictive abilities of these dosimetric/functional parameters on 2 cohorts of patients from Duke and The Netherlands Cancer Institute (NKI).

METHODS AND MATERIALS

For the Duke cohort, 55 eligible patients treated between 1999 and 2005 on a prospective IRB-approved study to monitor RT-induced lung injury were analyzed. A similar group of patients treated at the NKI between 1996 and 2002 were identified. Patients believed to be at high and low risk for pneumonitis were defined based on: (1) MLD; (2) OpRP (sum of predicted perfusion reduction based on regional dose-response curve); and (3) pre-RT DLCO. All doses reflected tissue density heterogeneity. The rates of grade > or =2 pneumonitis in the "presumed" high and low risk groups were compared using Fisher's exact test.

RESULTS

In the Duke group, pneumonitis rates in patients prospectively deemed to be at "high" vs. "low" risk are 7 of 20 and 9 of 35, respectively; p = 0.33 one-tailed Fisher's. Similarly, comparable rates for the NKI group are 4 of 21 and 6 of 44, respectively, p = 0.41 one-tailed Fisher's.

CONCLUSION

The prospective model appears unable to accurately segregate patients into high vs. low risk groups. However, considered retrospectively, these data are consistent with prior studies suggesting that dosimetric (e.g., MLD) and functional (e.g., PFTs or SPECT) parameters are predictive for RT-induced pneumonitis. Additional work is needed to better identify, and prospectively assess, predictors of RT-induced lung injury.

NTCP modelling and pulmonary function tests evaluation for the prediction of radiation induced pneumonitis in non-small-cell lung cancer radiotherapy

This work aims to evaluate the predictive strength of the relative seriality, parallel and Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) models regarding the incidence of radiation pneumonitis (RP), in a group of patients following lung cancer radiotherapy and also to examine their correlation with pulmonary function tests (PFTs). The study was based on 47 patients who received radiation therapy for stage III non-small-cell lung cancer. For each patient, lung dose volume histograms (DVHs) and the clinical treatment outcome were available. Clinical symptoms, radiological findings and pulmonary function tests incorporated in a post-treatment follow-up period of 18 months were used to assess the manifestation of radiation induced complications. Thirteen of the 47 patients were scored as having radiation induced pneumonitis, with RTOG criteria grade 3 and 28 of the 47 with RTOG criteria grade 2. Using this material, different methods of estimating the likelihood of radiation effects were evaluated, by analysing patient data based on their full dose distributions and associating the calculated complication rates with the clinical follow-up records. Lungs were evaluated as a paired organ as well as individual lungs. Of the NTCP models examined in the overall group considering the dose distribution in the ipsilateral lung, all models were able to predict radiation induced pneumonitis only in the case of grade 2 radiation pneumonitis score, with the LKB model giving the best results (chi2-test: probability of agreement between the observed and predicted results Pchi(chi2)=0.524 using the 0.05 significance level). The NTCP modelling considering lungs as a paired organ did not give statistically acceptable results. In the case of lung cancer radiotherapy, the application of different published radiobiological parameters alters the NTCP results, but not excessively as in the case of breast cancer radiotherapy. In this relatively small group of lung cancer patients, no positive statistical correlation could be established between the incidence of radiation pneumonitis as estimated by NTCP models and the pulmonary function test evaluation. However, the use of PFTs as markers or predictors for the incidence or severity of radiation induced pneumonitis must be investigated further.

Early clinical and radiological pulmonary complications following breast cancer radiation therapy: NTCP fit with four different models

OBJECTIVE

To fit four different NTCP (Normal Tissue Complication Probability) models to prospectively collected data on short-term pulmonary complications following breast cancer radiotherapy (RT).

MATERIALS/METHODS

Four hundred and seventy-five breast cancer patients, referred to the Radiotherapy Department at Stockholm Söder Hospital (1994-1998) for adjuvant post-operative RT were prospectively followed for pulmonary complications 1, 4 and 7 months after the completion of RT. Eighty-seven patients with complete dose-volume histogram (DVH) of the ipsilateral lung were selected for the present analysis. Mean dose to the ipsilateral lateral lung ranged from 2.5 to 18Gy (median 12Gy). Three different endpoints were considered: (1) clinical pneumonitis scored according to CTC-NCIC criteria: asymptomatic (grade 0) vs grade 1 and grade 2; (2) radiological changes assessed with diagnostic chest X-ray: no/slight radiological changes vs moderate/severe; (3) radiological changes assessed with CT: no/slight vs moderate/severe. Four NTCP models were used: the Lyman model with DVH reduced to the equivalent uniform dose (LEUD), the Logit model with DVH reduced to EUD, the Mean Lung Dose (MLD) model and the Relative Seriality (RS) model. The data fitting procedure was done using the maximum likelihood analysis. The analysis was done on the entire population (n=87) and on a subgroup of patients treated with loco-regional RT (n=44).

RESULTS

24/87 patients (28%) developed clinical pneumonitis; 28/81 patients (35%) had radiological side effects on chest X-rays and 11/75 patients (15%) showed radiological density changes on Computed Tomography (CT). The analysis showed that the risk of clinical pneumonitis was a smooth function of EUD (calculated from DVH using n=0.86+/-0.10, best fit result). With LEUD, the relationship between EUD and NTCP could be described with a D(50) of 16.4Gy+/-1.1Gy and a steepness parameter m of 0.36+/-0.7. The results found in the overall population were substantially confirmed in the subgroup of patients treated with loco-regional RT.

CONCLUSIONS

A large group of prospective patient data (87 pts), including grade 1 pneumonitis, were analysed. The four NTCP models fit quite accurately the considered endpoints. EUD or the mean lung dose are robust and simple parameters correlated with the risk of pneumonitis. For all endpoints the D(50) values ranged in an interval between 10 and 20Gy.

The atlas of complication incidence: a proposal for a new standard for reporting the results of radiotherapy protocols

We present a new method of reporting the results of radiotherapy protocols. The dose-volume atlas of complication incidence is a comprehensive and unbiased summary of the dose-volume exposures and complications occurring in patients after treatment. This new tool provides clear and systematic information about the safety of regions of dose-volume exposure previously treated that can be used when considering new treatments. Actuarial and model-dependent versions of the atlas are described. By using the raw data in the appropriate forms of the atlas, logistic regression, Kaplan-Meier, and Cox proportional hazards analysis can be performed, allowing for the independent calculation of dose-volume response. The data required are simple enough that provided compatible definitions of dose, volume, and complications are used, atlases from different protocols are potentially additive, facilitating the meta-analysis of inter-interinstitutional data. If this method were adopted as a standard for reporting the outcome of treatment protocols, a potentially synergistic increase in the utility of each protocol could result.

Complications aiguës et tardives des irradiations thoraciques

Radiotherapy of lung cancer is often complexe because of several organs at risk such as lung, heart, esophagus or spinal cord. An accurate balance needs to be defined between the necessities to reach a local control and to limit the risks of toxicities for each of them. Several parameters were significantly correlated with radiation induced lung toxicities (V13, V20, 30, Mean Dose) and esophagitis (V40 to V60). However no parameters are clearly defined for heart toxicities. The large number of parameters described for lung and esophagus highlights the necessity to perform an overall analysis of the DVHs that could be possible by using Nomal Tissue Complication Probability models.

The management of respiratory motion in radiation oncology report of AAPM Task Group 76

This document is the report of a task group of the AAPM and has been prepared primarily to advise medical physicists involved in the external-beam radiation therapy of patients with thoracic, abdominal, and pelvic tumors affected by respiratory motion. This report describes the magnitude of respiratory motion, discusses radiotherapy specific problems caused by respiratory motion, explains techniques that explicitly manage respiratory motion during radiotherapy and gives recommendations in the application of these techniques for patient care, including quality assurance (QA) guidelines for these devices and their use with conformal and intensity modulated radiotherapy. The technologies covered by this report are motion-encompassing methods, respiratory gated techniques, breath-hold techniques, forced shallow-breathing methods, and respiration-synchronized techniques. The main outcome of this report is a clinical process guide for managing respiratory motion. Included in this guide is the recommendation that tumor motion should be measured (when possible) for each patient for whom respiratory motion is a concern. If target motion is greater than 5 mm, a method of respiratory motion management is available, and if the patient can tolerate the procedure, respiratory motion management technology is appropriate. Respiratory motion management is also appropriate when the procedure will increase normal tissue sparing. Respiratory motion management involves further resources, education and the development of and adherence to QA procedures.

Comments on ‘Reconsidering the definition of a dose–volume histogram’—dose–mass histogram (DMH) versus dose–volume histogram (DVH) for predicting radiation-induced pneumonitis

In a recently published paper (Nioutsikou et al 2005 Phys. Med. Biol. 50 L17) the authors showed that the use of the dose-mass histogram (DMH) concept is a more accurate descriptor of the dose delivered to lung than the traditionally used dose-volume histogram (DVH) concept. Furthermore, they state that if a functional imaging modality could also be registered to the anatomical imaging modality providing a functional weighting across the organ (functional mass) then the more general and realistic concept of the dose-functioning mass histogram (D[F]MH) could be an even more appropriate descriptor. The comments of the present letter to the editor are in line with the basic arguments of that work since their general conclusions appear to be supported by the comparison of the DMH and DVH concepts using radiobiological measures. In this study, it is examined whether the dose-mass histogram (DMH) concept deviated significantly from the widely used dose-volume histogram (DVH) concept regarding the expected lung complications and if there are clinical indications supporting these results. The problem was investigated theoretically by applying two hypothetical dose distributions (Gaussian and semi-Gaussian shaped) on two lungs of uniform and varying densities. The influence of the deviation between DVHs and DMHs on the treatment outcome was estimated by using the relative seriality and LKB models using the Gagliardi et al (2000 Int. J. Radiat. Oncol. Biol. Phys. 46 373) and Seppenwoolde et al (2003 Int. J. Radiat. Oncol. Biol. Phys. 55 724) parameter sets for radiation pneumonitis, respectively. Furthermore, the biological equivalent of their difference was estimated by the biologically effective uniform dose (D) and equivalent uniform dose (EUD) concepts, respectively. It is shown that the relation between the DVHs and DMHs varies depending on the underlying cell density distribution and the applied dose distribution. However, the range of their deviation in terms of the expected clinical outcome was proven to be very large. Concluding, the effectiveness of the dose distribution delivered to the patients seems to be more closely related to the radiation effects when using the DMH concept.

Dose–volume modeling of the risk of postoperative pulmonary complications among esophageal cancer patients treated with concurrent chemoradiotherapy followed by surgery

PURPOSE

The aim of this study was to investigate the effect of radiation dose distribution in the lung on the risk of postoperative pulmonary complications among esophageal cancer patients.

METHODS AND MATERIALS

We analyzed data from 110 patients with esophageal cancer treated with concurrent chemoradiotherapy followed by surgery at our institution from 1998 to 2003. The endpoint for analysis was postsurgical pneumonia or acute respiratory distress syndrome. Dose-volume histograms (DVHs) and dose-mass histograms (DMHs) for the whole lung were used to fit normal-tissue complication probability (NTCP) models, and the quality of fits were compared using bootstrap analysis.

RESULTS

Normal-tissue complication probability modeling identified that the risk of postoperative pulmonary complications was most significantly associated with small absolute volumes of lung spared from doses > or = 5 Gy (VS5), that is, exposed to doses < 5 Gy. However, bootstrap analysis found no significant difference between the quality of this model and fits based on other dosimetric parameters, including mean lung dose, effective dose, and relative volume of lung receiving > or = 5 Gy, probably because of correlations among these factors. The choice of DVH vs. DMH or the use of fractionation correction did not significantly affect the results of the NTCP modeling. The parameter values estimated for the Lyman NTCP model were as follows (with 95% confidence intervals in parentheses): n = 1.85 (0.04, infinity), m = 0.55 (0.22, 1.02), and D50 = 17.5 Gy (9.4 Gy, 102 Gy).

CONCLUSIONS

In this cohort of esophageal cancer patients, several dosimetric parameters including mean lung dose, effective dose, and absolute volume of lung receiving < 5 Gy provided similar descriptions of the risk of postoperative pulmonary complications as a function of the radiation dose distribution in the lung.

Analysis of clinical and dosimetric factors associated with treatment-related pneumonitis (TRP) in patients with non-small-cell lung cancer (NSCLC) treated with concurrent chemotherapy and three-dimensional conformal radiotherapy (3D-CRT)

PURPOSE

To investigate factors associated with treatment-related pneumonitis in non-small-cell lung cancer patients treated with concurrent chemoradiotherapy.

PATIENTS AND METHODS

We retrospectively analyzed data from 223 patients treated with definitive concurrent chemoradiotherapy. Treatment-related pneumonitis was graded according to Common Terminology Criteria for Adverse Events version 3.0. Univariate and multivariate analyses were performed to identify predictive factors.

RESULTS

Median follow-up was 10.5 months (range, 1.4-58 months). The actuarial incidence of Grade > or =3 pneumonitis was 22% at 6 months and 32% at 1 year. By univariate analyses, lung volume, gross tumor volume, mean lung dose, and relative V5 through V65, in increments of 5 Gy, were all found to be significantly associated with treatment-related pneumonitis. The mean lung dose and rV5-rV65 were highly correlated (p < 0.0001). By multivariate analysis, relative V5 was the most significant factor associated with treatment-related pneumonitis; the 1-year actuarial incidences of Grade > or =3 pneumonitis in the group with V5 < or =42% and V5 >42% were 3% and 38%, respectively (p = 0.001).

CONCLUSIONS

In this study, a number of clinical and dosimetric factors were found to be significantly associated with treatment-related pneumonitis. However, rV5 was the only significant factor associated with this toxicity. Until it is better understood which dose range is most relevant, multiple clinical and dosimetric factors should be considered in treatment planning for non-small-cell lung cancer patients receiving concurrent chemoradiotherapy.

Clinical dose-volume histogram analysis in predicting radiation pneumonitis in Hodgkin's lymphoma

PURPOSE

To quantify the incidence of radiation pneumonitis (RP) in a modern Hodgkin's lymphoma (HL) cohort, and to identify any clinically relevant parameters that may influence the risk of RP.

METHODS AND MATERIALS

Between January 2003 and February 2005, 64 consecutive HL patients aged 18 years or older receiving radical mediastinal radiation therapy (RT) were retrospectively reviewed. Symptomatic cases of radiation pneumonitis were identified. Dose-volume histogram parameters, including V(13), V(20), V(30), and mean lung dose (MLD), were quantified.

RESULTS

At a median follow-up of 2.1 years, the actuarial survival for all patients was 91% at 3 years. There were 2 (2/64) cases of Radiation Therapy Oncology Group (RTOG) Grade 2 RP (incidence 3.1%). Both index cases with corresponding V(20) values of 47.0% and 40.7% were located in the upper quartile (2/16 cases), defined by a V(20) value of > or =36%, an incidence of 12.5% (p = 0.03). Similarly for total MLD, both index cases with values of 17.6 Gy and 16.4 Gy, respectively, were located in the upper quartile defined by MLD > or =14.2 Gy, an incidence of 11.8% (2/17 cases, p = 0.02).

CONCLUSIONS

Despite relatively high V(20) values in this study of HL patients, the incidence of RP was only 3%, lower compared with the lung cancer literature. We suggest the following clinically relevant parameters be considered in treatment plan assessment: a V(20) greater than 36% and an MLD greater than 14 Gy, over and above which the risk of RTOG Grade 2 or greater RP would be considered clinically significant.

Does transforming growth factor-beta1 predict for radiation-induced pneumonitis in patients treated for lung cancer?

The purpose of the study was to reassess the utility of transforming growth factor-beta-1 (TGF-beta1) together with dosimetric and tumor parameters as a predictor for radiation pneumonitis (RP). Of the 121 patients studied, 32 (26.4%) developed grade > or =1 RP, and 27 (22.3%) developed grade > or =2 RP. For the endpoint of grade > or =1 RP, those with V30>30% and an end-RT/baseline TGF-beta1 ratio> or =1 had a significantly higher incidence of RP than did those with V30>30% and an end-RT/baseline TGF-beta1 ratio<1. For most other patient groups, there were no clear associations between TGF-beta1 values and rates of RP. These findings suggest that TGF-beta1 is generally not predictive for RP except for the group of patients with a high V30.

Potential benefits of using non coplanar field and intensity modulated radiation therapy to preserve the heart in irradiation of lung tumors in the middle and lower lobes

PURPOSE

Investigate whether the use of non coplanar fields and intensity modulated radiation therapy (IMRT) reduces the dose to the heart, in irradiation of middle and lower lung tumors.

MATERIALS AND METHODS

Four plans are compared on 10 CT scans: (1) a reference plan, corresponding to coplanar plan of 3D conformal radiotherapy (3DCRT); (2) a 3DCRT(noncopl) plan, differing from reference plan by the change of one field in non coplanar fields; (3) an IMRT(copl) plan optimized on the same coplanar plan as reference plan; and (4) an IMRT(noncopl) plan optimized on the same non coplanar beams as 3DCRT(noncopl) plan. The equivalent uniform dose (EUD) for PTV was 74 Gy in 37 fractions.

RESULTS

In all plans, the 95% isodose surface covers at least 99% of the PTV with very similar conformity index values. A significant reduction in EUD, V30, V40 and V50 is observed for heart when either non coplanar fields or IMRT is used. IMRT also reduces the lung NTCP, V5, V13, V20 and V30 values and esophagus NTCP.

CONCLUSION

Both the use of non coplanar fields and IMRT dramatically reduces the dose received by the heart. The largest benefit is seen when the two techniques are combined.

A theoretical approach to the problem of dose-volume constraint estimation and their impact on the dose-volume histogram selection

This paper outlines a theoretical approach to the problem of estimating and choosing dose-volume constraints. Following this approach, a method of choosing dose-volume constraints based on biological criteria is proposed. This method is called "reverse normal tissue complication probability (NTCP) mapping into dose-volume space" and may be used as a general guidance to the problem of dose-volume constraint estimation. Dose-volume histograms (DVHs) are randomly simulated, and those resulting in clinically acceptable levels of complication, such as NTCP of 5 +/- 0.5%, are selected and averaged producing a mean DVH that is proven to result in the same level of NTCP. The points from the averaged DVH are proposed to serve as physical dose-volume constraints. The population-based critical volume and Lyman NTCP models with parameter sets taken from literature sources were used for the NTCP estimation. The impact of the prescribed value of the maximum dose to the organ, D(max), on the averaged DVH and the dose-volume constraint points is investigated. Constraint points for 16 organs are calculated. The impact of the number of constraints to be fulfilled based on the likelihood that a DVH satisfying them will result in an acceptable NTCP is also investigated. It is theoretically proven that the radiation treatment optimization based on physical objective functions can sufficiently well restrict the dose to the organs at risk, resulting in sufficiently low NTCP values through the employment of several appropriate dose-volume constraints. At the same time, the pure physical approach to optimization is self-restrictive due to the preassignment of acceptable NTCP levels thus excluding possible better solutions to the problem.

Reporting and analyzing statistical uncertainties in Monte Carlo-based treatment planning

PURPOSE

To investigate methods of reporting and analyzing statistical uncertainties in doses to targets and normal tissues in Monte Carlo (MC)-based treatment planning.

METHODS AND MATERIALS

Methods for quantifying statistical uncertainties in dose, such as uncertainty specification to specific dose points, or to volume-based regions, were analyzed in MC-based treatment planning for 5 lung cancer patients. The effect of statistical uncertainties on target and normal tissue dose indices was evaluated. The concept of uncertainty volume histograms for targets and organs at risk was examined, along with its utility, in conjunction with dose volume histograms, in assessing the acceptability of the statistical precision in dose distributions. The uncertainty evaluation tools were extended to four-dimensional planning for application on multiple instances of the patient geometry. All calculations were performed using the Dose Planning Method MC code.

RESULTS

For targets, generalized equivalent uniform doses and mean target doses converged at 150 million simulated histories, corresponding to relative uncertainties of less than 2% in the mean target doses. For the normal lung tissue (a volume-effect organ), mean lung dose and normal tissue complication probability converged at 150 million histories despite the large range in the relative organ uncertainty volume histograms. For "serial" normal tissues such as the spinal cord, large fluctuations exist in point dose relative uncertainties.

CONCLUSIONS

The tools presented here provide useful means for evaluating statistical precision in MC-based dose distributions. Tradeoffs between uncertainties in doses to targets, volume-effect organs, and "serial" normal tissues must be considered carefully in determining acceptable levels of statistical precision in MC-computed dose distributions.