Search

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

BackgroundHigh-dimensional clinical data are becoming more accessible in biobank-scale datasets. However, accurately phenotyping high-dimensional clinical data remains a major impediment to genetic discovery.MethodsWe introduce a general deep learning framework, RE presentation learning for Genetic discovery on Low-dimensional Embeddings (REGLE), for discovering associations between genetic variants and high-dimensional clinical data. REGLE uses convolutional variational autoencoders to compute anon-linear, low-dimensional, disentangled embeddingof the data and can also incorporate expert clinical metrics. We demonstrate the utility of REGLE by application to spirograms, which measure lung function. We generate two types of synthetic representations of pulmonary functions we call spirogram encodings (SPINCs) and residual spirogram encodings (RSPINCs).FindingsGenome-wide association studies on (R)SPINCs identify more genome-wide significant loci than existing methods while replicating most known lung function loci. Furthermore, (R)SPINCs are associated with overall survival and, under the latent causal variable model, they exhibit significantly high genetic causality proportion with asthma, chronic obstructive pulmonary disease (COPD), and inflammatory diseases. Finally, we construct a set of polygenic risk scores (PRS) that are generally predictive of pulmonary traits and diseases. We demonstrate superior performance predicting asthma and COPD, in multiple ancestries and across four biobanks, compared to PRSs constructed using expert-defined pulmonary function measurements.InterpretationREGLE is a method for generating low-dimensional, disentangled representations of high-dimensional clinical data that does not require labels, and improves upon expert-defined phenotypes for genetic discovery and disease prediction. It can flexibly incorporate expert-defined or clinical features and provides a framework to create accurate disease-specific PRS in datasets which have minimal expert phenotyping. (R)SPINCs are quantifying clinically relevant features that are not currently captured in a standardized or automated way.FundingGoogle LLC.

BackgroundSpirometry measures lung function by selecting the best of multiple efforts meeting pre-specified quality control (QC), and reporting two key metrics: forced expiratory volume in 1 second (FEV1) and forced vital capacity (FVC). We hypothesize that discarded submaximal and QC-failing data meaningfully contribute to the prediction of airflow obstruction and all-cause mortality.MethodsWe evaluated volume-time spirometry data from the UK Biobank. We identified “best” spirometry efforts as those passing QC with the maximum FVC. “Discarded” efforts were either submaximal or failed QC. To create a combined representation of lung function we implemented a contrastive learning approach,Spirogram-basedContrastiveLearningFramework (Spiro-CLF), which utilized all recorded volume-time curves per participant and applied different transformations (e.g. flow-volume, flow-time). In a held-out 20% testing subset we applied the Spiro-CLF representation of a participant’s overall lung function to 1) binary predictions of FEV1/FVC < 0.7 and FEV1Percent Predicted (FEV1PP) < 80%, indicative of airflow obstruction, and 2) Cox regression for all-cause mortality.FindingsWe included 940,705 volume-time curves from 352,684 UK Biobank participants with 2-3 spirometry efforts per individual (66.7% with 3 efforts) and at least one QC-passing spirometry effort. Of all spirometry efforts, 24.1% failed QC and 37.5% were submaximal. Spiro-CLF prediction of FEV1/FVC < 0.7 utilizing discarded spirometry efforts had an Area under the Receiver Operating Characteristics (AUROC) of 0.981 (0.863 for FEV1PP prediction). Incorporating discarded spirometry efforts in all-cause mortality prediction was associated with a concordance index (c-index) of 0.654, which exceeded the c-indices from FEV1(0.590), FVC (0.559), or FEV1/FVC (0.599) from each participant’s single best effort.InterpretationA contrastive learning model using raw spirometry curves can accurately predict lung function using submaximal and QC-failing efforts. This model also has superior prediction of all-cause mortality compared to standard lung function measurements.FundingMHC is supported by NIH R01HL137927, R01HL135142, HL147148, and HL089856.BDH is supported by NIH K08HL136928, U01 HL089856, and an Alpha-1 Foundation Research Grant.DH is supported by NIH 2T32HL007427-41EKS is supported by NIH R01 HL152728, R01 HL147148, U01 HL089856, R01 HL133135, P01 HL132825, and P01 HL114501.PJC is supported by NIH R01HL124233 and R01HL147326.SPB is supported by NIH R01HL151421 and UH3HL155806.TY, FH, and CYM are employees of Google LLC

