ORIGINAL RESEARCH article

Front. Med., 23 September 2021

Sec. Pulmonary Medicine

Volume 8 - 2021 | https://doi.org/10.3389/fmed.2021.705184

Pulmonary Emphysema Regional Distribution and Extent Assessed by Chest Computed Tomography Is Associated With Pulmonary Function Impairment in Patients With COPD

  • 1. Faculty of Medicine, Universidade do Porto, Porto, Portugal

  • 2. Serviço de Pneumologia, Centro Hospitalar de São João EPE, Porto, Portugal

  • 3. Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal

  • 4. Instituto de Biologia Molecular e Celular, Universidade do Porto, Porto, Portugal

  • 5. Serviço de Radiologia, Centro Hospitalar de São João EPE, Porto, Portugal

  • 6. Centro Hospitalar Vila Nova de Gaia/Espinho, Porto, Portugal

  • 7. Laboratory of Pulmonary Engineering, Biomedical Engineering Program, Alberto Luiz Coimbra Institute of Post-Graduation and Research in Engineering, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

  • 8. Department of Radiology, Universidade Federal Do Rio de Janeiro, Rio de Janeiro, Brazil

  • 9. IDOR–D'Or Institute for Research and Education, Rio de Janeiro, Brazil

  • 10. Laboratory of Respiration Physiology, Carlos Chagas Filho Institute of Biophysics, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

  • 11. Cardiovascular R&D Center, Department of Surgery and Physiology, Faculty of Medicine of the University of Porto, Porto, Portugal

Abstract

Objective: This study aimed to evaluate how emphysema extent and its regional distribution quantified by chest CT are associated with clinical and functional severity in patients with chronic obstructive pulmonary disease (COPD).

Methods/Design: Patients with a post-bronchodilator forced expiratory volume in one second (FEV1)/forced vital capacity (FVC) < 0.70, without any other obstructive airway disease, who presented radiological evidence of emphysema on visual CT inspection were retrospectively enrolled. A Quantitative Lung Imaging (QUALI) system automatically quantified the volume of pulmonary emphysema and adjusted this volume to the measured (EmphCTLV) or predicted total lung volume (TLV) (EmphPLV) and assessed its regional distribution based on an artificial neural network (ANN) trained for this purpose. Additionally, the percentage of lung volume occupied by low-attenuation areas (LAA) was computed by dividing the total volume of regions with attenuation lower or equal to −950 Hounsfield units (HU) by the predicted [LAA (%PLV)] or measured CT lung volume [LAA (%CTLV)]. The LAA was then compared with the QUALI emphysema estimations. The association between emphysema extension and its regional distribution with pulmonary function impairment was then assessed.

Results: In this study, 86 patients fulfilled the inclusion criteria. Both EmphCTLV and EmphPLV were significantly lower than the LAA indices independently of emphysema severity. CT-derived TLV significantly increased with emphysema severity (from 6,143 ± 1,295 up to 7,659 ± 1,264 ml from mild to very severe emphysema, p < 0.005) and thus, both EmphCTLV and LAA significantly underestimated emphysema extent when compared with those values adjusted to the predicted lung volume. All CT-derived emphysema indices presented moderate to strong correlations with residual volume (RV) (with correlations ranging from 0.61 to 0.66), total lung capacity (TLC) (from 0.51 to 0.59), and FEV1 (~0.6) and diffusing capacity for carbon monoxide DLCO (~0.6). The values of FEV1 and DLCO were significantly lower, and RV (p < 0.001) and TLC (p < 0.001) were significantly higher with the increasing emphysema extent and when emphysematous areas homogeneously affected the lungs.

Conclusions: Emphysema volume must be referred to the predicted and not to the measured lung volume when assessing the CT-derived emphysema extension. Pulmonary function impairment was greater in patients with higher emphysema volumes and with a more homogeneous emphysema distribution. Further studies are still necessary to assess the significance of CTpLV in the clinical and research fields.

Introduction

Chronic obstructive pulmonary disease (COPD) presents important morbidity and mortality worldwide with a prevalence as high as 14.2% in some areas of Portugal (13). COPD is characterized by persistent respiratory symptoms and airflow limitation caused by the pathophysiological changes ranging from chronic inflammation and narrowing of peripheral airways to emphysema (47).

The patients with predominant emphysema phenotype exhibit a more severe lung function impairment, such as decreased gas exchange efficiency, lower forced expiratory volume in one second (FEV1), and FEV1/forced vital capacity (FVC) ratio (810).

The correlation between the CT-derived emphysema indices and lung function impairment has been well-established (11, 12). However, only recent advances allowed the quantification of emphysema distribution by CT and its association with different degrees of clinical severity (4, 1315). None of these reports considered that hyperinflation is a major characteristic of emphysematous patients, which would increase total lung volume (TLV), possibly skewing the measurement of emphysema extent (16). The adjustment of emphysema volume to the predicted TLV might be of potential interest when one wishes to assess disease progression, since the increase in TLV may lead to a progressive underestimation of the extent of emphysema (13, 17, 18).

