Histogram reduction method for calculating complication probabilities for three-dimensional treatment planning evaluations

Radiobiological impact of reduced margins and treatment technique for prostate cancer in terms of tumor control probability (TCP) and normal tissue complication probability (NTCP)

Dose escalation in prostate radiotherapy is limited by normal tissue toxicities. The aim of this study was to assess the impact of margin size on tumor control and side effects for intensity-modulated radiation therapy (IMRT) and 3D conformal radiotherapy (3DCRT) treatment plans with increased dose. Eighteen patients with localized prostate cancer were enrolled. 3DCRT and IMRT plans were compared for a variety of margin sizes. A marker detectable on daily portal images was presupposed for narrow margins. Prescribed dose was 82 Gy within 41 fractions to the prostate clinical target volume (CTV). Tumor control probability (TCP) calculations based on the Poisson model including the linear quadratic approach were performed. Normal tissue complication probability (NTCP) was calculated for bladder, rectum and femoral heads according to the Lyman-Kutcher-Burman method. All plan types presented essentially identical TCP values and very low NTCP for bladder and femoral heads. Mean doses for these critical structures reached a minimum for IMRT with reduced margins. Two endpoints for rectal complications were analyzed. A marked decrease in NTCP for IMRT plans with narrow margins was seen for mild RTOG grade 2/3 as well as for proctitis/necrosis/stenosis/fistula, for which NTCP <7% was obtained. For equivalent TCP values, sparing of normal tissue was demonstrated with the narrow margin approach. The effect was more pronounced for IMRT than 3DCRT, with respect to NTCP for mild, as well as severe, rectal complications.

Focal dose escalation using FDG-PET-guided intensity-modulated radiation therapy boost for postoperative local recurrent rectal cancer: a planning study with comparison of DVH and NTCP

Background

To evaluate the safety of focal dose escalation to regions with standardized uptake value (SUV) >2.0 using intensity-modulated radiation therapy (IMRT) by comparison of radiotherapy plans using dose-volume histograms (DVHs) and normal tissue complication probability (NTCP) for postoperative local recurrent rectal cancer

Methods

First, we performed conventional radiotherapy with 40 Gy/20 fr. (CRT 40 Gy) for 12 patients with postoperative local recurrent rectal cancer, and then we performed FDG-PET/CT radiotherapy planning for those patients. We defined the regions with SUV > 2.0 as biological target volume (BTV) and made three boost plans for each patient: 1) CRT boost plan, 2) IMRT without dose-painting boost plan, and 3) IMRT with dose-painting boost plan. The total boost dose was 20 Gy. In IMRT with dose-painting boost plan, we increased the dose for BTV+5 mm by 30% of the prescribed dose. We added CRT boost plan to CRT 40 Gy (summed plan 1), IMRT without dose-painting boost plan to CRT 40 Gy (summed plan 2) and IMRT with dose-painting boost plan to CRT 40 Gy (summed plan 3), and we compared those plans using DVHs and NTCP.

Results

