Citations

Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.

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 Rationale High-attenuation area (HAA) is a computed tomography (CT) tool that correlates with lung inflammation and fibrosis. Systemic molecular correlates of HAA (e.g., plasma proteins) may inform biological processes involved in interstitial lung disease. Objectives To identify plasma proteins that associate with HAA and correlate with a higher probability of developing new-onset fibrotic or subpleural interstitial lung abnormalities (ILAs). Methods Plasma protein levels were measured using a semiquantitative aptamer-based platform in MESA (the Multi-Ethnic Study of Atherosclerosis; N = 5,486) and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study; N = 1,781). Linear regression models identified HAA-associated proteins after adjustment for demographic and socioeconomic factors, CT scanner parameters, study center, and batch. Associations of HAA-related proteins with new-onset fibrotic or subpleural ILAs were examined in MESA participants with ILA assessments on full-lung CT 10 years later. Immunohistochemical staining of select proteins was performed in lung tissue from pulmonary fibrosis cases. Measurements and Main Results There were 75 proteins detected that were significantly associated with HAA in MESA and SPIROMICS. Gene Ontology analysis of these proteins identified processes involved in immune cell chemotaxis and cellular growth and apoptosis. Seven proteins were associated with a higher probability of new-onset fibrotic or subpleural ILAs in MESA, and two of these, junctional adhesion molecule-like protein and GTP cyclohydrolase 1 feedback regulatory protein, stained in areas of fibrosis in lung tissue from patients with interstitial lung disease. Conclusions Plasma proteins associated with more HAA are involved in immune and cellular processes and associate with new-onset fibrotic-subpleural ILA.

ImportanceIdentifying factors associated with resilience during the COVID-19 pandemic can inform targeted interventions and resource allocation for groups disproportionately affected by systemic inequities.ObjectiveTo examine factors associated with self-reported resilience during the COVID-19 pandemic in racially and ethnically diverse, community-dwelling US adults.Design, Setting, and ParticipantsThis cross-sectional study was conducted as part of the Collaborative Cohort of Cohorts for COVID-19 Research (C4R) study, which assessed the associations of the pandemic with self-reported resilience of participants from 14 established US prospective cohorts since January 2021. This report includes participants who responded to the self-reported resilience question on C4R questionnaires. Data was initially analyzed from October 2023 to May 2024, with updated analyses performed from August 2024 to April 2025.ExposureRace and ethnicity, behavior factors, health conditions, and social determinants of health measurements accessed before and during the COVID-19 pandemic through cohort visits and C4R questionnaires.Main Outcomes and MeasuresSelf-reported resilience was collected via 1 question (from the Brief Resilience Scale) in C4R questionnaires, “I tend to bounce back quickly after hard times.” Participants who answered agree or strongly agree were classified as resilient, and those who reported neutral, disagree, or strongly disagree were classified as not resilient. Modified Poisson regression was performed to estimate prevalence ratios (PRs) and access multivariable-adjusted associations with resilience.ResultsOf 31 045 participants (18 672 [60%] women; 10 746 [34.6%] aged &amp;lt;65 years), 1185 (3.8%) identified as American Indian, 6728 (21.7%) as Black, 293 (0.9%) as East Asian, 6311 (20.3%) as Hispanic, 565 (1.8%) as South Asian, and 15 961 (51.3%) as White; a total of 23 103 participants (74.4%) self-identified as resilient. Compared with White participants, Black and Hispanic participants had higher prevalence of self-reported resilience (adjusted PR [aPR], 1.04; 95% CI, 1.02-1.06; aPR, 1.08; 95% CI, 1.06-1.11; respectively) and American Indian and East Asian participants had lower prevalence (aPR, 0.90; 95% CI, 0.86-0.94; aPR, 0.76; 95% CI, 0.68-0.84; respectively). Higher education, being married or living as married, higher income, and overweight were also associated with higher prevalence of resilience. Being female, having diabetes, and being unemployed were associated with lower prevalence of self-reported resilience. Compared with participants with public insurance only, participants with private insurance had higher prevalence of resilience (aPR, 1.07; 95% CI, 1.03-1.10). COVID-19 vaccination and infection statuses were not significantly associated with resilience. Modification analyses showed important racial and ethnic differences in how factors such as hypertension, marital status, and insurance status were associated with resilience.Conclusions and RelevanceIn this cross-sectional study of 31 045 adults, self-reported resilience varied by race, ethnicity, and sociodemographic factors. These findings highlight the complex interplay of individual and social factors in shaping the perception of resilience.

