Cho, Michael H.

Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.

Abstract Rationale Idiopathic Pulmonary Fibrosis (IPF) is characterized by chronic progressive pulmonary fibrosis and high mortality. Genetic markers, summarized into a polygenic risk score (PRS), associate with IPF in well-phenotyped research cohorts. Objectives To evaluate the performance of the PRS using real-world data from routinely captured electronic healthcare records. Methods We conducted an observational study evaluating the association of a PRS for IPF with electronic healthcare record diagnosis of IPF as well as lung transplant-free survival in four independent cohorts; the Mass General Brigham Biobank (MGBB), Mayo Clinic Biobank (MCBB), Mayo Clinic Tapestry Cohort (Tapestry), and U.K. Biobank (UKBB). We used multivariable logistic regression and multivariable Cox proportional hazards models adjusting for age, gender, and principal components of ancestry. The cohorts then underwent fixed and random effects meta-analysis. Measurements and Main Results Of 37,709; 44,195; 43,202; and 447,422 participants from MGBB, MCBB, Tapestry and UKBB respectively, 1,015 (2.7%) 2,879 (6.5%), 1,310 (3.0%), and 2,742 (0.6%) participants had an IPF diagnosis. Meta-analysis demonstrated a high-risk PRS associated with IPF diagnosis, OR 2.88 (95%CI 2.41-3.44) compared to all other individuals. A high-risk PRS also associated with the composite endpoint of mortality or lung transplant among those with an IPF diagnosis, HR 1.23(95%CI 1.11-1.35) compared to all other individuals. Conclusions A PRS can identify those at risk for an IPF diagnosis and mortality in biobank-scale data, which may have implications for clinical decisions. Further work is necessary to evaluate the utility of adding genetics in clinical settings.

Abstract Rationale Interstitial lung abnormalities (ILA) are abnormal findings on chest CT scans that may represent early-stages of interstitial lung disease. The presence of fibrosis (traction bronchiectasis and/or honeycombing) identifies those with ILA at an increased risk for mortality, however, determining who will progress remains a challenge, and the role of blood-based biomarkers in ILA risk-stratification is unknown. Methods Data from four prospective cohort studies were used to identify associations between plasma proteins measured with DNA aptamers (SomaScan) and ILA, fibrosis, and mortality. The Genetic Epidemiology of COPD Study (COPDGene) was used to identify protein associations, and findings were replicated by meta-analysis of the Age, Gene/Environment Susceptibility-Reykjavik study (AGES-Reykjavik), SubPopulations and InteRmediate Outcome Measures in COPD Study (SPIROMICS), and Framingham Heart Study (FHS). Protein associations with ILA and with fibrotic ILA compared to non-fibrotic ILA were performed using multivariable logistic regression; among those with ILA, Cox proportional hazards models were used to identify protein associations with mortality. Models were adjusted for age, gender, race, BMI, smoking status, pack-year smoking history, and for mortality, COPD. Benjamini-Hochberg false discovery rate p-values <0.05 were considered significant. Results 524 (10.7%) of 4,899 COPDGene participants had ILA, of whom 57 (10.9%) had fibrosis; of 8,146 AGES-Reykjavik, SPIROMICS and FHS participants included in the meta-analysis, 666 (8.2%) had ILA, of whom 198 (29.7%) had fibrosis. In COPDGene, 230 proteins were associated with ILA; the most significant of which were SFTPD (adjusted odds ratio per standard deviation [SD] of protein [OR]=2.21, p = 3.7E-41) and SFTPB (OR = 1.82, p = 7.9E-37). In the meta-analysis, 159 of the ILA-associated proteins were replicated. Of the 230 ILA-associated proteins in COPDGene, 9 were associated with fibrosis, and 5 of these 9 proteins were associated with mortality among those with ILA. The 5 proteins associated with all three outcomes were WFDC2 (ORILA=1.90, p = 2.6E-34; ORFibrosis=2.00, p = 0.0001; adjusted hazards of mortality per SD of protein [HR]=1.59, p = 3.4E-05), GDF15 (ORILA=1.67, p = 3.6E-23; ORFibrosis=1.78, p = 0.0006; HR = 1.63, p = 1.4E-06), SCGB3A1 (ORILA=1.52, p = 4.5e-14; ORFibrosis=1.80, p = 0.0002; HR = 1.25, p = 0.02), TREM1 (ORILA=1.44, p = 1.4e-14; ORFibrosis=1.62, p = 0.001; HR = 1.25, p = 0.02) ICAM5 (ORILA=1.42, p = 1.2e-11; ORFibrosis=1.79, p = 0.0001; HR = 1.31, p = 0.002). In the meta-analysis, 6 of the 9 fibrosis-associated proteins and all 5 ILA mortality-associated proteins replicated. Conclusion This represents the first plasma proteomic screen demonstrating reproducible biomarker associations not only with ILA presence, but also with measures of severity and adverse outcomes. These results provide compelling targets for further study, including their prospective evaluation as clinical biomarkers for early disease identification and prognostication. This abstract is funded by: Multiple Sources