D_mean of PTV-PET and that of PTV-CT were 26.5 Gy and 21.3 Gy, respectively. V₅₀ of small bowel PRV in summed plan 1 was significantly higher than those in other plans ((summed plan 1 vs. summed plan 2 vs. summed plan 3: 47.11 ± 45.33 cm³ vs. 40.63 ± 39.13 cm³ vs. 41.25 ± 39.96 cm³(p < 0.01, respectively)). There were no significant differences in V₃₀, V₄₀, V₆₀, D_mean or NTCP of small bowel PRV.

Conclusions

FDG-PET-guided IMRT can facilitate focal dose-escalation to regions with SUV above 2.0 for postoperative local recurrent rectal cancer.

Slice-based plan evaluation methods for three dimensional conformal radiotherapy treatment planning

Dose volume histograms (DVHs) play a vital role in determining the optimal plan for radiotherapy treatment delivery. The current concepts of conformality index (CI), equivalent uniform dose (EUD) derived from dose volume histogram (DVH) does not provide any spatial information. In this study, slice-based evaluation methods have been proposed for spatially analyzing the radiotherapy treatment plans. A case of prostate cancer has been selected for demonstrating the proposed tools for evaluating the dose distribution. Three dimensional conformal radiotherapy treatment planning (3D-CRT) was performed to a dose of 27 Gy/15# with three fields (6 MV anteroposterior and two 15 MV lateral fields) employing multileaf collimator after delivering 45 Gy/25#. The dose was normalized to isocenter and the treatment plan was evaluated with DVH. The dose maximum point, conformality index, planning target volume coverage index (PCI), planning target volume overdose index (POI) and equivalent uniform dose (EUD) were evaluated for every single slice along the cranio-caudal direction for all the planning target volume (PTV) contours and plotted against the slice location. The dose maximum point plotted against the slice position helps in identifying the slices where the dose maximum point is outside the target volume. The plot of conformality index gives the information about the location of those slices where excess of surrounding normal tissues is encompassed inside the prescription isodose. POI quantifies the high dose regions inside the PTV slices that receive doses above 107% of the prescription dose. Similarly, the plot of PCI and EUD with slice position gives the information about those slices where the tumor is not covered adequately. The proposed methods in this study forms as a simpler way to assess the spatial distribution of the dose inside the target volume. It could be used in combination with the current plan evaluation tools and will be very helpful in analyzing the treatment plans.

Patient-specific internal radionuclide dosimetry

The development of patient-specific treatment planning systems is of outmost importance in the development of radionuclide dosimetry, taking into account that quantitative three-dimensional nuclear medical imaging can be used in this regard. At present, the established method for dosimetry is based on the measurement of the biokinetics by serial gamma-camera scans, followed by calculations of the administered activity and the residence times, resulting in the radiation-absorbed doses of critical organs. However, the quantification of the activity in different organs from planar data is hampered by inaccurate attenuation and scatter correction as well as because of background and organ overlay. In contrast, dosimetry based on quantitative three-dimensional data can be more accurate and allows an individualized approach, provided that all effects that degrade the quantitative content of the images have been corrected for. In addition, inhomogeneous organ accumulation of the radionuclide can be detected and possibly taken into account. The aim of this work is to provide adequate information on internal emitter dosimetry and a state-of-the-art review of the current methodology and future trends.

Quantitative Analyses of Normal Tissue Effects in the Clinic (QUANTEC): an introduction to the scientific issues

Advances in dose-volume/outcome (or normal tissue complication probability, NTCP) modeling since the seminal Emami paper from 1991 are reviewed. There has been some progress with an increasing number of studies on large patient samples with three-dimensional dosimetry. Nevertheless, NTCP models are not ideal. Issues related to the grading of side effects, selection of appropriate statistical methods, testing of internal and external model validity, and quantification of predictive power and statistical uncertainty, all limit the usefulness of much of the published literature. Synthesis (meta-analysis) of data from multiple studies is often impossible because of suboptimal primary analysis, insufficient reporting and variations in the models and predictors analyzed. Clinical limitations to the current knowledge base include the need for more data on the effect of patient-related cofactors, interactions between dose distribution and cytotoxic or molecular targeted agents, and the effect of dose fractions and overall treatment time in relation to nonuniform dose distributions. Research priorities for the next 5-10 years are proposed.

Equivalence in dose fall-off for isocentric and nonisocentric intracranial treatment modalities and its impact on dose fractionation schemes

PURPOSE

To investigate whether dose fall-off characteristics would be significantly different among intracranial radiosurgery modalities and the influence of these characteristics on fractionation schemes in terms of normal tissue sparing.

METHODS AND MATERIALS

An analytic model was developed to measure dose fall-off characteristics near the target independent of treatment modalities. Variations in the peripheral dose fall-off characteristics were then examined and compared for intracranial tumors treated with Gamma Knife, Cyberknife, or Novalis LINAC-based system. Equivalent uniform biologic effective dose (EUBED) for the normal brain tissue was calculated. Functional dependence of the normal brain EUBED on varying numbers of fractions (1 to 30) was studied for the three modalities.

RESULTS

The derived model fitted remarkably well for all the cases (R(2) > 0.99). No statistically significant differences in the dose fall-off relationships were found between the three modalities. Based on the extent of variations in the dose fall-off curves, normal brain EUBED was found to decrease with increasing number of fractions for the targets, with alpha/beta ranging from 10 to 20. This decrease was most pronounced for hypofractionated treatments with fewer than 10 fractions. Additionally, EUBED was found to increase slightly with increasing number of fractions for targets with alpha/beta ranging from 2 to 5.

CONCLUSION

Nearly identical dose fall-off characteristics were found for the Gamma Knife, Cyberknife, and Novalis systems. Based on EUBED calculations, normal brain sparing was found to favor hypofractionated treatments for fast-growing tumors with alpha/beta ranging from 10 to 20 and single fraction treatment for abnormal tissues with low alpha/beta values such as alpha/beta = 2.

Use of normal tissue complication probability models in the clinic

The Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) review summarizes the currently available three-dimensional dose/volume/outcome data to update and refine the normal tissue dose/volume tolerance guidelines provided by the classic Emami et al. paper published in 1991. A "clinician's view" on using the QUANTEC information in a responsible manner is presented along with a description of the most commonly used normal tissue complication probability (NTCP) models. A summary of organ-specific dose/volume/outcome data, based on the QUANTEC reviews, is included.