Abstract Rationale: Individuals with chronic obstructive pulmonary disease (COPD) have different clinical trajectories, some with minimal symptom burden, others experiencing recurrent exacerbations, the most severe of which require hospitalization. Existing diagnostics, including spirometry and symptom assessments, lack the precision to identify high-risk patients who may benefit from timely preventive therapies. Nuclear cell-free DNA (ncfDNA) is a well-established marker of systemic tissue injury, implicated in amplifying inflammation by triggering type 1 interferon release by activating DNA-sensing pathways. Because COPD is accompanied by inflammation and tissue injury systemically, we hypothesized that ncfDNA is associated with increased risk of severe exacerbations requiring hospitalization. Methods: Design: We utilized baseline plasma samples randomly selected from 672 individuals with COPD from the Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS) cohort, which prospectively collected exacerbation numbers and type. The primary endpoint was exacerbation requiring hospitalization. Secondary endpoints included total number of exacerbations (defined as use of oral steroids, antibiotics, or healthcare utilization) and severe exacerbations (defined as requiring emergency room care or hospitalization). Measure: We measured plasma ncfDNA using digital droplet PCR. Analysis: We assessed the association of ncfDNA with exacerbation outcomes using negative binomial regression or Cox regression, adjusting for age, percent post-bronchodilator forced expiratory volume in one second (FEV1), supplemental oxygen use, and smoking status. Results: In total, 57% (n=383) of COPD participants experienced exacerbations, 29% (n=197) at least one severe, and 22% (n=145) requiring hospitalization, the latter at a median of 622 days from enrollment (IQR=287-997). Higher ncfDNA levels at enrollment correlated with measures of COPD severity, including baseline FEV1 (beta = -0.12, 95% CI -0.18 to -0.049, p &lt;0.001). Higher ncfDNA associated with increased risk of COPD exacerbation requiring hospitalization (incidence rate ratio (IRR) 1.40, 95% CI 1.12-1.75, p =0.003) and severe exacerbations (IRR 1.36, 95% CI 1.11-1.67, p=0.002) and a trend towards total exacerbations (IRR 1.11, 95% CI 0.98-1.26, p= 0.09). The primary outcome, risk of time-to-first hospitalization, correlated significantly with higher ncfDNA with hazard ratio (HR) 1.23, 95% CI 1.01-1.49, p=0.037. Conclusion: These results suggest that ncfDNA may serve as a predictor for severe and hospitalization-associated COPD exacerbations. To better understand the relationship between ncfDNA and COPD pathophysiology, follow-up studies should explore longitudinal trends and cfDNA tissue sources, and their link to COPD exacerbations.

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.

Background: Despite the availability of effective vaccines and a recent decrease in annual deaths, coronavirus disease-2019 (COVID-19) remains a leading cause of death. Serological studies of anti-nucleocapsid (N) antibody response to infection provide insights into host immunobiology of adaptive immune response, which holds promise for identifying high-risk individuals for adverse acute and chronic COVID-19 outcomes. Hypothesis: Among participants previously infected with SARS-CoV-2, traditional vascular disease risk factors are associated with higher anti-N antibodies while vaccination is associated with lower anti-N antibodies. Methods: Among previously infected participants, anti-N antibodies were measured from dried blood spots collected between February 2021-February 2023 among 1,419 Collaborative Cohort of Cohorts for COVID-19 Research (C4R) participants with prior SARS-CoV-2 infection. We measured anti-N IgG antibodies reported as median fluorescence intensity (MFI) units; reactivity (i.e., seropositivity) was defined based on MFI above a validated threshold. Vascular disease risk factors were assessed via pre-pandemic in-person examination, questionnaires, and medical record review. Multivariable generalized linear models regressed anti-N reactivity (modified Poisson) or LN transformed MFI (linear). Results: Among 1,419 participants with prior infection, mean age (standard deviation) was 65.8(12.1) years, 61% were women, and 42.8% were self-reported from a race/ethnicity minority group. Participant reactivity to nucleocapsid peaked at 69% by 4 months post-infection and waned to only 44% ≥12 months after infection. After multivariable adjustment, higher anti-N antibody response was associated with older age, Hispanic (vs. White) or American Indian (vs. White) race/ethnicity, lower income and education, former smoking, and higher anti-spike antibody levels. Asian race (vs. White) and vaccination ( even after infection ) were associated with lower nucleocapsid reactivity; common cardiometabolic co-morbidities were not associated with anti-N reactivity or levels (Figures 1 and 2). Conclusions: Among participants previously infected with SARS-CoV-2, select sociodemographic and behavioral risk factors were associated with higher anti-N antibody levels while vaccination was associated with lower anti-N antibody levels. The observation that vaccination, even after infection, is related to lower anti-N antibody levels merits further investigation.