A new method to determine dose-effect relations for local lung-function changes using correlated SPECT and CT data

Quantifying Regional Radiation-Induced Lung Injury in Patients Using Hyperpolarized 129Xe Gas Exchange Magnetic Resonance Imaging

PURPOSE

Radiation-induced lung injury has been shown to alter regional ventilation and perfusion in the lung. However, changes in regional pulmonary gas exchange have not previously been measured.

METHODS AND MATERIALS

Ten patients receiving conventional radiation therapy (RT) for lung cancer underwent pre-RT and 3-month post-RT magnetic resonance imaging (MRI) using an established hyperpolarized 129Xe gas exchange technique to map lung function. Four patients underwent an additional 8-month post-RT MRI. The MR signal from inhaled xenon was measured in the following 3 pulmonary compartments: the lung airspaces, the alveolar membrane tissue, and the pulmonary capillaries (interacting with red blood cells [RBCs]). Thoracic 1H MRI scans were acquired, and deformable registration was used to transfer 129Xe functional maps to the RT planning computed tomography scan. The RT-associated changes in ventilation, membrane uptake, and RBC transfer were computed as a function of regional lung dose (equivalent dose in 2-Gy fractions). Pearson correlations and t tests were used to determine statistical significance, and weighted sum of squares linear regression subsequently characterized the dose dependence of each functional component. The pulmonary function testing metrics of forced vital capacity and diffusing capacity for carbon monoxide were also acquired at each time point.

RESULTS

Compared with pre-RT baseline, 3-month post-RT ventilation decreased by an average of -0.24 ± 0.05%/Gy (ρ = -0.88; P < .001), membrane uptake increased by 0.69 ± 0.14%/Gy (ρ = 0.94; P < .001), and RBC transfer decreased by -0.41 ± 0.06%/Gy (ρ = -0.92; P < .001). Membrane uptake maintained a strong positive correlation with regional dose at 8 months post-RT, demonstrating an increase of 0.73 ± 0.11%/Gy (ρ = 0.92; P = .006). Changes in membrane uptake and RBC transfer appeared greater in magnitude (%/Gy) for individuals with low heterogeneity in their baseline lung function. An increase in whole-lung membrane uptake showed moderate correlation with decreases in forced vital capacity (ρ = -0.50; P = .17) and diffusing capacity for carbon monoxide (ρ = -0.44; P = .23), with neither correlation reaching statistical significance.

CONCLUSIONS

Hyperpolarized 129Xe MRI measured and quantified regional, RT-associated, dose-dependent changes in pulmonary gas exchange. This tool could enable future work to improve our understanding and management of radiation-induced lung injury.

Is individual perfusion dose-response different than ventilation dose-response for lung cancer patients treated with radiotherapy?

OBJECTIVE

Current ventilation and perfusion dose-response studies focus on single-modalities (ventilation or perfusion) and perform pulmonary-toxicity assessment related to radiotherapy on a population-based basis. This study aims at quantitative and clinical evaluation of intrapatient differences between ventilation and perfusion dose-responses among lung cancer patients treated with radiotherapy.

METHODS

20 patients enrolled on a prospective functional avoidance protocol underwent single photon emission computed tomography-CT ventilation and perfusion scans pre- and post-radiotherapy. Relative changes in pre- to post-treatment ventilation and perfusion in lung regions receiving ≥20 Gy were calculated. In addition, the slopes of the linear fit to the relative ventilation and perfusion changes in regions receiving 0-60 Gy were calculated. A radiologist read and assigned a functional defect score to pre- and post-treatment ventilation/perfusion scans.

RESULTS

25% of patients had a difference >35% between ventilation and perfusion pre- to post-treatment changes and 20-30% of patients had opposite directions for ventilation and perfusion pre- to post-treatment changes. Using a semi-quantitative scale, radiologist assessment showed that 20% of patients had different pre- to post-treatment ventilation changes when compared to pre- to post-treatment perfusion changes.

CONCLUSION

Our data showed that ventilation dose-response can differ from perfusion dose-response for 20-30% of patients. Therefore, when performing thoracic dose-response in cancer patients, it is insufficient to look at ventilation or perfusion alone; but rather both modes of functional imaging may be needed when predicting for clinical outcomes.

ADVANCES IN KNOWLEDGE

The significance of this study can be highlighted by the differences between the intrapatient dose-response assessments of this analysis compared to existing population-based dose-response analyses. Elucidating intrapatient ventilation and perfusion dose-response differences may be valuable in predicting pulmonary toxicity in lung cancer patients post-radiotherapy.

Radiotherapy for primary lung cancer

Herein are presented the recommendations from the Société française de radiothérapie oncologique regarding indications and modalities of lung cancer radiotherapy. The recommendations for delineation of the target volumes and organs at risk are detailed.

Synthetic pulmonary perfusion images from 4DCT for functional avoidance using deep learning

Purpose.To develop and evaluate the performance of a deep learning model to generate synthetic pulmonary perfusion images from clinical 4DCT images for patients undergoing radiotherapy for lung cancer.Methods. A clinical data set of 58 pre- and post-radiotherapy^99mTc-labeled MAA-SPECT perfusion studies (32 patients) each with contemporaneous 4DCT studies was collected. Using the inhale and exhale phases of the 4DCT, a 3D-residual network was trained to create synthetic perfusion images utilizing the MAA-SPECT as ground truth. The training process was repeated for a 50-imaging study, five-fold validation with twenty model instances trained per fold. The highest performing model instance from each fold was selected for inference upon the eight-study test set. A manual lung segmentation was used to compute correlation metrics constrained to the voxels within the lungs. From the pre-treatment test cases (N = 5), 50th percentile contours of well-perfused lung were generated from both the clinical and synthetic perfusion images and the agreement was quantified.Results. Across the hold-out test set, our deep learning model predicted perfusion with a Spearman correlation coefficient of 0.70 (IQR: 0.61-0.76) and a Pearson correlation coefficient of 0.66 (IQR: 0.49-0.73). The agreement of the functional avoidance contour pairs was Dice of 0.803 (IQR: 0.750-0.810) and average surface distance of 5.92 mm (IQR: 5.68-7.55).Conclusion. We demonstrate that from 4DCT alone, a deep learning model can generate synthetic perfusion images with potential application in functional avoidance treatment planning.

Investigating the SPECT Dose-Function Metrics Associated With Radiation-Induced Lung Toxicity Risk in Patients With Non-small Cell Lung Cancer Undergoing Radiation Therapy

Purpose

Dose to normal lung has commonly been linked with radiation-induced lung toxicity (RILT) risk, but incorporating functional lung metrics in treatment planning may help further optimize dose delivery and reduce RILT incidence. The purpose of this study was to investigate the impact of the dose delivered to functional lung regions by analyzing perfusion (Q), ventilation (V), and combined V/Q single-photon-emission computed tomography (SPECT) dose-function metrics with regard to RILT risk in patients with non-small cell lung cancer (NSCLC) patients who received radiation therapy (RT).

Methods and Materials

SPECT images acquired from 88 patients with locally advanced NSCLC before undergoing conventionally fractionated RT were retrospectively analyzed. Dose was converted to the nominal dose equivalent per 2 Gy fraction, and SPECT intensities were normalized. Regional lung segments were defined, and the average dose delivered to each lung region was quantified. Three functional categorizations were defined to represent low-, normal-, and high-functioning lungs. The percent of functional lung category receiving ≥20 Gy and mean functional intensity receiving ≥20 Gy (iV₂₀) were calculated. RILT was defined as grade 2+ radiation pneumonitis and/or clinical radiation fibrosis. A logistic regression was used to evaluate the association between dose-function metrics and risk of RILT.

Results

By analyzing V/Q normalized intensities and functional distributions across the population, a wide range in functional capability (especially in the ipsilateral lung) was observed in patients with NSCLC before RT. Through multivariable regression models, global lung average dose to the lower lung was found to be significantly associated with RILT, and Q and V iV₂₀ were correlated with RILT when using ipsilateral lung metrics. Through a receiver operating characteristic analysis, combined V/Q low-function receiving ≥20 Gy (low-functioning V/Q₂₀) in the ipsilateral lung was found to be the best predictor (area under the curce: 0.79) of RILT risk.

Conclusions

Irradiation of the inferior lung appears to be a locational sensitivity for RILT risk. The multivariable correlation between ipsilateral lung iV₂₀ and RILT, as well as the association of low-functioning V/Q₂₀ and RILT, suggest that irradiating low-functioning regions in the lung may lead to higher toxicity rates.

Functional liver-image guided hepatic therapy (FLIGHT): A technique to maximize hepatic functional reserve

INTRODUCTION