The effect of interfraction prostate motion on IMRT plans: a dose-volume histogram analysis using a Gaussian error function model

The Gaussian error function model, containing pairs of error and complementary error functions, was used to carry out cumulative dose-volume histogram (cDVH) analysis on prostate intensity modulated radiation therapy (IMRT) plans with interfraction prostate motion. Cumulative DVHs for clinical target volumes (CTVs) shifted in the anterior-posterior directions based on a 7-beam IMRT plan were calculated and modeled using the Pinnacle3 treatment planning system and a Gaussian error function, respectively. As the parameters in the error function model, namely, a, b and c were related to the shape of the cDVH curve, evaluation of cDVHs corresponding to the prostate motion based on the model parameters becomes possible as demonstrated in this study. It was found that deviations of the cDVH for the CTV were significant, when the CTV-planning target volume (PTV) margin was underestimated in the anterior-posterior directions, particularly in the posterior direction for a patient with relatively small prostate volume (39 cm3). Analysis of the cDVH for the CTV shifting in the anterior-posterior directions using the error function model showed that parameters a1,2, which were related to the maximum relative volume of the cDVH, changed symmetrically when the prostate was shifted in the anterior and posterior directions. This change was more significant for the larger prostate. For parameters b related to the slope of the cDVH, b1,2 changed symmetrically from the isocenter, when the CTV was within the PTV. This was different from parameters c (c1,2 are related to the maximum dose of the cDVH), which did not vary significantly with the prostate motion in the anterior-posterior directions and prostate volume. Using the patient data, this analysis validates the error function model, and further verified the clinical application of this mathematical model on treatment plan evaluations.

A comparison of HDR brachytherapy and IMRT techniques for dose escalation in prostate cancer: a radiobiological modeling study

A course of one to three large fractions of high dose rate (HDR) interstitial brachytherapy is an attractive alternative to intensity modulated radiation therapy (IMRT) for delivering boost doses to the prostate in combination with additional external beam irradiation for intermediate risk disease. The purpose of this work is to quantitatively compare single-fraction HDR boosts to biologically equivalent fractionated IMRT boosts, assuming idealized image guided delivery (igIMRT) and conventional delivery (cIMRT). For nine prostate patients, both seven-field IMRT and HDR boosts were planned. The linear-quadratic model was used to compute biologically equivalent dose prescriptions. The cIMRT plan was evaluated as a static plan and with simulated random and setup errors. The authors conclude that HDR delivery produces a therapeutic ratio which is significantly better than the conventional IMRT and comparable to or better than the igIMRT delivery. For the HDR, the rectal gBEUD analysis is strongly influenced by high dose DVH tails. A saturation BED, beyond which no further injury can occur, must be assumed. Modeling of organ motion uncertainties yields mean outcomes similar to static plan outcomes.

Modeling of alpha/beta for late rectal toxicity from a randomized phase II study: conventional versus hypofractionated scheme for localized prostate cancer

BACKGROUND

Recently, the use of hypo-fractionated treatment schemes for the prostate cancer has been encouraged due to the fact that alpha/beta ratio for prostate cancer should be low. However a major concern on the use of hypofractionation is the late rectal toxicity, it is important to be able to predict the risk of toxicity for alternative treatment schemes, with the best accuracy. The main purpose of this study is to evaluate the response of rectum wall to changes in fractionation and to quantify the alpha/beta ratio for late rectal toxicity

METHODS

162 patients with localized prostate cancer, treated with conformal radiotherapy, were enrolled in a phase II randomized trial. The patients were randomly assigned to 80 Gy in 40 fractions over 8 weeks (arm A) or 62 Gy in 20 fractions over 5 weeks (arm B). The median follow-up was 30 months. The late rectal toxicity was evaluated using the Radiation Therapy Oncology Group (RTOG) scale. It was assumed >or= Grade 2 (G2) toxicity incidence as primary end point. Fit of toxicity incidence by the Lyman-Burman-Kutcher (LKB) model was performed.

RESULTS

The crude incidence of late rectal toxicity >or= G2 was 14% and 12% for the standard arm and the hypofractionated arm, respectively. The crude incidence of late rectal toxicity >or= G2 was 14.0% and 12.3% for the arm A and B, respectively. For the arm A, volumes receiving >or= 50 Gy (V50) and 70 Gy (V70) were 38.3 +/- 7.5% and 23.4 +/- 5.5%; for arm B, V38 and V54 were 40.9 +/- 6.8% and 24.5 +/- 4.4%. An alpha/beta ratio for late rectal toxicity very close to 3 Gy was found.

CONCLUSION