In the present study, we employed a QUAntitative Lung Imaging (QUALI) system as a computer-aided diagnosis (CAD) system that applies a previously proposed artificial neural network (ANN) to automatically classify and locate the lung patterns commonly used in clinical practice: emphysema, normal parenchyma, ground-glass opacity (GGO), crazy paving and linear opacities (CP/LO), and consolidation from the chest CT images (19). Additionally, low attenuation areas (LAA) were also calculated as the percentage of lung volume occupied by voxels with attenuation of −950 Hounsfield Units (HU) or less (6, 20, 21).

Thereafter, we aimed to evaluate whether the extension of emphysema quantified by several indices derived from CT-scan images and its regional distribution are associated with the clinical and functional features in emphysematous patients.

Methods

Study Design and Study Subjects

Patients from the Department of Pulmonology of Centro Hospitalar e Universitário de São João, Porto, Portugal, with a post-bronchodilator FEV1/FVC <0.70 were retrospectively selected. The patients with <18 years old and a history of asthma, bronchiectasis, tuberculosis sequelae, cystic fibrosis, lung fibrosis, asthma/COPD overlap syndrome, thoracic surgery, or other confounding diseases were excluded (Figure 1).

Figure 1

Chest CT Acquisition and Image Processing

The CT scans were done with either a 32-channel SOMATOM go. Up or a 16-channel SOMATOM go. Now (Siemens Healthineers, Germany). The acquisitions were gathered in the supine position with 120 kV and 120–300 mA, slice thickness ranging from 3 to 5 mm without superposition, and a 512 × 512, 768 × 768, or 1,024 × 1,024 voxels matrix.

First, we visually inspected the thoracic CT scans of the selected patients in images reconstructed with a soft tissue (standard) convolutional kernel. Only patients showing evidence of emphysema, and tested for pulmonary function [spirometry, body plethysmography and diffusing capacity for carbon monoxide (DLCO)] within a year prior to the CT-scan were included.

Thereafter, the CT images series reconstructed with a sharp (bone) convolution kernel were selected and lung parenchyma and airways were automatically segmented by applying the Region Growing algorithm using the module Chest Imaging Platform (CIP) and segmentation tool with the Generate Partial Lung Map Label (3D Slicer version 4.8.1 software) (22). The right and left lungs were vertically divided into three equal segments (Figure 2A). All the images were visually inspected, manually edited, and exported to an in-house developed software (QUALI) written in MATLAB® (MathWorks®, Natick, MA, USA) (19).

Figure 2

After lung segmentation, CT-derived TLV was computed as:

The predicted TLV for a healthy person lying supine in deep inspiration was computed as previously reported:

for males 800.6 ml should be added to PLV (23).

Development of the Supervised Neural Network Architecture

Two chest radiologists blinded to patient identification, clinical data, and outcomes independently selected regions of interests (ROIs) from the 228 chest CT-scan images of patients with SARS-CoV-2 and 30 patients with pulmonary emphysema to visually classify the ROIs as emphysema regions (Emph, n = 14,162), normal parenchyma (n = 31,785), ground-glass opacities (GGO, n = 2,369), crazy paving/linear opacities (CP/LO, n = 7,337), and consolidation (n = 16,053), totaling 71,706 ROIs.

This ROI dataset comprises images acquired with several CT equipment, such as a 64-channel multislice (Brilliance 40 scanner, Philips Medical Systems, Cleveland, OH, USA and General Electrics Lightspeed VCT, Chicago, IL, USA), a 128-channel multislice dual-source CT system (Somatom Definition Flash, Siemens, Forchheim, Germany), and a 16-channel multislice (Emotion 16 CT, Siemens, Erlangen, Germany). Moreover, the images were reconstructed with several convolution kernels (B50f, B60f, B70s, C, FC13, FC86, I50f2, L, LUNG, and SOFT), depending on the manufacturer of CT.

The ROIs were then balanced by the lowest number of ROIs (2,369 ROIs). Thus, the percentiles 2.5, 25, 50, 75, and 97.5% derived from the density histograms of 11,845 ROIs were used for ANN training (8,291 ROIs), validation (1,777 ROIs), and test (1,777 ROIs). Each ROI consisted of a circle with a fixed radius of 4 mm with a spanning area of about 30 voxels in each CT section.

The whole lung CT scans were processed by the proposed QUALI system (19), resulting in a 3D scan renderization that includes all the classes (Figure 2B) and a summary glyph (Figure 2C) to the underlying 3D scan data.

Determination of Emphysema Extent

From the ANN-based classification procedure, the extent of emphysema was assessed as the cumulative sum of the QUALI-referred emphysema areas adjusted to CTLV (EmphyCTLV) or PLV (EmphPLV).

Additionally, the low attenuation areas (LAA) were calculated accounting for the volume occupied by voxels with attenuation of −950 HU or less (21). LAA index was then adjusted to CTLV [LAA (%CTLV)] or PLV [LAA (%PLV)].

For the assessment of the regional distribution of emphysema, we verified that at least 5% of each third of the lung was classified as emphysematous regions (24).

Clinical and Pulmonary Function Assessment

The medical records were reviewed retrospectively regarding demographic and anthropometric data (sex, age, height, weight, and body mass index [BMI]), smoking status (current smoker, former smoker, or non-smoker), inhaled corticosteroid therapy, dyspnea severity as determined by the modified Medical Research Council (mMRC) scale (25), and number of exacerbations in the year before enrolment (26). The patients with pulmonary function tests (PFTs) taken within a year prior to enrollment and at least 6 months without any clinical exacerbation were accepted and their air flow, lung volumes, and DLCO were recorded, in accordance with the international guidelines (27, 28).