IntroductionIn the personalized risk quantification of chronic obstructive pulmonary disease (COPD), genome-wide association studies and polygenic risk scores (PRS) complement traditional risk factors, such as age and cigarette smoking. However, despite being at considerable levels of risk, some individuals do not develop COPD. Research on COPD resilience remains largely unexplored.MethodsWe applied the previously published COPD PRS to whole genome sequencing data from non-Hispanic white and African American individuals in the COPDGene study. We defined genetic resilience as individuals unaffected by COPD with a polygenic risk score above the 90thpercentile. We defined risk-matched case individuals as those with COPD (i.e., FEV1/FVC < 0.70) and a PRS above the 90thpercentile. We defined low risk individuals without COPD (i.e., FEV1/FVC > 0.70) as a polygenic risk score below the 10thpercentile. We compared genetically resilient individuals to risk-matched individuals with COPD and low risk individuals by demographics, lung function, respiratory symptoms, co-morbidities, and chest CT scan measurements. We also performed survival analyses, differential expression analysis, and matching for sensitivity analyses.ResultsWe identified 211 resilient individuals without COPD, 605 genetic risk-matched individuals with COPD, and 527 low-risk individuals without COPD. Resilient individuals had higher FEV1% predicted and lower percent emphysema. In contrast, resilient individuals had higher airway wall thickness compared to low-risk unaffected individuals. While there was no difference in survival between low-risk and resilient individuals, resilient individuals had higher survival compared to risk matched cases. We also identified two genes that were differentially expressed between low-risk unaffected individuals and resilient individuals.ConclusionGenetically resilient individuals had a reduced burden of COPD disease-related measures compared to risk-matched cases but had subtly increased measures compared to low-risk unaffected individuals. Further genetic studies will be needed to illuminate the underlying pathobiology of our observations.

Among human organs, the lung harbors one of the highest rates of somatic mutations. However, the relationship of these mutations to lung disease and function is not known. We analyzed the somatic mutational pattern from 1,251 samples of normal and diseased non-cancerous lung tissue from the Lung Tissue Research Consortium using RNA-seq. In two of the most common diseases represented in our dataset, chronic obstructive pulmonary disease (COPD, 29%) and idiopathic pulmonary fibrosis (IPF, 13%), we found a significantly increased burden of somatic mutations compared to normal. Using deconvoluted cell type proportions, we found that a major predictor of somatic mutations was the airway to alveolar cell proportion and pathogenic cell types. We also found that mutational burden was associated with reduced lung function. This relationship remained even after adjustment for age, sex, smoking, and cell type proportion and in COPD and IPF. Our identification of an increased prevalence of somatic mutation in diseased lung that correlates with cell type proportion and disease severity highlights for the first time the role of somatic mutational processes in lung disease genetics.

Rationale: The ability of peripheral blood biomarkers to assess COPD risk and progression is unknown. Genetics and gene expression may capture important aspects of COPD-related biology that predict...