The >or= G2 late toxicities in both arms were comparable, indicating the feasibility of hypofractionated regimes in prostate cancer. An alpha/beta ratio for late rectal toxicity very close to 3 Gy was found.

SORS: a new software for the simulation of radiotherapy schedule

We present a software for choosing the best radiotherapy treatment schedule for head and neck cancers as a beginning radiotherapy plan or a temporarily interrupted plan. Its application occurs according to two modalities: the first adopts the best estimates for model parameters; the second takes into account the parameters' uncertainty too. In both cases, the choice becomes the schedule with the highest uncomplicated tumor control probability (UTCP). In the UTCP valuation, the normal tissue complication probability (NTCP) of each organ is related to the gravity of its possible late injury. For NTCP calculation, it has been adopted the empirical LKB (Lyman-Kutcher-Burman) model corrected for dose/fraction via linear-quadratic model and the incomplete repair effect. The tumor control probability (TCP) model is Poisson based and contains corrections for dose/fraction and regrowth effect; optionally, it can be accounted for the incomplete repair effect as well. At the end of processing, a detailed file with all informations about UTCP, TCP and single organ NTCP is furnished for every examined schedule. Moreover, a useful 3-D graphic representation of the schedule's UTCP is available, allowing the physician to easily understand the schedules with the highest radiotherapeutic efficacy. The open source characteristic allows the program to adapt to the individual clinical case as well as to be a valid support in radiobiological research.

Expected clinical impact of the differences between planned and delivered dose distributions in helical tomotherapy for treating head and neck cancer using helical megavoltage CT images

Helical Tomotherapy (HT) has become increasingly popular over the past few years. However, its clinical efficacy and effectiveness continues to be investigated. Pre-treatment patient repositioning in highly conformal image-guided radiation therapy modalities is a prerequisite for reducing setup uncertainties. A MVCT image set has to be acquired to account for daily changes in the patient's internal anatomy and setup position. Furthermore, a comparison should be performed to the kVCT study used for dosimetric planning, by a registration process which results in repositioning the patient according to specific transitional and rotational shifts. Different image registration techniques may lead to different repositioning of the patient and, as a result, to varying delivered doses. This study aims to investigate the expected effect of patient setup correction using the Hi-Art tomotherapy system by employing radiobiological measures such as the biologically effective uniform dose (BEUD) and the complication-free tumor control probability (P+). In this study, a typical case of lung cancer with metastatic head & neck disease was investigated by developing a Helical Tomotherapy plan. For the Tomotherapy HiArt plan, the dedicated Tomotherapy treatment planning station was used. Three dose distributions (planned and delivered with and without patient setup correction) were compared based on radiobiological measures by using the P+ index and the BEUD concept as the common prescription point of the plans and plotting the tissue response probabilities against the mean target dose for a range of prescription doses. The applied plan evaluation method shows that in this cancer case the planned and delivered dose distributions with and without patient setup correction give a P+ of 81.6%, 80.9% and 72.2%, for a BEUD to the planning target volume (PTV) of 78.0Gy, 77.7Gy and 75.4Gy, respectively. The corresponding tumor control probabilities are 86.3%, 85.1% and 75.1%, whereas the total complication probabilities are 4.64%, 4.20% and 2.89%, respectively. HT can encompass the often large PTV required while minimizing the volume of the organs at risk receiving high dose. However, the effectiveness of a HT treatment plan can be considerably deteriorated if an accurate patient setup system is not available. Taking into account the dose-response relations of the irradiated tumors and normal tissues, a radiobiological treatment plan evaluation can be performed, which may provide a closer association of the delivered treatment with the clinical outcome. In such situations, for effective evaluation and comparison of different treatment plans, traditional dose based evaluation tools can be complemented by the use of P+,BEUD diagrams.

Analysis of outcomes in radiation oncology: an integrated computational platform

Radiotherapy research and outcome analyses are essential for evaluating new methods of radiation delivery and for assessing the benefits of a given technology on locoregional control and overall survival. In this article, a computational platform is presented to facilitate radiotherapy research and outcome studies in radiation oncology. This computational platform consists of (1) an infrastructural database that stores patient diagnosis, IMRT treatment details, and follow-up information, (2) an interface tool that is used to import and export IMRT plans in DICOM RT and AAPM/RTOG formats from a wide range of planning systems to facilitate reproducible research, (3) a graphical data analysis and programming tool that visualizes all aspects of an IMRT plan including dose, contour, and image data to aid the analysis of treatment plans, and (4) a software package that calculates radiobiological models to evaluate IMRT treatment plans. Given the limited number of general-purpose computational environments for radiotherapy research and outcome studies, this computational platform represents a powerful and convenient tool that is well suited for analyzing dose distributions biologically and correlating them with the delivered radiation dose distributions and other patient-related clinical factors. In addition the database is web-based and accessible by multiple users, facilitating its convenient application and use.

