Gene regulatory network reconstruction using single-cell RNA sequencing of barcoded genotypes in diverse environments

Abstract

Understanding how gene expression programs are controlled requires identifying regulatory relationships between transcription factors and target genes. Gene regulatory networks are typically constructed from gene expression data acquired following genetic perturbation or environmental stimulus. Single-cell RNA sequencing (scRNAseq) captures the gene expression state of thousands of individual cells in a single experiment, offering advantages in combinatorial experimental design, large numbers of independent measurements, and accessing the interaction between the cell cycle and environmental responses that is hidden by population-level analysis of gene expression. To leverage these advantages, we developed a method for scRNAseq in budding yeast (Saccharomyces cerevisiae). We pooled diverse transcriptionally barcoded gene deletion mutants in 11 different environmental conditions and determined their expression state by sequencing 38,285 individual cells. We benchmarked a framework for learning gene regulatory networks from scRNAseq data that incorporates multitask learning and constructed a global gene regulatory network comprising 12,228 interactions.

Data availability

Sequencing data has been deposited in GEO: GSE125162Figures 2-7 (& Supplemental Figures 2-7) are generated from a single R markdown document. The scripts and all data necessary to do this analysis are provided as Supplemental Data 1. The raw output (knit HTML file) is provided as Supplemental Data 6.Interactive versions of several figures are available have been made available with the Shiny library in R: http://shiny.bio.nyu.edu/cj59/YeastSingleCell2019/The Inferelator package is available on GitHub and through python package managers (i.e. pip) under an open source license (BSD).

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

Article and author information

Author details

  1. Christopher A Jackson

    Center For Genomics and Systems Biology, New York University, New York, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8769-2710
  2. Dayanne M Castro

    Department of Biology, New York University, New York, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Giuseppe-Antonio Saldi

    Department of Biology, New York University, New York, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Richard Bonneau

    Center For Genomics and Systems Biology, New York University, New York, United States
    For correspondence
    bonneau@nyu.edu
    Competing interests
    The authors declare that no competing interests exist.
  5. David Gresham

    Center For Genomics and Systems Biology, New York University, New York, United States
    For correspondence
    dgresham@nyu.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4028-0364

Funding

National Institute of Diabetes and Digestive and Kidney Diseases (R01DK103358)

  • Richard Bonneau

National Institute of General Medical Sciences (R01GM107466)

  • David Gresham

National Science Foundation (MCB1818234)

  • David Gresham

Eunice Kennedy Shriver National Institute of Child Health and Human Development (R01HD096770)

  • Richard Bonneau

National Science Foundation (IOS1546218)

  • Richard Bonneau

National Cancer Institute (R01CA229235)

  • Richard Bonneau

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

Reviewing Editor

  1. Naama Barkai, Weizmann Institute of Science, Israel

Version history

  1. Received: August 21, 2019
  2. Accepted: January 10, 2020
  3. Accepted Manuscript published: January 27, 2020 (version 1)
  4. Version of Record published: February 6, 2020 (version 2)
  5. Version of Record updated: February 27, 2020 (version 3)

Copyright

© 2020, Jackson 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

  • 25,960
    Page views
  • 2,020
    Downloads
  • 72
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, 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. Christopher A Jackson
  2. Dayanne M Castro
  3. Giuseppe-Antonio Saldi
  4. Richard Bonneau
  5. David Gresham
(2020)
Gene regulatory network reconstruction using single-cell RNA sequencing of barcoded genotypes in diverse environments
eLife 9:e51254.
https://doi.org/10.7554/eLife.51254

Further reading

    1. Computational and Systems Biology
    2. Neuroscience
    George Zacharopoulos, Francesco Sella ... Roi Cohen Kadosh
    Research Article

    Previous research has highlighted the role of glutamate and gamma-aminobutyric acid (GABA) in perceptual, cognitive, and motor tasks. However, the exact involvement of these neurochemical mechanisms in the chain of information processing, and across human development, is unclear. In a cross-sectional longitudinal design, we used a computational approach to dissociate cognitive, decision, and visuomotor processing in 293 individuals spanning early childhood to adulthood. We found that glutamate and GABA within the intraparietal sulcus (IPS) explained unique variance in visuomotor processing, with higher glutamate predicting poorer visuomotor processing in younger participants but better visuomotor processing in mature participants, while GABA showed the opposite pattern. These findings, which were neurochemically, neuroanatomically and functionally specific, were replicated ~21 mo later and were generalized in two further different behavioral tasks. Using resting functional MRI, we revealed that the relationship between IPS neurochemicals and visuomotor processing is mediated by functional connectivity in the visuomotor network. We then extended our findings to high-level cognitive behavior by predicting fluid intelligence performance. We present evidence that fluid intelligence performance is explained by IPS GABA and glutamate and is mediated by visuomotor processing. However, this evidence was obtained using an uncorrected alpha and needs to be replicated in future studies. These results provide an integrative biological and psychological mechanistic explanation that links cognitive processes and neurotransmitters across human development and establishes their potential involvement in intelligent behavior.

    1. Computational and Systems Biology
    2. Neuroscience
    Huu Hoang, Shinichiro Tsutsumi ... Keisuke Toyama
    Research Article Updated

    Cerebellar climbing fibers convey diverse signals, but how they are organized in the compartmental structure of the cerebellar cortex during learning remains largely unclear. We analyzed a large amount of coordinate-localized two-photon imaging data from cerebellar Crus II in mice undergoing ‘Go/No-go’ reinforcement learning. Tensor component analysis revealed that a majority of climbing fiber inputs to Purkinje cells were reduced to only four functional components, corresponding to accurate timing control of motor initiation related to a Go cue, cognitive error-based learning, reward processing, and inhibition of erroneous behaviors after a No-go cue. Changes in neural activities during learning of the first two components were correlated with corresponding changes in timing control and error learning across animals, indirectly suggesting causal relationships. Spatial distribution of these components coincided well with boundaries of Aldolase-C/zebrin II expression in Purkinje cells, whereas several components are mixed in single neurons. Synchronization within individual components was bidirectionally regulated according to specific task contexts and learning stages. These findings suggest that, in close collaborations with other brain regions including the inferior olive nucleus, the cerebellum, based on anatomical compartments, reduces dimensions of the learning space by dynamically organizing multiple functional components, a feature that may inspire new-generation AI designs.