Statistical Analysis

The normality of the data (Kolmogorov–Smirnov test with Lilliefors' correction) and the homogeneity of variances (Levene's median test) were tested. Since all the variables had normal distribution and homogeneous variances, they were expressed as means and SD; categorical variables as absolute values and proportion.

A paired t-test was used to compare CT-derived emphysema areas adjusted or not to PLV. The Bland–Altman plot was used to evaluate the agreement among EmphPLV, EmphCTLV, LAA(%PLV), and LAA(%CTLV).

The correlations among EmphPLV, EmphCTLV, LAA(%PLV), and LAA(%CTLV) with pulmonary functional variables were assessed using the Spearman's correlation analysis (very weak, ρ = 0.00–0.19; weak, ρ = 0.20–0.39; moderate, ρ = 0.40–0.59; strong, ρ = 0.60–0.79; very strong, ρ ≥ 0.80).

One-way ANOVA was used to compare the variables among the Global Initiative for Chronic Obstructive Lung Disease (GOLD) severity subgroups. Bonferroni post-hoc analysis was applied for multiple comparisons among the subgroups. All the statistical analyses were performed with the SPSS Statistics software package, version 26.0 (IBM Corporation, Armonk, NY, USA).

Data Management/Ethical Considerations

All data were fairly and lawfully collected by the authors for the specified, explicit, and legitimate purposes and not further processed. The anonymity and confidentiality of collected information and its protection were guaranteed according to the EU Directive 95/46/CE. All data were initially stored on an external disk and later in the cloud with restricted access by user password.

Results

In this study, 389 subjects had FEV1/FVC <0.7 being 86 COPD patients with pulmonary emphysema (Figure 1). Their mean age was 67 ± 11 (SD) years, most of them were men (88.4%), all but two of the patients were current or former smokers. Less than half of the patients were under inhaled corticosteroid therapy (47.7%). Only 26.7% of them had frequent exacerbations (≥2 exacerbations in the last year). Severe symptoms (mMRC dyspnea scale score ≥2) were found in 20.9% of them.

The characteristics of the patients (Table 1) revealed a wide range of airflow limitation according to the GOLD severity classification (29), with a slight predominance of severe airflow limitation: mild, in 14 (16.3%); moderate, in 22 (25.6%); severe, in 36 (41.9%); and very severe, in 14 (16.3%).

Table 1

GOLD severity classification
Clinical featuresAll patients
(n = 86)
Mild
(n = 14)
Moderate
(n = 22)
Severe
(n = 36)
Very severe
(n = 14)
P-value
Clinical variables
Age (years)67 ± 1159.7 ± 13.369.8 ± 9.170.1 ± 10.361.3 ± 8.80.002
Male (gender)76 (88.4%)12 (85.7%)17 (77.3%)34 (97.1%)13 (92.9%)
BMI (kg/m2)23.6 ± 4.123.3 ± 324.7 ± 3.823.8 ± 3.821.1 ± 5.10.07
Smoking status
   Previous smokers53 (61.6%)6 (42.9%)17 (77.3%)22 (62.9%)7 (50%)
   Current smokers31 (36%)8 (57.1%)4 (18.2%)12 (34.3%)7 (50%)
Exacerbations ≥ 2 previous years23 (26.7%)3 (21.4%)3 (13.6%)10 (28.6%)6 (42.9%%)
Inhaled corticoid therapy41 (47.7%)2 (14.3%)13 (59.1%)23 (67.5%)3 (21.4%)
mMRC dyspnoea score ≥ 218 (20.9%)2 (14.3%)3 (13.6%)9 (25.7%)3 (21.4%)
Pulmonary function test variables
RV, % of predicted179.8 ± 59124.3 ± 32.5156.6 ± 48.2192.2 ± 40.7268.2 ± 34 <0.001
TLC, % of predicted120.3 ± 23108.4 ± 15.9114.7 ± 24.5124.4 ± 23.4136.3 ± 16.30.011
RV/TLC (%)61.6 ± 2246.4 ± 24.261.7 ± 27.965.4 ± 17.271.6 ± 7.40.021
FEV1, % of predicted51.3 ± 2393.1 ± 27.961.4 ± 10.339.4 ± 5.523.6 ± 3.7 <0.001
DLCO, % of predicted50 ± 22.671.3 ± 15.362.7 ± 15.445.2 ± 17.116.9 ± 7.8 <0.001
DLCO/VA, % of predicted61.6 ± 2475.1 ± 14.873.1 ± 14.660.1 ± 22.928 ± 13.2 <0.001
CT-derived variables
TLV (mL)6680 ± 13986143 ± 12956200 ± 13916856 ± 12877659 ± 12640.005
Predicted TLV (mL)5101 ± 5065275 ± 5114976 ± 6185126 ± 3975146 ± 4580.350
LAA (%PLV)17.2 ± 18.45.3 ± 6.88.5 ± 7.418.8 ± 18.738.4 ± 19.8 <0.001
LAA (%CTLV)12 ± 11.64.4 ± 5.46.4 ± 4.913 ± 11.425.4 ± 12.8 <0.001
QUALI-derived variables
EmphPLV (%)15.7 ± 18*4.6 ± 6.4*6.8 ± 7.0*17.3 ± 18.3*36.5 ± 20.2* <0.001
EmphCTLV (%)10.9 ± 113.9 ± 55.1 ± 4.611.9 ± 11.224.2 ± 13.1 <0.001
Emphysema (mL)804 ± 939249 ± 347345 ± 366871 ± 9101918 ± 1104 <0.001
Well-Aerated (mL)4839 ± 12484478 ± 16954923 ± 11684991 ± 11944801 ± 9680.623
GGO (mL)487 ± 346817 ± 573398 ± 233449 ± 280390 ± 1300.001
CP/LO (mL)394 ± 97433 ± 114383 ± 86391 ± 100383 ± 930.436
Consolidation (mL)156 ± 62165 ± 49150 ± 73153 ± 66166 ± 500.830

