Cho, Michael H.

Abstract Rationale Acute exacerbations (AEs) of chronic obstructive pulmonary disease (COPD) are associated with significant morbidity and mortality. Whether mucus plugs are associated with prospective exacerbations has not been examined extensively. Objectives To examine associations between mucus plugs on chest computed tomography (CT) and future moderate-to-severe AEs in two independent cohorts with spirometrically-confirmed COPD. Methods Mucus plugs were visually identified on baseline chest computed tomography scans from smokers with Global Initiative for Chronic Obstructive Lung Disease grade 2–4 COPD enrolled in two multicenter cohort studies: ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) and COPDGene (Genetic Epidemiology of COPD). Associations between ordinal mucus plug score categories (0, 1–2, and ≥3) and prospectively ascertained AEs, defined as worsening respiratory symptoms requiring systemic steroids and/or antibiotics (moderate to severe) and/or emergency room visit or hospitalization (severe), were assessed using multivariable-adjusted zero-inflated Poisson regression; subjects were exacerbation-free at enrollment. Measurements and Main Results Among 3,250 participants in COPDGene (mean age ± SD, 63.7 ± 8.4 yr; FEV1, 50.6 ± 17.8% predicted; 45.1% female) and 1,716 participants in ECLIPSE (age, 63.3 ± 7.1 yr; FEV1, 48.3 ± 15.8% predicted; 36.2% female), 44.4% and 46.0% had mucus plugs, respectively. The incidence rates of AEs were 61.0 (COPDGene) and 125.7 (ECLIPSE) per 100 person-years. Relative to those without mucus plugs, the presence of 1–2 and ≥3 mucus plugs was associated with increased risk (adjusted rate ratio [aRR], 1.07 [95% CI, 1.05–1.09] and 1.15 [1.1–1.2] in COPDGene; aRR, 1.06 [95% CI, 1.02–1.09] and 1.12 [1.04–1.2] in ECLIPSE, respectively) for prospective moderate to severe AEs. The presence of 1–2 and ≥3 mucus plugs was also associated with increased risk for severe AEs during follow-up (aRR, 1.05 [95% CI, 1.01–1.08] and 1.09 [1.02–1.18] in COPDGene; aRR, 1.17 [95% CI, 1.07–1.27] and 1.37 [1.15–1.62] in ECLIPSE, respectively). Conclusions Computed tomography–detected mucus plugs are associated with an increased risk for future COPD AEs.

Abstract Rationale Dysanapsis refers to a mismatch between airway tree caliber and lung size arising early in life. Dysanapsis assessed by computed tomography (CT) is evident by early adulthood and associated with chronic obstructive pulmonary disease (COPD) risk later in life. Objectives By examining the genetic factors associated with CT-assessed dysanapsis, we aimed to elucidate its molecular underpinnings and physiological significance across the lifespan. Methods We performed a genome-wide association study of CT-assessed dysanapsis in 11,951 adults, including individuals from two population-based and two COPD-enriched studies. We applied colocalization analysis to integrate genome-wide association study and gene expression data from whole blood and lung. Genetic variants associated with dysanapsis were combined into a genetic risk score that was applied to examine association with lung function in children from a population-based birth cohort (n = 1,278) and adults from the UKBiobank (n = 369,157). Measurements and Main Results CT-assessed dysanapsis was associated with genetic variants from 21 independent signals in 19 gene regions, implicating HHIP (hedgehog interacting protein), DSP, and NPNT as potential molecular targets based on colocalization of their expression. A higher dysanapsis genetic risk score was associated with obstructive spirometry among 5-year-old children and among adults in the fifth, sixth, and seventh decades of life. Conclusions CT-assessed dysanapsis is associated with variation in genes previously implicated in lung development, and dysanapsis genetic risk is associated with obstructive lung function from early life through older adulthood. Dysanapsis may represent an endophenotype link between the genetic variations associated with lung function and COPD.

Abstract Rationale Constantly exposed to the external environment and mutagens such as tobacco smoke, human lungs have one of the highest somatic mutation rates among all human organs. However, the relationship of these mutations to lung disease and function is not known. Objectives To identify the prevalence and significance of clonal somatic mutations in chronic lung diseases. Methods We analyzed the clonal somatic mutations from 1,251 samples of normal and diseased noncancerous lung tissue RNA sequencing with paired whole-genome sequencing from the Lung Tissue Research Consortium. We examined the associations of somatic mutations with lung function, disease status, and computationally deconvoluted cell types in two of the most common diseases represented in our dataset, chronic obstructive pulmonary disease (COPD; 29%) and idiopathic pulmonary fibrosis (IPF; 13%). Measurements and Main Results Clonal somatic mutational burden was associated with reduced lung function in both COPD and IPF. We identified an increased prevalence of clonal somatic mutations in individuals with IPF compared with normal control subjects and individuals with COPD independent of age and smoking status. IPF clonal somatic mutations were enriched in disease-related and airway epithelial–expressed genes such as MUC5B in IPF. Patients who were MUC5B risk variant carriers had increased odds of developing somatic mutations of MUC5B that were explained by increased expression of MUC5B. Conclusions Our identification of an increased prevalence of clonal somatic mutation in diseased lung that correlates with airway epithelial gene expression and disease severity highlights for the first time the role of somatic mutational processes in lung disease genetics.