Chronic Obstructive Pulmonary Disease (COPD) has a simple physiological diagnostic criterion but a wide range of clinical characteristics. The mechanisms underlying this variability in COPD phenotypes are unclear. To investigate the potential contribution of genetic variants to phenotypic heterogeneity, we examined the association of genome-wide associated lung function, COPD, and asthma variants with other phenotypes using phenome-wide association results derived in the UK Biobank. Our clustering analysis of the variants-phenotypes association matrix identified three clusters of genetic variants with different effects on white blood cell counts, height, and body mass index (BMI). To assess the potential clinical and molecular effects of these groups of variants, we investigated the association between cluster-specific genetic risk scores and phenotypes in the COPDGene cohort. We observed differences in steroid use, BMI, lymphocyte counts, chronic bronchitis, and differential gene and protein expression across the three genetic risk scores. Our results suggest that multi-phenotype analysis of obstructive lung disease-related risk variants may identify genetically driven phenotypic patterns in COPD.

RATIONALE: Current guidelines do not sufficiently capture the heterogeneous nature of asthma; a more detailed molecular classification is needed. Metabolomics represents a novel and compelling approach to derive asthma endotypes (i.e., subtypes defined by functional and/or pathobiological mechanisms). OBJECTIVES: To validate metabolomic-driven endotypes of asthma and explore their underlying biology. METHODS: In the Genetics of Asthma in Costa Rica Study (GACRS), untargeted metabolomic profiling, similarity network fusion, and spectral clustering was used to identify metabo-endotypes of asthma, and differences in asthma-relevant phenotypes across these metabo-endotypes were explored. The metabo-endotypes were recapitulated in the Childhood Asthma Management Program (CAMP), and clinical differences were determined. Metabolomic drivers of metabo-endotype membership were investigated by meta-analyzing findings from GACRS and CAMP. MEASUREMENTS AND MAIN RESULTS: Five metabo-endotypes were identified in GACRS with significant differences in asthma-relevant phenotypes, including prebronchodilator (p-ANOVA = 8.3 × 10(−5)) and postbronchodilator (p-ANOVA = 1.8 × 10(−5)) FEV(1)/FVC. These differences were validated in the recapitulated metabo-endotypes in CAMP. Cholesterol esters, trigylcerides, and fatty acids were among the most important drivers of metabo-endotype membership. The findings suggest dysregulation of pulmonary surfactant homeostasis may play a role in asthma severity. CONCLUSIONS: Clinically meaningful endotypes may be derived and validated using metabolomic data. Interrogating the drivers of these metabo-endotypes has the potential to help understand their pathophysiology.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI’s Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

BACKGROUND: Smokers manifest varied phenotypes of pulmonary impairment. RESEARCH QUESTION: Which pulmonary phenotypes are associated with coronary artery disease (CAD) in smokers? STUDY DESIGN AND METHODS: We analyzed data from the University of Pittsburgh COPD Specialized Center for Clinically Oriented Research (SCCOR) cohort (n = 481) and the Genetic Epidemiology of COPD (COPDGene) cohort (n = 2,580). Participants were current and former smokers with > 10 pack-years of tobacco exposure. Data from the two cohorts were analyzed separately because of methodologic differences. Lung hyperinflation was assessed by plethysmography in the SCCOR cohort and by inspiratory and expiratory CT scan lung volumes in the COPDGene cohort. Subclinical CAD was assessed as the coronary artery calcium score, whereas clinical CAD was defined as a self-reported history of CAD or myocardial infarction (MI). Analyses were performed in all smokers and then repeated in those with airflow obstruction (FEV(1) to FVC ratio, < 0.70). RESULTS: Pulmonary phenotypes, including airflow limitation, emphysema, lung hyperinflation, diffusion capacity, and radiographic measures of airway remodeling, showed weak to moderate correlations (r < 0.7) with each other. In multivariate models adjusted for pulmonary phenotypes and CAD risk factors, lung hyperinflation was the only phenotype associated with calcium score, history of clinical CAD, or history of MI (per 0.2 higher expiratory and inspiratory CT scan lung volume; coronary calcium: OR, 1.2; 95% CI, 1.1-1.5; P = .02; clinical CAD: OR, 1.6; 95% CI, 1.1-2.3; P = .01; and MI in COPDGene: OR, 1.7; 95% CI, 1.0-2.8; P = .05). FEV(1) and emphysema were associated with increased risk of CAD (P < .05) in models adjusted for CAD risk factors; however, these associations were attenuated on adjusting for lung hyperinflation. Results were the same in those with airflow obstruction and were present in both cohorts. INTERPRETATION: Lung hyperinflation is associated strongly with clinical and subclinical CAD in smokers, including those with airflow obstruction. After lung hyperinflation was accounted for, FEV(1) and emphysema no longer were associated with CAD. Subsequent studies should consider measuring lung hyperinflation and examining its mechanistic role in CAD in current and former smokers.

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.