Demographic, clinical, and imaging characteristics of enrolled patients.

Values are presented as mean ± SD or n (%). GOLD, the Global Initiative for Chronic Obstructive Lung Disease; mMRC, modified Medical Research Council; RV, Residual volume; TLC, Total lung capacity; FEV1, Forced expiratory volume measured during the first second; DLCO, Diffusing capacity for carbon monoxide; VA, Alveolar volume; LAA (%PLV), Low attenuation areas expressed as a percentage of the predicted lung volume; LAA (%CTLV), Low attenuation areas expressed as a percentage of CT measured lung volume; EmphPLV (%), QUALI-derived emphysema volume expressed as a percentage of the predicted lung volume; EmphCTLV (%), QUALI-derived emphysema volume expressed as a percentage of CT measured lung volume. GGO, Ground-glass opacities; CP/LO, Crazy paving/linear opacities.

*

significantly different from LAA (%PLV);

*significantly different from LAA(%CTLV). Bold numbers highlight statistical signficance.

An overall tendency for hyperinflation was clear, with a mean (± SD) RV of 179.8 ± 59% of the predicted value, TLC of 120.3 ± 23% of predicted and a RV/TLC ratio of 61.6 ± 22%. The patients also showed a reduced FEV1 of 51.3 ± 23% and DLCO of 50 ± 22.6% of predicted. A significant impairment in pulmonary function was observed from mild to very severe GOLD emphysema severity classification (Table 1).

The overall agreement between QUALI and the classifications of radiologists was 95%, being 99.4% for emphysema, 96.3% for well-aerated regions, 91.1% for GGO, 88.9% for CP/LO, and 98.7% for consolidation in the test set (Figure 3A, second row and first column). A receiver operating characteristic (ROC) curve from each radiological pattern is presented in Figure 3B. The ANN classifier performance was higher than 0.99 (0.97–1.0, 95% CI) in all classes. The best validation performance occurred at epoch 75 with a minimal entropy of 0.032 (Figure 3C).

Figure 3

The CT density histogram of all ROIs visually assigned as emphysema, well-aerated regions, GGO, CP/LO, and consolidation is depicted in Figure 4. Note that there are several overlapping bins between the densities assigned to normal parenchyma and emphysema, as well as between normal parenchyma and GGO. Additionally, there is an important overlap among GGO, CP/LO, and consolidation, as previously reported (19).

Figure 4

Even though EmphCTLV and EmphPLV were strongly correlated (ρ2 = 0.98; ρ = 0.98 and p < 0.001), a significant bias of −4.8% (±2 SD ranging from −19 to 8.9%, p < 0.001) was observed especially in the subjects with higher emphysema content (EmphCTLV > 20%) (Figures 5A,B). In fact, in all but mild GOLD severity classification, EmphPLV was always significantly higher than EmphCTLV (Figure 5C).

Figure 5

A very strong correlation was observed between LAA(%PLV) and EmphPLV (Figure 6A) as well as between LAA(%CTLV) and EmphCTLV (Figure 6B). Although the LAA estimations were always significantly higher than EmphCTLV or EmphPLV (Table 1), the bias ranged ~1.1–1.5% of the measured or predicted lung volume (Figure 6).

Figure 6

The Spearman's correlation coefficients between EmphPLV, EmphCTLV, LAA (%PLV), or LAA (%CTLV) and RV, TLC, FEV1, DLCO, DLCO/VA, and FVC are provided in Table 2. In general, all the CT-derived indices moderately to strongly correlated with RV, TLC, FEV1, and DLCO and very weakly with FVC (Table 2).

Table 2

Pulmonary function test variables
RVTLCFEV1DLCODLCO/VAFVCRV/TLV
EmphPLV0.6560.586−0.603−0.581−0.594−0.1390.527
EmphCTLV0.6110.517−0.602−0.603−0.606−0.1750.518
LAA (%PLV)0.6560.580−0.598−0.594−0.610−0.1430.520
LAA (%)0.6020.508−0.589−0.617−0.608−0.1750.514

Spearman's correlation coefficients between QUALI-derived emphysema indices (EmphPLV and EmphCTLV), LAA index obtained directed from CT densitometry and pulmonary function test variables.

