Deep learning models predict regulatory variants in pancreatic islets and refine type 2 diabetes association signals

Abstract

Genome-wide association analyses have uncovered multiple genomic regions associated with T2D, but identification of the causal variants at these remains a challenge. There is growing interest in the potential of deep learning models - which predict epigenome features from DNA sequence - to support inference concerning the regulatory effects of disease-associated variants. Here, we evaluate the advantages of training convolutional neural network (CNN) models on a broad set of epigenomic features collected in a single disease-relevant tissue – pancreatic islets in the case of type 2 diabetes (T2D) - as opposed to models trained on multiple human tissues. We report convergence of CNN-based metrics of regulatory function with conventional approaches to variant prioritization – genetic fine-mapping and regulatory annotation enrichment. We demonstrate that CNN-based analyses can refine association signals at T2D-associated loci and provide experimental validation for one such signal. We anticipate that these approaches will become routine in downstream analyses of GWAS.

Data availability

The datasets analysed during the current study are available in the public repositories under accessions listed in STable 1.

The following previously published data sets were used

Article and author information

Author details

  1. Agata Wesolowska-Andersen

    The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8688-2814
  2. Grace Zhuo Yu

    Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
  3. Vibe Nylander

    Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
  4. Fernando Abaitua

    The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
  5. Matthias Thurner

    The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7329-9769
  6. Jason M Torres

    The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7537-7035
  7. Anubha Mahajan

    The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
  8. Anna L Gloyn

    The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1205-1844
  9. Mark I McCarthy

    Department of Human Genetics, Genentech, South San Francisco, United States
    For correspondence
    mark.mccarthy@drl.ox.ac.uk
    Competing interests
    Mark I McCarthy, Senior editor, eLife; has served on advisory panels for Pfizer, NovoNordisk, Zoe Global; has received honoraria from Merck, Pfizer, NovoNordisk and Eli Lilly; has stock options in Zoe Global and has received research funding from Abbvie, Astra Zeneca, Boehringer Ingelheim, Eli Lilly, Janssen, Merck, NovoNordisk, Pfizer, Roche, Sanofi Aventis, Servier & Takeda. As of June 2019, is an employee of Genentech, and holds stock in Roche.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4393-0510

Funding

Wellcome (099673/Z/12/Z)

  • Matthias Thurner

Horizon 2020 Framework Programme (T2D Systems)

  • Anna L Gloyn

NIH Clinical Center (U01-DK105535)

  • Anna L Gloyn
  • Mark I McCarthy

NIH Clinical Center (U01-DK085545)

  • Anna L Gloyn

National Institute for Health Research (NF-SI-0617-10090)

  • Anna L Gloyn

Wellcome (090532)

  • Mark I McCarthy

Wellcome (106130)

  • Anna L Gloyn
  • Mark I McCarthy

Wellcome (098381)

  • Mark I McCarthy

Wellcome (203141)

  • Anna L Gloyn
  • Mark I McCarthy

Wellcome (212259)

  • Mark I McCarthy

Wellcome (095101)

  • Anna L Gloyn

Wellcome (200837)

  • Anna L Gloyn

Medical Research Council (MR/L020149/1)

  • Anna L Gloyn

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

Copyright

© 2020, Wesolowska-Andersen 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,166
    views
  • 631
    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. Agata Wesolowska-Andersen
  2. Grace Zhuo Yu
  3. Vibe Nylander
  4. Fernando Abaitua
  5. Matthias Thurner
  6. Jason M Torres
  7. Anubha Mahajan
  8. Anna L Gloyn
  9. Mark I McCarthy
(2020)
Deep learning models predict regulatory variants in pancreatic islets and refine type 2 diabetes association signals
eLife 9:e51503.
https://doi.org/10.7554/eLife.51503

Share this article

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

Further reading

    1. Computational and Systems Biology
    2. Genetics and Genomics
    Jia-Ying Su, Yun-Lin Wang ... Chien-Ling Lin
    Research Article

    Untranslated regions (UTRs) contain crucial regulatory elements for RNA stability, translation and localization, so their integrity is indispensable for gene expression. Approximately 3.7% of genetic variants associated with diseases occur in UTRs, yet a comprehensive understanding of UTR variant functions remains limited due to inefficient experimental and computational assessment methods. To systematically evaluate the effects of UTR variants on RNA stability, we established a massively parallel reporter assay on 6555 UTR variants reported in human disease databases. We examined the RNA degradation patterns mediated by the UTR library in two cell lines, and then applied LASSO regression to model the influential regulators of RNA stability. We found that UA dinucleotides and UA-rich motifs are the most prominent destabilizing element. Gain of UA dinucleotide outlined mutant UTRs with reduced stability. Studies on endogenous transcripts indicate that high UA-dinucleotide ratios in UTRs promote RNA degradation. Conversely, elevated GC content and protein binding on UA dinucleotides protect high-UA RNA from degradation. Further analysis reveals polarized roles of UA-dinucleotide-binding proteins in RNA protection and degradation. Furthermore, the UA-dinucleotide ratio of both UTRs is a common characteristic of genes in innate immune response pathways, implying a coordinated stability regulation through UTRs at the transcriptomic level. We also demonstrate that stability-altering UTRs are associated with changes in biobank-based health indices, underscoring the importance of precise UTR regulation for wellness. Our study highlights the importance of RNA stability regulation through UTR primary sequences, paving the way for further exploration of their implications in gene networks and precision medicine.

    1. Computational and Systems Biology
    2. Medicine
    Hong Yang, Cheng Zhang ... Adil Mardinoglu
    Research Article

    Excessive consumption of sucrose, in the form of sugar-sweetened beverages, has been implicated in the pathogenesis of metabolic dysfunction‐associated fatty liver disease (MAFLD) and other related metabolic syndromes. The c-Jun N-terminal kinase (JNK) pathway plays a crucial role in response to dietary stressors, and it was demonstrated that the inhibition of the JNK pathway could potentially be used in the treatment of MAFLD. However, the intricate mechanisms underlying these interventions remain incompletely understood given their multifaceted effects across multiple tissues. In this study, we challenged rats with sucrose-sweetened water and investigated the potential effects of JNK inhibition by employing network analysis based on the transcriptome profiling obtained from hepatic and extrahepatic tissues, including visceral white adipose tissue, skeletal muscle, and brain. Our data demonstrate that JNK inhibition by JNK-IN-5A effectively reduces the circulating triglyceride accumulation and inflammation in rats subjected to sucrose consumption. Coexpression analysis and genome-scale metabolic modeling reveal that sucrose overconsumption primarily induces transcriptional dysfunction related to fatty acid and oxidative metabolism in the liver and adipose tissues, which are largely rectified after JNK inhibition at a clinically relevant dose. Skeletal muscle exhibited minimal transcriptional changes to sucrose overconsumption but underwent substantial metabolic adaptation following the JNK inhibition. Overall, our data provides novel insights into the molecular basis by which JNK inhibition exerts its metabolic effect in the metabolically active tissues. Furthermore, our findings underpin the critical role of extrahepatic metabolism in the development of diet-induced steatosis, offering valuable guidance for future studies focused on JNK-targeting for effective treatment of MAFLD.