|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Thoracic Imaging |
1 From the Departments of Radiology (A.C.B., D.M.), Obstetrics and Gynecology (A.M.L.), and Preventive Medicine and Biometrics (G.K.G.), University of Colorado Denver, 4200 E 9th Ave, Denver, CO 80262; Pediatric Immunology Research Department, Children's Mercy Hospital, Kansas City, Mo (J.M.); Biogen, Cambridge, Mass (C.M.B.); and Division of Radiology, National Jewish Medical and Research Center, Denver, Colo (D.A.L.). Received December 29, 2006; revision requested March 1, 2007; revision received May 14; final version accepted August 13. Supported by Biogen. Address correspondence to A.C.B. (e-mail: alan.best{at}uchsc.edu).
| ABSTRACT |
|---|
|
|
|---|
Materials and Methods: Institutional review board approval and informed consent were obtained at all participating institutions. One hundred sixty-seven patients (110 men, 57 women; mean age, 63 years ± 9 [standard deviation]) with IPF were enrolled in a clinical trial. Patients underwent thin-section CT in the supine position at full inspiration at enrollment (baseline) and at 12-month follow-up. After segmentation of the lungs, mean lung attenuation (MLA), skewness, and kurtosis were measured. Extent of ground glass opacity and lung fibrosis were assessed visually. Forced vital capacity (FVC) and total lung capacity (TLC) were measured. Median duration of follow-up for mortality was 1.5 years. Univariate and multivariate survival analyses were used to determine the predictive value of baseline variables for survival.
Results: At univariate analysis, baseline variables predictive of death included TLC, fibrosis, skewness, and kurtosis. At multivariate analysis, FVC (P = .006) and fibrosis (P = .002) were predictors of short-term mortality. In 95 patients who had both baseline and follow-up CT scans, fibrosis (P = .030), MLA (P = .003), skewness (P < .001), and kurtosis (P < .001) all showed change indicating disease progression.
Conclusion: Visually determined disease extent on CT images is a strong independent predictor of mortality in IPF. Serial evaluation of quantitative CT measures can show disease progression in these patients.
© RSNA, 2008
| INTRODUCTION |
|---|
|
|
|---|
Quantitative CT indexes have been proposed as a standardized method for scoring the extent of interstitial lung disease. The validity of such indexes was previously examined in a population of patients with IPF by a multi-institutional study (2). That study showed that measurements of skewness, kurtosis, and mean lung attenuation (MLA) (computer-derived values that describe the shape of thin-section CT frequency histograms) correlate with the degree of pulmonary function abnormality. Quantitative CT measurements are attractive for use as indexes of IPF because of their noninvasive and quantitative nature and minimal requirements for user intervention. However, the pattern of change in quantitative CT measurements over time in patients with IPF has not been evaluated.
The two purposes of our study were (a) to determine which baseline variables (CT visual scores, quantitative CT indexes, or pulmonary function test results) are the best predictors of mortality during follow-up and (b) to describe serial changes in quantitative CT indexes over a 12-month period in a population with IPF. We hypothesized that CT visual scores, quantitative CT indexes, and pulmonary function test results could be directly compared for their prognostic value in a population with IPF and that serial quantitative CT measurements could be used to document progression of the disease over time.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Patients
Our study included 167 patients (110 men, 57 women; mean age, 63 years ± 9 [standard deviation]) enrolled in a double-blind placebo-controlled clinical trial of interferon β-1a (Avonex; Biogen) for treatment of IPF. Patient enrollment in the original study began in March 1998 and was completed in March 2001. In the study IPF was diagnosed by using the criteria developed by the American Thoracic Society (3). Diagnosis of IPF was confirmed in each patient by a panel of clinicians and was based on a review of clinical history, occupational and environmental exposure, pulmonary function test results, thin-section CT images of the lungs, and, if available, transbronchial or surgical lung biopsy slides. All patients were required to have had progression of IPF, defined as meeting at least one of the following criteria: (a) a greater than 10% relative decrease in total lung capacity (TLC) or FVC or a greater than 15% relative decrease in single-breath diffusing capacity of the lung for carbon monoxide; (b) a greater than 3% decrease in the resting oxygen saturation level, a 3 mm Hg increase in the resting gradient between the partial pressure of oxygen in the artery and that in the alveoli, or a 5% decrease in oxygen saturation with exercise; or (c) radiologic progression of disease (increase in reticular and/or ground glass opacity [GGO]) as assessed on chest radiographs or thin-section CT images.
Patients were excluded if they had end-stage IPF, which was defined as the presence of at least two of the following: (a) TLC less than 45% of the predicted volume, (b) hemoglobin-corrected single-breath diffusing capacity of the lung for carbon monoxide less than 25% of the predicted capacity, (c) resting gradient between the partial pressure of oxygen in the artery and that in the alveoli of more than 40 mm Hg, (d) oxygen saturation of less than 80% with exercise (walking on level ground on room air at own pace for 6 minutes), or (e) New York Heart Association class III or IV heart disease. Other exclusion criteria included (a) environmental or drug exposures likely to cause interstitial lung disease, (b) connective tissue disease, and (c) emphysema occupying more than 50% of the lung.
After the conclusion of this double-blind study, treatment assignment data were revealed. Twenty-nine patients (17.4%) had received a daily placebo, 54 patients (32.3%) had received 15 mg of interferon β-1a daily, 41 patients (24.6%) had received 30 mg of the drug daily, and 38 patients (22.8%) had received 60 mg of the drug daily. Treatment data for five patients (3.0%) were not available for analysis.
Our retrospective study used preexisting data from the prospective clinical trial, including image analyses. Images were not reinterpreted for our study. Pulmonary function test and thin-section CT histogram data were obtained at 30 participating hospitals in the United States and Canada at baseline and at 12-month follow-up. FVC and TLC were measured by using commonly accepted techniques and were expressed as a percentage of predicted performance for each patient, by using the formula developed by Crapo et al (4). Of the 167 enrolled patients, 95 had complete data sets for both baseline and follow-up. Data from these 95 patients were analyzed for assessment of serial change in visual and digital CT variables.
Patients were followed up from the time of enrollment until the end of the entire study to record mortality. Therefore, depending on the date of enrollment for each individual, the duration of follow-up ranged from 0.9 to 2.7 years (median follow-up, 1.5 years).
CT Image Analysis
Inspiratory thin-section (1-mm-thick) CT images were obtained at 2-cm intervals with the patient in the supine position. Images in Digital Imaging and Communications in Medicine, or DICOM, format were sent to a central imaging facility. Visual analysis was performed by two of three thoracic radiologists (including D.A.L., 16 years thoracic CT experience) who independently recorded the extent of pulmonary abnormality. The mean extent of GGO, reticular abnormality, honeycombing, and emphysema was scored to the nearest 5% in three zones in each lung. The upper zone was defined as at or above the aortic arch, the middle zone was defined as between the aortic arch and pulmonary veins, and the lower zone was defined as at or below the pulmonary veins. For the purposes of this study, the mean extent of lung fibrosis was calculated as the mean of the mean extent of reticular abnormality and honeycombing as scored by two observers in each of the six zones. Observers received training sessions prior to the study, and interobserver variability was evaluated and determined to be satisfactory for consistency in visual scoring. Baseline and follow-up images were scored during the trial, and observers were blinded to the temporal sequence of the examinations. Image analysis to determine skewness, kurtosis, and MLA was performed by using semiautomated image segmentation and standard image-analysis techniques (2).
Statistical Analysis
Statistical analyses were carried out by using statistical software (SAS, version 9.1.3; SAS Institute, Cary, NC). Means and standard deviations for measures of pulmonary function (FVC and TLC), quantitative CT indexes (MLA, skewness, and kurtosis), and CT visual scores (fibrosis, GGO, and emphysema) were determined at baseline and at 12 months. The magnitude of change in each of these variables was calculated, as were the number and percentage of patients with increasing or decreasing values. A paired t test was performed to evaluate the change in variables from baseline to 12 months. Mean and standard deviation values for each variable at baseline were also calculated separately for the patients who were alive and for those who were dead at the end of follow-up. The difference between the two groups was assessed with a t test after equal variance was established with the F test.
Univariate and multivariate logistic regression analyses were used to determine the ability of each variable to predict mortality. In the univariate logistic analysis, each variable was examined separately for its association with death. The odds ratio was defined as a ratio of odds in the death group divided by odds in the control group, and 95% confidence intervals, the Wald
2 statistic value, and P values were calculated. A P value of less than .05 was considered to indicate a significant difference. The odds ratio reflects an increase (>1) or decrease (<1) in the odds of death within 12 months for each unit of increase in the explanatory variable. Multivariate logistic regression analysis was applied to estimate the odds ratio of each potential predictor of mortality while other possible predictors were taken into account. Time-varying covariates were not used in the multivariate model because we wished to examine the predictive value of the baseline measurements, not of the follow-up measurements.
Survival analysis was performed to assess the value of each variable in predicting mortality while accounting for differences in duration of follow-up. To find the potential predictors we used the Cox proportional hazards model for univariate analysis. Proportional hazard assumptions were tested by using supremum tests and the exact method was used for ties.
| RESULTS |
|---|
|
|
|---|
|
A greater extent of fibrosis at baseline (Table 2) was the strongest independent predictor of death (P = .017). Kurtosis was the second-best predictor of death during follow-up (P = .072). As Table 2 included all available data points in the analysis, the logistic regression analysis was repeated by using only patients (n = 145) who had data for all eight baselines variables, and results were essentially unchanged; visual extent of lung fibrosis was again the most significant predictor (P = .003), followed by kurtosis.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Because of the documented accuracy of using CT images to establish a confident diagnosis of IPF (12), CT is often used in place of biopsy as an entry criterion in clinical trials (11). There is increasing evidence that the extent of lung disease on CT images is an important determinant of prognosis in patients with IPF (1).
While the role of quantitative CT indexes appears to be quite well established in patients with obstructive lung diseases, particularly emphysema (13–16), quantitative analysis of infiltrative lung diseases is more challenging, largely because it requires methods to identify and quantify the increase in lung attenuation that occurs in these conditions. Although automated analysis of multiple features (17) and fractal analysis (18) have been used, histogram-based methods have been used more extensively (19) and have shown reasonable correlation with pulmonary function test results in patients with asbestosis (20) and IPF (2,20). While the serial change in quantitative CT indexes observed in our study does not by itself confer clinical relevance, this change over 12 months suggests that these techniques may be used to identify progressive fibrosis. The greatest degree of change was found in skewness and kurtosis (P < .001 for each).
The 21.0% mortality rate in our study (median follow-up, 1.5 years) is comparable to published survival figures in patients with IPF (6,7). Our results from univariate logistic regression analysis show that baseline lung fibrosis, kurtosis, skewness, and TLC were statistically significant predictors of short-term mortality. The directionality of the values was as expected: Greater extent of fibrosis and lower TLC, skewness, and kurtosis values predicted mortality during follow-up. Survival analysis demonstrated that these trends did not change when the duration of follow-up was taken into account. These findings support the hypothesis that quantitative CT indexes, specifically kurtosis and skewness, are predictive of short-term mortality. However, multivariate logistic regression analysis showed that the single best predictor of mortality during follow-up was the visual extent of fibrosis on CT images, and there was little additional predictive power to be gained from the more quantitative measures. These results are very similar to those of a recent study by Lynch et al (1) that also found the visual extent of lung fibrosis to be the strongest independent predictor of mortality. That study evaluated physiologic features but did not include quantitative CT measurements. Given the subjectivity of visual assessment of disease extent, it is a little surprising that quantitative CT indexes, which correlate with extent of disease, were not stronger predictors of mortality. Our findings suggest that texture-sensitive image-analysis algorithms with the ability to segment fibrotic lung might be more useful than the relatively crude global quantitative CT indexes used in our study.
Early CT studies of patients defined as having IPF or fibrosing alveolitis suggested that the extent of GGO was a predictor of treatment responsiveness and suggested a better prognosis. However, these studies were published before the wide recognition of the entity of nonspecific interstitial pneumonia, which has a better prognosis than usual interstitial pneumonia and tends to be associated with GGO on CT images. The better prognosis associated with GGO in those studies was likely due to the inclusion of some patients with nonspecific interstitial pneumonia. Our results, like those of Lynch et al (1), indicate that in patients with IPF diagnosed according to current definitions (3), the extent of GGO does not have prognostic value.
Limitations of our study included the lack of spirometric gating in the acquisition of chest CT images. Differing levels of inspiration at CT could lead to substantial differences in quantitative CT indexes, including MLA (21–23). However, spirometric gating or control of CT acquisition remains technically cumbersome and is not widely available on clinical CT scanners. Our findings that quantitative CT can show progression of lung fibrosis even without spirometric control suggest that spirometric standardization may not be necessary for routine use. Another limitation was our specific study group, which included only patients with moderate disease who were eligible for treatment. Further, patient deaths likely had an effect on serial change measurements, as patients who died before follow-up likely had declines in physiologic measurements and progression of CT index values. The relatively short duration of follow-up was a further limitation. The fact that patients were randomized to treatment and placebo groups is a potential concern, but we believe that the results still have implications for an untreated IPF population, since treatment group was not a significant factor in outcome.
We conclude that, in patients with IPF, visual assessment of baseline disease extent on CT images is the strongest available independent predictor of short-term mortality and that serial evaluation of quantitative CT indexes can show disease progression.
| ADVANCES IN KNOWLEDGE |
|---|
|
|
|---|
| IMPLICATIONS FOR PATIENT CARE |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Abbreviations: FVC = forced vital capacity GGO = ground glass opacity IPF = idiopathic pulmonary fibrosis MLA = mean lung attenuation TLC = total lung capacity
See Materials and Methods for pertinent disclosures.
Author contributions: Guarantors of integrity of entire study, A.C.B., D.A.L.; study concepts/study design or data acquisition or data analysis/interpretation, all authors; manuscript drafting or manuscript revision for important intellectual content, all authors; approval of final version of submitted manuscript, all authors; literature research, A.C.B., D.M.; clinical studies, A.C.B., A.M.L., C.M.B., D.A.L.; statistical analysis, J.M., G.K.G.; and manuscript editing, A.C.B., J.M., C.M.B., D.M., G.K.G., D.A.L.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. M. Groves, T. Win, N. J. Screaton, M. Berovic, R. Endozo, H. Booth, I. Kayani, L. J. Menezes, J. C. Dickson, and P. J. Ell Idiopathic Pulmonary Fibrosis and Diffuse Parenchymal Lung Disease: Implications from Initial Experience with 18F-FDG PET/CT J. Nucl. Med., April 1, 2009; 50(4): 538 - 545. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Behr and V. J. Thannickal Update in Diffuse Parenchymal Lung Disease 2008 Am. J. Respir. Crit. Care Med., March 15, 2009; 179(6): 439 - 444. [Full Text] [PDF] |
||||
![]() |
R. Kim and K. C. Meyer Review: Therapies for interstitial lung disease: past, present and future Therapeutic Advances in Respiratory Disease, October 1, 2008; 2(5): 319 - 338. [Abstract] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |