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

[Comparison of biological functional assessment in intensity-modulated radiotherapy: two-dimensional study]

Modulated intensity dose distributions are obtained by inverse planning. It requires an inversion algorithm and an objective function that can be physical or biological. The biological objective functions aim at quantifying the probability of the favourable end of the treatment. The inversion algorithm used is analytical and is based on the mathematical analysis of the singular values decomposition. It proposes as many solutions as there are elementary beams. From the Tumour Control Probability, Normal Tissue Complication Probabilities and complication free tumour control, three biological assessment functions of the proposed solutions are compared with the least square difference between the prescribed and obtained dose distributions. We used a simplified irradiation configuration: Brahme's dose prescription (2D modelling of a prostate) and 9 beams (1D). The choice by mean of biological criterion of the optimal solution makes it possible to increase the average dose in the tumour, so as its homogeneity compared to physical optimisation. Conversely, the organs at risk are then less protected. The laying down of relevant constraints makes it possible to obtain satisfactory dose distributions. Concerning the validity of the models and data used, some limitations appear. At present time, it seems to exclude the use in clinical routine of an only biological optimisation. The future availability of new biological data will allow the development and in particular the clinical use of biological optimisation.

Induction chemotherapy plus three-dimensional conformal radiation therapy in the definitive treatment of locally advanced non-small-cell lung cancer

PURPOSE

To evaluate our institution's experience using chemotherapy in conjunction with three-dimensional conformal radiation therapy (3D-CRT).

METHODS AND MATERIALS

From 1991 to 1998, 152 patients with Stage III non-small-cell lung cancer (NSCLC) were treated with 3D-CRT at Memorial Sloan-Kettering Cancer Center. A total of 137 patients (90%) were surgically staged with either thoracotomy or mediastinoscopy. The remainder were staged radiographically. Seventy patients were treated with radiation therapy alone, and 82 patients received induction chemotherapy before radiation. The majority of chemotherapy-treated patients received a platinum-containing regimen. Radiation was delivered with a 3D conformal technique using CT-based treatment planning. The median dose in the radiation alone group was 70.2 Gy, while in the combined modality group, it was 64.8 Gy.

RESULTS

The median follow-up time was 30.5 months among survivors. Stage IIIB disease was present in 36 patients (51%) in the radiation-alone group and 57 patients (70%) in the combined-modality group. Thirty-nine patients had poor prognostic factors (KPS < 70 or weight loss > 5%), and they were equally distributed between the two groups. The median survival times for the radiation-alone and the combined-modality groups were 11.7 months and 18.1 months, respectively (p = 0.001). The 2-year rates of local control in the radiation-alone and combined-modality groups were 35.4% and 43.1%, respectively (p = 0.1). Grade 3 or worse nonhematologic toxicity occurred in 20% of the patients receiving radiation alone and in 16% of those receiving chemotherapy and radiation. Overall, there were only 4 cases of Grade 3 or worse esophagitis.

CONCLUSION

Despite more Stage IIIB patients in the combined-modality group, the addition of chemotherapy to 3D-CRT produced a survival advantage over 3D-CRT alone in Stage III NSCLC without a concomitant increase in toxicity. Chemotherapy thus appears to be beneficial, even in patients who are receiving higher doses of radiation therapy than are typically given with conventional techniques. Because locoregional failure remains a major challenge in patients with advanced disease, 3D-CRT in conjunction with chemotherapy may allow safe treatment to the dose levels required to further enhance local control.

Three-Dimensional Conformal Radiation Therapy (3D-CRT) for Early-Stage Non–Small-Cell Lung Cancer

The standard treatment for early-stage non-small-cell lung cancer is surgical resection. However, many patients are inoperable due to medical comorbidities. Thirty-two medically inoperable patients with early-stage non-small-cell lung cancer were treated with 3-dimensional conformal radiation therapy between January 1991 and December 2000. The median dose was 70.2 Gy, and the median follow-up time in survivors was 30 months. The 2-year actuarial local control, overall survival, and cancer-specific survival rates were 43%, 54%, and 57%, respectively. The 5-year actuarial local control, overall survival, and cancer-specific survival rates were 43%, 33%, and 39%, respectively. This report suggests that local control is improved with high-dose conformal radiation therapy when compared to other institutions' retrospective experiences.

Isotropic beam bouquets for shaped beam linear accelerator radiosurgery

In stereotactic radiosurgery and radiotherapy treatment planning, the steepest dose gradient is obtained by using beam arrangements with maximal beam separation. We propose a treatment plan optimization method that optimizes beam directions from the starting point of a set of isotropically convergent beams, as suggested by Webb. The optimization process then individually steers each beam to the best position, based on beam's-eye-view (BEV) critical structure overlaps with the target projection and the target's projected cross sectional area at each beam position. This final optimized beam arrangement maintains a large angular separation between adjacent beams while conformally avoiding critical structures. As shown by a radiosurgery plan, this optimization method improves the critical structure sparing properties of an unoptimized isotropic beam bouquet, while maintaining the same degree of dose conformity and dose gradient. This method provides a simple means of designing static beam radiosurgery plans with conformality indices that are within established guidelines for radiosurgery planning, and with dose gradients that approach those achieved in conventional radiosurgery planning.

Pulmonary complications following different radiotherapy techniques for breast cancer, and the association to irradiated lung volume and dose

PURPOSE

This study investigates the incidence of short-term pulmonary complications following radiotherapy (RT) for breast cancer (BC) with different treatment techniques/incidentally irradiated lung volumes and the importance of confounding factors on RT-induced pulmonary complications.

PATIENTS AND METHODS

Prospectively, 475 patients with BC were followed for pulmonary complications 1, 4 and 7 months post-RT. Mean lung dose volume histograms (MDVH) were constructed and compared for the different RT-techniques. Among a subset of the mastectomized patients treated with loco-regional (LR-) RT, who had undergone complete three-dimensional (3-D) dose planning (n = 43), MDVH for asymptomatic patients was compared with MDVH for patients experiencing both radiological and clinical pulmonary side-effects.

RESULTS

Moderate pulmonary complications, that is requiring treatment with corticosteroids, were rare following local RT (< 1%), but were diagnosed among 11% of the patients treated with LR-RT. A correlation between increasing irradiated lung volumes at the >20 Gy-level (V20), based on MDVH for the RT-techniques, and pulmonary complications was found (P < 0.001). Furthermore, increasing age and reduced pre-RT functional level were independently associated with a higher rate of pulmonary complications (P = 0.005 and P = 0.018). Among the subgroup of mastectomized patients treated with LR-RT, who had undergone complete 3-D dose planning, a difference in mean V20 was found between patients experiencing both clinical and radiological pulmonary side-effects compared to patients experiencing neither of the two side-effects (P = 0.007).

CONCLUSION

Moderate pulmonary complications following local RT for BC are rare. The incidence of short-term moderate pulmonary complications in LR-RT is, however, clinically significant and to define quality assurance guidelines for these RT-techniques, 3-D RT planning can be used.

The controversies and pitfalls in modeling normal tissue radiation injury/damage

Highly conformal fields have become achievable in routine clinical practice. The optimal shape of the resultant dose distributions depends on information that is not currently available. This missing information is the dose-volume response of the normal tissues at risk. These functions are now the subject of aggressive research. The research involves collecting the dose-response data, modeling the dose-response function, and fitting the models to the data. The controversies addressed here influence the selection of the biomathematical model that one might use to describe such a function. The form that the dose-volume response function takes depends on the nature of the volume effect. The nature of the volume effect for a given radiation response is the subject of considerable debate. Related to this debate, this report addresses the existence of the volume effect, the existence of a threshold volume, and the existence of functional subunits. The pitfalls relate to the problems in accurate determination and application of the dose-response functions.

Can dose-response models predict reliable normal tissue complication probabilities in radical radiotherapy of urinary bladder cancer? The impact of alternative radiation tolerance models and parameters

PURPOSE

To analyze the consequences of selecting alternative normal tissue complication probability (NTCP) models and parameters for evaluation of radiotherapy of urinary bladder cancer.

METHODS AND MATERIALS

Treatment plans of 24 bladder cancer patients referred to radical 4-field conformal radiotherapy were analyzed. Small intestinal and rectal NTCPs were determined using both the probit and relative seriality models with several sets of published radiation tolerance parameters. Various combinations of NTCP models and parameters were applied to find the prescription dose in individual patients as well as to estimate the benefit of the conformal radiotherapy setup.

RESULTS

Different risk estimates were predicted from the two NTCP models, even when the same clinical radiation tolerance doses were fitted into the two models. The demonstrated variability translated into significant deviations (7-10 Gy) in the recommended prescription doses. Even if it was possible to discriminate between a 2-field plan and the 4-field conformal setup using a given complication model and set of tolerance parameters, the estimated benefit of the conformal treatment in terms of permitted dose escalation varied with as much as 10-12 Gy between the different NTCP models/parameters used.

CONCLUSION

Different NTCP models and tolerance parameters might propose different answers to important clinical questions in radiotherapy treatment of bladder cancer, such as dose prescription and scoring of rival treatment plans. We therefore recommend that the variability caused by tolerance parameter uncertainty and model selection should be taken into account in dose-response modeling of radiotherapy treatment.

A prospective study of salivary function sparing in patients with head-and-neck cancers receiving intensity-modulated or three-dimensional radiation therapy: initial results

OBJECTIVES

In a prospective clinical study, we tested the hypothesis that sparing the parotid glands may result in significant objective and subjective improvement of xerostomia in patients with head-and-neck cancers. The functional outcome 6 months after the completion of radiation therapy is presented.

METHODS AND MATERIALS

