Hastie, Annette T.

Abstract Rationale Eosinophilic airway inflammation in COPD is associated with more frequent exacerbations,lower lung function,and corticosteroid responsiveness.DNA methylation captures information about genetic factors and environmental exposures, but is under-studied in eosinophilic COPD.We hypothesized that sputum cell methylation patterns reflect genomic regulation of inflammatory endotypes in COPD by DNA sequence variation. Methods We analyzed DNA previously extracted from induced sputum cells for a prior study(PMID38261629) for which sufficient quantity DNA remained for methylation arrays.For this analysis,we focused on n = 132 SPIROMICS participants with ≥20 pack-years smoking history.Samples were stratified into tertiles based on sputum eosinophil proportions,previously determined from sputum cell differential analysis by the centralized sputum analysis core for SPIROMICS(Univ of North Carolina).We assessed proportion of DNA methylation (β values) spanning 936,866 CpG sites using the Illumina Infinium Methylation EPICv2.0 BeadChip Array,which we transformed into M values(meffil R package). To identify differentially methylated regions (DMRs)(FDR<0.05 in minimum of 2 CpG sites per gene region),we contrasted CpG β values between the high-eosinophil and low-eosinophil tertiles(N = 44 each),using the DMRcate package.We tested for enrichment of genes overlapping with significant DMRs using Enrichr (GWAS Catalog 2025 gene set library).To identify genetic variants associated with DNA methylation levels,we performed Methylation Quantitative Trait Loci (meQTL) analysis using TensorQTL.For each CpG site,we tested associations with SNPs within a ± 500 kb window (cis-meQTL).The linear regression model included age,sex,BMI,smoking status,smoking pack-years,study site,inhaled corticosteroid use, the first two genetic principal components, cell-type proportions (macrophages/monocytes, neutrophils, squamous epithelial cells, lymphocytes), and eosinophil tertile . Significant meQTLs were identified at an FDR<0.05 using Benjamini-Hochberg correction. Results A total of 532 eosinophil-associated DMRs (FDR < 0.05) were identified between the high- and low-eosinophil groups,encompassing 446 unique genes.The top enriched pathways related to counts of eosinophils (adjusted p = 0. 00002464), lymphocytes (adjusted p = 0.0001860), and the sum of eosinophils-plus-basophils (adjusted p = 0.0001860).We identified 18,432 significant SNP-CpG/meQTL pairs (FDR < 0.05),suggesting extensive genetic influence on airway DNA methylation patterns.Among the n = 521 eosinophil-associated CpGs that overlapped with meQTLs,rs55687040 (TENT2,a gene involved in mRNA and miRNA processing) is a GWAS locus for FEV1/FVC ratio (PMID: 26634245), suggesting possible genetic control of airway eosinophilic methylation signatures. Conclusions Our findings disclose distinct DNA methylation signatures associated with sputum eosinophil proportion in COPD.A subset of these CpGs are regulated by genetic variants,suggesting an integrated genetic-epigenetic mechanism underlying airway eosinophilic inflammation with the potential to influence COPD risk and severity. This abstract is funded by: SPIROMICS was supported by contracts from the NIH/NHLBI (HHSN268200900013C, HHSN268200900014C, HHSN268200900015C, HHSN268200900016C, HHSN268200900017C, HHSN268200900018C, HHSN268200900019C, HHSN268200900020C), grants from the NIH/NHLBI (U01 HL137880, U24 HL141762, R01 HL182622, and R01 HL144718), and supplemented by contributions made through the Foundation for the NIH and the COPD Foundation from Amgen; AstraZeneca/MedImmune; Bayer; Bellerophon Therapeutics; Boehringer-Ingelheim Pharmaceuticals, Inc.; Chiesi Farmaceutici S.p.A.; Forest Research Institute, Inc.; Genentech; GlaxoSmithKline; Grifols Therapeutics, Inc.; Ikaria, Inc.; MGC Diagnostics; Novartis Pharmaceuticals Corporation; Nycomed GmbH; Polarean; ProterixBio; Regeneron Pharmaceuticals, Inc.; Sanofi; Sunovion; Takeda Pharmaceutical Company; and Theravance Biopharma and Mylan/Viatris.

