Search

Rationale Emphysema progression is heterogeneous. Predicting temporal changes in lung density and detecting rapid progressors may facilitate selection of individuals for targeted therapies. Objective To test whether computed tomography (CT) radiomics can be used to predict changes in lung density and detect rapid progressors. Methods We extracted radiomics features from inspiratory chest CT in 4,575 subjects with and without airflow obstruction at enrollment, who completed a follow-up visit at approximately 5 years. We quantified emphysema using adjusted lung density (ALD) and estimated emphysema progression as the annualized change in ALD (∆ALD/year) between visits. We categorized participants into rapid progressors (>1% ∆ALD/year) and stable disease (≤1% ∆ALD/year). A gradient boosting model was used (1) to predict ALD at 5-years and (2) to identify rapid progressors. Four models using demographics (base clinical model); CT density; radiomics; and combined features (clinical, radiomics, and CT density) were evaluated and tested. Results There were 1,773 (38.8%) rapid progressors. For predicting ALD at 5-years in the 20% held-out data, the base model explained 31% of the variance (adjusted R2 = 0.31) whereas R2 was 0.74 for the CT density model, 0.66 for the radiomics-only model, and 0.77 for the combined features model. For detecting rapid progressors, the base model (AUC = 0.57, 95%CI 0.53-0.61) was outperformed by the radiomics-only model (AUC = 0.73, 95%CI 0.69-0.76, ∆ =0.0003, p < 0.001) and the combined model (AUC = 0.74, 95%CI 0.71-0.77, ∆ = 0.0003, p < 0.001). Conclusions Parenchymal and airway radiomics features derived from inspiratory scans can be used to predict temporal changes in lung density and help identify rapid progressors.

Abstract: Bronchiolitis obliterans syndrome (BOS) after hematopoietic cell transplantation (HCT) is associated with substantial morbidity and mortality. Quantitative computed tomography (qCT) can help diagnose advanced BOS meeting National Institutes of Health (NIH) criteria (NIH-BOS) but has not been used to diagnose early, often asymptomatic BOS (early BOS), limiting the potential for early intervention and improved outcomes. Using pulmonary function tests (PFTs) to define NIH-BOS, early BOS, and mixed BOS (NIH-BOS with restrictive lung disease) in patients from 2 large cancer centers, we applied qCT to identify early BOS and distinguish between types of BOS. Patients with transient impairment or healthy lungs were included for comparison. PFTs were done at month 0, 6, and 12. Analysis was performed with association statistics, principal component analysis, conditional inference trees (CITs), and machine learning (ML) classifier models. Our cohort included 84 allogeneic HCT recipients, 66 with BOS (NIH-defined, early, or mixed) and 18 without BOS. All qCT metrics had moderate correlation with forced expiratory volume in 1 second, and each qCT metric differentiated BOS from those without BOS (non-BOS; P < .0001). CITs distinguished 94% of participants with BOS vs non-BOS, 85% of early BOS vs non-BOS, 92% of early BOS vs NIH-BOS. ML models diagnosed BOS with area under the curve (AUC) of 0.84 (95% confidence interval [CI], 0.74-0.94) and early BOS with AUC of 0.84 (95% CI, 0.69-0.97). qCT metrics can identify individuals with early BOS, paving the way for closer monitoring and earlier treatment in this vulnerable population., (© 2024 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

Background: Luminal narrowing is a hallmark feature of airway remodeling in COPD, but current measures focus on airway wall remodeling. Quantification of the natural increase in cumulative cross-sectional area along the length of the human airway tree can facilitate assessment of airway narrowing., Methods: We analysed the airway trees of 7641 subjects enrolled in the multicenter COPDGene cohort. Airway luminal tapering was assessed by estimating the slope of the change in cumulative cross-sectional area along the length of the airway tree over successive generations (T-Slope). We performed multivariable regression analyses to test the associations between T-Slope and lung function, St. George's Respiratory Questionnaire (SGRQ), modified Medical Research Council (mMRC) dyspnea score, 6-minute walk distance (6 MWD), FEV 1 change, exacerbations, and all-cause mortality after adjusting for demographics, %CT emphysema, and total airway count., Results: The T-Slope decreased with increasing COPD severity: 2.69 (0.70) in nonsmokers and 2.33 (0.70), 2.11 (0.65), 1.78 (0.58), 1.60 (0.53), and 1.57 (0.52) in GOLD stages 0 through 4 respectively (Jonckheere-Terpstra p=0.04). On multivariable analyses, the T-Slope was independently associated with FEV 1 (β=0.13 L, 95% CI 0.10 to 0.15, p<0.001), 6MWD (β=15.0 m, 95%CI 10.8 to 19.2, p<0.001), change in FEV 1 (β=-4.50 ml·year -1 , 95% CI -7.32 to -1.67; p=0.001), exacerbations (IRR=0.78, 95% CI 0.73 to 0.83, p<0.001), and mortality (HR=0.79, 95% CI 0.72 to 0.86, p<0.001)., Conclusion: T-Slope is a measure of airway luminal remodeling and is associated with respiratory morbidity and mortality., (Copyright ©The authors 2024. For reproduction rights and permissions contact permissions@ersnet.org.)

Rationale: Quantifying functional small airways disease (fSAD) requires additional expiratory computed tomography (CT) scan, limiting clinical applicability. Artificial intelligence (AI) could enable fSAD quantification from chest CT scan at total lung capacity (TLC) alone (fSAD TLC )., Objectives: To evaluate an AI model for estimating fSAD TLC and study its clinical associations in chronic obstructive pulmonary disease (COPD)., Methods: We analyzed 2513 participants from the SubPopulations and InteRmediate Outcome Measures in COPD Study (SPIROMICS). Using a subset ( n = 1055), we developed a generative model to produce virtual expiratory CTs for estimating fSAD TLC in the remaining 1458 SPIROMICS participants. We compared fSAD TLC with dual volume, parametric response mapping fSAD PRM . We investigated univariate and multivariable associations of fSAD TLC with FEV 1 , FEV 1 /FVC, six-minute walk distance (6MWD), St. George's Respiratory Questionnaire (SGRQ), and FEV 1 decline. The results were validated in a subset ( n = 458) from COPDGene study. Multivariable models were adjusted for age, race, sex, BMI, baseline FEV 1 , smoking pack years, smoking status, and percent emphysema., Measurements and Main Results: Inspiratory fSAD TLC was highly correlated with fSAD PRM in SPIROMICS (Pearson's R = 0.895) and COPDGene (R = 0.897) cohorts. In SPIROMICS, fSAD TLC was associated with FEV 1 (L) (adj.β = -0.034, P < 0.001), FEV 1 /FVC (adj.β = -0.008, P < 0.001), SGRQ (adj.β = 0.243, P < 0.001), and FEV 1 decline (mL / year) (adj.β = -1.156, P < 0.001). fSAD TLC was also associated with FEV 1 (L) (adj.β = -0.032, P < 0.001), FEV 1 /FVC (adj.β = -0.007, P < 0.001), SGRQ (adj.β = 0.190, P = 0.02), and FEV 1 decline (mL / year) (adj.β = -0.866, P = 0.001) in COPDGene. We found fSAD TLC to be more repeatable than fSAD PRM with intraclass correlation of 0.99 (95% CI: 0.98, 0.99) vs. 0.83 (95% CI: 0.76, 0.88)., Conclusions: Inspiratory fSAD TLC captures small airways disease as reliably as fSAD PRM and is associated with FEV 1 decline.

Mechanical ventilation exposes the lung to injurious stresses and strains that can negatively affect clinical outcomes in acute respiratory distress syndrome or cause pulmonary complications after general anesthesia. Excess global lung strain, estimated as increased respiratory system driving pressure, is associated with mortality related to mechanical ventilation. The role of small-dimension biomechanical factors underlying this association and their spatial heterogeneity within the lung are currently unknown. Using four-dimensional computed tomography with a voxel resolution of 2.4 cubic millimeters and a multiresolution convolutional neural network for whole-lung image segmentation, we dynamically measured voxel-wise lung inflation and tidal parenchymal strains. Healthy or injured ovine lungs were evaluated as the mechanical ventilation positive end-expiratory pressure (PEEP) was titrated from 20 to 2 centimeters of water. The PEEP of minimal driving pressure (PEEP DP ) optimized local lung biomechanics. We observed a greater rate of change in nonaerated lung mass with respect to PEEP below PEEP DP compared with PEEP values above this threshold. PEEP DP similarly characterized a breaking point in the relationships between PEEP and SD of local tidal parenchymal strain, the 95th percentile of local strains, and the magnitude of tidal overdistension. These findings advance the understanding of lung collapse, tidal overdistension, and strain heterogeneity as local triggers of ventilator-induced lung injury in large-animal lungs similar to those of humans and could inform the clinical management of mechanical ventilation to improve local lung biomechanics.

Purpose: To determine whether 4-dimensional computed tomography (4DCT) ventilation-based functional lung avoidance radiation therapy preserves pulmonary function compared with standard radiation therapy for non-small cell lung cancer (NSCLC)., Methods and Materials: This single center, randomized, phase 2 trial enrolled patients with NSCLC receiving curative intent radiation therapy with either stereotactic body radiation therapy or conventionally fractionated radiation therapy between 2016 and 2022. Patients were randomized 1:1 to standard of care radiation therapy or functional lung avoidance radiation therapy. The primary endpoint was the change in Jacobian-based ventilation as measured on 4DCT from baseline to 3 months postradiation. Secondary endpoints included changes in volume of high- and low-ventilating lung, pulmonary toxicity, and changes in pulmonary function tests (PFTs)., Results: A total of 122 patients were randomized and 116 were available for analysis. Median follow up was 29.9 months. Functional avoidance plans significantly (P < .05) reduced dose to high-functioning lung without compromising target coverage or organs at risk constraints. When analyzing all patients, there was no difference in the amount of lung showing a reduction in ventilation from baseline to 3 months between the 2 arms (1.91% vs 1.87%; P = .90). Overall grade ≥2 and grade ≥3 pulmonary toxicities for all patients were 24.1% and 8.6%, respectively. There was no significant difference in pulmonary toxicity or changes in PFTs between the 2 study arms. In the conventionally fractionated cohort, there was a lower rate of grade ≥2 pneumonitis (8.2% vs 32.3%; P = .049) and less of a decline in change in forced expiratory volume in 1 second (-3 vs -5; P = .042) and forced vital capacity (1.5 vs -6; P = .005) at 3 months, favoring the functional avoidance arm., Conclusions: There was no difference in posttreatment ventilation as measured by 4DCT between the arms. In the cohort of patients treated with conventionally fractionated radiation therapy with functional lung avoidance, there was reduced pulmonary toxicity, and less decline in PFTs suggesting a clinical benefit in patients with locally advanced NSCLC., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Chest computed tomography (CT) at inspiration is often complemented by an expiratory CT to identify peripheral airways disease. Additionally, co-registered inspiratory-expiratory volumes can be used to derive various markers of lung function. Expiratory CT scans, however, may not be acquired due to dose or scan time considerations or may be inadequate due to motion or insufficient exhale; leading to a missed opportunity to evaluate underlying small airways disease. Here, we propose LungViT- a generative adversarial learning approach using hierarchical vision transformers for translating inspiratory CT intensities to corresponding expiratory CT intensities. LungViT addresses several limitations of the traditional generative models including slicewise discontinuities, limited size of generated volumes, and their inability to model texture transfer at volumetric level. We propose a shifted-window hierarchical vision transformer architecture with squeeze-and-excitation decoder blocks for modeling dependencies between features. We also propose a multiview texture similarity distance metric for texture and style transfer in 3D. To incorporate global information into the training process and refine the output of our model, we use ensemble cascading. LungViT is able to generate large 3D volumes of size 320×320×320 . We train and validate our model using a diverse cohort of 1500 subjects with varying disease severity. To assess model generalizability beyond the development set biases, we evaluate our model on an out-of-distribution external validation set of 200 subjects. Clinical validation on internal and external testing sets shows that synthetic volumes could be reliably adopted for deriving clinical endpoints of chronic obstructive pulmonary disease.