RV, Residual volume; TLC, Total lung capacity; FEV1, Forced expiratory volume measured during the first second; DLCO, Diffusing capacity for carbon monoxide; VA, Alveolar volume; LAA (%PLV), Low attenuation areas expressed as a percentage of the predicted lung volume; LAA (%CTLV), Low attenuation areas expressed as a percentage of CT measured lung volume; EmphPLV(%), QUALI-derived emphysema volume expressed as a percentage of the predicted lung volume; EmphCTLV (%), QUALI-derived emphysema volume expressed as a percentage of CT measured lung volume; Orange squares marked strong (r = 0.60–0.79), yellow moderate (r = 0.40–0.59), and gray very weak (r = 0.00–0.19) Spearman's correlations.

After dividing the left and right lungs in 3 thirds, EmphPLV was assessed for each segmented lung region. Thus, when emphysema was present in more than 4 lung thirds, the patients presented a significantly higher RV, TLC, and RV/TLC ratio and a significantly lower FEV1, DLCO, DLCO/VA than the patients with 1–4 emphysematous segmented regions of the lungs (Figure 7).

Figure 7

Discussion

Recently, we reported the performance of a software (QUALI) for automated classification and quantification of coronavirus disease-2019 (COVID-19) pneumonia extent from chest CT (19). In the present study, we included emphysema as another radiological pattern in QUALI and aimed to quantify emphysema extent and its regional distribution in the lungs. The updated QUALI system presented an excellent performance in classifying emphysema, with 99.4% agreement with consensual visually assigned emphysema areas by the two independent radiologists (Figure 3).

Even though the emphysematous changes recognized on density-based CT quantified emphysema (21, 30) have been proven to present strong associations with airflow obstruction, worsening the GOLD severity and air trapping (31, 32), the extent of pulmonary emphysema is usually calculated by considering the volume of areas identified as being emphysematous divided by the TLV measured by CT.

However, hyperinflation is a major characteristic of COPD patients with emphysema and correlates with disease severity (16). Thus, worsening of emphysema might imply an increased TLV leading to a progressive underestimation of emphysema extent by CT-derived indices. This fact could be particularly important in the follow-up studies to assess the rate of emphysema progression, and several approaches have been advanced to overcome such limitation (18, 33, 34).

Hence, we proposed to express all CT-derived emphysema indices as a percentage of the predicted TLV using a mathematical model that estimates the TLV of healthy subjects under identical conditions to those existing during the acquisition of chest-CT images (23).

In fact, all CT-derived emphysema estimates calculated both from an ANN-based classifier (EmphyCT) or from densitometric CT-derived analysis (LAA) consistently tend to underestimate the amount of emphysema when compared with the same indices adjusted to the predicted lung volume (Table 1 and Figure 5). Although the correlations between the emphysema indices adjusted by the predicted lung volume and pulmonary function did not differ from those adjusted by measured lung volume, the present data suggest that importance of recording the height of patient in the DICOM header file as a routine procedure in radiology department to calculate the predicted lung volume and avoid the underestimation of emphysema severity.

In addition, we should be aware that PFTs are conventionally performed with seated individuals, volumes are higher in the sitting position (35) and that mild and locally restricted emphysema can be identified on CT, but these radiologic changes may not affect the pulmonary function (5, 36).

Thus, we can probably highlight that the CT-derived emphysema indices should be adjusted by the predicted lung volume to better estimate emphysema extent especially in the patients with more severe GOLD emphysema severity classification. Future studies must also evaluate the influence of such a strategy on the evaluation of the progression of emphysema. Hence, we emphasize the importance of collecting the height of the patient and inserting it in the DICOM header file as a routine procedure in the radiology departments.

A very strong correlation between LAA and QUALI-derived emphysema indices (Figure 6) is quite expected as both use density characteristics in their computation. In LAA index, all voxels with densities equal to or below −950 HU are automatically assigned as emphysema areas. In fact, this strategy seems to overestimate the emphysema extent by the inclusion of low density but still, normal parenchyma voxels due to some density overlap between highly inflated normal areas and emphysema regions (Figure 4). Several studies also report the dependence of image acquisition and reconstruction parameters as potential limitations for lung densitometry analysis (17, 18, 33, 34, 37). This drawback might also affect the QUALI-derived emphysema indices. However, we emphasize that the same convolution sharp kernel filter was used in all the calculations, as recommended (37).

A potential advantage of the QUALI-derived emphysema indices compared to LAA densitometry method lies in its unique capability of identifying other parenchymal radiological patterns in addition to emphysema, since the overlap between COPD and other lung parenchymal diseases is common (30, 38, 39). Moreover, more complex computational methods using convolutional neural networks (30) that account not only to voxel density but also for voxel texture can be easily incorporated into QUALI soon.

Regarding emphysema distribution, our data showed that not just total emphysema extent but also a more homogenous emphysema distribution occurring in almost all lung regions seems to be present in patients with worse airflow obstruction, diffusing capacity, and hyperinflation (Figure 7). However, this finding must be cautiously interpreted because several studies report that the occurrence of emphysema in the middle and lower lung thirds seems to be the main contributor to lung function decline (4042). Conversely, it has also been reported that independently of the emphysema extent, heavy (former) smokers with upper lobe-predominant CT-quantified emphysema exhibit a faster decrease in lung function than those with lower lobe-predominant CT-quantified emphysema (4).

