CT Imaging-Based Low-Attenuation Super Clusters in Three Dimensions and the Progression of Emphysema

CT Chest Imaging Using Normalized Join-Count: Predicting Emphysema Progression in the CanCOLD Study

Background Pre-existing emphysema is recognized as an indicator of future worsening in patients with chronic obstructive pulmonary disease (COPD) when observed through CT imaging. However, it remains uncertain whether additional factors, such as the spatial compactness of CT emphysema, might also serve as predictors of disease progression. Purpose To evaluate the relationship between the compactness of CT emphysema voxels and emphysema progression. Materials and Methods This secondary analysis uses data from the prospective Canadian Cohort Obstructive Lung Disease (CanCOLD) study, examining CT images obtained in participants with and without COPD at baseline and a 3-year follow-up time point (November 2009 to November 2018). Measurements of forced expiratory volume in first second of expiration (FEV1) and diffusing capacity of lung for carbon monoxide (DLco) were collected. The normalized join-count (NJC) measurement from baseline CT images and lung density (LD) changes were analyzed. Emphysema progression was defined as an annualized LD change of less than half an SD below the mean of the participants without COPD with no smoking history. Multivariable linear and logistic regression models were used to assess the association between baseline CT NJC measurements and the annualized change in LD, FEV1, DLco, and emphysema progression versus nonprogression. Results A total of 524 participants (mean age, 66 years ± 10 [SD]; 293 male) (FEV1 percent predicted, 88% ± 19; FEV1/FVC, 67% ± 9; DLco percent predicted, 105% ± 25) were analyzed, 187 (36%) of whom had COPD. CT NJC was associated with the annualized change in LD (P < .001), FEV1 (P = .02), and DLco (P = .01). Additionally, CT NJC predicted emphysema progression versus nonprogression (odds ratio, 2.24; 95% CI: 1.37, 3.50; P < .001). Conclusion The spatial distribution, or "compactness," of CT emphysema voxels predicted emphysema progression in individuals with and without COPD. ClinicalTrials.gov Identifier: NCT00920348 © RSNA, 2024 Supplemental material is available for this article.

A scale-free model of acute and ventilator-induced lung injury: a network theory approach inspired by seismology

Introduction

Acute respiratory distress syndrome (ARDS) presents a significant clinical challenge, with ventilator-induced lung injury (VILI) being a critical complication arising from life-saving mechanical ventilation. Understanding the spatial and temporal dynamics of VILI can inform therapeutic strategies to mitigate lung damage and improve outcomes.

Methods

Histological sections from initially healthy mice and pulmonary lavage-injured mice subjected to a second hit of VILI were segmented with Ilastik to define regions of lung injury. A scale-free network approach was applied to assess the correlation between injury regions, with regions of injury represented as 'nodes' in the network and 'edges' quantifying the degree of correlation between nodes. A simulated time series analysis was conducted to emulate the temporal sequence of injury events.

Results

Automated segmentation identified different lung regions in good agreement with manual scoring, achieving a sensitivity of 78% and a specificity of 85% across 'injury' pixels. Overall accuracy across 'injury', 'air', and 'other' pixels was 81%. The size of injured regions followed a power-law distribution, suggesting a 'rich-get-richer' phenomenon in the distribution of lung injury. Network analysis revealed a scale-free distribution of injury correlations, highlighting hubs of injury that could serve as focal points for therapeutic intervention. Simulated time series analysis further supported the concept of secondary injury events following an initial insult, with patterns resembling those observed in seismological studies of aftershocks.

Conclusion

The size distribution of injured regions underscores the spatially heterogeneous nature of acute and ventilator-induced lung injury. The application of network theory demonstrates the emergence of injury 'hubs' that are consistent with a 'rich-get-richer' dynamic. Simulated time series analysis demonstrates that the progression of injury events in the lung could follow spatiotemporal patterns similar to the progression of aftershocks in seismology, providing new insights into the mechanisms of injury distribution and propagation. Both phenomena suggest a potential for interventions targeting these injury 'hubs' to reduce the impact of VILI in ARDS management.

