Comellas, Alejandro P.

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 Thymic stromal lymphopoietin (TSLP) is an epithelial alarmin chiefly implicated in T2 airway inflammation, though it may also participate in broader (non-T2) inflammatory processes. The utility of airway epithelial gene expression of TSLP and its specific receptor subunit, cytokine receptor-like factor 2 (CRLF2), as biomarkers of downstream inflammation and clinical outcomes remains unclear. Methods We analyzed RNA sequencing data from bronchial airway brushing obtained in phase 1 of the SPIROMICS bronchoscopy sub-study (N = 161: 23 never-smokers, 67 ever-smokers with COPD, 71 ever-smokers without COPD). We related airway TSLP and CRLF2 expression to clinical, hematologic/sputum measures, and gene-expression signatures of T2/eosinophilic and neutrophilic inflammation, including to our published T2 signature and to eosinophil and neutrophil signatures derived from single-cell atlas data. Associations used Spearman correlations and linear models adjusted for age, sex, race/ethnicity, COPD status, smoking status, and pack-years. Results CRLF2 expression correlated with gene-expression evidence of both T2/eosinophilic and neutrophilic inflammation: T2 signature (ρ = 0.26, p < 0.001), eosinophil signature (ρ = 0.44, p < 0.001), and neutrophil signature (ρ = 0.39, p < 0.001). CRLF2 was positively associated with log sputum neutrophil count (ρ = 0.24, p = 0.024) and blood neutrophil count (ρ = 0.036, p < 0.001), but not with sputum or blood eosinophil counts. Higher CRLF2 related to worse baseline lung function, particularly among current smokers, (lower post-bronchodilator FEV1: ρ = −0.21, p = 0.008), as well as greater 5-year FEV1 decline across current and former smokers (ρ = −0.22, p = 0.007). Significant findings persisted after covariate adjustment. In contrast, TSLP expression, although strongly correlated with expression of the epithelial alarmin IL33 (ρ = 0.53, p < 0.001), was negatively correlated with CRLF2 (ρ = −0.24, p = 0.002) and with eosinophil and neutrophil gene-expression signatures (ρ = −0.26 and −0.33, both p < 0.001). TSLP expression showed no association with sputum or blood cell counts. Conclusions In bronchial epithelium from SPIROMICS, gene expression of CRLF2, but not of TSLP, tracked with downstream alarmin-relevant inflammatory programs, consistent with a regulatory role in TSLP signaling through TSLPR, as well as with clinically meaningful outcomes. CRLF2 expression was associated with evidence of both T2/eosinophilic and neutrophilic inflammation, and with lower baseline FEV1 and accelerated FEV1 decline, suggesting engagement of a broader, clinically relevant inflammatory milieu. By contrast, epithelial TSLP transcript levels mirrored IL33 gene expression but were not associated with cellular or transcriptomic inflammation. These findings position CRLF2, rather than TSLP expression itself, as a potential airway biomarker of TSLP-axis activity and disease impact with potential implications for patient stratification as well as strategic approaches for TSLP-pathway-targeted therapies. This abstract is funded by: SPIROMICS was supported by contracts from the NIH/NHLBI (HHSN268200900013C, HHSN268200900014C, HHSN268200900015C, HHSN268200900016C, HHSN268200900017C, HHSN268200900018C, HHSN268200900019C, HHSN268200900020C, 75N92024D00012), grants from the NIH/NHLBI (U01HL137880, U24HL141762, R01HL182622, R01HL144718, and R01HL093081), 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.; Bristol Myers Squibb; 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; Theravance Biopharma; Verona; and Mylan/Viatris.

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.

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. 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 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.

Abstract Rationale Ground-glass opacities (GGOs) in the absence of interstitial lung disease are understudied. Objectives To assess the association of GGOs with white blood cells (WBCs) and progression of quantified chest computed tomography emphysema. Methods We analyzed data of participants in the SPIROMICS study (Subpopulations and Intermediate Outcome Measures in COPD Study). Chest radiologists and pulmonologists labeled regions of the lung as GGOs, and the adaptive multiple feature method (AMFM) trained the computer to assign those labels to image voxels and quantify the volume of the lung with GGOs (%GGOAMFM). We used multivariable linear regression, zero-inflated negative binomial, and proportional hazards regression models to assess the association of %GGOAMFM with WBCs, changes in percentage emphysema, and clinical outcomes. Measurements and Main Results Among 2,714 participants, 1,680 had chronic obstructive pulmonary disease (COPD) and 1,034 had normal spirometry. Among participants with COPD, on the basis of multivariable analysis, current smoking and chronic productive cough were associated with higher %GGOAMFM. Higher %GGOAMFM was cross-sectionally associated with higher WBC and neutrophil concentrations. Higher %GGOAMFM per interquartile range at visit 1 (baseline) was associated with an increase in emphysema at 1-year follow-up visit by 11.7% (relative increase; 95% confidence interval, 7.5–16.1%; P &lt; 0.001). We found no association between %GGOAMFM and 1-year FEV1 decline, but %GGOAMFM was associated with exacerbations and all-cause mortality during a median follow-up of 1,544 days (interquartile interval, 1,118–2,059). Among normal spirometry participants, we found similar results, except that %GGOAMFM was associated with progression to COPD at 1-year follow-up. Conclusions Our findings suggest that GGOAMFM is associated with increased systemic inflammation and emphysema progression.