Radiation planning approaches for liver radiation often do not consider the regional variation that can exist in liver function. This study dosimetrically compares functional liver image-guided hepatic therapy (FLIGHT) to standard stereotactic body radiation therapy (SBRT) plans. In the FLIGHT plans, functional data from hepatobiliary iminodiacetic acid (HIDA) single photon emission computed tomography (SPECT) scans serve as a road map to guide beam arrangement. While meeting the same target volume coverage, plans are optimized to reduce dose to high-functioning liver.

MATERIALS AND METHODS

The study included 10 patients with hepatocellular carcinoma (HCC) with baseline HIDA SPECT imaging. Standard SBRT plans which did not systematically incorporate these scans had previously been completed on all 10 plans. Retrospectively, FLIGHT plans were created based on the use of contours of relative liver function from the HIDA SPECT as avoidance structures. Resulting dose to each relative functional liver structure was examined and compared qualitatively and using Wilcoxin rank-sum tests. Target coverage, doses to organs at risk (OARs), conformity index (CI), and gradient index (GI) were also evaluated.

RESULTS

While maintaining the same target coverage, FLIGHT plans reduced the mean dose to the high functioning liver by a median of 3.0 Gy (range 0.7 to 4.6 Gy), which represented a 31.4% mean reduction compared to standard planning. FLIGHT plans reduced the volume of high functioning liver receiving 15 Gy by a mean of 59.3 cc (range 7 to 170 cc), for a mean reduction of 41.9%. The mean dose to areas of liver function defined by 25% to 100% and 50% to 100% maximum was reduced with FLIGHT from 10.5 Gy to 8.5 Gy and from 10.5 Gy to 7.5 Gy, respectively (p < 0.005 for both comparisons). The FLIGHT plans' mean CI and GI did not differ significantly from the standard plans' (p = 0.721 and 0.169, respectively).

CONCLUSION

FLIGHT SBRT allows for field design and plan optimization individualized to a patient's baseline regional liver function to maximize hepatic functional reserve. This personalized approach is achieved without compromising target coverage or OAR sparing.

The Role of Mesenchymal Stem Cells in Radiation-Induced Lung Fibrosis

Radiation therapy is one of the most important treatment modalities for thoracic tumors. Despite significant advances in radiation techniques, radiation-induced lung injury (RILI) still occurs in up to 30% of patients undergoing thoracic radiotherapy, and therefore remains the main dose-limiting obstacle. RILI is a potentially lethal clinical complication of radiotherapy that has 2 main stages: an acute stage defined as radiation pneumonitis, and a late stage defined as radiation-induced lung fibrosis. Patients who develop lung fibrosis have a reduced quality of life with progressive and irreversible organ malfunction. Currently, the most effective intervention for the treatment of lung fibrosis is lung transplantation, but the lack of available lungs and transplantation-related complications severely limits the success of this procedure. Over the last few decades, advances have been reported in the use of mesenchymal stem cells (MSCs) for lung tissue repair and regeneration. MSCs not only replace damaged lung epithelial cells but also promote tissue repair through the secretion of anti-inflammatory and anti-fibrotic factors. Here, we present an overview of MSC-based therapy for radiation-induced lung fibrosis, focusing in particular on the molecular mechanisms involved and describing the most recent preclinical and clinical studies carried out in the field.

Quantifying ventilation change due to radiation therapy using 4DCT Jacobian calculations

PURPOSE

Regional ventilation and its response to radiation dose can be estimated using four-dimensional computed tomography (4DCT) and image registration. This study investigated the impact of radiation therapy (RT) on ventilation and the dependence of radiation-induced ventilation change on pre-RT ventilation derived from 4DCT.

METHODS AND MATERIALS

Three 4DCT scans were acquired from each of 12 subjects: two scans before RT and one scan 3 months after RT. The 4DCT datasets were used to generate the pre-RT and post-RT ventilation maps by registering the inhale phase image to the exhale phase image and computing the Jacobian determinant of the resulting transformation. The ventilation change between pre-RT and post-RT was calculated by taking a ratio of the post-RT Jacobian map to the pre-RT Jacobian map. The voxel-wise ventilation change between pre- and post-RT was investigated as a function of dose and pre-RT ventilation.

RESULTS

Lung regions receiving over 20 Gy exhibited a significant decrease in function (3.3%, P < 0.01) compared to those receiving less than 20 Gy. When the voxels were stratified into high and low pre-RT function by thresholding the Jacobian map at 10% volume expansion (Jacobian = 1.1), high-function voxels exhibited 4.8% reduction in function for voxels receiving over 20 Gy, a significantly greater decline (P = 0.037) than the 2.4% reduction in function for low-function voxels. Ventilation decreased linearly with dose in both high-function and low-function regions. High-function regions showed a significantly larger decline in ventilation (P ≪ 0.001) as dose increased (1.4% ventilation reduction/10 Gy) compared to low-function regions (0.3% ventilation reduction/10 Gy). With further stratification of pre-RT ventilation, voxels exhibited increasing dose-dependent ventilation reduction with increasing pre-RT ventilation, with the largest pre-RT Jacobian bin (pre-RT Jacobian between 1.5 and 1.6) exhibiting a ventilation reduction of 4.8% per 10 Gy.

CONCLUSIONS

Significant ventilation reductions were measured after radiation therapy treatments, and were dependent on the dose delivered to the tissue and the pre-RT ventilation of the tissue. For a fixed radiation dose, lung tissue with high pre-RT ventilation experienced larger decreases in post-RT ventilation than lung tissue with low pre-RT ventilation.

Correlation of Regional Lung Ventilation and Gas Transfer to Red Blood Cells: Implications for Functional-Avoidance Radiation Therapy Planning

PURPOSE

To investigate the degree to which lung ventilation and gas exchange are regionally correlated, using the emerging technology of hyperpolarized (HP)-129Xe magnetic resonance imaging (MRI).

METHODS AND MATERIALS

Hyperpolarized-129Xe MRI studies were performed on 17 institutional review board-approved human subjects, including 13 healthy volunteers, 1 emphysema patient, and 3 non-small cell lung cancer patients imaged before and approximately 11 weeks after radiation therapy (RT). Subjects inhaled 1 L of HP-129Xe mixture, followed by the acquisition of interleaved ventilation and gas exchange images, from which maps were obtained of the relative HP-129Xe distribution in three states: (1) gaseous, in lung airspaces; (2) dissolved interstitially, in alveolar barrier tissue; and (3) transferred to red blood cells (RBCs), in the capillary vasculature. The relative spatial distributions of HP-129Xe in airspaces (regional ventilation) and RBCs (regional gas transfer) were compared. Further, we investigated the degree to which ventilation and RBC transfer images identified similar functional regions of interest (ROIs) suitable for functionally guided RT. For the RT patients, both ventilation and RBC functional images were used to calculate differences in the lung dose-function histogram and functional effective uniform dose.

RESULTS

The correlation of ventilation and RBC transfer was ρ = 0.39 ± 0.15 in healthy volunteers. For the RT patients, this correlation was ρ = 0.53 ± 0.02 before treatment and ρ = 0.39 ± 0.07 after treatment; for the emphysema patient it was ρ = 0.24. Comparing functional ROIs, ventilation and RBC transfer demonstrated poor spatial agreement: Dice similarity coefficient = 0.50 ± 0.07 and 0.26 ± 0.12 for the highest-33%- and highest-10%-function ROIs in healthy volunteers, and in RT patients (before treatment) these were 0.58 ± 0.04 and 0.40 ± 0.04. The average magnitude of the differences between RBC- and ventilation-derived functional effective uniform dose, fV20Gy, fV10Gy, and fV5Gy were 1.5 ± 1.4 Gy, 4.1% ± 3.8%, 5.0% ± 3.8%, and 5.3% ± 3.9%, respectively.

CONCLUSION

Ventilation may not be an effective surrogate for true regional lung function for all patients.

Functional liver image guided hepatic therapy (FLIGHT) with hepatobiliary iminodiacetic acid (HIDA) scans

PURPOSE

Hepatobiliary iminodiacetic acid (HIDA) scans provide global and regional assessments of liver function that can serve as a road map for functional avoidance in stereotactic body radiation therapy (SBRT) planning. Functional liver image guided hepatic therapy (FLIGHT), an innovative planning technique, is described and compared with standard planning using functional dose-volume histograms. Thresholds predicting for decompensation during follow up are evaluated.

METHODS AND MATERIALS

We studied 17 patients who underwent HIDA scans before SBRT. All SBRT cases were replanned using FLIGHT. The following dosimetric endpoints were compared for FLIGHT versus standard SBRT planning: functional residual capacity <15 Gy (FRC₁₅HIDA), mean liver dose (MLD), equivalent uniform dose (EUD), and functional EUD (FEUD). Receiver operating characteristics curves were used to evaluate whether baseline HIDA values, standard cirrhosis scoring, and/or dosimetric data predicted clinical decompensation.