Accrediting radiation technique in a multicentre trial of chemoradiation for pancreatic cancer

Before a multicentre trial of 3-D conformal radiotherapy to treat cancer of the pancreas, participating clinicians were asked to complete an accreditation exercise. This involved planning two test cases according to the study protocol, then returning hard copies of the plans and dosimetric data for review. Any radiation technique that achieved the specified constraints was allowed. Eighteen treatment plans were assessed. Seven plans were prescribed incorrect doses and two of the planning target volumes did not comply with protocol guidelines. All plans met predefined normal tissue dose constraints. The identified errors were attributable to unforeseen ambiguities in protocol documentation. They were addressed by feedback and corresponding amendments to protocol documentation. Summary radiobiological measures including total weighted normal tissue equivalent uniform dose varied significantly between centres. This accreditation exercise successfully identified significant potential sources of protocol violations, which were then easily corrected. We believe that this process should be applied to all clinical trials involving radiotherapy. Due to the limitations of data analysis with hard-copy information only, it is recommended that complete planning datasets from treatment-planning systems be collected through a digital submission process.

Analysis of salivary flow and dose-volume modeling of complication incidence in patients with head-and-neck cancer receiving intensity-modulated radiotherapy

PURPOSE

To investigate dose-volume effects of salivary flow and the functional recovery over time, using salivary function data and different models of normal tissue complication probability (NTCP).

METHODS AND MATERIALS

A total of 59 patients with head-and-neck cancer treated with intensity-modulated radiotherapy (IMRT) were analyzed in the present study. The toxicity was evaluated using the Radiation Therapy Oncology Group (RTOG) scale and salivary flows, both unstimulated (USF) and stimulated (SSF). The assessments were done before radiotherapy (RT) and at 3, 6, 12, 18, and 24 months after RT. Grade 3 toxicity was the primary endpoint. Analyses of toxicity incidence at 3, 6, and 12 months after RT were performed by both the Lyman-Kutcher-Burman (LKB) and relative seriality (S) models.

RESULTS

A significant correlation was found between the incidence of Grade 3 toxicity and the incidence of patients with a reduction in SSF to <25% of the pre-RT value. Better correlations resulted between the RTOG toxicity score and the dosimetric parameters, compared with USF/SSF. The TD(50), assessed by the LKB model, was 21.4, 27.8, and 41.6 Gy at 3, 6, and 12 months after RT, respectively. The TD(50), assessed by the S model, was 20.0, 26.3, and 40.0 Gy at 3, 6, and 12 months after RT, respectively.

CONCLUSION

Recovery of salivary gland function vs. time after RT took place mostly within 1 year after RT. The RTOG Grade 3 was a reliable score to perform the NTCP modeling. The choice of NTCP model had no influence on the accuracy of predictions.

Comparison of the helical tomotherapy and MLC-based IMRT radiation modalities in treating brain and cranio-spinal tumors

The investigation of the clinical efficacy and effectiveness of Intensity Modulated Radiotherapy (IMRT) using Multileaf Collimators (MLC) and Helical Tomotherapy (HT) has been an issue of increasing interest over the past few years. In order to assess the suitability of a treatment plan, dosimetric criteria such as dose-volume histograms (DVH), maximum, minimum, mean, and standard deviation of the dose distribution are typically used. Nevertheless, the radiobiological parameters of the different tumors and normal tissues are often not taken into account. The use of the biologically effective uniform dose (D=) together with the complication-free tumor control probability (P(+)) were applied to evaluate the two radiation modalities. Two different clinical cases of brain and cranio-spinal axis cancers have been investigated by developing a linac MLC-based step-and-shoot IMRT plan and a Helical Tomotherapy plan. The treatment plans of the MLC-based IMRT were developed on the Philips treatment planning station using the Pinnacle 7.6 software release while the dedicated Tomotherapy treatment planning station was used for the HT plan. With the use of the P(+) index and the D(=) concept as the common prescription point, the different treatment plans were compared based on radiobiological measures. The tissue response probabilities were plotted against D(=) for a range of prescription doses. The applied plan evaluation method shows that in the brain cancer, the HT treatment gives slightly better results than the MLC-based IMRT in terms of optimum expected clinical outcome (P(+) of 66.1% and 63.5% for a D(=) to the PTV of 63.0 Gy and 62.0 Gy, respectively). In the cranio-spinal axis cancer, the HT plan is significantly better compared to the MLC-based IMRT plan over the clinically useful dose prescription range (P(+) of 84.1% and 28.3% for a D(=) to the PTV of 50.6 Gy and 44.0 Gy, respectively). If a higher than 5% risk for complications could be allowed, the complication-free tumor control could be increased by almost 30% compared to the initial dose prescription. In comparison to MLC based-IMRT, HT can better encompass the often large PTV while minimizing the volume of the OARs receiving high dose. A radiobiological treatment plan evaluation can provide a closer association of the delivered treatment with the clinical outcome by taking into account the dose-response relations of the irradiated tumors and normal tissues. The use of P - (D=) diagrams can complement the traditional tools of evaluation such as DVHs, in order to compare and effectively evaluate different treatment plans.

