Computed tomography-based visual assessment of chronic obstructive pulmonary disease: comparison with pulmonary function test and quantitative computed tomography.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) has variable subtypes involving mixture of large airway inflammation, small airway disease, and emphysema. This study evaluated the relationship between visually assessed computed tomography (CT) subtypes and clinical/imaging characteristics.
METHODS: In total, 452 participants were enrolled in this study between 2012 and 2017. Seven subtypes were defined by visual evaluation of CT images using Fleischner Society classification: normal, paraseptal emphysema (PSE), bronchial disease, and centrilobular emphysema (trace, mild, moderate and confluent/advanced destructive). The differences in several variables, including clinical, laboratory, spirometric, and quantitative CT features among CT-based visual subtypes, were compared using the chi-square tests and one-way analysis of variance.
RESULTS: Subjects who had PSE had better forced expiratory volume in 1 second (FEV1) (P=0.03) percentage and higher lung density (P<0.05) than those with moderate to confluent/advanced destructive centrilobular emphysema. As the visual grade of centrilobular emphysema worsened, pulmonary function declined and modified Medical Research Council, COPD assessment test (CAT) score, and quantitative assessment (emphysema index and air trapping) increased. The bronchial subtype was associated with higher body mass index (BMI), better lung function and higher lung density. Participants with trace emphysema showed a rapid increase in functional small airway disease.
CONCLUSIONS: Classifying subtypes using visual CT imaging features can reflect heterogeneity and pathological processes of COPD. 2021 Journal of Thoracic Disease. All rights reserved.

Introduction

Chronic obstructive pulmonary disease (COPD) is characterized by nonreversible airflow limitation caused by the inhalation of irritants such as tobacco smoke (1). Approximately 15‒20% of smokers are likely to develop COPD, and there is substantial heterogeneity in the clinical manifestation, physiological characteristics, imaging characteristics, disease progression, and prognosis of the disease (1,2). COPD is routinely diagnosed based on pulmonary function test (PFT) results, respiratory symptoms, and history of exposure (3). Nevertheless, the technique has several limitations. PFT results are difficult to evaluate in patients with respiratory problem, and tests tends to be relatively insensitive at the early stage of the disease, and cannot distinguish the subtypes of COPD, including emphysema, small airway disease, and bronchial inflammation (2,4).

Computed tomography (CT) is a useful imaging modality for evaluating structural abnormalities and for classifying the COPD subtype. Quantitative CT evaluation is useful for assessing the severity of emphysema, airway wall thickening, and small airway disease (5). Emphysema quantification by means of CT correlated well with histopathology (6,7). However, CT densitometry has limited ability to provide information about the emphysema phenotype. Recently, visual categorization based on CT was proposed by the Fleischner Society, and emphysema has been divided into centrilobular emphysema, paraseptal emphysema (PSE), and panlobular emphysema (8). A study reported that the visual presence and severity of emphysema was an independent predictor of mortality (9). However, the study analyzed only centrilobular emphysema, and did not analyze other COPD phenotypes, such as PSE and bronchial disease. Recently, Park et al. (2) reported that the phenotypes of COPD could classified by combining visual and quantitative CT imaging features, presenting a useful tool for characterizing the heterogeneity of COPD. Nevertheless, quantitative CT evaluation of COPD is not yet widely used in clinical practice. Therefore, we here classified COPD phenotypes by means of visual characterization only, in order to evaluate the relationship of visually assessed CT subtypes with clinical presentation and quantitative CT characteristics.

We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/jtd-20-3041).

Methods

Study population

