1. Genetics and Genomics

Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation

  1. Courtney J Smith  Is a corresponding author
  2. Nasa Sinnott-Armstrong  Is a corresponding author
  3. Anna Cichońska
  4. Heli Julkunen
  5. Eric B Fauman
  6. Peter Würtz
  7. Jonathan K Pritchard  Is a corresponding author
  1. Stanford University, United States
  2. Fred Hutchinson Cancer Research Center, United States
  3. Nightingale Health Plc, Finland
  4. Pfizer, United States
https://doi.org/10.7554/eLife.79348
Share this article
Cite this article
  1. Courtney J Smith
  2. Nasa Sinnott-Armstrong
  3. Anna Cichońska
  4. Heli Julkunen
  5. Eric B Fauman
  6. Peter Würtz
  7. Jonathan K Pritchard
(2022)
Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation
eLife 11:e79348.
https://doi.org/10.7554/eLife.79348
  2. Download BibTeX
  3. Download .RIS

Abstract

Pleiotropy and genetic correlation are widespread features in GWAS, but they are often difficult to interpret at the molecular level. Here, we perform GWAS of 16 metabolites clustered at the intersection of amino acid catabolism, glycolysis, and ketone body metabolism in a subset of UK Biobank. We utilize the well-documented biochemistry jointly impacting these metabolites to analyze pleiotropic effects in the context of their pathways. Among the 213 lead GWAS hits, we find a strong enrichment for genes encoding pathway-relevant enzymes and transporters. We demonstrate that the effect directions of variants acting on biology between metabolite pairs often contrast with those of upstream or downstream variants as well as the polygenic background. Thus, we find that these outlier variants often reflect biology local to the traits. Finally, we explore the implications for interpreting disease GWAS, underscoring the potential of unifying biochemistry with dense metabolomics data to understand the molecular basis of pleiotropy in complex traits and diseases.

Data availability