Proton radiotherapy for liver tumors: dosimetric advantages over photon plans

The purpose of the study is to dosimetrically investigate the advantages of proton radiotherapy over photon radiotherapy for liver tumors. The proton plan and the photon plan were designed using commercial treatment planning systems. The treatment target dose conformity and heterogeneity and dose-volume analyses of normal structures were compared between proton and photon radiotherapy for 9 patients with liver tumors. Proton radiotherapy delivered a more conformal target dose with slightly less homogeneity when compared with photon radiotherapy. Protons significantly reduced the fractional volume of liver receiving dose greater or equal to 30 Gy (V(30)) and the mean liver dose. The stomach and duodenal V(45) were significantly lower with the use of proton radiotherapy. The V(40) and V(50) of the heart and the maximum spinal cord dose were also significantly lower with the use of proton radiotherapy. Protons were better able to spare one kidney completely and deliver less dose to one (generally the left) kidney than photons. The mean dose to the total body and most critical structures was significantly decreased using protons when compared to corresponding photon plans. In conclusion, our study suggests the dosimetric benefits of proton radiotherapy over photon radiotherapy. These dosimetric advantages of proton plans may permit further dose escalation with lower risk of complications.

Lacrimal gland radiosensitivity in uveal melanoma patients

PURPOSE

To find a dose-volume effect for inhomogeneous irradiated lacrimal glands.

METHODS AND MATERIALS

Between 1999 and 2006, 72 patients (42 men and 30 women) were treated with fractionated stereotactic radiotherapy in a prospective, nonrandomized clinical trial (median follow-up, 32 months). A total dose of 50 Gy was given on 5 consecutive days. The mean of all Schirmer test results obtained > or =6 months after treatment was correlated with the radiation dose delivered to the lacrimal gland. Also, the appearance of dry eye syndrome (DES) was related to the lacrimal gland dose distribution.

RESULTS

Of the 72 patients, 17 developed a late Schirmer value <10 mm; 9 patients developed DES. A statistically significant relationship was found between the received median dose in the lacrimal gland vs. reduced tear production (p = 0.000) and vs. the appearance of DES (p = 0.003), respectively. A median dose of 7 Gy/fraction to the lacrimal gland caused a 50% risk of low Schirmer results. A median dose of 10 Gy resulted in a 50% probability of DES.

CONCLUSION

We found a clear dose-volume relationship for irradiated lacrimal glands with regard to reduced tear production and the appearance of DES.

Normal tissue dose-effect models in biological dose optimisation

Sophisticated radiotherapy techniques like intensity modulated radiotherapy with photons and protons rely on numerical dose optimisation. The evaluation of normal tissue dose distributions that deviate significantly from the common clinical routine and also the mathematical expression of desirable properties of a dose distribution is difficult. In essence, a dose evaluation model for normal tissues has to express the tissue specific volume effect. A formalism of local dose effect measures is presented, which can be applied to serial and parallel responding tissues as well as target volumes and physical dose penalties. These models allow a transparent description of the volume effect and an efficient control over the optimum dose distribution. They can be linked to normal tissue complication probability models and the equivalent uniform dose concept. In clinical applications, they provide a means to standardize normal tissue doses in the face of inevitable anatomical differences between patients and a vastly increased freedom to shape the dose, without being overly limiting like sets of dose-volume constraints.

Methods to calculate normal tissue complication and tumour control probabilities for fractionated inhomogeneous dose distribution of intensity-modulated radiation therapy

Objectives: This study is designed to present and evaluate radiobiological-based dose–volume histogram (DVH) reduction schemes to calculate normal tissue complication probability (NTCP) and tumour control probability (TCP) for intensity-modulated radiation therapy (IMRT).

Methods: The proposed DVH reduction schemes were derived for 2 Gy per fraction and prescribed dose per fraction for critical organs and tumours, respectively. Sample computed tomography scans were used to generate two IMRT plans to deliver 54 Gy to PTV1 and 24 Gy to PTV2 via sequential IMRT boost (SqIB) and simultaneous integrated IMRT boost (SIB) plans. Differential DVHs were used to calculate effective volumes using published values of related parameters of critical organs and prostate.

