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.
Neutrophil-to-Lymphocyte Ratio as a Biomarker in Clinically Stable Chronic Obstructive Pulmonary Disease: SPIROMICS Cohort
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 (>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.
Epigenomic and Transcriptomic Regulation of the MUC5B/MUC5AC Locus in Bronchial Airway Epithelium and Association With Early Small Airways Disease in Spiromics
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.
Phenotypes and Trajectories of Tobacco-Exposed Persons with Preserved Spirometry: Insights from Lung Volumes
Abstract Rationale Among tobacco-exposed persons with preserved spirometry (TEPSs), we previously demonstrated that different lung volume indices—specifically, elevated total lung capacity (TLC) versus elevated ratio of functional residual capacity to TLC (FRC/TLC)—identify different lung disease characteristics in the COPDGene cohort. Objective We sought to determine differential disease characteristics and trajectories associated with lung volume indices among TEPSs in the SPIROMICS cohort. Methods We categorized TEPSs (n = 814) by tertiles (low, intermediate, and high) of TLC or residual volume–to-TLC ratio (RV/TLC) derived from baseline computed tomography images and then examined clinical and spirometric disease trajectories in mutually exclusive categories of participants with high TLC without high RV/TLC ([TLC]high) versus high RV/TLC without high TLC ([RV/TLC]high). We examined differences in computed tomography–measured emphysema (Hounsfield units [HU] ⩽−950; parametric response mapping [PRM] of emphysema), air trapping (HU⩽−856; PRM of functional small airway disease; a disease probability measure for non-emphysematous gas trapping), airway geometry (the mean square root of wall area of a hypothetical airway with 10 mm internal perimeter), respiratory symptoms (on the modified Medical Research Council Dyspnea Scale; COPD Assessment Test [CAT]; St. George’s Respiratory Questionnaire [SGRQ]; and Short Form-12 [SF12]), and outcomes (annualized exacerbation rate) between the two categories at baseline and over follow-up time up to 8.5 years, using regression modeling adjusted for age, sex, height, weight, and smoking status (current vs. former smoker) and burden (pack-years). Results In TEPSs, the pattern of spirometric disease progression differed between participants with [TLC]high and those with [RV/TLC]high: There was increased forced vital capacity with stable forced expiratory volume in 1 second in participants with [TLC]high, versus unchanged forced vital capacity but nominally decreased forced expiratory volume in 1 second in those with [RV/TLC]high. Compared with participants with [TLC]high, TEPSs with [RV/TLC]high had less emphysema (by HU ⩽−950) but more airway disease (by HU ⩽−856; PRM of functional small airway disease; disease probability measure for gas trapping, and mean square root of wall area of a hypothetical airway with 10 mm internal perimeter), more respiratory symptoms (on the modified Medical Research Council Dyspnea Scale, CAT, SGRQ, and SF12), and more severe exacerbations at baseline. Over an average follow-up of 4.1 ± 2.4 years (range = 0.5–8.5 yr), TEPSs with [RV/TLC]high also had a higher likelihood of developing more severe spirometric disease (preserved ratio impaired spirometry or Global Initiative for Chronic Obstructive Lung Disease Classification 2) and worsening of their respiratory symptoms (on the CAT and SGRQ). Although the incidence rates of respiratory exacerbations, hospitalizations, and mortality were not significantly different between the two categories over the follow-up period, TEPSs with [RV/TLC]high were more likely to have been prescribed a respiratory inhaler at their last follow-up visit. Conclusions In these TEPSs from the SPIROMICS cohort, lung volume stratification by TLC versus RV/TLC identifies two pre-COPD phenotypes with distinct respiratory symptoms, radiographic features, and clinical trajectories. The characteristics of these pre-COPD phenotypes match those previously described in the COPDGene cohort using TLC versus FRC/TLC stratification.
Multi-omic Studies of Lower Airway Epithelium in Spiromics Reveals COPD-associated Methylation Sites and Co-methylation Networks Enriched for Lung Function GWAS Loci and Linked to Microbiota
Importance People who smoked cigarettes may experience respiratory symptoms without spirometric airflow obstruction. These individuals are typically excluded from chronic obstructive pulmonary disease (COPD) trials and lack evidence-based therapies. Objective To define the natural history of persons with tobacco exposure and preserved spirometry (TEPS) and symptoms (symptomatic TEPS). Design, Setting, and Participants SPIROMICS II was an extension of SPIROMICS I, a multicenter study of persons aged 40 to 80 years who smoked cigarettes (&gt;20 pack-years) with or without COPD and controls without tobacco exposure or airflow obstruction. Participants were enrolled in SPIROMICS I and II from November 10, 2010, through July 31, 2015, and followed up through July 31, 2021. Exposures Participants in SPIROMICS I underwent spirometry, 6-minute walk distance testing, assessment of respiratory symptoms, and computed tomography of the chest at yearly visits for 3 to 4 years. Participants in SPIROMICS II had 1 additional in-person visit 5 to 7 years after enrollment in SPIROMICS I. Respiratory symptoms were assessed with the COPD Assessment Test (range, 0 to 40; higher scores indicate more severe symptoms). Participants with symptomatic TEPS had normal spirometry (postbronchodilator ratio of forced expiratory volume in the first second [FEV 1 ] to forced vital capacity &gt;0.70) and COPD Assessment Test scores of 10 or greater. Participants with asymptomatic TEPS had normal spirometry and COPD Assessment Test scores of less than 10. Patient-reported respiratory symptoms and exacerbations were assessed every 4 months via phone calls. Main Outcomes and Measures The primary outcome was assessment for accelerated decline in lung function (FEV 1 ) in participants with symptomatic TEPS vs asymptomatic TEPS. Secondary outcomes included development of COPD defined by spirometry, respiratory symptoms, rates of respiratory exacerbations, and progression of computed tomographic–defined airway wall thickening or emphysema. Results Of 1397 study participants, 226 had symptomatic TEPS (mean age, 60.1 [SD, 9.8] years; 134 were women [59%]) and 269 had asymptomatic TEPS (mean age, 63.1 [SD, 9.1] years; 134 were women [50%]). At a median follow-up of 5.76 years, the decline in FEV 1 was −31.3 mL/y for participants with symptomatic TEPS vs −38.8 mL/y for those with asymptomatic TEPS (between-group difference, −7.5 mL/y [95% CI, −16.6 to 1.6 mL/y]). The cumulative incidence of COPD was 33.0% among participants with symptomatic TEPS vs 31.6% among those with asymptomatic TEPS (hazard ratio, 1.05 [95% CI, 0.76 to 1.46]). Participants with symptomatic TEPS had significantly more respiratory exacerbations than those with asymptomatic TEPS (0.23 vs 0.08 exacerbations per person-year, respectively; rate ratio, 2.38 [95% CI, 1.71 to 3.31], P &lt; .001). Conclusions and Relevance Participants with symptomatic TEPS did not have accelerated rates of decline in FEV 1 or increased incidence of COPD vs those with asymptomatic TEPS, but participants with symptomatic TEPS did experience significantly more respiratory exacerbations over a median follow-up of 5.8 years.
Loss of Bone Density with Inhaled Triamcinolone in Lung Health Study II
Abstract Inhaled glucocorticosteroids (ICS) are commonly prescribed for chronic obstructive pulmonary disease. No adverse effect on bone mineral density (BMD) has been proven. In a randomized double-blind, placebo-controlled trial at seven centers in North America, we recruited 412 current smokers or recent quitters with mild to moderate chronic obstructive pulmonary disease. They used inhaled triamcinolone acetonide, 600 mcg, or placebo, twice daily. We measured femoral neck and lumbar spine BMD at baseline and after 1 and 3 years, and serum osteocalcin at baseline, 3 months, 1 year, and 3 years. After 3 years, BMD at the femoral neck decreased 1.78% more with ICS than with placebo (p < 0.001). More participants in the ICS group experienced 6% or more loss of femoral neck BMD (p = 0.002). Lumbar spine BMD increased in the placebo group by 0.98% but decreased by 0.35% in the ICS group (a difference of 1.33%, p = 0.007). Changes in osteocalcin did not correlate with changes in BMD. Fractures, lost height, or osteoporosis diagnoses were not increased among ICS users compared with placebo users. In summary, the use of inhaled triamcinolone acetonide was associated with loss of BMD at the femoral neck and lumbar spine after 3 years of treatment.