Identifying gene expression programs of cell-type identity and cellular activity with single-cell RNA-Seq
Abstract
Identifying gene expression programs underlying both cell-type identity and cellular activities (e.g. life-cycle processes, responses to environmental cues) is crucial for understanding the organization of cells and tissues. Although single-cell RNA-Seq (scRNA-Seq) can quantify transcripts in individual cells, each cell's expression profile may be a mixture of both types of programs, making them difficult to disentangle. Here we benchmark and enhance the use of matrix factorization to solve this problem. We show with simulations that a method we call consensus non-negative matrix factorization (cNMF) accurately infers identity and activity programs, including their relative contributions in each cell. To illustrate the insights this approach enables, we apply it to published brain organoid and visual cortex scRNA-Seq datasets; cNMF refines cell types and identifies both expected (e.g. cell cycle and hypoxia) and novel activity programs, including programs that may underlie a neurosecretory phenotype and synaptogenesis.
Data availability
All of the analyzed real datasets are publicly available and the relevant GEO accession codes are included in the manuscript. All of the simulated and real data can be accessed through Code Ocean at the following URL: https://doi.org/10.24433/CO.9044782e-cb96-4733-8a4f-bf42c21399e6
-
Cell diversity and network dynamics in photosensitive human brain organoids.Gene Expression Omnibus, GSE86153.
-
Adult mouse cortical cell taxonomy by single cell transcriptomicsGene Expression Omnibus, GSE71585.
Article and author information
Author details
Funding
National Institute of General Medical Sciences (T32GM007753)
- Dylan Kotliar
- Adrian Veres
- M Aurel Nagy
- Eran Hodis
National Institute of Allergy and Infectious Diseases (R01AI099210)
- Pardis C Sabeti
U.S. Food and Drug Administration (HHSF223201810172C)
- Dylan Kotliar
- Pardis C Sabeti
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Alfonso Valencia, Barcelona Supercomputing Center - BSC, Spain
Publication history
- Received: November 21, 2018
- Accepted: July 7, 2019
- Accepted Manuscript published: July 8, 2019 (version 1)
- Version of Record published: July 18, 2019 (version 2)
Copyright
© 2019, Kotliar 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
-
- 23,336
- Page views
-
- 2,138
- Downloads
-
- 59
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Further reading
-
- Computational and Systems Biology
- Evolutionary Biology
Studies of protein fitness landscapes reveal biophysical constraints guiding protein evolution and empower prediction of functional proteins. However, generalisation of these findings is limited due to scarceness of systematic data on fitness landscapes of proteins with a defined evolutionary relationship. We characterized the fitness peaks of four orthologous fluorescent proteins with a broad range of sequence divergence. While two of the four studied fitness peaks were sharp, the other two were considerably flatter, being almost entirely free of epistatic interactions. Mutationally robust proteins, characterized by a flat fitness peak, were not optimal templates for machine-learning-driven protein design – instead, predictions were more accurate for fragile proteins with epistatic landscapes. Our work paves insights for practical application of fitness landscape heterogeneity in protein engineering.
-
- Computational and Systems Biology
- Evolutionary Biology
Using a neural network to predict how green fluorescent proteins respond to genetic mutations illuminates properties that could help design new proteins.