Quantifying the spatial clustering characteristics of radiographic emphysema explains variability in pulmonary function

Quantitative assessment of emphysema in CT scans has mostly focused on calculating the percentage of lung tissue that is deemed abnormal based on a density thresholding strategy. However, this overall measure of disease burden discards virtually all the spatial information encoded in the scan that is implicitly utilized in a visual assessment. This simplification is likely grouping heterogenous disease patterns and is potentially obscuring clinical phenotypes and variable disease outcomes. To overcome this, several methods that attempt to quantify heterogeneity in emphysema distribution have been proposed. Here, we compare three of those: one based on estimating a power law for the size distribution of contiguous emphysema clusters, a second that looks at the number of emphysema-to-emphysema voxel adjacencies, and a third that applies a parametric spatial point process model to the emphysema voxel locations. This was done using data from 587 individuals from Phase 1 of COPDGene that had an inspiratory CT scan and plasma protein abundance measurements. The associations between these imaging metrics and visual assessment with clinical measures (FEV[Formula: see text], FEV[Formula: see text]-FVC ratio, etc.) and plasma protein biomarker levels were evaluated using a variety of regression models. Our results showed that a selection of spatial measures had the ability to discern heterogeneous patterns among CTs that had similar emphysema burdens. The most informative quantitative measure, average cluster size from the point process model, showed much stronger associations with nearly every clinical outcome examined than existing CT-derived emphysema metrics and visual assessment. Moreover, approximately 75% more plasma biomarkers were found to be associated with an emphysema heterogeneity phenotype when accounting for spatial clustering measures than when they were excluded.

Climate-mediated air pollution associated with COPD severity

Air pollution has been reported to be associated with chronic obstructive pulmonary disease (COPD). Our study aim was to examine the mediating effects of air pollution on climate-associated health outcomes of COPD patients. A cross-sectional study of 117 COPD patients was conducted in a hospital in Taiwan. We measured the lung function, 6-min walking distance, oxygen desaturation, white blood cell count, and percent emphysema (low attenuation area, LAA) and linked these to 0-1-, 0-3-, and 0-5-year lags of individual-level exposure to relative humidity (RH), temperature, and air pollution. Linear regression models were conducted to examine associations of temperature, RH, and air pollution with severity of health outcomes. A mediation analysis was conducted to examine the mediating effects of air pollution on the associations of RH and temperature with health outcomes. We observed that a 1 % increase in the RH was associated with increases in forced expiratory volume in 1 s (FEV₁), eosinophils, and lymphocytes, and a decrease in the total-lobe LAA. A 1 °C increase in temperature was associated with decreases in oxygen desaturation, and right-, left-, and upper-lobe LAA values. Also, a 1 μg/m³ increase in PM_2.5 was associated with a decrease in the FEV₁ and an increase in oxygen desaturation. A 1 μg/m³ increases in PM₁₀ and PM_2.5 was associated with increases in the total-, right-, left, upper-, and lower-lobe (PM_2.5 only) LAA. A one part per billion increase in NO₂ was associated with a decrease in the FEV₁ and an increase in the upper-lobe LAA. Next, we found that NO₂ fully mediated the association between RH and FEV₁. We found PM_2.5 fully mediated associations of temperature with oxygen saturation and total-, right-, left-, and upper-lobe LAA. In conclusion, climate-mediated air pollution increased the risk of decreasing FEV₁ and oxygen saturation and increasing emphysema severity among COPD patients. Climate change-related air pollution is an important public health issue, especially with regards to respiratory disease.

Spatiotemporal distribution of cellular injury and leukocytes during the progression of ventilator-induced lung injury