Rational: Ground glass opacities (GGO) in the absence of interstitial lung disease are understudied., Objective: To assess the association of GGO with white blood cells (WBCs) and progression of quantified chest CT emphysema., Methods: We analyzed data of participants in the Subpopulations and Intermediate Outcome Measures In COPD Study (SPIROMICS). Chest radiologists and pulmonologists labeled regions of the lung as GGO and adaptive multiple feature method (AMFM) trained the computer to assign those labels to image voxels and quantify the volume of the lung with GGO (%GGO AMFM ). We used multivariable linear regression, zero-inflated negative binomial, and proportional hazards regression models to assess the association of %GGO AMFM with WBC, changes in %emphysema, and clinical outcomes., Measurements and Main Results: Among 2,714 participants, 1,680 had COPD and 1,034 had normal spirometry. Among COPD participants, based on the multivariable analysis, current smoking and chronic productive cough was associated with higher %GGO AMFM . Higher %GGO AMFM was cross-sectionally associated with higher WBCs and neutrophils levels. Higher %GGO AMFM per interquartile range at visit 1 (baseline) was associated with an increase in emphysema at one-year follow visit by 11.7% (Relative increase; 95%CI 7.5-16.1%;P<0.001). We found no association between %GGO AMFM and one-year FEV 1 decline but %GGO AMFM was associated with exacerbations and all-cause mortality during a median follow-up time of 1,544 days (Interquartile Interval=1,118-2,059). Among normal spirometry participants, we found similar results except that %GGO AMFM was associated with progression to COPD at one-year follow-up., Conclusions: Our findings suggest that GGO AMFM is associated with increased systemic inflammation and emphysema progression.

Sliding between lung lobes along lobar fissures is a poorly understood aspect of lung mechanics. The objective of this study was to test the hypothesis that lobar sliding helps reduce distortion in the lung parenchyma during breathing. Finite element models of left lungs with geometries and boundary conditions derived from medical images of human subjects were developed. Effect of lobar sliding was studied by comparing nonlinear finite elastic contact mechanics simulations that allowed and disallowed lobar sliding. Lung parenchymal distortion during simulated breath-holds and tidal breathing was quantified with the model's spatial mean anisotropic deformation index (ADI), a measure of directional preference in volume change that varies spatially in the lung. Models that allowed lobar sliding had significantly lower mean ADI (i.e., lesser parenchymal distortion) than models that disallowed lobar sliding under simulations of both tidal breathing (5.3% median difference, P = 0.008, n = 8) and lung deformation between breath-holds at total lung capacity and functional residual capacity (3.2% median difference, P = 0.03, n = 6). This effect was most pronounced in the lower lobe where lobar sliding reduced parenchymal distortion with statistical significance, but not in the upper lobe. In addition, more lobar sliding was correlated with greater reduction in distortion between sliding and nonsliding models in our study cohorts (Pearson's correlation coefficient of 0.95 for tidal breathing, 0.87 for breath-holds, and 0.91 for the combined dataset). These findings are consistent with the hypothesis that lung lobar sliding reduces parenchymal distortion during breathing. NEW & NOTEWORTHY The role of lobar sliding in lung mechanics is poorly understood. Delineating this role could help explain how breathing is affected by anatomical differences between subjects such as incomplete and missing lobar fissures. We used computational contact mechanics models of lungs from human subjects to delineate the effect of lobar sliding by comparing simulations that allowed and disallowed sliding. We found evidence consistent with the hypothesis that lung lobar sliding reduces parenchymal distortion during breathing.