The complexity between the emphysema progression and the compensatory pulmonary mechanisms to optimize gas-exchange might explain the relevance of regional distribution of emphysematous areas and its impact on the pulmonary. In fact, both emphysema extent and the regional occurrence of emphysema take place in the face of disease progression, thus rendering the isolation of their individual contribution or even to determine which one prevails was quite difficult (13, 41).

The present study presents some limitations. First, the female gender is not well-represented, but this lower proportion reflects daily clinical practice. Second, the disparity between scanners, scanner models, and image reconstruction kernels limit the application of both LAA and QUALI derived indices application in the clinical scenario. As a mitigation procedure, we used the same convolution kernel (sharp) and with the same slice thickness in all the studied images as previously recommended (37). Third, at the beginning of our study, we excluded 28 patients because the segmentation tool was not capable of performing segmentation due to their very hyperinflated lungs. This is an undesirable limitation since it prevents the use of our ANN in patients that would putatively benefit from it to a greater extent. Additionally, the absence of these patients could cause an overestimation of the correlation between lung function and the extent of emphysema as a percentage of the CTLV. Additionally, the division of the right and left lung into thirds represents only an approximation of the actual anatomy of the lung lobes. Certainly, segmentation into lobes considering each respective fissure would be very laborious and would require more computational resources.

To overcome such problems, we are working to apply a more complex convolutional neural network that does not need previous segmentation of the CT scans and works independently of the image reconstruction kernels used and will be able to automatically segment all the pulmonary lobes.

Apart from these limitations, the CT-derived emphysema indices can be adjusted to the predicted lung volume, being a simple and fast estimate of emphysema extent that can be obtained with minimal user intervention. Moreover, since COPD patients with emphysema benefit from long-term follow-up (29, 43), these are considerable advantages over visual evaluations of CT scans that are time-consuming, encompass interobserver variability, and are unable to quantify emphysema (44). Certainly, CTpLV is a potential index for the proper assessment of emphysematous lungs, however, the significance of CTpLV in clinical and research fields should be better elucidated in future studies.

In summary, in this group of COPD patients with pulmonary emphysema, the automatic quantification of the CT-derived emphysema indices adjusted to the predicted TLV enabled the quantification of emphysema extension and its regional distribution. The latter presented a moderate to strong correlation with pulmonary function tests parameters, especially with those associated with air trapping and hyperinflation. Pulmonary function impairment seemed to be more important in the patients with higher emphysema volumes and with a more homogeneous emphysema distribution among upper, middle, and lower lungs.

Funding

This research was supported by the Brazilian Council for Scientific and Technological Development (Conselho Nacional de Desenvolvimento Científico e Tecnológico, Grants Nos. 302702/2017-2 and 302839/2017-8) and the Rio de Janeiro State Research Supporting Foundation (Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro, Grants Nos. E-26/211.867/2016, E-26/202.785/2017, E-26/203.001/2018), and by national funds through FCT, Cardiovascular R&D Center – UnIC (UIDB/00051/2020 and UIDP/00051/2020).

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The paper was submitted to the Centro Hospitalar de São João (CHSJ) Research Ethics Committee, who approved the present study (CES 56/21), and that it complied with the current national and international standards. As this is a retrospective study without intervention the consent form was waived by the respective ethics committees.

Author contributions

