The single-cell eQTLGen consortium

  1. Monique GP van der Wijst  Is a corresponding author
  2. Dylan H de Vries
  3. Hilde E Groot
  4. Gosia Trynka
  5. Chung-Chau Hon
  6. Marc-Jan Bonder
  7. Oliver Stegle
  8. Martijn Nawijn
  9. Youssef Idaghdour
  10. Pim van der Harst
  11. Chun J Ye
  12. Joseph Powell
  13. Fabian J Theis
  14. Ahmed Mahfouz
  15. Matthias Heinig
  16. Lude Franke
  1. University of Groningen, University Medical Center Groningen, Netherlands
  2. Wellcome Sanger Institute, United Kingdom
  3. RIKEN Center for Integrative Medical Sciences, Japan
  4. European Molecular Biology Laboratory, European Bioinformatics Institute, United Kingdom
  5. DKFZ, Germany
  6. New York University Abu Dhabi, United Arab Emirates
  7. University of California, San Francisco, United States
  8. Garvan Institute, Australia
  9. Helmholtz Zentrum München, Germany
  10. Leiden University Medical Center, Netherlands
  11. Institute of Computational Biology, Helmholtz Zentrum München, Technical University of Munich, Germany

Abstract

In recent years, functional genomics approaches combining genetic information with bulk RNA-sequencing data have identified the downstream expression effects of disease-associated genetic risk factors through so-called expression quantitative trait locus (eQTL) analysis. Single-cell RNA-sequencing creates enormous opportunities for mapping eQTLs across different cell types and in dynamic processes, many of which are obscured when using bulk methods. Rapid increase in throughput and reduction in cost per cell now allow this technology to be applied to large-scale population genetics studies. To fully leverage these emerging data resources, we have founded the single-cell eQTLGen consortium (sc-eQTLGen), aimed at pinpointing the cellular contexts in which disease-causing genetic variants affect gene expression. Here, we outline the goals, approach and potential utility of the sc-eQTLGen consortium. We also provide a set of study design considerations for future single-cell eQTL studies.

Data availability

Not applicable

Article and author information

Author details

  1. Monique GP van der Wijst

    Genetics, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
    For correspondence
    m.g.p.van.der.wijst@umcg.nl
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1520-3970
  2. Dylan H de Vries

    Genetics, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
  3. Hilde E Groot

    Cardiology, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8265-3085
  4. Gosia Trynka

    Cellular Genetics, Wellcome Sanger Institute, Hinxton, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6955-9529
  5. Chung-Chau Hon

    Genome Information Analysis, RIKEN Center for Integrative Medical Sciences, Yokahama, Japan
    Competing interests
    The authors declare that no competing interests exist.
  6. Marc-Jan Bonder

    Wellcome Trust Genome Campus, European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8431-3180
  7. Oliver Stegle

    DKFZ, Heidelberg, Germany
    Competing interests
    The authors declare that no competing interests exist.
  8. Martijn Nawijn

    Department of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3372-6521
  9. Youssef Idaghdour

    Program in Biology, Public Health Research Center, New York University Abu Dhabi, Abu Dhabi, United Arab Emirates
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2768-9376
  10. Pim van der Harst

    Department of Cardiology, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2713-686X
  11. Chun J Ye

    Division of Rheumatology, Department of Medicine, Department of Bioengineering and Therapeutic Sciences, Department of Epidemiology and Biostatistics, University of California, San Francisco, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
  12. Joseph Powell

    Garvan Institute, Sydney, Australia
    Competing interests
    The authors declare that no competing interests exist.
  13. Fabian J Theis

    Institute of Computational Biology, German Research Center for Environmental Health, Helmholtz Zentrum München, Neuherberg, Germany
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2419-1943
  14. Ahmed Mahfouz

    Single cell analysis, Leiden University Medical Center, Leiden, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8601-2149
  15. Matthias Heinig

    Germany Department of Informatics, Institute of Computational Biology, Helmholtz Zentrum München, Technical University of Munich, Neuherberg, München, Germany
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5612-1720
  16. Lude Franke

    Genetics, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
    Competing interests
    The authors declare that no competing interests exist.