Background: Chest computed tomography (CT) enables characterization of pulmonary diseases by producing high-resolution and high-contrast images of the intricate lung structures. Deformable image registration is used to align chest CT scans at different lung volumes, yielding estimates of local tissue expansion and contraction., Purpose: We investigated the utility of deep generative models for directly predicting local tissue volume change from lung CT images, bypassing computationally expensive iterative image registration and providing a method that can be utilized in scenarios where either one or two CT scans are available., Methods: A residual regression convolutional neural network, called Reg3DNet+, is proposed for directly regressing high-resolution images of local tissue volume change (i.e., Jacobian) from CT images. Image registration was performed between lung volumes at total lung capacity (TLC) and functional residual capacity (FRC) using a tissue mass- and structure-preserving registration algorithm. The Jacobian image was calculated from the registration-derived displacement field and used as the ground truth for local tissue volume change. Four separate Reg3DNet+ models were trained to predict Jacobian images using a multifactorial study design to compare the effects of network input (i.e., single image vs. paired images) and output space (i.e., FRC vs. TLC). The models were trained and evaluated on image datasets from the COPDGene study. Models were evaluated against the registration-derived Jacobian images using local, regional, and global evaluation metrics., Results: Statistical analysis revealed that both factors - network input and output space - were significant determinants for change in evaluation metrics. Paired-input models performed better than single-input models, and model performance was better in the output space of FRC rather than TLC. Mean structural similarity index for paired-input models was 0.959 and 0.956 for FRC and TLC output spaces, respectively, and for single-input models was 0.951 and 0.937. Global evaluation metrics demonstrated correlation between registration-derived Jacobian mean and predicted Jacobian mean: coefficient of determination (r 2 ) for paired-input models was 0.974 and 0.938 for FRC and TLC output spaces, respectively, and for single-input models was 0.598 and 0.346. After correcting for effort, registration-derived lobar volume change was strongly correlated with the predicted lobar volume change: for paired-input models r 2 was 0.899 for both FRC and TLC output spaces, and for single-input models r 2 was 0.803 and 0.862, respectively., Conclusions: Convolutional neural networks can be used to directly predict local tissue mechanics, eliminating the need for computationally expensive image registration. Networks that use paired CT images acquired at TLC and FRC allow for more accurate prediction of local tissue expansion compared to networks that use a single image. Networks that only require a single input image still show promising results, particularly after correcting for effort, and allow for local tissue expansion estimation in cases where multiple CT scans are not available. For single-input networks, the FRC image is more predictive of local tissue volume change compared to the TLC image., (© 2023 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Computed Tomography (CT) imaging is routinely used for imaging of the lungs. Deep learning can effectively automate complex and laborious tasks in medical imaging. In this work, a deep learning technique is utilized to assess lobar fissure completeness (also known as fissure integrity) from pulmonary CT images. The human lungs are divided into five separate lobes, divided by the lobar fissures. Fissure integrity assessment is important to endobronchial valve treatment screening. Fissure integrity is known to be a biomarker of collateral ventilation between lobes impacting the efficacy of valves designed to block airflow to diseased lung regions. Fissure integrity is also likely to impact lobar sliding which has recently been shown to affect lung biomechanics. Further widescale study of fissure integrity's impact on disease susceptibility and progression requires rapid, reproducible, and noninvasive fissure integrity assessment. In this paper we describe IntegrityNet, an attention U-Net based automatic fissure integrity analysis tool. IntegrityNet is able to predict fissure integrity with an accuracy of 95.8%, 96.1%, and 89.8% for left oblique, right oblique, and right horizontal fissures, compared to manual analysis on a dataset of 82 subjects. We also show that our method is robust to COPD severity and reproducible across subject scans acquired at different time points., (© 2023. Springer Nature Limited.)

Importance: People who smoked cigarettes may experience respiratory symptoms without spirometric airflow obstruction. These individuals are typically excluded from chronic obstructive pulmonary disease (COPD) trials and lack evidence-based therapies., Objective: To define the natural history of persons with tobacco exposure and preserved spirometry (TEPS) and symptoms (symptomatic TEPS)., Design, Setting, and Participants: SPIROMICS II was an extension of SPIROMICS I, a multicenter study of persons aged 40 to 80 years who smoked cigarettes (>20 pack-years) with or without COPD and controls without tobacco exposure or airflow obstruction. Participants were enrolled in SPIROMICS I and II from November 10, 2010, through July 31, 2015, and followed up through July 31, 2021., Exposures: Participants in SPIROMICS I underwent spirometry, 6-minute walk distance testing, assessment of respiratory symptoms, and computed tomography of the chest at yearly visits for 3 to 4 years. Participants in SPIROMICS II had 1 additional in-person visit 5 to 7 years after enrollment in SPIROMICS I. Respiratory symptoms were assessed with the COPD Assessment Test (range, 0 to 40; higher scores indicate more severe symptoms). Participants with symptomatic TEPS had normal spirometry (postbronchodilator ratio of forced expiratory volume in the first second [FEV1] to forced vital capacity >0.70) and COPD Assessment Test scores of 10 or greater. Participants with asymptomatic TEPS had normal spirometry and COPD Assessment Test scores of less than 10. Patient-reported respiratory symptoms and exacerbations were assessed every 4 months via phone calls., Main Outcomes and Measures: The primary outcome was assessment for accelerated decline in lung function (FEV1) in participants with symptomatic TEPS vs asymptomatic TEPS. Secondary outcomes included development of COPD defined by spirometry, respiratory symptoms, rates of respiratory exacerbations, and progression of computed tomographic-defined airway wall thickening or emphysema., Results: Of 1397 study participants, 226 had symptomatic TEPS (mean age, 60.1 [SD, 9.8] years; 134 were women [59%]) and 269 had asymptomatic TEPS (mean age, 63.1 [SD, 9.1] years; 134 were women [50%]). At a median follow-up of 5.76 years, the decline in FEV1 was -31.3 mL/y for participants with symptomatic TEPS vs -38.8 mL/y for those with asymptomatic TEPS (between-group difference, -7.5 mL/y [95% CI, -16.6 to 1.6 mL/y]). The cumulative incidence of COPD was 33.0% among participants with symptomatic TEPS vs 31.6% among those with asymptomatic TEPS (hazard ratio, 1.05 [95% CI, 0.76 to 1.46]). Participants with symptomatic TEPS had significantly more respiratory exacerbations than those with asymptomatic TEPS (0.23 vs 0.08 exacerbations per person-year, respectively; rate ratio, 2.38 [95% CI, 1.71 to 3.31], P < .001)., Conclusions and Relevance: Participants with symptomatic TEPS did not have accelerated rates of decline in FEV1 or increased incidence of COPD vs those with asymptomatic TEPS, but participants with symptomatic TEPS did experience significantly more respiratory exacerbations over a median follow-up of 5.8 years.

Imaging biomarkers can assess disease progression or prognoses and are valuable tools to help guide interventions. Particularly in lung imaging, biomarkers present an opportunity to extract regional information that is more robust to the patient's condition prior to intervention than current gold standard pulmonary function tests (PFTs). This regional aspect has particular use in functional avoidance radiation therapy (RT) in which treatment planning is optimized to avoid regions of high function with the goal of sparing functional lung and improving patient quality of life post-RT. To execute functional avoidance, detailed dose-response models need to be developed to identify regions which should be protected. Previous studies have begun to do this, but for these models to be clinically translated, they need to be validated. This work validates two metrics that encompass the main components of lung function (ventilation and perfusion) through post-mortem histopathology performed in a novel porcine model. With these methods validated, we can use them to study the nuanced radiation-induced changes in lung function and develop more advanced models., (© 2023. The Author(s).)

Background Approximately half of adults with chronic obstructive pulmonary disease (COPD) remain undiagnosed. Chest CT scans are frequently acquired in clinical practice and present an opportunity to detect COPD. Purpose To assess the performance of radiomics features in COPD diagnosis using standard-dose and low-dose CT models. Materials and Methods This secondary analysis included participants enrolled in the Genetic Epidemiology of COPD, or COPDGene, study at baseline (visit 1) and 10 years after baseline (visit 3). COPD was defined by a forced expiratory volume in the 1st second of expiration to forced vital capacity ratio less than 0.70 at spirometry. The performance of demographics, CT emphysema percentage, radiomics features, and a combined feature set derived from inspiratory CT alone was evaluated. CatBoost (Yandex), a gradient boosting algorithm, was used to perform two classification experiments to detect COPD; the two models were trained and tested on standard-dose CT data from visit 1 (model I) and low-dose CT data from visit 3 (model II). Classification performance of the models was evaluated using area under the receiver operating characteristic curve (AUC) and precision-recall curve analysis. Results A total of 8878 participants (mean age, 57 years ± 9 [SD]; 4180 female, 4698 male) were evaluated. Radiomics features in model I achieved an AUC of 0.90 (95% CI: 0.88, 0.91) in the standard-dose CT test cohort versus demographics (AUC, 0.73; 95% CI: 0.71, 0.76; P < .001), emphysema percentage (AUC, 0.82; 95% CI 0.80, 0.84; P < .001), and combined features (AUC, 0.90; 95% CI: 0.89, 0.92; P = .16). Model II, trained on low-dose CT scans, achieved an AUC of 0.87 (95% CI: 0.83, 0.91) on the 20% held-out test set for radiomics features compared with demographics (AUC, 0.70; 95% CI: 0.64, 0.75; P = .001), emphysema percentage (AUC, 0.74; 95% CI: 0.69, 0.79; P = .002), and combined features (AUC, 0.88; 95% CI: 0.85, 0.92; P = .32). Density and texture features were the majority of the top 10 features in the standard-dose model, whereas shape features of lungs and airways were significant contributors in the low-dose CT model. Conclusion A combination of features representing parenchymal texture and lung and airway shape on inspiratory CT scans can be used to accurately detect COPD. ClinicalTrials.gov registration no. NCT00608764 © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Vliegenthart in this issue.

Purpose: To identify metrics of radiation dose delivered to highly ventilated lung that are predictive of radiation-induced pneumonitis., Methods and Materials: A cohort of 90 patients with locally advanced non-small cell lung cancer treated with standard fractionated radiation therapy (RT) (60-66 Gy in 30-33 fractions) were evaluated. Regional lung ventilation was determined from pre-RT 4-dimensional computed tomography (4DCT) using the Jacobian determinant of a B-spline deformable image registration to estimate lung tissue expansion during respiration. Multiple voxel-wise population- and individual-based thresholds for defining high functioning lung were considered. Mean dose and volumes receiving dose ≥ 5-60 Gy were analyzed for both total lung-ITV (MLD,V5-V60) and highly ventilated functional lung-ITV (fMLD,fV5-fV60). The primary endpoint was symptomatic grade 2+ (G2+) pneumonitis. Receiver operator curve (ROC) analyses were used to identify predictors of pneumonitis., Results: G2+ pneumonitis occurred in 22.2% of patients, with no differences between stage, smoking status, COPD, or chemo/immunotherapy use between G<2 and G2+ patients (P≥ 0.18). Highly ventilated lung was defined as voxels exceeding the population-wide median of 18% voxel-level expansion. All total and functional metrics were significantly different between patients with and without pneumonitis (P≤ 0.039). Optimal ROC points predicting pneumonitis from functional lung dose were fMLD ≤ 12.3 Gy, fV5 ≤ 54% and fV20 ≤ 19 %. Patients with fMLD ≤ 12.3 Gy had a 14% risk of developing G2+ pneumonitis whereas risk significantly increased to 35% for those with fMLD > 12.3 Gy (P = 0.035)., Conclusions: Dose to highly ventilated lung is associated with symptomatic pneumonitis and treatment planning strategies should focus on limiting dose to functional regions. These findings provide important metrics to be used in functional lung avoidance RT planning and designing clinical trials., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Background: Functional lung avoidance radiation therapy (RT) is a technique being investigated to preferentially avoid specific regions of the lung that are predicted to be more susceptible to radiation-induced damage. Reducing the dose delivered to high functioning regions may reduce the occurrence radiation-induced lung injuries (RILIs) and toxicities. However, in order to develop effective lung function-sparing plans, accurate predictions of post-RT ventilation change are needed to determine which regions of the lung should be spared., Purpose: To predict pulmonary ventilation change following RT for nonsmall cell lung cancer using machine learning., Methods: A conditional generative adversarial network (cGAN) was developed with data from 82 human subjects enrolled in a randomized clinical trial approved by the institution's IRB to predict post-RT pulmonary ventilation change. The inputs to the network were the pre-RT pulmonary ventilation map and radiation dose distribution. The loss function was a combination of the binary cross-entropy loss and an asymmetrical structural similarity index measure (aSSIM) function designed to increase penalization of under-prediction of ventilation damage. Network performance was evaluated against a previously developed polynomial regression model using a paired sample t-test for comparison. Evaluation was performed using eight-fold cross-validation., Results: From the eight-fold cross-validation, we found that relative to the polynomial model, the cGAN model significantly improved predicting regions of ventilation damage following radiotherapy based on true positive rate (TPR), 0.14±0.15 to 0.72±0.21, and Dice similarity coefficient (DSC), 0.19±0.16 to 0.46±0.14, but significantly declined in true negative rate, 0.97±0.05 to 0.62±0.21, and accuracy, 0.79±0.08 to 0.65±0.14. Additionally, the average true positive volume increased from 104±119 cc in the POLY model to 565±332 cc in the cGAN model, and the average false negative volume decreased from 654±361 cc in the POLY model to 193±163 cc in the cGAN model., Conclusions: The proposed cGAN model demonstrated significant improvement in TPR and DSC. The higher sensitivity of the cGAN model can improve the clinical utility of functional lung avoidance RT by identifying larger volumes of functional lung that can be spared and thus decrease the probability of the patient developing RILIs., (© 2023 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Background: Biomarkers estimating local lung ventilation have been derived from computed tomography (CT) imaging using various image acquisition and post-processing techniques. CT-ventilation biomarkers have potential clinical use in functional avoidance radiation therapy (RT), in which RT treatment plans are optimized to reduce dose delivered to highly-ventilated lung. Widespread clinical implementation of CT-ventilation biomarkers necessitates understanding of biomarker repeatability. Performing imaging within a highly controlled experimental design enables quantification of error associated with remaining variables., Purpose: To characterize CT-ventilation biomarker repeatability and dependence on image acquisition and post-processing methodology in anesthetized and mechanically ventilated pigs., Methods: Five mechanically ventilated Wisconsin Miniature Swine (WMS) received multiple consecutive four-dimensional CT (4DCT) and maximum inhale and exhale breath-hold CT (BH-CT) scans on five dates to generate CT-ventilation biomarkers. Breathing maneuvers were controlled with an average tidal volume difference <200 cc. As surrogates for ventilation, multiple local expansion ratios (LERs) were calculated from the acquired CT scans using Jacobian-based post-processing techniques. LER 2 measured local expansion between an image pair using either inhale and exhale BH-CT images or two 4DCT breathing phase images. LER N measured the maximum local expansion across the 4DCT breathing phase images. Breathing maneuver consistency, intra- and inter-day biomarker repeatability, image acquisition and post-processing technique dependence were quantitatively analyzed., Results: Biomarkers showed strong agreement with voxel-wise Spearman correlation ρ > 0.9 for intra-day repeatability and ρ > 0.8 for all other comparisons, including between image acquisition techniques. Intra- and inter-day repeatability were significantly different (p<0.01). LER 2 and LER N post-processing did not significantly affect intra-day repeatability., Conclusions: 4DCT and BH-CT ventilation biomarkers derived from consecutive scans show strong agreement in controlled experiments with non-human subjects. This article is protected by copyright. All rights reserved., (This article is protected by copyright. All rights reserved.)

Purpose: To quantify the impact of image noise on CT-based lung ventilation biomarkers calculated using Jacobian determinant techniques. Methods: Five mechanically ventilated swine were imaged on a multi-row CT scanner with acquisition parameters of 120 kVp and 0.6 mm slice thickness in static and 4-dimensional CT (4DCT) modes with respective pitches of 1 and 0.09. A range of tube current time product (mAs) values were used to vary image dose. On two dates, subjects received two 4DCTs: one with 10 mAs/rotation (low-dose, high-noise) and one with CT simulation standard of care 100 mAs/rotation (high-dose, low-noise). Additionally, 10 intermediate noise level breath-hold (BHCT) scans were acquired with inspiratory and expiratory lung volumes. Images were reconstructed with and without iterative reconstruction (IR) using 1 mm slice thickness. The Jacobian determinant of an estimated transformation from a B-spline deformable image registration was used to create CT-ventilation biomarkers estimating lung tissue expansion. 24 CT-ventilation maps were generated per subject per scan date: four 4DCT ventilation maps (two noise levels each with and without IR) and 20 BHCT ventilation maps (10 noise levels each with and without IR). Biomarkers derived from reduced dose scans were registered to the reference full dose scan for comparison. Evaluation metrics were gamma pass rate (Γ) with 2 mm distance-to-agreement and 6% intensity criterion, voxel-wise Spearman correlation ( ρ ) and Jacobian ratio coefficient of variation ( CoV JR ). Results: Comparing biomarkers derived from low (CTDI vol = 6.07 mGy) and high (CTDI vol = 60.7 mGy) dose 4DCT scans, mean Γ, ρ and CoV JR values were 93% ± 3%, 0.88 ± 0.03 and 0.04 ± 0.009, respectively. With IR applied, those values were 93% ± 4%, 0.90 ± 0.04 and 0.03 ± 0.003. Similarly, comparisons between BHCT-based biomarkers with variable dose (CTDI vol = 1.35-7.95 mGy) had mean Γ, ρ and CoV JR of 93% ± 4%, 0.97 ± 0.02 and 0.03 ± 0.006 without IR and 93% ± 4%, 0.97 ± 0.03 and 0.03 ± 0.007 with IR. Applying IR did not significantly change any metrics ( p > 0.05). Discussion: This work demonstrated that CT-ventilation, calculated using the Jacobian determinant of an estimated transformation from a B-spline deformable image registration, is invariant to Hounsfield Unit (HU) variation caused by image noise. This advantageous finding may be leveraged clinically with potential applications including dose reduction and/or acquiring repeated low-dose acquisitions for improved ventilation characterization., Competing Interests: JR is a shareholder in VIDA Diagnostics, Inc. GC receives licensing fees from VIDA Diagnostics, Inc. and is a consultant and owns stock in PowerPollen, Inc., Ames IA. JB has ownership interest in MR Guidance, LLC (MRG), which has performed business with radiation oncology venders. The remaining 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., (Copyright © 2023 Flakus, Wuschner, Wallat, Shao, Shanmuganayagam, Christensen, Reinhardt, Li and Bayouth.)

Background: Quantitative CT is becoming increasingly common for the characterisation of lung disease; however, its added potential as a clinical tool for predicting severe exacerbations remains understudied. We aimed to develop and validate quantitative CT-based models for predicting severe chronic obstructive pulmonary disease (COPD) exacerbations., Methods: We analysed the Subpopulations and Intermediate Outcome Measures In COPD Study (SPIROMICS) cohort, a multicentre study done at 12 clinical sites across the USA, of individuals aged 40-80 years from four strata: individuals who never smoked, individuals who smoked but had normal spirometry, individuals who smoked and had mild to moderate COPD, and individuals who smoked and had severe COPD. We used 3-year follow-up data to develop logistic regression classifiers for predicting severe exacerbations. Predictors included age, sex, race, BMI, pulmonary function, exacerbation history, smoking status, respiratory quality of life, and CT-based measures of density gradient texture and airway structure. We externally validated our models in a subset from the Genetic Epidemiology of COPD (COPDGene) cohort. Discriminative model performance was assessed using the area under the receiver operating characteristic curve (AUC), which was also compared with other predictors, including exacerbation history and the BMI, airflow obstruction, dyspnoea, and exercise capacity (BODE) index. We evaluated model calibration using calibration plots and Brier scores., Findings: Participants in SPIROMICS were enrolled between Nov 12, 2010, and July 31, 2015. Participants in COPDGene were enrolled between Jan 10, 2008, and April 15, 2011. We included 1956 participants from the SPIROMICS cohort who had complete 3-year follow-up data: the mean age of the cohort was 63·1 years (SD 9·2) and 1017 (52%) were men and 939 (48%) were women. Among the 1956 participants, 434 (22%) had a history of at least one severe exacerbation. For the CT-based models, the AUC was 0·854 (95% CI 0·852-0·855) for at least one severe exacerbation within 3 years and 0·931 (0·930-0·933) for consistent exacerbations (defined as ≥1 acute episode in each of the 3 years). Models were well calibrated with low Brier scores (0·121 for at least one severe exacerbation; 0·039 for consistent exacerbations). For the prediction of at least one severe event during 3-year follow-up, AUCs were significantly higher with CT biomarkers (0·854 [0·852-0·855]) than exacerbation history (0·823 [0·822-0·825]) and BODE index 0·812 [0·811-0·814]). 6965 participants were included in the external validation cohort, with a mean age of 60·5 years (SD 8·9). In this cohort, AUC for at least one severe exacerbation was 0·768 (0·767-0·769; Brier score 0·088)., Interpretation: CT-based prediction models can be used for identification of patients with COPD who are at high risk of severe exacerbations. The newly identified CT biomarkers could potentially enable investigation into underlying disease mechanisms responsible for exacerbations., Funding: National Institutes of Health and the National Heart, Lung, and Blood Institute., Competing Interests: Declaration of interests EAH has received grants from the National Institutes of Health (NIH) and American Lung Association (ALA); has received royalties from VIDA Diagnostics; is a participant on Siemens photon counting CT advisory board; and is founder and shareholder of VIDA Diagnostics. JG has received grants from NIH and is a shareholder of VIDA Diagnostics. APC has received grants from NIH and is a paid consultant for GlaxoSmithKline and AstraZeneca. JDN has received grants from NIH and VIDA Diagnostics; has received royalties from Elsevier; and has received consulting fees, honoraria for lectures, travel expenses, fees for leadership roles, is also shareholder for, shares multiple patents with, and has received computer equipment from VIDA Diagnostics. PN has received grants from the NIH and honoraria from the GE Medical–University of Washington Imaging Symposium. SF has received grants from the American Thoracic Society, Fisher, and Paykel; and has served as a consultant for Genentech. GEC has received grants from the NIH; royalties from VIDA Diagnostics; fees for consultancy work from PowerPollen; and holds stocks or stock options in PowerPollen. FA has received grants from the NIH. IZB has received grants from Theravance & Viatris; fees for consultancy work from Theravance & Viatris, GlaxoSmithKline, AstraZeneca, and GE Healthcare; payment or honoraria for lectures, presentations, speaking, bureaus, manuscript writing, or educational events from Theravance & Viatris, AstraZeneca, and Grifols; is a participant on a data safety monitoring board for Theravance & Viatris, GlaxoSmithKline, and AstraZeneca; and is a member of the American Thoracic Society pulmonary function test committee. RGB has received grants from the NIH, the National Heart, Lung, and Blood Institute (NHLBI), The COPD Foundation, ALA, and the Foundation for the NIH; has received travel expenses from The COPD Foundation for an NIH-funded study; and has served The COPD Foundation in an unpaid leadership or fiduciary role. SPB has received grants from NIH; royalties from Springer Humana; fees for consultancy work from Boehringer Ingelheim, Sanofi/Regeneron, Sunovion, and GlaxoSmithKline; and holds stock options for Vigor Medical Systems. CBC has received grants from the NIH, The COPD Foundation, and the Foundation for the NIH; royalties from the Cambridge University Press; fees for consultancy work from Nuvaira and MGC Diagnostics; and personal fees from GlaxoSmithKline, and medicolegal personal fees from various law firms. MKH has received grants from the NHLBI, Sanofi, Novartis, Nuvaira, Sunovion; royalties from UpToDate and Norton Publishing; fees for consultancy work from AstraZeneca, Boehringer Ingelheim, GlaxoSmithKline, Novartis, Pulmonx, Teva, Verona, Merck, Sanofi, DevPro, Aerogen, and United Therapeutics; payments or honoraria from Cipla, Chiesi, Astra Zeneca, Boehringer Ingelheim, and GlaxoSmithKline; has participated on a data safety monitoring board or advisory board for Novartis and MedTronic; and is a member of The COPD Foundation Board, The COPD Foundation Scientific Advisory Committee, and ALA advisory committee. FJM has received grants from NHLBI, AstraZeneca, Chiesi, GlaxoSmithKline, and Sanofi/Regeneron; fees for consultancy work from Astra Zeneca, Boehringer Ingelheim, Chiesi, CsL Behring, Gala, GlaxoSmithKline, Novartis, Polarean, PulmonX, Sanofi/Regeneron, Sunovion, Teva, Theravance & Viatris, and Virona; payment or honoraria for lectures, presentations, speaking, bureaus, manuscript writing, or educational events from UpToDate; and is a member of the data safety monitoring board of MedTronic. MGM has received grants from the NIH and NHLBI. VEO is member of an independent data safety monitoring board for Sanofi and Regeneron. RP has received grants from NHLBI, The COPD Foundation, and Department of Veteran Affairs; and has received fees for consultancy work from Partner Therapeutics. PGW has received grants from The COPD Foundation, and Genentech; fees for consultancy work from Glenmark Pharmaceuticals, the University of Wisconsin, NGM Pharma, GlaxoSmithKline, Theravance, Sanofi, and AstraZeneca; and honoraria from the Western Society of Allergy, Asthma and Immunology. JMR has received grants from the NHLBI, and The Roy J Carver Charitable Trust; royalties from VIDA Diagnostics; personal fees from Boehringer Ingelheim; payment for expert testimony from Desmarais LLP; and is a shareholder of VIDA Diagnostics. All other authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY-NC-ND 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Background The prevalence of chronic obstructive pulmonary disease (COPD) in women is fast approaching that in men, and women experience greater symptom burden. Although sex differences in emphysema have been reported, differences in airways have not been systematically characterized. Purpose To evaluate whether structural differences in airways may underlie some of the sex differences in COPD prevalence and clinical outcomes. Materials and Methods In a secondary analyses of a multicenter study of never-, current-, and former-smokers enrolled from January 2008 to June 2011 and followed up longitudinally until November 2020, airway disease on CT images was quantified using seven metrics: airway wall thickness, wall area percent, and square root of the wall thickness of a hypothetical airway with internal perimeter of 10 mm (referred to as Pi10) for airway wall; and lumen diameter, airway volume, total airway count, and airway fractal dimension for airway lumen. Least-squares mean values for each airway metric were calculated and adjusted for age, height, ethnicity, body mass index, pack-years of smoking, current smoking status, total lung capacity, display field of view, and scanner type. In ever-smokers, associations were tested between each airway metric and postbronchodilator forced expiratory volume in 1 second (FEV 1 )-to-forced vital capacity (FVC) ratio, modified Medical Research Council dyspnea scale, St George's Respiratory Questionnaire score, and 6-minute walk distance. Multivariable Cox proportional hazards models were created to evaluate the sex-specific association between each airway metric and mortality. Results In never-smokers ( n = 420), men had thicker airway walls than women as quantified on CT images for segmental airway wall area percentage (least-squares mean, 47.68 ± 0.61 [standard error] vs 45.78 ± 0.55; difference, -1.90; P = .02), whereas airway lumen dimensions were lower in women than men after accounting for height and total lung capacity (segmental lumen diameter, 8.05 mm ± 0.14 vs 9.05 mm ± 0.16; difference, -1.00 mm; P < .001). In ever-smokers ( n = 9363), men had greater segmental airway wall area percentage (least-squares mean, 52.19 ± 0.16 vs 48.89 ± 0.18; difference, -3.30; P < .001), whereas women had narrower segmental lumen diameter (7.80 mm ± 0.05 vs 8.69 mm ± 0.04; difference, -0.89; P < .001). A unit change in each of the airway metrics (higher wall or lower lumen measure) resulted in lower FEV 1 -to-FVC ratio, more dyspnea, poorer respiratory quality of life, lower 6-minute walk distance, and worse survival in women compared with men (all P < .01). Conclusion Airway lumen sizes quantified at chest CT were smaller in women than in men after accounting for height and lung size, and these lower baseline values in women conferred lower reserves against respiratory morbidity and mortality for equivalent changes compared with men. © RSNA, 2022 Online supplemental material is available for this article.