The source data and analyzed data have been deposited in Dryad. Code are available at the github link (https://github.com/courtrun/Pleiotropy-of-UKB-Metabolites). The raw individual level data are available through application to UK Biobank.

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Courtney J Smith

    Department of Genetics, Stanford University, Stanford, United States
    For correspondence
    courtrun@stanford.edu
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7812-0083

  2. Nasa Sinnott-Armstrong

    Herbold Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, United States
    For correspondence
    nasa@fredhutch.org
    Competing interests
    No competing interests declared.

  3. Anna Cichońska

    Nightingale Health Plc, Helsinki, Finland
    Competing interests
    Anna Cichońska, is a former employee and holds stock options with Nightingale Health Plc..

  4. Heli Julkunen

    Nightingale Health Plc, Helsinki, Finland
    Competing interests
    Heli Julkunen, is an employee and holds stock options with Nightingale Health Plc..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4282-0248

  5. Eric B Fauman

    Internal Medicine Research Unit, Pfizer, Cambridge, United States
    Competing interests
    Eric B Fauman, is affiliated with Pfizer Worldwide Research, has no financial interests to declare, contributed as an individual and the work was not part of a Pfizer collaboration nor was it funded by Pfizer..

  6. Peter Würtz

    Nightingale Health Plc, Helsinki, Finland
    Competing interests
    Peter Würtz, is an employee and shareholder of Nightingale Health Plc..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5832-0221

  7. Jonathan K Pritchard

    Departments of Genetics and Biology, Stanford University, Stanford, United States
    For correspondence
    pritch@stanford.edu
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8828-5236

Funding

Stanford Knight-Hennessy Scholars Program (Graduate Student Fellowship)

  • Courtney J Smith

National Science Foundation (Graduate Student Fellowship)

  • Courtney J Smith

National Institute of Health (5R01HG011432 and 5R01AG066490)

  • Jonathan K Pritchard

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Human subjects: All participants provided written informed consent and ethical approval was obtained from the North West Multi-Center Research Ethics Committee (11/NW/0382). The current analysis was approved under UK Biobank Project 24983 and 30418.

Copyright

© 2022, Smith et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,224
    views
  • 723
    downloads
  • 27
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Courtney J Smith
  2. Nasa Sinnott-Armstrong
  3. Anna Cichońska
  4. Heli Julkunen
  5. Eric B Fauman
  6. Peter Würtz
  7. Jonathan K Pritchard
(2022)
Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation
eLife 11:e79348.
https://doi.org/10.7554/eLife.79348

Share this article

https://doi.org/10.7554/eLife.79348

Categories and tags

Research organism

Further reading

    1. Genetics and Genomics

    A transcription network underlies the dual genomic coordination of mitochondrial biogenesis

    Fan Zhang, Annie Lee ... Hong Xu
    Research Article

    Mitochondrial biogenesis requires the expression of genes encoded by both the nuclear and mitochondrial genomes. However, aside from a handful transcription factors regulating specific subsets of mitochondrial genes, the overall architecture of the transcriptional control of mitochondrial biogenesis remains to be elucidated. The mechanisms coordinating these two genomes are largely unknown. We performed a targeted RNAi screen in developing eyes with reduced mitochondrial DNA content, anticipating a synergistic disruption of tissue development due to impaired mitochondrial biogenesis and mitochondrial DNA (mtDNA) deficiency. Among 638 transcription factors annotated in the Drosophila genome, 77 were identified as potential regulators of mitochondrial biogenesis. Utilizing published ChIP-seq data of positive hits, we constructed a regulatory network revealing the logic of the transcription regulation of mitochondrial biogenesis. Multiple transcription factors in core layers had extensive connections, collectively governing the expression of nearly all mitochondrial genes, whereas factors sitting on the top layer may respond to cellular cues to modulate mitochondrial biogenesis through the underlying network. CG1603, a core component of the network, was found to be indispensable for the expression of most nuclear mitochondrial genes, including those required for mtDNA maintenance and gene expression, thus coordinating nuclear genome and mtDNA activities in mitochondrial biogenesis. Additional genetic analyses validated YL-1, a transcription factor upstream of CG1603 in the network, as a regulator controlling CG1603 expression and mitochondrial biogenesis.

    1. Genetics and Genomics

    Accurate predictions of SARS-CoV-2 infectivity from comprehensive analysis

    Jongkeun Park, WonJong Choi ... Dongwan Hong
    Research Article

    An unprecedented amount of SARS-CoV-2 data has been accumulated compared with previous infectious diseases, enabling insights into its evolutionary process and more thorough analyses. This study investigates SARS-CoV-2 features as it evolved to evaluate its infectivity. We examined viral sequences and identified the polarity of amino acids in the receptor binding motif (RBM) region. We detected an increased frequency of amino acid substitutions to lysine (K) and arginine (R) in variants of concern (VOCs). As the virus evolved to Omicron, commonly occurring mutations became fixed components of the new viral sequence. Furthermore, at specific positions of VOCs, only one type of amino acid substitution and a notable absence of mutations at D467 were detected. We found that the binding affinity of SARS-CoV-2 lineages to the ACE2 receptor was impacted by amino acid substitutions. Based on our discoveries, we developed APESS, an evaluation model evaluating infectivity from biochemical and mutational properties. In silico evaluation using real-world sequences and in vitro viral entry assays validated the accuracy of APESS and our discoveries. Using Machine Learning, we predicted mutations that had the potential to become more prominent. We created AIVE, a web-based system, accessible at https://ai-ve.org to provide infectivity measurements of mutations entered by users. Ultimately, we established a clear link between specific viral properties and increased infectivity, enhancing our understanding of SARS-CoV-2 and enabling more accurate predictions of the virus.

    1. Cell Biology
    2. Genetics and Genomics

    Full-length direct RNA sequencing uncovers stress granule-dependent RNA decay upon cellular stress

    Showkat Ahmad Dar, Sulochan Malla ... Manolis Maragkakis
    Research Article

    Cells react to stress by triggering response pathways, leading to extensive alterations in the transcriptome to restore cellular homeostasis. The role of RNA metabolism in shaping the cellular response to stress is vital, yet the global changes in RNA stability under these conditions remain unclear. In this work, we employ direct RNA sequencing with nanopores, enhanced by 5ʹ end adapter ligation, to comprehensively interrogate the human transcriptome at single-molecule and -nucleotide resolution. By developing a statistical framework to identify robust RNA length variations in nanopore data, we find that cellular stress induces prevalent 5ʹ end RNA decay that is coupled to translation and ribosome occupancy. Unlike typical RNA decay models in normal conditions, we show that stress-induced RNA decay is dependent on XRN1 but does not depend on deadenylation or decapping. We observed that RNAs undergoing decay are predominantly enriched in the stress granule transcriptome while inhibition of stress granule formation via genetic ablation of G3BP1 and G3BP2 rescues RNA length. Our findings reveal RNA decay as a key component of RNA metabolism upon cellular stress that is dependent on stress granule formation.