Self-assembling manifolds in single-cell RNA sequencing data
- Stanford University, United States
Abstract
Single-cell RNA sequencing has spurred the development of computational methods that enable researchers to classify cell types, delineate developmental trajectories, and measure molecular responses to external perturbations. Many of these technologies rely on their ability to detect genes whose cell-to-cell variations arise from the biological processes of interest rather than transcriptional or technical noise. However, for datasets in which the biologically relevant differences between cells are subtle, identifying these genes is challenging. We present the self-assembling manifold (SAM) algorithm, an iterative soft feature selection strategy to quantify gene relevance and improve dimensionality reduction. We demonstrate its advantages over other state-of-the-art methods with experimental validation in identifying novel stem cell populations of Schistosoma mansoni, a prevalent parasite that infects hundreds of millions of people. Extending our analysis to a total of 56 datasets, we show that SAM is generalizable and consistently outperforms other methods in a variety of biological and quantitative benchmarks.
Data availability
The schistosome stem cell scRNAseq data generated in this study is available through the Gene Expression Omnibus (GEO) under accession number GSE116920.
-
Single-cell RNA sequencing of proliferative stem cell population from juvenile Schistosoma mansoniNCBI Gene Expression Omnibus, GSE116920.
-
Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cellsNCBI Gene Expression Omnibus, GSE36552.
-
Adult mouse cortical cell taxonomy revealed by single cell transcriptomicsNCBI Gene Expression Omnibus, GSE71585-GPL17021.
-
The transcriptome and DNA methylome landscapes of human primordial germ cellsNCBI Gene Expression Omnibus, GSE63818.
-
Single-cell analysis uncovers clonal acinar cell heterogeneity in the adult pancreasNCBI Gene Expression Omnibus, GSE80032.
-
Single-cell RNA-seq reveals dynamic, random monoallelic gene expression in mammalian cellsNCBI Gene Expression Omnibus, GSE45719.
-
Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastomaNCBI Gene Expression Omnibus, GSE57872.
-
Single-cell topological RNA-seq analysis reveals insights into cellular differentiation and developmentNCBI Gene Expression Omnibus, GSE94883.
-
Innate-like functions of natural killer T cell subsets result from highly divergent gene programsNCBI Gene Expression Omnibus, GSE74596.
-
Single-cell RNA-Seq resolves cellular complexity in sensory organs from the neonatal inner earNCBI Gene Expression Omnibus, GSE71982.
-
Cell fate inclination within 2-cell and 4-cell mouse embryos revealed by single-cell RNA sequencingNCBI Gene Expression Omnibus, GSE57249.
-
The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cellsNCBI Gene Expression Omnibus, GSE52529-GPL16791.
-
Single-cell RNA-seq reveals dynamic paracrine control of cellular variationNCBI Gene Expression Omnibus, GSE48968-GPL13112.
-
Oscope identifies oscillatory genes in unsynchronized single-cell RNAseq experimentsNCBI Gene Expression Omnibus, GSE64016.
-
Reconstructing lineage hierarchies of the distal lung epithelium using single-cell RNA-seqNCBI Gene Expression Omnibus, GSE52583-GPL13112.
-
Single-cell analysis of mixed-lineage states leading to a binary cell fate choiceNCBI Gene Expression Omnibus, GSE70245.
-
Single-cell RNA-seq with waterfall reveals molecular cascades underlying adult neurogenesisNCBI Gene Expression Omnibus, GSE71485.
-
Deciphering cell lineage specification during male sex determination with single-cell RNA sequencingNCBI Gene Expression Omnibus, GSE97519.
-
A molecular atlas of cell types and zonation in the brain vasculatureNCBI Gene Expression Omnibus, GSE99235.
-
Multipotent peripheral glial cells generate neuroendocrine cells of the adrenal medullaNCBI Gene Expression Omnibus, GSE99933.
-
Defining the earliest step of cardiovascular lineage segregation by single-cell RNA-seqNCBI Gene Expression Omnibus, GSE100471.
-
Temporal tracking of microglia activation in neurodegeneration at single-cell resolutionNCBI Gene Expression Omnibus, GSE103334.
-
Early emergence of cortical interneuron diversity in the mouse embryoNCBI Gene Expression Omnibus, GSE109796.
-
Single-cell transcriptional dynamics of flavivirus infectionNCBI Gene Expression Omnibus, GSE110496.
-
Single-cell RNA-seq supports a developmental hierarchy in human oligodendrogliomaNCBI Gene Expression Omnibus, GSE70630.
-
Transcriptional heterogeneity and lineage commitment in myeloid progenitorsNCBI Gene Expression Omnibus, GSE72857.
-
A survey of human brain transcriptome diversity at the single cell levelNCBI Gene Expression Omnibus, GSE67835.
-
Single-cell transcriptomics of the human endocrine pancreasNCBI Gene Expression Omnibus, GSE83139.
-
Cell type transcriptome atlas for the planarian Schmidtea mediterraneaNCBI Gene Expression Omnibus, GSE111764.
-
A Single-cell transcriptome atlas of the human pancreasNCBI Gene Expression Omnibus, GSE85241.
Article and author information
Author details
Funding
Burroughs Wellcome Fund
- Bo Wang
Arnold and Mabel Beckman Foundation
- Bo Wang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Animal experimentation: In adherence to the Animal Welfare Act and the Public Health Service Policy on Humane Care and Use of Laboratory Animals, all experiments with and care of mice were performed in accordance with protocols approved by the Institutional Animal Care and Use Committees (IACUC) of Stanford University (protocol approval number 30366).
Copyright
© 2019, Tarashansky 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
-
- 8,695
- views
-
- 1,083
- downloads
-
- 66
- 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)
Self-assembling manifolds in single-cell RNA sequencing data
eLife 8:e48994.
https://doi.org/10.7554/eLife.48994
Further reading
-
- Computational and Systems Biology
Plasmid construction is central to life science research, and sequence verification is arguably its costliest step. Long-read sequencing has emerged as a competitor to Sanger sequencing, with the principal benefit that whole plasmids can be sequenced in a single run. Nevertheless, the current cost of nanopore sequencing is still prohibitive for routine sequencing during plasmid construction. We develop a computational approach termed Simple Algorithm for Very Efficient Multiplexing of Oxford Nanopore Experiments for You (SAVEMONEY) that guides researchers to mix multiple plasmids and subsequently computationally de-mixes the resultant sequences. SAVEMONEY defines optimal mixtures in a pre-survey step, and following sequencing, executes a post-analysis workflow involving sequence classification, alignment, and consensus determination. By using Bayesian analysis with prior probability of expected plasmid construction error rate, high-confidence sequences can be obtained for each plasmid in the mixture. Plasmids differing by as little as two bases can be mixed as a single sample for nanopore sequencing, and routine multiplexing of even six plasmids per 180 reads can still maintain high accuracy of consensus sequencing. SAVEMONEY should further democratize whole-plasmid sequencing by nanopore and related technologies, driving down the effective cost of whole-plasmid sequencing to lower than that of a single Sanger sequencing run.
-
- Biochemistry and Chemical Biology
- Computational and Systems Biology
Protein–protein interactions are fundamental to understanding the molecular functions and regulation of proteins. Despite the availability of extensive databases, many interactions remain uncharacterized due to the labor-intensive nature of experimental validation. In this study, we utilized the AlphaFold2 program to predict interactions among proteins localized in the nuage, a germline-specific non-membrane organelle essential for piRNA biogenesis in Drosophila. We screened 20 nuage proteins for 1:1 interactions and predicted dimer structures. Among these, five represented novel interaction candidates. Three pairs, including Spn-E_Squ, were verified by co-immunoprecipitation. Disruption of the salt bridges at the Spn-E_Squ interface confirmed their functional importance, underscoring the predictive model’s accuracy. We extended our analysis to include interactions between three representative nuage components—Vas, Squ, and Tej—and approximately 430 oogenesis-related proteins. Co-immunoprecipitation verified interactions for three pairs: Mei-W68_Squ, CSN3_Squ, and Pka-C1_Tej. Furthermore, we screened the majority of Drosophila proteins (~12,000) for potential interaction with the Piwi protein, a central player in the piRNA pathway, identifying 164 pairs as potential binding partners. This in silico approach not only efficiently identifies potential interaction partners but also significantly bridges the gap by facilitating the integration of bioinformatics and experimental biology.