RESULTS

Compared with standard planning, FLIGHT significantly improved FRC₁₅HIDA (mean improvement: 5.3%) as well as MLD, EUD, and FEUD (P < .05). Considerable interindividual variations in the extent of benefit were noted. Decompensation during follow-up was associated with baseline global HIDA <2.915%/min/m², FRC₁₅HIDA <2.11%/min/m², and MELD ≥11 (P < .05).

CONCLUSIONS

FLIGHT with HIDA-based parameters may complement blood chemistry-based assessments of liver function and facilitate individualized, adaptive liver SBRT planning.

Modeling Patient-Specific Dose-Function Response for Enhanced Characterization of Personalized Functional Damage

PURPOSE

Functional-guided radiation therapy (RT) plans have the potential to limit damage to normal tissue and reduce toxicity. Although functional imaging modalities have continued to improve, a limited understanding of the functional response to radiation and its application to personalized therapy has hindered clinical implementation. The purpose of this study was to retrospectively model the longitudinal, patient-specific dose-function response in non-small cell lung cancer patients treated with RT to better characterize the expected functional damage in future, unknown patients.

METHODS AND MATERIALS

Perfusion single-photon emission computed tomography/computed tomography scans were obtained at baseline (n = 81), midtreatment (n = 74), 3 months post-treatment (n = 51), and 1 year post-treatment (n = 26) and retrospectively analyzed. Patients were treated with conventionally fractionated RT or stereotactic body RT. Normalized perfusion single-photon emission computed tomography voxel intensity was used as a surrogate for local lung function. A patient-specific logistic model was applied to each individual patient's dose-function response to characterize functional reduction at each imaging time point. Patient-specific model parameters were averaged to create a population-level logistic dose-response model.

RESULTS

A significant longitudinal decrease in lung function was observed after RT by analyzing the voxelwise change in normalized perfusion intensity. Generated dose-function response models represent the expected voxelwise reduction in function, and the associated uncertainty, for an unknown patient receiving conventionally fractionated RT or stereotactic body RT. Differential treatment responses based on the functional status of the voxel at baseline suggest that initially higher functioning voxels are damaged at a higher rate than lower functioning voxels.

CONCLUSIONS

This study modeled the patient-specific dose-function response in patients with non-small cell lung cancer during and after radiation treatment. The generated population-level dose-function response models were derived from individual patient assessment and have the potential to inform functional-guided treatment plans regarding the expected functional lung damage. This type of patient-specific modeling approach can be applied broadly to other functional response analyses to better capture intrapatient dependencies and characterize personalized functional damage.

Pulmonary injury associated with radiation therapy - Assessment, complications and therapeutic targets

Pulmonary injury is more common in patients undergoing radiation therapy for lungs and other thoracic malignancies. Recently with the use of most-advanced technologies powerful doses of radiation can be delivered directly to tumor site with exquisite precision. The awareness of technical and clinical parameters that influence the chance of radiation induced lung injury is important to guide patient selection and toxicity minimization strategies. At the cellular level, radiation activates free radical production, leading to DNA damage, apoptosis, cell cycle changes, and reduced cell survival. Preclinical research shows the potential for therapies targeting transforming growth factor-β (TGF-B), Toll like receptor (TLRs), Tumour necrosis factor-alpha (TNF-alpha), Interferon gamma (IFN-γ) and so on that may restore lung function. At present Amifostine (WR-2721) is the only approved broad spectrum radioprotector in use for patients undergoing radiation therapy. Newer techniques also offer the opportunity to identify new biomarkers and new targets for interventions to prevent or ameliorate these late effects of lung damage.

Decreased Lung Perfusion After Breast/Chest Wall Irradiation: Quantitative Results From a Prospective Clinical Trial

PURPOSE

To quantify lung perfusion changes after breast/chest wall radiation therapy (RT) using pre- and post-RT single photon emission computed tomography/computed tomography (SPECT/CT) attenuation-corrected perfusion scans; and correlate decreased perfusion with adjuvant RT dose for breast cancer in a prospective clinical trial.

METHODS AND MATERIALS

As part of an institutional review board-approved trial studying the impact of RT technique on lung function in node-positive breast cancer, patients received breast/chest wall and regional nodal irradiation including superior internal mammary node RT to 50 to 52.2 Gy with a boost to the tumor bed/mastectomy scar. All patients underwent quantitative SPECT/CT lung perfusion scanning before RT and 1 year after RT. The SPECT/CT scans were co-registered, and the ratio of decreased perfusion after RT relative to the pre-RT perfusion scan was calculated to allow for direct comparison of SPECT/CT perfusion changes with delivered RT dose. The average ratio of decreased perfusion was calculated in 10-Gy dose increments from 0 to 60 Gy.

RESULTS

Fifty patients had complete lung SPECT/CT perfusion data available. No patient developed symptoms consistent with pulmonary toxicity. Nearly all patients demonstrated decreased perfusion in the left lung according to voxel-based analyses. The average ratio of lung perfusion deficits increased for each 10-Gy increment in radiation dose to the lung, with the largest changes in regions of lung that received 50 to 60 Gy (ratio 0.72 [95% confidence interval 0.64-0.79], P<.001) compared with the 0- to 10-Gy region. For each increase in 10 Gy to the left lung, the lung perfusion ratio decreased by 0.06 (P<.001).

CONCLUSIONS

In the assessment of 50 patients with node-positive breast cancer treated with RT in a prospective clinical trial, decreased lung perfusion by SPECT/CT was demonstrated. Our study allowed for quantification of lung perfusion defects in a prospective cohort of breast cancer patients for whom attenuation-corrected SPECT/CT scans could be registered directly to RT treatment fields for precise dose estimates.

Loss of lung function after chemo-radiotherapy for NSCLC measured by perfusion SPECT/CT: Correlation with radiation dose and clinical morbidity

BACKGROUND

The purpose of the study was to assess dose and time dependence of radiotherapy (RT)-induced changes in regional lung function measured with single photon emission computed tomography (SPECT) of the lung and relate these changes to the symptomatic endpoint of radiation pneumonitis (RP) in patients treated for non-small cell lung cancer (NSCLC).

MATERIAL AND METHODS

NSCLC patients scheduled to receive curative RT of minimum 60 Gy were included prospectively in the study. Lung perfusion SPECT/CT was performed before and three months after RT. Reconstructed SPECT/CT data were registered to treatment planning CT. Dose to the lung was segmented into regions corresponding to 0-5, 6-20, 21-40, 41-60 and > 60 Gy. Changes (%) in regional lung perfusion before and after RT were correlated with regional dose and symptomatic RP (CTC grade 2-5) outcome.

RESULTS

A total of 58 patients were included, of which 45 had three-month follow-up SPECT/CT scans. Analysis showed a statistically significant dose-dependent reduction in regional perfusion at three-month follow-up. The largest population composite perfusion loss was in 41-60 Gy (42.2%) and > 60 Gy (41.7%) dose bins. Lung regions receiving low dose of 0-5 Gy and 6-20 Gy had corresponding perfusion increase (-7.2% and -6.1%, respectively). Regional perfusion reduction was different in patients with and without RP with the largest difference in 21-40 Gy bin (p = 0.02), while for other bins the difference did not reach statistical significance. The risk of symptomatic RP was higher for the patients with perfusion reduction after RT (p = 0.02), with the relative risk estimate of 3.6 (95% CI 1.1-12).

CONCLUSION

Perfusion lung function changes in a dose-dependent manner after RT. The severity of radiation-induced lung symptoms is significantly correlated with SPECT perfusion changes. Perfusion reduction early after RT is associated with a high risk of later development of symptomatic RP.

Reproducibility of four-dimensional computed tomography-based lung ventilation imaging

RATIONALE AND OBJECTIVES

A novel ventilation imaging method based on four-dimensional (4D) computed tomography (CT) has been applied to the field of radiation oncology. Understanding its reproducibility is a prerequisite for clinical applications. The purpose of this study was to quantify the reproducibility of 4D CT ventilation imaging over different days and the same session.

MATERIALS AND METHODS

Two ventilation images were created from repeat 4D CT scans acquired over the average time frames of 15 days for 6 lung cancer patients and 5 minutes for another 6 patients. The reproducibility was quantified using the voxel-based Spearman rank correlation coefficients for all lung voxels and Dice similarity coefficients (DSC) for the spatial overlap of segmented high-, moderate-, and low-functional lung volumes. Furthermore, the relationship between the variation in abdominal motion range as a measure of the depth of breathing and variation in ventilation was evaluated using linear regression.

