Topologic Parametric Response Mapping Identifies Tissue Subtypes Associated with Emphysema Progression

RATIONALE AND OBJECTIVES

Small airways disease (SAD) and emphysema are significant components of chronic obstructive pulmonary disease (COPD), a heterogenous disease where predicting progression is difficult. SAD, a principal cause of airflow obstruction in mild COPD, has been identified as a precursor to emphysema. Parametric Response Mapping (PRM) of chest computed tomography (CT) can help distinguish SAD from emphysema. Specifically, topologic PRM can define local patterns of both diseases to characterize how and in whom COPD progresses. We aimed to determine if distribution of CT-based PRM of functional SAD (fSAD) is associated with emphysema progression.

MATERIALS AND METHODS

We analyzed paired inspiratory-expiratory chest CT scans at baseline and 5-year follow up in 1495 COPDGene subjects using topological analyses of PRM classifications. By spatially aligning temporal scans, we mapped local emphysema at year five to baseline lobar PRM-derived topological readouts. K-means clustering was applied to all observations. Subjects were subtyped based on predominant PRM cluster assignments and assessed using non-parametric statistical tests to determine differences in PRM values, pulmonary function metrics, and clinical measures.

RESULTS

We identified distinct lobar imaging patterns and classified subjects into three radiologic subtypes: emphysema-dominant (ED), fSAD-dominant (FD), and fSAD-transition (FT: transition from healthy lung to fSAD). Relative to year five emphysema, FT showed rapid local emphysema progression (-57.5% ± 1.1) compared to FD (-49.9% ± 0.5) and ED (-33.1% ± 0.4). FT consisted primarily of at-risk subjects (roughly 60%) with normal spirometry.

CONCLUSION

The FT subtype of COPD may allow earlier identification of individuals without spirometrically-defined COPD at-risk for developing emphysema.

Local heterogeneity of normal lung parenchyma and small airways disease are associated with COPD severity and progression

BACKGROUND

Small airways disease (SAD) is a major cause of airflow obstruction in COPD patients and has been identified as a precursor to emphysema. Although the amount of SAD in the lungs can be quantified using our Parametric Response Mapping (PRM) approach, the full breadth of this readout as a measure of emphysema and COPD progression has yet to be explored. We evaluated topological features of PRM-derived normal parenchyma and SAD as surrogates of emphysema and predictors of spirometric decline.

METHODS

PRM metrics of normal lung (PRMNorm) and functional SAD (PRMfSAD) were generated from CT scans collected as part of the COPDGene study (n = 8956). Volume density (V) and Euler-Poincaré Characteristic (χ) image maps, measures of the extent and coalescence of pocket formations (i.e., topologies), respectively, were determined for both PRMNorm and PRMfSAD. Association with COPD severity, emphysema, and spirometric measures were assessed via multivariable regression models. Readouts were evaluated as inputs for predicting FEV1 decline using a machine learning model.

RESULTS

Multivariable cross-sectional analysis of COPD subjects showed that V and χ measures for PRMfSAD and PRMNorm were independently associated with the amount of emphysema. Readouts χfSAD (β of 0.106, p < 0.001) and VfSAD (β of 0.065, p = 0.004) were also independently associated with FEV1% predicted. The machine learning model using PRM topologies as inputs predicted FEV1 decline over five years with an AUC of 0.69.

CONCLUSIONS

We demonstrated that V and χ of fSAD and Norm have independent value when associated with lung function and emphysema. In addition, we demonstrated that these readouts are predictive of spirometric decline when used as inputs in a ML model. Our topological PRM approach using PRMfSAD and PRMNorm may show promise as an early indicator of emphysema onset and COPD progression.

Quantitative accuracy of lung function measurement using parametric response mapping: A virtual imaging study

Parametric response mapping (PRM) is a voxel-based quantitative CT imaging biomarker that measures the severity of chronic obstructive pulmonary disease (COPD) by analyzing both inspiratory and expiratory CT scans. Although PRM-derived measurements have been shown to predict disease severity and phenotyping, their quantitative accuracy is impacted by the variability of scanner settings and patient conditions. The aim of this study was to evaluate the variability of PRM-based measurements due to the changes in the scanner types and configurations. We developed 10 human chest models with emphysema and air-trapping at end-inspiration and end-expiration states. These models were virtually imaged using a scanner-specific CT simulator (DukeSim) to create CT images at different acquisition settings for energy-integrating and photon-counting CT systems. The CT images were used to estimate PRM maps. The quantified measurements were compared with ground truth values to evaluate the deviations in the measurements. Results showed that PRM measurements varied with scanner type and configurations. The emphysema volume was overestimated by 3 ± 9.5 % (mean ± standard deviation) of the lung volume, and the functional small airway disease (fSAD) volume was underestimated by 7.5±19 % of the lung volume. PRM measurements were more accurate and precise when the acquired settings were photon-counting CT, higher dose, smoother kernel, and larger pixel size. This study demonstrates the development and utility of virtual imaging tools for systematic assessment of a quantitative biomarker accuracy.

Deep learning parametric response mapping from inspiratory chest CT scans: a new approach for small airway disease screening

OBJECTIVES

