Cho, Michael H

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 Idiopathic Pulmonary Fibrosis (IPF) is characterized by chronic progressive pulmonary fibrosis and high mortality. Genetic markers, summarized into a polygenic risk score (PRS), associate with IPF in well-phenotyped research cohorts. Objectives To evaluate the performance of the PRS using real-world data from routinely captured electronic healthcare records. Methods We conducted an observational study evaluating the association of a PRS for IPF with electronic healthcare record diagnosis of IPF as well as lung transplant-free survival in four independent cohorts; the Mass General Brigham Biobank (MGBB), Mayo Clinic Biobank (MCBB), Mayo Clinic Tapestry Cohort (Tapestry), and U.K. Biobank (UKBB). We used multivariable logistic regression and multivariable Cox proportional hazards models adjusting for age, gender, and principal components of ancestry. The cohorts then underwent fixed and random effects meta-analysis. Measurements and Main Results Of 37,709; 44,195; 43,202; and 447,422 participants from MGBB, MCBB, Tapestry and UKBB respectively, 1,015 (2.7%) 2,879 (6.5%), 1,310 (3.0%), and 2,742 (0.6%) participants had an IPF diagnosis. Meta-analysis demonstrated a high-risk PRS associated with IPF diagnosis, OR 2.88 (95%CI 2.41-3.44) compared to all other individuals. A high-risk PRS also associated with the composite endpoint of mortality or lung transplant among those with an IPF diagnosis, HR 1.23(95%CI 1.11-1.35) compared to all other individuals. Conclusions A PRS can identify those at risk for an IPF diagnosis and mortality in biobank-scale data, which may have implications for clinical decisions. Further work is necessary to evaluate the utility of adding genetics in clinical settings.

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.

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.

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program ( n = 63,302) and UK Biobank ( n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Abstract Rationale Multiple studies have demonstrated an increased risk of chronic obstructive pulmonary disease (COPD) in heterozygous carriers of the AAT (alpha-1 antitrypsin) Z allele. However, it is not known if MZ subjects with COPD are phenotypically different from noncarriers (MM genotype) with COPD. Objectives To assess if MZ subjects with COPD have different clinical features compared with MM subjects with COPD. Methods Genotypes of SERPINA1were ascertained by using whole-genome sequencing data in three independent studies. We compared outcomes between MM subjects with COPD and MZ subjects with COPD in each study and combined the results in a meta-analysis. We performed longitudinal and survival analyses to compare outcomes in MM and MZ subjects with COPD over time. Measurements and Main Results We included 290 MZ subjects with COPD and 6,184 MM subjects with COPD across the three studies. MZ subjects had a lower FEV1% predicted and greater quantitative emphysema on chest computed tomography scans compared with MM subjects. In a meta-analysis, the FEV1was 3.9% lower (95% confidence interval [CI], −6.55% to −1.26%) and emphysema (the percentage of lung attenuation areas <−950 HU) was 4.14% greater (95% CI, 1.44% to 6.84%) in MZ subjects. We found one gene, PGF(placental growth factor), to be differentially expressed in lung tissue from one study between MZ subjects and MM subjects. Conclusions Carriers of the AAT Z allele (those who were MZ heterozygous) with COPD had lower lung function and more emphysema than MM subjects with COPD. Taken with the subtle differences in gene expression between the two groups, our findings suggest that MZ subjects represent an endotype of COPD.