Identification of cell types in a mouse brain single-cell atlas using low sampling coverage

BackgroundHigh throughput methods for profiling the transcriptomes of single cells have recently emerged as transformative approaches for large-scale population surveys of cellular diversity in heterogeneous primary tissues. However, the efficient generation of such atlases will depend on sufficient sampling of diverse cell types while remaining cost-effective to enable a comprehensive examination of organs, developmental stages, and individuals.ResultsTo examine the relationship between sampled cell numbers and transcriptional heterogeneity in the context of unbiased cell type classification, we explored the population structure of a publicly available 1.3 million cell dataset from E18.5 mouse brain and validated our findings in published data from adult mice. We propose a computational framework for inferring the saturation point of cluster discovery in a single-cell mRNA-seq experiment, centered around cluster preservation in downsampled datasets. In addition, we introduce a “complexity index,” which characterizes the heterogeneity of cells in a given dataset. Using Cajal-Retzius cells as an example of a limited complexity dataset, we explored whether the detected biological distinctions relate to technical clustering. Surprisingly, we found that clustering distinctions carrying biologically interpretable meaning are achieved with far fewer cells than the originally sampled, though technical saturation of rare populations such as Cajal-Retzius cells is not achieved. We additionally validated these findings with a recently published atlas of cell types across mouse organs and again find using subsampling that a much smaller number of cells recapitulates the cluster distinctions of the complete dataset.ConclusionsTogether, these findings suggest that most of the biologically interpretable cell types from the 1.3 million cell database can be recapitulated by analyzing 50,000 randomly selected cells, indicating that instead of profiling few individuals at high “cellular coverage,” cell atlas studies may instead benefit from profiling more individuals, or many time points at lower cellular coverage and then further enriching for populations of interest. This strategy is ideal for scenarios where cost and time are limited, though extremely rare populations of interest (< 1%) may be identifiable only with much higher cell numbers.

[1]  Aleksandra A. Kolodziejczyk,et al.  Classification of low quality cells from single-cell RNA-seq data , 2016, Genome Biology.

[2]  Vilas Menon,et al.  Clustering single cells: a review of approaches on high-and low-depth single-cell RNA-seq data. , 2018, Briefings in functional genomics.

[3]  S. Orkin,et al.  Mapping the Mouse Cell Atlas by Microwell-Seq , 2018, Cell.

[4]  Christof Koch,et al.  Adult Mouse Cortical Cell Taxonomy by Single Cell Transcriptomics , 2016, Nature Neuroscience.

[5]  A. Pierani,et al.  A Novel Role for Dbx1-Derived Cajal-Retzius Cells in Early Regionalization of the Cerebral Cortical Neuroepithelium , 2010, PLoS biology.

[6]  Vilas Menon,et al.  Clustering single cells: a review of approaches on high-and low-depth single-cell RNA-seq data. , 2018, Briefings in functional genomics.

[7]  R. Irizarry,et al.  Missing data and technical variability in single‐cell RNA‐sequencing experiments , 2018, Biostatistics.

[8]  Ian R. Wickersham,et al.  The BRAIN Initiative Cell Census Consortium: Lessons Learned toward Generating a Comprehensive Brain Cell Atlas , 2017, Neuron.

[9]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

[10]  M. Cugmas,et al.  On comparing partitions , 2015 .

[11]  S. Orkin,et al.  Mapping the Mouse Cell Atlas by Microwell-Seq , 2018, Cell.

[12]  Luca Muzio,et al.  Foxg1 Confines Cajal-Retzius Neuronogenesis and Hippocampal Morphogenesis to the Dorsomedial Pallium , 2005, The Journal of Neuroscience.

[13]  Sara Ballouz,et al.  Characterizing the replicability of cell types defined by single cell RNA-sequencing data using MetaNeighbor , 2018, Nature Communications.

[14]  Michael J. T. Stubbington,et al.  The Human Cell Atlas: from vision to reality , 2017, Nature.

[15]  Grace X. Y. Zheng,et al.  Massively parallel digital transcriptional profiling of single cells , 2016, Nature Communications.

[16]  C. Englund,et al.  Cajal-Retzius cells in the mouse: transcription factors, neurotransmitters, and birthdays suggest a pallial origin. , 2003, Brain research. Developmental brain research.

[17]  Sébastien Vigneau,et al.  Multiple origins of Cajal-Retzius cells at the borders of the developing pallium , 2005, Nature Neuroscience.

[18]  M. Schaub,et al.  SC3 - consensus clustering of single-cell RNA-Seq data , 2016, Nature Methods.

[19]  Mauricio Barahona,et al.  SC3 - consensus clustering of single-cell RNA-Seq data , 2016, Nature Methods.

[20]  T. Curran,et al.  A protein related to extracellular matrix proteins deleted in the mouse mutant reeler , 1995, Nature.

[21]  Evan Z. Macosko,et al.  Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics , 2016, Cell.

[22]  Allon M. Klein,et al.  Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells , 2015, Cell.

[23]  L. Hubert,et al.  Comparing partitions , 1985 .

[24]  S. Rétaux,et al.  Differential expression of LIM-homeodomain factors in Cajal-Retzius cells of primates, rodents, and birds. , 2010, Cerebral cortex.