Predicting mechanisms of action at genetic loci associated with discordant effects on type 2 diabetes and abdominal fat accumulation

  1. Yonathan Tamrat Aberra  Is a corresponding author
  2. Lijiang Ma
  3. Johan LM Björkegren
  4. Mete Civelek  Is a corresponding author
  1. University of Virginia, United States
  2. Icahn School of Medicine at Mount Sinai, United States

Abstract

Metabolic syndrome (MetSyn) is a cluster of dysregulated metabolic conditions that occur together to increase the risk for cardiometabolic disorders such as type 2 diabetes (T2D). One key condition associated with MetSyn, abdominal obesity, is measured by computing the ratio of waist-to-hip circumference adjusted for the body-mass index (WHRadjBMI). WHRadjBMI and T2D are complex traits with genetic and environmental components, which has enabled genome-wide association studies (GWAS) to identify hundreds of loci associated with both. Statistical genetics analyses of these GWAS have predicted that WHRadjBMI is a strong causal risk factor of T2D and that these traits share genetic architecture at many loci. To date, no variants have been described that are simultaneously associated with protection from T2D but with increased abdominal obesity. Here, we used colocalization analysis to identify genetic variants with a shared association for T2D and abdominal obesity. This analysis revealed the presence of five loci associated with discordant effects on T2D and abdominal obesity. The alleles of the lead genetic variants in these loci that were protective against T2D were also associated with increased abdominal obesity. We further used publicly available expression, epigenomic, and genetic regulatory data to predict the effector genes (eGenes) and functional tissues at the 2p21, 5q21.1, and 19q13.11 loci. We also computed the correlation between the subcutaneous adipose tissue (SAT) expression of predicted effector genes (eGenes) with metabolic phenotypes and adipogenesis. We proposed a model to resolve the discordant effects at the 5q21.1 locus. We find that eGenes gypsy retrotransposon integrase 1 (GIN1), diphosphoinositol pentakisphosphate kinase 2 (PPIP5K2), and peptidylglycine alpha-amidating monooxygenase (PAM) represent the likely causal eGenes at the 5q21.1 locus. Taken together, these results are the first to describe a potential mechanism through which a genetic variant can confer increased abdominal obesity but protection from T2D risk. Understanding precisely how and which genetic variants confer increased risk for MetSyn will develop the basic science needed to design novel therapeutics for metabolic syndrome.

Data availability

The current manuscript is a computational investigation using publicly available data, so no data have been generated for this manuscript. All publicly obtained data sets are included in Supplementary Table 1. All analysis and figure-generating code uploaded to the following Github repository: https://github.com/aberrations/predicting-functional-mechanisms-discordant-loci.

The following previously published data sets were used
    1. Raulerson CK
    2. Ko A
    3. Kidd JC
    4. et al.
    (2019) METSIM eQTL
    FTP, https://doi.org/10.1016/j.ajhg.2019.09.001.

Article and author information

Author details

  1. Yonathan Tamrat Aberra

    Department of Biomedical Engineering, University of Virginia, Charlottesville, United States
    For correspondence
    ya8eb@virginia.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6055-2291
  2. Lijiang Ma

    Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Johan LM Björkegren

    Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Mete Civelek

    Department of Biomedical Engineering, University of Virginia, Charlottesville, United States
    For correspondence
    mete@virginia.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8141-0284

Funding

National Heart, Lung, and Blood Institute (2T32HL007284-46)

  • Yonathan Tamrat Aberra

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

Reviewing Editor

  1. Tune H Pers, University of Copenhagen, Denmark

Version history

  1. Received: April 28, 2022
  2. Preprint posted: April 29, 2022 (view preprint)
  3. Accepted: May 31, 2023
  4. Accepted Manuscript published: June 30, 2023 (version 2)
  5. Version of Record published: June 16, 2023 (version 1)

Copyright

© 2023, Aberra 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

  • 995
    views
  • 171
    downloads
  • 2
    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. Yonathan Tamrat Aberra
  2. Lijiang Ma
  3. Johan LM Björkegren
  4. Mete Civelek
(2023)
Predicting mechanisms of action at genetic loci associated with discordant effects on type 2 diabetes and abdominal fat accumulation
eLife 12:e79834.
https://doi.org/10.7554/eLife.79834

Share this article

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

Further reading

    1. Computational and Systems Biology
    Maksim Kleverov, Daria Zenkova ... Alexey A Sergushichev
    Research Article

    Transcriptomic profiling became a standard approach to quantify a cell state, which led to accumulation of huge amount of public gene expression datasets. However, both reuse of these datasets or analysis of newly generated ones requires significant technical expertise. Here we present Phantasus - a user-friendly web-application for interactive gene expression analysis which provides a streamlined access to more than 96000 public gene expression datasets, as well as allows analysis of user-uploaded datasets. Phantasus integrates an intuitive and highly interactive JavaScript-based heatmap interface with an ability to run sophisticated R-based analysis methods. Overall Phantasus allows users to go all the way from loading, normalizing and filtering data to doing differential gene expression and downstream analysis. Phantasus can be accessed on-line at https://alserglab.wustl.edu/phantasus or can be installed locally from Bioconductor (https://bioconductor.org/packages/phantasus). Phantasus source code is available at https://github.com/ctlab/phantasus under MIT license.

    1. Computational and Systems Biology
    2. Evolutionary Biology
    Ryan T Bell, Harutyun Sahakyan ... Eugene V Koonin
    Research Article

    A comprehensive census of McrBC systems, among the most common forms of prokaryotic Type IV restriction systems, followed by phylogenetic analysis, reveals their enormous abundance in diverse prokaryotes and a plethora of genomic associations. We focus on a previously uncharacterized branch, which we denote coiled-coil nuclease tandems (CoCoNuTs) for their salient features: the presence of extensive coiled-coil structures and tandem nucleases. The CoCoNuTs alone show extraordinary variety, with three distinct types and multiple subtypes. All CoCoNuTs contain domains predicted to interact with translation system components, such as OB-folds resembling the SmpB protein that binds bacterial transfer-messenger RNA (tmRNA), YTH-like domains that might recognize methylated tmRNA, tRNA, or rRNA, and RNA-binding Hsp70 chaperone homologs, along with RNases, such as HEPN domains, all suggesting that the CoCoNuTs target RNA. Many CoCoNuTs might additionally target DNA, via McrC nuclease homologs. Additional restriction systems, such as Type I RM, BREX, and Druantia Type III, are frequently encoded in the same predicted superoperons. In many of these superoperons, CoCoNuTs are likely regulated by cyclic nucleotides, possibly, RNA fragments with cyclic termini, that bind associated CARF (CRISPR-Associated Rossmann Fold) domains. We hypothesize that the CoCoNuTs, together with the ancillary restriction factors, employ an echeloned defense strategy analogous to that of Type III CRISPR-Cas systems, in which an immune response eliminating virus DNA and/or RNA is launched first, but then, if it fails, an abortive infection response leading to PCD/dormancy via host RNA cleavage takes over.