Computed tomographic measurements of airway dimensions and emphysema in smokers. Correlation with lung function

Role of medical and molecular imaging in COPD

Chronic obstructive pulmonary disease (COPD) is expected to climb on the podium of the leading causes of mortality worldwide in the upcoming decade. Clinical diagnosis of COPD has classically relied upon detecting irreversible airflow obstruction on pulmonary function testing as a global assessment of pulmonary physiology. However, the outcome is still not favorable to decrease mortality due to COPD. Progress made in both medical and molecular imaging fields are beginning to offer additional tools to address this clinical problem. This review aims to describe medical and molecular imaging modalities used to diagnose COPD and to select patients for appropriate treatments and to monitor response to therapy.

Quantitative CT analysis of honeycombing area predicts mortality in idiopathic pulmonary fibrosis with definite usual interstitial pneumonia pattern: A retrospective cohort study

Background

Honeycombing on high-resolution computed tomography (HRCT) images is a key finding in idiopathic pulmonary fibrosis (IPF). In IPF, honeycombing area determined by quantitative CT analysis is correlated with pulmonary function test findings. We hypothesized that quantitative CT-derived honeycombing area (HA) might predict mortality in patients with IPF.

Materials and methods

Chest HRCT images of 52 IPF patients with definite usual interstitial pneumonia (UIP) pattern were retrospectively evaluated. Mortality data up to July 31, 2016, were recorded. Using a computer-aided system, HA and percentage of HA (%HA) were measured quantitatively. Predictors of 3-year mortality were evaluated using logistic regression models.

Results

The median %HA, %predicted forced vital capacity (FVC) and composite physiologic index (CPI) were 3.8%, 83.6%, and 33.6, respectively. According to GAP (gender, age, and physiology) stage, 20, 14, and 5 patients were classified under stages I-II-III, respectively. Percentage of HA was significantly correlated with %FVC, CPI, and GAP stage (all, p < 0.001). In univariate analysis, %HA, %FVC, and CPI were statistically significant predictors of mortality. In multivariate analysis using the stepwise regression method, only %HA (odds ratio [OR], 1.27; p = 0.011) was a significant independent predictors of mortality. Patients with %HA ≥ 4.8% had significantly lower survival rates than those with lesser %HA (median survival time, 1.3 vs 5.0 years; log-rank test; p < 0.001).

Conclusion

Quantitative CT-derived HA might be an important and independent predictor of mortality in IPF patients with definite UIP pattern.

Predictive and modifying factors of bone mineral density decline in patients with COPD

RATIONALE

Various determinants of osteoporosis have been previously identified. However, only a few longitudinal studies have examined related factors. We aimed to investigate factors predicting and modifying rapid decline of bone mineral density in patients with chronic obstructive pulmonary disease.

METHODS

We analyzed patients with chronic obstructive pulmonary disease whose bone mineral density were measured at least three times over three years (n = 111). We divided annual per cent changes of bone mineral density in different body parts into tertiles. Rapid decliners (n = 33) were defined as those with the largest decline in at least two parts; all other participants were defined as non-rapid decliners (n = 78).

RESULTS

At enrollment, bone mineral density did not differ between the two groups. However, rapid decliners had a significantly greater rate of new vertebral fractures over 3 years compared with non-rapid decliners. On multivariate logistic regression analysis, age, moderate to severe emphysema, no daily exercise habits, and anemia increased the likelihood of rapid decliners. Furthermore, patients who newly started and continued bisphosphonate exhibited higher annual per cent changes of bone mineral density than did those without bisphosphonate use.

CONCLUSIONS

A rapid decline in bone mineral density correlates to a higher likelihood of vertebral fracture. We clarified the predictors of bone mineral density decline and demonstrated that bisphosphonate use might modify bone mineral density in patients with chronic obstructive pulmonary disease.

Uncontrolled asthma phenotypes defined from parameters using quantitative CT analysis

OBJECTIVE

Asthma is a heterogeneous disease with diverse clinical phenotypes that have been identified via cluster analyses. However, the classification of phenotypes based on quantitative CT (qCT) is poorly understood. The study was conducted to investigate CT determination of uncontrolled asthma phenotypes.

METHODS

Sixty-five patients with uncontrolled asthma (37 with severe asthma, 28 with non-severe asthma) underwent detailed clinical, laboratory, and pulmonary function tests, as well as qCT analysis. Twenty-five healthy subjects were also included in this study and underwent clinical physical examinations, pulmonary function tests, and low-dose CT scans.

RESULTS

The mean lumen area/body surface area ratio was smaller in patients with severe uncontrolled asthma compared with that in healthy subjects (9.84 mm² [SD, 2.57 mm²], 11.96 mm² [SD, 3.09 mm²]; p = 0.026). However, the percentage of mean wall area (WA) was greater (64.39% [SD, 2.55%], 62.09% [SD, 3.81%], p = 0.011). Air trapping (measured based on mean lung density and VI_-856 [%] on expiratory scan) was greater in patients with severe uncontrolled asthma than in those with non-severe uncontrolled asthma and was higher in all patients with uncontrolled asthma than that in healthy subjects (all p < 0.001). Three CT-determined uncontrolled asthma phenotypes were identified. Cluster 1 had mild air trapping with or without proximal airway remodeling. Cluster 2 had moderate air trapping with or without proximal airway remodeling. Cluster 3 had severe air trapping with proximal airway remodeling.

CONCLUSIONS

There was obvious air trapping and proximal airway remodeling in patients with severe uncontrolled asthma. The three CT-determined uncontrolled asthma phenotypes might reflect underlying mechanisms of disease in patient stratification and in the different stages of disease development.

KEY POINTS

• Obvious air trapping and proximal airway remodeling were present in patients with severe uncontrolled asthma. • CT air trapping indices showed a good correlation with disease duration, total IgE, atopy, and OCS and ICS doses, and were even more strongly correlated with clinical lung function. • Three CT-determined uncontrolled asthma phenotypes were identified, which might reflect underlying mechanisms of disease in patient stratification and in the different stages of disease development.

Oscillometry and pulmonary magnetic resonance imaging in asthma and COPD

Developed over six decades ago, pulmonary oscillometry has re-emerged as a noninvasive and effort-independent method for evaluating respiratory-system impedance in patients with obstructive lung disease. Here, we evaluated the relationships between hyperpolarized 3 He ventilation-defect-percent (VDP) and respiratory-system resistance, reactance and reactance area (AX ) measurements in 175 participants including 42 never-smokers without respiratory disease, 56 ex-smokers with chronic-obstructive-pulmonary-disease (COPD), 28 ex-smokers without COPD and 49 asthmatic never-smokers. COPD participants were dichotomized based on x-ray computed-tomography (CT) evidence of emphysema (relative-area CT-density-histogram ≤ 950HU (RA950 ) ≥ 6.8%). In asthma and COPD subgroups, MRI VDP was significantly related to the frequency-dependence of resistance (R5-19 ; asthma: ρ = 0.48, P = 0.0005; COPD: ρ = 0.45, P = 0.0004), reactance at 5 Hz (X5 : asthma, ρ = -0.41, P = 0.004; COPD: ρ = -0.38, P = 0.004) and AX (asthma: ρ = 0.47, P = 0.0007; COPD: ρ = 0.43, P = 0.0009). MRI VDP was also significantly related to R5-19 in COPD participants without emphysema (ρ = 0.54, P = 0.008), and to X5 in COPD participants with emphysema (ρ = -0.36, P = 0.04). AX was weakly related to VDP in asthma (ρ = 0.47, P = 0.0007) and COPD participants with (ρ = 0.39, P = 0.02) and without (ρ = 0.43, P = 0.04) emphysema. AX is sensitive to obstruction but not specific to the type of obstruction, whereas the different relationships for MRI VDP with R5-19 and X5 may reflect the different airway and parenchymal disease-specific biomechanical abnormalities that lead to ventilation defects.

Emphysema extent on computed tomography is a highly specific index in diagnosing persistent airflow limitation: a real-world study in China

Objective

The diagnostic value of emphysema extent in consistent air flow limitation remains controversial. Therefore, we aimed to assess the value of emphysema extent on computed tomography (CT) on the diagnosis of persistent airflow limitation. Furthermore, we developed a diagnostic criterion for further verification.

Materials and methods

We retrospectively enrolled patients who underwent chest CT and lung function test. To be specific, 671 patients were enrolled in the derivation group (Group 1.1), while 479 patients were in the internal validation group (Group 1.2). The percentage of lung volume occupied by low attenuation areas (LAA%) and the percentile of the histogram of attenuation values were calculated.

Results

In patients with persistent airflow limitation, the LAA% was higher and the percentile of the histogram of attenuation values was lower, compared with patients without persistent airflow limitation. Using LAA% with a threshold of −950 HU >1.4% as the criterion, the sensitivity was 44.3% and 47.2%, and the specificity was 95.2% and 95.7%, in Group 1.1 and Group 1.2, respectively. The specificity was influenced by the coexistence of interstitial lung disease, pneumothorax, and post-surgery, rather than the coexistence of pneumonia, nodule, or mass. Multivariable models were also developed.

Conclusion

The emphysema extent on CT is a highly specific marker in the diagnosis of persistent airflow limitation.