RESULTS

The voxel-based correlation between the two ventilation images was moderate on average (0.50 ± 0.15). The DSCs were also moderate for the high- (0.60 ± 0.08), moderate- (0.46 ± 0.06), and low-functional lung (0.58 ± 0.09). No patients demonstrated strong correlations. The relationship between the motion range variation and ventilation variation was found to be moderate and significant.

CONCLUSIONS

We investigated the reproducibility of 4D CT ventilation imaging over the time frames of 15 days and 5 minutes and found that it was only moderately reproducible. Respiratory variation during 4D CT scans was found to deteriorate the reproducibility. Improvement of 4D CT imaging is necessary to increase the reproducibility of 4D CT ventilation imaging.

Functional dosimetric metrics for predicting radiation-induced lung injury in non-small cell lung cancer patients treated with chemoradiotherapy

Background

Radiation-induced lung injury (RILI) is an important dose-limiting toxicity during thoracic radiotherapy. The purpose of this study is to investigate single photon emission computed tomography (SPECT) perfusion-weighted functional dose-volume histogram (FDVH) for predicting RILI in non-small cell lung cancer (NSCLC) patients treated with definitive chemoradiotherapy.

Methods

Fifty-seven locally advanced NSCLC patients receiving chemoradiotherapy were enrolled prospectively. Patients had treatment scans and dose calculations to provide a standard dose-volume histogram (DVH). Fusion of SPECT and computed tomography scans provided perfusion-weighted FDVH and associated functional dosimetric parameters (relative volumes of functional lung receiving more than a threshold dose of 5 – 60 Gy at increments of 5 Gy [FV5 – FV60]). The predictive abilities of FDVH and DVH were calculated and compared based on the area under receiver operating characteristic (ROC) curve (AUC).

Results

The accumulative incidence of ≥ 2 grade RILI was 19.3% with a median follow-up of 12 months. Univariate analysis showed that the functional (FV5 – FV60) and standard (V5 – V40) parameters were associated with RILI (all value of p < 0.05). Close correlations between a variety of functional and standard parameters were found. By ROC curve analysis, functional metrics (AUCs were 0.784 – 0.869) provided similarly (p value 0.233 – 1.000) predictive outcome to standard metrics (AUCs were 0.695 – 0.902) in lower – median dose level parameters (FV5 – FV40). However, FDVH seemed to add some predictive value in higher dose level, the best statistical significance for comparing FV60 with V60 was 0.693 vs. 0.511 (p = 0.055).

Conclusions

Functional metrics are identified as reliable predictors for RILI, however, this observation still needs to be further verified using a larger sample size.

ICRP publication 118: ICRP statement on tissue reactions and early and late effects of radiation in normal tissues and organs--threshold doses for tissue reactions in a radiation protection context

This report provides a review of early and late effects of radiation in normal tissues and organs with respect to radiation protection. It was instigated following a recommendation in Publication 103 (ICRP, 2007), and it provides updated estimates of 'practical' threshold doses for tissue injury defined at the level of 1% incidence. Estimates are given for morbidity and mortality endpoints in all organ systems following acute, fractionated, or chronic exposure. The organ systems comprise the haematopoietic, immune, reproductive, circulatory, respiratory, musculoskeletal, endocrine, and nervous systems; the digestive and urinary tracts; the skin; and the eye. Particular attention is paid to circulatory disease and cataracts because of recent evidence of higher incidences of injury than expected after lower doses; hence, threshold doses appear to be lower than previously considered. This is largely because of the increasing incidences with increasing times after exposure. In the context of protection, it is the threshold doses for very long follow-up times that are the most relevant for workers and the public; for example, the atomic bomb survivors with 40-50years of follow-up. Radiotherapy data generally apply for shorter follow-up times because of competing causes of death in cancer patients, and hence the risks of radiation-induced circulatory disease at those earlier times are lower. A variety of biological response modifiers have been used to help reduce late reactions in many tissues. These include antioxidants, radical scavengers, inhibitors of apoptosis, anti-inflammatory drugs, angiotensin-converting enzyme inhibitors, growth factors, and cytokines. In many cases, these give dose modification factors of 1.1-1.2, and in a few cases 1.5-2, indicating the potential for increasing threshold doses in known exposure cases. In contrast, there are agents that enhance radiation responses, notably other cytotoxic agents such as antimetabolites, alkylating agents, anti-angiogenic drugs, and antibiotics, as well as genetic and comorbidity factors. Most tissues show a sparing effect of dose fractionation, so that total doses for a given endpoint are higher if the dose is fractionated rather than when given as a single dose. However, for reactions manifesting very late after low total doses, particularly for cataracts and circulatory disease, it appears that the rate of dose delivery does not modify the low incidence. This implies that the injury in these cases and at these low dose levels is caused by single-hit irreparable-type events. For these two tissues, a threshold dose of 0.5Gy is proposed herein for practical purposes, irrespective of the rate of dose delivery, and future studies may elucidate this judgement further.

Functional dose-volume histograms for predicting radiation pneumonitis in locally advanced non-small cell lung cancer treated with late-course accelerated hyperfractionated radiotherapy

The aim of this study was to determine whether functional dose-volume histograms (FDVHs) are valuable for predicting radiation pneumonitis (RP), and to identify whether FDVHs have advantages over conventional dose-volume histograms (DVHs) for the prediction of RP in patients with locally advanced non-small cell lung cancer (LANSCLC). Fifty-seven patients with LANSCLC undergoing functional image-guided late-course accelerated hyperfractionated radiotherapy were enrolled. The grade of RP was evaluated according to the Common Toxicity Criteria 3.0. To identify predictive factors of RP, the FDVHs, including the volume of the functional lung receiving 5 Gy (FV(5)) through 50 Gy (FV(50)), mean perfusion-weighted lung dose (MPWLD) and functional normal tissue complication probability (FNTCP), were analyzed and compared to their counterparts [total lung receiving 5 Gy (V(5)) through 50 Gy (V(50)), mean lung dose (MLD) and normal tissue complication probability (NTCP)] derived from conventional DVHs. Univariate analysis revealed that V(5)-V(40), MLD, NTCP and FV(5)-FV(50), MPWLD, FNTCP were all statistically significant relative to the development of RP (all p<0.05). Multivariate analysis showed that only MLD and FV(15) were associated with RP (p=0.001 and 0.044, respectively). Receiver operator characteristic curve anaysis indicated that almost all of the FDVHs had larger areas under the curve compared to the DVHs, although no statistically significant difference was observed (p-value ranged from 0.066 to 0.951). FDVHs are valuable for predicting RP with the predictive efficiency equivalent to or slightly advantageous over conventional DVHs. More homogeneous studies involving larger numbers of patients are required to further assess the value of FDVHs for predicting RP.

Lung perfusion imaging can risk stratify lung cancer patients for the development of pulmonary complications after chemoradiation

INTRODUCTION

We investigated the value of lung perfusion imaging in predicting the risk of developing pulmonary complications after chemoradiation (CRT) or radiation therapy (RT) for lung cancer.

METHODS

Fifty patients who underwent lung perfusion imaging before RT for lung cancer were included. Planar and single photon emission computed tomography/computed tomography images of the lungs were obtained. Lung perfusion score (LPS) was developed to visually grade localized perfusion defect per lung on a scale of 0 to 4 and perfusion pattern in the remaining lungs on a scale of 1 to 4. The LPS is the sum of the score for the localized perfusion defect in each lung plus the score for the remaining lungs perfusion. LPSs were correlated with pulmonary function tests and the patients were followed for 8 months after therapy to determine the incidence of grade 2 to 5 symptomatic therapy related pulmonary complications according to the common terminology criteria for adverse events (CTCAE 3.0).

RESULTS

Thirty-four patients underwent CRT and 16 underwent RT. The mean total radiation dose delivered was 56.1 +/- 10.4 Gy. Eighteen patients (36%) suffered from pulmonary complications at a mean interval of 3.4 months after therapy. Nine patients had grade 2, 7 had grade 3, 1 had grade 4, and 1 had grade 5 pulmonary complications. The mean LPS was 4.9 in patients who developed pulmonary complications versus 3.5 in patients who did not (p = 0.01). There were no significant difference between pulmonary function tests in the patients with pulmonary complications and the patient without. In addition, there were no significant differences between the mean lung radiation dose, the volume of lung irradiated or the percentage of lung receiving greater than 20 Gy between the two groups.

CONCLUSIONS

LPS using lung perfusion imaging is useful for predicting possible pulmonary complications after CRT or RT in lung cancer patients.

Liver function after irradiation based on computed tomographic portal vein perfusion imaging