Parametric response mapping (PRM) enables the evaluation of small airway disease (SAD) at the voxel level, but requires both inspiratory and expiratory chest CT scans. We hypothesize that deep learning PRM from inspiratory chest CT scans can effectively evaluate SAD in individuals with normal spirometry.

METHODS

We included 537 participants with normal spirometry, a history of smoking or secondhand smoke exposure, and divided them into training, tuning, and test sets. A cascaded generative adversarial network generated expiratory CT from inspiratory CT, followed by a UNet-like network predicting PRM using real inspiratory CT and generated expiratory CT. The performance of the prediction is evaluated using SSIM, RMSE and dice coefficients. Pearson correlation evaluated the correlation between predicted and ground truth PRM. ROC curves evaluated predicted PRMfSAD (the volume percentage of functional small airway disease, fSAD) performance in stratifying SAD.

RESULTS

Our method can generate expiratory CT of good quality (SSIM 0.86, RMSE 80.13 HU). The predicted PRM dice coefficients for normal lung, emphysema, and fSAD regions are 0.85, 0.63, and 0.51, respectively. The volume percentages of emphysema and fSAD showed good correlation between predicted and ground truth PRM (|r| were 0.97 and 0.64, respectively, p < 0.05). Predicted PRMfSAD showed good SAD stratification performance with ground truth PRMfSAD at thresholds of 15%, 20% and 25% (AUCs were 0.84, 0.78, and 0.84, respectively, p < 0.001).

CONCLUSION

Our deep learning method generates high-quality PRM using inspiratory chest CT and effectively stratifies SAD in individuals with normal spirometry.

Parametric response mapping cut-off values that predict survival of hepatocellular carcinoma patients after TACE

PURPOSE

Parametric response mapping (PRM) is a novel image-analysis technique applicable to assess tumor viability and predict intrahepatic recurrence of hepatocellular carcinoma (HCC) patients treated with transarterial chemoembolization (TACE). However, to date, the prognostic value of PRM for prediction of overall survival in HCC patients undergoing TACE is unclear. The objective of this explorative, single-center study was to identify cut-off values for voxel-specific PRM parameters that predict the post TACE overall survival in HCC patients.

METHODS

PRM was applied to biphasic CT data obtained at baseline and following 3 TACE treatments of 20 patients with HCC tumors ≥ 2 cm. The individual portal venous phases were registered to the arterial phases followed by segmentation of the largest lesion, i.e., the region of interest (ROI). Segmented voxels with their respective arterial and portal venous phase density values were displayed as a scatter plot. Voxel-specific PRM parameters were calculated and compared to patients' survival at 1, 2, and 3 years post treatment to identify the maximal predictive parameters.

RESULTS

The hypervascularized tissue portion of the ROI was found to represent an independent predictor of the post TACE overall survival. For this parameter, cut-off values of 3650, 2057, and 2057 voxels, respectively, were determined to be optimal to predict overall survival at 1, 2, and 3 years after TACE. Using these cut points, patients were correctly classified as having died with a sensitivity of 80, 92, and 86% and as still being alive with a specificity of 60, 75, and 83%, respectively. The prognostic accuracy measured by area under the curve (AUC) values ranged from 0.73 to 0.87.

CONCLUSION

PRM may have prognostic value to predict post TACE overall survival in HCC patients.

Parametric response mapping on chest computed tomography associates with clinical and functional parameters in chronic obstructive pulmonary disease

BACKGROUND

In the search for specific phenotypes of chronic obstructive pulmonary disease (COPD) computed tomography (CT) derived Parametric Response Mapping (PRM) has been introduced. This study evaluates the association between PRM and currently available biomarkers of disease severity in COPD.

METHODS

Smokers with and without COPD were characterized based on questionnaires, pulmonary function tests, body plethysmography, and low-dose chest CT scanning. PRM was used to calculate the amount of emphysema (PRM^Emph) and non-emphysematous air trapping (i.e. functional small airway disease, PRM^fSAD). PRM was first compared with other biomarkers for emphysema (Perc15) and air trapping (E/I-ratio_MLD). Consequently, linear regression models were utilized to study associations of PRM measurements with clinical parameters.

RESULTS

166 participants were included with a mean ± SD age of 50.5 ± 17.7 years. Both PRM^Emph and PRM^fSAD were more strongly correlated with lung function parameters as compared to Perc15 and E/I-ratio_MLD. PRM^Emph and PRM^fSAD were higher in COPD participants than non-COPD participants (14.0% vs. 1.1%, and 31.6% vs. 8.2%, respectively, both p < 0.001) and increased with increasing GOLD stage (all p < 0.001). Multivariate analysis showed that PRM^fSAD was mainly associated with total lung capacity (TLC) (β = -7.90, p < 0.001), alveolar volume (VA) (β = 7.79, p < 0.001), and residual volume (β = 6.78, p < 0.001), whilst PRM^Emph was primarily associated with Kco (β = 8.95, p < 0.001), VA (β = -6.21, p < 0.001), and TLC (β = 6.20, p < 0.001).

CONCLUSIONS

PRM strongly associates with the presence and severity of COPD. PRM therefore appears to be a valuable tool in differentiating COPD phenotypes.