From February 1997 to February 1999, 41 patients with head-and-neck cancers were enrolled in a prospective salivary function study. Inverse-planning intensity-modulated radiation therapy (IMRT) was used to treat 27 patients, and forward-planning three-dimensional radiation therapy in 14. To avoid potential bias in data interpretation, only patients whose submandibular glands received greater than 50 Gy were eligible. Attempts were made to spare the superficial lobe of the parotid glands to avoid underdosing tumor targets in the parapharyngeal space; however, the entire parotid volume was used to compute dose-volume histograms (DVHs) for this analysis. DVHs were computed for each gland separately. Parotid function was assessed objectively by measuring stimulated and unstimulated saliva flow before and 6 months after the completion of radiation therapy. Measurements were converted to flow rate (mL/min) and normalized relative to that before treatment. The corresponding quality-of-life (QOL) outcome was assessed by five questions regarding the patient's oral discomfort and eating/speaking problems.

RESULTS

We observed a correlation between parotid mean dose and the fractional reduction of stimulated saliva output at 6 months after the completion of radiation therapy. We further examined whether the functional outcome could be modeled as a function of dose. Two models were found to describe the dose-response data well. The first model assumed that each parotid gland is comprised of multiple independent parallel functional subunits (corresponding to computed tomography voxels) and that each gland contributes equally to overall flow, and that saliva output decreases exponentially as a quadratic function of irradiation dose to each voxel. The second approach uses the equivalent uniform dose (EUD) metrics, which assumes loss of salivary function with increase in EUD for each parotid gland independently. The analysis suggested that the mean dose to each parotid gland is a reasonable indicator for the functional outcome of each gland. The corresponding exponential coefficient was 0.0428/Gy (95% confidence interval: 0.01, 0.09). The QOL questions on eating/speaking function were significantly correlated with stimulated and unstimulated saliva flow at 6 months. In a multivariate analysis, a toxicity score derived from the model based on radiation dose to the parotid gland was found to be the sole significant predictive factor for xerostomia. Neither radiation technique (IMRT vs. non-IMRT) nor chemotherapy (yes or no) independently influenced the functional outcome of the salivary glands.

CONCLUSION

Sparing of the parotid glands translates into objective and subjective improvement of both xerostomia and QOL scores in patients with head-and-neck cancers receiving radiation therapy. Modeling results suggest an exponential relationship between saliva flow reduction and mean parotid dose for each gland. We found that the stimulated saliva flow at 6 months after treatment is reduced exponentially, for each gland independently, at a rate of approximately 4% per Gy of mean parotid dose.

Dose escalation in non-small-cell lung cancer using three-dimensional conformal radiation therapy: update of a phase I trial

PURPOSE

High-dose radiation may improve outcomes in non-small-cell lung cancer (NSCLC). By using three-dimensional conformal radiation therapy and limiting the target volume, we hypothesized that the dose could be safely escalated.

MATERIALS AND METHODS

A standard phase I design was used. Five bins were created based on the volume of normal lung irradiated, and dose levels within bins were chosen based on the estimated risk of radiation pneumonitis. Starting doses ranged from 63 to 84 Gy given in 2.1-Gy fractions. Target volumes included the primary tumor and any nodes >or= 1 cm on computed tomography. Clinically uninvolved nodal regions were not included purposely. More recently, selected patients received neoadjuvant cisplatin and vinorelbine.

RESULTS

At the time of this writing, 104 patients had been enrolled. Twenty-four had stage I, four had stage II, 43 had stage IIIA, 26 had stage IIIB, and seven had locally recurrent disease. Twenty-five received chemotherapy, and 63 were assessable for escalation. All bins were escalated at least twice. Although grade 2 radiation pneumonitis occurred in five patients, grade 3 radiation pneumonitis occurred in only two. The maximum-tolerated dose was only established for the largest bin, at 65.1 Gy. Dose levels for the four remaining bins were 102.9, 102.9, 84 and 75.6 Gy. The majority of patients failed distantly, though a significant proportion also failed in the target volume. There were no isolated failures in clinically uninvolved nodal regions.

CONCLUSION

Dose escalation in NSCLC has been accomplished safely in most patients using three-dimensional conformal radiation therapy, limiting target volumes, and segregating patients by the volume of normal lung irradiated.

Potential role of intensity modulated proton beams in prostate cancer radiotherapy

PURPOSE

The present study was undertaken to assess the potential benefit of intensity modulated (IM) proton beams in optimizing the dose distribution to safely escalate the tumor dose in prostate cancer radiotherapy.

METHODS AND MATERIALS

Four treatment plans were compared in a prostate cancer patient aiming to deliver 81 Gy to the target: 1) conformal 18 MV X-rays, 6-fields; 2) 214 MeV protons, 2-fields; 3) IM 15 MV X-rays, 5-fields; and 4) 177-200 Mev IM protons, 5-fields as in Plan 3. In addition, IM methods were used to further escalate the tumor dose to 99 Gy. Dose-volume histograms (DVH) were used to physically compare the treatment plans. DVH data were also used to obtain normal tissue complication probabilities (NTCP) for the rectum, bladder, femoral heads, and tumor control probabilities.

RESULTS

Although the planning target volume dose distribution was satisfactory with the four treatment plans, the homogeneity was slightly reduced in both X-ray plans (IM and standard) and the low-to-medium doses delivered to all organs at risk, and other normal tissues were significantly reduced by both proton plans. For a prescribed dose of 81 Gy, only the IM X-ray and IM proton plans both succeeded in predicting an acceptably low NTCP for the rectum (<5%, Grade 3). The integral nontarget dose was significantly reduced with IM proton beams (i.e., 3.1, 1.3, and 1.7 times less than Plans 1, 2, and 3, respectively). When escalating the dose to 99 Gy, no additional improvement between IM protons and IM X-ray beams was observed.

CONCLUSION

Both IM X-ray and proton beams were able to optimize the dose distribution and comply with the goal of delivering the highest dose to the target while reducing the risk of severe morbidity to acceptable levels. The main advantage compared to IM X-rays was that IM protons succeeded in significantly reducing the low-to-medium dose to the nontarget tissues and achieved a small improvement in planning target volume (PTV) dose heterogeneity.

Loss of heterozygosity at the mannose 6-phosphate insulin-like growth factor 2 receptor (M6P/IGF2R) locus predisposes patients to radiation-induced lung injury

PURPOSE

To investigate the relationship between loss of heterozygosity (LOH) at the mannose 6-phosphate/insulin-like growth factor 2 receptor (M6P/IGF2R) gene locus and the development of radiation-induced lung injury.

MATERIAL AND METHODS

Thirty-five lung cancer patients with both stored plasma for Transforming Growth Factor beta1 (TGFbeta1) analysis and sufficient quantities of archival pathology tissue to screen for LOH were studied. All patients had been treated with thoracic radiotherapy for their malignancy and had radiographically detectable tumor present before beginning radiotherapy. Tumor and normal cells were microdissected from archival lung cancer pathology specimens. Two polymorphisms in the 3' untranslated region of the M6P/IGF2R were used to screen for LOH. Plasma TGFbeta1 levels were measured using acid-ethanol extraction and an ELISA. TGFbeta1 and M6P/IGF2R protein expression was estimated by immunofluorescence and immunohistochemical staining. Symptomatic radiation pneumonitis was scored according to National Cancer Institute Common Toxicity Criteria without knowledge of the results of TGFbeta or LOH analyses.

RESULTS

Of the 35 patients, 10 were homozygous for this polymorphism (noninformative) and were excluded. Of the 25 informative patients, 13 had LOH. Twelve of 13 patients with LOH had increased pretreatment plasma TGFbeta1 levels, vs. 3/12 patients without LOH (p < 0.01). A decrease or loss of M6P/IGF2R protein in the malignant cell accompanied by increased latent TGFbeta1 protein in extracellular matrix and tumor stroma was found in tumors with LOH, suggesting that this mutation resulted in loss of function of the receptor. Seven of 13 (54%) LOH patients developed symptomatic radiation-induced lung injury vs. 1/12 (8%) of patients without LOH (p = 0.05).

CONCLUSION

Loss of the M6P/IGF2R gene strongly correlates with the development of radiation pneumonitis after thoracic radiotherapy (RT). Furthermore, patients with LOH (in the setting of measurable tumor) are much more likely to have elevated plasma TGFbeta, suggesting an inability to normally process this cytokine. Thus, loss of the M6P/IGF2R gene may predispose patients to the development of radiation-induced lung injury.

Technical aspects of the deep inspiration breath-hold technique in the treatment of thoracic cancer

PURPOSE

The goal of this paper is to describe our initial experience with the deep inspiration breath-hold (DIBH) technique in conformal treatment of non-small-cell lung cancer with particular emphasis on the technical aspects required for implementation.

METHODS AND MATERIALS

In the DIBH technique, the patient is verbally coached through a modified slow vital capacity maneuver and brought to a reproducible deep inspiration breath-hold level. The goal is to immobilize the tumor and to expand normal lung out of the high-dose region. A physicist or therapist monitors and records patient breathing during simulation, verification, and treatment using a spirometer with a custom computer interface. Examination of internal anatomy during fluoroscopy over multiple breath holds establishes the reproducibility of the DIBH maneuver for each patient. A reference free-breathing CT scan and DIBH planning scan are obtained. To provide an estimate of tumor motion during normal tidal breathing, additional scan sets are obtained at end inspiration and end expiration. These are also used to set the spirometer action levels for treatment. Patient lung inflation is independently verified over the course of treatment by comparing the distance from the isocenter to the diaphragm measured from the DIBH digitally reconstructed radiographs to the distance measured on the portal films. Patient breathing traces obtained during treatment were examined retrospectively to assess the reproducibility of the technique.

RESULTS

Data from the first 7 patients, encompassing over 250 treatments, were analyzed. The inferred displacement of the centroid of gross tumor volume from its position in the planning scan, as calculated from the spirometer records in over 350 breath holds was 0.02 +/- 0.14 cm (mean and standard deviation). These data are consistent with the displacements of the diaphragm (-0.1 +/- 0.4 cm; range, from -1.2 to 1.1 cm) relative to the isocenter, as measured on the (92) portal films. The latter measurements include the patient setup error. The patient averaged displacement of the tumor during free breathing, determined from the tumor displacement between end inspiration and end expiration, was 0.8 +/- 0.5 cm in both the superior-inferior and anterior-posterior directions and 0.1 cm (+/- 0.1 cm) medial-laterally.

CONCLUSION

Treatment of patients with the DIBH technique is feasible in a clinical setting. With this technique, consistent lung inflation levels are achieved in patients, as judged by both spirometry and verification films. Breathing-induced tumor motion is significantly reduced using DIBH compared to free breathing, enabling better target coverage.

[Conformal radiotherapy: tomotherapy]

Conformal radiation therapy allows the possibility of delivering high doses at the tumor volume whilst limiting the dose to the surrounding tissues and diminishing the secondary effects. With the example of the conformal radiation therapy used at the AZ VUB (3DCRT and tomotherapy), two treatment plans of a left ethmoid carcinoma will be evaluated and discussed in detail. The treatment of ethmoid cancer is technically difficult for both radiation therapy and surgery because of the anatomic constraints and patterns of local spread. A radiation therapy is scheduled to be delivered after surgical resection of the tumor. The treatment plan for the radiation therapy was calculated on a three-dimensional (3D) treatment planning system based on virtual simulation with a beam's eye view: George Sherouse's Gratis. An effort was made to make the plan as conformal and as homogeneous as possible to deliver a dose of 66 Gy in 33 fractions at the tumor bed with a maximum dose of 56 Gy to the right optic nerve and the chiasma. To establish the clinical utility and potential advantages of tomotherapy over 3DCRT for ethmoid carcinoma, the treatment of this patient was also planned with Peacock Plan. For both treatment plans the isodose distributions and cumulative dose volume histograms (CDVH) were computed. Superimposing the CDVHs yielded similar curves for the target and an obvious improvement for organs at risk such as the chiasma, brainstem and the left eye when applying tomotherapy. These results have also been reflected in the tumor control probabilities (equal for both plans) and the normal tissue complication probabilities (NTCP), yielded significant reductions in NTCP for tomotherapy. The probability of uncomplicated tumor control was 52.7% for tomotherapy against 38.3% for 3DCRT.

The physical parameters and molecular events associated with radiation-induced lung toxicity

Radiation therapy (RT) is frequently used to treat patients with tumors in and around the thorax. Clinical radiation pneumonitis is a common side effect, occurring in 5% to 20% of patients. Efforts to identify patients at risk for pneumonitis have focused on physical factors, such as dose and volume. Recently, the underlying molecular biological mechanisms behind RT-induced lung injury have come under study. Improved knowledge of the molecular events associated with RT-induced lung injury may translate into a better ability to individualized therapy. This review discusses our current understanding of the physical and molecular factors contributing to RT-induced pulmonary injury.

The deep inspiration breath-hold technique in the treatment of inoperable non-small-cell lung cancer

PURPOSE

Conventional radiotherapeutic techniques are associated with lung toxicity that limits the treatment dose. Motion of the tumor during treatment requires the use of large safety margins that affect the feasibility of treatment. To address the control of tumor motion and decrease the volume of normal lung irradiated, we investigated the use of three-dimensional conformal radiation therapy (3D-CRT) in conjunction with the deep inspiration breath-hold (DIBH) technique.

METHODS AND MATERIALS

In the DIBH technique, the patient is initially maintained at quiet tidal breathing, followed by a deep inspiration, a deep expiration, a second deep inspiration, and breath-hold. At this point the patient is at approximately 100% vital capacity, and simulation, verification, and treatment take place during this phase of breath-holding.

RESULTS

Seven patients have received a total of 164 treatment sessions and have tolerated the technique well. The estimated normal tissue complication probabilities decreased in all patients at their prescribed dose when compared to free breathing. The dose to which patients could be treated with DIBH increased on average from 69.4 Gy to 87.9 Gy, without increasing the risk of toxicity.

CONCLUSIONS

The DIBH technique provides an advantage to conventional free-breathing treatment by decreasing lung density, reducing normal safety margins, and enabling more accurate treatment. These improvements contribute to the effective exclusion of normal lung tissue from the high-dose region and permit the use of higher treatment doses without increased risks of toxicity.

Correlation of radiation tolerance dose of normal human organs with organ weight, blood, and water content

The concerted effort to minimize the radiation exposure to normal human tissues while delivering a high radiation dose to the tumor often results in complications. This limits the efficacy of radiation treatment. Analysis of radiation tolerance dose with organ weight in 15 human organs yields a correlation coefficient of 0.62, whereas the correlation of radiation tolerance dose with blood and water content yields correlation coefficients of 0.82 and 0.60, respectively. Results indicate that as the organ weight and/or blood and water content increases, radiation tolerance dose decreases.

Beam intensity modulation using tissue compensators or dynamic multileaf collimation in three-dimensional conformal radiotherapy of primary cancers of the oropharynx and larynx, including the elective neck

INTRODUCTION

The treatment of midline tumors in the head and neck by conventional radiotherapy almost invariably results in xerostomia. This study analyzes whether a simple three-dimensional conformal radiotherapy (3D-CRT) technique with beam intensity modulation (BIM) (using a 10-MV beam of the MM50 Racetrack Microtron) can spare parotid and submandibular glands without compromising the dose distribution in the planning target volume (PTV).

METHODS

For 15 T2 tumors of the tonsillar fossa with extension into the soft palate (To) and 15 T3 tumors of the supraglottic larynx (SgL), conventional treatment plans, consisting of lateral parallel opposed beams, were used for irradiation of both the primary tumor (70 Gy) and the elective neck regions (46 Gy). Separately, for each tumor a 3-D conformal treatment plan was developed using the 3-D computer planning system, CadPlan, and Optimize, a noncommercial program to compute optimal beam profiles. Beam angles were selected with the intention of optimal sparing of the salivary glands. The intensity of the beams was then modulated to achieve a homogeneous dose distribution in the target for the given 3D-CRT techniques. The dose distributions, dose-volume histograms (DVHs) of target and salivary glands, tumor control probabilities (TCPs), salivary gland volumes absorbing a biologically equivalent dose of greater than 40 or 50 Gy, and normal tissue complication probabilities (NTCPs) of each treatment plan were computed. The parameters of the 3D-CRT plans were compared with those of the conventional plans.

RESULTS

In comparison with the conventional technique, the dose homogeneity in the target volume was improved by the conformal technique for both tumor sites. In addition, for the SgL conformal technique, the average volumes of the parotid glands absorbing a BED of greater than 40 Gy (V40) decreased by 23%, and of the submandibular glands by 7% (V40) and 6% (V50). Consequently, the average NTCPs for the parotid and submandibular glands were reduced by 7% and 6%, respectively. For the To conformal techniques, the V40 of the parotid glands was decreased on average by 31%, resulting in an average reduction of the NTCP by 49%. Both the average V50 and the NTCP of the submandibular glands were decreased by 7%.

CONCLUSION

For primary tumors of the oropharynx, the parotid glands could be spared to a considerable degree with the 3D-CRT technique. However, particularly the ipsilateral submandibular gland could not be spared. For primary tumors of the larynx, the 3D-CRT technique allows sparing of all salivary glands to a considerable and probably clinically relevant degree. Moreover, the conformal techniques resulted in an increased dose homogeneity in the PTV of both tumor sites.

Calculation of cranial nerve complication probability for acoustic neuroma radiosurgery

PURPOSE

Estimations of complications from stereotactic radiosurgery usually rely simply on dose-volume or dose-diameter isoeffect curves. Due to the sparse clinical data available, these curves have typically not considered the target location in the brain, target histology, or treatment plan conformality as parameters in the calculation. In this study, a predictive model was generated to estimate the probability of cranial neuropathies as a result of acoustic schwannoma radiosurgery.

METHODS AND MATERIALS

The dose-volume histogram reduction scheme was used to calculate the normal tissue complication probability (NTCP) from brainstem dose-volume histograms. The model's fitting parameters were optimized to provide the best fit to the observed complication data for acoustic neuroma patients treated with stereotactic radiosurgery at the University of Florida. The calculation was then applied to the remainder of the patients in the database.

RESULTS

The best fit to our clinical data was obtained using n = 0.04, m = 0.15, and alpha/beta = 2.1 Gy(-1). Although the fitting parameter m is relatively consistent with ranges found in the literature, both the volume parameter, n, and alpha/beta are much smaller than the values quoted in the literature. The fit to our clinical data indicates that brainstem, or possibly a specific portion of the brainstem, is more radiosensitive than the parameters in the literature indicate, and that there is very little volume effect; in other words, irradiation of a small fraction of the brainstem yields NTCPs that are nearly as high as those calculated for entire volume irradiation. These new fitting parameters are specific to acoustic neuroma radiosurgery, and the small volume effect that we observe may be an artifact of the fixed relationship of acoustic tumors to specific regions of the brainstem. Applying the model to our patient database, we calculate an average NTCP of 7.2% for patients who had no cranial nerve complications, and the average NTCP for was 66% for patients who sustained a cranial neuropathy. For the entire patient population, the actual percentage of patients suffering either facial or trigeminal neuropathy was 14.7%, whereas the calculated average NTCP was 14.8%.

DISCUSSION

NTCP calculations using brainstem dose-volume histograms can be used to estimate the rate of cranial neuropathies from acoustic neuroma radiosurgery. More clinical data and further study will lead to refinement of the model with time.

