Hastie, A 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