Results: NTCP values for bladder were almost zero for both IMRT plans. The plots between k and NTCP for rectum and femurs (k = 0.1–1.0) show higher NTCP for SqIB than that for SIB. The TCP decreases with increasing clonogenic cell density and is higher for SIB than that for SqIB for all clonogenic cell densities. The value of α proposed by Brenner and Hall shows very low radio sensitivity of clonogens of the prostate, which gives very low TCP for conventional doses of 70–80 Gy delivered in 7–8 weeks, even for very low clonogenic cell density in the prostate.

Conclusion: The presented DVH reduction schemes have radiobiological bearing and therefore seem to be effective in calculating fairly accurate NTCP and TCP.

Treatment planning comparison between conformal radiotherapy and helical tomotherapy in the case of locally advanced-stage NSCLC

BACKGROUND AND PURPOSE

To investigate the impact of Helical Tomotherapy (HT) upon the dose distribution when compared to our routinely delivered 3D conformal radiotherapy (CRT) in the case of patients affected by stage III non-small-cell lung cancer (NSCLC).

MATERIAL AND METHODS

Thirteen stage III inoperable NSCLC patients were scheduled to receive 61.2-70.2Gy, 1.8Gy/fraction. Two treatment techniques (HT and CRT) were considered, and in the case of CRT the dose calculation was performed using both the pencil beam (PB) and Anisotropic Analytical Algorithm (AAA) available on the Varian Eclipse planning system. Dose volume constraints for PTV coverage and OAR sparing were assessed for the HT inverse planning with the highest priority upon PTV coverage and spinal cord sparing. The three plans were compared in terms of dose-volume histograms (DVHs) and normal tissue complication probability (NTCP). A statistical analysis was performed using non-parametric Wilcoxon matched pairs tests.

RESULTS

In CRT the use of a less accurate algorithm (PB) decreased the monitor unit number by 2.4%. HT significantly improved dose homogeneity within PTV compared with CRT_AAA. For lung parenchyma V20-V40 were lower with HT, corresponding to a decrease of 7% in the risk of radiation pneumonitis. The volume of the heart and esophagus irradiated to >45-60Gy were reduced using HT plans. For eight PTs with an esophagus-PTV overlap >5%, HT significantly reduced both late and acute esophageal complication probability.

CONCLUSIONS

Our findings obtained in stage III NSCLC patients underline that HT guarantees an important sparing of lungs and esophagus, thus HT has the potential to improve therapeutic ratio, when compared with CRT, by means of dose escalation and/or combined treatment strategy. In CRT of locally advanced lung cancers, the use of a more advanced algorithm would give significantly better modeling of target dose and coverage.

Dosimetric and anatomic indicators of late rectal toxicity after high-dose intensity modulated radiation therapy for prostate cancer

We seek to identify dosimetric and anatomic indicators of late rectal toxicity in prostate cancer patients treated with intensity modulated radiation therapy (IMRT). Data from 49 patients sampled from 698 patients treated for clinically localized prostate cancer at the Memorial Sloan-Kettering Cancer Center with IMRT to a dose of 81 Gy were analyzed. The end point of the study was late Grade 2 or worse rectal toxicity within 30 months of treatment. Dosimetric analysis was performed on the rectum surface in three dimensions and on two-dimensional dose maps obtained by flattening the rectum surface using a conformal mapping procedure. Several parameters including the percentage and absolute surface area of the rectum irradiated, mean dose as a function of location on the rectum, planning target volume (PTV) size and rectum size were analyzed for correlation to toxicity. Significance was set at p < 0.05 for a two-sided t-test. Correlation between absolute areas irradiated and toxicity was observed on both the rectum surface and flattened rectum. Patients with toxicity also received a significantly higher mean dose to the superior 25% of the rectum surface and 15% of the flattened rectum. PTV volume, PTV height, rectum surface area and average cross-sectional area were significantly larger in patients with toxicity. The conformal mapping procedure has potential utility for evaluating dose to the rectum and risk of toxicity. Late rectal toxicity was related to the irradiation of the upper part of the rectum and also to the absolute area irradiated, PTV size, and rectum size on the planning computed tomography (CT) scan.

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.

A constrained non-rigid registration algorithm for use in prostate image-guided radiotherapy