Chronic obstructive pulmonary disease (COPD) is an umbrella term used to define a collection of inflammatory lung diseases that cause airflow obstruction and severe damage to the lung parenchyma. This study investigated the robustness of image-registration-based local biomechanical properties of the lung in individuals with COPD as a function of Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage. Image registration was used to estimate the pointwise correspondences between the inspiration (total lung capacity) and expiration (residual volume) computed tomography (CT) images of the lung for each subject. In total, three biomechanical measures were computed from the correspondence map: the Jacobian determinant; the anisotropic deformation index (ADI); and the slab-rod index (SRI). CT scans from 245 subjects with varying GOLD stages were analyzed from the SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS). Results show monotonic increasing or decreasing trends in the three biomechanical measures as a function of GOLD stage for the entire lung and on a lobe-by-lobe basis. Furthermore, these trends held across all five image registration algorithms. The consistency of the five image registration algorithms on a per individual basis is shown using Bland-Altman plots.

Vessel segmentation in the lung is an ongoing challenge. While many methods have been able to successfully identify vessels in normal, healthy, lungs, these methods struggle in the presence of abnormalities. Following radiotherapy, these methods tend to identify regions of radiographic change due to post-radiation therapytoxicities as vasculature falsely. By combining texture analysis and existing vasculature and masking techniques, we have developed a novel vasculature segmentation workflow that improves specificity in irradiated lung while preserving the sensitivity of detection in the rest of the lung. Furthermore, radiation dose has been shown to cause vascular injury as well as reduce pulmonary function post-RT. This work shows the improvements our novel vascular segmentation method provides relative to existing methods. Additionally, we use this workflow to show a dose dependent radiation-induced change in vasculature which is correlated with previously measured perfusion changes ( R 2 = 0.72) in both directly irradiated and indirectly damaged regions of perfusion. These results present an opportunity to extend non-contrast CT-derived models of functional change following radiation therapy., Competing Interests: JR is a shareholder in VIDA Diagnostics, Inc., GC receives licensing fees from VIDA Diagnostics, Inc., and JB has ownership interest in MR Guidance, LLC. MR Guidance has business activity with ViewRay, Inc., and while this project was not sponsored in any way by ViewRay, data were collected on the ViewRay MRIdian system. Data were collected on a Radixact system (Accuray, Inc.) provided to UW-Madison under a research agreement (JB, PI) The remaining 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., (Copyright © 2022 Wuschner, Flakus, Wallat, Reinhardt, Shanmuganayagam, Christensen, Gerard and Bayouth.)

