Pratte, Katherine

BACKGROUND We constructed multi-trait polygenic risk scores (PRSs) predicting chronic obstructive pulmonary disease (COPD) and exacerbations, validated their performance in diverse cohorts, and identified PRS-related proteins for potential therapeutic targeting. METHODS PRSmix + , a multi-trait PRS framework, is used to train a composite PRS (PRS multi ) in COPDGene non-Hispanic White participants ( n = 6,647). Associations of PRS multi with COPD status (GOLD 2–4 vs. GOLD 0 or ICD) and exacerbation frequency were tested in COPDGene African American ( n = 2,466), ECLIPSE ( n = 1,858), Mass General Brigham Biobank ( n = 15,152), and All of Us ( n = 118,566). Protein prediction models were applied to GWAS summary statistics from traits contributing to PRS multi and were validated with proteomic data in COPDGene ( n = 5,173) and UK Biobank ( n = 5,012). RESULTS PRSmix + selected 7 traits for PRS multi . In multivariable models, PRS multi was associated with COPD status (meta-analysis random effects [RE] OR 1.58 [95% CI: 1.28–1.94]) and exacerbation frequency (meta-analysis RE β 0.21 [95% CI: 0.11–0.31]), with higher effect sizes observed in smoking-enriched cohorts. PRS multi outperformed traditional single-trait PRS in all tested cohorts. Using protein prediction models, we identified 73 proteins associated with the PRSs that were also validated with measured protein levels in COPDGene and UK Biobank. Of these proteins, 25 were linked to approved or investigational drugs. Notable targets include RAGE/sRAGE, IL1RL1, and SCARF2, all implicated in COPD pathogenesis and exacerbations. CONCLUSIONS Multi-trait PRS improves prediction of COPD and exacerbation risk. Integration with proteomic data identifies druggable protein targets, offering a promising avenue for precision medicine in COPD management. TRIAL REGISTRATION COPDGene: ClinicalTrials.gov NCT00608764; ECLIPSE: ClinicalTrials.gov NCT00292552.

Introduction There is an increased risk of chronic obstructive pulmonary disease (COPD) in heterozygotes for the alpha-1 antitrypsin (AAT) Z allele (PI*MZ), but there is significant variation in outcomes. Matrix metalloproteinases (MMPs) contribute to worse airway disease and emphysema. Given that the AAT protein is an antiprotease, we hypothesised that MMPs play a modifying role among AAT-deficient individuals. Methods We studied PI*MZ (n=39) and PI*ZZ (n=27) individuals from the Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis Study cohort. MMP (MMP3, MMP9, MMP10, MMP12 and MMP15) gene expression was assessed in bronchoalveolar lavage cells and peripheral blood mononuclear cells. Plasma protein MMP (MMP3, MMP9, MMP10 and MMP12) expression was assessed using SomaScan and BAL MMP15 protein using ELISA. MMP values were log-transformed, and linear regression analyses were used to assess lung function, patient-reported outcomes and CT emphysema adjusted for age, body mass index, sex and duration of smoking. Results Alveolar MMP3 and MMP15 were uniquely associated with worse lung function, patient-reported outcomes and emphysema among PI*MZ individuals. MMP10 was associated with less disease severity among PI*MZ individuals. Discussion Our study shows that alveolar MMP3 and MMP15 may uniquely contribute to disease among PI*MZ individuals, raising the possibility that they could serve as biomarkers or therapeutic targets in this highly prevalent COPD phenotype. Furthermore, we show that MMP10 may serve a protective effect. In aggregate, these data support further studies of the MMP pathway among lung-affected PI*MZ individuals, as they may serve as biomarkers or pharmacological targets in future studies.

Abstract Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to match proteomes to genomes for 2,812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We were able to correctly match 90%-95% of proteomes to their correct genome and for 95%-99% we could match the proteome to the 1% most likely genome. The accuracy of matching in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. When serial proteomes are available, the matching algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1,000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Chronic obstructive pulmonary disease (COPD) is a disease with marked metabolic disturbance. Previous studies have shown the association between single metabolites and lung function for COPD, but whether a combination of metabolites could predict phenotype is unknown. We developed metabolomic severity scores using plasma metabolomics from the Metabolon platform from two US cohorts of ever-smokers: the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) (n = 648; training/testing cohort; 72% non-Hispanic, white; average age 63 years) and the COPDGene Study (n = 1120; validation cohort; 92% non-Hispanic, white; average age 67 years). Separate adaptive LASSO (adaLASSO) models were used to model forced expiratory volume at one second (FEV1) and MESA-adjusted lung density using 762 metabolites common between studies. Metabolite coefficients selected by the adaLASSO procedure were used to create a metabolomic severity score (metSS) for each outcome. A total of 132 metabolites were selected to create a metSS for FEV1. The metSS-only models explained 64.8% and 31.7% of the variability in FEV1 in the training and validation cohorts, respectively. For MESA-adjusted lung density, 129 metabolites were selected, and metSS-only models explained 59.0% of the variability in the training cohort and 17.4% in the validation cohort. Regression models including both clinical covariates and the metSS explained more variability than either the clinical covariate or metSS-only models (53.4% vs. 46.4% and 31.6%) in the validation dataset. The metabolomic pathways for arginine biosynthesis; aminoacyl-tRNA biosynthesis; and glycine, serine, and threonine pathway were enriched by adaLASSO metabolites for FEV1. This is the first demonstration of a respiratory metabolomic severity score, which shows how a metSS can add explanation of variance to clinical predictors of FEV1 and MESA-adjusted lung density. The advantage of a comprehensive metSS is that it explains more disease than individual metabolites and can account for substantial collinearity among classes of metabolites. Future studies should be performed to determine whether metSSs are similar in younger, and more racially and ethnically diverse populations as well as whether a metabolomic severity score can predict disease development in individuals who do not yet have COPD.