Chronic obstructive pulmonary disease (COPD), the third leading cause of death worldwide, is highly heritable. While COPD is clinically defined by applying thresholds to summary measures of lung function, a quantitative liability score has more power to identify new genetic signals. Here we train a deep convolutional neural network on noisy self-reported and ICD-based labels to predict COPD case/control status from high-dimensional raw spirograms and use the model predictions as a liability score. The machine-learning-based (ML-based) liability score accurately discriminates COPD cases and controls (AUROC = 0.82 ± 0.01) and COPD-related hospitalization (AUROC = 0.89 ± 0.01) without any domain-specific knowledge. Moreover, the ML-based liability score is associated with overall survival (Hazard ratio = 1.22 ± 0.01; P ≤ 2 × 10−16) and exacerbation events (R2 = 0.10 ± 0.01; P ≤ 4 × 10−101). A genome-wide association study on the ML-based liability score replicates existing COPD and lung function loci, but also identifies 67 new loci. Thirty-eight of these have supportive evidence in independent datasets, including a locus near LTBR. We demonstrate the biological plausibility of the novel variants through enrichment analyses, phenome-wide association studies, and generalizability of COPD prediction in multiple datasets. These results provide an example of the potential to improve genetic discovery of disease-relevant variants by training deep neural networks to predict noisy labels from high-dimensional raw data.

The serpin family A member 1 (SERPINA1) Z allele is present in approximately one in 25 individuals of European ancestry. Z allele homozygosity (Pi*ZZ) is the most common cause of alpha 1‐antitrypsin deficiency and is a proven risk factor for cirrhosis. We examined whether heterozygous Z allele (Pi*Z) carriers in United Kingdom (UK) Biobank, a population‐based cohort, are at increased risk of liver disease. We replicated findings in Massachusetts General Brigham Biobank, a hospital‐based cohort. We also examined variants associated with liver disease and assessed for gene–gene and gene–environment interactions. In UK Biobank, we identified 1,493 cases of cirrhosis, 12,603 Z allele heterozygotes, and 129 Z allele homozygotes among 312,671 unrelated white British participants. Heterozygous carriage of the Z allele was associated with cirrhosis compared to noncarriage (odds ratio [OR], 1.53; P = 1.1×10−04); homozygosity of the Z allele also increased the risk of cirrhosis (OR, 11.8; P = 1.8 × 10−09). The OR for cirrhosis of the Z allele was comparable to that of well‐established genetic variants, including patatin‐like phospholipase domain containing 3 (PNPLA3) I148M (OR, 1.48; P = 1.1 × 10−22) and transmembrane 6 superfamily member 2 (TM6SF2) E167K (OR, 1.34; P = 2.6 × 10−06). In heterozygotes compared to noncarriers, the Z allele was associated with higher alanine aminotransferase (ALT; P = = 4.6 × 10−46), aspartate aminotransferase (AST; P = 2.2 × 10−27), alkaline phosphatase (P = 3.3 × 10−43), gamma‐glutamyltransferase (P = 1.2 × 10−05), and total bilirubin (P = 6.4 × 10−06); Z allele homozygotes had even greater elevations in liver biochemistries. Body mass index (BMI) amplified the association of the Z allele for ALT (P interaction = 0.021) and AST (P interaction = 0.0040), suggesting a gene–environment interaction. Finally, we demonstrated genetic interactions between variants in PNPLA3, TM6SF2, and hydroxysteroid 17‐beta dehydrogenase 13 (HSD17B13); there was no evidence of epistasis between the Z allele and these variants. Conclusion: SERPINA1 Z allele heterozygosity is an important risk factor for liver disease; this risk is amplified by increasing BMI.