PURPOSE

To determine whether individual and regional liver sensitivity to radiation could be assessed by measuring liver perfusion during a course of treatment using dynamic contrast-enhanced computed tomography scanning.

METHODS AND MATERIALS

Patients with intrahepatic cancer undergoing conformal radiotherapy underwent dynamic contrast-enhanced computed tomography (to measure perfusion distribution) and an indocyanine extraction study (to measure liver function) before, during, and 1 month after treatment. We hoped to determine whether the residual functioning liver (i.e., those regions showing portal vein perfusion) could be used to predict overall liver function after irradiation.

RESULTS

Radiation doses from 45 to 84 Gy resulted in undetectable regional portal vein perfusion 1 month after treatment. The volume of each liver with undetectable portal vein perfusion ranged from 0 to 39% and depended both on the patient's sensitivity and on dose distribution. There was a significant correlation between indocyanine green clearance and the mean of the estimated portal vein perfusion in the functional liver parenchyma (p < 0.001).

CONCLUSION

This study reveals substantial individual variability in the sensitivity of the liver to irradiation. In addition, these findings suggest that hepatic perfusion imaging may be a marker for liver function and has the potential to be a tool for individualizing therapy.

The prediction of radiation-induced liver dysfunction using a local dose and regional venous perfusion model

We have shown that high dose conformal radiation combined with chemotherapy appears to prolong the survival of patients with unresectable intrahepatic cancers. The ability to safely deliver higher doses is primarily limited by the development of radiation-induced liver disease, characterized by venous occlusion. In this study, we investigated whether portal venous perfusion measured prior to the end of radiation therapy (RT) together with dose could predict liver venous perfusion dysfunction after treatment. Ten patients with unresectable intrahepatic cancer participated in an IRB-approved computer tomography (CT) perfusion study. Hepatic arterial and portal vein perfusion distributions were estimated by using dynamic contrast enhanced CT and the single compartmental model. Scans were obtained at four time points: prior to treatment, after 15 and 30 fractions of 1.5 Gy treatments, and one month following the completion of RT. Multivariant linear regression was used to determine covariances among the first three time point measurements plus dose for prediction of the post RT measurement. The reduction in the regional venous perfusion one month following RT was predicted by the local accumulated dose and the change in the regional venous perfusion after -30 fractions (F=90.6,p <0.000 01). Each Gy produced an approximately 1.2% of reduction in the venous perfusion. This local dose and venous perfusion model has the potential to predict individual sensitivity to radiation. This is the first step toward developing a method to deliver higher and potentially more curative radiation doses to the patients who can safely receive these higher doses.

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

PURPOSE

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

METHODS AND MATERIALS

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

RESULTS

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

CONCLUSION

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

Individual radiation response of parotid glands investigated by dynamic 11C-methionine PET

BACKGROUND AND PURPOSE

Previously, we showed that the net metabolic clearance of 11C-methionine of the parotid gland, K, calculated from dynamic 11C-methionine PET, can be used as a measure of parotid gland function. The aim of this study was to investigate by dynamic 11C-methionine PET the individual radiation dose response relationship of parotid glands in head and neck cancer patients.

PATIENTS AND METHODS

Twelve head and neck cancer patients were examined by dynamic 11C-methionine PET after radiotherapy. Parametric images of K were generated, co-registered and compared voxel-by-voxel with the 3D radiation dose plan within the parotid gland to assess the individual radiation dose-function relationship.

RESULTS

In each patient, voxel-values of K decreased with increasing radiation dose. Population based analysis showed a sigmoid dose response relationship of parotid gland, from which we estimated a threshold radiation dose of 16 Gy and a mean TD50 of 30 Gy. TD50 ranged from 7 to 50 Gy in the group of patients.

CONCLUSIONS

Individual radiation dose response of parotid glands can be measured by dynamic 11C-methionine PET. The dose response analysis revealed a sigmoid relationship, a threshold radiation dose of 16 Gy, and a mean TD50 of 30 Gy.

Incorporation of functional imaging data in the evaluation of dose distributions using the generalized concept of equivalent uniform dose

Advances in the fields of IMRT and functional imaging have greatly increased the prospect of escalating the dose to highly active or hypoxic tumour sub-volumes and steering the dose away from highly functional critical structure regions. However, current clinical treatment planning and evaluation tools assume homogeneous activity/function status in the tumour/critical structures. A method was developed to incorporate tumour/critical structure heterogeneous functionality in the generalized concept of equivalent uniform dose (EUD). The tumour and critical structures functional EUD (FEUD) values were calculated from the dose-function histogram (DFH), which relates dose to the fraction of total function value at that dose. The DFH incorporates flouro-deoxyglucose positron emission tomography (FDG-PET) functional data for tumour, which describes the distribution of metabolically active tumour clonogens, and single photon emission computed tomography (SPECT) perfusion data for critical structures. To demonstrate the utility of the method, the lung dose distributions of two non-small cell lung cancer patients, who received 3D conformal external beam radiotherapy treatment with curative intent, were evaluated. Differences between the calculated lungs EUD and FEUD values of up to 50% were observed in the 3D conformal plans. In addition, a non-small cell lung cancer patient was inversely planned with a target dose prescription of 76 Gy. Two IMRT plans (plan-A and plan-B) were generated for the patient based on the CT, FDG-PET and SPECT treatment planning images using dose-volume objective functions. The IMRT plans were generated with the goal of achieving more critical structures sparing in plan-B than plan-A. Results show the target volume EUD in plan-B is lower than plan-A by 5% with a value of 73.31 Gy, and the FEUD in plan-B is lower than plan-A by 2.6% with a value of 75.77 Gy. The FEUD plan-B values for heart and lungs were lower than plan-A by 22% and 18%, respectively. While EUD values show plan-A is marginally better than plan-B in terms of target volumetric coverage, the FEUD plan-B values show adequate target function coverage with significant critical structure function sparing. In conclusion, incorporating functional data in the calculation of EUD is important in evaluating the biological merit of treatment plans.

Advanced radiation treatment planning and delivery approaches for treatment of lung cancer

Great technologic progress has been made in the last decade in the radiation treatment of lung cancer. In particular, the widespread use of 3D conformal therapy has the potential to escalate the dose to the tumor while sparing dose to normal tissue. Current technology, however, has yet to impact local control and survival. It could be hypothesized that this is due to geographic misses because of poor target definition, movement of the tumor due to respiration, and dose/ fractionation levels. Several emerging technologies that are described in this article have the potential to address these problems, with results expected in the near future. The technical delivery of radiation has not reached its limit.

Incorporating an improved dose-calculation algorithm in conformal radiotherapy of lung cancer: re-evaluation of dose in normal lung tissue

BACKGROUND AND PURPOSE

The low density of lung tissue causes a reduced attenuation of photons and an increased range of secondary electrons, which is inaccurately predicted by the algorithms incorporated in some commonly available treatment planning systems (TPSs). This study evaluates the differences in dose in normal lung tissue computed using a simple and a more correct algorithm. We also studied the consequences of these differences on the dose-effect relations for radiation-induced lung injury.

MATERIALS AND METHODS

The treatment plans of 68 lung cancer patients initially produced in a TPS using a calculation model that incorporates the equivalent-path length (EPL) inhomogeneity-correction algorithm, were recalculated in a TPS with the convolution-superposition (CS) algorithm. The higher accuracy of the CS algorithm is well-established. Dose distributions in lung were compared using isodoses, dose-volume histograms (DVHs), the mean lung dose (MLD) and the percentage of lung receiving >20 Gy (V20). Published dose-effect relations for local perfusion changes and radiation pneumonitis were re-evaluated.

RESULTS

Evaluation of isodoses showed a consistent overestimation of the dose at the lung/tumor boundary by the EPL algorithm of about 10%. This overprediction of dose was also reflected in a consistent shift of the EPL DVHs for the lungs towards higher doses. The MLD, as determined by the EPL and CS algorithm, differed on average by 17+/-4.5% (+/-1SD). For V20, the average difference was 12+/-5.7% (+/-1SD). For both parameters, a strong correlation was found between the EPL and CS algorithms yielding a straightforward conversion procedure. Re-evaluation of the dose-effect relations showed that lung complications occur at a 12-14% lower dose. The values of the TD(50) parameter for local perfusion reduction and radiation pneumonitis changed from 60.5 and 34.1 Gy to 51.1 and 29.2 Gy, respectively.

CONCLUSIONS

A simple tissue inhomogeneity-correction algorithm like the EPL overestimates the dose to normal lung tissue. Dosimetric parameters for lung injury (e.g. MLD, V20) computed using both algorithms are strongly correlated making an easy conversion feasible. Dose-effect relations should be refitted when more accurate dose data is available.