Supportive mechanical ventilation is a necessary lifesaving treatment for acute respiratory distress syndrome (ARDS). This intervention often leads to injury exacerbation by ventilator-induced lung injury (VILI). Patterns of injury in ARDS and VILI are recognized to be heterogeneous; however, quantification of these injury distributions remains incomplete. Developing a more detailed understanding of injury heterogeneity, particularly how it varies in space and time, can help elucidate the mechanisms of VILI pathogenesis. Ultimately, this knowledge can be used to develop protective ventilation strategies that slow disease progression. To expand existing knowledge of VILI heterogeneity, we document the spatial evolution of cellular injury distribution and leukocyte infiltration, on the micro- and macroscales, during protective and injurious mechanical ventilation. We ventilated naïve mice using either high inspiratory pressure and zero positive end-expiratory pressure ventilation or low tidal volume with positive end-expiratory pressure. Distributions of cellular injury, identified with propidium iodide staining, were microscopically analyzed at three levels of injury severity. Cellular injury initiated in diffuse, quasi-random patterns, and progressed through expansion of high-density regions of injured cells termed "injury clusters." The density profile of the expanding injury regions suggests that stress shielding occurs, protecting the already injured regions from further damage. Spatial distribution of leukocytes did not correlate with that of cellular injury or ventilation-induced changes in lung function. These results suggest that protective ventilation protocols should protect the interface between healthy and injured regions to stymie injury propagation.

Differentiating COPD and Asthma using Quantitative CT Imaging and Machine Learning

There are similarities and differences between chronic obstructive pulmonary disease (COPD) and asthma patients in terms of computed tomography (CT) disease-related features. Our objective was to determine the optimal subset of CT imaging features for differentiating COPD and asthma using machine learning.COPD and asthma patients were recruited from Heidelberg University Hospital. CT was acquired and 93 features were extracted (VIDA Diagnostics): percentage of low-attenuating-areas below -950HU (LAA₉₅₀), LAA₉₅₀ hole count, estimated airway-wall-thickness for a 10 mm internal perimeter airway (Pi10), total-airway-count (TAC), as well as inner/outer perimeter/areas and wall thickness for each of five segmental airways, and the average of those five airways. Hybrid feature selection was used to select the optimum number of features, and support vector machine was used to classify COPD and asthma.Ninety-five participants were included (n=48 COPD; n=47 asthma); there were no differences between COPD and asthma for age (p=0.25) or FEV₁ (p=0.31). In a model including all CT features, the accuracy and F1-score was 80% and 81%, respectively. The top features were: LAA_950, LAA₉₅₀ hole count, average outer and inner airway perimeter, outer and inner airway area RB1, and TAC. In the model with only airway features, the accuracy and F1-score were 66% and 68%, respectively. The top features were: inner area RB1, wall thickness RB1, outer area LB1, TAC LB10, average outer/inner perimeter, Pi10, and TAC.In conclusions, COPD and asthma can be differentiated using machine learning with moderate-high accuracy by a subset of only 7 CT features.

A Personalized Spring Network Representation of Emphysematous Lungs From CT Images

Emphysema is a progressive disease characterized by irreversible tissue destruction and airspace enlargement, which manifest as low attenuation area (LAA) on CT images. Previous studies have shown that inflammation, protease imbalance, extracellular matrix remodeling and mechanical forces collectively influence the progression of emphysema. Elastic spring network models incorporating force-based mechanical failure have been applied to investigate the pathogenesis and progression of emphysema. However, these models were general without considering the patient-specific information on lung structure available in CT images. The aim of this work was to develop a novel approach that provides an optimal spring network representation of emphysematous lungs based on the apparent density in CT images, allowing the construction of personalized networks. The proposed method takes into account the size and curvature of LAA clusters on the CT images that correspond to a pre-stressed condition of the lung as opposed to a naïve method that excludes the effects of pre-stress. The main findings of this study are that networks constructed by the new method 1) better preserve LAA cluster sizes and their distribution than the naïve method; and 2) predict different course of emphysema progression compared to the naïve method. We conclude that our new method has the potential to predict patient-specific emphysema progression which needs verification using clinical data.

A homological approach to a mathematical definition of pulmonary fibrosis and emphysema on computed tomography