Centrilobular emphysema and coronary artery calcification: mediation analysis in the SPIROMICS cohort

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is associated with a two-to-five fold increase in the risk of coronary artery disease independent of shared risk factors. This association is hypothesized to be mediated by systemic inflammation but this link has not been established.

METHODS

We included 300 participants enrolled in the SPIROMICS cohort, 75 each of lifetime non-smokers, smokers without airflow obstruction, mild-moderate COPD, and severe-very severe COPD. We quantified emphysema and airway disease on computed tomography, characterized visual emphysema subtypes (centrilobular and paraseptal) and airway disease, and used the Weston visual score to quantify coronary artery calcification (CAC). We used the Sobel test to determine whether markers of systemic inflammation mediated a link between spirometric and radiographic features of COPD and CAC.

RESULTS

FEV1/FVC but not quantitative emphysema or airway wall thickening was associated with CAC (p = 0.036), after adjustment for demographics, diabetes mellitus, hypertension, statin use, and CT scanner type. To explain this discordance, we examined visual subtypes of emphysema and airway disease, and found that centrilobular emphysema but not paraseptal emphysema or bronchial thickening was independently associated with CAC (p = 0.019). MMP3, VCAM1, CXCL5 and CXCL9 mediated 8, 8, 7 and 16% of the association between FEV1/FVC and CAC, respectively. Similar biomarkers partially mediated the association between centrilobular emphysema and CAC.

CONCLUSIONS

The association between airflow obstruction and coronary calcification is driven primarily by the centrilobular subtype of emphysema, and is linked through bioactive molecules implicated in the pathogenesis of atherosclerosis.

TRIAL REGISTRATION

ClinicalTrials.gov: Identifier: NCT01969344 .

Longitudinal airway remodeling in active and past smokers in a lung cancer screening population

OBJECTIVES

To longitudinally investigate smoking cessation-related changes of quantitative computed tomography (QCT)-based airway metrics in a group of heavy smokers.

METHODS

CT scans were acquired in a lung cancer screening population over 4 years at 12-month intervals in 284 long-term ex-smokers (ES), 405 continuously active smokers (CS), and 31 subjects who quitted smoking within 2 years after baseline CT (recent quitters, RQ). Total diameter (TD), lumen area (LA), and wall percentage (WP) of 1st-8th generation airways were computed using airway analysis software. Inter-group comparison was performed using Mann-Whitney U test or Student's t test (two groups), and ANOVA or ANOVA on ranks with Dunn's multiple comparison test (more than two groups), while Fisher's exact test or chi-squared test was used for categorical data. Multiple linear regression was used for multivariable analysis.

RESULTS

At any time, TD and LA were significantly higher in ES than CS, for example, in 5th-8th generation airways at baseline with 6.24 mm vs. 5.93 mm (p < 0.001) and 15.23 mm² vs. 13.51 mm² (p < 0.001), respectively. RQ showed higher TD (6.15 mm vs. 5.93 mm, n.s.) and significantly higher LA (14.77 mm² vs. 13.51 mm², p < 0.001) than CS after 3 years, and after 4 years. In multivariate analyses, smoking status independently predicted TD, LA, and WP at baseline, at 3 years and 4 years (p < 0.01-0.001), with stronger impact than pack years.

CONCLUSIONS

Bronchial dimensions depend on the smoking status. Smoking-induced airway remodeling can be partially reversible after smoking cessation even in long-term heavy smokers. Therefore, QCT-based airway metrics in clinical trials should consider the current smoking status besides pack years.

KEY POINTS

• Airway lumen and diameter are decreased in active smokers compared to ex-smokers, and there is a trend towards increased airway wall thickness in active smokers. • Smoking-related airway changes improve within 2 years after smoking cessation. • Smoking status is an independent predictor of airway dimensions.

Variability of CT Airways Measurements in COPD Patients Between Morning and Afternoon: Comparisons to Variability of Spirometric Measurements

RATIONALE AND OBJECTIVES

Computed tomography (CT) airways measurements can be used as surrogates to spirometric measurements for assessing bronchodilation in a particular patient with chronic obstructive pulmonary disease. Although spirometric measurements show variations within the opening hours of a hospital department, we aimed to compare the variability of CT airways measurements between morning and afternoon in patients with chronic obstructive pulmonary disease to that of spirometric measurements.

MATERIALS AND METHODS

Twenty patients had pulmonary function tests and CT around 8 am and 4 pm. Luminal area (LA) and wall thickness (WT) of third and fourth generation airways were measured twice by three readers. The percentage of airway area occupied by the wall (WA%) and the square root of wall area at an internal perimeter of 10 mm (√WAPi10) were calculated. The effects of examination time, reader, and measurement session on CT airways measurements were assessed, and the variability of these measurements was compared to that of spirometric measurements.

RESULTS

Variability of LA^3rd and LA^4th was greater than that of spirometric measurements (P values ranging from <.001 to .033). There was no examination time effect on √WAPi10, WT^3rd, LA^4th, or WA%^4th (P values ranging from .102 to .712). There was a reader effect on all CT airways measurements (P values ranging from <.001 to .028), except in WT^3rd (P> .999). There was no effect of measurement session on any CT airway measurement (P values ranging from .535 to >.999).

CONCLUSION

As the variability of LA^3rd and LA^4th is greater than that of spirometric measurements, clinical studies should include cohorts with larger numbers of patients when considering LA than when considering spirometric measurements as end points.

Blood eosinophil count as a prognostic biomarker in COPD

Background

High blood eosinophil count is a predictive biomarker for response to inhaled corticosteroids in prevention of acute exacerbation of COPD, and low blood eosinophil count is associated with pneumonia risk in COPD patients taking inhaled corticosteroids. However, the prognostic role of blood eosinophil count remains underexplored. Therefore, we investigated the associated factors and mortality based on blood eosinophil count in COPD.

Methods

Patients with COPD were recruited from 16 hospitals of the Korean Obstructive Lung Disease cohort (n=395) and COPD in Dusty Area cohort (n=234) of Kangwon University Hospital. The two merged cohorts were divided based on blood eosinophil count into three groups: high (≥5%), middle (2%-5%), and low (<2%).

Results

The high group had longer six-minute walk distance (high =445.8±81.4, middle =428.5±88.0, and low =414.7±86.3 m), higher body mass index (23.3±3.1, 23.1±3.1, and 22.5±3.2 kg/m²), lower emphysema index (18.5±14.1, 22.2±15.3, and 23.7±16.3), and higher inspiratory capacity/total lung capacity ratio (32.6±7.4, 32.4±9.2, and 29.9% ± 8.9%) (P<0.05). The survival period increased with increasing blood eosinophil count (high =9.52±0.23, middle =8.47±1.94, and low =7.42±0.27 years, P<0.05). Multivariate linear regression analysis revealed that the emphysema index was independently and negatively correlated with blood eosinophil count (P<0.05).

Conclusion

In COPD, the severity of emphysema was independently linked with low blood eosinophil count and the longer survival period was associated with increased blood eosinophil count, though it was not proven in the multivariate analysis.

Airway fractal dimension predicts respiratory morbidity and mortality in COPD

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is characterized by airway remodeling. Characterization of airway changes on computed tomography has been challenging due to the complexity of the recurring branching patterns, and this can be better measured using fractal dimensions.

METHODS

We analyzed segmented airway trees of 8,135 participants enrolled in the COPDGene cohort. The fractal complexity of the segmented airway tree was measured by the Airway Fractal Dimension (AFD) using the Minkowski-Bougliand box-counting dimension. We examined associations between AFD and lung function and respiratory morbidity using multivariable regression analyses. We further estimated the extent of peribronchial emphysema (%) within 5 mm of the airway tree, as this is likely to affect AFD. We classified participants into 4 groups based on median AFD, percentage of peribronchial emphysema, and estimated survival.

RESULTS

AFD was significantly associated with forced expiratory volume in one second (FEV1; P < 0.001) and FEV1/forced vital capacity (FEV1/FVC; P < 0.001) after adjusting for age, race, sex, smoking status, pack-years of smoking, BMI, CT emphysema, air trapping, airway thickness, and CT scanner type. On multivariable analysis, AFD was also associated with respiratory quality of life and 6-minute walk distance, as well as exacerbations, lung function decline, and mortality on longitudinal follow-up. We identified a subset of participants with AFD below the median and peribronchial emphysema above the median who had worse survival compared with participants with high AFD and low peribronchial emphysema (adjusted hazards ratio [HR]: 2.72; 95% CI: 2.20-3.35; P < 0.001), a substantial number of whom were not identified by traditional spirometry severity grades.

CONCLUSION

Airway fractal dimension as a measure of airway branching complexity and remodeling in smokers is associated with respiratory morbidity and lung function change, offers prognostic information additional to traditional CT measures of airway wall thickness, and can be used to estimate mortality risk.

TRIAL REGISTRATION

ClinicalTrials.gov identifier: NCT00608764.

FUNDING

This study was supported by NIH K23 HL133438 (SPB) and the COPDGene study (NIH Grant Numbers R01 HL089897 and R01 HL089856). The COPDGene project is also supported by the COPD Foundation through contributions made to an Industry Advisory Board comprised of AstraZeneca, Boehringer Ingelheim, Novartis, Pfizer, Siemens, Sunovion and GlaxoSmithKline.