Methodological issues in radiation dose-volume outcome analyses: summary of a joint AAPM/NIH workshop

This report represents a summary of presentations at a joint workshop of the National Institutes of Health and the American Association of Physicists in Medicine (AAPM). Current methodological issues in dose-volume modeling are addressed here from several different perspectives. Areas of emphasis include (a) basic modeling issues including the equivalent uniform dose framework and the bootstrap method, (b) issues in the valid use of statistics, including the need for meta-analysis, (c) issues in dealing with organ deformation and its effects on treatment response, (d) evidence for volume effects for rectal complications, (e) the use of volume effect data in liver and lung as a basis for dose escalation studies, and (f) implications of uncertainties in volume effect knowledge on optimized treatment planning. Taken together, these approaches to studying volume effects describe many implications for the development and use of this information in radiation oncology practice. Areas of significant interest for further research include the meta-analysis of clinical data; interinstitutional pooled data analyses of volume effects; analyses of the uncertainties in outcome prediction models, minimal parameter number outcome models for ranking treatment plans (e.g., equivalent uniform dose); incorporation of the effect of motion in the outcome prediction; dose-escalation/isorisk protocols based on outcome models; the use of functional imaging to study radioresponse; and the need for further small animal tumor control probability/normal tissue complication probability studies.

Update on Results of Multifield Conformal Radiation Therapy of Non—Small-Cell Lung Cancer Using Multileaf Collimated Beams

We evaluated the treatment outcome for 5-field 3-dimensional conformal radiation therapy (3D-CRT) in 46 consecutive patients with unresectable, nonmetastatic non-small-cell lung cancer treated from 1993-2001. Four percent of the patients had stage I tumors, 6% had stage II, 44% had stage IIIA, and 46% had stage IIIB tumors. The median radiation therapy (RT) dose to the gross tumor volume with a median of 467.5 cc (range, 75.0-3073.0 cc) was 6120 cGy (range, 3000-6840 cGy). Thirty-one of 46 patients (67.4%) received combined chemoradiotherapy. Mean follow-up was 13.2 months (range, 3-159 months). Survival for stage III patients was 48.7% +/-9.1% at 1 year and 25.0% +/-8.4% at 2 years, with a median survival of 12.0 months+/-4.4 months. The local control rate for stage III patients was 66.8%+/- 9.4% at 1 year and 28.5%+/- 10.4% at 2 years. Patients who received chemotherapy had better survival (P = 0.0533) and local control (P = 0.0984) compared with patients receiving RT alone. Esophageal toxicity >or= grade 3 was significantly greater in combined chemoradiotherapy patients (29% early, 13% late) compared to the patients receiving RT alone (0% early and late). Pulmonary toxicity (early and late) was limited to grades 1/2 in 24% of patients and early grade 3 in 2% of patients. Chemotherapy appears to improve survival and local control when added to 3D-CRT in this series. The addition of concurrent chemotherapy to RT significantly increased esophageal toxicity (within acceptable levels) and did not effect pulmonary toxicity in this series.

The time course of radiation therapy-induced reductions in regional perfusion: a prospective study with >5 years of follow-up

PURPOSE

To assess the time-dependence of radiation therapy (RT)-induced reductions in regional lung perfusion, as measured by single photon emission computed tomography (SPECT) lung perfusion scans.

METHODS AND MATERIALS

Between 1991 and 1999, 79 patients had SPECT lung perfusion scans before and serially after RT. Changes in regional perfusion were correlated with regional dose using 3D planning tools and image fusion (PLUNC-Plan UNC). Multiple post-RT follow-up scans were evaluated to determine the temporal nature of RT-induced regional perfusion changes. To facilitate the comparison of dose-response curves (DRCs) at different post-RT intervals, each DRC was fit to a linear model and thus described by its slope.

RESULTS

There was a dose-dependent reduction in regional perfusion at nearly all time points post-RT (p = 0.0001). The slope of the DRCs for RT-induced reductions in regional perfusion became steeper at essentially each successive follow-up interval (p = 0.0001). However, the increases in slope became progressively smaller at later follow-up intervals. Overall, about 80% of the long-term RT-induced regional perfusion injury was manifest within 12 months post-RT.

CONCLUSION

There is a progression of RT-induced reductions in regional perfusion, with most of this injury manifest within 12 months post-RT. Additional regional injury appears to evolve for years.

Can we predict radiation-induced changes in pulmonary function based on the sum of predicted regional dysfunction?

PURPOSE

To determine whether changes in whole-lung pulmonary function test (PFT) values are related to the sum of predicted radiation therapy (RT)-induced changes in regional lung perfusion.

PATIENTS AND METHODS

Between 1991 and 1998, 96 patients (61% with lung cancer) who were receiving incidental partial lung irradiation were studied prospectively. The patients were assessed with pre- and post-RT PFTs (forced expiratory volume in one second [FEV1] and diffusion capacity for carbon monoxide [DLCO]) for at least a 6-month follow-up period, and patients were excluded if it was determined that intrathoracic recurrence had an impact on lung function. The maximal declines in PFT values were noted. A dose-response model based on RT-induced reduction in regional perfusion (function) was used to predict regional dysfunction. The predicted decline in pulmonary function was calculated as the weighted sum of the predicted regional injuries: equation [see text] where Vd is the volume of lung irradiated to dose d, and Rd is the reduction in regional perfusion anticipated at dose d.

RESULTS

The relationship between the predicted and measured reduction in PFT values was significant for uncorrected DLCO (P = .005) and borderline significant for DLCO (P = .06) and FEV1 (P = .08). However, the correlation coefficients were small (range,.18 to.30). In patients with lung cancer, the correlation coefficients improved as the number of follow-up evaluations increased (range,.43 to.60), especially when patients with hypoperfusion in the lung adjacent to a central mediastinal/hilar thoracic mass were excluded (range,.59 to.91).

CONCLUSION

The sum of predicted RT-induced changes in regional perfusion is related to RT-induced changes in pulmonary function. In many patients, however, the percentage of variation explained is small, which renders accurate predictions difficult.

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.

Changes in local pulmonary injury up to 48 months after irradiation for lymphoma and breast cancer

PURPOSE

To assess the recovery from early local pulmonary injury after irradiation and to determine whether regional differences exist.

METHODS

For 110 patients treated for breast cancer or malignant lymphoma, single photon emission computed tomography (SPECT) perfusion and ventilation scans and CT scans were made before, 3, 18, and 48 months after radiotherapy. Dose-effect relations for changes in local perfusion, ventilation, and density were determined for each individual patient using spatially correlated SPECT and CT data sets, for each follow-up period. Average dose-effect relations for both subgroups were determined, as well as dose-effect relations for different regions.

RESULTS

In general, partial improvement of local pulmonary injury was observed between 3 and 18 months for each of the three endpoints. After 18 months, no further improvement was seen. Patients with breast cancer and malignant lymphoma showed a similar improvement (except for the perfusion parameter), which was attributed to a recovery from the early radiation response and could not be explained by contraction effects of fibrosis of lung parenchyma. No regional differences in radiosensitivity 18 months after treatment were observed, except for the dorsal versus ventral region. This difference was attributed to a gravity-related effect in the measuring procedure.

CONCLUSION

For all patients, a partial recovery from early local perfusion, ventilation, and density changes, was seen between 3 and 18 months after radiotherapy. After 18 months, local lung function did not further improve (lymphoma patients).

Radiation dose-effect relations and local recovery in perfusion for patients with non-small-cell lung cancer

PURPOSE

To determine local dose-effect relations for lung perfusion and density changes due to irradiation for patients with non-small-cell lung cancer (NSCLC) and to quantify the effect of reperfusion.

METHODS AND MATERIALS

For 25 NSCLC patients and a reference group of 81 patients with healthy lungs, registered single photon emission computed tomography (SPECT) lung perfusion and CT scans were made, before and after radiotherapy. Average dose-effect relations for perfusion and CT-density changes were calculated and compared with the dose-effect relation of the reference group. On the basis of these dose-effect relations, the post-RT perfusion was predicted for each patient and compared to the measured post-RT perfusion.

RESULTS

Well-perfused lung regions of the NSCLC patients showed the same dose-effect relation as the reference patients. By comparing predicted and measured post-treatment perfusion scans, regions of reperfusion could be determined for 18 of 25 NSCLC patients but for none of the reference patients.

CONCLUSION

Well-perfused lung tissue of patients with NSCLC behaves like healthy lung tissue with respect to radiation. The dose-effect relation for perfusion and CT density was extended for doses up to 80 Gy. Radiation damage in poorly perfused lung regions was less than predicted as a consequence of local reperfusion.