Background Clustering key clinical characteristics of participants in the Severe Asthma Research Program (SARP), a large, multicenter prospective observational study of patients with asthma and healthy controls, has led to the identification of novel asthma phenotypes. Purpose To determine whether quantitative CT (qCT) could help distinguish between clinical asthma phenotypes. Materials and Methods A retrospective cross-sectional analysis was conducted with the use of qCT images (maximal bronchodilation at total lung capacity [TLC], or inspiration, and functional residual capacity [FRC], or expiration) from the cluster phenotypes of SARP participants (cluster 1: minimal disease; cluster 2: mild, reversible; cluster 3: obese asthma; cluster 4: severe, reversible; cluster 5: severe, irreversible) enrolled between September 2001 and December 2015. Airway morphometry was performed along standard paths (RB1, RB4, RB10, LB1, and LB10). Corresponding voxels from TLC and FRC images were mapped with use of deformable image registration to characterize disease probability maps (DPMs) of functional small airway disease (fSAD), voxel-level volume changes (Jacobian), and isotropy (anisotropic deformation index [ADI]). The association between cluster assignment and qCT measures was evaluated using linear mixed models. Results A total of 455 participants were evaluated with cluster assignments and CT (mean age ± SD, 42.1 years ± 14.7; 270 women). Airway morphometry had limited ability to help discern between clusters. DPM fSAD was highest in cluster 5 (cluster 1 in SARP III: 19.0% ± 20.6; cluster 2: 18.9% ± 13.3; cluster 3: 24.9% ± 13.1; cluster 4: 24.1% ± 8.4; cluster 5: 38.8% ± 14.4; P < .001). Lower whole-lung Jacobian and ADI values were associated with greater cluster severity. Compared to cluster 1, cluster 5 lung expansion was 31% smaller (Jacobian in SARP III cohort: 2.31 ± 0.6 vs 1.61 ± 0.3, respectively, P < .001) and 34% more isotropic (ADI in SARP III cohort: 0.40 ± 0.1 vs 0.61 ± 0.2, P < .001). Within-lung Jacobian and ADI SDs decreased as severity worsened (Jacobian SD in SARP III cohort: 0.90 ± 0.4 for cluster 1; 0.79 ± 0.3 for cluster 2; 0.62 ± 0.2 for cluster 3; 0.63 ± 0.2 for cluster 4; and 0.41 ± 0.2 for cluster 5; P < .001). Conclusion Quantitative CT assessments of the degree and intraindividual regional variability of lung expansion distinguished between well-established clinical phenotypes among participants with asthma from the Severe Asthma Research Program study. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Verschakelen in this issue.

Recent functional lung imaging studies have presented evidence of an “indirect effect” on perfusion damage, where regions that are unirradiated or lowly irradiated but that are supplied by highly irradiated regions observe perfusion damage post-radiation therapy (RT). The purpose of this work was to investigate this effect using a contrast-enhanced dynamic CT protocol to measure perfusion change in five novel swine subjects. A cohort of five Wisconsin Miniature Swine (WMS) were given a research course of 60 Gy in five fractions delivered locally to a vessel in the lung using an Accuray Radixact tomotherapy system with Synchrony motion tracking to increase delivery accuracy. Imaging was performed prior to delivering RT and 3 months post-RT to yield a 28−36 frame image series showing contrast flowing in and out of the vasculature. Using MIM software, contours were placed in six vessels on each animal to yield a contrast flow curve for each vessel. The contours were placed as follows: one at the point of max dose, one low-irradiated (5−20 Gy) branching from the max dose vessel, one low-irradiated (5−20 Gy) not branching from the max dose vessel, one unirradiated (<5 Gy) branching from the max dose vessel, one unirradiated (<5 Gy) not branching from the max dose vessel, and one in the contralateral lung. Seven measurements (baseline-to-baseline time and difference, slope up and down, max rise and value, and area under the curve) were acquired for each vessel’s contrast flow curve in each subject. Paired Student t-tests showed statistically significant (p < 0.05) reductions in the area under the curve in the max dose, and both fed contours indicating an overall reduction in contrast in these regions. Additionally, there were statistically significant reductions observed when comparing pre- and post-RT in slope up and down in the max dose, low-dose fed, and no-dose fed contours but not the low-dose not-fed, no-dose not-fed, or contralateral contours. These findings suggest an indirect damage effect where irradiation of the vasculature causes a reduction in perfusion in irradiated regions as well as regions fed by the irradiated vasculature.

This paper presents the Population Learning followed by One Shot Learning (PLOSL) pulmonary image registration method. PLOSL is a fast unsupervised learning-based framework for 3D-CT pulmonary image registration algorithm based on combining population learning (PL) and one-shot learning (OSL). The PLOSL image registration has the advantages of the PL and OSL approaches while reducing their respective drawbacks. The advantages of PLOSL include improved performance over PL, substantially reducing OSL training time and reducing the likelihood of OSL getting stuck in local minima. PLOSL pulmonary image registration uses tissue volume preserving and vesselness constraints for registration of inspiration-to-expiration and expiration-to-inspiration pulmonary CT images. A coarse-to-fine convolution encoder-decoder CNN architecture is used to register large and small shape features. During training, the sum of squared tissue volume difference (SSTVD) compensates for intensity differences between inspiration and expiration computed tomography (CT) images and the sum of squared vesselness measure difference (SSVMD) helps match the lung vessel tree. Results show that the PLOSL (SSTVD+SSVMD) algorithm achieved subvoxel landmark error while preserving pulmonary topology on the SPIROMICS data set, the public DIR-LAB COPDGene and 4DCT data sets., Competing Interests: Declaration of Competing Interest The following authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Di Wang, Yue Pan, Oguz C. Durumeric, Joyce D Schroeder The following authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Joseph M. Reinhardt: Joseph M. Reinhardt is a shareholder in VIDA Diagnostics, Inc. and has received personal fees from Desmarais LLP and Boehringer Ingelheim GmbH. Eric A. Hoffman: Eric A. Hoffman is a founder and shareholder of VIDA Diagnostics, a company commercializing lung image analysis software developed, in part, at the University of Iowa. Gary E. Christensen: Gary E Christensen receives licensing fees from VIDA Diagnostics, Inc. and is a consultant for PowerPollen, Ankeny, IA., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Purpose. To investigate indirect radiation-induced changes in airways as precursors to atelectasis post radiation therapy (RT). Methods. Three Wisconsin Miniature Swine (WMS TM ) underwent a research course of 60 Gy in 5 fractions delivered to a targeted airway/vessel in the inferior left lung. The right lung received a max point dose <5 Gy. Airway segmentation was performed on the pre- and three months post-RT maximum inhale phase of the four-dimensional (4D) computed tomography (CT) scans. Changes in luminal area (Ai) and square root of wall area (WA) for each airway were investigated. Changes in ventilation were assessed using the Jacobian ratio and were measured in three different regions: the inferior left lung <5 Gy (ILL), the superior left lung <5 Gy (SLL), and the contralateral right lung <5 Gy (RL). Results. Airways (n = 25) in the right lung for all swine showed no significant changes (p = 0.48) in Ai post-RT compared to pre-RT. Airways (n = 28) in the left lung of all swine were found to have a significant decrease (p < 0.001) in Ai post-RT compared to pre-RT, correlated (Pearson R = -0.97) with airway dose. Additionally,WAdecreased significantly (p < 0.001) with airway dose. Lastly, the Jacobian ratio of the ILL (0.883) was lower than that of the SLL (0.932) and the RL (0.955). Conclusions. This work shows that for the swine analyzed, there were significant correlations between Ai andWAchange with radiation dose. Additionally, there was a decrease in lung function in the regions of the lung supplied by the irradiated airways compared to the regions supplied by unirradiated airways. These results support the hypothesis that airway dose should be considered during treatment planning in order to potentially preserve functional lung and reduce lung toxicities., (© 2021 IOP Publishing Ltd.)

Enhanced intrapulmonary gas transport enables oscillatory ventilation modalities to support gas exchange using extremely low tidal volumes at high frequencies. However, it is unknown whether gas transport rates can be improved by combining multiple frequencies of oscillation simultaneously. The goal of this study was to investigate distributed gas transport in vivo during multi-frequency oscillatory ventilation (MFOV) as compared with conventional mechanical ventilation (CMV) or high-frequency oscillatory ventilation (HFOV). We hypothesized that MFOV would result in more uniform rates of gas transport compared to HFOV, measured using contrast-enhanced CT imaging during wash-in of xenon gas. In 13 pigs, xenon wash-in equilibration rates were comparable between CMV and MFOV, but 21 to 39% slower for HFOV. By contrast, the root-mean-square delivered volume was lowest for MFOV, increased by 70% during HFOV and 365% during CMV. Overall gas transport heterogeneity was similar across all modalities, but gravitational gradients and regional patchiness of specific ventilation contributed to regional ventilation heterogeneity, depending on ventilator modality. We conclude that MFOV combines benefits of low lung stretch, similar to HFOV, but with fast rates of gas transport, similar to CMV., (© 2021. Biomedical Engineering Society.)

Pulmonary respiratory motion artifacts are common in four-dimensional computed tomography (4DCT) of lungs and are caused by missing, duplicated, and misaligned image data. This paper presents a geodesic density regression (GDR) algorithm to correct motion artifacts in 4DCT by correcting artifacts in one breathing phase with artifact-free data from corresponding regions of other breathing phases. The GDR algorithm estimates an artifact-free lung template image and a smooth, dense, 4D (space plus time) vector field that deforms the template image to each breathing phase to produce an artifact-free 4DCT scan. Correspondences are estimated by accounting for the local tissue density change associated with air entering and leaving the lungs, and using binary artifact masks to exclude regions with artifacts from image regression. The artifact-free lung template image is generated by mapping the artifact-free regions of each phase volume to a common reference coordinate system using the estimated correspondences and then averaging. This procedure generates a fixed view of the lung with an improved signal-to-noise ratio. The GDR algorithm was evaluated and compared to a state-of-the-art geodesic intensity regression (GIR) algorithm using simulated CT time-series and 4DCT scans with clinically observed motion artifacts. The simulation shows that the GDR algorithm has achieved significantly more accurate Jacobian images and sharper template images, and is less sensitive to data dropout than the GIR algorithm. We also demonstrate that the GDR algorithm is more effective than the GIR algorithm for removing clinically observed motion artifacts in treatment planning 4DCT scans. Our code is freely available at https://github.com/Wei-Shao-Reg/GDR., Competing Interests: Declaration of Competing Interest John E. Bayouth has ownership interest in MR Guidance, LLC, which has business activity with a company that utilizes deformable image registration (ViewRay, Inc.). Gary E. Christensen receives licensing fees from VIDA Diagnostics, Inc. Other authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier B.V. All rights reserved.)