Abstract Rationale Higher African genetic ancestry has been associated with lower lung function, among individuals from specific racial-ethnic groups. The influence and interactions of genetic variation and ancestry in lung function is poorly understood. This study investigates ancestry-associated loci influencing lung function via admixture mapping and fine-scale allelic association. Methods 512 African American patients with asthma (≥12 years) from the Severe Asthma Research Program and Collaborative Study on the Genetics of Asthma (SARP/CSGA) cohorts were genotyped on the OmniExpress and HumanHap1M arrays or had whole-genome sequencing (WGS). For admixture mapping, global ancestry was inferred using ADMIXTURE and local ancestry estimated using RFMix v1.5.4 with CEU (N = 99) and YRI (N = 108) references across 15,095 local ancestry segments. Peak regions and contiguous segments (p < 10-3) were fine mapped in SARP using local ancestry estimated from WGS with MAF ≥0.03. The Local Ancestry Adjusted Allelic Association (LAAA) model was used to fine map local genetic effect, local ancestry effect, and the interaction between the two, adjusted for age, sex, height, BMI, cohort, and percentage global African ancestry. SNPs with local genetic or interaction p < 10⁻³ were tested for replication in 282 adolescent/adult participants from the Best African Response to Asthma Drugs (BARD) cohort genotyped with the MEGA array imputed to TOPMed (MAF ≥0.03). Results Mean African genetic ancestry was 81.3% in SARP and 80.1% in BARD. In SARP, African genetic ancestry was inversely associated with lung function (FEV₁:β = −7.9, p = 0.005; FVC:β = −9.0, p = 0.004), with stronger effects in females (FEV₁: β = −9.3, p = 0.003; FVC: β = −11.4, p = 0.001) than males (FEV₁:β = −6.4, p = 0.237; FVC:β = −5.7, p = 0.340). Similar total and sex-specific associations were identified in BARD. Admixture mapping in SARP/CSGA identified significant peaks for maximum FVC (p < 1.87 × 10⁻4) on chromosomes 8 and 17 near RP1 and RBFOX3, and near-significant peaks on chromosome 5 near IQGAP2, SNRPCP2, SV2C and for the maximum FEV₁/FVC ratio on chromosome 10 near PCDH15. Fine mapping identified the top SNP on chromosome 8 (rs73587686; p-allele-dose=9.0x10-5; RGS20; q11.23), with 62 additional SNPs overlapping RGS20. SNP rs7825816 was also identified in SARP (β = 2.35, p-allele-dose =1.03 × 10⁻4) and replicated in BARD (β = 1.15, p-allele-dose =0.03). In BARD, LAAA also identified eight additional SNPs (p ≤ 0.05) in the same LD block (chr8:53914976-53918297) overlapping with RGS20, ATP6V1H, and TCEA1. Conclusions We identified multiple chromosomal regions associated with lung function in asthma patients of predominantly African genetic ancestry. Lung function associations were strongest at genetic variation near a key that was replicated in an independent asthma cohort. This abstract is funded by: R01 K08 HL118128, R01 HL142992, R01 NR013700, R01 HL69167, R01 HL161674, HL69116, HL69167, HL69170, HL69174, U10 HL109164, U10 HL109257, U10 HL109146, U10 HL109172, U10 HL109250, U10 HL109250, U10 HL109250, U10 HL109168, U10 HL109152, U10 HL109086-04, U10 HL098115, HL098102, HL098096, HL098075, HL098090, HL098177, HL098098, HL098107, HL098112, HL098103, and K12 HL89992

Introduction Despite the identification of multiple susceptibility loci by genome-wide association studies (GWAS), considerable chronic obstructive pulmonary disease (COPD) heritability remains unexplained. Aim To identify interaction networks of airway epithelial cell DNA methylation in COPD and further explore potential correlations with airway bacterial composition, as potentially collective regulators of biological pathways influencing COPD severity. Methods Using DNA isolated from bronchial airway brushings of 67 ever-smokers (>20 pack-years) from the SubPopulations and InteRmediate Outcomes Measures in COPD Study (SPIROMICS), we assessed proportion of DNA methylation (β) by epigenome-wide association study (EWAS) and examined associations of differentially methylated CpG probes (DMPs) with risk for moderate-to-severe COPD (N=34) versus absent or mild COPD (N=33). We tested co-methylation modules generated by Weighted Correlation Network Analyses (WGCNA) for associations with moderate-to-severe COPD and with bacterial genus-level relative abundances (16S rRNA sequencing). Results EWAS-identified nominally significant DMPs enriched for lung function GWAS loci. Eigengenes in six WGCNA modules were associated with moderate-to-severe COPD (false discovery rate <0.05). Four of those modules were enriched for forced expiratory volume in 1 s/forced vital capacity GWAS loci, and five overlapped with DMPs from EWAS. Overlapping CpG loci in three COPD-associated modules were adjacent to mucin genes; one had 10 genes highly ranked by connectivity with MUC5B, including important pathway genes: B3GNT6 , DGKI and ITGA8 . CpGs in an independent COPD-associated module showed the most correlations with Rothia , with directionality suggestive of negative associations with moderate-severe COPD. Conclusions Bronchial epithelial DNA methylation modules enriched for lung function GWAS loci associate with COPD severity in SPIROMICS. Potential module relationships to bronchial bacterial composition require further validation.