Abstract Rationale Deficiency of alpha-1 antitrypsin (AAT), the protein encoded by SERPINA1, remains the strongest genetic risk factor for COPD. Individuals heterozygote for the type Z variant, termed the PI*MZ genotype (estimated global prevalence n = 35,000,000), experience varying respiratory morbidity beyond expected from tobacco smoking alone. AAT is a pleiotropic protein which inhibits not only neutrophil elastase but also neutrophil chemotaxis mediated by leukotriene B4 (LTB4) [by antagonizing the LTB4 receptor (LTB4R)] and interleukin-8 (IL-8). Although LTB-4 and IL-8 gene variants and protein levels have independently been implicated in COPD, their potential role in PI*MZ is unknown. Methods COPDGene and SPIROMICS are observational tobacco-smoke and COPD-enriched cohort studies. Participants were phenotyped by spirometry and CT imaging, with whole genome sequencing data available via TOPMed. We analyzed participants with PI*MM (COPDGene n = 9,732; SPIROMICS n = 2,105) and PI*MZ (COPDGene n = 300; SPIROMICS n = 83). LTB-4 (ALOX5AP, LTA4H, LTB4R) and IL-8 (CXCL8, CXCR1, CXCR2) associated gene variants (exonic and intronic) with a minor allele frequency of ≥ 0.05 were extracted. We assessed statistical interactions between MZ status and chest CT outcomes or FEV1% predicted using the SPAGxECCT package with adjustment for relatedness (genetic related matrix) and ancestry (top ten PCs). Baseline protein LTB4R (SPIROMICS only) and IL-8 (both cohorts) protein levels were measured via SomaScan, and statistical interactions between MZ status and AAT concentration were explored. Models were adjusted for age, sex, BMI, smoking status, and pack-years. Given the candidate gene discovery approach, significance was defined as nominal p-value<0.05 Results We identified the following number of variants in genes of the LTB-4 or IL-8 pathways: 87 shared between cohorts; 133 unique to COPDGene; 1 unique to SPIROMICS. LTB-4-related polymorphisms were associated (p < 0.05) with any radiographic outcomes in 22 variants (COPDGene) vs. 23 variants (SPIROMICS), with a single overlap (rs17074984 in ALOX5AP). Two LTB-4-pathway variants were associated with FEV1% predicted in SPIROMICS (both in ALOX5AP), but none in COPDGene. AAT plasma concentrations displayed a statistically significant interaction with FEV1% predicted (p = 0.046) and PRMemphysema (p = 0.013) for LTBR4, but not for IL-8. There was not an interaction by MM/MZ status and protein levels. Conclusions Predominantly LTB-4, and not IL-8, pathway-related genetic variants and protein levels may contribute to disease heterogeneity among PI*MZ individuals. The association of circulating LTB4R protein with reduced lung function, particularly among those with low AAT levels, further suggest interactions between AAT and the leukotriene pathway in PI*MZ individuals. Thus, PI*MZ individuals may uniquely derive benefit from leukotriene-targeting pharmacotherapy. This abstract is funded by: NIH NHLBI K23HL173570