Three-dimensional imaging is essential to evaluate local abnormalities and understand structure-function relationships in an organ. However, quantifiable and interpretable methods to localize abnormalities remain unestablished. Visual assessments are prone to bias, machine learning methods depend on training images, and the underlying decision principle is usually difficult to interpret. Here, we developed a homological approach to mathematically define emphysema and fibrosis in the lungs on computed tomography (CT). With the use of persistent homology, the density of homological features, including connected components, tunnels, and voids, was extracted from the volumetric CT scans of lung diseases. A pair of CT values at which each homological feature appeared (birth) and disappeared (death) was computed by sweeping the threshold levels from higher to lower CT values. Consequently, fibrosis and emphysema were defined as voxels with dense voids having a longer lifetime (birth-death difference) and voxels with dense connected components having a lower birth, respectively. In an independent dataset including subjects with idiopathic pulmonary fibrosis (IPF), chronic obstructive pulmonary disease (COPD), and combined pulmonary fibrosis and emphysema (CPFE), the proposed definition enabled accurate segmentation with comparable quality to deep learning in terms of Dice coefficients. Persistent homology-defined fibrosis was closely associated with physiological abnormalities such as impaired diffusion capacity and long-term mortality in subjects with IPF and CPFE, and persistent homology-defined emphysema was associated with impaired diffusion capacity in subjects with COPD. The present persistent homology-based evaluation of structural abnormalities could help explore the clinical and physiological impacts of structural changes and morphological mechanisms of disease progression.NEW & NOTEWORTHY This study proposes a homological approach to mathematically define a three-dimensional texture feature of emphysema and fibrosis on chest computed tomography using persistent homology. The proposed definition enabled accurate segmentation with comparable quality to deep learning while offering higher interpretability than deep learning-based methods.

Three-Dimensional CT for Quantification of Longitudinal Lung and Pneumonia Variations in COVID-19 Patients

Objectives: Visual chest CT is subjective with interobserver variability. We aimed to quantify the dynamic changes of lung and pneumonia on three-dimensional CT (3D-CT) images in coronavirus disease 2019 (COVID-19) patients during hospitalization.

Methods: A total of 110 laboratory-confirmed COVID-19 patients who underwent chest CT from January 3 to February 29, 2020 were retrospectively reviewed. Pneumonia lesions were classified as four stages: early, progressive, peak, and absorption stages on chest CT. A computer-aided diagnostic (CAD) system calculated the total lung volume (TLV), the percentage of low attenuation areas (LAA%), the volume of pneumonia, the volume of ground-glass opacities (GGO), the volume of consolidation plus the GGO/consolidation ratio. The CT score was visually assessed by radiologists. Comparisons of lung and pneumonia parameters among the four stages were performed by one-way ANOVA with post-hoc tests. The relationship between the CT score and the volume of pneumonia, and between LAA% and the volume of pneumonia in four stages was assessed by Spearman's rank correlation analysis.

Results: A total of 534 chest CT scans were performed with a median interval of 4 days. TLV, LAA%, and the GGO/consolidation ratio were significantly decreased, while the volume of pneumonia, GGO, and consolidation were significantly increased in the progressive and peak stages (for all, P < 0.05). The CT score was significantly correlated with the pneumonia volume in the four stages (r = 0.731, 0.761, 0.715, and 0.669, respectively, P < 0.001).

Conclusion: 3D-CT could be used as a useful quantification method in monitoring the dynamic changes of COVID-19 pneumonia.

Fractal Analysis of Lung Structure in Chronic Obstructive Pulmonary Disease

Chest CT is often used for localizing and quantitating pathologies associated with chronic obstructive pulmonary disease (COPD). While simple measurements of areas and volumes of emphysema and airway structure are common, these methods do not capture the structural complexity of the COPD lung. Since the concept of fractals has been successfully applied to evaluate complexity of the lung, this review is aimed at describing the fractal properties of airway disease, emphysema, and vascular abnormalities in COPD. An object forms a fractal if it exhibits the property of self-similarity at different length scales of evaluations. This fractal property is governed by power-law functions characterized by the fractal dimension (FD). Power-laws can also manifest in other statistical descriptors of structure such as the size distribution of emphysema clusters characterized by the power-law exponent D. Although D is not the same as FD of emphysematous clusters, it is a useful index to characterize the spatial pattern of disease progression and predict clinical outcomes in patients with COPD. The FD of the airway tree shape and the D of the size distribution of airway branches have been proposed indexes of structural assessment and clinical predictions. Simulations are also useful to understand the mechanism of disease progression. Therefore, the power-law and fractal analysis of the parenchyma and airways, especially when combined with computer simulations, could lead to a better understanding of the structural alterations during the progression of COPD and help identify subjects at a high risk of severe COPD.