Abstract Rationale Inflammation is central to chronic obstructive pulmonary disease (COPD) pathogenesis but incompletely represented in COPD prognostic models. The neutrophil-to-lymphocyte ratio (NLR) is a readily available inflammatory biomarker. Objectives To explore the associations of NLR with smoking status, clinical features of COPD, and future adverse outcomes. Methods We analyzed NLR calculated from the complete blood count of participants who currently or formerly smoked (n = 2,624) and tobacco-naive control subjects (n = 187) in the SPIROMICS multicenter observational cohort study. We assessed the stability of NLR at 6 weeks and 1 year, the association with select blood biomarkers, and the impact of smoking on NLR and cell counts. We stratified participants by NLR quartiles to compare cross-sectional clinical features at enrollment, prospectively observed exacerbations at 1 year, and mortality during longitudinal follow up. Results Higher NLR quartiles were broadly associated with more severe clinical features of COPD. NLR values were repeatable at 6 weeks (intraclass correlation coefficient, 0.74) and 1 year (intraclass correlation coefficient, 0.62). The impact of smoking on NLR varied with the severity of airflow limitation, mediated by an interaction between smoking, forced expiratory volume in 1 second percent predicted, and neutrophil counts but not lymphocyte counts. The highest NLR quartile (&gt;3.11) was associated with an increased risk of exacerbation over 1 year (adjusted odds ratio, 1.51; 95% confidence interval, 1.18, 1.92) and increased risk of mortality (adjusted hazard ratio, 1.41; 95% confidence interval, 1.20, 1.66) compared with quartiles 1–3. Conclusions Elevated NLR in stable COPD is a widely available biomarker associated with increased risk for exacerbation and death. The impact of cigarette smoking on NLR varies with disease severity.

