Martinez, F.J.

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 >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<0.70, FEV1<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 Chronic obstructive pulmonary disease (COPD) manifests later in life and is a leading cause of death worldwide. While smoking-related lung pathology is a long-recognized COPD etiotype, emerging evidence suggests that early-life growth conditions also play an important role. A recently-developed index of early-life growth conditions – the difference between attained height and genotype-predicted height (height-GaP) – has been validated against prospectevely ascertained adverse early-life growth conditions. This study used height-GaP to evaluate the association of early-life growth conditions with COPD phenotype and prognosis later in life. Methods Among current and former smokers with spirometry-defined COPD in the SubPopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS), we computed height-GaP as the difference between attained height and genotype-predicted height. Regression models were fit for respiratory mortality (defined by primary diagnostic code), dyspnea prevalence (Medical Research Council dyspnea rating >1), chronic bronchitis prevalence, airflow obstruction severity as quantified by FEV1Q, and functional capacity (six-minute walk distance), adjusting for age, age2, sex, race-ethnicity, genotype-predicted height, cigarette smoking status and pack-years. Sensitivity analyses accounted for competing risks of death (Fine-Gray subdistribution hazard model) and accounted for later-life height loss (restricting sample to participants under 50 years). Results Among 1,258 participants with COPD (mean±SD age: 65.4±8.0 years, 42.4% female, mean±SD pack-years: 53.6±28.2, mean±SD FEV1%: 63.3±23.9%), there were 202 respiratory deaths over a median 7.5 years of follow-up (Figure). Larger height-GaP deficit was associated with higher respiratory mortality (adjusted hazard ratio [aHR]: 1.30 [95%CI: 1.12 to 1.50] per 1-SD increment in height-GaP deficit). This association remained significant after accounting for competing events (p=0.001) and when restricted to participants under 50 years (p=0.032). Larger height-GaP deficit was associated with higher dyspnea prevalence (adjusted prevalence ratio: 1.13 [95%CI: 1.02 to 1.25] per 1-SD increment in height-GaP deficit), lower FEV1Q (adjusted mean difference: -0.37 [95%CI: -0.45 to -0.29] per 1-SD height-GaP deficit), and shorter 6-minute walk (adjusted mean difference: -16.6 m [95%CI: -23.1 m to -10.0 m] per 1-SD increment in height-GaP deficit), but not chronic bronchitis or CAT score ≥10 (p>0.05 for each). Conclusions Among current and former smokers with COPD, a quantitative index of early-life conditions was associated with higher respiratory mortality, in addition to dyspnea, airflow obstruction severity and reduced functional capacity. These results suggest that early-life growth conditions may influence COPD phenotype and prognosis later in life.

Abstract Rationale: We have previously shown reduced total IgG, IgG subclasses, and levels/function of pneumococcal-specific antibodies are associated with increased risk for exacerbations of COPD (ECOPD). The binding of IgG to its receptor, the Fc receptor for IgG (FcγR) on immune cells, has a critical role in determining antibody levels and function. In addition to our research demonstrating that IgG heavy chain gene (IGHG) variants are associated with ECOPD, single nucleotide variants (SNVs) in genes encoding FcγR (FCGR1-3, which encode receptor subtypes FcγRI, FcγRIIa/b/c, and FcγRIIIa/b) may also result in antibody deficiency and increased exacerbation risk. Due to technical challenges in accurately identifying SNVs in FCGR, this region has not been well-characterized. The limited data that exist are from people of European descent, and prior studies have not included African American participants. We therefore sought to investigate the hypothesis that distinct SNVs in genes encoding FcγR are associated with disease severity and history of exacerbations in African American people with COPD. Methods: We developed a novel pipeline to identify SNVs in the coding (CDS) regions of chromosome 1 using short-read whole genome sequence data. The pipeline was benchmarked using paired long-read samples from Genome in a Bottle. SNVs were called within CDS regions and analyzed in a meta-analysis of self-reported African American participants in SPIROMICS and COPDGene (n=3063). Genetic variants were tested for association with COPD severity (GOLD spirometry stage 1-2 vs 3-4) as well as COPD exacerbations, adjusting for age, sex and principal components controlling for population structure. Results: Eight SNVs were significantly associated with frequent ECOPD (≥2 in one year) in the meta-analysis of African American COPDGene and SPIROMICS participants (Table). Seven of these 8 SNVs were found in FCGR3; 4 in the activating FCGR3A subtype, and 3 in FCGR3B. Seven of 11 SNVs associated with COPD severity were in FCGR3 (4 in FCGR3A; 3 in FCGR3B). One SNV (rs371849) was associated with both frequent exacerbations (P=0.002) and COPD severity (P=0.044). Conclusions: We used a cutting-edge variant calling pipeline to identify multiple SNVs in FCGR, primarily involving FcγRIII receptor subtypes, associated with COPD severity and exacerbations in African American participants in SPIROMICS and COPDGene. FcγRIII variants are known to cause impaired antibody levels and function in people of European descent. The structural and functional importance of newly recognized variants in African American cohorts require further study, as they may represent immunogenetic risk factors for antibody deficiency and exacerbations.