Funding

Dutch Research Council (NWO-Veni 192.029)

  • Monique GP van der Wijst

Dutch Research Council (ZonMW-VIDI 917.14.374)

  • Lude Franke

European Research Council (ERC Starting grant Immrisk 637640)

  • Lude Franke

Oncode Institute

  • Lude Franke

National Health and Medical Research Council (Investigator grant 1175781)

  • Joseph Powell

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

Reviewing Editor

  1. Helena Pérez Valle, eLife, United Kingdom

Publication history

  1. Received: September 24, 2019
  2. Accepted: March 3, 2020
  3. Accepted Manuscript published: March 9, 2020 (version 1)
  4. Version of Record published: March 17, 2020 (version 2)

Copyright

© 2020, van der Wijst 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

  • 9,024
    Page views
  • 1,090
    Downloads
  • 67
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

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. Monique GP van der Wijst
  2. Dylan H de Vries
  3. Hilde E Groot
  4. Gosia Trynka
  5. Chung-Chau Hon
  6. Marc-Jan Bonder
  7. Oliver Stegle
  8. Martijn Nawijn
  9. Youssef Idaghdour
  10. Pim van der Harst
  11. Chun J Ye
  12. Joseph Powell
  13. Fabian J Theis
  14. Ahmed Mahfouz
  15. Matthias Heinig
  16. Lude Franke
(2020)
The single-cell eQTLGen consortium
eLife 9:e52155.
https://doi.org/10.7554/eLife.52155

Further reading

    1. Chromosomes and Gene Expression
    2. Genetics and Genomics
    Joseph V Geisberg, Zarmik Moqtaderi ... Kevin Struhl
    Research Advance

    Alternative polyadenylation yields many mRNA isoforms whose 3' termini occur disproportionately in clusters within 3' UTRs. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. Within yeast isoform clusters, these shifts occur steadily from one isoform to the next across nucleotide distances. In contrast, the shift between clusters from the last isoform of one cluster to the first isoform of the next - is much less pronounced, even over large distances. GC content in a region 13-30 nt downstream from isoform clusters correlates with their sensitivity to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at the nucleotide level within clusters, but not between them. Pol II occupancy increases just downstream of the most speed-sensitive poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that 1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, 2) poly(A) site clusters are linked to the local elongation rate and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, 3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and 4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.

    1. Genetics and Genomics
    Amel Lamri, Jayneel Limbachia ... Sonia S Anand
    Research Article Updated

    South Asian women are at increased risk of developing gestational diabetes mellitus (GDM). Few studies have investigated the genetic contributions to GDM risk. We investigated the association of a type 2 diabetes (T2D) polygenic risk score (PRS), on its own, and with GDM risk factors, on GDM-related traits using data from two birth cohorts in which South Asian women were enrolled during pregnancy. 837 and 4372 pregnant South Asian women from the SouTh Asian BiRth CohorT (START) and Born in Bradford (BiB) cohort studies underwent a 75-g glucose tolerance test. PRSs were derived using genome-wide association study results from an independent multi-ethnic study (~18% South Asians). Associations with fasting plasma glucose (FPG); 2 hr post-load glucose (2hG); area under the curve glucose; and GDM were tested using linear and logistic regressions. The population attributable fraction (PAF) of the PRS was calculated. Every 1 SD increase in the PRS was associated with a 0.085 mmol/L increase in FPG ([95% confidence interval, CI=0.07–0.10], p=2.85×10−20); 0.21 mmol/L increase in 2hG ([95% CI=0.16–0.26], p=5.49×10−16); and a 45% increase in the risk of GDM ([95% CI=32–60%], p=2.27×10−14), independent of parental history of diabetes and other GDM risk factors. PRS tertile 3 accounted for 12.5% of the population’s GDM alone, and 21.7% when combined with family history. A few weak PRS and GDM risk factors interactions modulating FPG and GDM were observed. Taken together, these results show that a T2D PRS and family history of diabetes are strongly and independently associated with multiple GDM-related traits in women of South Asian descent, an effect that could be modulated by other environmental factors.