Abstract Rationale Alpha-1 antitrypsin (AAT) deficiency, primarily caused by homozygous SERPINA1 Z mutation, increases the risk of COPD. Individuals carrying one Z and one normal M allele (MZ) have higher COPD risk than MM individuals, although many never develop disease. COPD is a polygenic condition influenced by numerous genetic factors across the genome and by environmental exposures such as smoking. We aimed to examine how polygenic factors and smoking modify risk of airflow obstruction in MZ compared with MM. Methods We included participants with the MZ and MM genotypes from the Framingham Heart Study. We first tested the association between AAT genotype (MZ coded as 1 and MM as 0) and FEV1/FVC, stratified by ever-smoking status, and further evaluated interactions between AAT genotype and smoking pack-years. Next, for each participant we calculated a polygenic risk score (PRS) for FEV1/FVC, which was derived from external genome-wide association studies in the UK Biobank and SpiroMeta Consortium. To reflect risk for airflow limitation, higher PRSs correspond to lower FEV1/FVC. We then tested the interaction between the PRS and AAT genotype, stratified by ever-smoking. Models were adjusted for age, sex, height, current smoking, pack-years (among ever-smokers), familial relatedness, and in PRS models, principal components of genetic ancestry. Results Among 4,220 participants (mean age 53.7 years, 54.2% female, 50.2% ever-smokers, mean FEV1/FVC 74.6%), 133 (3.2%) had MZ genotype. In never-smokers, MZ individuals had 1.72% higher FEV1/FVC compared with MM (P = 0.032), whereas no difference was observed among ever-smokers (β= -0.39%; P = 0.65). Cumulative smoking modified the effect of MZ genotype on FEV1/FVC (interaction β= -0.13% per pack-year; P = 0.003); this interaction persisted among ever-smokers (β= -0.11% per pack-year; P = 0.044). A one-standard-deviation increase in the FEV1/FVC PRS was associated with lower FEV1/FVC in both never- (β= -2.37%; P < 0.0001) and ever-smokers (β= -2.46%; P < 0.0001). A higher PRS was disproportionately associated with reduced FEV1/FVC in MZ compared to MM (interaction β= -1.87% per standard-deviation PRS increase, P = 0.027) among ever-smokers, but not in never-smokers (P = 0.15). Conclusions In this community-based cohort, the AAT MZ genotype was associated with higher FEV1/FVC among non-smokers compared to MM. The interaction between MZ genotype and pack-years suggested a dose-dependent modifying effect of smoking on the MZ-related COPD risk. Among ever-smokers, higher polygenic risk further amplified susceptibility to airflow obstruction in MZ compared with MM. These findings highlight the roles of both polygenic and environmental risk factors in COPD susceptibility in MZ. Replication of the results is needed. This abstract is funded by: This work is supported by an Alpha-1 Foundation Grant and a Boston University Department of Medicine Career Investment Award