This study reports systematic longitudinal pathophysiology of lung parenchymal and vascular effects of asymptomatic COVID-19 pneumonia in a young, healthy never-smoking male. Inspiratory and expiratory noncontrast along with contrast dual-energy computed tomography (DECT) scans of the chest were performed at baseline on the day of acute COVID-19 diagnosis ( day 0 ), and across a 90-day period. Despite normal vital signs and pulmonary function tests on the day of diagnosis, the CT scans and corresponding quantification metrics detected abnormalities in parenchymal expansion based on image registration, ground-glass (GGO) texture (inflammation) as well as DECT-derived pulmonary blood volume (PBV). Follow-up scans on day 30 showed improvement in the lung parenchymal mechanics as well as reduced GGO and improved PBV distribution. Improvements in lung PBV continued until day 90 . However, the heterogeneity of parenchymal mechanics and texture-derived GGO increased on days 60 and 90 . We highlight that even asymptomatic COVID-19 infection with unremarkable vital signs and pulmonary function tests can have measurable effects on lung parenchymal mechanics and vascular pathophysiology, which may follow apparently different clinical courses. For this asymptomatic subject, post COVID-19 regional mechanics demonstrated persistent increased heterogeneity concomitant with return of elevated GGOs, despite early improvements in vascular derangement. NEW & NOTEWORTHY We characterized the temporal changes of lung parenchyma and microvascular pathophysiology from COVID-19 infection in an asymptomatic young, healthy nonsmoking male using dual-energy CT. Lung parenchymal mechanics and microvascular disease followed different clinical courses. Heterogeneous perfused blood volume became more uniform on follow-up visits up to 90 days. However, post COVID-19 mechanical heterogeneity of the lung parenchyma increased after apparent improvements in vascular abnormalities, even with normal spirometric indices.

Rationale : Intratidal changes in regional lung aeration, as assessed with dynamic four-dimensional computed tomography (CT; 4DCT), may indicate the processes of recruitment and derecruitment, thus portending atelectrauma during mechanical ventilation. In this study, we characterized the time constants associated with deaeration during the expiratory phase of pressure-controlled ventilation in pigs before and after acute lung injury using respiratory-gated 4DCT and image registration. Methods : Eleven pigs were mechanically ventilated in pressure-controlled mode under baseline conditions and following an oleic acid model of acute lung injury. Dynamic 4DCT scans were acquired without interrupting ventilation. Automated segmentation of lung parenchyma was obtained by a convolutional neural network. Respiratory structures were aligned using 4D image registration. Exponential regression was performed on the time-varying CT density in each aligned voxel during exhalation, resulting in regional estimates of intratidal aeration change and deaeration time constants. Regressions were also performed for regional and total exhaled gas volume changes. Results : Normally and poorly aerated lung regions demonstrated the largest median intratidal aeration changes during exhalation, compared to minimal changes within hyper- and non-aerated regions. Following lung injury, median time constants throughout normally aerated regions within each subject were greater than respective values for poorly aerated regions. However, parametric response mapping revealed an association between larger intratidal aeration changes and slower time constants. Lower aeration and faster time constants were observed for the dependent lung regions in the supine position. Regional gas volume changes exhibited faster time constants compared to regional density time constants, as well as better correspondence to total exhaled volume time constants. Conclusion : Mechanical time constants based on exhaled gas volume underestimate regional aeration time constants. After lung injury, poorly aerated regions experience larger intratidal changes in aeration over shorter time scales compared to normally aerated regions. However, the largest intratidal aeration changes occur over the longest time scales within poorly aerated regions. These dynamic 4DCT imaging data provide supporting evidence for the susceptibility of poorly aerated regions to ventilator-induced lung injury, and for the functional benefits of short exhalation times during mechanical ventilation of injured lungs., Competing Interests: JH and DK are co-founders and shareholders of OscillaVent, Inc. and consultants for ZOLL Medical Corporation. JR and EH are co-founders and shareholders of VIDA Diagnostics, Inc., and GC is paid licensing fees from VIDA Diagnostics, Inc., (Copyright © 2021 Herrmann, Gerard, Shao, Xin, Cereda, Reinhardt, Christensen, Hoffman and Kaczka.)

To analyze radiation induced changes in Hounsfield units and determine their correlation with changes in perfusion and ventilation. Additionally, to compare the post-RT changes in human subjects to those measured in a swine model used to quantify perfusion changes and validate their use as a preclinical model. A cohort of 5 Wisconsin Miniature Swine (WMS) were studied. Additionally, 19 human subjects were recruited as part of an IRB approved clinical trial studying functional avoidance radiation therapy for lung cancer and were treated with SBRT. Imaging (a contrast enhanced dynamic perfusion CT in the swine and 4DCT in the humans) was performed prior to and post-RT. Jacobian elasticity maps were calculated on all 4DCT images. Contours were created from the isodose lines to discretize analysis into 10 Gy dose bins. B-spline deformable image registration allowed for voxel-by-voxel comparative analysis in these contours between timepoints. The WMS underwent a research course of 60 Gy in 5 fractions delivered locally to a target in the lung using an MRI-LINAC system. In the WMS subjects, the dose-bin contours were copied onto the contralateral lung, which received < 5 Gy for comparison. Changes in HU and changes in Jacobian were analyzed in these contours. Statistically significant (p < 0.05) changes in the mean HU value post-RT compared to pre-RT were observed in both the human and WMS groups at all timepoints analyzed. The HU increased linearly with dose for both groups. Strong linear correlation was observed between the changes seen in the swine and humans (Pearson coefficient > 0.97, p < 0.05) at all timepoints. Changes seen in the swine closely modeled the changes seen in the humans at 12 months post RT (slope = 0.95). Jacobian analysis showed between 30 and 60% of voxels were damaged post-RT. Perfusion analysis in the swine showed a statistically significant (p < 0.05) reduction in contrast inside the vasculature 3 months post-RT compared to pre-RT. The increases in contrast outside the vasculature was strongly correlated (Pearson Correlation 0.88) with the reduction in HU inside the vasculature but were not correlated with the changes in Jacobians. Radiation induces changes in pulmonary anatomy at 3 months post-RT, with a strong linear correlation with dose. The change in HU seen in the non-vessel lung parenchyma suggests this metric is a potential biomarker for change in perfusion. Finally, this work suggests that the WMS swine model is a promising pre-clinical model for analyzing radiation-induced changes in humans and poses several benefits over conventional swine models.

Rationale: Airway remodeling in chronic obstructive pulmonary disease (COPD) is due to luminal narrowing and/or loss of airways. Existing computed tomographic metrics of airway disease reflect only components of these processes. With progressive airway narrowing, the ratio of the airway luminal surface area to volume (SA/V) should increase, and with predominant airway loss, SA/V should decrease. Objectives: To phenotype airway remodeling in COPD. Methods: We analyzed the airway trees of 4,325 subjects with COPD Global Initiative for Chronic Obstructive Lung Disease stages 0 to 4 and 73 nonsmokers enrolled in the multicenter COPDGene (Genetic Epidemiology of COPD) cohort. Surface area and volume measurements were estimated for the subtracheal airway tree to derive SA/V. We performed multivariable regression analyses to test associations between SA/V and lung function, 6-minute-walk distance, St. George's Respiratory Questionnaire, change in FEV 1 , and mortality, adjusting for demographics, total airway count, airway wall thickness, and emphysema. On the basis of the change in SA/V over 5 years, we categorized subjects into predominant airway narrowing [positive ∆(SA/V) more than 0] and predominant airway loss [negative ∆(SA/V) less than 0] and compared survival between the two groups. Measurements and Main Results: Airway SA/V was independently associated with FEV 1 /FVC (β = 0.12; 95% confidence interval [CI], 0.09-0.14; P  < 0.001) and FEV 1 % predicted (β = 20.10; 95% CI, 15.13-25.08; P  < 0.001). Airway SA/V was also independently associated with 6-minute-walk distance, respiratory quality of life, and lung function decline. Compared with subjects with predominant airway narrowing ( n  = 2,914; 66.3%), those with predominant airway loss ( n  = 1,484; 33.7%) had worse survival (adjusted hazard ratio for all-cause mortality = 1.58; 95% CI, 1.18-2.13; P  = 0.002). Conclusions: Computed tomography-based airway SA/V is an imaging biomarker of airway remodeling and provides differential information on predominant airway narrowing and loss in COPD. SA/V is associated with respiratory morbidity, lung function decline, and survival.

The purpose of this study was to develop a fully-automated segmentation algorithm, robust to various density enhancing lung abnormalities, to facilitate rapid quantitative analysis of computed tomography images. A polymorphic training approach is proposed, in which both specifically labeled left and right lungs of humans with COPD, and nonspecifically labeled lungs of animals with acute lung injury, were incorporated into training a single neural network. The resulting network is intended for predicting left and right lung regions in humans with or without diffuse opacification and consolidation. Performance of the proposed lung segmentation algorithm was extensively evaluated on CT scans of subjects with COPD, confirmed COVID-19, lung cancer, and IPF, despite no labeled training data of the latter three diseases. Lobar segmentations were obtained using the left and right lung segmentation as input to the LobeNet algorithm. Regional lobar analysis was performed using hierarchical clustering to identify radiographic subtypes of COVID-19. The proposed lung segmentation algorithm was quantitatively evaluated using semi-automated and manually-corrected segmentations in 87 COVID-19 CT images, achieving an average symmetric surface distance of [Formula: see text] mm and Dice coefficient of [Formula: see text]. Hierarchical clustering identified four radiographical phenotypes of COVID-19 based on lobar fractions of consolidated and poorly aerated tissue. Lower left and lower right lobes were consistently more afflicted with poor aeration and consolidation. However, the most severe cases demonstrated involvement of all lobes. The polymorphic training approach was able to accurately segment COVID-19 cases with diffuse consolidation without requiring COVID-19 cases for training.