Results of multifield conformal radiation therapy of nonsmall-cell lung carcinoma using multileaf collimation beams

A five-field conformal technique with three-dimensional radiation therapy treatment planning (3-DRTP) has been shown to permit better definition of the target volume for lung cancer, while minimizing the normal tissue volume receiving greater than 50% of the target dose. In an initial study to confirm the safety of conventional doses, we used the five-field conformal 3-DRTP technique. We then used the technique in a second study, enhancing the therapeutic index in a series of 42 patients, as well as to evaluate feasibility, survival outcome, and treatment toxicity. Forty-two consecutive patients with nonsmall-cell lung carcinoma (NSCLC) were evaluated during the years 1993-1997. The median age was 60 years (range 34-80). The median radiation therapy (RT) dose to the gross tumor volume was 6,300 cGy (range 5,000-6,840 cGy) delivered over 6 to 6.5 weeks in 180-275 cGy daily fractions, 5 days per week. There were three patients who received a split course treatment of 5,500 cGy in 20 fractions, delivering 275 cGy daily with a 2-week break built into the treatment course after 10 fractions. The stages of disease were II in 2%, IIIA in 40%, IIIB in 42.9%, and recurrent disease in 14.3% of the patients. The mean tumor volume was 324.14 cc (range 88.3-773.7 cc); 57.1% of the patients received combined chemoradiotherapy, while the others were treated with radiation therapy alone. Of the 42 patients, 7 were excluded from the final analysis because of diagnosis of distant metastasis during treatment. Two of the patients had their histology reinterpreted as being other than NSCLC, 2 patients did not complete RT at the time of analysis, and 1 patient voluntarily discontinued treatment because of progressive deterioration. Median follow-up was 11.2 months (range 3-32.5 months). Survival for patients with Stage III disease was 70.2% at 1 year and 51.5% at 2 years, with median survival not yet reached. Local control for the entire series was 23.3+/-11.4% at 2 years. However, for Stage III patients, local control was 50% at 1 year and 30% at 2 years. Patients who received concurrent chemotherapy had significantly improved survival (P = 0.002) and local control (P = 0.004), compared with RT alone. Late esophageal toxicity of > or =Grade 3 occurred in 14.1+/-9.3% of patients (3 of 20) receiving combined chemoradiotherapy, but in none of the 15 patients treated with RT alone. Pulmonary toxicity limited to Grades 1-2 occurred in 6.8% of the patients, and none developed > or =Grade 3 pulmonary toxicity. Patients with locally advanced NSCLC, who commonly have tumor volumes in excess of 200 cc, presenta challenge for adequate dose delivery without significant toxicity. Our five-field conformal 3-DRTP technique, which incorporates treatment planning by dose/volume histogram (DVH) was associated with minimal toxicity and may facilitate dose escalation to the gross tumor.

Effect of radiotherapy and chemotherapy on pulmonary function after treatment for breast cancer and lymphoma: A follow-up study

PURPOSE

To determine the changes in pulmonary function tests (PFTs) 0 to 48 months after treatment for breast cancer and lymphoma.

PATIENTS AND METHODS

The alveolar volume (V(A)), vital capacity, forced expiratory volume in 1 second, and corrected transfer factor of carbon monoxide (T(L,COc)) were measured in 69 breast cancer and 41 lymphoma patients before treatment and 3, 18, and 48 months after treatment with radiotherapy alone or radiotherapy in combination with chemotherapy (mechlorethamine, vincristine, procarbazine, prednisone, doxorubicin, bleomycin, vinblastine; cyclophosphamide, epidoxorubicin, fluorouracil; cyclophosphamide, thiotepa, carboplatin; cyclophosphamide, methotrexate, fluorouracil). The three-dimensional dose distribution in the lung of each patient was converted to the mean lung dose. Statistical analysis was used to evaluate the changes in PFT values over time in relation to age, sex, smoking, chemotherapy, and the mean lung dose.

RESULTS

After an initial reduction in PFT values at 3 months, significant recovery was seen at 18 months for all patients. Thereafter, no further improvement could be demonstrated. Reductions in spirometry values and V(A) were related to the mean lung dose only (0.9% per Gy at 3 months and 0.4% per Gy mean dose at 18 months). T(L,COc) decreased 1. 1% per Gy mean dose and additionally decreased 6% when chemotherapy was given after radiotherapy. Chemotherapy administered before radiotherapy reduced baseline T(L,COc) values by 8% to 21%. All patients showed an improvement of 5% at 18 months.

CONCLUSION

On the basis of the mean lung dose and the chemotherapy regimen, the changes in PFT values can be estimated before treatment within 10% of the values actually observed in 72% to 85% of our patients with healthy lungs.

The effect of patient-specific factors on radiation-induced regional lung injury

PURPOSE

To assess the impact of patient-specific factors on radiation (RT)-induced reductions in regional lung perfusion.

METHODS