Five-hundred-and-four participants from the Chronic Obstructive Pulmonary Disease in Dusty Areas (CODA) cohort, a Korean study that observed the clinical outcomes of COPD in Korean patients residing near cement plants, were initially selected. The details of this cohort had been described previously (10). All enrolled participants underwent medical interview, physical examination, spirometry, and chest CT scan. We excluded 52 participants with quantification error (n=8), lobectomy state (n=4), and severe lung parenchymal distortion due to pneumoconiosis, bronchiectasis, and pulmonary tuberculosis (n=30), and accompanying lung disease (n=10) that could affect the quantitative CT measurement. Finally, 452 participants were enrolled (). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). This study was approved by the institutional review board of our hospital (approval #2012-06-007), and all participants provided written informed consent.

Figure 1

Flowchart of the study participants.

Clinical parameters and pulmonary function

The interview questionnaire included questions on demographic data, medical history, environmental exposure, and respiratory symptoms. Dyspnea was assessed using the modified Medical Research Council, and health-related quality of life was assessed using the patient-reported COPD assessment test (CAT). The peripheral levels of white blood cells (WBCs), neutrophil, eosinophils, interleukin (IL)-6, IL-8, and C-reactive protein were used as markers of inflammation.

PFTs were performed using the Easy One Kit (NDD Medizintechnik AG, Zurich, Switzerland), before and after inhalation of 400 µg salbutamol. The spirometric airflow limitation is defined as a postbronchodilator ratio of forced expiratory volume in 1 second (FEV1) to forced vital capacity (FVC) of <0.7. The severity of airflow limitation was classified as described by the 2017 Global Initiative for Chronic Obstructive Lung Disease (GOLD) grading system as grade 1 (≥80%), grade 2 (50‒79%), grade 3 (30‒49%), or grade 4 (<30%) (11). GOLD grade 0 was defined as normal spirometry, and the group with normal FEV1/FVC and FEV1 <80% predicted was defined as having preserved ratio impaired spirometry (PRISm) (3). Bronchodilator responsiveness (BDR) was defined as an increase of at least 12% or 200 mL in FEV1 or FVC after salbutamol administration (12).

CT acquisition

Participants underwent thin-section, chest CT at full inspiration and expiration in the supine position. All chest CT examinations were performed using a Somatom Definition dual-source CT scanner (Siemens Healthcare, Forchheim, Germany) with the following parameters: 140 kVp, 100 mA, 09-1 beam pitch and 0.6 mm slice thickness. CT images were reconstructed with B30f.

Visual analysis and subtypes based on CT images

CT scans were visually assessed on the basis of the Fleischner Society classification system (8). Centrilobular emphysema was graded as trace (ca. 0.5%), mild (0.5‒5%), moderate (>5%), confluent (coalescent centrilobular or lobular lucencies) and advanced destructive emphysema (panlobular lucencies with hyperexpansion and distortion). Substantial PSE was defined as large cyst-like lucencies located near the pleura and interlobar fissure. Visual assessment of CT was independently interpreted by 2 radiologists (with 3 and 11 years of experience in thoracic CT scan interpretation) who were unaware of clinical data. Different interpretations were resolved by consensus.

Participants were categorized into 7 subtypes based on visual assessment of CT images: (I) normal (subjects with no visual emphysema and bronchial wall thickening), (II) PSE (subjects with substantial PSE regardless of other imaging characteristics), (III) bronchial airway disease (subjects with bronchial wall thickening and absent or trace emphysema or mild PSE), (IV) trace centrilobular emphysema, (V) mild centrilobular emphysema, (VI) moderate centrilobular emphysema, and (VII) confluent and advanced destructive emphysema () (2).

Figure 2

Classification of 7 computed tomography (CT) subtypes. (A) Normal. (B) Paraseptal emphysema; CT scan shows numerous low attenuation areas in the subpleural regions. (C) Bronchial disease; CT scan shows diffuse bronchial wall thickening without emphysema. (D) Trace centrilobular emphysema; CT scan shows scattered centrilobular low attenuated lesions involving about 0.5‒5% of the upper lung zone. (E) Mild centrilobular emphysema; CT scans shows well-defined low attenuation lesions involving about 0.5‒5% of the upper lung zone. (F) Moderate centrilobular emphysema; CT scans shows well-defined low attenuated lesions involving more than 5% of the upper lung zone. (G) Confluent and advanced destructive emphysema; CT scans shows hyperexpanded secondary pulmonary lobules and multiple hyperlucencies in both lungs.

Quantitative analysis based on CT image

Lung segmentation and image co-registration of inspiratory and expiratory CT were performed using an Aview® system (Coreline Soft Inc., Seoul, South Korea). Parametric response mapping (PRM) analysis was classified into PRMfSAD (≥−950 HU at inspiration and <-856 HU at expiration), or PRMemph (<−950 HU at inspiration and <−856 HU at expiration, ) (13). PRM data were represented as a percentage of the total lung volume. Lung attenuation at the 15th percentile (Perc 15) was defined as CT attenuation at the 15th percentile of the lung CT histogram (14). Pi10 was defined as the airway wall area standardized to an internal airway perimeter of 10 mm.

Figure 3

Parametric response mapping in 3 different patients. Normal lung tissue is denoted in green (PRMnormal), functional small airway disease in yellow (PRMfSAD), and emphysematous lung in red (PRMemphy). (A) Image of a 56-year-old male with COPD GOLD grade 2. CT shows significant bronchial wall thickening without emphysema. CT subtype was assigned as bronchial disease subtype. PRMemphy was 0.7% and PRMfSAD was 19.2%. Pi-10 was 5.5 mm. (B) Image of a 78-year-old male with COPD GOLD grade 4. CT shows paraseptal emphysema in both upper lungs. CT subtypes was classified paraseptal emphysema. PRMemphy was 18.2% and PRMfSAD was 50.4%. (C) Image of a 78-year-old male with COPD GOLD grade 2. CT show diffuse lobular lucencies with hyperexpansion in both lungs. CT subtype was classified as advanced destructive emphysema. PRMemphy was 38.5% and PRMfSAD was 30.0%. CT, computed tomography; COPD, chronic obstructive pulmonary disease; GOLD, Global Initiative for Chronic Obstructive Lung Disease.

Statistical analysis

Parametric data are expressed as mean ± standard deviation (SD), whereas nonparametric data are expressed as numbers and percentages. The chi-square test and one-way analysis of variance were used to determine differences in clinical characteristics and quantitative CT variables among the subtypes. P values <0.05 were considered statistically significant. Data analysis was performed using SAS version 9.3 (SAS Institute Inc., Cary, NC, USA).

Results

For the 452 participants, the mean age at enrollment was 72.4±7.4 years, and 327 of these participants (72.4%) were male. Two-hundred-and-eighty-five participants (63.5%) were current or former smokers at enrollment, while 3 participants had no recorded smoking history. With the subjects group divided into GOLD grades, 142 participants (31.4%) were GOLD grade 0, 20 participants (4.4%) were PRISm, 148 participants (32.7%) were GOLD grade 1, 119 participants (26.3%) were GOLD grade 2, 21 participants (4.7%) were GOLD grade 3 and only 2 participants (0.5%) were GOLD grade 4. Detailed characteristics of the participants are shown in .

Table 1

Characteristics of 452 participants

Characteristics	Value
Age	72.4±7.4
Male*	327 (72.4)
BMI	23.6±3.2
Smoking status*
Current	98 (21.8)
Former	187 (41.7)
None	164 (36.5)
Pack-years	17.5±23.6
mMRC	1.36±1.1
CAT score	16.1±5.6
WBC (×103/μL)	6.8±1.9
Neutrophil (%)	58.5±9.5
Eosinophil (%)	2.7±2.8
CRP (mg/dL)	0.3±0.6
IL-6 (pg/mL)	2.5±3.5
IL-8 (pg/mL)	17.7±21.3
GOLD grades*
0	142 (31.4)
PRISm	20 (4.4)
1	148 (32.7)
2	119 (26.3)
3	21 (4.7)
4	2 (0.5)
Pulmonary function test
FVC (L)	3.0±0.8
FVC (%)	97.9±18.7
FEV1 (L)	1.96±0.6
FEV1 (%)	87.5±22.3
FEV1/FVC	65.2±11.5
BDR* (%)	70 (15.5)

Data are shown as mean ± standard deviation unless indicated otherwise. *Data are number ad data in parentheses are percentages. BDR, bronchodilator responsiveness; BMI, body mass index; CAT, chronic obstructive pulmonary disease assessment test; CRP, C-reactive protein; FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; GOLD, Global Initiative for Chronic Obstructive Lung Disease; mMRC, modified Medical Research Council.

Overall comparison of subtypes

Subjects were divided into 7 subtypes. There were 212 (46.9%) subjects with no abnormality in visual CT analysis. Thirty-one (6.9%) subjects had PSE type and 10 (2.2%) had bronchial type. One-hundred-and-ninety-nine (44.0%) subjects had centrilobular emphysema on visual CT analysis and trace, mild and moderate centrilobular emphysema (11.1%, 20.1% and 9.7% each respectively) were predominant type among the 7 subtypes ().

Table 2

Clinical and imaging characteristics in 7 computed tomography imaging subtypes

Characteristics	Normal	PSE	Bronchial	Trace	Mild	Moderate	Confluent and advanced destructive	P value
Number of subjects	212 (46.9)	31 (6.9)	10 (2.2)	50 (11.1)	91 (20.1)	44 (9.7)	14 (3.1)
Age	71.3±7.6	73.4±6.3	72.0±10.2	73.6±7.0	73.3±6.8	73.3±8.0	73.7±7.0	0.170
Male*	119 (56.1)	31 (100.0)	5 (50.0)	36 (72.0)	80 (87.9)	42 (95.5)	14 (100.0)	<0.001
BMI	24.6±3.0	22.5±3.6	24.2±3.2	24.3±3.2	22.5±2.8	21.8±3.1	20.4±2.3	<0.001
Smoking status*								<0.001
Current	26 (12.3)	10 (32.3)	3 (30.0)	8 (16.0)	27 (29.7)	20 (48.8)	4 (28.6)
Former	63 (29.7)	20 (64.5)	2 (20.0)	22 (44.0)	51 (56.0)	19 (46.3)	10 (71.4)
None	123 (58.0)	1 (3.2)	5 (50.0)	20 (40.0)	13 (14.3)	2 (4.9)	0 (0.0)
mMRC	1.1±1.1	1.5±1.1	1.3±1.3	1.6±1.2	1.6±1.2	1.5±1.2	1.8±1.3	0.006
CAT score	14.5±9.3	20.5±9.6	16.4±9.6	17.1±10.0	17.4±9.1	16.6±9.9	18.1±11.6	0.020
WBC (×103/μL)	6.6±1.8	7.5±1.7	7.7±1.8	6.5±2.0	6.8±2.0	6.7±2.1	7.2±1.5	0.116
Neutrophil (%)	57.3±10.1	62.0±8.1	60.0±10.4	57.2±9.8	60.0±8.7	59.0±8.4	59.3±7.9	0.087
Eosinophil (%)	2.6±2.3	2.7±4.5	2.4±1.5	3.3±4.8	2.4±2.1	2.8±2.0	2.0±1.5	0.717
CRP (mg/dL)	0.3±0.7	0.4±0.9	0.2±0.3	0.2±0.2	0.2±0.5	0.3±0.4	0.4±0.7	0.786
IL-6 (pg/mL)	2.0±3.0	4.9±7.0	3.2±4.5	2.3±2.5	2.5±2.5	2.9±3.5	2.8±3.3	0.004
IL-8 (pg/mL)	16.8±20.5	21.8±28.1	30.6±34.7	16.8±26.6	17.2±1.3	18.9±25.8	14.0±9.4
GOLD grade*								<0.001
0	12 (5.6)	1 (3.2)	2 (20.0)	0 (0.0)	5 (5.5)	0 (0.0)	0 (0.0)
PRISm	106 (50.0)	3 (9.7)	2 (20.0)	14 (28.0)	11 (12.1)	6 (13.6)	0 (0.0)
1	56 (26.4)	13 (41.9)	3 (30.0)	19 (38.0)	42 (46.1)	13 (29.6)	2 (14.3)
2	33 (15.6)	10 (32.3)	3 (30.0)	15 (30.0)	30 (33.0)	21 (47.7)	7 (50.0)
3	5 (2.4)	3 (9.7)	0 (0.0)	2 (4.0)	2 (2.2)	4 (9.1)	5 (35.7)
4	0 (0.0)	1 (3.2)	0 (0.0)	0 (0.0)	1 (1.1)	0 (0.0)	0 (0.0)
Pulmonary function
FVC (L)	2.9±0.8	3.3±0.9	2.6±0.9	3.0±0.9	3.1± 0.8	3.1±0.8	3.1±0.7	0.1370
FVC (%)	99.4±17.9	97.1±22.6	93.6±22.8	100.0±18.8	96.6±19.9	95.6±15.5	88.3±16.9	0.2860
FEV1 (L)	2.0±0.6	1.9±0.7	1.6±0.5	2.0±0.6	1.9±0.5	1.9±0.7	1.6±0.6	0.0270
FEV1 (%)	93.8±21.1	77.5±20.9	79.6±20.3	91.6±23.7	83.3±19.8	78.5±21.4	61.8±19.3	<0.001
FEV1/FVC	69.7±10.7	57.7±10.9	63.3±10.3	66.0±9.0	62.3±8.5	59.1±12.6	50.5±12.4	<0.001
BDR (%)	25 (11.8)	7 (22.6)	3 (30.0)	10 (20.0)	17 (18.7)	7 (15.9)	1 (7.1)	0.2980
Perc 15 (HU)	−902.0±31.7	−941.6±24.0	−886.6±46.4	−913.0±24.1	−922.3±22.0	−933.2±17.0	−962.5±15.1	<0.001
Pi10 (mm)	4.7±0.4	4.6±0.5	4.9±0.5	4.7±0.4	4.6±0.35	4.7±0.4	4.5±0.1	0.241
PRMemph	2.1±2.9	12.5±8.7	1.3±1.5	3.1±2.9	5.3±4.2	8.2±4.7	22.5±9.3	<0.001
PRMfSAD	19.2±15.3	29.8±14.8	17.2±17.0	29.4±19.2	28.3±15.4	33.7±15.8	35.0±8.9	<0.001

Data are shown as mean ± standard deviation unless indicated otherwise. *Data are number ad data in parentheses are percentages. BDR, bronchodilator responsiveness; BMI, body mass index; CAT, chronic obstructive pulmonary disease assessment test; CRP, C-reactive protein; FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; GOLD, Global Initiative for Chronic Obstructive Lung Disease; mMRC, modified Medical Research Council; Perc 15, lung attenuation at 15th percentile; PRM, parametric response mapping; PSE, paraseptal emphysema.

Body mass index (BMI) lowered statistically significantly as the degree of emphysema, based on CT findings, worsened. The bronchial group tended to have a higher BMI (24.2±3.2) than the emphysema groups, except for the normal group. The CAT score was highest in subjects with PSE (20.5±9.6). Although there was no statistically significant difference, neutrophil counts were the highest in subjects with PSE. In subjects with PSE, IL-6 was higher than in other subtypes (P=0.004). Pulmonary function lowered as the degree of emphysema worsened.

The Pi-10 was greatest in the bronchial disease subtype. The PRMemph and PRMfSAD raised and Perc 15 lowered as the visual grade of emphysema worsened (all P<0.001). The PRMemph gradually increased in the order of the normal, trace emphysema, and mild emphysema subtypes. However, PRMfSAD was markedly increased in the trace emphysema group compared to the normal group and decreased slightly in the mild emphysema group as compared to the trace emphysema group. The differences among the subtypes are shown in .

PSE compared with moderate and severe emphysema

Participants with confluent and advanced destructive emphysema had the lowest lung function and highest emphysema (P<0.05). Compared with participants with moderate centrilobular emphysema, participants with PSE had lower lung function and Perc 15, and higher PRMemph. Although there were no statistically significant differences, functional small airway disease was least common in participants with PSE ().

Table 3

Clinical and imaging characteristics in paraseptal emphysema, moderate centrilobular emphysema and severe emphysema subtypes

Characteristics	PSE	Moderate	Confluent and advanced destructive	P value
Age	73.4±6.3	73.3±8.0	73.7±7.0	0.984
Male*	31 (100.0)	42 (95.5)	14 (100.0)	0.351
BMI	22.5±3.6	21.8±3.1	20.4±2.3	0.126
Smoking status*				0.395
Current	10 (32.3)	20 (48.8)	4 (28.6)
Former	20 (64.5)	19 (46.3)	10 (71.4)
None	1 (3.2)	2 (4.9)	0 (0.0)
mMRC	1.5±1.1	1.5±1.2	1.8±1.3	0.700
CAT score	20.5±9.6	16.6±9.9	18.1±11.6	0.278
WBC (×103/μL)	7.5±1.7	6.7±2.1	7.2±1.5	0.183
Neutrophil (%)	62.0±8.1	59.0±8.4	59.3±7.9	0.280
Eosinophil (%)	2.7±4.5	2.8±2.0	2.0±1.5	0.738
CRP (mg/dL)	0.4±0.9	0.3±0.4	0.4±0.7	0.837
IL-6 (pg/mL)	4.9±7.0	2.9±3.5	2.8±3.3	0.219
IL-8 (pg/mL)	21.8±28.1	18.9±25.8	14.0±9.4	0.679
GOLD grade*				0.116
0	3 (9.7)	6 (13.6)	0 (0.0)
PRISm	1 (3.2)	0 (0.0)	0 (0.0)
1	13 (41.9)	13 (29.6)	2 (14.3)
2	10 (32.3)	21 (47.7)	7 (50.0)
3	3 (9.7)	4 (9.1)	5 (35.7)
4	1 (3.2)	0 (0.0)	0 (0.0)
Pulmonary function
FVC (L)	3.3±0.9	3.1±0.8	3.1±0.7	0.561
FVC (%)	97.1±22.6	95.6±15.5	88.3±16.9	0.323
FEV1 (L)	1.9±0.7	1.9±0.7	1.6±0.6	0.237
FEV1 (%)	77.5±20.9	78.5±21.4	61.8±19.3	0.031
FEV1 /FVC	57.7±10.9	59.1±12.6	50.5±12.4	0.072
BDR (%)	7 (22.6)	7 (15.9)	1 (7.1)	0.428
Perc 15 (HU)	−941.6±24.0	−933.2±17.0	−962.5±15.1	<0.001
Pi10 (mm)	4.6±0.5	4.7±0.4	4.5±0.1	0.481
PRMemph	12.5±8.7	8.2±4.7	22.5±9.3	<0.001
PRMfSAD	29.8±14.8	33.7±15.8	35.0±8.9	0.419

Data are shown as mean ± standard deviation unless indicated otherwise. *Data are number ad data in parentheses are percentages. BDR, bronchodilator responsiveness; BMI, body mass index; CAT, chronic obstructive pulmonary disease assessment test; CRP, C-reactive protein; FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; GOLD, Global Initiative for Chronic Obstructive Lung Disease; mMRC, modified Medical Research Council; Perc 15, lung attenuation at 15th percentile; PRM, parametric response mapping; PSE, paraseptal emphysema.

Discussion

In this study, we demonstrate that COPD phenotypes could be evaluated by using only visual assessment of CT images. We found that pulmonary function and CT quantitative characteristics of PSE were intermediate to those of moderate emphysema and confluent advanced destructive emphysema. Participants with PSE had the highest CAT score, neutrophil count, and IL-6 levels among the subtypes. In addition, participants with trace emphysema showed a higher degree in functional small airway disease among the emphysema subtypes. Pulmonary function and quantitative features worsened with the visual grade of emphysema.

The diagnosis, severity evaluation and longitudinal monitoring of COPD usually depend on PFTs, including FEV1/FVC and FEV1. However, PFT can underestimate lung function in younger patients and overestimate lung function in older (15). A previous study found that FEV1 did not significantly reflect the progression of emphysema and air trapping detected by CT (less than 10% of emphysema and up to 50% for air trapping) (16). In addition, PFT provides a global evaluation of lung damage, rather than regional information (17). CT could be used to evaluate the extent, subtype, distribution and progression of COPD (15). Quantitative CT analysis provides global and regional information regarding emphysema, air trapping, and airways (8), and in patients with COPD, quantitative CT analysis generally correlated with pathology-based quantification of emphysema and the results of PFT (17). A study showed that emphysema associated with the FEV1/FVC, Pi10 did with the FEV1, and the ratio of mean lung density at expiration and inspiration (E/I-ratioMLD) did with the residual volume (18). PRM for COPD quantification is a voxel-based image analysis technique between inspiration and expiration scans after image registration (19). Pompe et al. showed that PRMemph and PRMfSAD strongly associated with presence and severity of COPD (20). In addition, the baseline model (age, BMI, smoking status and pack-years) with PRM had higher diagnostic value compared with other CT-derived biomarkers (Pi10, Perc 15, E/I-ratioMLD) (20). Therefore, in our study, PRMemph and PRMfSAD were used as the COPD quantification biomarkers. However, the heterogeneous nature of COPD cannot be explained by quantitative CT features only (4). In addition to quantitative analysis, visual evaluation of CT images provides information on patterns of altered lung structure and distinct phenotypes in COPD that cannot currently be identified by quantitative CT (8). A previous study showed that visual assessment of COPD could provide reproducible, physiologically significant information (21).

Emphysema is classified into centrilobular emphysema, PSE, and panlobular emphysema. Centrilobular emphysema shows destruction of the center of the lung lobule, whereas PSE shows destruction of peripheral lung lobule located near the pleural surface and in the interlobular fissures (8). The clinical manifestation and pathologic features differ between PSE and centrilobular emphysema. Compared with subjects with PSE, subjects with centrilobular emphysema tended to have more dyspnea, reduced exercise capacity, and impaired lung function (2,22). The number of terminal bronchioles are relatively preserved in PSE-dominant regions compared to centrilobular emphysema dominant regions in patients with COPD, which implies that patients with the PSE had better physiological and functional characteristics (23). In our study, participants with PSE had the highest CAT score, particularly regarding items on sputum and dyspnea when climbing stairs. Our study classified PSE subtypes in subjects with substantial PSE regardless of other characteristics; most subjects with moderate to advanced destructive emphysema were included in the PSE subtype group. Therefore, this group may have slightly different physiological and functional characteristics to a pure PSE. As in other studies, compared with participants with confluent and advanced destructive emphysema, participants with PSE had better pulmonary function, and higher lung density.

COPD is associated with chronic inflammation, which is characterized by innate and adaptive immune responses (24). The infiltration of macrophages, CD4+-, CD8+- and B-cell lymphocytes were associated with the centrilobular phenotype of emphysematous destruction in COPD (25). On the other hand, histological evaluation showed that neutrophil infiltration is increased in PSE regions (26). In our study, the PSE subtype showed the highest neutrophil count among the subtypes. Neutrophils secrete neutrophil elastase, cathepsin G, and proteinase-3, which lead to alveolar destruction (24). Neutrophilia and sputum neutrophil count correlated with the rate of decline in lung function and peripheral airway dysfunction on CT (27). Tanabe et al. (23) explained that toxic gas spreads to the periphery of the lobule leading to PSE, which leads to an increase in neutrophils. In addition to neutrophils, pro-inflammatory cytokines also play an important role in the pathogenesis of COPD (28). In the present study, the PSE subtype showed significantly higher IL-6 levels; IL-6 increases alveolar cell apoptosis and promotes pulmonary emphysema (29).

Our study showed that functional small airway disease in the trace emphysema subtype had higher degree than in the mild emphysema subtype. Airflow limitation occurred in small airways less than 2 mm in diameter. Since these small airways are beyond the current CT imaging resolution, it is assumed that air trapping is caused by small airway disease (30). Several studies have shown that the narrowing and disappearance of small conducting airways precede emphysematous destruction (31,32). PRM is significant diagnostic value of CT quantification because of its ability to differentiation of emphysema from non-emphysematous air trapping (function small airway disease, PRMfSAD) within lung parenhcyma (19). Also, several studies have shown that PRM are valuable CT derived biomarker assessment of small airway disease and COPD regardless of subtypes of COPD (19,20). Our study using PRM demonstrate that functional small airway disease occurs in the trace emphysema subtype, in the early stage of the disease, and that emphysema increased subsequently, supporting previous studies that indicated that small airway disease is a precursor to emphysema development (30).

Our study showed that participants with bronchial disease tended to have a higher BMI, airway thickness, and lung density, and better lung function. BMI is associated with the severity of emphysema: as BMI lowered, emphysema worsened (33). In contrast, BMI and airway thickness were positively correlated (34). In our study, the bronchial subtype had the highest BMI among the groups, except for the normal group. Due to airway inflammation, the bronchial subtype had higher value of inflammatory markers, including WBC and neutrophil count, and IL-6 and IL-8 levels.

This study had several limitations. First, the number of participants was relatively small for each subtype. In particular, there were few subjects who had confluent centrilobular emphysema and advanced destructive emphysema, and the 2 subtypes had to be combined for analysis. Second, as we did not externally validate our results, a prospective study with a large population is needed. Third, we classified PSE subtype if substantial PSE was present, regardless of other CT features, and thus our PSE subtype groups may not necessarily reflect the characteristics of pure PSE. Fourth, in our study, COPD were defined according to the chronic exposure to noxious agents, chronic respiratory symptoms, postbronchodilator ratio of FEV1/FVC of <0.7, and no serious structural change. Smoking is a significant factor that causes COPD, but various factors other than smoking cause COPD. Our study included patients with non-smoker COPD, so there may be limitation in generalizing smoker COPD. Therefore, further validation is needed.

In conclusion, we demonstrate that visual evaluation of CT images can provide information on lung structure changes and COPD subtypes that cannot currently be identified by quantitative CT. The PSE subtype had better lung function than the severe centrilobular emphysema subtype, and had higher value of inflammatory markers, including neutrophil counts and IL-6 levels. Small airway disease occurs early in this disease, before the development of significant emphysema. Thus, CT-based, visually defined subtypes of COPD can be useful for distinguishing the heterogeneity and pathological processes of COPD.

The article’s supplementary files as