BACKGROUNDCurrently recommended traditional spirometry outputs do not reflect the relative contributions of emphysema and airway disease to airflow obstruction. We hypothesized that machine-learning algorithms can be trained on spirometry data to identify these structural phenotypes.METHODSParticipants enrolled in a large multicenter study (COPDGene) were included. The data points from expiratory flow-volume curves were trained using a deep-learning model to predict structural phenotypes of chronic obstructive pulmonary disease (COPD) on CT, and results were compared with traditional spirometry metrics and an optimized random forest classifier. Area under the receiver operating characteristic curve (AUC) and weighted F-score were used to measure the discriminative accuracy of a fully convolutional neural network, random forest, and traditional spirometry metrics to phenotype CT as normal, emphysema-predominant (>5% emphysema), airway-predominant (Pi10 > median), and mixed phenotypes. Similar comparisons were made for the detection of functional small airway disease phenotype (>20% on parametric response mapping).RESULTSAmong 8980 individuals, the neural network was more accurate in discriminating predominant emphysema/airway phenotypes (AUC 0.80, 95%CI 0.79-0.81) compared with traditional measures of spirometry, FEV1/FVC (AUC 0.71, 95%CI 0.69-0.71), FEV1% predicted (AUC 0.70, 95%CI 0.68-0.71), and random forest classifier (AUC 0.78, 95%CI 0.77-0.79). The neural network was also more accurate in discriminating predominant emphysema/small airway phenotypes (AUC 0.91, 95%CI 0.90-0.92) compared with FEV1/FVC (AUC 0.80, 95%CI 0.78-0.82), FEV1% predicted (AUC 0.83, 95%CI 0.80-0.84), and with comparable accuracy with random forest classifier (AUC 0.90, 95%CI 0.88-0.91).CONCLUSIONSStructural phenotypes of COPD can be identified from spirometry using deep-learning and machine-learning approaches, demonstrating their potential to identify individuals for targeted therapies.TRIAL REGISTRATIONClinicalTrials.gov NCT00608764.FUNDINGThis study was supported by NIH grants K23 HL133438 and R21EB027891 and an American Thoracic Foundation 2018 Unrestricted Research Grant. The COPDGene study is supported by NIH grants NHLBI U01 HL089897 and U01 HL089856. The COPDGene study (NCT00608764) is also supported by the COPD Foundation through contributions made to an Industry Advisory Committee comprising AstraZeneca, Boehringer-Ingelheim, GlaxoSmithKline, Novartis, and Sunovion.

Purpose: To create a dose-response model that predicts lung ventilation change following radiation therapy, and examine the effects of out-of-phase ventilation., Methods: The dose-response model was built using 27 human subjects who underwent radiation therapy (RT) from an IRB-approved trial. For each four-dimensional computed tomography, two ventilation maps were created by calculating the N-phase local expansion ratio (LER N ) using most or all breathing phases and the 2-phase LER (LER 2 ) using only the end inspiration and end expiration breathing phases. A polynomial regression model was created using the LER N ventilation maps pre-RT and post-RT and dose distributions for each subject, and crossvalidated with a leave-one-out method. Further validation of the model was performed using 15 additional human subjects using common statistical operating characteristics and gamma pass rates., Results: For voxels receiving 20 Gy or greater, there was a significant increase from 52% to 59% (P = 0.03) in the gamma pass rates of the LER N model predicted post-RT Jacobian maps to the actual post-RT Jacobian maps, relative to the LER 2 model. Additionally, accuracy significantly increased (P = 0.03) from 68% to 75% using the LER N model, relative to the LER 2 model., Conclusions: The LER N model was significantly more accurate than the LER 2 model at predicting post-RT ventilation maps. More accurate post-RT ventilation maps will aid in producing a higher quality functional avoidance treatment plan, allowing for potentially better normal tissue sparing., (© 2020 American Association of Physicists in Medicine.)

Out-of-phase ventilation occurs when local regions of the lung reach their maximum or minimum volumes at breathing phases other than the global end inhalation or exhalation phases. This paper presents the N-phase local expansion ratio (LER N ) as a surrogate for lung ventilation. A common approach to estimate lung ventilation is to use image registration to align the end exhalation and inhalation 3DCT images and then analyze the resulting correspondence map. This 2-phase local expansion ratio (LER 2 ) is limited because it ignores out-of-phase ventilation and thus may underestimate local lung ventilation. To overcome this limitation, LER N measures the maximum ratio of local expansion and contraction over the entire breathing cycle. Comparing LER 2 to LER N provides a means for detecting and characterizing locations of the lung that experience out-of-phase ventilation. We present a novel in-phase/out-of-phase ventilation (IOV) function plot to visualize and measure the amount of high-function IOV that occurs during a breathing cycle. Treatment planning 4DCT scans collected during coached breathing from 32 human subjects with lung cancer were analyzed in this study. Results show that out-of-phase breathing occurred in all subjects and that the spatial distribution of out-of-phase ventilation varied from subject to subject. For the 32 subjects analyzed, 50% of the out-of-phase regions on average were mislabeled as low-function by LER 2 (high-function threshold of 1.1, IOV threshold of 1.05). 4DCT and Xenon-enhanced CT of four sheep showed that LER 8 is more accurate than LER 2 for measuring lung ventilation.

Mechanical ventilation strategies that reduce the heterogeneity of regional lung stress and strain may reduce the risk of ventilator-induced lung injury (VILI). In this study, we used registration of four-dimensional computed tomographic (4DCT) images to assess regional lung aeration and deformation in 10 pigs under baseline conditions and following acute lung injury induced with oleic acid. CT images were obtained via dynamic axial imaging (Siemens SOMATOM Force) during conventional pressure-controlled mechanical ventilation (CMV), as well as high-frequency and multi-frequency oscillatory ventilation modalities (HFOV and MFOV, respectively). Our results demonstrate that oscillatory modalities reduce intratidal strain throughout the lung in comparison to conventional ventilation, as well as the spatial gradients of dynamic strain along the dorsal-ventral axis. Harmonic distortion of parenchymal deformation was observed during HFOV with a single discrete sinusoid delivered at the airway opening, suggesting inherent mechanical nonlinearity of the lung tissues. MFOV may therefore provide improved lung-protective ventilation by reducing strain magnitudes and spatial gradients of strain compared to either CMV or HFOV., (Copyright © 2020 Herrmann, Gerard, Shao, Hawley, Reinhardt, Christensen, Hoffman and Kaczka.)

Segmentation of lungs with acute respiratory distress syndrome (ARDS) is a challenging task due to diffuse opacification in dependent regions which results in little to no contrast at the lung boundary. For segmentation of severely injured lungs, local intensity and texture information, as well as global contextual information, are important factors for consistent inclusion of intrapulmonary structures. In this study, we propose a deep learning framework which uses a novel multi-resolution convolutional neural network (ConvNet) for automated segmentation of lungs in multiple mammalian species with injury models similar to ARDS. The multi-resolution model eliminates the need to tradeoff between high-resolution and global context by using a cascade of low-resolution to high-resolution networks. Transfer learning is used to accommodate the limited number of training datasets. The model was initially pre-trained on human CT images, and subsequently fine-tuned on canine, porcine, and ovine CT images with lung injuries similar to ARDS. The multi-resolution model was compared to both high-resolution and low-resolution networks alone. The multi-resolution model outperformed both the low- and high-resolution models, achieving an overall mean Jacaard index of 0.963 ± 0.025 compared to 0.919 ± 0.027 and 0.950 ± 0.036, respectively, for the animal dataset (N=287). The multi-resolution model achieves an overall average symmetric surface distance of 0.438 ± 0.315 mm, compared to 0.971 ± 0.368 mm and 0.657 ± 0.519 mm for the low-resolution and high-resolution models, respectively. We conclude that the multi-resolution model produces accurate segmentations in severely injured lungs, which is attributed to the inclusion of both local and global features., (Copyright © 2019. Published by Elsevier B.V.)

Background: Chronic obstructive pulmonary disease (COPD) remains a major cause of morbidity and mortality. Present-day diagnostic criteria are largely based solely on spirometric criteria. Accumulating evidence has identified a substantial number of individuals without spirometric evidence of COPD who suffer from respiratory symptoms and/or increased morbidity and mortality. There is a clear need for an expanded definition of COPD that is linked to physiologic, structural (computed tomography [CT]) and clinical evidence of disease. Using data from the COPD Genetic Epidemiology study (COPDGene ® ), we hypothesized that an integrated approach that includes environmental exposure, clinical symptoms, chest CT imaging and spirometry better defines disease and captures the likelihood of progression of respiratory obstruction and mortality., Methods: Four key disease characteristics - environmental exposure (cigarette smoking), clinical symptoms (dyspnea and/or chronic bronchitis), chest CT imaging abnormalities (emphysema, gas trapping and/or airway wall thickening), and abnormal spirometry - were evaluated in a group of 8784 current and former smokers who were participants in COPDGene ® Phase 1. Using these 4 disease characteristics, 8 categories of participants were identified and evaluated for odds of spirometric disease progression (FEV 1 > 350 ml loss over 5 years), and the hazard ratio for all-cause mortality was examined., Results: Using smokers without symptoms, CT imaging abnormalities or airflow obstruction as the reference population, individuals were classified as Possible COPD, Probable COPD and Definite COPD. Current Global initiative for obstructive Lung Disease (GOLD) criteria would diagnose 4062 (46%) of the 8784 study participants with COPD. The proposed COPDGene ® 2019 diagnostic criteria would add an additional 3144 participants. Under the new criteria, 82% of the 8784 study participants would be diagnosed with Possible, Probable or Definite COPD. These COPD groups showed increased risk of disease progression and mortality. Mortality increased in patients as the number of their COPD characteristics increased, with a maximum hazard ratio for all cause-mortality of 5.18 (95% confidence interval [CI]: 4.15-6.48) in those with all 4 disease characteristics., Conclusions: A substantial portion of smokers with respiratory symptoms and imaging abnormalities do not manifest spirometric obstruction as defined by population normals. These individuals are at significant risk of death and spirometric disease progression. We propose to redefine the diagnosis of COPD through an integrated approach using environmental exposure, clinical symptoms, CT imaging and spirometric criteria. These expanded criteria offer the potential to stimulate both current and future interventions that could slow or halt disease progression in patients before disability or irreversible lung structural changes develop., Competing Interests: The COPDGene® study is funded by National Heart, Lung, and Blood Institute grants U01 HL089897 and U01 HL089856. The COPDGene® study (NCT00608764) is also supported by the COPD Foundation through contributions made to an Industry Advisory Committee comprised of AstraZeneca, Boehringer-Ingelheim, Genentech, GlaxoSmithKline, Novartis, Pfizer, Siemens, and Sunovion. While some individual authors of this manuscript were employed by one of the listed funders at the time the work of this study was conducted, these employment relationships did not constitute undue influence by funders. These funders have had no official role in the collection, management, analysis and interpretation of the data or design and conduct of the study. All authors have completed a Conflict of Interest form, disclosing any real or apparent financial relationships including receiving royalties, honoraria or fees for consulting, lectures, speakers’ bureaus, continuing education, medical advisory boards or expert testimony; receipt of grants; travel reimbursement; direct employment compensation. These disclosure forms have been filed with the Chronic Obstructive Pulmonary Diseases: Journal of the COPD Foundation Editorial Office and are available for review, upon request, at COPDC@njhealth.org., (JCOPDF © 2019.)