PG: capture and organization of clinical data, image processing and segmentation, analysis of results and statistical evaluation, and writing of the text and submission of the article. HB: capture and organization of clinical data, results discussion, and draft review. AC: capture of clinical data. AL: capture and organization of clinical data. AG: image processing and analysis of results and theoretical development of the neural network. RR and WZ: results discussion and draft review. ASC: image processing and analysis of results, theoretical development of the neural network and the computation method of voxel-to-voxel analysis, results discussion, and draft review. All authors approved the final version of the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    LozanoRNaghaviMForemanKLimSShibuyaKAboyansVet al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. (2013) 380:2095128. 10.1016/S0140-6736(12)61728-0

  • 2.

    VosTFlaxmanADNaghaviMLozanoRMichaudCEzzatiMet al. Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. (2013) 380:216396. 10.1016/S0140-6736(12)61729-2

  • 3.

    BárbaraCRodriguesFDiasHCardosoJAlmeidaJMatosMJet al. Prevalência da doença pulmonar obstrutiva crónica em Lisboa, Portugal: estudo Burden of Obstructive Lung Disease. Rev Port Pneumol. (2013) 19:96105. 10.1016/j.rppneu.2012.11.004

  • 4.

    HoeseinFAAMRikxoortE vanGinnekenB vanJongPA deProkopMLammersJ-WJet al. Computed tomography-quantified emphysema distribution is associated with lung function decline. Eur Respir J. (2012) 40:84450. 10.1183/09031936.00186311

  • 5.

    PrattPCKilburnKH. A modern concept of the emphysemas based on correlations of structure and function. Hum Pathol. (1970) 1:44363. 10.1016/S0046-8177(70)80077-6

  • 6.

    ThurlbeckWMMüllerNL. Emphysema: definition, imaging, and quantification. Am J Roentgenol. (1994) 163:101725. 10.2214/ajr.163.5.7976869

  • 7.

    SniderGLKleinermanJThurlbeckWMBengaliZH. The definition of emphysema. Am Rev Respir Dis. (1985) 132:1825.

  • 8.

    BoschettoPQuintavalleSZeniELeprottiSPotenaABallerinLet al. Association between markers of emphysema and more severe chronic obstructive pulmonary disease. Thorax. (2006) 61:1037. 10.1136/thx.2006.058321

  • 9.

    BoschettoPMiniatiMMiottoDBraccioniFRosaEDBononiIet al. Predominant emphysema phenotype in chronic obstructive pulmonary disease patients. Eur Respir J. (2003) 21:4504. 10.1183/09031936.03.00048703

  • 10.

    ElbehairyAFCiavagliaCEWebbKAGuenetteJAJensenDMouradSMet al. Pulmonary gas exchange abnormalities in mild chronic obstructive pulmonary disease. Implications for dyspnea and exercise intolerance. Am J Resp Crit Care. (2015) 191:138494. 10.1164/rccm.201501-0157OC

  • 11.

    NakanoYMuroSSakaiHHiraiTChinKTsukinoMet al. Computed tomographic measurements of airway dimensions and emphysema in smokers. Am J Resp Crit Care. (2000) 162:11028. 10.1164/ajrccm.162.3.9907120

  • 12.

    ParrDG. Validation of computed tomographic lung densitometry for monitoring emphysema in 1-antitrypsin deficiency. Thorax. (2006) 61:48590. 10.1136/thx.2005.054890

  • 13.

    ParrDGSevenoaksMDengCStoelBCStockleyRA. Detection of emphysema progression in alpha 1-antitrypsin deficiency using CT densitometry; methodological advances. Respir Res. (2008) 9:21. 10.1186/1465-9921-9-21

  • 14.

    MartiniKFrauenfelderT. Advances in imaging for lung emphysema. Ann Transl Med. (2020) 8:1467. 10.21037/atm.2020.04.44

  • 15.

    SørensenLNielsenMPetersenJPedersenJHDirksenABruijneM. Chronic obstructive pulmonary disease quantification using CT texture analysis and densitometry: results from the Danish lung cancer screening trial. Am J Roentgenol. (2020) 214:126979. 10.2214/AJR.19.22300

  • 16.

    DubéB-PGuerderAMorelot-PanziniCLavenezianaP. The clinical relevance of the emphysema-hyperinflated phenotype in COPD. Copd Res Pract. (2016) 2:1. 10.1186/s40749-015-0017-7

  • 17.

    StoelBCVroomanHStolkJReiberJ. Sources of error in lung densitometry with CT. Investig Radiol. (1999) 230:3039. 10.1148/radiol.2303031314

  • 18.

    StoelBCPutterHBakkerMEDirksenAStockleyRAPiitulainenEet al. Volume correction in computed tomography densitometry for follow-up studies on pulmonary emphysema. Proc Am Thoracic Soc. (2008) 5:91924. 10.1513/pats.200804-040QC

  • 19.

    CarvalhoARSGuimarãesAWerberichGMCastroSN dePintoJSFSchmittWRet al. COVID-19 chest computed tomography to stratify severity and disease extension by artificial neural network computer-aided diagnosis. Front Med. (2020) 7:577609. 10.3389/fmed.2020.577609

  • 20.

    WangZGuSLeaderJKKunduSTedrowJRSciurbaFCet al. Optimal threshold in CT quantification of emphysema. Eur Radiol. (2013) 23:97584. 10.1007/s00330-012-2683-z

  • 21.

    MüllerNLStaplesCAMillerRRAbboudRT. “Density mask” an objective method to quantitate emphysema using computed tomography. Chest. (1988) 94:7827. 10.1378/chest.94.4.782

  • 22.

    FedorovABeichelRKalpathy-CramerJFinetJFillion-RobinJ-CPujolSet al. 3D Slicer as an image computing platform for the Quantitative Imaging Network. Magnetic Resonance Imaging. (2012) 30:132341. 10.1016/j.mri.2012.05.001

  • 23.

    CarvalhoARSGuimarãesAGarciaT de SOWerberichGMCeottoVFBozzaFAet al. Estimating COVID-19 pneumonia extent and severity from chest computed tomography. Front Physiol. (2021) 12:617657. 10.3389/fphys.2021.617657

  • 24.

    SmithBMBarrRG. Establishing normal reference values in quantitative computed tomography of emphysema. J Thoracic Imaging. (2013) 28:2803. 10.1097/RTI.0b013e3182a0d805

  • 25.

    MahlerDAWellsCK. Evaluation of clinical methods for rating dyspnea. Chest. (1988) 93:5806. 10.1378/chest.93.3.580

  • 26.

    Rodriguez-RoisinR. Toward a consensus definition for COPD exacerbations. Chest. (2000) 117:398S401. 10.1378/chest.117.5_suppl_2.398S

  • 27.

    RocaJBurgosFBarberàJASunyerJRodriguez-RoisinRCastellsaguéJet al. Prediction equations for plethysmographic lung volumes. Respiratory Med. (1998) 92:45460. 10.1016/S0954-6111(98)90291-8

  • 28.

    RocaJBurgosFSunyerJSaezMChinnSAntóJMet al. References values for forced spirometry. Eur Respir J. (1998) 11:135462. 10.1183/09031936.98.11061354

  • 29.

    HalpinDMGCrinerGJPapiASinghDAnzuetoAMartinezFJet al. Global initiative for the diagnosis, management, and prevention of chronic obstructive lung disease. The 2020 GOLD science committee report on COVID-19 and chronic obstructive pulmonary disease. Am J Resp Crit Care. (2021) 203:2436. 10.1164/rccm.202009-3533SO

  • 30.

    Bermejo-PeláezDAshSYWashkoGREstéparRSJLedesma-CarbayoMJ. Classification of interstitial lung abnormality patterns with an ensemble of deep convolutional neural networks. Sci Rep-uk. (2020) 10:338. 10.1038/s41598-019-56989-5

  • 31.

    RambodMPorszaszJMakeBJCrapoJDCasaburiR. Six-minute walk distance predictors, including CT scan measures, in the COPDGene cohort. Chest. (2012) 141:86775. 10.1378/chest.11-0870

  • 32.

    SchroederJDMcKenzieASZachJAWilsonCGCurran-EverettDStinsonDSet al. Relationships between airflow obstruction and quantitative CT measurements of emphysema, air trapping, and airways in subjects with and without chronic obstructive pulmonary disease. Am J Roentgenol. (2013) 201:W46070. 10.2214/AJR.12.10102

  • 33.

    ParrDGStoelBCStolkJNightingalePGStockleyRA. Influence of calibration on densitometric studies of emphysema progression using computed tomography. Am J Respiratory Critical Care Med. (2004) 170:88390. 10.1164/rccm.200403-326OC

  • 34.

    StaringMBakkerMEStolkJShamoninDPReiberJHCStoelBC. Towards local progression estimation of pulmonary emphysema using CT. Med Phys. (2016) 41:02190514. 10.1118/1.4851535

  • 35.

    CressoniMGallazziEChiurazziCMarinoABrioniMMengaFet al. Limits of normality of quantitative thoracic CT analysis. Critical Care. (2013) 17:R93. 10.1186/cc12738

  • 36.

    MeiLanKHNabihahTVictorKSusanMGrahamBStephenIRet al. CT emphysema score (EMPHysema, Age, Smoking, SIZE - EMPHASIZE) identifies clinically significant COPD. in C37. COPD: WHAT IS NEW IN IMAGING? American Thoracic Society International Conference Abstracts. American Thoracic Society (2014).

  • 37.

    MascalchiMCamiciottoliGDiciottiS. Lung densitometry: why, how and when. J Thoracic Dis. (2017) 9:331945. 10.21037/jtd.2017.08.17

  • 38.

    RocaJVargasCCanoISelivanovVBarreiroEMaierDet al. Chronic Obstructive Pulmonary Disease heterogeneity: challenges for health risk assessment, stratification and management. J Transl Med. (2014) 12:S3. 10.1186/1479-5876-12-S2-S3

  • 39.

    NemotoMNeiYBartholmaiBYoshidaKMatsuiHNakashitaTet al. Automated computed tomography quantification of fibrosis predicts prognosis in combined pulmonary fibrosis and emphysema in a real-world setting: a single-centre, retrospective study. Respir Res. (2020) 21:275. 10.1186/s12931-020-01545-3

  • 40.

    HaraguchiMShimuraSHidaWShiratoK. Pulmonary function and regional distribution of emphysema as determined by high-resolution computed tomography. Respiration. (1998) 65:1259. 10.1159/000029243

  • 41.

    TanabeNMuroSTanakaSSatoSOgumaTKiyokawaHet al. Emphysema distribution and annual changes in pulmonary function in male patients with chronic obstructive pulmonary disease. Respir Res. (2012) 13:31. 10.1186/1465-9921-13-31

  • 42.

    BastosHN eNevesIRedondoMCunhaRPereiraJMMagalhãesAet al. Influence of emphysema distribution on pulmonary function parameters in COPD patients. J Bras Pneumol. (2015) 41:48995. 10.1590/s1806-37562015000000136

  • 43.

    TorresJP deBastarrikaGWisniveskyJPAlcaideABCampoASeijoLMet al. Assessing the relationship between lung cancer risk and emphysema detected on low-dose CT of the chest. Chest. (2007) 132:19328. 10.1378/chest.07-1490

  • 44.

    UppaluriRMitsaTSonkaMHoffmanEAMcLennanG. Quantification of pulmonary emphysema from lung computed tomography images. Am J Resp Crit Care. (1997) 156:24854. 10.1164/ajrccm.156.1.9606093

Summary

Keywords

chronic obstructive pulmonary disease, pulmonary emphysema, computed tomography, CT-estimated emphysema, quantitative chest CT analysis

Citation

Gomes P, Bastos HN, Carvalho A, Lobo A, Guimarães A, Rodrigues RS, Zin WA and Carvalho ARS (2021) Pulmonary Emphysema Regional Distribution and Extent Assessed by Chest Computed Tomography Is Associated With Pulmonary Function Impairment in Patients With COPD. Front. Med. 8:705184. doi: 10.3389/fmed.2021.705184

Received

04 May 2021

Accepted

24 August 2021

Published

23 September 2021

Volume

8 - 2021

Edited by

Naoya Tanabe, Kyoto University, Japan

Reviewed by

Chin Kook Rhee, The Catholic University of Korea, South Korea; Kaoruko Shimizu, Hokkaido University, Japan

Updates

Copyright

*Correspondence: Alysson Roncally S. Carvalho ;

This article was submitted to Pulmonary Medicine, a section of the journal Frontiers in Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics