Randall K. Ten Haken

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.

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.

Quantitative Analyses of Normal Tissue Effects in the Clinic (QUANTEC): An Introduction to the Scientific Issues

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

Comparing different NTCP models that predict the incidence of radiation pneumonitis

Purpose: To compare different normal tissue complication probability (NTCP) models to predict the incidence of radiation pneumonitis on the basis of the dose distribution in the lung. Methods and Materials: The data from 382 breast cancer, malignant lymphoma, and inoperable non-small-cell lung cancer patients from two centers were studied. Radiation pneumonitis was scored using the Southwestern Oncology Group criteria. Dose-volume histograms of the lungs were calculated from the dose distributions that were corrected for dose per fraction effects. The dose-volume histogram of each patient was reduced to a single parameter using different local dose-effect relationships. Examples of single parameters were the mean lung dose (MLD) and the volume of lung receiving more than a threshold dose (VDth). The parameters for the different NTCP models were fit to patient data using a maximum likelihood analysis. Results: The best fit resulted in a linear local dose-effect relationship, with the MLD as the resulting single parameter. The relationship between the MLD and NTCP could be described with a median toxic dose (TD50) of 30.8 Gy and a steepness parameter m of 0.37. The best fit for the relationship between the VDth and the NTCP was obtained with a Dth of 13 Gy. The MLD model was found to be significantly better than the VDth model (p 35%. For arbitrary dose distributions, an estimate of the uncertainty in the NTCP could be determined using the probability distribution of the parameter values of the Lyman-Kutcher-Burman model. Conclusion: The maximum likelihood method revealed that the underlying local dose-effect relation for radiation pneumonitis was linear (the MLD model), rather than a step function (the VDth model). Thus, for the studied patient population, the MLD was the most accurate predictor for the incidence of radiation pneumonitis.PURPOSE To compare different normal tissue complication probability (NTCP) models to predict the incidence of radiation pneumonitis on the basis of the dose distribution in the lung. METHODS AND MATERIALS The data from 382 breast cancer, malignant lymphoma, and inoperable non-small-cell lung cancer patients from two centers were studied. Radiation pneumonitis was scored using the Southwestern Oncology Group criteria. Dose-volume histograms of the lungs were calculated from the dose distributions that were corrected for dose per fraction effects. The dose-volume histogram of each patient was reduced to a single parameter using different local dose-effect relationships. Examples of single parameters were the mean lung dose (MLD) and the volume of lung receiving more than a threshold dose (V(Dth)). The parameters for the different NTCP models were fit to patient data using a maximum likelihood analysis. RESULTS The best fit resulted in a linear local dose-effect relationship, with the MLD as the resulting single parameter. The relationship between the MLD and NTCP could be described with a median toxic dose (TD(50)) of 30.8 Gy and a steepness parameter m of 0.37. The best fit for the relationship between the V(Dth) and the NTCP was obtained with a D(th) of 13 Gy. The MLD model was found to be significantly better than the V(Dth) model (p <0.03). However, for 85% of the studied patients, the difference in NTCP calculated with both models was <10%, because of the high correlation between the two parameters. For dose distributions outside the range of the studied dose-volume histograms, the difference in NTCP, using the two models could be >35%. For arbitrary dose distributions, an estimate of the uncertainty in the NTCP could be determined using the probability distribution of the parameter values of the Lyman-Kutcher-Burman model. CONCLUSION The maximum likelihood method revealed that the underlying local dose-effect relation for radiation pneumonitis was linear (the MLD model), rather than a step function (the V(Dth) model). Thus, for the studied patient population, the MLD was the most accurate predictor for the incidence of radiation pneumonitis.

Measurement of prostate movement over the course of routine radiotherapy using implanted markers

PURPOSE To measure the range and frequency of occurrence of intertreatment movement of the prostate gland over the course of radiotherapy, and to demonstrate that the prostate may move independently of the surrounding bones of the pelvis. METHODS AND MATERIALS Ten patients underwent implantation of radiopaque markers around the prostate. Orthogonal portal films were taken at multiple stages during the course of treatment and digitized. An image registration tool was used to solve for film detector placement and, subsequently, to determine positional changes between structures on a reference portal image pair and all subsequent pairs for each patient. Transformations describing prostate movement were measured independently of those describing setup variations of the pelvic girdle. RESULTS Translation and/or rotation of the prostate was detected in 70% of the treatments for which films were taken. The maximum measured displacement was 7.5 mm along a major axis. Typical translations of the prostate were between 0-4 mm. The translation and rotation had a predominant direction, suggesting a natural axis for prostate movement. CONCLUSION Although significant prostate displacement can occur between treatments, the typical range of movement seen along a major axis was less than 5 mm. Proper treatment planning should consider the movement of the target independent of surrounding bony anatomy. Advances in online portal imaging, image registration, and dynamic field shaping may permit shaped fields that encompass the prostate gland in its position at the time of treatment, allowing for the use of smaller fields while ensuring proper target coverage.

Uncertainties in CT-based radiation therapy treatment planning associated with patient breathing

PURPOSE To evaluate uncertainties associated with treatment-planning computed tomography (CT) data obtained with the patient breathing freely. METHODS AND MATERIALS Patients with thoracic or abdominal tumors underwent a standard treatment-planning CT study while breathing quietly and freely, followed by CT scans while holding their breath at normal inhalation and normal exhalation. Identical treatment plans on all three CT data sets for each patient pointed out differences in: (a) radiation path lengths; (b) positions of the organs; (c) physical volumes of the lung, liver, and kidneys; (d) the interpretation of plan evaluation tools such as dose-volume histograms and normal tissue complication probability (NTCP) models; and (e) how well the planning CT data set represented the average of the inhalation and exhalation studies. RESULTS Inhalation and exhalation data differ in terms of radiation path length (nearly one quarter of the cases had path-length differences > 1 cm), although the free breathing and average path lengths do not exhibit large differences (0-9 mm). Liver and kidney movements averaged 2 cm, whereas differences between the free breathing and average positions averaged 0.6 cm. The physical volume of the liver between the free breathing and static studies varied by as much as 12%. The NTCP calculations on exhale and inhale studies varied from 3 to 43% for doses that resulted in a 15% NTCP on the free-breathing studies. CONCLUSION Free-breathing CT studies may improperly estimate the position and volume of critical structures, and thus may mislead evaluation of plans based on such volume dependent criteria such as dose-volume histograms and NTCP calculations.

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.

The use of 3-D dose volume analysis to predict radiation hepatitis

Although it is well known that the tolerance of the liver to external beam irradiation depends on the volume of liver irradiated, few data exist which quantify this dependence. Therefore, a review was carried out of our clinical trial for the treatment of intrahepatic malignancies in which the dose of radiation delivered depended on the volume of normal liver treated. Three dimensional treatment planning using dose-volume histogram analysis of the normal liver was used for all patients. Nine of the 79 patients treated developed clinical radiation hepatitis. None of the patient related variables assessed were associated with radiation hepatitis. All patients who developed radiation hepatitis received whole liver irradiation, as all or part of their treatment, which produced a mean dose ≥ 37 Gy. Dose volume histograms were used to calculate normal tissue complication probabilities based on parameters derived from the literature. The risk of complication was greatly overestimated among patients receiving a high dose of radiation to part of the liver without whole liver treatment. An estimation of model parameters based on the clinical results indicated a larger magnitude for the “volume effect parameter” than the literature estimate (n = 0.69 ± 0.05 vs 0.32; p < 0.001). Computation of the normal tissue complication probabilities using the larger value of n produced a good description of the observed risk of radiation hepatitis. These findings suggest that dose volume histogram analysis can be used to quantify the tolerance of the liver to radiation. The predictive value of this parameterization of the normal tissue complication probability model will need to be tested with liver tolerance and dose volume histogram data from an independent clinical trial.

Radiation-Associated Liver Injury

The liver is a critically important organ that has numerous functions including the production of bile, metabolism of ingested nutrients, elimination of many waste products, glycogen storage, and plasma protein synthesis. The liver is often incidentally irradiated during radiation therapy (RT) for tumors in the upper- abdomen, right lower lung, distal esophagus, or during whole abdomen or whole body RT. This article describes the endpoints, time-course, and dose-volume effect of radiation on the liver.

Escalated Focal Liver Radiation and Concurrent Hepatic Artery Fluorodeoxyuridine for Unresectable Intrahepatic Malignancies

PURPOSE To evaluate the response, time to progression, survival, and impact of radiation (RT) dose on survival in patients with intrahepatic malignancies treated on a phase I trial of escalated focal liver RT. PATIENTS AND METHODS From April 1996 to January 1998, 43 patients with unresectable intrahepatic hepatobiliary cancer (HB; 27 patients) and colorectal liver metastases (LM; 16 patients) were treated with high-dose conformal RT. The median tumor size was 10 x 10 x 8 cm. The median RT dose was 58.5 Gy (range, 28.5 to 90 Gy), 1.5 Gy twice daily, with concurrent continuous-infusion hepatic arterial fluorodeoxyuridine (0.2 mg/kg/d) during the first 4 weeks of RT. RESULTS The response rate in 25 assessable patients was 68% (16 partial and one complete response). With a median potential follow-up period of 26.5 months, the median times to progression for all tumors, LM, and HB were 6, 8, and 3 months, respectively. The median survival times of all patients, patients with LM, and patients with HB were 16, 18, and 11 months, respectively. On multivariate analyses, escalated RT dose was independently associated with improved progression-free and overall survival. The median survival of patients treated with 70 Gy or more has not yet been reached (16.4+ months), compared with 11.6 months in patients treated with lower RT doses (P =.0003). CONCLUSION The excellent response rate, prolonged intrahepatic control, and improved survival in patients treated with RT doses of 70 Gy or more motivate continuation of dose-escalation studies for patients with intrahepatic malignancies.