Barjaktarević, Igor

Abstract Rationale Mucus plugs occur in chronic obstructive pulmonary disease (COPD), but little is known about their size and airway location or whether these mucus plug features change as emphysema worsens. Methods Computed tomographic lung scans from 31 participants with COPD from the SubPopulations and InteRmediate Outcomes In COPD Study (SPIROMICS) cohort were analyzed to quantify mucus plug size and airway location and extent of emphysema. Radiologist annotations of mucus plugs were incorporated in an image processing pipeline to generate size and location information. Emphysema was quantified as the percentage of voxels less than -950 Hounsfield units. Measurements and Main Results The length distribution of 563 annotated mucus plugs varied from 2 to 50 mm. The distribution was multimodal, and an 11-mm length defined short (“stubby,” ≤11 mm) and long (“stringy,” >11 mm) plug phenotypes, with stubby plugs being most common (64%). Mucus plugs localized predominantly to airway generations 6 to 9, and their prevalence was highest in lower lobes. On a lobar basis, the mucus plug number decreased as emphysema percentage increased, whereas average mucus plug volume tended to increase. Conclusion COPD is characterized by mucus plugs of varying size and shape, and emphysematous lung regions may be differentially susceptible to formation of these plugs. The location of mucus plugs in segmental and subsegmental airways in COPD makes them amenable to treatment with inhaled medications or bronchoscopy.

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.

BACKGROUND: Although chronic obstructive pulmonary disease (COPD) and heart failure with preserved ejection fraction often coexist with overlapping clinical features, they are usually studied separately. The SPIROMICS HF (Subpopulations and Intermediate Outcome Measures in COPD and Heart Failure Study) is testing hypotheses that new computed tomography emphysema subtypes are associated with specific cardiovascular phenotypes (eg, cor pulmonale , cor pulmonale parvus ), common airway branch variants are associated with right heart dysfunction, and symptomatic tobacco-exposed persons with preserved spirometry have signs of increased left ventricular afterload. METHODS: SPIROMICS is a multicenter observational study of COPD with extensive pulmonary phenotyping of participants with ≥20 pack-years smoking and nonsmoking controls. COPD and COPD severity were defined by standard spirometric criteria and symptomatic tobacco-exposed persons with preserved spirometry by ≥20 pack-years, normal spirometry, and COPD Assessment Test score >10. SPIROMICS HF selected all participants in SPIROMICS visit 5 at major sites. Its comprehensive speckle-tracking echocardiography, which included physiological perturbations of leg raise and low-intensity exercise, was harmonized prospectively with the Multi-Ethnic Study of Atherosclerosis Early Heart Failure and HeartSHARE (Combining Omics, Deep Phenotyping, and Electronic Health Records for Heart Failure Subtypes and Treatment Targets) studies. The cardiopulmonary magnetic resonance imaging protocol with gadolinium administration included myocardial fibrosis sequences, pulmonary angiography, time-resolved 3-dimensional cine magnetic resonance imaging (4-dimensional flow) of venous return, and metronome-paced tachypnea to induce dynamic hyperinflation. Coronary artery calcium was assessed on computed tomography scans. The Kansas City Cardiomyopathy Questionnaire was administered. RESULTS: Of the final sample of 753 participants, 57% had COPD (15% mild, 27% moderate, and 15% severe), 18% had symptomatic tobacco-exposed persons with preserved spirometry, 16% were smoking controls, and 8% were nonsmoking controls. Reproducibility of the main measures from speckle-tracking echocardiography (intraclass correlation coefficient, 0.83–0.99), exercise echocardiography (intraclass correlation coefficient, 0.71–0.99) and magnetic resonance imaging (intraclass correlation coefficient, 0.57–0.99) were good-to-excellent, including in severe COPD. CONCLUSIONS: SPIROMICS HF aims to characterize and understand cardiopulmonary interactions in COPD and COPD-related phenotypes to inform targeted treatments for combined cardiopulmonary failure.

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 Introduction Resilient smokers, who maintain respiratory health despite smoking, remain poorly characterized. Airway epithelial cell gene expression in resilient smokers can provide insight into mechanisms and pathways that may protect smokers from developing smoking-related lung disease. The objective of this study is to identify genes and potential pathways associated with lung-related resilience. Methods Resilient smokers were defined based on a previously published multi-dimensional definition which includes the absence of: 1) spirometric airflow limitation, 2) symptoms (cough/sputum, dyspnea), 3) radiographic (HRCT) abnormalities, 4) exacerbations, and 5) an accelerated FEV1 decline. “Non-resilient” smokers had preserved spirometry but showed abnormalities in at least one domain of resilience. We analyzed airway epithelial brush samples from 21 resilient and 34 non-resilient smokers from SPIROMICS. Differential gene expression was performed using DESeq2 package in R, which estimates differential expression using generalized linear models. We adjusted for sequencing batch, age (continuous), race, gender, and current smoking status. While controlling for these covariates we evaluated differential gene expression between resilient and non-resilient smokers, the interaction between resilience and smoking status, and in each sub-category of resilient (absence of symptoms, radiographic abnormalities, exacerbations, and accelerated decline in FEV1). Results We identified 98 genes that were differentially expressed (FDR ≤0.1) when comparing resilient to non-resilient smokers, 20 of those genes were up regulated in resilient smokers (Figure 1). Notable genes included cadherin-12 (CDH12), aspartate beta-hydroxylase (ASPH), and protein kinase AMP-activated alpha 2 (PRKAA2). CDH12 is involved in cell adhesion and epithelial integrity. ASPH influences cell adhesion and signaling. PRKAA2 plays role in cellular homeostasis and stress response, which may be altered in non-resilient smokers, affecting cellular metabolism and inflammation. BPI folding containing family A member (BPIFA1) was significantly associated with smoking status in the interaction analysis, suggesting its role in modulating the impact of smoking on resilience. BPIFA1 has anti-microbial properties, regulates immune responses, and contributes mucosal surface protection. Sub-category analysis identified 17 genes associated with resilience to symptoms and 2 genes related to physiologic resilience (FEV1 decline), though no significant genes were associated with radiographic or exacerbation resilience. Conclusion Differential gene expression of airway epithelial cells identified distinct expression patterns between resilient and non-resilient smokers. Specifically, genes associated with epithelial cell adhesion, signaling, and cellular homeostasis were up regulated in resilient smokers. These findings suggest that resilient smokers may have transcriptomic profiles and potential mechanisms that mitigate the effects of cigarette smoke.

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 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 &gt;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&gt;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. Bronchoscopy can directly assess cells in airway epithelium in participants with COPD, but it is invasive, expensive, difficult to perform at scale, and contraindicated in severe disease. Blood is much more easily available but contains only a limited subset of cell types, many of which are less relevant to lung biology. Recent advances in RNA-sequencing technology have allowed for reliable measurement of transcriptomics in a wider variety of more accessible tissues, including hair follicle cells. We collected and performed bulk RNA-seq on hair follicles as an ancillary study within the NIH-funded SubPopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) and investigated their cell type composition to better understand the opportunities and limitations of this tissue for the study of lung disease. Methods. SPIROMICS recruited individuals with and without COPD aged 40-80 years with ≥ 20 pack-years smoking history. Participant at several sites underwent collection of approximately 10 hair follicles from the head, which were immediately frozen at -80ºC in 2018-21. mRNA in hair follicles was measured by bulk RNA-seq, and cell type deconvolution was performed using GEDIT v1.7 with the provided BlueCodeV1.0.tsv reference matrix. Cell types were collapsed into broader categories by summing the estimated cell type proportions. Results. The cell type deconvolution of the hair follicle bulk RNA-seq from 80 participants with and without COPD estimated that the tissue was largely composed of immune cells (median 36%, including B cells, T cells, neutrophils, macrophages, NK cells, and monocytes) and epithelial cell types (median 31%), with the largest subgroup of epithelial cell types identified as hair follicular keratinocytes (Figure 1). The next largest component was fibroblasts (median 9%), followed by endothelial cells (median 8%). Signals for myocytes (median 5%), melanocytes (median 5%), and chondrocytes (median 5%) were also present. Conclusion. This initial study demonstrated promising signals for a variety of lung-relevant cell types in hair follicles in SPIROMICS, including immune and epithelial cells in addition to fibroblasts and endothelial cells that are not otherwise easily accessible. Hair follicle gene expression profiling shows promise for understanding numerous facets of lung biology and disease that both complements and exceeds that of other noninvasive tissues such as blood.