Quantitative measurement of airway dimensions using ultra-high resolution computed tomography

BACKGROUND

Quantitative measurement of airway dimensions using computed tomography (CT) is performed in relatively larger airways due to the limited resolution of CT scans. Nevertheless, the small airway is an important pathological lesion in lung diseases such as chronic obstructive pulmonary disease (COPD) and asthma. Ultra-high resolution scanning may resolve the smaller airway, but its accuracy and limitations are unclear.

METHODS

Phantom tubes were imaged using conventional (512 × 512) and ultra-high resolution (1024 × 1024 and 2048 × 2048) scans. Reconstructions were performed using the forward-projected model-based iterative reconstruction solution (FIRST) algorithm in 512 × 512 and 1024 × 1024 matrix scans and the adaptive iterative dose reduction 3D (AIDR-3D) algorithm for all scans. In seven subjects with COPD, the airway dimensions were measured using the 1024 × 1024 and 512 × 512 matrix scans.

RESULTS

Compared to the conventional 512 × 512 scan, variations in the CT values for air were increased in the ultra-high resolution scans, except in the 1024×1024 scan reconstructed through FIRST. The measurement error of the lumen area of the tube with 2-mm diameter and 0.5-mm wall thickness (WT) was minimal in the ultra-high resolution scans, but not in the conventional 512 × 512 scan. In contrast to the conventional scans, the ultra-high resolution scans resolved the phantom tube with ≥ 0.6-mm WT at an error rate of < 11%. In seven subjects with COPD, the WT showed a lower value with the 1024 × 1024 scans versus the 512 × 512 scans.

CONCLUSIONS

The ultra-high resolution scan may allow more accurate measurement of the bronchioles with smaller dimensions compared with the conventional scan.

Automatic emphysema detection using weakly labeled HRCT lung images

PURPOSE

A method for automatically quantifying emphysema regions using High-Resolution Computed Tomography (HRCT) scans of patients with chronic obstructive pulmonary disease (COPD) that does not require manually annotated scans for training is presented.

METHODS

HRCT scans of controls and of COPD patients with diverse disease severity are acquired at two different centers. Textural features from co-occurrence matrices and Gaussian filter banks are used to characterize the lung parenchyma in the scans. Two robust versions of multiple instance learning (MIL) classifiers that can handle weakly labeled data, miSVM and MILES, are investigated. Weak labels give information relative to the emphysema without indicating the location of the lesions. The classifiers are trained with the weak labels extracted from the forced expiratory volume in one minute (FEV1) and diffusing capacity of the lungs for carbon monoxide (DLCO). At test time, the classifiers output a patient label indicating overall COPD diagnosis and local labels indicating the presence of emphysema. The classifier performance is compared with manual annotations made by two radiologists, a classical density based method, and pulmonary function tests (PFTs).

RESULTS

The miSVM classifier performed better than MILES on both patient and emphysema classification. The classifier has a stronger correlation with PFT than the density based method, the percentage of emphysema in the intersection of annotations from both radiologists, and the percentage of emphysema annotated by one of the radiologists. The correlation between the classifier and the PFT is only outperformed by the second radiologist.

CONCLUSIONS

The presented method uses MIL classifiers to automatically identify emphysema regions in HRCT scans. Furthermore, this approach has been demonstrated to correlate better with DLCO than a classical density based method or a radiologist, which is known to be affected in emphysema. Therefore, it is relevant to facilitate assessment of emphysema and to reduce inter-observer variability.

Pre-existing pulmonary fibrosis is a risk factor for anti-PD-1-related pneumonitis in patients with non-small cell lung cancer: A retrospective analysis

OBJECTIVES

Pneumonitis related to the use of anti-programmed death 1 (PD-1) antibodies is a common immune-related adverse event that can be life-threatening. However, the relationship between pulmonary fibrosis/emphysema and the incidence of anti-PD-1-related pneumonitis is unclear.

MATERIALS AND METHODS

We retrospectively reviewed data from 123 patients who were diagnosed with non-small-cell lung cancer and treated with anti-PD-1 antibodies (nivolumab or pembrolizumab) at the Aichi Cancer Center Hospital, Japan, between December 17, 2015, and November 30, 2017. Patients who previously received thoracic radiotherapy to the primary lesion, mediastinum, spinal, or rib metastases were excluded from the analysis. Fibrosis score (0-5) and emphysema score (0-4) on baseline chest computed tomography (CT) were determined by two diagnostic radiologists.

RESULTS

Eighteen patients (14.6%) experienced anti-PD-1-related pneumonitis, of which four (3.3%) were grade ≥3. The median onset time of pneumonitis after starting anti-PD-1 therapy was 60 days (range, 6-634 days). According to the analysis of fibrosis score, pneumonitis occurred in 13 (35.1%) of the 37 patients with a fibrosis score ≥1 and in 5 (5.8%) of 86 patients with a fibrosis score of 0. Univariate and multivariate logistic regression analysis revealed that fibrosis score ≥1 was the only risk factor for anti-PD-1-related pneumonitis (p = 0.0008).

CONCLUSION

Our results indicate that pre-existing pulmonary fibrosis significantly increases the risk of anti-PD-1-related pneumonitis. Further studies are needed to identify predictive factors of anti-PD-1-related pneumonitis in patients with fibrotic changes on CT findings.

Comparison between montelukast and tiotropium as add-on therapy to inhaled corticosteroids plus a long-acting β2-agonist in for patients with asthma

Objective: Asthma often remains uncontrolled despite treatment with inhaled corticosteroids (ICS) alone or with ICS plus a long-acting β₂-agonist (LABA). The recommended alternative is the addition of either montelukast or tiotropium. The aim of this study was to compare the effects of montelukast and tiotropium on airway inflammation and remodeling in persistent asthma. Methods: Eighty-seven patients with asthma were treated with budesonide and formoterol (640/18 μg); then, the patients were randomly allocated to three groups to receive oral montelukast (10 mg/day), inhaled tiotropium (5 μg/day), or no add-on to the maintenance therapy for 48 weeks. Fractional exhaled nitric oxide (FeNO) and pulmonary function were measured, and quantitative computed tomography was performed. Results: Compared to the maintenance therapy, add-on montelukast significantly decreased FeNO (p < 0.05) and improved airflow obstruction (p < 0.05), whereas airway dimensions remained unchanged. Changes in FeNO were significantly correlated with changes in FEV₁ (r = -0.71, p < 0.001). In contrast, the addition of tiotropium significantly decreased airway wall area corrected for body surface area (WA/BSA) (p < 0.05), decreased wall thickness (T/√BSA) (p < 0.05) and improved airflow obstruction (p < 0.05) with no change in FeNO. Changes in WA/BSA and T/√BSA were significantly correlated with the change in percentage predicted FEV₁ (r = -0.84, p < 0.001 and r = -0.59, p < 0.01, respectively). Conclusions:Adding either montelukast or tiotropium to ICS/LABA may provide additive benefits with respect to the pulmonary function and airway inflammation or remodeling in patients with asthma.

Accurate Measurement of Airway Morphology on Chest CT Images

In recent years, the ability to accurately measuring and analyzing the morphology of small pulmonary structures on chest CT images, such as airways, is becoming of great interest in the scientific community. As an example, in COPD the smaller conducting airways are the primary site of increased resistance in COPD, while small changes in airway segments can identify early stages of bronchiectasis. To date, different methods have been proposed to measure airway wall thickness and airway lumen, but traditional algorithms are often limited due to resolution and artifacts in the CT image. In this work, we propose a Convolutional Neural Regressor (CNR) to perform cross-sectional measurements of airways, considering wall thickness and airway lumen at once. To train the networks, we developed a generative synthetic model of airways that we refined using a Simulated and Unsupervised Generative Adversarial Network (SimGAN). We evaluated the proposed method by first computing the relative error on a dataset of synthetic images refined with SimGAN, in comparison with other methods. Then, due to the high complexity to create an in-vivo ground-truth, we performed a validation on an airway phantom constructed to have airways of different sizes. Finally, we carried out an indirect validation analyzing the correlation between the percentage of the predicted forced expiratory volume in one second (FEV1%) and the value of the Pi10 parameter. As shown by the results, the proposed approach paves the way for the use of CNNs to precisely and accurately measure small lung airways with high accuracy.

Complementary regional heterogeneity information from COPD patients obtained using oxygen-enhanced MRI and chest CT

Background

The heterogeneous distribution of emphysema is a key feature of chronic obstructive pulmonary disease (COPD) patients that typically is evaluated using high-resolution chest computed tomography (HRCT). Oxygen-enhanced pulmonary magnetic resonance imaging (OEMRI) is a new method to obtain information regarding regional ventilation, diffusion, and perfusion in the lung without radiation exposure. We aimed to compare OEMRI with HRCT for the assessment of heterogeneity in COPD patients.

Methods

