Background: Chronic obstructive pulmonary disease (COPD) exhibits considerable progression heterogeneity. We hypothesized that elastic principal graph analysis (EPGA) would identify distinct clinical phenotypes and their longitudinal relationships.
Methods: Cross-sectional data from 8,972 tobacco-exposed COPDGene participants, with and without COPD, were used to train a model with EPGA, using thirty clinical, physiologic and CT features. Principal component analysis (PCA) was used to reduce data dimensionality to six principal components. An elastic principal tree was fitted to the reduced space. 4,585 participants from COPDGene Phase 2 were used to test longitudinal trajectories. 2,652 participants from SPIROMICS tested external reproducibility.
Results: Our analysis used cross-sectional data to create an elastic principal tree, where the concept of time is represented by distance on the tree. Six clinically distinct tree segments were identified that differed by lung function, symptoms, and CT features: 1) Subclinical (SC); 2) Parenchymal Abnormality (PA); 3) Chronic Bronchitis (CB); 4) Emphysema Male (EM); 5) Emphysema Female (EF); and 6) Severe Airways (SA) disease. Cross-sectional SPIROMICS data confirmed similar groupings. 5-year data from COPDGene mapped longitudinal changes onto the tree. 29% of patients changed segment during follow-up; longitudinal trajectories confirmed a net flow of patients along the tree, from SC towards Emphysema, although alternative trajectories were noted, through airway disease predominant phenotypes, CB and SA.
Conclusion: This novel analytic methodology provides an approach to defining longitudinal phenotypic trajectories using cross sectional data. These insights are clinically relevant and could facilitate precision therapy and future trials to modify disease progression.
Keywords: Airway Disease; COPD; Disease Progression; Emphysema; Phenotypes.