Genome-wide association studies (GWASs) have identified regions associated with chronic obstructive pulmonary disease (COPD). GWASs of other diseases have shown an approximately 10-fold overrepresentation of nonsynonymous variants, despite limited exonic coverage on genotyping arrays. We hypothesized that a large-scale analysis of coding variants could discover novel genetic associations with COPD, including rare variants with large effect sizes. We performed a meta-analysis of exome arrays from 218,399 controls and 33,851 moderate-to-severe COPD cases. All exome-wide significant associations were present in regions previously identified by GWAS. We did not identify any novel rare coding variants with large effect sizes. Within GWAS regions on chromosomes 5q, 6p, and 15q, four coding variants were conditionally significant ( P < 0.00015) when adjusting for lead GWAS single-nucleotide polymorphisms A common gasdermin B ( GSDMB) splice variant (rs11078928) previously associated with a decreased risk for asthma was nominally associated with a decreased risk for COPD [minor allele frequency (MAF) = 0.46, P = 1.8e-4]. Two stop variants in coiled-coil α-helical rod protein 1 ( CCHCR1), a gene involved in regulating cell proliferation, were associated with COPD (both P < 0.0001). The SERPINA1 Z allele was associated with a random-effects odds ratio of 1.43 for COPD (95% confidence interval = 1.17–1.74), though with marked heterogeneity across studies. Overall, COPD-associated exonic variants were identified in genes involved in DNA methylation, cell-matrix interactions, cell proliferation, and cell death. In conclusion, we performed the largest exome array meta-analysis of COPD to date and identified potential functional coding variants. Future studies are needed to identify rarer variants and further define the role of coding variants in COPD pathogenesis.

Background Secondary polycythemia is associated with cigarette smoking and chronic obstructive pulmonary disease (COPD). However, the prevalence of polycythemia in COPD and the contributing risk factors for polycythemia in COPD have not been extensively studied. Methods We analyzed the presence of secondary polycythemia in current and former smokers with moderate to very severe COPD at the five-year follow-up visit in the observational COPDGene study. We used logistic regression to evaluate the association of polycythemia with age, sex, race, altitude, current smoking status, spirometry, diffusing capacity for carbon monoxide (DLCO), quantitative chest CT measurements (including emphysema, airway wall thickness, and pulmonary artery to aorta diameter ratio), resting hypoxemia, exercise-induced hypoxemia, and long-term oxygen therapy. Results In a total of 1928 COPDGene participants with moderate to very severe COPD, secondary polycythemia was found in 97 (9.2%) male and 31 (3.5%) female participants. In a multivariable logistic model, severe resting hypoxemia (OR 3.50, 95% CI 1.41–8.66), impaired DLCO (OR 1.28 for each 10-percent decrease in DLCO % predicted, CI 1.09–1.49), male sex (OR 3.60, CI 2.20–5.90), non-Hispanic white race (OR 3.33, CI 1.71–6.50), current smoking (OR 2.55, CI 1.49–4.38), and enrollment in the Denver clinical center (OR 4.42, CI 2.38–8.21) were associated with higher risk for polycythemia. In addition, continuous (OR 0.13, CI 0.05–0.35) and nocturnal (OR 0.46, CI 0.21–0.97) supplemental oxygen were associated with lower risk for polycythemia. Results were similar after excluding participants with anemia and participants enrolled at the Denver clinical center. Conclusions In a large cohort of individuals with moderate to very severe COPD, male sex, current smoking, enrollment at the Denver clinical center, impaired DLCO, and severe hypoxemia were associated with increased risk for secondary polycythemia. Continuous or nocturnal supplemental oxygen use were associated with decreased risk for polycythemia.