A semi-empirical treatment planning model for optimization of multiprobe cryosurgery

A model is presented for treatment planning of multiprobe cryosurgery. In this model a thermal simulation algorithm is used to generate temperature distribution from cryoprobes, visualize isotherms in the anatomical region of interest (ROI) and provide tools to assist estimation of the amount of freezing damage to the target and surrounding normal structures. Calculations may be performed for any given freezing time for the selected set of operation parameters. The thermal simulation is based on solving the transient heat conduction equation using finite element methods for a multiprobe geometry. As an example, a semi-empirical optimization of 2D placement of six cryoprobes and their thermal protocol for the first freeze cycle is presented. The effectiveness of the optimized treatment protocol was estimated by generating temperature-volume histograms and calculating the objective function for the anatomy of interest. Two phantom experiments were performed to verify isotherm locations predicted by calculations. A comparison of the predicted 0 degrees C isotherm with the actual iceball boundary imaged by x-ray CT demonstrated a spatial agreement within +/-2 mm.

Late rectal toxicity after conformal radiotherapy of prostate cancer (I): multivariate analysis and dose-response

PURPOSE

The purpose of this paper is to use the outcome of a dose escalation protocol for three-dimensional conformal radiation therapy (3D-CRT) of prostate cancer to study the dose-response for late rectal toxicity and to identify anatomic, dosimetric, and clinical factors that correlate with late rectal bleeding in multivariate analysis.

METHODS AND MATERIALS

Seven hundred forty-three patients with T1c-T3 prostate cancer were treated with 3D-CRT with prescribed doses of 64.8 to 81.0 Gy. The 5-year actuarial rate of late rectal toxicity was assessed using Kaplan-Meier statistics. A retrospective dosimetric analysis was performed for patients treated to 70.2 Gy (52 patients) or 75.6 Gy (119 patients) who either exhibited late rectal bleeding (RTOG Grade 2/3) within 30 months after treatment (i.e., 70.2 Gy-13 patients, 75. 6 Gy-36 patients) or were nonbleeding for at least 30 months (i.e., 70.2 Gy-39 patients, 75.6 Gy-83 patients). Univariate and multivariate logistic regression was performed to correlate late rectal bleeding with several anatomic, dosimetric, and clinical variables.

RESULTS

A dose response for >/= Grade 2 late rectal toxicity was observed. By multivariate analysis, the following factors were significantly correlated with >/= Grade 2 late rectal bleeding for patients prescribed 70.2 Gy: 1) enclosure of the outer rectal contour by the 50% isodose on the isocenter slice (i.e., Iso50) (p < 0.02), and 2) smaller anatomically defined rectal wall volume (p < 0.05). After 75.6 Gy, the following factors were significant: 1) smaller anatomically defined rectal wall volume (p < 0.01), 2) higher rectal D(max) (p < 0.01), 3) enclosure of rectal contour by Iso50 (p < 0.01), 4) patient age (p = 0.02), and 5) history of diabetes mellitus (p = 0.04). In addition to these five factors, acute rectal toxicity was also significantly correlated (p = 0.05) with late rectal bleeding when patients from both dose groups were combined in multivariate analysis.

CONCLUSION

A multivariate logistic regression model is presented which describes the probability of developing late rectal bleeding after conformal irradiation of prostate cancer. Late rectal bleeding correlated with factors which may indicate that a greater fractional volume of rectal wall was exposed to high dose, such as smaller rectal wall volume, inclusion of the rectum within the 50% isodose on the isocenter slice, and higher rectal D(max).

Influence of patient positioning on dose-volume histogram and normal tissue complication probability for small bowel and bladder in patients receiving pelvic irradiation: a prospective study using a 3D planning system and a radiobiological model

PURPOSE

A prospective study was undertaken to evaluate the influence of patient positioning (prone position using a belly board vs. supine position) on the dose-volume histograms (DVHs) of organs of risk, and to analyze its possible clinical relevance using radiobiological models.

METHODS AND MATERIALS