Purpose: CT ventilation imaging (CTVI) is being used to achieve functional avoidance lung cancer radiation therapy in three clinical trials (NCT02528942, NCT02308709, NCT02843568). To address the need for common CTVI validation tools, we have built the Ventilation And Medical Pulmonary Image Registration Evaluation (VAMPIRE) Dataset, and present the results of the first VAMPIRE Challenge to compare relative ventilation distributions between different CTVI algorithms and other established ventilation imaging modalities., Methods: The VAMPIRE Dataset includes 50 pairs of 4DCT scans and corresponding clinical or experimental ventilation scans, referred to as reference ventilation images (RefVIs). The dataset includes 25 humans imaged with Galligas 4DPET/CT, 21 humans imaged with DTPA-SPECT, and 4 sheep imaged with Xenon-CT. For the VAMPIRE Challenge, 16 subjects were allocated to a training group (with RefVI provided) and 34 subjects were allocated to a validation group (with RefVI blinded). Seven research groups downloaded the Challenge dataset and uploaded CTVIs based on deformable image registration (DIR) between the 4DCT inhale/exhale phases. Participants used DIR methods broadly classified into B-splines, Free-form, Diffeomorphisms, or Biomechanical modeling, with CT ventilation metrics based on the DIR evaluation of volume change, Hounsfield Unit change, or various hybrid approaches. All CTVIs were evaluated against the corresponding RefVI using the voxel-wise Spearman coefficient r S , and Dice similarity coefficients evaluated for low function lung ( DSC low ) and high function lung ( DSC high )., Results: A total of 37 unique combinations of DIR method and CT ventilation metric were either submitted by participants directly or derived from participant-submitted DIR motion fields using the in-house software, VESPIR. The r S and DSC results reveal a high degree of inter-algorithm and intersubject variability among the validation subjects, with algorithm rankings changing by up to ten positions depending on the choice of evaluation metric. The algorithm with the highest overall cross-modality correlations used a biomechanical model-based DIR with a hybrid ventilation metric, achieving a median (range) of 0.49 (0.27-0.73) for r S , 0.52 (0.36-0.67) for DSC low , and 0.45 (0.28-0.62) for DSC high . All other algorithms exhibited at least one negative r S value, and/or one DSC value less than 0.5., Conclusions: The VAMPIRE Challenge results demonstrate that the cross-modality correlation between CTVIs and the RefVIs varies not only with the choice of CTVI algorithm but also with the choice of RefVI modality, imaging subject, and the evaluation metric used to compare relative ventilation distributions. This variability may arise from the fact that each of the different CTVI algorithms and RefVI modalities provides a distinct physiologic measurement. Ultimately this variability, coupled with the lack of a "gold standard," highlights the ongoing importance of further validation studies before CTVI can be widely translated from academic centers to the clinic. It is hoped that the information gleaned from the VAMPIRE Challenge can help inform future validation efforts., (© 2018 American Association of Physicists in Medicine.)

Sputum deposition blocks the airways of patients and leads to blood oxygen desaturation. Medical staff must periodically check the breathing state of intubated patients. This process increases staff workload. In this paper, we describe a system designed to acquire respiratory sounds from intubated subjects, extract the audio features, and classify these sounds to detect the presence of sputum. Our method uses 13 features extracted from the time-frequency spectrum of the respiratory sounds. To test our system, 220 respiratory sound samples were collected. Half of the samples were collected from patients with sputum present, and the remainder were collected from patients with no sputum present. Testing was performed based on ten-fold cross-validation. In the ten-fold cross-validation experiment, the logistic classifier identified breath sounds with sputum present with a sensitivity of 93.36% and a specificity of 93.36%. The feature extraction and classification methods are useful and reliable for sputum detection. This approach differs from waveform research and can provide a better visualization of sputum conditions. The proposed system can be used in the ICU to inform medical staff when sputum is present in a patient's trachea.

Pulmonary fissure detection in computed tomography (CT) is a critical component for automatic lobar segmentation. The majority of fissure detection methods use feature descriptors that are hand-crafted, low-level, and have local spatial extent. The design of such feature detectors is typically targeted toward normal fissure anatomy, yielding low sensitivity to weak, and abnormal fissures that are common in clinical data sets. Furthermore, local features commonly suffer from low specificity, as the complex textures in the lung can be indistinguishable from the fissure when the global context is not considered. We propose a supervised discriminative learning framework for simultaneous feature extraction and classification. The proposed framework, called FissureNet, is a coarse-to-fine cascade of two convolutional neural networks. The coarse-to-fine strategy alleviates the challenges associated with training a network to segment a thin structure that represents a small fraction of the image voxels. FissureNet was evaluated on a cohort of 3706 subjects with inspiration and expiration 3DCT scans from the COPDGene clinical trial and a cohort of 20 subjects with 4DCT scans from a lung cancer clinical trial. On both data sets, FissureNet showed superior performance compared with a deep learning approach using the U-Net architecture and a Hessian-based fissure detection method in terms of area under the precision-recall curve (PR-AUC). The overall PR-AUC for FissureNet, U-Net, and Hessian on the COPDGene (lung cancer) data set was 0.980 (0.966), 0.963 (0.937), and 0.158 (0.182), respectively. On a subset of 30 COPDGene scans, FissureNet was compared with a recently proposed advanced fissure detection method called derivative of sticks (DoS) and showed superior performance with a PR-AUC of 0.991 compared with 0.668 for DoS.

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.

Purpose: Regional ventilation and its response to radiation dose can be estimated using four-dimensional computed tomography (4DCT) and image registration. This study investigated the impact of radiation therapy (RT) on ventilation and the dependence of radiation-induced ventilation change on pre-RT ventilation derived from 4DCT., Methods and Materials: Three 4DCT scans were acquired from each of 12 subjects: two scans before RT and one scan 3 months after RT. The 4DCT datasets were used to generate the pre-RT and post-RT ventilation maps by registering the inhale phase image to the exhale phase image and computing the Jacobian determinant of the resulting transformation. The ventilation change between pre-RT and post-RT was calculated by taking a ratio of the post-RT Jacobian map to the pre-RT Jacobian map. The voxel-wise ventilation change between pre- and post-RT was investigated as a function of dose and pre-RT ventilation., Results: Lung regions receiving over 20 Gy exhibited a significant decrease in function (3.3%, P < 0.01) compared to those receiving less than 20 Gy. When the voxels were stratified into high and low pre-RT function by thresholding the Jacobian map at 10% volume expansion (Jacobian = 1.1), high-function voxels exhibited 4.8% reduction in function for voxels receiving over 20 Gy, a significantly greater decline (P = 0.037) than the 2.4% reduction in function for low-function voxels. Ventilation decreased linearly with dose in both high-function and low-function regions. High-function regions showed a significantly larger decline in ventilation (P ≪ 0.001) as dose increased (1.4% ventilation reduction/10 Gy) compared to low-function regions (0.3% ventilation reduction/10 Gy). With further stratification of pre-RT ventilation, voxels exhibited increasing dose-dependent ventilation reduction with increasing pre-RT ventilation, with the largest pre-RT Jacobian bin (pre-RT Jacobian between 1.5 and 1.6) exhibiting a ventilation reduction of 4.8% per 10 Gy., Conclusions: Significant ventilation reductions were measured after radiation therapy treatments, and were dependent on the dose delivered to the tissue and the pre-RT ventilation of the tissue. For a fixed radiation dose, lung tissue with high pre-RT ventilation experienced larger decreases in post-RT ventilation than lung tissue with low pre-RT ventilation., (© 2018 The Authors Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine.)

Lung imaging is increasingly being used to diagnose, quantify, and phenotype chronic obstructive pulmonary disease (COPD). Although spirometry is the gold standard for the diagnosis of COPD and for severity staging, the role of computed tomography (CT) imaging has expanded in both clinical practice and research. COPD is a heterogeneous disease with considerable variability in clinical features, radiographic disease, progression, and outcomes. Recent studies have examined the utility of CT imaging in enhancing diagnostic certainty, improving phenotyping, predicting disease progression and prognostication, selecting patients for intervention, and also in furthering our understanding of the complex pathophysiology of this disease. Multiple CT metrics show promise for use as imaging biomarkers in COPD.

Rationale: A substantial proportion of subjects without overt airflow obstruction have significant respiratory morbidity and structural abnormalities as visualized by computed tomography. Whether regions of the lung that appear normal using traditional computed tomography criteria have mild disease is not known., Objectives: To identify subthreshold structural disease in normal-appearing lung regions in smokers., Methods: We analyzed 8,034 subjects with complete inspiratory and expiratory computed tomographic data participating in the COPDGene Study, including 103 lifetime nonsmokers. The ratio of the mean lung density at end expiration (E) to end inspiration (I) was calculated in lung regions with normal density (ND) by traditional thresholds for mild emphysema (-910 Hounsfield units) and gas trapping (-856 Hounsfield units) to derive the ND-E/I ratio. Multivariable regression analysis was used to measure the associations between ND-E/I, lung function, and respiratory morbidity., Measurements and Main Results: The ND-E/I ratio was greater in smokers than in nonsmokers, and it progressively increased from mild to severe chronic obstructive pulmonary disease severity. A proportion of 26.3% of smokers without airflow obstruction had ND-E/I greater than the 90th percentile of normal. ND-E/I was independently associated with FEV 1 (adjusted β = -0.020; 95% confidence interval [CI], -0.032 to -0.007; P = 0.001), St. George's Respiratory Questionnaire scores (adjusted β = 0.952; 95% CI, 0.529 to 1.374; P < 0.001), 6-minute-walk distance (adjusted β = -10.412; 95% CI, -12.267 to -8.556; P < 0.001), and body mass index, airflow obstruction, dyspnea, and exercise capacity index (adjusted β = 0.169; 95% CI, 0.148 to 0.190; P < 0.001), and also with FEV 1 change at follow-up (adjusted β = -3.013; 95% CI, -4.478 to -1.548; P = 0.001)., Conclusions: Subthreshold gas trapping representing mild small airway disease is prevalent in normal-appearing lung regions in smokers without airflow obstruction, and it is associated with respiratory morbidity. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Rationale: The rate of decline of lung function is greater than age-related change in a substantial proportion of patients with chronic obstructive pulmonary disease, even after smoking cessation. Regions of the lung adjacent to emphysematous areas are subject to abnormal stretch during respiration, and this biomechanical stress likely influences emphysema initiation and progression., Objectives: To assess whether quantifying this penumbra of lung at risk would predict FEV 1 decline., Methods: We analyzed paired inspiratory-expiratory computed tomography images at baseline of 680 subjects participating in a large multicenter study (COPDGene) over approximately 5 years. By matching inspiratory and expiratory images voxel by voxel using image registration, we calculated the Jacobian determinant, a measure of local lung expansion and contraction with respiration. We measured the distance between each normal voxel to the nearest emphysematous voxel, and quantified the percentage of normal voxels within each millimeter distance from emphysematous voxels as mechanically affected lung (MAL). Multivariable regression analyses were performed to assess the relationship between the Jacobian determinant, MAL, and FEV 1 decline., Measurements and Main Results: The mean (SD) rate of decline in FEV 1 was 39.0 (58.6) ml/yr. There was a progressive decrease in the mean Jacobian determinant of both emphysematous and normal voxels with increasing disease stage (P < 0.001). On multivariable analyses, the mean Jacobian determinant of normal voxels within 2 mm of emphysematous voxels (MAL 2 ) was significantly associated with FEV 1 decline. In mild-moderate disease, for participants at or above the median MAL 2 (threshold, 36.9%), the mean decline in FEV 1 was 56.4 (68.0) ml/yr versus 43.2 (59.9) ml/yr for those below the median (P = 0.044)., Conclusions: Areas of normal-appearing lung are mechanically influenced by emphysematous areas and this lung at risk is associated with lung function decline. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).