Spontaneous Breathing Through Increased Airway Resistance Augments Elastase-Induced Pulmonary Emphysema

Introduction

Resistive breathing (RB), the pathophysiologic hallmark of chronic obstructive pulmonary disease (COPD), especially during exacerbations, is associated with significant inflammation and mechanical stress on the lung. Mechanical forces are implicated in the progression of emphysema that is a major pathologic feature of COPD. We hypothesized that resistive breathing exacerbates emphysema.

Methods

C57BL/6 mice were exposed to 0.75 units of pancreatic porcine elastase intratracheally to develop emphysema. Resistive breathing was applied by suturing a nylon band around the trachea to reduce surface area to half for the last 24 or 72 hours of a 21-day time period after elastase treatment in total. Following RB (24 or 72 hours), lung mechanics were measured and bronchoalveolar lavage (BAL) was performed. Emphysema was quantified by the mean linear intercept (Lm) and the destructive index (DI) in lung tissue sections.

Results

Following 21 days of intratracheal elastase exposure, Lm and DI increased in lung tissue sections [Lm (μm), control 39.09±0.76, elastase 62.05±2.19, p=0.003 and DI, ctr 30.95±2.75, elastase 73.12±1.75, p<0.001]. RB for 72 hours further increased Lm by 64% and DI by 19%, compared to elastase alone (p<0.001 and p=0.02, respectively). RB induced BAL neutrophilia in elastase-treated mice. Static compliance (C_st) increased in elastase-treated mice [C_st (mL/cmH2O), control 0.067±0.001, elastase 0.109±0.006, p<0.001], but superimposed RB decreased C_st, compared to elastase alone [C_st (mL/cmH2O), elastase+RB24h 0.090±0.004, p=0.006 to elastase, elastase+RB72h 0.090±0.005, p=0.006 to elastase].

Conclusion

Resistive breathing augments pulmonary inflammation and emphysema in an elastase-induced emphysema mouse model.

Quantitative CT analysis in patients with pulmonary emphysema: is lung function influenced by concomitant unspecific pulmonary fibrosis?

Purpose

Quantitative analysis of CT scans has proven to be a reproducible technique, which might help to understand the pathophysiology of chronic obstructive pulmonary disease (COPD) and combined pulmonary fibrosis and emphysema. The aim of this retrospective study was to find out if the lung function of patients with COPD with Global Initiative for Chronic Obstructive Lung Disease (GOLD) stages III or IV and pulmonary emphysema is measurably influenced by high attenuation areas as a correlate of concomitant unspecific fibrotic changes of lung parenchyma.

Patients and methods

Eighty-eight patients with COPD GOLD stage III or IV underwent CT and pulmonary function tests. Quantitative CT analysis was performed to determine low attenuation volume (LAV) and high attenuation volume (HAV), which are considered to be equivalents of fibrotic (HAV) and emphysematous (LAV) changes of lung parenchyma. Both parameters were determined for the whole lung, as well as peripheral and central lung areas only. Multivariate regression analysis was used to correlate HAV with different parameters of lung function.

Results

Unlike LAV, HAV did not show significant correlation with parameters of lung function. Even in patients with a relatively high HAV of more than 10%, in contrast to HAV (p=0.786) only LAV showed a significantly negative correlation with forced expiratory volume in 1 second (r=-0.309, R2=0.096, p=0.003). A severe decrease of DLCO% was associated with both larger HAV (p=0.045) and larger LAV (p=0.001). Residual volume and FVC were not influenced by LAV or HAV.

Conclusion

In patients with COPD GOLD stage III-IV, emphysematous changes of lung parenchyma seem to have such a strong influence on lung function, which is a possible effect of concomitant unspecific fibrosis is overwhelmed.