Abstract Rationale: MUC5B and MUC5AC are neighboring genes on chromosome 11p15 that encode the two predominant mucins secreted in the respiratory airways. Increased airway mucin concentrations and an imbalance between these two mucins drive COPD symptoms and the development of airflow obstruction or COPD, but the epigenomic regulation of this loci has not been characterized in the airways. Methods: We isolated DNA from bronchial airway brushings of 72 participants with &gt;20 pack-years smoking history from the Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS). 36 participants had GOLD-based moderate-to-severe COPD (FEV1/FVC&lt;0.70, FEV1&lt;80%) and 36 were without or with mild COPD. We used the Illumina Infinium Methylation EPIC BeadChip Array to assess proportion of DNA methylation (β) spanning 689,937 CpG sites. We evaluated proportional DNA methylation (β) at 96 CpG loci within a 103kB region using matrixQTL for cis-expression quantitative trait methylation (eQTM: 10kB flanks relative to MUC5AC/MUC5B) integrating bulk RNA-seq data and mQTL models integrating gene sequencing data to identify DNA variation associated with proportion eQTM methylation. Models included age, sex, BMI, clinical site, current smoking vs past smoking, total pack-years, three ancestral components of ancestry, chipset, and position. Results: We identified a top eQTM, cg00524486, marginally associated with MUC5B expression (coefficient estimate=8.75, punadjusted=0.0014, pFDR=0.065, Figure 1B). This eQTM was associated with CT scan parametric response mapping-based functional small airways disease (PRM-fSAD, beta coefficient[β]=3.64, p=0.011). We tested for DNA variation associated with methylation of this top eQTM which identified two methylation quantitative trait loci (mQTL: rs885454 [β=0.044, minor allele frequency[MAF]=0.09, pFDR=0.01], rs868902 [β=0.020, MAF=0.43, pFDR=0.08, Figure 1A) associated with cg00524486 methylation, of which the more frequent mQTL, rs868902, was marginally associated with PRM-fSAD (p=0.07) in the larger SPIROMICS cohort with MUC5B DNA sequencing data. Conclusions: This preliminary multi-omic study of the lower airway epithelium demonstrates that the expression of MUC5B is epigenetically regulated by neighboring CpG probe methylation (eQTM) represented by DNA variation (mQTL). These studies require further confirmation in a larger dataset of airway samples and support a role for epigenetic variation as a biologic contributor and biomarker of early COPD pathogenesis.

Abstract Rationale: The club cell secretory protein (CC16), encoded by the SCGB1A1 gene, has an anti-inflammatory role in airways diseases, including asthma (Li, AJRCCM 2023; Voraphani, AJRCCM 2023). We hypothesize that additional pathway genes regulated or co-expressed with CC16 have cumulative effects on asthma severity outcomes by modulating airway inflammation and maintaining lung homeostasis. Methods: As published previously, four genes were downregulated in CC16 knockout experiments and correlated with CC16 levels: BPIFA1, SFTPD, LTF, and LYZ (Iannuzo, Front. Immunol 2023). We developed a four-gene biomarker score in 94 patients with asthma from the NHLBI-sponsored Severe Asthma Research Program (SARP3) who underwent research bronchoscopy for lower airway epithelial brushings, RNA isolation, and bulk RNA sequencing. We log-transformed the raw values of each transcript; obtained standardized residuals from regression models that included sex, age, and batch effects; and summed the standardized residuals of all four proteins to estimate a standardized sum score (Zhai, JACI Allergy 2024). Logistic regression models tested the four-gene score, three-gene score without LTF, and individual transcripts for associations with exacerbations during 12-month follow-up and zero-inflated negative binomial regression for number of exacerbations. Results: 30.85% percent of 94 patients experienced at least one exacerbation at the end of follow-up. The CC16-associated four-gene biomarker score showed protective effects on asthma exacerbation risk with a 54% percent reduction in exacerbation risk for every one-fold higher score (OR=0.46, 95% CI: 0.27-0.78, p =0.0040, Figure 1). Transcript expression of BPIFA1 (OR=0.87, 95% CI:0.66-1.16, p=0.337), SFTPD (OR=0.33, 95% CI:0.04-2.52, p=0.283), LYZ (OR=0.62, 95%CI:0.28-1.39, p=0.248), and LTF (OR=0.36, 95% CI:0.18-0.74, p=0.0053) all indicated protective directionality with respect to exacerbation risk. However, only LTF demonstrated a statistically significant association. When we excluded LTF, the remaining three-gene score showed a trend toward association with risk (OR=0.64, 95% CI:0.40-1.03, p=0.068) and was significantly associated with number of exacerbations (RR=0.49, 95% CI:0.26-0.95, p=0.034). Compared to LTF alone (OR=0.48, PPV=0.50, NPV=0.71, Sensitivity=0.14, AUC=0.65), the combined four-gene score (OR=0.46, PPV=0.60, NPV= 0.73, Sensitivity=0.21, AUC=0.68) showed improved predictive performance. Conclusions: We identified a CC16-associated transcriptomic predictive score for asthma exacerbations, with LTF showing the strongest effect. Although the other three genes showed individually weak and non-significant associations with exacerbations, the three gene score was significant for number of exacerbations.We demonstrate the potential of this transcriptomic score as a predictive biomarker for asthma severity-related outcomes. Further research investigating these genes and their interaction with other CC16-associated immune pathways could provide insight into therapeutic targets.

Abstract Rationale It is unknown whether air pollution is associated with radiographic features of interstitial lung disease in individuals with chronic obstructive pulmonary disease (COPD). Objectives To determine whether air pollution increases the prevalence of interstitial lung abnormalities (ILA) or percent high-attenuation areas (HAA) on computed tomography (CT) in individuals with a heavy smoking history and COPD. Methods We performed a cross-sectional study of SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study), focused on current or former smokers with COPD. Ten-year exposure to particulate matter ⩽2.5 μm in aerodynamic diameter (PM2.5), nitrogen oxides (NOx), nitrogen dioxide (NO2), and ozone before enrollment CT (completed between 2010 and 2015) were estimated with validated spatiotemporal models at residential addresses. We applied adjusted multivariable modified Poisson regression and linear regression to investigate associations between pollution exposure and relative risk (RR) of ILA or increased percent HAA (between −600 and −250 Hounsfield units), respectively. We assessed for effect modification by MUC5B-promoter polymorphism (variant allele carriers GT or TT vs. GG at rs3705950), smoking status, sex, and percent emphysema. Results Among 1,272 participants with COPD assessed for HAA, 424 were current smokers, and 249 were carriers of the variant MUC5B allele. A total of 519 participants were assessed for ILA. We found no association between pollution exposure and ILA or HAA. Associations between pollutant exposures and risk of ILA were modified by the presence of MUC5B polymorphism (P value interaction term for NOx = 0.04 and PM2.5 = 0.05) and smoking status (P value interaction term for NOx = 0.05; NO2 = 0.01; and ozone = 0.05). With higher exposure to NOx and PM2.5, MUC5B variant carriers had an increased risk of ILA (RR per 26 ppb NOx, 2.41; 95% confidence interval [CI], 0.97–6.0; and RR per 4 μg ⋅ m−3 PM2.5, 1.43; 95% CI, 0.93–2.2, respectively). With higher exposure to NO2, former smokers had an increased risk of ILA (RR per 10 ppb, 1.64; 95% CI, 1.0–2.7). Conclusions Exposure to ambient air pollution was not associated with interstitial features on CT in this population of heavy smokers with COPD. MUC5B modified the association between pollution and ILA, suggesting that gene–environment interactions may influence prevalence of interstitial lung features in COPD.