Forty patients with stable COPD underwent quantitative HRCT, OEMRI, and pulmonary function tests, including arterial blood gas analysis. OEMRI was also performed on nine healthy control subjects. We measured the severity of emphysema (percent low attenuation volume; LAV%) in whole lungs and the standard deviations (SDs) of the LAV% values of 10 isovolumetric partitions (SD-LAV) as an index of cranial-caudal heterogeneity. Similarly, relative enhancement ratios of oxygen (RERs) in whole lungs from OEMRI and SD-RER were analyzed.

Results

COPD patients showed a lower mean RER than control subjects (12.6% vs 22.0%, p<0.01). The regional heterogeneity of the RERs was not always consistent with the LAV distribution. Both the HRCT (LAV% and SD-LAV) and the OEMRI (RER and SD-RER) indices were significantly associated with the diffusion capacity (DL_CO) and partial pressure of oxygen in arterial blood (PaO₂). The PaO₂ was associated only with the heterogeneity index of HRCT (SD-LAV) (R² = 0.39); however, the PaO₂ was associated with both the mean RER and heterogeneity (SD-RER) in the multivariate analysis (R² = 0.38).

Conclusions

OEMRI-derived parameters were directly associated with oxygen uptake in COPD patients. Although the OEMRI-derived parameters were not identical to the HRCT-derived parameters, the cranial-caudal heterogeneity in HRCT or OEMRI was complementary to that in evaluations of oxygen uptake in the lungs. Functional imaging seems to provide new insights into COPD pathophysiology without radiation exposure.

Fractal analysis of low attenuation clusters on computed tomography in chronic obstructive pulmonary disease

Background

The fractal dimension characterizing the cumulative size distribution of low attenuation area (LAA) clusters, identified with a fixed threshold such as − 950 Hounsfield Units (HU), on computed tomography (CT) sensitively detects parenchymal destruction in chronic obstructive pulmonary disease (COPD) even when the percent LAA (LAA%), a standard emphysema index, is unchanged. This study examines whether the cumulative size distribution of LAA clusters, defined with thresholds of the 15th, 25th, and 35th percentiles of a CT density histogram instead of the fixed-threshold of − 950 HU, exhibits a fractal property and whether its fractal dimension (D’15, D’25, and D’35, respectively) provides additional structural information in emphysematous lungs that is difficult to detect with the conventional − 950-HU-based fractal dimension (D950).

Methods

Chest inspiratory CT scans and pulmonary functions were cross-sectionally examined in 170 COPD subjects. A proxy for the inspiration level at CT scan was obtained by dividing CT-measured total lung volume (CT-TLV) by physiologically measured total lung capacity. Moreover, long-term (> 5 years) changes in D950 and the new fractal dimensions were longitudinally evaluated in 17 current and 42 former smokers with COPD.

Results

D950, but not D’15, D’25, or D’35 was weakly correlated with the proxy for the inspiration. D950, D’25, and D’35 but not D’15 correlated with LAA% and diffusion capacity. In the long-term longitudinal study, LAA% was increased and D950 and D’35 were decreased in both current and former smokers, while D’25 was decreased only in current smokers and D’15 was not changed in either group. The longitudinal changes in D’25 but not those in LAA%, D950, D’15, and D’35 were greater in current smokers than in former smokers. This greater change in D’25 in current smokers was confirmed after adjusting the change in CT-TLV and the baseline D’25.

Conclusions

D’25 reflects diffusion capacity in emphysematous lungs and is robust against inspiration levels during CT scans. This new fractal dimension might provide additional structural information that is difficult to detect with the conventional D950 and LAA% and allow for more sensitive evaluation of emphysema progression over time.

Electronic supplementary material

The online version of this article (10.1186/s12890-018-0714-5) contains supplementary material, which is available to authorized users.

Quantitative CT detects changes in airway dimensions and air-trapping after bronchial thermoplasty for severe asthma

OBJECTIVES

Bronchial thermoplasty (BT) can be considered in the treatment of severe asthma to reduce airway smooth muscle mass and bronchoconstriction. We hypothesized that BT may thus have long-term effects on airway dimensions and air-trapping detectable by quantitative computed tomography (QCT).

METHODS

Paired in- and expiratory CT and inspiratory CT were acquired in 17 patients with severe asthma before and up to two years after bronchial thermoplasty and in 11 additional conservatively treated patients with serve asthma, respectively. A fully automatic software calculated the airways metrics for wall thickness (WT), wall percentage (WP), lumen area (LA) and total diameter (TD). Furthermore, lung air-trapping was quantified by determining the quotient of mean lung attenuation in expiration vs. inspiration (E/I MLA) and relative volume change in the Hounsfield interval -950 to -856 in expiration to inspiration (RVC_856-950) in a generation- and lobe-based approach, respectively.

RESULTS

BT reduced WT for the combined analysis of the 2nd-7th airway generation significantly by 0.06 mm (p = 0.026) and WP by 2.05% (p < 0.001), whereas LA and TD did not change significantly (p = 0.147, p = 0.706). No significant changes were found in the control group. Furthermore, E/I MLA and RVC_856-950 decreased significantly after BT by 12.65% and 1.77% (p < 0.001), respectively.

CONCLUSION

BT significantly reduced airway narrowing and air-trapping in patients with severe asthma. This can be interpreted as direct therapeutic effects caused by a reduction in airway-smooth muscle mass and changes in innervation. A reduction in air-trapping indicates an influence on more peripheral airways not directly treated by the BT procedure.

Prognostic Significance of Large Airway Dimensions on Computed Tomography in the General Population. The Multi-Ethnic Study of Atherosclerosis (MESA) Lung Study

RATIONALE

Large airway dimensions on computed tomography (CT) have been associated with lung function, symptoms, and exacerbations in chronic obstructive pulmonary disease (COPD), as well as with symptoms in smokers with preserved spirometry. Their prognostic significance in persons without lung disease remains undefined.

OBJECTIVES

To examine associations between large airway dimensions on CT and respiratory outcomes in a population-based cohort of adults without prevalent lung disease.

METHODS

The Multi-Ethnic Study of Atherosclerosis recruited participants ages 45-84 years without cardiovascular disease in 2000-2002; we excluded participants with prevalent chronic lower respiratory disease (CLRD). Spirometry was measured in 2004-2006 and 2010-2012. CLRD hospitalizations and deaths were classified by validated criteria through 2014. The average wall thickness for a hypothetical airway of 10-mm lumen perimeter on CT (Pi10) was calculated using measures of airway wall thickness and lumen diameter. Models were adjusted for age, sex, principal components of ancestry, body mass index, smoking, pack-years, scanner, percent emphysema, genetic risk score, and initial forced expiratory volume in 1 second (FEV1) percent predicted.

RESULTS

Greater Pi10 was associated with 9% faster FEV1 decline (95% confidence interval [CI], 2 to 15%; P = 0.012) and increased incident COPD (odds ratio, 2.22; 95% CI, 1.43-3.45; P = 0.0004) per standard deviation among 1,830 participants. Over 78,147 person-years, higher Pi10 was associated with a 57% higher risk of first CLRD hospitalization or mortality (P = 0.0496) per standard deviation. Of Pi10's component measures, both greater airway wall thickness and narrower lumen predicted incident COPD and CLRD clinical events.

CONCLUSIONS

In adults without CLRD, large airway dimensions on CT were prospectively associated with accelerated lung function decline and increased risks of COPD and CLRD hospitalization and mortality.

Relationship of spirometric, body plethysmographic, and diffusing capacity parameters to emphysema scores derived from CT scans

Phenotyping of chronic obstructive pulmonary disease (COPD) with computed tomography (CT) is used to distinguish between emphysema- and airway-dominated type. The phenotype is reflected in correlations with lung function measures. Among these, the relative value of body plethysmography has not been quantified. We addressed this question using CT scans retrospectively collected from clinical routine in a large COPD cohort. Three hundred and thirty five patients with baseline data of the German COPD cohort COPD and Systemic Consequences-Comorbidities Network were included. CT scans were primarily evaluated using a qualitative binary emphysema score. The binary score was positive for emphysema in 52.5% of patients, and there were significant differences between the positive/negative groups regarding forced expiratory volume in 1 second (FEV₁), FEV₁/forced vital capacity (FVC), intrathoracic gas volume (ITGV), residual volume (RV), specific airway resistance (sRaw), transfer coefficient (KCO), transfer factor for carbon monoxide (TLCO), age, pack-years, and body mass index (BMI). Stepwise discriminant analyses revealed the combination of FEV₁/FVC, RV, sRaw, and KCO to be significantly related to the binary emphysema score. The additional positive predictive value of body plethysmography, however, was only slightly higher than that of the conventional combination of spirometry and diffusing capacity, which if taken alone also achieved high predictive values, in contrast to body plethysmography. The additional information on the presence of CT-diagnosed emphysema as conferred by body plethysmography appeared to be minor compared to the well-known combination of spirometry and CO diffusing capacity.

Semi-automatic Methods for Airway and Adjacent Vessel Measurement in Bronchiectasis Patterns in Lung HRCT Images of Cystic Fibrosis Patients