Abstract Rationale Reliable data about the natural history of lung function decline in alpha-1 antitrypsin (AAT)-deficient Pi*MZ heterozygotes is largely missing. Objectives We hypothesized that, in adults with a tobacco smoking history, lung function deteriorates faster in Pi*MZ compared with the Pi*MM genotype. Methods We identified 1,856 Pi*MM and 79 Pi*MZ participants with ⩾20 pack-years tobacco smoking history from the SPIROMICS (Subpopulations and Intermediate Outcomes Measures in COPD Study) cohort by DNA sequencing and followed them over a median of 4.8 years, comparing radiographic and clinical characteristics between the two groups over time using regression models. Results Adjusted for age, sex, race, smoking pack-year history, and smoking status, Pi*MZ participants had a lower baseline percent predicted forced expiratory volume in 1 second (FEV1) (65.4% vs. 75.1%; difference 95% confidence interval [CI], −15.4% to −3.9%), more radiographic emphysema (&lt;−950 Hounsfield units%, 12.9% vs. 7.8%; difference 95% CI, 2.8% to 7.5%), and nonsignificantly lower lung density. In the longitudinal analysis, the FEV1 annual rate of decline was similar in both groups over the course of the study (−34.5 ml/yr vs. −34.6 ml/yr for Pi*MZ and Pi*MM; difference 95% CI, −16.9 to 17.1 ml/yr). There were no significant differences between Pi*MZ and Pi*MM individuals in the annualized change in lung density, emphysema, patient-reported outcomes, exacerbations, or survival. The proportion with faster FEV1 decline (annual loss ⩾40 ml) was similar in Pi*MZ and Pi*MM groups. In both groups, faster FEV1 decline was associated with more air trapping and small airway disease at baseline. In Pi*MZ only, faster decline was associated with higher blood eosinophil counts (310 vs. 220 cells/μl; difference 95% CI, 30 to 140 cells/μl). In the subgroup analysis limited to a small number of currently smoking participants, no significant differences in longitudinal outcomes were found. Conclusions Despite a lower FEV1 and more emphysema at enrollment, the longitudinal analysis did not demonstrate significantly greater lung function decline or lung density loss in Pi*MZ compared with Pi*MM participants with tobacco smoking history. Limited sample size and duration of longitudinal follow-up constrain generalizability of our findings, thus prohibiting the conclusion that longitudinal trajectories did not differ between these groups. However, our results may suggest that earlier life events could be responsible for more extensive lung disease at enrollment in Pi*MZ compared with Pi*MM tobacco-exposed individuals.