Abstract
BACKGROUND: The purpose of this study was to correlate airway parameters of COPD determined by low-dose high-resolution computed tomography (HRCT) with pulmonary function testing (PFT) results.
METHODS: PFT data were collected for subjects with COPD and healthy controls. All subjects received inspiratory and expiratory phase low-dose HRCT. Bronchi in the apical segment of the right upper lobe (RB1), posterior segment of the right lower lobe (RB6), and lower lingual segment of the left upper lobe (LB5) were the target bronchi. Software automatically calculated airway wall area, inner area, and airway wall area percentage (percentage wall area for bronchial external area).
RESULTS: A total of 75 COPD and 20 control subjects were included. The subjects with COPD were classified according to COPD stage, with 20 grade I, II, and III subjects, respectively, and 15 grade IV subjects. In COPD grade II, residual volume/total lung capacity was negatively correlated with airway wall area in LB5 (r = −0.51). In COPD grade III, FVC was negatively correlated with airway wall area percentage in LB5 (r = −0.49) but positively correlated with airway wall area in RB6 (r = 0.52); percent-of-predicted FEV1 was negatively correlated with airway wall area percentage in RB1 (r = −0.49); residual volume was negatively correlated with airway wall area (r = −0.47), and total lung capacity was negatively correlated with airway wall area in RB1 (r = −0.52) (all, P < .001).
CONCLUSIONS: The results of this study suggest that airway parameters in different COPD grades have no uniform tendency of correlation with PFT, but some HRCT parameters are correlated to some PFT parameters.
Introduction
COPD can be simply defined as partial air-flow limitation. However, COPD is in fact a complex systemic disease and chiefly consists of chronic bronchitis and emphysema.1,2 Although small-airway pathology and emphysema can induce air-flow obstruction in patients with COPD, small-airway pathology is the most important factor. Micro-computed tomography (CT) studies have shown that the diameter of terminal bronchioles in subjects with severe COPD is 100 times smaller than that of normal healthy controls, and this is before the occurrence of emphysema.3,4 Another study5 has shown that airway wall thickness in subjects with COPD is related to external airway diameter. As COPD develops, large and small airways both undergo reconfiguration, thickening of airway walls, narrowing of airway lumen, and airway fibrosis after inflammation. Airway remodeling in subjects with COPD is the chief factor inducing air-flow obstruction.6
At present, although pulmonary function testing (PFT) is still the accepted standard for the confirmation and clinical grading of COPD, the results can be affected by a patient's age and other health factors, and PFT is associated with a relatively large margin of error.7 Furthermore, PFT may be normal in some patients in whom there are already abnormalities in small-airway function (diameter < 2 mm), especially in the early stage of the disease.8,9
High-resolution computed tomography (HRCT) allows the assessment of airways as small as 1 mm in diameter and is useful for examining changes of lung tissue and airways in pulmonary diseases, including asthma and COPD.10,11 CT is also useful for volume measurements, examination of tracheal morphology, quantification of emphysema, and analysis of perfusion in patients with COPD.12–15 A number of studies have sought to compare CT data, such as airway diameter and airway wall thickness, in subjects with and without lung disease (asthma, COPD, emphysema) and to correlate these parameters with disease severity.16–20 Much of the prior research, however, has consisted of quantitative studies of emphysema in subjects with COPD.7,21–23 The purpose of this study was to correlate airway parameters of subjects with COPD determined by low-dose HRCT with PFT to determine the feasibility of using CT to determine the severity of COPD.
QUICK LOOK
Current knowledge
High-resolution computed tomography (HRCT) is useful for examining changes of lung tissue and airways in pulmonary diseases, including asthma and COPD. Much of the prior research, however, has consisted of quantitative studies of emphysema in subjects with COPD.
What this paper contributes to our knowledge
The results showed that the wall area of the apical segment of the right upper lobe of subjects with COPD was greater than that of control subjects. The apical segment of the right upper lobe, posterior segment of the right lower lobe, and lower lingual segment of the left upper lobe wall area proportions of subjects with COPD were uniformly greater than those of the control group, and the higher values in the apical segment of the right upper lobe and posterior segment of the right lower lobe were significant. The percentage of wall area was best able to reflect small-airway pathology.
Methods
Subjects
This study recruited subjects with COPD seen at the First Affiliated Hospital of Guangzhou Medical University (Guangzhou Respiratory Disease Research Institute). The study was approved by the institutional review board of the hospital, and all subjects provided written informed consent.
Subjects were eligible for inclusion if they met the 2007 Global Initiative for Chronic Obstructive Lung Disease (GOLD) diagnostic standards.24 Briefly, the standards consist of evident chronic coughing, expectoration, FEV1/FVC < .70, and FEV1 < 80% of predicted value after inhaling a bronchodilator; no cough or expectoration, FEV1/FVC < .70, and FEV1 ≥ 80% of predicted value, and excluding other diseases; stable COPD (stable or mild symptoms of expectoration, coughing, and shortness of breath). Exclusion criteria included α1-antitrypsin deficiency, pneumonia, thoracic deformity, asthma, pulmonary fibrosis, lung cancer, tuberculosis, and history of thoracic surgery or other diseases that might affect pulmonary function. A normal control group consisted of subjects with normal PFT; no chronic coughing, expectoration, or shortness of breath; no evident respiratory pathologies before inclusion; no history of contact with dust or smoking history; no deformity of the chest profile caused by thoracic disease; and no other thoracic pathologies as determined by chest CT.
PFT and low-dose HRCT were performed within 3 d of each other. In all subjects, drugs that might affect the test results were suspended before testing because the subjects included in the study were considered stable. Before PFT and HRCT, subjects were administered 200 μg of salbutamol sulfate by inhalation.
Low-Dose HRCT
A Toshiba Aquilion 16-slice spiral CT machine was used. Scanning was performed after subjects inhaled deeply and held their breath, and the scanning time was 6–8 s. Scanning parameters were: tube voltage, 120 kV; tube current, 40 mA; pitch, 16; x-ray tube rotation velocity, 0.5 slice/rotation; field of view, 400 mm × 400 mm; acquisition array, 512 × 512; scanning acquisition thickness, 0.6 mm; image reconstructed slice thickness, 1 mm. A high-resolution reconstructive algorithm was used. Both inspiratory and expiratory studies were performed, but only inspiratory views were included in this study. This was because the airway was very constricted during the expiratory phase, especially due to collapse of the posterior wall, so the software could not accurately assess the condition of the airway during the expiratory phase. In addition, the majority of studies on this topic analyzed only the inspiratory phase, and expiratory views are primarily for visualizing the parenchyma.
Data Analysis
Vida Apollo 1.2 software developed by University of Iowa was employed for data analysis. After the collected DICOM data were entered into the computer system, the software automatically removed the mediastinum, major blood vessels, trachea, and tissue of the thoracic wall. This study opted to perform measurement of the apical segmental bronchus of the right upper lobe (RB1), posterior segmental bronchus of the right lower lobe (RB6), and lower lingual segmental bronchus of the left upper lobe (LB5) because these 3 bronchi have little variability and can reflect the overall condition of the airways of the lower lobes to a certain degree because of their location, although different approaches have been described.5,25 Lumen volume representation technology was used to perform post-processing of the images, and the software automatically generated CMPR images of the bronchial tree. In accordance with images of the bronchial trees, axial point array diagrams and virtual endoscopic images of the bronchi were used to track the route of the target bronchi. The target bronchi were located from different angles of rotation. Despite the fact that the planar scale was perpendicular to the long axis of the target bronchi, short-axis images were obtained of the target bronchi, measurements were made at a distance of 3 mm from the distal bronchial openings (Fig. 1), and the software automatically calculated the airway wall area, inner area, and airway wall area percentage. The mean values of these 3 measurements were calculated and used in the analysis.
Pulmonary Function Testing
A Quark PFT series (Cosmed, Rome, Italy) PFT instrument meeting American Thoracic Society/European Respiratory Society instrument quality control standards was used. After 10 min, the subject was instructed to breathe calmly for 3–4 cycles and then, after deep inhalation, suddenly, continuously, and steadily exhale completely with maximum effort to completion via an interface device. After complete exhalation with no interruption or coughing, the patient inhaled completely with maximum effort. Requirements of an adequate test consisted of a peak expiratory flow that quickly appeared on the maximum expiratory flow-volume curve, extrapolated capacity <5% FVC or 150 mL, smooth expiratory phase descending curve, and time capacity curve displaying an expiratory phase platform and volume change of <25 mL/s. At least 3 measurements were performed for FEV1, and FVC should be <150 mL. The maximum FVC and FEV1 values were reported, and the best of the other values were taken as the maximum FVC + FEV1 and measurement value with a smooth curve. PFT values obtained included FVC, percent-of-predicted FEV1, FEV1/FVC, residual volume (RV), and total lung capacity (TLC). TLC and RV were measured using body plethysmography, and the RV/TLC was calculated.
Statistical Analysis
Subject demographic characteristics, lung function data, and bronchial size were summarized as mean ± SD for continuous variables and n (%) for categorical variables. Differences among groups were compared using one-way analysis of variance with a post hoc pair-wise comparison. The Bonferroni test was used for continuous variables with normal distribution; the Kruskal-Wallis test was used with a post hoc pair-wise comparison, the Mann-Whitney U test was used for continuous variable without normal distribution; and the Pearson chi-square test was used for categorical variables. Two-way analysis of variance was also performed with a post hoc pair-wise comparison, the Bonferroni test, for continuous variables by considering some covariates (sex and body mass index [BMI]). A partial Pearson correlation analysis was performed to identify the correlation of lung function with bronchial size with adjustment for BMI and COPD grade, and a coefficient of correlation was reported. All statistical assessments were 2-tailed and considered significant at P < .05. An adjusted significance level of .01 (.05/5) was considered for the post hoc pair-wise comparisons. All statistical analyses were performed with IBM SPSS 22 for Windows (IBM Corp, Armonk, New York).
Results
A total of 75 subjects with COPD and 20 control subjects were included in the analysis, and demographic characteristics, lung function data, and bronchial sizes are summarized in Table 1. The subjects with COPD were classified into 4 groups according to COPD stage, with 20 grade I, II, and III subjects, respectively, and 15 grade IV subjects. The majority of the subjects with COPD were male, and the mean age of grade I, III, and IV subjects was greater than that of the control group (sex: P < .001 and < .001 for grade I and III subjects vs control group; age: all, P < .001 for grade I, III, and IV subjects vs control group). Lung function parameters were significantly different among the groups (all, P < .001, except for percent-of-predicted TLC, P = .03).
Bronchial sizes were significantly different among the groups, except for the inner area and airway wall area in RB6 (RB1: P < .001, P < .001, and P < .001 for inner area, airway wall area, and airway wall area percentage, respectively; RB6: P = .20, P = .48, and P < .001 for inner area, airway wall area, and airway wall area percentage, respectively; LB5: P < .001, P = .01, and P = .032 for inner area, airway wall area, and airway wall area percentage, respectively). COPD grade IV subjects had lower inner area in RB1 and LB5 than grade I subjects (P = .01 in RB1 and P < .001 in LB5). COPD grade IV subjects had a lower airway wall area in RB1 than grade II subjects (P = .01 in RB1). COPD grade IV subjects had a higher airway wall area percentage in RB1 and RB6 than control subjects (both, P < .001 in RB1 and RB6) and a higher airway wall area percentage in RB1 than grade I subjects (P < .001) (see Table 1).
The correlations of lung function with bronchial size adjusted for BMI and COPD grade are shown in Table 2. FVC was positively correlated with inner area in RB1, LB5, and RB6 (r = 0.25, 0.32, and 0.33, P = .02, P < .001, P < .001, respectively). FVC was also positively correlated with airway wall area in LB5 and RB6 (r = 0.26 and 0.25, P = .01 and P = .02, respectively). Percent-of-predicted FEV1 was positively correlated with airway wall area in LB5 (r = 0.22, P = .031), and FEV1/FVC was negatively correlated with airway wall area in RB6 (r = −0.24, P = .02). RV was negatively correlated with inner area in RB1 (r = −0.21, P = .047). The RV/TLC was negatively correlated with inner area and airway wall area in RB1 (r = −0.27 and −0.26, P = .01 and P = .01, respectively). There was no significant correlation between TLC and bronchial size.
Correlation analysis between lung function and bronchial size was also performed with adjustment for BMI and COPD grade separately (Table 3). In the control group, FVC was positively correlated with inner area in LB5 (r = 0.47, P = .043), FEV1/FVC was negatively correlated with airway wall area in RB6 (r = −0.59, P = .01), and TLC was positively correlated with airway wall area percentage in LB5 (r = 0.51, P = .03). In subjects with COPD grade I, there were no significant correlations. In subjects with COPD grade II, the RV/TLC was negatively correlated with airway wall area in LB5 (r = −0.51, P = .03). In subjects with COPD grade III, FVC was negatively correlated with airway wall area percentage in LB5 (r = −0.49, P = .033) but positively correlated with airway wall area in RB6 (r = 0.52, P = .02); percent-of-predicted FEV1 was negatively correlated with airway wall area percentage in RB1 (r = −0.49, P = .035); RV was negatively correlated with airway wall area in RB1 (r = −0.47, P = .040); and TLC was negatively correlated with airway wall area in RB1 (r = −0.52, P = .02).
Because the mean age of subjects with COPD grade I, III, and IV was greater than that of the control group, the correlation between bronchial size and lung function was analyzed with adjustment for age, BMI, and COPD grade (see supplementary Table 1 in the supplementary materials at http://www.rcjournal.com). The results were the same as those presented in Table 2, except for the correlation between airway wall area and FVC in RB6; however, the significance level was borderline (P = .050). Also, considering the age effect, correlation between bronchial size and lung function was analyzed with adjustment for age and BMI with respect to COPD grade (supplementary Table 2). In grade I COPD, an additional significant correlation was noted in RB1 (airway wall area percentage vs percent-of-predicted RV: r = −0.51, P = .03). In grade II COPD, 2 additional significant correlations were noted in LB5 (inner area vs percent-of-predicted FEV1: r = 0.53, P = .02; airway wall area percentage vs percent-of-predicted FEV1: r = −0.48, P = .045). In addition, in grade III COPD, 2 correlations that were significant in the prior analysis became nonsignificant (RB1: airway wall area vs percent-of-predicted RV: r = −0.46, P = .056; RB6: airway wall area vs FVC: r = 0.47, P = .050).
Lung function was compared between smokers and nonsmokers, and the results showed that smokers had poorer lung function than nonsmokers, except for the ratio of percentage-of-predicted TLC (supplementary Table 3). We then performed a subgroup analysis for the correlation of lung function with bronchial size that considered the effect of smoking on lung function (ie, correlation analysis was performed for smokers and nonsmokers) (supplementary Table 4). The analysis showed that there was no significant correlation between lung function and bronchial size in RB1. In LB5 in nonsmokers, there was a significant correlation between inner area and FVC (r = 0.42, P = .02); inner area and percent-of-predicted RV (r = −0.39, P = .032); inner area and percent-of-predicted RV/TLC (r = −0.51, P < .01); airway wall area and percent-of-predicted FEV1 (r = 0.39, P = .034); and airway wall area and percent-of-predicted RV/TLC (r = −0.48, P = .01). There was no significant correlation between lung function and bronchial size in smokers. In RB6, a significant correlation was found between inner area and FEV (r = 0.28, P = .031) in smokers and between airway wall area and FEV1/FVC (r = −0.488, P = .01) in non-smokers. Subgroup analysis with respect to COPD grade was not performed because the data were not suitable for subgroup analysis.
Discussion
This study measured parameters of the apical segmental bronchus of the right upper lobe (RB1), posterior segmental bronchus of the right lower lobe (RB6), and lower lingual segmental bronchus of the left upper lobe (LB5), which can reveal pathology of the large airways. The primary result was that the airway wall area percentage values of subjects with COPD were uniformly greater than those of healthy control subjects, and the airway wall area percentage values were greater than those in subjects with COPD grade I, II, and III disease. However, we must emphasize that statistical significance was not observed.
Pairwise comparisons of control group data and the data of subjects with different grades of COPD showed that although there was no statistically significant difference between some groups and although early subjects with COPD may have slightly enlarged airway lumens, the airway lumen was significantly smaller in grade III and grade IV subjects. In addition, whereas the airway wall area may be relatively large in grade I and II COPD, it begins to shrink by grade III and has shrunk significantly by grade IV. Airway wall area percentage gradually increases with the severity of COPD and is significantly increased in grade III and IV disease. These findings indicate that although airway lumen and airway wall area have both decreased by late-stage COPD, the decrease in airway lumen area tends to be more evident than the decrease in the airway wall area. Nevertheless, previous studies have suggested that the airway walls of subjects with COPD have a greater thickness than those of normal persons, and as the condition worsens, airway walls will gradually thicken, the lumen area will gradually decrease, and the wall area will gradually increase.26–28 This will result in a significant increase of airway wall area percentage, giving meaning to this parameter. On the other hand, Kosciuch et al18 found no significant correlations between airway wall thickness and spirometric parameters in subjects with asthma or COPD; however, in COPD, airway wall area percentage was positively correlated with airway resistance (r = 0.72, P < .001). Copley et al29 studied asymptomatic elderly subjects and reported that wall thickening was statistically more frequent in subjects who were >75 y old than in those <55 y of age.
Although the sample size in our study was relatively small, a control group was included, and measurements were made along a plane as perpendicular as possible to the long axis of the airway. In addition, in contrast to manual or semi-automatic measurements, the relatively mature Vida fully automated measurement software, which has been validated in other studies,30–32 was used. We thus believe that discrepancies between the results of other studies and our results are not just a variation between the methods but reflect the real airway status. It appears that airway remodeling caused by chronic inflammation and fibrosis appearing at late stages will cause the airway lumen to shrink, and the airway wall will become correspondingly thinner. However, this conclusion must be verified by pathological examination.
Past studies of quantitative measurements of bronchial size in subjects with asthma and COPD have suggested a correlation between bronchial size and pulmonary function,33,34 although this depends on the parameters examined and has not been shown in some studies.18 The current study also found correlations between the airway size of subjects with COPD and pulmonary function, including RB1 in subjects with grade IV COPD and pulmonary function indices. However, at an advanced stage (GOLD grade IV), changes in RB6 and LB5 not were observed. Although we cannot know for certain, we believe that these results can be explained as follows. Bronchial plasticity was not evident in grade I subjects, whereas in grade II subjects, bronchial plasticity began to occur, and there simultaneously existed inflammatory edema. Although grade III and IV subjects have more severe disease, their bronchial plasticity was evident, a high degree of fibrosis was present, and the fibrosis caused even greater contraction of the bronchial wall. As for whether different grades have effects in different locations, clinical experience indicates that inflammation significantly changes the parallel routing of the bronchi to a greater degree, and RB1 is vertical or oblique. Because of this, this bronchus (RB1) may be influenced even more significantly in patients with a higher COPD grade. Plasticity should be even more evident in the case of RB6 and LB5, and inflammatory edema of the bronchial wall might not be very severe.
Among airway inner area, airway wall area, and airway wall area percentage values, inner area and airway wall area percentage had relatively high correlations with pulmonary function. Inner area was positively correlated with FVC and negatively correlated with RV/TLC. This indicates that the smaller the airway lumen, the more pulmonary ventilation and small-airway function will deteriorate, and RV will worsen even more. Airway wall area percentage was negatively correlated with FVC, FEV1, and FEV1/FVC and positively correlated with RV/TLC, which indicates that the larger the airway wall area and the smaller the airway lumen, the greater the degree of air-flow obstruction. This further indicates that bronchial changes not only occur in peripheral small airways, but also affect central, relatively large airways. Previous studies have reported a negative correlation between quantitative airway wall thickness determined by CT and pulmonary function, a negative correlation between airway wall area percentage and percent-of-predicted FEV1 and FEV1/FVC, and a positive correlation between airway wall area percentage and RV/TLC.27,35 Patel et al36 have suggested that the best correlation exists between airway wall thickness and RV/TLC.
It should be noted that although the inner areas for subjects with COPD and controls were not significantly different, the data in Table 2 appear to show a significant correlation between inner area and FVC for all 3 segments examined. This is because the subjects included for the analysis in the tables are different. The correlation noted between inner area and FVC was probably caused by grade (Table 1), and when all of the subjects were pooled (Table 3), the correlations differed.
There are certain limitations, as well as strengths of the current study. The numbers of subjects with the different grades of COPD were relatively small. The study, however, evaluated a sufficient number of bronchi to select optimal sites for diagnosis, and this study focused on airway parameters, whereas prior similar studies focused on emphysema. Given the pathogenesis of COPD, detailed airway measurement probably provides more meaningful information. It would be ideal if it were possible to measure all 18 bronchopulmonary segments; however, the work involved in this would be beyond the scope of most researchers. RB1 represents the upper lobe and has been used in most prior studies, and because RB6 and LB5 represent the largest lobes of the right and left lungs, analysis of them makes it relatively easy to assess changes. We believe that assessment of these 3 segments is representative of the condition of the lobes they represent. It may have been beneficial to analyze emphysema and airway size together; however, other studies have examined airway dimensions and emphysema, and the purpose of this study was to focus on airway size and simplify the analysis. Both inspiratory and expiratory studies were performed, but only inspiratory views were included in this study. This was because in expiratory views, the airway was very constricted during the expiratory phase, especially due to collapse of the posterior wall, so the software could not accurately assess the condition of the airway. Coupled inspiratory-expiratory studies, however, have proven to be informative in some reports, because assessment of changes in inspiration-expiration is important for correlating bronchial size with PFT.35–37 TLC and RV were measured using body plethysmography, and results obtained with this method may be discrepant with those of other methods, such as helium dilution. COPD affects the airway wall and lung parenchyma to varying degrees, and air trapping and air wall thickening are associated. The results showed no correlation between the degree of bronchial obstruction and air trapping. Air trapping is more related to lung parenchyma changes, and the aim of the present study is to evaluate the airway wall changes among different degrees of COPD. For this reason, we did not analyze parenchyma changes.
Conclusions
The results of this study suggest that airway parameters in different COPD grades have no uniform tendency of correlation with PFT, but some HRCT parameters are correlated to some PFT parameters. This may allow the use of HRCT as an easier and more objective assessment of patients with COPD, as compared with PFT, in routine clinical practice. Our results indicated that in COPD grade II, RV/TLC was negatively correlated with airway wall area in LB5. In COPD grade III, FVC was negatively correlated with airway wall area percentage in LB5 but positively correlated with airway wall area in RB6; percent-of-predicted FEV1 was negatively correlated with airway wall area percentage in RB1; RV and TLC were negatively correlated with airway wall area in RB1. A complete HRCT assessment system for patients with different grades of COPD could be developed after more studies and more detailed analysis of this topic.
Footnotes
- Correspondence: Shi-jun Qiu MD, Medical Imaging Center, Nanfang Hospital, Southern Medical University, No. 1838, Guangzhou Dadaobei Road, Baiyun District, Guangzhou City, Guangdong Province 510515, China. E-mail: shijunqiu1963{at}163.com.
This study was supported by Science and Technology Planning Project of Guangdong Province, China (Grant 2014A020212608) and Public Science and Technology Research Funds Projects of NHFPC of the People's Republic of China (Grant 201402013). The authors have disclosed no conflicts of interest.
Supplementary material related to this paper is available at http://www.rcjournal.com.
- Copyright © 2017 by Daedalus Enterprises