Airway and vessel characterization of bronchiectasis patterns in lung high-resolution computed tomography (HRCT) images of cystic fibrosis (CF) patients is very important to compute the score of disease severity. We propose a hybrid and evolutionary optimized threshold and model-based method for characterization of airway and vessel in lung HRCT images of CF patients. First, the initial model of airway and vessel is obtained using the enhanced threshold-based method. Then, the model is fitted to the actual image by optimizing its parameters using particle swarm optimization (PSO) evolutionary algorithm. The experimental results demonstrated the outperformance of the proposed method over its counterpart in R-squared, mean and variance of error, and run time. Moreover, the proposed method outperformed its counterpart for airway inner diameter/vessel diameter (AID/VD) and airway wall thickness/vessel diameter (AWT/VD) biomarkers in R-squared and slope of regression analysis.

Definition and diagnosis of asthma-COPD overlap (ACO)

It is now widely recognized that asthma and COPD can coexist as asthma-COPD overlap (ACO), but the preliminary attempts at providing universal guidelines for the diagnosis of ACO still need to be improved. We believe that a case can be made for devising guidelines for the diagnosis of this increasingly common disease that are specific to Japan. In this paper, we present our consensus-based description of ACO which we believe is realistic for use in our country. In addition, we cite the scientific evidence for our own "objective" features used to develop the criteria for COPD and asthma diagnosis. We acknowledge that they will need to be validated and updated over time, but hope the results will encourage further research on the characteristics and treatment of this commonly encountered clinical problem.

Characteristics of COPD phenotypes classified according to the findings of HRCT and spirometric indices and its correlation to clinical characteristics

INTRODUCTION

In recent years, there has been increasing interest in diagnosing various components of chronic obstructive pulmonary disease (COPD) using high-resolution computed tomography (HRCT). The present study was undertaken to evaluate HRCT features in patients with COPD.

MATERIALS AND METHODS

Fifty patients of COPD (confirmed on Spirometry as per the GOLD guidelines 2014 guidelines) were enrolled, out of which 35 patients got a HRCT done. The Philips computer program for lung densitometry was used with these limits (-800/-1, 024 Hounsfield unit [HU]) to calculate densities, after validating densitometry values with phantoms. We established the area with a free hand drawing of the region of interest, then we established limits (in HUs) and the computer program calculated the attenuation as mean lung density (MLD) of the lower and upper lobes.

RESULTS

There was a significant correlation between smoking index and anteroposterior tracheal diameter (P = 0.036). Tracheal index was found to be decreasing with increasing disease severity which was statistically significant (P = 0.037). A mild linear correlation of pre-forced expiratory volume in the first second (FEV1) was observed with lower lobe and total average MLD while a mild linear correlation of post-FEV1 was observed with both coronal (P = 0.042) and sagittal (P = 0.001) lower lobes MLD. In addition, there was a linear correlation between both pre (P = 0.050) and post (P = 0.024) FEV1/forced vital capacity with sagittal lower lobe MLD.

CONCLUSION

HRCT may be an important additional tool in the holistic evaluation of COPD.

Heterogeneity in pulmonary emphysema: Analysis of CT attenuation using Gaussian mixture model

PURPOSE

To utilize Gaussian mixture model (GMM) for the quantification of chronic obstructive pulmonary disease (COPD) and to evaluate the combined use of multiple types of quantification.

MATERIALS AND METHODS

Eighty-seven patients (67 men, 20 women; age, 67.4 ± 11.0 years) who had undergone computed tomography (CT) and pulmonary function test (PFT) were included. The heterogeneity of CT attenuation in emphysema (HC) was obtained by analyzing a distribution of CT attenuation with GMM. The percentages of low-attenuation volume in the lungs (LAV), wall area of bronchi (WA), and the cross-sectional area of small pulmonary vessels (CSA) were also calculated. The relationships between COPD quantifications and the PFT results were evaluated by Pearson's correlation coefficients and through linear models, with the best models selected using Akaike information criterion (AIC).

RESULTS

The correlation coefficients with FEV1 were as follows: LAV, -0.505; HC, -0.277; CSA, 0.384; WA, -0.196. The correlation coefficients with FEV1/FVC were: LAV, -0.640; HC, -0.136; CSA, 0.288; WA, -0.131. For predicting FEV1, the smallest AIC values were obtained in the model with LAV, HC, CSA, and WA. For predicting FEV1/FVC, the smallest AIC values were obtained in the model with LAV and HC. In both models, the coefficient of HC was statistically significant (P-values = 0.000880 and 0.0441 for FEV1 and FEV1/FVC, respectively).

CONCLUSION

GMM was applied to COPD quantification. The results of this study show that COPD severity was associated with HC. In addition, it is shown that the combined use of multiple types of quantification made the evaluation of COPD severity more reliable.

Respiratory impedance is correlated with airway narrowing in asthma using three-dimensional computed tomography

BACKGROUND

Respiratory impedance comprises the resistance and reactance of the respiratory system and can provide detailed information on respiratory function. However, details of the relationship between impedance and morphological airway changes in asthma are unknown.

OBJECTIVE

We aimed to evaluate the correlation between imaging-based airway changes and respiratory impedance in patients with asthma.

METHODS

Respiratory impedance and spirometric data were evaluated in 72 patients with asthma and 29 reference subjects. We measured the intraluminal area (Ai) and wall thickness (WT) of third- to sixth-generation bronchi using three-dimensional computed tomographic analyses, and values were adjusted by body surface area (BSA, Ai/BSA, and WT/the square root (√) of BSA).

RESULTS

Asthma patients had significantly increased respiratory impedance, decreased Ai/BSA, and increased WT/√BSA, as was the case in those without airflow limitation as assessed by spirometry. Ai/BSA was inversely correlated with respiratory resistance at 5 Hz (R5) and 20 Hz (R20). R20 had a stronger correlation with Ai/BSA than did R5. Ai/BSA was positively correlated with forced expiratory volume in 1 second/forced vital capacity ratio, percentage predicted forced expiratory volume in 1 second, and percentage predicted mid-expiratory flow. WT/√BSA had no significant correlation with spirometry or respiratory impedance.

CONCLUSIONS & CLINICAL RELEVANCE

Respiratory resistance is associated with airway narrowing.

Morphological and functional properties of the conducting human airways investigated by in vivo computed tomography and in vitro MRI

The accurate representation of the human airway anatomy is crucial for understanding and modeling the structure-function relationship in both healthy and diseased lungs. The present knowledge in this area is based on morphometric studies of excised lung casts, partially complemented by in vivo studies in which computed tomography (CT) was used on a small number of subjects. In the present study, we analyzed CT scans of a cohort of healthy subjects and obtained comprehensive morphometric information down to the seventh generation of bronchial branching, including airway diameter, length, branching angle, and rotation angle. Although some of the geometric parameters (such as the child-to-parent branch diameter ratio) are found to be in line with accepted values, for others (such as the branch length-to-diameter ratio) our findings challenge the common assumptions. We also evaluated several metrics of self-similarity, including the fractal dimension of the airway tree. Additionally, we used phase-contrast magnetic resonance imaging (MRI) to obtain the volumetric flow field in the three-dimensional-printed airway model of one of the subjects during steady inhalation. This is used to relate structural and functional parameters and, in particular, to close the power-law relationship between branch flow rate and diameter. The diameter exponent is found to be significantly lower than in the usually assumed Poiseuille regime, which we attribute to the strong secondary (i.e., transverse) velocity component. The strength of the secondary velocity with respect to the axial component exceeds the levels found in idealized airway models and persists within the first seven generations. NEW & NOTEWORTHY We performed a comprehensive computed tomography-based study of the conductive airway morphology in normal human subjects, including branch diameter, length, and mutual angles. We found significant departure from classic homothetic relationships. We also carried out MRI measurements of the three-dimensional inspiratory flow in an anatomy-based model and directly assessed structure-function relationships that have so far been assumed. We found that strong secondary flows (i.e., transverse velocity components) persist through the first seven generations of bronchial branching.

Emphysematous and Nonemphysematous Gas Trapping in Chronic Obstructive Pulmonary Disease: Quantitative CT Findings and Pulmonary Function

Purpose To identify a prevalent computed tomography (CT) subtype in patients with chronic obstructive pulmonary disease (COPD) by separating emphysematous from nonemphysematous contributions to total gas trapping and to attempt to predict and grade the emphysematous gas trapping by using clinical and functional data. Materials and Methods Two-hundred and two consecutive eligible patients (159 men and 43 women; mean age, 70 years [age range, 41-85 years]) were prospectively studied. Pulmonary function and CT data were acquired by pulmonologists and radiologists. Noncontrast agent-enhanced thoracic CT scans were acquired at full inspiration and expiration, and were quantitatively analyzed by using two software programs. CT parameters were set as follows: 120 kVp; 200 mAs; rotation time, 0.5 second; pitch, 1.1; section thickness, 0.75 mm; and reconstruction kernels, b31f and b70f. Gas trapping obtained by difference of inspiratory and expiratory CT density thresholds (percentage area with CT attenuation values less than -950 HU at inspiration and percentage area with CT attenuation values less than -856 HU at expiration) was compared with that obtained by coregistration analysis. A logistic regression model on the basis of anthropometric and functional data was cross-validated and trained to classify patients with COPD according to the relative contribution of emphysema to total gas trapping, as assessed at CT. Results Gas trapping obtained by difference of inspiratory and expiratory CT density thresholds was highly correlated (r = 0.99) with that obtained by coregistration analysis. Four groups of patients were distinguished according to the prevalent CT subtype: prevalent emphysematous gas trapping, prevalent functional gas trapping, mixed severe, and mixed mild. The predictive model included predicted forced expiratory volume in 1 second/vital capacity, percentage of predicted forced expiratory volume in 1 second, percentage of diffusing capacity for carbon monoxide, and body mass index as emphysema regressors at CT, with 81% overall accuracy in classifying patients according to its extent. Conclusion The relative contribution of emphysematous and nonemphysematous gas trapping obtained by coregistration of inspiratory and expiratory CT scanning can be determined accurately by difference of CT inspiratory and expiratory density thresholds. CT extent of emphysema can be predicted with accuracy suitable for clinical purposes by pulmonary function data and body mass index. ^© RSNA, 2018 Online supplemental material is available for this article.

Relationships between pulmonary function test parameters and quantitative computed tomography measurements of emphysema in subjects with chronic obstructive pulmonary disease

Objective.

CT is able to precisely define the pathological process in COPD. There are a number of previous articles discussing the distribution of emphysema and its connection with pulmonary function tests. However, the results obtained by the researchers are not identical.

Purpose.

To assess relationships between emphysema and pulmonary function test parameters in COPD patients.

Materials and methods.

Fifty-nine patients diagnosed to have COPD underwent chest CT examinations and pulmonary function tests.

For the quantitative assessment, percentages of low attenuation volume LAV _{950 HU} (%) of a both lungs, the right lung, the left lung, and each lobe were obtained. Quantitative CT measurements were compared with forced expiratory volume in 1 s (FEV₁), the ratio of FEV₁ to forced vital capacity (FEV₁/FVC), the diffusing capacity for carbon monoxide (DLco) and total lung capacity (TLC).

Results.

Except for the right middle lobe and the right upper lobe, respectively, all the quantitative CT measurements showed weak to moderate negative correlations with diffusing capacity (DLco) (r = –0.35 to –0.61, p < 0.05) and weak positive correlations with TLC (r = 0.34 to 0.44, p < 0.05). Group analysis indicated that LAV_{–950 HU} (%) values of both lungs, right lung, left lung, and each lobe, except for right middle lobe, were increased in patients with GOLD stages 3 and 4 of COPD compared to GOLD stages 1 and 2 (p < 0.05).

Conclusion.

CT measurements of emphysema are significantly related to pulmonary function tests results, particularly DLco.

Disease Staging and Prognosis in Smokers Using Deep Learning in Chest Computed Tomography

RATIONALE

Deep learning is a powerful tool that may allow for improved outcome prediction.

OBJECTIVES

To determine if deep learning, specifically convolutional neural network (CNN) analysis, could detect and stage chronic obstructive pulmonary disease (COPD) and predict acute respiratory disease (ARD) events and mortality in smokers.

METHODS