A constrained non-rigid registration (CNRR) algorithm for use in updating prostate external beam image-guided radiotherapy treatment plans is presented in this paper. The developed algorithm is based on a multi-resolution cubic B-spline FFD transformation and has been tested and verified using 3D CT images from 10 sets of real patient data acquired from 4 different patients on different treatment days. The registration can be constrained to any combination of the prostate, rectum, bladder, pelvis, left femur, and right femur. The CNRR was tested with 5 different combinations of constraints and each test significantly outperformed both rigid and non-rigid registration at aligning constrained bones and critical organs. The CNRR was then used to update the treatment plans to account for articulated, rigid bone motion and non-rigid organ deformation. Each updated treatment plan outperformed the original treatment plan by increasing radiation dosage to the prostate and lowering radiation dosage to the rectum and bladder.

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.

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.

Involved-field radiation therapy for inoperable non small-cell lung cancer

PURPOSE

Dose escalation has been shown to improve local control in non-small-cell lung cancer (NSCLC) treated with definitive radiation therapy, but with increased complications. We implemented the use of involved-field radiotherapy (IFRT) in an effort to reduce toxicity while treating the gross tumor to higher doses. This analysis reports failure rates in uninvolved nodal regions with the use of IFRT.

PATIENTS AND METHODS

A total of 524 patients with NSCLC treated with three-dimensional conformal radiotherapy at Memorial Sloan-Kettering Cancer Center between 1991 and 2005 were reviewed. Only lymph node regions initially involved with tumor by either biopsy or radiographic criteria were included in the clinical target volume. Elective nodal failure (ENF) was defined as a recurrence in an initially uninvolved lymph node in the absence of local failure.

RESULTS

Only 32 patients (6.1%) with ENF were identified. The 2-year actuarial rates of elective nodal control and primary tumor control were 92.4% and 51%, respectively, with a median follow-up of 41 months in survivors. In patients who achieved local disease control, the 2-year elective nodal control rate was 91%. The median time to ENF was 6 months (range, 0 to 56 months). Many patients experienced treatment failure in multiple lymph node regions simultaneously.

CONCLUSION

The use of IFRT did not cause a significant amount of failure in lymph node regions not included in the tumor volume. Therefore, IFRT remains an acceptable method of treatment that allows for dose escalation while minimizing toxicity.

Potential for intensity-modulated radiation therapy to permit dose escalation for canine nasal cancer

We evaluated the impact of inverse planned intensity-modulated radiation therapy (IMRT) on the dose-volume histograms (DVHs) and on the normal tissue complication probabilities (NTCPs) of brain and eyes in dogs with nasal tumors. Nine dogs with large, caudally located nasal tumors were planned using conventional techniques and inverse planned IMRT for a total prescribed dose of 52.5 Gy in 3.5 Gy fractions. The equivalent uniform dose for brain and eyes was calculated to estimate the normal tissue complication probability (NTCP) of these organs. The NTCP values as well as the DVHs were used to compare the treatment plans. The dose distribution in IMRT plans was more conformal than in conventional plans. The average dose delivered to one-third of the brain was 10 Gy lower with the IMRT plan compared with conventional planning. The mean partial brain volume receiving 43.6 Gy or more was reduced by 25.6% with IMRT. As a consequence, the NTCPs were also significantly lower in the IMRT plans. The mean NTCP of brain was two times lower and at least one eye could be saved in all patients planed with IMRT. Another possibility with IMRT is dose escalation in the target to improve tumor control while keeping the NTCPs at the same level as for conventional planning. Veterinary

High-Dose Radiotherapy for the Treatment of Inoperable Non–Small Cell Lung Cancer

PURPOSE

Local failure continues to be a major cause of mortality in patients with inoperable non-small cell lung cancer (NSCLC) treated with radiation therapy. Dose escalation is one method of improving local control. We investigated whether high-dose radiotherapy enhances outcomes in patients with inoperable NSCLC.

MATERIALS AND METHODS

Eighty-two patients with inoperable NSCLC stage I-IIIB were treated with three-dimensional conformal radiotherapy to doses of > or =80 Gy. Patients were divided into 2 groups based on stage: those with stage I/II disease (55 patients) and those with stage III (IIIA or IIIB) disease (27 patients).

RESULTS

The 5-year local control and overall survival rates for the patients with stage I/II disease were 67% and 36%, respectively, with a median survival time of 41 months. For the patients with stage III disease, 5-year local control and overall survival rates were observed to be 39% and 31%, respectively, with a median survival time of 32 months.

CONCLUSIONS

Our data show a favorable 5-year overall survival rate (36%) with an acceptable toxicity profile in patients with early-stage NSCLC treated to doses of > or =80 Gy using three-dimensional conformal radiotherapy. Sequential chemotherapy combined with high-dose radiation gave survival rates equivalent to those seen with concurrent chemoradiation therapy in locally advanced disease. The overall survival and local control rates observed among patients with all stages of disease are consistent with and comparable to results from other dose-escalation studies reported in the literature.

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.

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.