Background Chronic obstructive pulmonary disease (COPD) is a highly morbid and heterogenous disease. While COPD is defined by spirometry, many COPD characteristics are seen in cigarette smokers with normal spirometry. The extent to which COPD and COPD heterogeneity is captured in omics of lung tissue is not known. Methods We clustered gene expression and methylation data in 78 lung tissue samples from former smokers with normal lung function or severe COPD. We applied two integrative omics clustering methods: (1) Similarity Network Fusion (SNF) and (2) Entropy-Based Consensus Clustering (ECC). Results SNF clusters were not significantly different by the percentage of COPD cases (48.8% vs. 68.6%, p = 0.13), though were different according to median forced expiratory volume in one second (FEV1) % predicted (82 vs. 31, p = 0.017). In contrast, the ECC clusters showed stronger evidence of separation by COPD case status (48.2% vs. 81.8%, p = 0.013) and similar stratification by median FEV1% predicted (82 vs. 30.5, p = 0.0059). ECC clusters using both gene expression and methylation were identical to the ECC clustering solution generated using methylation data alone. Both methods selected clusters with differentially expressed transcripts enriched for interleukin signaling and immunoregulatory interactions between lymphoid and non-lymphoid cells. Conclusions Unsupervised clustering analysis from integrated gene expression and methylation data in lung tissue resulted in clusters with modest concordance with COPD, though were enriched in pathways potentially contributing to COPD-related pathology and heterogeneity.

Summary Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380 ), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Background: Small studies have recently suggested that there are specific plasma metabolic signatures in chronic obstructive pulmonary disease (COPD), but there have been no large comprehensive study of metabolomic signatures in COPD that also integrate genetic variants. Materials and Methods: Fresh frozen plasma from 957 non-Hispanic white subjects in COPDGene was used to quantify 995 metabolites with Metabolon's global metabolomics platform. Metabolite associations with five COPD phenotypes (chronic bronchitis, exacerbation frequency, percent emphysema, post-bronchodilator forced expiratory volume at one second [FEV1]/forced vital capacity [FVC], and FEV1 percent predicted) were assessed. A metabolome-wide association study was performed to find genetic associations with metabolite levels. Significantly associated single-nucleotide polymorphisms were tested for replication with independent metabolomic platforms and independent cohorts. COPD phenotype-driven modules were identified in network analysis integrated with genetic associations to assess gene-metabolite-phenotype interactions. Results: Of metabolites tested, 147 (14.8%) were significantly associated with at least 1 COPD phenotype. Associations with airflow obstruction were enriched for diacylglycerols and branched chain amino acids. Genetic associations were observed with 109 (11%) metabolites, 72 (66%) of which replicated in an independent cohort. For 20 metabolites, more than 20% of variance was explained by genetics. A sparse network of COPD phenotype-driven modules was identified, often containing metabolites missed in previous testing. Of the 26 COPD phenotype-driven modules, 6 contained metabolites with significant met-QTLs, although little module variance was explained by genetics. Conclusion: A dysregulation of systemic metabolism was predominantly found in COPD phenotypes characterized by airflow obstruction, where we identified robust heritable effects on individual metabolite abundances. However, network analysis, which increased the statistical power to detect associations missed previously in classic regression analyses, revealed that the genetic influence on COPD phenotype-driven metabolomic modules was modest when compared with clinical and environmental factors.

Lung function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry GWAS meta-analysis of lung function to date, comprising 580,869 participants, 1020 independent association signals identified 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score (GRS) showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies (PheWAS) for selected associated variants, and trait and pathway-specific GRS to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008–2010) and SpiroMeta Consortium (multiple countries, 1947–2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic β4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.