From November 1996 to August 1997 a computed tomography (CT) scan was done in the prone position using a belly board and in supine position in 20 consecutive patients receiving postoperative pelvic irradiation because of rectal cancer. Using a three-dimensional (3D) planning system (Helax, TMS) the DVH for small bowel, bladder, a standard planning target volume (PTV) of postoperative irradiation of rectal cancer, the intersection of volume of PTV and small bowel (PTV intersection V(SB), respectively, of PTV and bladder (PTV intersection V(B)) were defined in each axial CT slice. The normal tissue complication probability (NTCP) was determined by the radiobiological model of Lyman and Kutcher using the tolerance data of Emami. For evaluation of late toxicity alpha/beta ratio was 2.5; for evaluation of acute toxicity, it was 10. Total dose was 50.4 Gy (1.8 Gy/fraction) (ICRU 50).

RESULTS

Using the prone position compared to the supine position, the median volume of PTV intersection V(B) was reduced by 18.5 cm3 (62%). Median dose (related to the reference dose) to the bladder was 44.5% (22.4 Gy) in prone and 66.05% (33.3 Gy) in supine position (p<0.001). Median V(B) within the 90% (45.4 Gy), 80% (40.3 Gy), 60% (30.2 Gy), and 40% (20.2 Gy) isodose was significantly lower in the prone position when compared to the supine position. Using the radiobiological models, however, there was no difference of NTCP between prone position or supine position. In the prone position, median volume of PTV intersection V(SB) was reduced by 32.5 cm3 (54%). The median dose to small bowel was 30.85% (15.4 Gy) in the prone position and 47.35% (23.9Gy) in the supine position (p<0.001). Significant differences between prone and supine position were found for median V(SB) within the 90% (45.4 Gy), 80% (40.3 Gy), 60% (30.2 Gy), and 40% (20.2 Gy) isodose. According to the method of Lyman, median NTCP of small bowel was significant lower in prone than in supine position.

CONCLUSION

The prone position with a standard belly board should be the standard positioning technique for patients receiving adjuvant postoperative radiation therapy following surgery of rectal cancer. Both irradiated volume and total dose to the organs of risk can be reduced significantly. As a consequence of this, radiation induced toxicity will be minimized.

The rationale and use of three-dimensional radiation treatment planning for lung cancer

Treatment of lung cancer with conventional radiation therapy is associated with suboptimal local tumor control and poor long-term survival. Poor local tumor control may result from inaccurate tumor targeting, failure to satisfactorily conform to dose distribution with the target volume, and/or inadequate radiation doses. Three-dimensional treatment planning is a radiotherapy technique that provides more accurate dose targeting via the direct transfer of three-dimensional anatomic information from diagnostic scans into the planning process. This technology can assist treatment planning by providing dose-volume histograms, an estimation of normal tissue complication probabilities, and facilitate dose escalation. Preliminary clinical studies suggest that this is a feasible approach worthy of additional study. The three-dimensional tools provide new opportunities to better understand radiation-induced changes in pulmonary function.

Bone marrow doses and leukaemia risk in radiotherapy of prostate cancer

BACKGROUND AND PURPOSE

As more and more patients with prostate cancer are cured and survive with only minor chronic morbidity, other potentially treatment related morbidity, in particular second cancers, becomes an urgent problem which may influence decisions on treatment strategy and treatment plan optimisation. Epidemiological data suggest a radiotherapy associated risk of AML in prostate cancer patients of approximately 0.1% in 10 years. The aim of the study was to determine the range of bone marrow doses from different treatment plans and in different patients in order to develop criteria for optimisation of treatment plans in conformal radiotherapy of prostate cancer to further minimise the small risk of secondary leukaemia.

MATERIALS AND METHODS

Doses to the pelvic bone marrow were calculated for eight different plans used in radiotherapy of prostate cancer to determine the variability of bone marrow doses in radiotherapy of prostate cancer. Computer tomography (CT) slices of the entire pelvic region of an Alderson phantom were acquired and transferred to the TPS. Critical bone marrow structures were outlined in each slice. Different treatment plans were evaluated on this phantom and dose-volume histograms (DVH) for the pelvic bone marrow were obtained. Similarly, the DVH for the bone marrow of 14 patients who received conformal radiotherapy for prostate cancer was determined.

RESULTS

Mean total bone marrow doses ranged from 3.4 to 5.6 Gy in the phantom study. Approximately 99% of the mean dose to the total bone marrow comes from the dose to bone marrow located in the pelvic bones and lumbar vertebrae. Mean bone marrow doses of 14 patients given the same conformal radiotherapy plan ranged from 3.5 to 7.7 Gy.

CONCLUSIONS

No correlation was found between the rectum normal tissue complication probability (NTCP) and the mean bone marrow dose. This means that in the process of treatment planning, exposure to both critical organs, the rectum as well as the bone marrow, should be minimised independently to arrive at the optimal treatment plan.

Dose-wall histograms and normalized dose-surface histograms for the rectum: a new method to analyze the dose distribution over the rectum in conformal radiotherapy

PURPOSE

To develop an accurate method to generate a dose-volume histogram (DVH) of the rectum wall, solely based on the outer contours of the rectum wall.

METHODS AND MATERIALS

A mathematical model for the rectum wall is developed, incorporating the stretching of the rectum wall due to variable rectal filling and neighboring structures. The model is based on the assumption that the amount of intersected rectum wall tissue normal to the central axis of the rectum is constant. The main objective of the model is to determine the thickness of the rectum wall in each wall element. Two approaches are described, each yielding a DVH of the rectum wall, based only on the delineated outer contours of the rectum. In the first approach, the model is used to create a set of inner contours out of the axial outer contours. Both sets of contours are used to derive a dose-wall histogram (DWH) of the rectum. In the second approach, the model is used to generate a normalized 2D sampling space, which is subsequently binned into a normalized dose-surface histogram (NDSH). The model is verified using 20 sets of CT data (5 patients x 4 scans) in which both outer and inner contours of the rectum are carefully delineated. The DWHs and NDSHs are compared with DVHs of the rectum wall, which require contouring of the outer and inner surfaces of the rectum wall, and with DVHs of the total rectum (including rectal filling). The variation between DWHs, NDSHs, and DVHs is investigated using normal tissue complication probability (NTCP) calculations.

RESULTS

The local wall thickness of the rectum as outlined on CT data was in conformity with the described rectum model. The amount of rectum wall tissue per unit length rectum varied considerably between patients (27%, 1 SD). In all analyzed patients, the DWHs and NDSHs corresponded well to the DVHs of the rectum wall. Much more discrepancies were observed between the DVHs of the total rectum and the DVHs of the rectum wall.

CONCLUSION

The applied methods yield accurate dose distributions of the rectum wall, without delineating the inner surface of the rectum. This reduces both the workload and variations due to inaccurate delineation of the rectum wall. The DWH and NDSH are effective tools to evaluate 3D dose distributions of the rectum wall and to estimate the complication probability of the rectum in high-dose conformal radiotherapy.

Deep inspiration breath-hold technique for lung tumors: the potential value of target immobilization and reduced lung density in dose escalation

PURPOSE/OBJECTIVE

This study evaluates the dosimetric benefits and feasibility of a deep inspiration breath-hold (DIBH) technique in the treatment of lung tumors. The technique has two distinct features--deep inspiration, which reduces lung density, and breath-hold, which immobilizes lung tumors, thereby allowing for reduced margins. Both of these properties can potentially reduce the amount of normal lung tissue in the high-dose region, thus reducing morbidity and improving the possibility of dose escalation.

METHODS AND MATERIALS

Five patients treated for non-small cell lung carcinoma (Stage IIA-IIIB) received computed tomography (CT) scans under 4 respiration conditions: free-breathing, DIBH, shallow inspiration breath-hold, and shallow expiration breath-hold. The free-breathing and DIBH scans were used to generate 3-dimensional conformal treatment plans for comparison, while the shallow inspiration and expiration scans determined the extent of tumor motion under free-breathing conditions. To acquire the breath-hold scans, the patients are brought to reproducible respiration levels using spirometry, and for DIBH, modified slow vital capacity maneuvers. Planning target volumes (PTVs) for free-breathing plans included a margin for setup error (0.75 cm) plus a margin equal to the extent of tumor motion due to respiration (1-2 cm). Planning target volumes for DIBH plans included the same margin for setup error, with a reduced margin for residual uncertainty in tumor position (0.2-0.5 cm) as determined from repeat fluoroscopic movies. To simulate the effects of respiration-gated treatments and estimate the role of target immobilization alone (i.e., without the benefit of reduced lung density), a third plan is generated from the free-breathing scan using a PTV with the same margins as for DIBH plans.

RESULTS

The treatment plan comparison suggests that, on average, the DIBH technique can reduce the volume of lung receiving more than 25 Gy by 30% compared to free-breathing plans, while respiration gating can reduce the volume by 18%. The DIBH maneuver was found to be highly reproducible, with intra breath-hold reproducibility of 1.0 (+/- 0.9) mm and inter breath-hold reproducibility of 2.5 (+/- 1.6) mm, as determined from diaphragm position. Patients were able to perform 10-13 breath-holds in one session, with a comfortable breath-hold duration of 12-16 s.

CONCLUSION

Patients tolerate DIBH maneuvers well and can perform them in a highly reproducible fashion. Compared to conventional free-breathing treatment, the DIBH technique benefits from reduced margins, as a result of the suppressed target motion, as well as a decreased lung density; both contribute to moving normal lung tissue out of the high-dose region. Because less normal lung tissue is irradiated to high dose, the possibility for dose escalation is significantly improved.

The local management of soft tissue sarcoma

Soft tissue sarcomas (STS) are rare tumors arising from the connective tissues. STS can arise at any anatomic site, can demonstrate varied behavior and prognosis, and therefore present a formidable challenge in management. The local treatment of STS demands technical complexity in the application of diagnostic tools, including pathology and imaging, as well as treatment approaches, including surgical ablation and reconstruction, radiotherapy, and, in defined cases, chemotherapy. The understanding of the management of these lesions is profoundly dependent on the multidisciplinary setting, where experience has been gained and skills are available to increase the likelihood of a successful result. Several proven options are available for optimal local management, and the choice of approach depends on the prevailing practice and resource profile of the treating center. With modern approaches, the local control rate can be expected to be at least 90% for extremity lesions, which constitute the most common STS. The experience in other anatomic sites is less favorable as a result of a combination of late diagnosis, technically difficult access sites, and possibly less familiarity with these less common presentations. The disappointing results make it all the more important for patients to be referred to a multidisciplinary setting with experience in sarcoma management to maximize the chance of successful local outcome.

Clinical dose-volume histogram analysis for pneumonitis after 3D treatment for non-small cell lung cancer (NSCLC)

PURPOSE

To identify a clinically relevant and available parameter upon which to identify non-small cell lung cancer (NSCLC) patients at risk for pneumonitis when treated with three-dimensional (3D) radiation therapy.

METHODS AND MATERIALS

Between January 1991 and October 1995, 99 patients were treated definitively for inoperable NSCLC. Patients were selected for good performance status (96%) and absence of weight loss (82%). All patients had full 3D treatment planning (including total lung dose-volume histograms [DVHs]) prior to treatment delivery. The total lung DVH parameters were compared with the incidence and grade of pneumonitis after treatment.

RESULTS

Univariate analysis revealed the percent of the total lung volume exceeding 20 Gy (V20), the effective volume (Veff) and the total lung volume mean dose, and location of the tumor primary (upper versus lower lobes) to be statistically significant relative to the development of > or = Grade 2 pneumonitis. Multivariate analysis revealed the V20 to be the single independent predictor of pneumonitis.

CONCLUSIONS

The V20 from the total lung DVH is a useful parameter easily obtained from most 3D treatment planning systems. The V20 may be useful in comparing competing treatment plans to evaluate the risk of pneumonitis for our individual patient treatment and may also be a useful parameter upon which to stratify patients or prospective dose escalation trials.

Dose escalation of chart in non-small cell lung cancer: is three-dimensional conformal radiation therapy really necessary?

PURPOSE

To evaluate, preclinically, the potential for dose escalation of continuous, hyperfractionated, accelerated radiation therapy (CHART) for non small-cell lung cancer (NSCLC), we examined the strategy of omission of elective nodal irradiation with and without the application of three-dimensional conformal radiation technology (3DCRT).

METHODS AND MATERIALS

2D, conventional therapy plans were designed according to the specifications of CHART for 18 patients with NSCLC (Stages Ib, IIb, IIIa, and IIIb). Further plans were generated with the omission of elective nodal irradiation (ENI) from the treatment portals (2D minus ENI plans [2D-ENI plans]). Both sets were inserted in the patient's planning computed tomographies (CTs). These reconstructed plans were then compared to alternative, three-dimensional treatment plans which had been generated de novo, with the omission of ENI: 3D minus elective nodal irradiation (3D-ENI plans). Dose delivery to the planning target volumes (PTVs) and to the organs at risk were compared between the 3 sets of corresponding plans. The potential for dose escalation of each patient's 2D-ENI and 3D-ENI plan beyond 54 Gy, standard to CHART, was also determined.

RESULTS

PTV coverage was suboptimal in the 2D CHART and the 2D-ENI plans. Only in the 3D-ENI plans did 100% of the PTV get > or = 95% of the dose prescribed (i.e., 51.5 Gy [51.3-52.2]). Using 3D-ENI plans significantly reduced the dose received by the spinal cord, the mean and median doses to the esophagus and the heart. It did not significantly reduce the lung dose when compared to 2D-ENI plans. Escalation of the dose (minimum > or = 1 Gy) with optimal PTV coverage was possible in 55.5% of patients using 3D-ENI, but was possible only in 16.6% when using the 2D-ENI planning strategy.

CONCLUSIONS

3DCRT is fundamental to achieving optimal PTV coverage in NSCLC. A policy of omission of elective nodal irradiation alone (and using 2D technology) will not achieve optimal PTV coverage or dose escalation. 3DCRT with omission of ENI can achieve true escalation of CHART in 55.5% of tumors, depending on their site and N-stage.

[Conformal radiotherapy: principles and classification]

'Conformal radiotherapy' is the name fixed by usage and given to a new form of radiotherapy resulting from the technological improvements observed during, the last ten years. While this terminology is now widely used, no precise definition can be found in the literature. Conformal radiotherapy refers to an approach in which the dose distribution is more closely 'conformed' or adapted to the actual shape of the target volume. However, the achievement of a consensus on a more specific definition is hampered by various difficulties, namely in characterizing the degree of 'conformality'. We have therefore suggested a classification scheme be established on the basis of the tools and the procedures actually used for all steps of the process, i.e., from prescription to treatment completion. Our classification consists of four levels: schematically, at level 0, there is no conformation (rectangular fields); at level 1, a simple conformation takes place, on the basis of conventional 2D imaging; at level 2, a 3D reconstruction of the structures is used for a more accurate conformation; and level 3 includes research and advanced dynamic techniques. We have used our personal experience, contacts with colleagues and data from the literature to analyze all the steps of the planning process, and to define the tools and procedures relevant to a given level. The corresponding tables have been discussed and approved at the European level within the Dynarad concerted action. It is proposed that the term 'conformal radiotherapy' be restricted to procedures where all steps are at least at level 2.

A neural network to predict symptomatic lung injury

A nonlinear neural network that simultaneously uses pre-radiotherapy (RT) biological and physical data was developed to predict symptomatic lung injury. The input data were pre-RT pulmonary function, three-dimensional treatment plan doses and demographics. The output was a single value between 0 (asymptomatic) and 1 (symptomatic) to predict the likelihood that a particular patient would become symptomatic. The network was trained on data from 97 patients for 400 iterations with the goal to minimize the mean-squared error. Statistical analysis was performed on the resulting network to determine the model's accuracy. Results from the neural network were compared with those given by traditional linear discriminate analysis and the dose-volume histogram reduction (DVHR) scheme of Kutcher. Receiver-operator characteristic (ROC) analysis was performed on the resulting network which had Az = 0.833 +/- 0.04. (Az is the area under the ROC curve.) Linear discriminate multivariate analysis yielded an Az = 0.813 +/- 0.06. The DVHR method had Az = 0.521 +/- 0.08. The network was also used to rank the significance of the input variables. Future studies will be conducted to improve network accuracy and to include functional imaging data.

Conformal radiation therapy treatment planning reduces the dose to the optic structures for patients with tumors of the paranasal sinuses

PURPOSE

Compare dose distributions of traditional versus conformal beam orientations in paranasal sinus malignancies.

MATERIALS AND METHODS

Maximum normal tissue doses, dose volume histograms (DVH), normal tissue complication probabilities (NTCP), and the percentage of each normal tissue receiving >80% of the average target dose (V80) were calculated.

RESULTS/CONCLUSIONS

Conformal planning reduced the V80 to the optic nerves and chiasm as well as the normal tissue maximum doses to the ipsilateral and contralateral optic nerves and chiasm, and mean NTCPs.

A method for incorporating organ motion due to breathing into 3D dose calculations

A method is proposed that incorporates the effects of intratreatment organ motion due to breathing on the dose calculations for the treatment of liver disease. Our method is based on the convolution of a static dose distribution with a probability distribution function (PDF) which describes the nature of the motion. The organ motion due to breathing is assumed here to be one-dimensional (in the superior-inferior direction), and is modeled using a periodic but asymmetric function (more time spent at exhale versus inhale). The dose distribution calculated using convolution-based methods is compared to the static dose distribution using dose difference displays and the effective volume (Veff) of the uninvolved liver, as per a liver dose escalation protocol in use at our institution. The convolution-based calculation is also compared to direct simulations that model individual fractions of a treatment. Analysis shows that incorporation of the organ motion could lead to changes in the dose prescribed for a treatment based on the Veff of the uninvolved liver. Comparison of convolution-based calculations and direct simulation of various worst-case scenarios indicates that a single convolution-based calculation is sufficient to predict the dose distribution for the example treatment plan given.

[Inverse radiotherapy planning]

BACKGROUND

In clinical practice it sometimes happens that with currently available conformal radiotherapy techniques no satisfactory dose distribution can be achieved. In these cases inverse radiotherapy planning and intensity modulated radiotherapy may give better solutions.

METHOD

Inverse planning is a technique using a computer program to automatically achieve a treatment plan which has an optimal merit. This merit may either depend on dose or dose-volume constraints like minimum and maximum doses in the target region or critical organs, respectively, or biological indices like the complication free tumor control rate. As the result of inverse planning the inhomogeneous intensity fluence of the beams is calculated. These fluence distributions may be generated by beam compensators or multi-leaf collimation.

RESULTS

Clinical studies to prove the advantage of inverse planning are already on the way. It has been shown that this technology is safe and that the dose distributions which can be achieved are superior to conventional methods.

CONCLUSIONS

Inverse treatment planning and intensity modulated radiation therapy will almost certainly come to be the technique of choice for selected clinical cases.

Beam intensity modulation for penumbra enhancement in the treatment of lung cancer

PURPOSE

A treatment planning study was performed for patients with lung cancer in order to investigate the extent to which doses to critical structures could be reduced by penumbra enhancement at the superior and inferior field edges, using beam intensity modulation (BIM) with a multileaf collimator. By applying two independent published models for the prediction of the incidence of normal tissue complications, the potential for dose escalation without increasing the incidence of pneumonitis was estimated.

METHODS AND MATERIALS

For 12 patients, the standard treatment technique was compared with the BIM technique using the Cadplan 3D planning system (Varian-Dosetek). Dose distributions in the healthy lung tissue were evaluated by considering both lungs minus the tumor as one functional unit. The following parameters were compared: (i) the average normalized total dose (NTD), (ii) the lung volume receiving an NTD of more than 20 Gy, and (iii) the calculated normal tissue complication probability (NTCP).

RESULTS

Due to the applied BIM technique, the field lengths could be reduced by 1.4 cm for all patients, while achieving a minimum dose at the superior and inferior parts of the target of 95% of the isocenter dose. Compared to the standard technique, BIM reduced the patient mean of the average NTD for the healthy lung tissue from 16.5 to 15.3 Gy. The volume of healthy lung tissue receiving an NTD of 20 Gy or more was reduced by 9.7% (range 2.2 to 23.1%). The calculated NTCP reduced from 10.7% to 7.6% on average. The length of the esophagus that received a dose of 60 Gy or more could be reduced for 5 of the 6 stage III patients in this study. Based on equal lung NTCPs for the standard technique and the BIM technique, a mean dose escalation of 5.7 Gy (range 1.1 to 16.0 Gy) was possible for the 12 patients in this study. Based on equal average NTDs for the two techniques, the patient mean of the allowed dose escalation was 6.5 Gy (range 1.1 to 18.2 Gy). All dose escalations would be possible without exceeding the spinal cord tolerance dose.

CONCLUSIONS

The BIM technique reduced the dose delivery to critical tissues. Two published methods for estimating the incidence of pneumonitis both pointed to a potential for dose escalation of 6 to 7 Gy on average with the BIM technique, without increasing the incidence of pneumonitis. For 2 of the 12 patients in this study the estimated allowed dose escalation even exceeded 15 Gy.

Treatment plan evaluation using dose-volume histogram (DVH) and spatial dose-volume histogram (zDVH)

OBJECTIVE

The dose-volume histogram (DVH) has been accepted as a tool for treatment-plan evaluation. However, DVH lacks spatial information. A new concept, the z-dependent dose-volume histogram (zDVH), is presented as a supplement to the DVH in three-dimensional (3D) treatment planning to provide the spatial variation, as well as the size and magnitude of the different dose regions within a region of interest.

MATERIALS AND METHODS

Three-dimensional dose calculations were carried out with various plans for three disease sites: lung, breast, and prostate. DVHs were calculated for the entire volume. A zDVH is defined as a differential dose-volume histogram with respect to a computed tomographic (CT) slice position. In this study, zDVHs were calculated for each CT slice in the treatment field. DVHs and zDVHs were compared.

RESULTS

In the irradiation of lung, DVH calculation indicated that the treatment plan satisfied the dose-volume constraint placed on the lung and zDVH of the lung revealed that a sizable fraction of the lung centered about the central axis (CAX) received a significant dose, a situation that warranted a modification of the treatment plan due to the removal of one lung. In the irradiation of breast with tangential fields, the DVH showed that about 7% of the breast volume received at least 110% of the prescribed dose (PD) and about 11% of the breast received less than 98% PD. However, the zDVHs of the breast volume in each of seven planes showed the existence of high-dose regions of 34% and 15%, respectively, of the volume in the two caudal-most planes and cold spots of about 40% in the two cephalic planes. In the treatment planning of prostate, DVHs showed that about 15% of the bladder and 40% of the rectum received 102% PD, whereas about 30% of the bladder and 50% of the rectum received the full dose. Taking into account the hollow structure of both the bladder and the rectum, the dose-surface histograms (DSH) showed larger hot-spot volume, about 37% of the bladder wall and 43% of the rectal wall. The zDVHs of the bladder revealed that the hot-spot region was superior to the central axis. The zDVHs of the rectum showed that the high-dose region was an 8-cm segment mostly superior to the central axis. The serial array-like of the rectum warrants a closer attention with regard to the complication probability of the organ.

CONCLUSIONS

Although DVH provides an averaged dose-volume information, zDVH provides differential dose-volume information with respect to the CT slice position. zDVH is a 2D analog of a 3D DVH and, in some situations, more superior. It provides additional information on plan evaluation that otherwise could not be appreciated. The zDVH may be used along with DVH for plan evaluation and for the correlation of radiation outcome.

An objective function for radiation treatment optimization based on local biological measures

The implementation of biological optimization of radiation treatment plans is impeded by both computational and modelling problems. We derive an objective function from basic model assumptions which includes the normal tissue constraints as interior penalty functions. For organs that are composed of parallel subunits, a mean response model is proposed which leads to constraints similar to dose-volume constraints. This objective function is convex in the case when no parallel organs lie in the treatment volume. Otherwise, an argument is given to show that a number of local minima may exist which are near degenerate to the global minimum. Thus, together with the measure quality of the objective function, highly efficient gradient algorithms can be used. The number of essential biological model parameters could be reduced to a minimum. However, if the optimization constraints are given as TCP/NTCP values, Lagrange multiplier updates have to be performed by invoking comprehensive biological models.

Optimization and clinical use of multisegment intensity-modulated radiation therapy for high-dose conformal therapy

Intensity-modulated radiation therapy (IMRT) may be performed with many different treatment delivery techniques. This article summarizes the clinical use and optimization of multisegment IMRT plans that have been used to treat more than 350 patients with IMRT over the last 4.5 years. More than 475 separate clinical IMRT plans are reviewed, including treatments of brain, head and neck, thorax, breast and chest wall, abdomen, pelvis, prostate, and other sites. Clinical planning, plan optimization, and treatment delivery are summarized, including efforts to minimize the number of additional intensity-modulated segments needed for particular planning protocols. Interactive and automated optimization of segmental and full IMRT approaches are illustrated, and automation of the segmental IMRT planning process is discussed.

Radiation pneumonitis as a function of mean lung dose: an analysis of pooled data of 540 patients

PURPOSE

To determine the relation between the incidence of radiation pneumonitis and the three-dimensional dose distribution in the lung.

METHODS AND MATERIALS

In five institutions, the incidence of radiation pneumonitis was evaluated in 540 patients. The patients were divided into two groups: a Lung group, consisting of 399 patients with lung cancer and 1 esophagus cancer patient and a Lymph./Breast group with 78 patients treated for malignant lymphoma, 59 for breast cancer, and 3 for other tumor types. The dose per fraction varied between 1.0 and 2.7 Gy and the prescribed total dose between 20 and 92 Gy. Three-dimensional dose calculations were performed with tissue density inhomogeneity correction. The physical dose distribution was converted into the biologically equivalent dose distribution given in fractions of 2 Gy, the normalized total dose (NTD) distribution, by using the linear quadratic model with an alpha/beta ratio of 2.5 and 3.0 Gy. Dose-volume histograms (DVHs) were calculated considering both lungs as one organ and from these DVHs the mean (biological) lung dose, NTDmean, was obtained. Radiation pneumonitis was scored as a complication when the pneumonitis grade was grade 2 (steroids needed for medical treatment) or higher. For statistical analysis the conventional normal tissue complication probability (NTCP) model of Lyman (with n=1) was applied along with an institutional-dependent offset parameter to account for systematic differences in scoring patients at different institutions.

RESULTS

The mean lung dose, NTDmean, ranged from 0 to 34 Gy and 73 of the 540 patients experienced pneumonitis, grade 2 or higher. In all centers, an increasing pneumonitis rate was observed with increasing NTDmean. The data were fitted to the Lyman model with NTD50=31.8 Gy and m=0.43, assuming that for all patients the same parameter values could be used. However, in the low dose range at an NTDmean between 4 and 16 Gy, the observed pneumonitis incidence in the Lung group (10%) was significantly (p=0.02) higher than in the Lymph./Breast group (1.4%). Moreover, between the Lung groups of different institutions, also significant (p=0.04) differences were present: for centers 2, 3, and 4, the pneumonitis incidence was about 13%, whereas for center 5 only 3%. Explicitly accounting for these differences by adding center-dependent offset values for the Lung group, improved the data fit significantly (p < 10(-5)) with NTD50=30.5+/-1.4 Gy and m=0.30+/-0.02 (+/-1 SE) for all patients, and an offset of 0-11% for the Lung group, depending on the center.

CONCLUSIONS

The mean lung dose, NTDmean, is relatively easy to calculate, and is a useful predictor of the risk of radiation pneumonitis. The observed dose-effect relation between the NTDmean and the incidence of radiation pneumonitis, based on a large clinical data set, might be of value in dose-escalating studies for lung cancer. The validity of the obtained dose-effect relation will have to be tested in future studies, regarding the influence of confounding factors and dose distributions different from the ones in this study.

Plasma transforming growth factor beta1 as a predictor of radiation pneumonitis

PURPOSE

To investigate prospectively the utility of plasma transforming growth factor beta1 (TGFbeta1) as a marker for the development of symptomatic radiation pneumonitis.

MATERIALS AND METHODS

Seventy-three patients with lung cancer treated with curative intent are reported herein. Plasma TGFbeta1 samples were obtained before, weekly during, and at each follow-up after radiation therapy (RT). TGFbeta1 was extracted using an acid/ethanol method. An enzyme-linked immunosorbent assay was used to quantify plasma TGFbeta1 concentrations. The TGFbeta1 level at the end of RT was considered "normal" if it was both < or = 7.5 ng/ml and less than the pretreatment value. All patients were followed for at least 6 months, unless symptomatic pneumonitis developed sooner. Pneumonitis was defined by National Cancer Institute (NCI) common toxicity criteria.

RESULTS

Fifteen of the 73 patients (21%) developed symptomatic pneumonitis and the remaining 58 (79%) did not. A normal plasma TGFbeta1 by the end of RT, as defined above, was more common in patients who did not develop pneumonitis. A return of the plasma TGFbeta1 to normal accurately identified patients who would not develop pneumonitis with both a sensitivity and positive predictive value of 90%.

CONCLUSION

Plasma TGFbeta1 levels appear to be a useful means to identify patients at low risk for the development of pneumonitis from thoracic RT. Thus, monitoring of plasma TGFbeta1 levels may identify candidates for dose escalation studies in the treatment of lung cancer.

Evaluation of two dose-volume histogram reduction models for the prediction of radiation pneumonitis

PURPOSE

To evaluate the similarities between the mean lung dose and two dose-volume histogram (DVH) reduction techniques of 3D dose distributions of the lung.

PATIENTS AND METHODS

DVHs of the lungs were calculated from 3D dose distributions of patients treated for malignant lymphoma (44), breast cancer (42) and lung cancer (20). With a DVH reduction technique, a DVH is summarized by the equivalent uniform dose (EUD), a quantity which is directly related to the normal tissue complication probability (NTCP). Two DVH reduction techniques were used. The first was based on an empirical model proposed by Kutcher et al. (Kutcher, G.J., Burman, C., Brewster, M.S., Goitein, M. and Mohan, R. Histogram reduction method for calculating complication probabilities for three-dimensional treatment planning evaluations. Int. J. Radiat. Oncol. Biol. Phys. 21: 137-146, 1991), which needs a volume exponent n. Several values for n were tested. The second technique was based on a radiobiological model, the parallel functional subunit model developed by Niemierko et al. (Niemierko, A. and Goitein, M. Modeling of normal tissue response to radiation: the critical volume model. Int. J. Radiat. Oncol. Biol. Phys. 25: 135-145, 1993) and Jackson et al. (Jackson, A., Kutcher, G.J. and Yorke, E.D. Probability of radiation-induced complications for normal tissues with parallel architecture subject to non-uniform irradiation. Med. Phys. 20: 613-625, 1993), for which a local dose-effect relation needed to be specified. This relation was obtained from an analysis of perfusion and ventilation SPECT data.

RESULTS

It can be shown analytically that the two DVH reduction techniques are identical, if the local dose-effect relation obeys a power-law relationship in the clinical dose range. Local dose-effect relations based on perfusion and ventilation SPECT data can indeed be fitted with a power-law relationship in the range 0-80 Gy, from which values of n = 0.8-0.9 were deduced. These correspond to the commonly used value of n = 0.87 for lung tissue and yielded EUDn=0.87 values which were almost identical to the mean lung doses. For other n values, for which no experimental data are present, differences exist between EUD and mean dose values. Six patients with malignant lymphoma (6/44) and none of the breast cancer patients (0/42) developed radiation pneumonitis. These cases occurred only at high values for the mean lung dose.

CONCLUSION

The two DVH reduction techniques are identical for lung and are very similar to mean dose calculations. The two techniques are also relatively similar for other model parameter values.

A test of the claim that plan rankings are determined by relative complication and tumor-control probabilities

PURPOSE

This study tests an accepted claim regarding tumor control (TCP) and normal tissue complication (NTCP) probability functions. The claim is that treatment plans can be ranked using relative probabilities, even when the absolute probabilities are unknown. The assumption supports the use of probability models for plan optimization and the comparison of treatment techniques.

METHODS

The claim was tested using a hypothetical model consisting of two tissues, and illustrated with clinical data. Plans were scored using the probability of uncomplicated tumor control. The scores of different plans were compared by fixing their relative risks for an individual tissue complication, but adjusting the absolute probability levels up or down. The tested claim is that the plan rankings should not change.

RESULTS

In the two-tissue model, the rankings of competing plans were reversed by doubling all the probabilities. The preference ordering of lung cancer plans changed after the risk of pulmonary complication was reduced by 3-fold. In another site, the ranking of plans by overall complication-free probability was disturbed by errors that preserved the ordering of plans with respect to any individual complication. An adjustment of +/- 2.5% in the initial NTCP values for two tissues changed the direction in which a plan score moved in response to a fixed tradeoff in complication risk in an optimization search.

CONCLUSIONS

Contrary to claims, plan rankings are not determined by the relative probabilities of adverse events. The effect on plan scores of trading one complication for another depends on the absolute levels of risk. Absolute errors in NTCP and TCP functions result in the wrong ranking of plans, even when relative probabilities are correct. An optimization routine based on TCP and NTCP calculations may be forced in the wrong direction by small errors in the probability estimates.

The relevance of transforming growth factor beta 1 in pulmonary injury after radiation therapy

The maximum dose of radiation which can be delivered to a tumor is limited by the tolerance of the surrounding normal tissues. The ability to identify patients at high or low risk of injury from radiation therapy might enable the clinician to tailor radiation doses in order to maximize efficacy and minimize risk. The cytokine transforming growth factor beta 1 (TGF beta 1) has been implicated in the development of normal tissue injury after irradiation in several organs, including the lung. Herein, the evidence supporting the role of TGF beta 1 in radiation-induced lung injury is reviewed. Using the treatment of non-small cell lung cancer as a model, we also discuss how it may be possible to identify patients at risk for this complication using measurements of plasma TGF beta 1, and how this information may be used in the future to adjust doses of radiation in the treatment of lung cancer.

[Conformational radiotherapy with multi-leaf collimators: one year experience at the Leon-Berard Centre]

Taking advantage of the renewal of a linear accelerator, the Radiation Therapy Department of the Centre Léon Bérard implemented, in collaboration with Philips Systèmes Médicaux, a conformal therapy set-up procedure using CT-scan for 3D treatment planning and a multileaf collimator that allows achievement of numerous irregular-shaped beams via the multileaf preparation system. The various elements of this equipment make possible well defined and structured procedures for treatment planning with different steps and essential tools used by this technique. We describe the means used and indicate future improvements that will lead to automation in order to provide good quality assurance, better security and substantial time saving. During the first year, 115 patients were treated with this new technique. They presented with central nervous system tumors (32 patients), lung cancer (29 patients), prostate cancer (20 patients), paranasal sinus tumors (14 patients) and tumors located in other sites (13 patients with soft sarcoma, hepato-bilary tumor, etc).

Comparison of radiosurgery treatment modalities based on complication and control probabilities

PURPOSE

The relative efficacy of Gamma Knife, Linac, and Proton treatment modalities for stereotactic radiosurgery (SRS) was investigated on the basis of normal tissue complication probability (NTCP) and tumor control probability (TCP), calculated for representative test cases.

METHODS AND MATERIALS

Five radiosurgery patient cases were selected to cover a range of treatment-planning situations from small spherical volumes to large irregular volumes. A target volume consisting of contours drawn on CT transverse slices was prepared for each case. Plans were developed using the three treatment modalities for each case, with the objective of encompassing the target as closely as possible with a prescription isodose line and minimizing dose to normal tissue, within the constraints of current clinical practice. Dose-volume histograms (DVH) were calculated for the target and for normal tissue, and these histograms were used to calculate NTCP and TCP values for each plan.

RESULTS AND CONCLUSIONS

Differences in NTCP and TCP values were found to depend on treatment modality, size, shape, and location of the target, the amount of effort devoted to treatment planning, and the complexity of the plan.

Potential clinical efficacy of intensity-modulated conformal therapy

PURPOSE

The purpose of this study was to examine the potential benefit of using intensity-modulated conformal therapy for a variety of lesions currently treated with stereotactic radiosurgery or conventional radiotherapy.

METHODS AND MATERIALS

Intensity-modulated conformal treatment plans were generated for small intracranial lesions, as well as head and neck, lung, breast, and prostate cases, using the Peacock Plan treatment-planning system (Nomos Corporation). For small intracranial lesions, intensity-modulated conformal treatment plans were compared with stereotactic radiosurgery treatment plans generated for patient treatment at the University of Florida Shands Cancer Center. For other sites (head and neck, lung, breast, and prostate), plans generated using the Peacock Plan were compared with conventional treatment plans, as well as beam's-eye-view conformal treatment plans. Plan comparisons were accomplished through conventional qualitative review of two-dimensional (2D) dose distributions in conjunction with quantitative techniques, such as dose-volume histograms, dosimetric statistics, normal tissue complication probabilities, tumor control probabilities, and objective numerical scoring.

RESULTS

For small intracranial lesions, there is little difference between intensity-modulated conformal treatment planning and radiosurgery treatment planning in the conformation of high isodose lines with the target volume. However, stereotactic treatment planning provides a steeper dose gradient outside the target volume and, hence, a lower normal tissue toxicity index. For extracranial sites, objective numerical scores for beam's-eye-view and intensity-modulated conformal planning techniques are superior to scores for conventional treatment plans. The beam's-eye-view planning technique prevents geographic target misses and better excludes healthy tissues from the treatment portal. Compared with scores for the beam's-eye-view planning technique, scores for intensity-modulated conformal plans using the Peacock Plan were significantly better for the lung and head and neck cases studied, equivalent for prostate cases, and inferior for breast cases.

CONCLUSION

Using the entire 3D data set to construct radiotherapy plans through virtual simulation is always advantageous, whether done for stereotactic radiosurgery, beam's-eye-view conformal therapy, or intensity-modulated conformal treatment. Intensity modulation of the photon beam further enhances treatment planning under specific conditions. In general, the intensity-modulated technique is advantageous for large, irregular targets with critical structures in close proximity. Intensity-modulated treatment planning does not appear advantageous for stereotactic radiosurgery or treatment of the intact breast.

Electron wedges for radiation therapy

PURPOSE

Brain tumors can be advantageously treated with electron over photon radiation, by exploiting the rapid fall-off in dose with depth. This advantage could be further enhanced by utilizing multiple electron beams. However, in some beam configurations, wedged dose profiles would be necessary for the dose uniformity. Unlike photons, shaped pieces of material placed in electron beam severely degrade the energy, give additional scattering and, therefore, are suboptimal. The purpose of this study was to create wedged electron fields, using intensity modulation. The combination of electron wedges enables a more uniform coverage of brain tumors with a reduced dose to normal tissue.

METHODS AND MATERIALS

Intensity modulation was performed for 10 to 50 MeV electrons using a narrow scanning elementary beam of a racetrack Microtron accelerator, delivering radiation pulses with coordinates and intensities prescribed by a custom scan matrix. Dispensing more pulses (or longer pulses) within the field to increase the local dose, one can sharpen the penumbra at depth and generate wedged dose distributions of arbitrary angle as well as many other desired profiles. We modulated the electron beams, measured dose distributions using film in an anthropomorphic phantom, and compared the results with conventional techniques.

RESULTS

Intensity modulation of electron beams decreases the 50-90% penumbra at depth by 40% and increases the flatness by 80%. Wedged profiles at depth can be created for any angle up to about 70 degrees, depending on the beam energy. Multiple modulated electron beams give smaller 20-70% but larger 70-100% isodose regions than photon beams.

CONCLUSIONS

Electron beams can improve dose distributions in brain compared to the same number of photon beams, reducing the 20-70% isodoses region in normal tissue by 30%. Intensity modulation significantly improves the dose distribution from combined electron beams providing a sharper penumbra, better conformity, and reduced margin.

[3-Dimensional irradiation planning in brain tumors. The advantages of the method and the clinical results]

PURPOSE

Radiotherapy became an important component in the treatment of brain gliomas. The aim of this study is to analyse several advantages of the three-dimensional conformal radiation therapy in comparison with a two-dimensional conventional technique and to present the clinical results of 43 patients with brain gliomas treated according to a three-dimensional planning.

PATIENTS AND METHOD

Between January 1994 and December 1995, 43 patients with malignant brain gliomas (WHO III and IV) were treated in our department according to a three-dimensional treatment planning. The patients received a total irradiation dose of 60 Gy, 2 Gy/day, 5 days/week. The rate of survival was analysed in relation with the known prognostical factors: histology, Karnofsky index, age, resection status. In 10 patients a three-dimensional treatment planning was compared with a conventional two-dimensional planning: the volume of the normal brain tissue irradiated to high dose levels (95% isodose) and the normal tissue complication probability (NTCP) for the brain by Kutcher and Lyman were comparatively analysed.

RESULTS

The survival rate for the whole group was 14 months. The histology of the tumor, age, Karnofsky index and resection status were important prognostical factors. The three-dimensional planning allows a 15 to 20% reduction in the volume of normal brain tissue irradiated to high dose levels (95% isodose). The NTCP is significantly lower using the three-dimensional technique (range 0.03% to 13%), in comparison with the two-dimensional conventional technique (range 0.1% to 26%). The value of NTCP increases with tumor volume.

CONCLUSIONS

Concerning the tumor control and survival rate, the three-dimensional treatment planning shows no advantages compared to the standard conventional methods. The main advantage of the three-dimensional treatment planning is the possibility to spare normal brain tissue. The possibility to integrate mathematical models in the evaluation of the therapy could give this technique new dimensions.

Physical and biological predictors of changes in whole-lung function following thoracic irradiation

PURPOSE

To develop methods of predicting the pulmonary consequences of thoracic irradiation (RT) by prospectively studying changes in pulmonary function following RT.

METHODS AND MATERIALS

100 patients receiving incidental partial-lung irradiation during treatment of tumors in or adjacent to the thorax had whole-lung function assessed via symptoms and pulmonary function tests (PFTs: FEV1-forced expiratory volume 1 s; DLCO-diffusion capacity) before and repeatedly 6-48 months following RT. All had computed tomography-based three-dimensional (3D) dose calculations with lung density heterogeneity corrections for dose-volume histogram (DVH) and normal tissue complication probability (NTCP) calculations. Functional DVHs (DVfH) based on SPECT (single photon emission computed tomography) lung perfusion scans, and serial transforming growth factor-beta (TGF-beta1) levels were available in 50 and 48 patients, respectively. The incidence and severity of changes in whole-lung function were correlated with clinical, physical, and biological factors. Exploratory statistical analyses were performed using chi-square, Pearson correlations, logistic regression, and multiple linear regression.

RESULTS

RT-induced symptoms developed in 21 patients. In the overall group, the single best predictor for the development of symptoms was the NTCP (p < 0.05). Pre-RT PFTs alone were less predictive (p = 0.1 for FEV1, p = 0.08 for DLCO). A multivariate model based on pre-RT DLCO and CT-based NTCP was strongly predictive for the development of symptoms (p < 0.001). NTCPs based on SPECT-derived DVf Hs and TGF-beta1 levels did not appear to provide additional predictive value. The presence or absence of pulmonary symptoms was correlated with the decline in PFT 6 months following RT (p < 0.05). In the overall group, the degree of decline in PFTs was not well correlated with any of the dose-volume variables considered. In patients with "good" pre-RT PFTs, there was a relationship between the percent reduction in PFT and dose-volume parameters such as the percent of lung volume receiving > 30 Gy (p < 0.05).

CONCLUSION

The extent of alteration in whole-lung function (symptoms or PFT changes) appears to be related to both dose-volume and pre-RT PFT parameters. The data suggest that no one variable is likely to be an adequate predictor and that multivariate predictive models will be needed. Additional studies are underway to develop better predictive models that consider physical factors such as the DVH and regional perfusion, as well as biological/clinical factors such as pre-RT PFTs and TGF-beta1.