ABSTRACT Background Chronic obstructive pulmonary disease (COPD) is associated with musculoskeletal comorbidities, including cachexia. Weight loss (WL) is the major criterion for cachexia and increases risk for mortality in COPD. Risk factors for WL in COPD are incompletely understood. We performed this whole genome sequencing (WGS) analysis to identify genetic risk variants for WL in COPD. Methods We studied 16 972 participants from the Trans‐Omics for Precision Medicine (TOPMed) Initiative and All of Us Research Program. COPD was diagnosed using spirometry in TOPMed, while diagnosis codes were used in All of Us. WL was defined as WL ≥ 5% or a final body mass index (BMI) < 20 kg/m 2 . WGS data came from white blood cells in all cohorts. Single‐variant testing was conducted on both race‐ and study‐stratified cohorts and in a cosmopolitan, ancestry‐independent manner using GENESIS in TOPMed and SAIGE in All of Us. SAIGE‐GENE+ gene‐based analyses were performed on race‐stratified and cosmopolitan cohorts. Single variant meta‐analyses were conducted using METAL within (B/AA and NHW analyses of Black/African–American and non‐Hispanic white participants, respectively) and across racial groups (COSMO). Rare variant gene‐based results were combined using Fisher's method. Transcriptomic effects were predicted using MetaXcan. We used the GWAS Catalogue to analyse for colocalization with other related traits. Results Two single variants were associated with WL in COPD among All of Us participants: one intronic variant in HCN1 in Black/African–American participants (chr5:45271359:TACACAC:T, odds ratio with 95% confidence interval (OR (CI 95 )) = 2.43(1.78–3.31), p = 1.95 × 10 −8 ) and one intergenic variant between PPP4R2 and PDZRN3 in the cosmopolitan and NHW cohorts (chr3:73345901:A:G, OR (CI 95 ) = 0.21(0.12–0.35), p = 8.84 × 10 −9 in cosmopolitan and OR (CI 95 ) = 0.18(0.10–0.33), p = 9.44 × 10 −9 in NHW). Single‐variant meta‐analysis identified two loci associated with WL in COPD: five variants within DRAIC in the B/AA meta‐analysis (lead variant chr15:69571341:A:G, OR (CI 95 ) = 1.37(1.23–1.51), p = 1.29 × 10 −9 ) and two intronic variants within RFX3 in the cosmopolitan meta‐analysis (lead variant chr9:3390983:T:C, OR (CI 95 ) = 1.50(1.31–1.73), p = 1.06 × 10 −8 ). Rare variants within RNU6‐565P (NHW analysis in All of Us; p = 2.83 × 10 −7 ) and LOC339298 (B/AA analysis in All of Us, p = 1.85 × 10 −6 ) were associated with WL in COPD. The RNU‐565P signal remained significant after combination with TOPMed results (p combined = 1.23 × 10 −6 ). MetaXcan predicted differential expression of LNC00959 in visceral adipose tissue ( p = 1.16 × 10 −6 in COSMO analysis). Colocalization analyses identified genomic associations between BMI and variants in or near DRAIC , RFX3 , PDZRN3 , and LINC00959 . Conclusions In the first WGS analysis of WL in COPD, we have identified seven novel loci. Further characterization of these loci will validate our findings and improve our understanding of the molecular pathophysiology of this condition.

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.

ImportanceChronic obstructive pulmonary disease (COPD) is often undiagnosed. Although genetic risk plays a significant role in COPD susceptibility, its utility in guiding spirometry testing and identifying undiagnosed cases is unclear.ObjectiveTo determine whether a COPD polygenic risk score (PRS) enhances the identification of undiagnosed COPD beyond a case-finding questionnaire (eg, the Lung Function Questionnaire) using conventional risk factors and respiratory symptoms.Design, Setting, and ParticipantsThis cross-sectional analysis of participants 35 years or older who reported no history of physician-diagnosed COPD was conducted using data from 2 observational studies: the community-based Framingham Heart Study (FHS) and the COPD-enriched Genetic Epidemiology of COPD (COPDGene) study.ExposuresModified Lung Function Questionnaire (mLFQ) scores and COPD PRS.Main Outcomes and MeasuresThe primary outcome was spirometry-defined moderate to severe COPD (forced expiratory volume in the first second of expiration/forced vital capacity [FEV1/FVC] &amp;lt;0.7 and FEV1 [percent predicted] &amp;lt;80%). The performance of logistic models was assessed using the PRS, mLFQ score, and PRS plus mLFQ score for predicting spirometry-defined COPD.ResultsAmong 3385 FHS participants (median age, 52.0 years; 45.9% male) and 4095 COPDGene participants (median age, 56.8 years; 55.5% male) who reported no history of COPD, 160 (4.7%) FHS and 775 (18.9%) COPDGene participants had spirometry-defined COPD. Adding the PRS to the mLFQ score significantly improved the area under the curve from 0.78 to 0.84 (P &amp;lt; .001) in FHS, 0.69 to 0.72 (P = .04) in COPDGene non-Hispanic African American, and 0.75 to 0.78 (P &amp;lt; .001) in COPDGene non-Hispanic White participants. At a risk threshold for spirometry referral of 10%, the addition of the PRS to the mLFQ score correctly reclassified 13.8% (95% CI, 6.6%-21.0%) of COPD cases in FHS, but not in COPDGene.Conclusions and RelevanceA COPD PRS enhances the identification of undiagnosed COPD beyond a conventional case-finding approach in the general population. Further research is needed to assess its impact on COPD diagnosis and outcomes.