A CNN was trained using computed tomography scans from 7,983 COPDGene participants and evaluated using 1,000 nonoverlapping COPDGene participants and 1,672 ECLIPSE participants. Logistic regression (C statistic and the Hosmer-Lemeshow test) was used to assess COPD diagnosis and ARD prediction. Cox regression (C index and the Greenwood-Nam-D'Agnostino test) was used to assess mortality.

MEASUREMENTS AND MAIN RESULTS

In COPDGene, the C statistic for the detection of COPD was 0.856. A total of 51.1% of participants in COPDGene were accurately staged and 74.95% were within one stage. In ECLIPSE, 29.4% were accurately staged and 74.6% were within one stage. In COPDGene and ECLIPSE, the C statistics for ARD events were 0.64 and 0.55, respectively, and the Hosmer-Lemeshow P values were 0.502 and 0.380, respectively, suggesting no evidence of poor calibration. In COPDGene and ECLIPSE, CNN predicted mortality with fair discrimination (C indices, 0.72 and 0.60, respectively), and without evidence of poor calibration (Greenwood-Nam-D'Agnostino P values, 0.307 and 0.331, respectively).

CONCLUSIONS

A deep-learning approach that uses only computed tomography imaging data can identify those smokers who have COPD and predict who are most likely to have ARD events and those with the highest mortality. At a population level CNN analysis may be a powerful tool for risk assessment.

Significance of Low-Attenuation Cluster Analysis on Quantitative CT in the Evaluation of Chronic Obstructive Pulmonary Disease

Objective

To assess clinical feasibility of low-attenuation cluster analysis in evaluation of chronic obstructive pulmonary disease (COPD).

Materials and Methods

Subjects were 199 current and former cigarette smokers that underwent CT for quantification of COPD and had physiological measurements. Quantitative CT (QCT) measurements included low-attenuation area percent (LAA%) (voxels ≤ -950 Hounsfield unit [HU]), and two-dimensional (2D) and three-dimensional D values of cluster analysis at three different thresholds of CT value (-856, -910, and -950 HU). Correlation coefficients between QCT measurements and physiological indices were calculated. Multivariable analyses for percentage of predicted forced expiratory volume at one second (%FEV1) was performed including sex, age, body mass index, LAA%, and D value had the highest correlation coefficient with %FEV1 as independent variables. These analyses were conducted in subjects including those with mild COPD (global initiative of chronic obstructive lung disease stage = 0-II).

Results

LAA% had a higher correlation coefficient (-0.549, p < 0.001) with %FEV1 than D values in subjects while 2D D_-910HU (-0.350, p < 0.001) revealed slightly higher correlation coefficient than LAA% (-0.343, p < 0.001) in subjects with mild COPD. Multivariable analyses revealed that LAA% and 2D D value_-910HU were significant independent predictors of %FEV1 in subjects and that only 2D D value_-910HU revealed a marginal p value (0.05) among independent variables in subjects with mild COPD.

Conclusion

Low-attenuation cluster analysis provides incremental information regarding physiologic severity of COPD, independent of LAA%, especially with mild COPD.

Performance of quantitative CT parameters in assessment of disease severity in COPD: A prospective study

BACKGROUND

Both emphysematous destruction of lung parenchyma and airway remodeling is thought to contribute to airflow limitation in cases of chronic obstructive pulmonary disease (COPD).

OBJECTIVE

To evaluate the value of quantitative computed tomography (QCT) parameters of emphysema and airway disease with disease severity in patients with COPD.

MATERIALS AND METHODS

We prospectively studied 50 patients with COPD, which included nonsmokers and patients with different degrees of cumulative smoking exposure. Three QCT parameters namely LAA% (low attenuation area percentage), WA% (Wall area percentage), and pi10 were calculated as per the standard technique. Forced expiratory volume in 1 s (FEV1), BODE score, and MMRC dyspnea scale were used as measures of disease severity.

RESULTS

FEV1 was inversely and significantly associated with all three QCT parameters. Receiver operated characteristic curves in prediction of GOLD class 3 COPD yielded cut-off values of 12.2, 61.45, and 3.5 for LAA%, WA%, and pi10, respectively, with high sensitivities and specificities. In multiple linear regression model, however, only LAA% proved to be significantly associated with FEV1, BODE, and dyspnea.

CONCLUSION

QCT indices of both emphysema and airway disease influence FEV1, dyspnea, and BODE score in patients with COPD. Emphysema, however, appears to be more closely related to disease severity.

A cluster analysis of chronic obstructive pulmonary disease in dusty areas cohort identified three subgroups

Background

Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease with variable clinical manifestations, structural changes, and treatment responses. In a cohort study, we performed a baseline cluster analysis to identify the subgroups of COPD and to assess the clinical outcomes of each subgroup during a 1-year follow-up.

Methods

We analyzed dusty areas cohort comprising 272 patients with COPD. The main factors with the highest loading in 15 variables were selected using principal component analysis (PCA) at baseline. The COPD patients were classified by hierarchical cluster analysis using clinical, physiological, and imaging data based on PCA-transformed data. The clinical parameters and outcomes during the 1-year follow-up were evaluated among the subgroups.

Results

PCA revealed that six independent components accounted for 77.3% of variance. Three distinct subgroups were identified through the cluster analysis. Subgroup 1 included younger subjects with fewer symptoms and mild airflow obstruction, and they had fewer exacerbations during the 1-year follow-up. Subgroup 2 comprised subjects with additional symptoms and moderate airflow obstruction, and they most frequently experienced exacerbations requiring hospitalization during the 1-year follow-up. Subgroup 3 included subjects with additional symptoms and mild airflow obstruction; this group had more female patients and a modest frequency of exacerbations requiring hospitalization.

Conclusions

Cluster analysis using the baseline data of a COPD cohort identified three distinct subgroups with different clinical parameters and outcomes. These findings suggest that the identified subgroups represent clinically meaningful subtypes of COPD.

Electronic supplementary material

The online version of this article (10.1186/s12890-017-0553-9) contains supplementary material, which is available to authorized users.

Anemia as a clinical marker of stable chronic obstructive pulmonary disease in the Korean obstructive lung disease cohort

Background

Anemia is a major co-morbidity in chronic obstructive pulmonary disease (COPD). However, the mechanism of development for anemia and the impact of anemia on the prognosis of COPD remain poorly understood. Therefore, this study attempted to evaluate the prognostic role of anemia on the clinical course of COPD and investigate the factors linked with the serum hemoglobin level in COPD.

Methods

We analyzed 407 COPD patients enrolled in the Korean obstructive lung disease (KOLD) cohort at 16 hospitals in Korea recruited over 9 years. Multivariate Cox regression analysis was performed to find independent predictors of survival and multivariate logistic regression analyses were done to find independent factors.

Results

Anemic COPD were older with lower body mass index (BMI) (P<0.001), lower serum cholesterol level (P=0.001), lower serum albumin level (P<0.001), and shorter 6-minute walking distance (P=0.046) compared to non-anemic COPD. A multivariate Cox regression analysis revealed that age (P=0.002), BMI (P=0.001), post-bronchodilator forced expiratory volume in 1 second (FEV₁) (P=0.007), 6-minute walk distance (P=0.008), anemia (P=0.025) were significant predictors for all-cause mortality. In multivariate regression analysis, older age (P<0.001), female gender (P=0.001), lower BMI (P=0.016), and lower serum albumin level (P<0.001) were independent factors associated with lower serum hemoglobin level.

Conclusions

Our data showed that anemia was an independent risk factor for mortality in COPD, and aging, lower serum albumin level, and lower BMI were independent factors associated with lower serum hemoglobin level.

Inflammatory biomarkers and radiologic measurements in never-smokers with COPD: A cross-sectional study from the CODA cohort

Various biomarkers have emerged as potential surrogates to represent various subgroups of chronic obstructive pulmonary disease (COPD), which manifest with different phenotypes. However, the biomarkers representing never-smokers with COPD have not yet been well elucidated. The aim of this study was to evaluate the associations of certain serum and radiological biomarkers with the presence of COPD in never-smokers. To explore the associations of serum and radiological biomarkers with the presence of COPD in never-smokers, we conducted a cross-sectional patient cohort study composed of never-smokers from the COPD in Dusty Areas (CODA) cohort, consisting of subjects living in dusty areas near cement plants in South Korea. Of the 131 never-smokers in the cohort, 77 (58.8%) had COPD. There were no significant differences in the number of subjects with high levels of inflammatory biomarkers (>90th percentile of never-smokers without COPD), including white blood cell count, total bilirubin, interleukin (IL)-6, IL-8, and C-reactive protein, or radiologic measurements (including emphysema index and mean wall area percentage) between never-smokers with COPD and those without COPD. However, the number of subjects with high uric acid was significantly higher in never-smokers with COPD than never-smokers without COPD (31.2% (24/77) vs. 11.1% (6/54); p = 0.013). In addition, multivariate analysis revealed that high uric acid was significantly associated with the presence of COPD in never-smokers (adjusted relative risk: 1.63; 95% confidence interval: 1.21, 2.18; p = 0.001). Our study suggests that high serum levels of uric acid might be a potential biomarker for assessing the presence of COPD in never-smokers.

Disease Severity Dependence of the Longitudinal Association Between CT Lung Density and Lung Function in Smokers

BACKGROUND

In smokers, the lung parenchyma is characterized by inflammation and emphysema, processes that can result in local gain and loss of lung tissue. CT measures of lung density might reflect lung tissue changes; however, longitudinal data regarding the effects of CT lung tissue on FEV₁ in smokers with and without COPD are scarce.

METHODS

The 15th percentile of CT lung density was obtained from the scans of 3,390 smokers who completed baseline and 5-year follow-up of the Genetic Epidemiology of COPD (COPDGene) study visits. The longitudinal relationship between total lung capacity-adjusted lung density (TLC-PD15) and FEV₁ was assessed by using multivariable mixed models. Separate models were performed in smokers at risk, smokers with preserved ratio and impaired spirometry (PRISm), and smokers with COPD according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) staging system.

RESULTS

The direction of the relationship between lung density and lung function was GOLD stage dependent. In smokers with PRISm, a 1-g/L decrease in TLC-PD15 was associated with an increase of 2.8 mL FEV₁ (P = .02). In contrast, among smokers with GOLD III to IV COPD, a 1-g/L decrease in TLC-PD15 was associated with a decrease of 4.1 mL FEV₁ (P = .002).

CONCLUSIONS

A decline in TLC-PD15 was associated with an increase or decrease in FEV₁ depending on disease severity. The associations are GOLD stage specific, and their presence might influence the interpretation of future studies that use CT lung density as an intermediate study end point for a decline in lung function.

TRIAL REGISTRY

ClinicalTrials.gov; No.: NCT00608764; URL: www.clinicaltrials.gov.

CT imaging of chronic obstructive pulmonary disease: insights, disappointments, and promise

CT imaging is a readily quantifiable tool that can provide in-vivo assessments of lung structure in conditions such as chronic obstructive pulmonary disease (COPD). The information extracted from these data has been used in many clinical, epidemiological, and genetic investigations for patient stratification and prognostication, and to determine intermediate endpoints for clinical trials. Although these efforts have informed our understanding of the heterogeneity of pulmonary disease in smokers, they have not yet translated into new treatments for COPD or the personalisation of patient care. There are a multitude of potential reasons for this, including the lack of insight that static imaging provides for lung function and dysfunction, the limited resolution of clinical CT scanning for microscopic changes to the lung architecture, and the challenges that the biomedical community faces when trying to translate discovery to therapy. Such limitations might be addressed through novel image analysis techniques, up-and-coming CT-based and MRI-based technologies, closer ties between academia and industry, and an expanded endeavour to share data across the biomedical community.

Evaluating Qualitative and Quantitative Computerized Tomography Indicators of Chronic Obstructive Pulmonary Disease and Their Correlation with Pulmonary Function Tests

BACKGROUND

With increasingly aging populations, chronic obstructive pulmonary disease (COPD) is the fourth leading cause of death today. Emphysematous changes, an important component of the disease, must be determined on HRCT, either qualitatively or quantitatively. The purpose of this study was to evaluate features that help determine emphysematous changes and correlate them with respiratory function tests (RFTs).

MATERIAL/METHODS

A total of thirty COPD patients and a control group of the same size, matched for age, were included in the study. The mean lung parenchyma density values on inspiration and expiration, visual HRCT scores, and pulmonary function tests were obtained. IBM SPSS statistical software (version 22) was used to perform correlation analysis (Pearson's coefficient) and the Mann-Whitney U test.

RESULTS

The most valuable RFTs for determining emphysematous changes were DLCO, FEV1, and FEV1/FVC, in that order. Quantitative measures of the mean lung density had the highest correlation with coefficient on expiration.

CONCLUSIONS

As regards the comparison between objective and subjective density values, the HRCT-based visual density values are satisfactory. On the other hand, the best assessment can be performed with the use of mean density values on expiration. DLCO, FEV1, and FEV1/FVC were found to be valuable parameters in determining parenchymal changes.

Relationship between fraction of exhaled nitric oxide and airway morphology assessed by three-dimensional CT analysis in asthma

Fraction of exhaled nitric oxide (FeNO) provides information about chronic inflammation in asthma. However, its relationship with structural changes in the airways is unknown. We aimed to evaluate the correlation between computer-based airway changes and FeNO in patients with asthma. The wall area (WA) and airway inner luminal area (Ai) of the third- to sixth-generation bronchi were measured using three-dimensional computed tomography in asthmatic patients. Each value was corrected by body surface area (BSA). Relationships between FeNO and WA/BSA and Ai/BSA were evaluated. Forty-one clinically stable patients with asthma were evaluated. FeNO was significantly correlated with WA/BSA of the third-, fourth-, fifth- and sixth-generation bronchi (Spearman correlation coefficient (ρ) = 0.326, p = 0.041; ρ = 0.356, p = 0.025; ρ = 0.496, p = 0.002; and ρ = 0.529, p < 0.001, respectively). The correlation with sixth-generation bronchi was significantly greater than with the third-generation bronchi (p = 0.047). Partial rank correlation analysis indicated FeNO was significantly correlated with WA/BSA of the sixth-generation bronchi, independent from confounding factors of Ai/BSA, age, duration of asthma, dose of inhaled corticosteroid, blood eosinophil percentage, and blood IgE (ρ = 0.360, p = 0.034). In contrast, there was no correlation between FeNO and Ai/BSA. FeNO correlates with bronchial wall thickening in asthma patients. Measurement of FeNO may be useful to detect airway remodeling in asthma.

Imaging approaches to understand disease complexity: chronic obstructive pulmonary disease as a clinical model

The clinical manifestations of chronic obstructive pulmonary disease (COPD) reflect an aggregate of multiple pulmonary and extrapulmonary processes. It is increasingly clear that full assessment of these processes is essential to characterize disease burden and to tailor therapy. Medical imaging has advanced such that it is now possible to obtain in vivo insight in the presence and severity of lung disease-associated features. In this review, we have assembled data from multiple disciplines of medical imaging research to review the role of imaging in characterization of COPD. Topics include imaging of the lungs, body composition, and extrapulmonary tissue metabolism. The primary focus is on imaging modalities that are widely available in clinical care settings and that potentially contribute to describing COPD heterogeneity and enhance our insight in underlying pathophysiological processes and their structural and functional effects.

Micro-Computed Tomography Comparison of Preterminal Bronchioles in Centrilobular and Panlobular Emphysema

RATIONALE

Very little is known about airways that are too small to be visible on thoracic multidetector computed tomography but larger than the terminal bronchioles.

OBJECTIVES

To examine the structure of preterminal bronchioles located one generation proximal to terminal bronchioles in centrilobular and panlobular emphysema.

METHODS

Preterminal bronchioles were identified by backtracking from the terminal bronchioles, and their centerlines were established along the entire length of their lumens. Multiple cross-sectional images perpendicular to the centerline were reconstructed to evaluate the bronchiolar wall and lumen, and the alveolar attachments to the outer airway walls in relation to emphysematous destruction in 28 lung samples from six patients with centrilobular emphysema, 20 lung samples from seven patients with panlobular emphysema associated with alpha-1 antitrypsin deficiency, and 47 samples from seven control (donor) lungs.

MEASUREMENTS AND MAIN RESULTS

The preterminal bronchiolar length, wall volume, total volume (wall + lumen), lumen circularity, and number of alveolar attachments were reduced in both centrilobular and panlobular emphysema compared with control lungs. In contrast, thickening of the wall and narrowing of the lumen were more severe and heterogeneous in centrilobular than in panlobular emphysema. The bronchiolar lumen was narrower in the middle than at both ends, and the decreased number of alveolar attachments was associated with increased wall thickness in centrilobular emphysema.

CONCLUSIONS

These results provide new information about small airways pathology in centrilobular and panlobular emphysema and show that these changes affect airways that are not visible with thoracic multidetector computed tomography scans but located proximal to the terminal bronchioles in chronic obstructive pulmonary disease.

Automated prediction of emphysema visual score using homology-based quantification of low-attenuation lung region

Objective

The purpose of this study was to investigate the relationship between visual score of emphysema and homology-based emphysema quantification (HEQ) and evaluate whether visual score was accurately predicted by machine learning and HEQ.

Materials and methods

A total of 115 anonymized computed tomography images from 39 patients were obtained from a public database. Emphysema quantification of these images was performed by measuring the percentage of low-attenuation lung area (LAA%). The following values related to HEQ were obtained: nb₀ and nb₁. LAA% and HEQ were calculated at various threshold levels ranging from −1000 HU to −700 HU. Spearman’s correlation coefficients between emphysema quantification and visual score were calculated at the various threshold levels. Visual score was predicted by machine learning and emphysema quantification (LAA% or HEQ). Random Forest was used as a machine learning algorithm, and accuracy of prediction was evaluated by leave-one-patient-out cross validation. The difference in the accuracy was assessed using McNemar’s test.

Results

The correlation coefficients between emphysema quantification and visual score were as follows: LAA% (−950 HU), 0.567; LAA% (−910 HU), 0.654; LAA% (−875 HU), 0.704; nb₀ (−950 HU), 0.552; nb₀ (−910 HU), 0.629; nb₀ (−875 HU), 0.473; nb₁ (−950 HU), 0.149; nb₁ (−910 HU), 0.519; and nb₁ (−875 HU), 0.716. The accuracy of prediction was as follows: LAA%, 55.7% and HEQ, 66.1%. The difference in accuracy was statistically significant (p = 0.0290).

Conclusion

LAA% and HEQ at −875 HU showed a stronger correlation with visual score than those at −910 or −950 HU. HEQ was more useful than LAA% for predicting visual score.

Phenotyping emphysema and airways disease: Clinical value of quantitative radiological techniques

The pathophysiology of chronic obstructive pulmonary disease (COPD) and Alpha one antitrypsin deficiency is increasingly recognised as complex such that lung function alone is insufficient for early detection, clinical categorisation and dictating management. Quantitative imaging techniques can detect disease earlier and more accurately, and provide an objective tool to help phenotype patients into predominant airways disease or emphysema. Computed tomography provides detailed information relating to structural and anatomical changes seen in COPD, and magnetic resonance imaging/nuclear imaging gives functional and regional information with regards to ventilation and perfusion. It is likely imaging will become part of routine clinical practice, and an understanding of the implications of the data is essential. This review discusses technical and clinical aspects of quantitative imaging in obstructive airways disease.

Chronic Obstructive Pulmonary Disease: Lobar Analysis with Hyperpolarized 129Xe MR Imaging

Purpose To compare lobar ventilation and apparent diffusion coefficient (ADC) values obtained with hyperpolarized xenon 129 (129Xe) magnetic resonance (MR) imaging to quantitative computed tomography (CT) metrics on a lobar basis and pulmonary function test (PFT) results on a whole-lung basis in patients with chronic obstructive pulmonary disease (COPD). Materials and Methods The study was approved by the National Research Ethics Service Committee; written informed consent was obtained from all patients. Twenty-two patients with COPD (Global Initiative for Chronic Obstructive Lung Disease stage II-IV) underwent hyperpolarized 129Xe MR imaging at 1.5 T, quantitative CT, and PFTs. Whole-lung and lobar 129Xe MR imaging parameters were obtained by using automated segmentation of multisection hyperpolarized 129Xe MR ventilation images and hyperpolarized 129Xe MR diffusion-weighted images after coregistration to CT scans. Whole-lung and lobar quantitative CT-derived metrics for emphysema and bronchial wall thickness were calculated. Pearson correlation coefficients were used to evaluate the relationship between imaging measures and PFT results. Results Percentage ventilated volume and average ADC at lobar 129Xe MR imaging showed correlation with percentage emphysema at lobar quantitative CT (r = -0.32, P < .001 and r = 0.75, P < .0001, respectively). The average ADC at whole-lung 129Xe MR imaging showed moderate correlation with PFT results (percentage predicted transfer factor of the lung for carbon monoxide [Tlco]: r = -0.61, P < .005) and percentage predicted functional residual capacity (r = 0.47, P < .05). Whole-lung quantitative CT percentage emphysema also showed statistically significant correlation with percentage predicted Tlco (r = -0.65, P < .005). Conclusion Lobar ventilation and ADC values obtained from hyperpolarized 129Xe MR imaging demonstrated correlation with quantitative CT percentage emphysema on a lobar basis and with PFT results on a whole-lung basis. © RSNA, 2016.

Longitudinal changes in structural abnormalities using MDCT in COPD: do the CT measurements of airway wall thickness and small pulmonary vessels change in parallel with emphysematous progression?

Background

Recent advances in multidetector computed tomography (MDCT) facilitate acquiring important clinical information for managing patients with COPD. MDCT can detect the loss of lung tissue associated with emphysema as a low-attenuation area (LAA) and the thickness of airways as the wall area percentage (WA%). The percentage of small pulmonary vessels <5 mm² (% cross-sectional area [CSA] <5) has been recently recognized as a parameter for expressing pulmonary perfusion. We aimed to analyze the longitudinal changes in structural abnormalities using these CT parameters and analyze the effect of exacerbation and smoking cessation on structural changes in COPD patients.

Methods

We performed pulmonary function tests (PFTs), an MDCT, and a COPD assessment test (CAT) in 58 patients with COPD at the time of their enrollment at the hospital and 2 years later. We analyzed the change in clinical parameters including CT indices and examined the effect of exacerbations and smoking cessation on the structural changes.

Results

The CAT score and forced expiratory volume in 1 second (FEV₁) did not significantly change during the follow-up period. The parameters of emphysematous changes significantly increased. On the other hand, the WA% at the distal airways significantly decreased or tended to decrease, and the %CSA <5 slightly but significantly increased over the same period, especially in ex-smokers. The parameters of emphysematous change were greater in patients with exacerbations and continued to progress even after smoking cessation. In contrast, the WA% and %CSA <5 did not change in proportion to emphysema progression.

Conclusion

The WA% at the distal bronchi and the %CSA <5 did not change in parallel with parameters of LAA over the same period. We propose that airway disease and vascular remodeling may be reversible to some extent by smoking cessation and appropriate treatment. Optimal management may have a greater effect on pulmonary vascularity and airway disease than parenchymal deconstruction in the early stage of COPD.

Respiratory impedance is correlated with morphological changes in the lungs on three-dimensional CT in patients with COPD

The forced oscillation technique provides information concerning respiratory impedance, which comprises resistance and reactance of the respiratory system. However, its relationship with morphological changes of the lungs in chronic obstructive pulmonary disease (COPD) remains unclear. Respiratory impedance and spirometric data were evaluated in 98 patients with COPD and 49 reference subjects. Wall thickness (WT) and airway intraluminal area (Ai) of third- to sixth-generation bronchi, and percentage low-attenuation area with less than -950 HU (%LAA) of lungs were measured using three-dimensional computed tomography. COPD patients had higher respiratory impedance, decreased Ai, and increased %LAA compared with reference subjects. Indices of respiratory resistance and reactance and forced expiratory volume in 1 second (FEV₁) were correlated with Ai, and the association between percent predicted FEV₁ and Ai was predominant in distal bronchi. The difference in respiratory resistance between 5 Hz and 20 Hz (R5-R20) and FEV₁/forced vital capacity ratio (FEV₁/FVC) were correlated with WT. The %LAA was correlated with the FEV₁/FVC ratio and respiratory reactance. Airway function measurements with the forced oscillation technique provide complementary information to spirometry in COPD.