Fifty patients (32 lung carcinoma, 7 Hodgkin's disease, 9 breast carcinoma and 2 other thoracic tumors) had pre-RT and > or = 24-week post-RT single photon emission computed tomography (SPECT) perfusion images to assess the dose dependence of RT-induced reductions in regional lung perfusion. The SPECT data were analyzed using a normalized and non-normalized approach. Furthermore, two different mathematical methods were used to assess the impact of patient-specific factors on the dose-response curve (DRC). First, DRCs for different patient subgroups were generated and compared. Second, in a more formal statistical approach, individual DRCs for regional lung injury for each patient were fit to a linear-quadratic model (reduction = coefficient 1 x dose + coefficient 2 x dose2). Multiple patient-specific factors including tobacco history, pre-RT diffusion capacity to carbon monoxide (DLCO), transforming growth factor-beta (TGF-beta), chemotherapy exposure, disease type, and mean lung dose were explored in a multivariate analysis to assess their impact on the coefficients.

RESULTS

None of the variables tested had a consistent impact on the radiation sensitivity of regional lung (i.e., the slope of the DRC). In the formal statistical analysis, there was a suggestion of a slight increase in radiation sensitivity in the dose range >40 Gy for nonsmokers (vs. smokers) and in those receiving chemotherapy (vs. no chemotherapy). However, this finding was very dependent on the specific statistical and normalization method used.

CONCLUSION

Patient-specific factors do not have a dramatic effect on RT-induced reduction in regional lung perfusion. Additional studies are underway to better clarify this issue. We continue to postulate that patient-specific factors will impact on how the summation of regional injury translates into whole organ injury. Refinements in our methods to generate and compare SPECT scans are needed.

Quantitative pulmonary single photon emission computed tomography for radiotherapy applications

Pulmonary imaging using single photon emission computed tomography (SPECT) is the focus of current radiotherapy research, including dose-response analysis and three-dimensional (3D) radiation treatment planning. Improvement in the quantitative capability of SPECT may help establish its potential role in this application as well as others requiring accurate knowledge of pulmonary blood flow. The purposes of this study were to quantitatively evaluate SPECT filtered backprojection (FBP) and ordered subset-expectation maximization (OS-EM) reconstruction implementations for measuring absolute activity concentration in lung phantom experiments, and to incorporate quantitative SPECT techniques in 3D-RTP for lung cancer. Quantitative FBP (nonuniform iterative Chang attenuation compensation, scatter correction, and 3D postreconstruction Metz filtering) and OS-EM implementations were compared with a "clinical" implementation of FBP (uniform multiplicative Chang attenuation compensation and post-reconstruction von Hann filtering), for their ability to improve quantification of inactive and active spherical defects in the lungs of an anthropomorphic torso phantom. Activity concentration estimates were found to depend on many factors, such as region of interest size, scatter subtraction constant (k), postreconstruction deconvolution filtering and, in the case of OS-EM, total number of iterations. In general, reconstruction implementations incorporating compensation for nonuniform attenuation and scatter provided reduced bias relative to the clinical implementation. Potential applications to lung radiotherapy, including dose-functional histograms and treatment planning are also discussed. SPECT has the potential to provide accurate estimates of lung activity distributions that, together with improved image quality, may be useful for the study and prediction of therapeutic response.

CT and PET lung image registration and fusion in radiotherapy treatment planning using the chamfer-matching method

PURPOSE

We present a validation study of CT and PET lung image registration and fusion based on the chamfer-matching method.

METHODS AND MATERIALS

The contours of the lung surfaces from CT and PET transmission images were automatically segmented by the thresholding technique. The chamfer-matching technique was then used to register the extracted lung surfaces. Arithmetic means of distance between the two data sets of the pleural surfaces were used as the cost function. Matching was then achieved by iteratively minimizing the cost function through three-dimensional (3D) translation and rotation with an optimization method.

RESULTS

Both anatomic thoracic phantom images and clinical patient images were used to evaluate the performance of our registration system. Quantitative analysis from five patients indicates that the registration error in translation was 2-3 mm in the transverse plane, 3-4 mm in the longitudinal direction, and about 1.5 degree in rotation. Typical computing time for chamfer matching is about 1 min. The total time required to register a set of CT and PET lung images, including contour extraction, was generally less than 30 min.

CONCLUSION

We have implemented and validated the chamfer-matching method for CT and PET lung image registration and fusion. Our preliminary results show that the chamfer-matching method for CT and PET images in the lung area is feasible. The described registration system has been used to facilitate target definition and treatment planning in radiotherapy.

Incorporation of functional status into dose-volume analysis

The dose-volume histogram (DVH) has gained wide acceptance as a mechanism for reducing the voluminous data of a three-dimensional dose distribution into a two-dimensional graph. These graphs are often converted to a single figure of merit. This data reduction technique is used both for clinical treatment plan evaluation and as part of proposed systems for estimating control and complication probabilities. It has long been recognized that a major shortcoming of the DVH as an analysis tool is that all spatial information is discarded. A subtler problem, which is addressed in this work, is that the DVH also implies homogeneity of biological consequence of irradiation in what may be a functionally heterogeneous volume of tissue. An extension to the DVH, the functional dose-volume histogram, or dose-function histogram (DFH), is proposed, that explicitly includes quantitative three-dimensional functional information. The concept is illustrated by the use of SPECT imaging to assess the functional status of irradiated lung.

Prediction of overall pulmonary function loss in relation to the 3-D dose distribution for patients with breast cancer and malignant lymphoma

PURPOSE

To predict the changes in pulmonary function tests (PFTs) 3-4 months after radiotherapy based on the three-dimensional (3-D) dose distribution and taking into account patient- and treatment-related factors.

METHODS

For 81 patients with malignant lymphoma and breast cancer, PFTs (VA, VC, FEV1 and TL,COc) were performed prior to and 3-4 months after irradiation and dose-effect relations for early changes in local perfusion, ventilation and air-filled fraction were determined using correlated CT and SPECT data. The 3-D dose distribution of each patient was converted into four different dose-volume parameters, i.e. the mean dose in the lung and three overall response parameters (ORPs, which represent the average local injury over the complete lung). ORPs were determined using the dose-effect relations for early changes in local perfusion, ventilation and air-filled fraction. Correlation coefficients were calculated between these dose-volume parameters and the changes in PFTs. In addition, the impact of the variables chemotherapy (MOPP/ABV and CMF), tamoxifen, smoking, age and gender on the relation between the mean lung dose and the relative changes in PFTs following radiotherapy was studied using multiple regression analysis.

RESULTS

The mean lung dose proved to be the easiest parameter to predict the reduction in PFTs 3-4 months following radiotherapy. For all patients the relation between the mean lung dose and the changes in PFTs could be described with one regression line through the origin and a slope of 1% reduction in PFT for each increase of 1 Gy in mean lung dose. Smoking and CMF chemotherapy influenced the reduction in PFTs significantly for VA and TL,COc, respectively. Patients treated with MOPP/ABV prior to radiotherapy had lower pre-radiotherapy PFTs than other patient groups, but did not show further deterioration after radiotherapy (at 3-4 months).

CONCLUSIONS

The relative reduction in VA, VC, FEV1 and TL,COc 3-4 months after radiotherapy for breast cancer and malignant lymphoma can be estimated before radiotherapy based on the mean lung dose of each individual patient and taking into account the use of chemotherapy and smoking habits of the patient.

Regional dose response to pulmonary irradiation using a manual method

PURPOSE

To better understand the dose dependence of radiation therapy (RT)-induced changes in regional lung perfusion and tissue density, using a manual method to reduce inaccuracies that might be present in previously described automated methods.

MATERIALS AND METHODS

Patients who were to receive RT for tumors in and around the thorax, wherein portions of healthy lung would be incidentally irradiated, were prospectively studied. Changes in regional perfusion and tissue density were assessed by comparison of pre- and post-RT single photon emission computed tomography (SPECT), lung perfusion scans and computed tomography (CT) scans, respectively. The three-dimensional dose distribution was calculated on the pre-RT CT scan and correlated to the other scans via image registration. Study volumes were defined by hand and individually visualized on pre- and post-RT scans. The manually generated dose response data were compared to data generated using automated methods. The relationship between CT density and SPECT perfusion was also determined.

RESULTS

Thirteen patients with lung cancer were evaluated for changes in tissue density and 11 patients were evaluated for changes in regional perfusion at 12 months post-RT. In general, density increases with increasing regional dose, with marked changes at >60 Gy. Regional perfusion decreases with increasing regional dose. In the low dose regions, relative perfusion increases by 35% on average. Manually measured dose responses correlated well with those determined automatically. The relationship between regional perfusion and CT density indicates a wide range of perfusion over a narrow range of CT density, with markedly reduced perfusion at CT densities of > -600 and < -900 H.

CONCLUSIONS

The manually generated CT density dose response data broadly agree with data previously generated using automated methods. The manually generated perfusion dose response data are in fairly good agreement with automated data, lending credibility to the accuracy of the automated methods. Regional perfusion is markedly diminished where CT density is outside the range of normal lung tissue.

Dose-effect relations for early local pulmonary injury after irradiation for malignant lymphoma and breast cancer

PURPOSE

To quantify the influence of treatment- and patient-related factors on the severity of early local pulmonary injury and to establish whether regional differences are present for local dose-effect relations for early radiation-induced pulmonary injury.

METHODS

Forty-two patients with malignant lymphoma and 40 breast cancer patients were examined prior to and 3 months after radiotherapy. The lymphoma patients were irradiated with mantle fields to an average dose of 38 Gy and the breast cancer patients were irradiated with internal mammary node fields with or without tangential breast fields to an average dose of 50 Gy. Dose-effect relations for local perfusion, ventilation and density changes were determined using correlated single photon emission computed tomography (SPECT) and CT data. A multivariate analysis was performed to study the influence of irradiated volume, chemotherapy (CMF and MOPP/ABV), smoking, age and gender. In addition, dose-effect relations for different regions in the lung were determined.

RESULTS

A similar and almost linear increase of early functional changes as a function of radiation dose was observed for perfusion and ventilation, whereas the shape of the dose-effect relation and the magnitude of early structural changes were different for density. For the three end-points studied, regional differences in radiosensitivity could not be demonstrated. For the posterior lung region compared to the anterior lung region, however, a difference was observed, which could be attributed to a gravity-related effect in the measuring procedure. Local structural changes (density) were significantly smaller for smokers (P = 0.002) and young patients (P = 0.007), whereas the CMF chemotherapy regimen given after radiotherapy (P = 0.017) significantly increased the amount of functional changes (perfusion). The magnitude of local pulmonary changes was independent of the irradiated volume, the MOPP/ABV chemotherapy regimen and gender.

CONCLUSION

The dose-effect relations for early radiation-induced local pulmonary changes were independent of the irradiated volume, MOPP/ABV, gender and lung region. CMF, smoking and age influenced the magnitude of early pulmonary changes and should be taken into account in dose-escalation protocols.

Registration of MR and SPECT without using external fiducial markers

The aim of our work is to present, test and validate an automated registration method used for matching brain SPECT scans with corresponding MR scans. The method was applied on a data set consisting of ten brain IDEX SPECT scans and ten T1- and T2-weighted MR scans of the same subjects. Of two subjects a CT scan was also made. (Semi-) automated algorithms were used to extract the brain from the MR, CT and SPECT images. Next, a surface registration technique called chamfer matching was used to match the segmented brains. A perturbation study was performed to determine the sensitivity of the matching results to the choice of the starting values. Furthermore, the SPECT segmentation threshold was varied to study its effect on the resulting parameters and a comparison between the use of MR T1- and T2-weighted images was made. Finally, the two sets of CT scans were used to estimate the accuracy by matching MR to CT and comparing the MR-SPECT match to the SPECT-CT match. The perturbation study showed that for initial perturbations up to 6 cm the algorithm fails in less than 4% of the cases. A variation of the SPECT segmentation threshold over a realistic range (25%) caused an average variation in the optimal match of 0.28 cm vector length. When T2 is used instead of T1 the stability of the algorithm is comparable but the results are less realistic due the large deformations. Finally, a comparison of the direct SPECT-MR match and the indirect match with CT as intermediate yields a discrepancy of 0.4 cm vector length. We conclude that the accuracy of our automatic matching algorithm for SPECT and MR, in which no external markers were used, is comparable to the accuracies reported in the literature for non-automatic methods or methods based on external markers. The proposed method is efficient and insensitive to small variations in SPECT segmentation.