Huang, Yvonne J

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 Chronic obstructive pulmonary disease (COPD) remains a major public health challenge that contributes greatly to mortality and morbidity worldwide. Although it has long been recognized that the epithelium is altered in COPD, there has been little focus on targeting it to modify the disease course. Therefore, mechanisms that disrupt epithelial cell function in patients with COPD are poorly understood. In this study, we sought to determine whether epigenetic reprogramming of the cell-cell adhesion molecule E-cadherin, encoded by the CDH1 gene, disrupts epithelial integrity. By reducing these epigenetic marks, we can restore epithelial integrity and rescue alveolar airspace destruction. We used differentiated normal and COPD-derived primary human airway epithelial cells, genetically manipulated mouse tracheal epithelial cells, and mouse and human precision-cut lung slices to assess the effects of epigenetic reprogramming. We show that the loss of CDH1 in COPD is due to increased DNA methylation site at the CDH1 enhancer D through the downregulation of the ten-eleven translocase methylcytosine dioxygenase (TET) enzyme TET1. Increased DNA methylation at the enhancer D region decreases the enrichment of RNA polymerase II binding. Remarkably, treatment of human precision-cut slices derived from patients with COPD with the DNA demethylation agent 5-aza-2′-deoxycytidine decreased cell damage and reduced air space enlargement in the diseased tissue. Here, we present a novel mechanism that targets epigenetic modifications to reverse the tissue remodeling in human COPD lungs and serves as a proof of concept for developing a disease-modifying target.