Abstract Rationale In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case–control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10−95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10−87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Importance Airway mucus plugs are common in patients with chronic obstructive pulmonary disease (COPD); however, the association of airway mucus plugging and mortality in patients with COPD is unknown. Objective To determine whether airway mucus plugs identified on chest computed tomography (CT) were associated with increased all-cause mortality. Design, Setting, and Participants Observational retrospective analysis of prospectively collected data of patients with a diagnosis of COPD in the Genetic Epidemiology of COPD cohort. Participants were non-Hispanic Black or White individuals, aged 45 to 80 years, who smoked at least 10 pack-years. Participants were enrolled at 21 centers across the US between November 2007 and April 2011 and were followed up through August 31, 2022. Exposures Mucus plugs that completely occluded airways on chest CT scans, identified in medium- to large-sized airways (ie, approximately 2- to 10-mm lumen diameter) and categorized as affecting 0, 1 to 2, or 3 or more lung segments. Main Outcomes and Measures The primary outcome was all-cause mortality, assessed with proportional hazard regression analysis. Models were adjusted for age, sex, race and ethnicity, body mass index, pack-years smoked, current smoking status, forced expiratory volume in the first second of expiration, and CT measures of emphysema and airway disease. Results Among the 4483 participants with COPD, 4363 were included in the primary analysis (median age, 63 years [IQR, 57-70 years]; 44% were women). A total of 2585 (59.3%), 953 (21.8%), and 825 (18.9%) participants had mucus plugs in 0, 1 to 2, and 3 or more lung segments, respectively. During a median 9.5-year follow-up, 1769 participants (40.6%) died. The mortality rates were 34.0% (95% CI, 32.2%-35.8%), 46.7% (95% CI, 43.5%-49.9%), and 54.1% (95% CI, 50.7%-57.4%) in participants who had mucus plugs in 0, 1 to 2, and 3 or more lung segments, respectively. The presence of mucus plugs in 1 to 2 vs 0 and 3 or more vs 0 lung segments was associated with an adjusted hazard ratio of death of 1.15 (95% CI, 1.02-1.29) and 1.24 (95% CI, 1.10-1.41), respectively. Conclusions and Relevance In participants with COPD, the presence of mucus plugs that obstructed medium- to large-sized airways was associated with higher all-cause mortality compared with patients without mucus plugging on chest CT scans.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Abstract Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to match proteomes to genomes for 2,812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We were able to correctly match 90%-95% of proteomes to their correct genome and for 95%-99% we could match the proteome to the 1% most likely genome. The accuracy of matching in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. When serial proteomes are available, the matching algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1,000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Abstract With aging, somatic mutations in hematopoietic stem and progenitor cells (HSPC) can give rise to clonal hematopoiesis of indeterminate potential (CHIP), a premalignant state associated with diverse age-related diseases. Here we report the largest multi-ancestry genome-wide analysis of CHIP to date (N = 1,018,305), including individuals of African (N = 85,978), Admixed American (N = 191,371), East Asian (N = 13,532), European (N = 694,015), and South Asian (N = 13,193) ancestry. Multi-ancestry meta-analyses identified 72 genome-wide significant loci, including 44 novel associations implicating genes such as AFF1 , ATF7IP , ATP8B4 , BCL2 , CEBPA , CYRIA , DNM2 , ELF1 , NKX2-3 , PIK3CB , PRDM16 , RPN1 , TERC , and TRIM4 . Notably, variants at MECOM and PHF20L1 showed opposite allelic effects between DNMT3A - and non- DNMT3A -driven CHIP, highlighting driver-specific germline influences. Implicated loci converge on pathways regulating telomere maintenance, cell-cycle control, hematopoietic transcription, DNA damage response and immune signaling. These findings support a model in which germline variation both expands the HSPC pool and biases clonal selection in a driver-dependent manner, providing a mechanistic basis for inter-individual heterogeneity in CHIP. Phenome-wide association analyses further linked CHIP to hematologic, neoplastic, and circulatory traits, with enrichment across hematopoietic and non-hematopoietic cell types. Together, this work expands the genomic and phenomic landscape of CHIP and reveals germline–somatic interactions that shape clonal evolution during aging.

Chronic obstructive pulmonary disease (COPD) is a disease with marked metabolic disturbance. Previous studies have shown the association between single metabolites and lung function for COPD, but whether a combination of metabolites could predict phenotype is unknown. We developed metabolomic severity scores using plasma metabolomics from the Metabolon platform from two US cohorts of ever-smokers: the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) (n = 648; training/testing cohort; 72% non-Hispanic, white; average age 63 years) and the COPDGene Study (n = 1120; validation cohort; 92% non-Hispanic, white; average age 67 years). Separate adaptive LASSO (adaLASSO) models were used to model forced expiratory volume at one second (FEV1) and MESA-adjusted lung density using 762 metabolites common between studies. Metabolite coefficients selected by the adaLASSO procedure were used to create a metabolomic severity score (metSS) for each outcome. A total of 132 metabolites were selected to create a metSS for FEV1. The metSS-only models explained 64.8% and 31.7% of the variability in FEV1 in the training and validation cohorts, respectively. For MESA-adjusted lung density, 129 metabolites were selected, and metSS-only models explained 59.0% of the variability in the training cohort and 17.4% in the validation cohort. Regression models including both clinical covariates and the metSS explained more variability than either the clinical covariate or metSS-only models (53.4% vs. 46.4% and 31.6%) in the validation dataset. The metabolomic pathways for arginine biosynthesis; aminoacyl-tRNA biosynthesis; and glycine, serine, and threonine pathway were enriched by adaLASSO metabolites for FEV1. This is the first demonstration of a respiratory metabolomic severity score, which shows how a metSS can add explanation of variance to clinical predictors of FEV1 and MESA-adjusted lung density. The advantage of a comprehensive metSS is that it explains more disease than individual metabolites and can account for substantial collinearity among classes of metabolites. Future studies should be performed to determine whether metSSs are similar in younger, and more racially and ethnically diverse populations as well as whether a metabolomic severity score can predict disease development in individuals who do not yet have COPD.