Reproducibility across single-cell RNA-seq protocols for spatial ordering analysis

As newer single-cell protocols generate increasingly more cells at reduced sequencing depths, the value of a higher read depth may be overlooked. Using data from three different single-cell RNA-seq protocols that lend themselves to having either higher read depth (Smart-seq) or many cells (MARS-seq and 10X), we evaluate their ability to recapitulate biological signals in the context of pseudo-spatial reconstruction. Overall, we find gene expression profiles after spatial-reconstruction analysis are highly reproducible between datasets despite being generated by different protocols and using different computational algorithms. While UMI based protocols such as 10X and MARS-seq allow for capturing more cells, Smart-seq’s higher sensitivity and read-depth allows for analysis of lower expressed genes and isoforms. Additionally, we evaluate trade-offs for each protocol by performing subsampling analyses, and find that optimizing the balance between sequencing depth and number of cells within a protocol is important for efficient use of resources. Our analysis emphasizes the importance of selecting a protocol based on the biological questions and features of interest.

[1]  Christoph Ziegenhain,et al.  The impact of amplification on differential expression analyses by RNA-seq , 2016, Scientific Reports.

[2]  A. Oshlack,et al.  Gene length and detection bias in single cell RNA sequencing protocols , 2017, bioRxiv.

[3]  Gene W. Yeo,et al.  Single-Cell Alternative Splicing Analysis with Expedition Reveals Splicing Dynamics during Neuron Differentiation. , 2017, Molecular cell.

[4]  P. Moldéus,et al.  Zonation of acetaminophen metabolism and cytochrome P450 2E1-mediated toxicity studied in isolated periportal and perivenous hepatocytes. , 1993, Biochemical pharmacology.

[5]  R. Gebhardt,et al.  Different capacities for amino acid transport in periportal and perivenous hepatocytes isolated by digitonin/collagenase perfusion , 1989, Hepatology.

[6]  M. Newton,et al.  SCnorm: robust normalization of single-cell RNA-seq data , 2017, Nature Methods.

[7]  M. Guzmán,et al.  Zonation of fatty acid metabolism in rat liver. , 1989, The Biochemical journal.

[8]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

[9]  Hongkai Ji,et al.  TSCAN: Pseudo-time reconstruction and evaluation in single-cell RNA-seq analysis , 2016, Nucleic acids research.

[10]  Samriddha Ray,et al.  Single cell RNA analysis identifies cellular heterogeneity and adaptive responses of the lung at birth , 2019, Nature Communications.

[11]  I. Hellmann,et al.  Comparative Analysis of Single-Cell RNA Sequencing Methods , 2016, bioRxiv.

[12]  Ahmed Mahfouz,et al.  Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells , 2018, Nature Biotechnology.

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

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

[15]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[16]  Ana Conesa,et al.  Single-cell RNAseq for the study of isoforms—how is that possible? , 2018, Genome Biology.

[17]  Eric T. Wang,et al.  Alternative Isoform Regulation in Human Tissue Transcriptomes , 2008, Nature.

[18]  R. Stewart,et al.  Single-cell RNA-seq reveals novel regulators of human embryonic stem cell differentiation to definitive endoderm , 2016, Genome Biology.

[19]  J. Reddy,et al.  Biochemical characterization of two functional human liver acyl-CoA oxidase isoforms 1a and 1b encoded by a single gene. , 2007, Biochemical and biophysical research communications.

[20]  James T. Webber,et al.  Single-cell transcriptomics of 20 mouse organs creates a Tabula Muris , 2018, Nature.

[21]  R. Sandberg,et al.  Full-Length mRNA-Seq from single cell levels of RNA and individual circulating tumor cells , 2012, Nature Biotechnology.

[22]  M. Bonin,et al.  Differential gene expression in periportal and perivenous mouse hepatocytes , 2006, The FEBS journal.

[23]  S. Linnarsson,et al.  Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. , 2011, Genome research.

[24]  Yvan Saeys,et al.  A comparison of single-cell trajectory inference methods , 2019, Nature Biotechnology.

[25]  Cole Trapnell,et al.  The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells , 2014, Nature Biotechnology.

[26]  Sten Linnarsson,et al.  Alternative TSSs are co‐regulated in single cells in the mouse brain , 2017, Molecular systems biology.

[27]  D. Tosh,et al.  Glucagon regulation of gluconeogenesis and ketogenesis in periportal and perivenous rat hepatocytes. Heterogeneity of hormone action and of the mitochondrial redox state. , 1988, The Biochemical journal.

[28]  I. Amit,et al.  Single-cell spatial reconstruction reveals global division of labor in the mammalian liver , 2016, Nature.

[29]  Ning Leng,et al.  Oscope identifies oscillatory genes in unsynchronized single cell RNA-seq experiments , 2015, Nature Methods.

[30]  Shu-Whei Tsai,et al.  Hepatic translation control in the late-gestation fetal rat. , 2008, American journal of physiology. Regulatory, integrative and comparative physiology.

[31]  K. Lindros,et al.  Urea synthesis in freshly isolated and in cultured periportal and perivenous hepatocytes. , 1986, The Biochemical journal.

[32]  S. Bhatia,et al.  Zonal liver cell heterogeneity : effects of oxygen on metabolic functions of hepatocytes , 2011 .

[33]  Gioele La Manno,et al.  Quantitative single-cell RNA-seq with unique molecular identifiers , 2013, Nature Methods.

[34]  Valentine Svensson,et al.  Power Analysis of Single Cell RNA-Sequencing Experiments , 2016, Nature Methods.

[35]  Principal Investigators,et al.  Single-cell transcriptomics of 20 mouse organs creates a Tabula Muris , 2018 .

[36]  Aleksandra A. Kolodziejczyk,et al.  Accounting for technical noise in single-cell RNA-seq experiments , 2013, Nature Methods.

[37]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[38]  Catalin C. Barbacioru,et al.  mRNA-Seq whole-transcriptome analysis of a single cell , 2009, Nature Methods.

[39]  S. Teichmann,et al.  A practical guide to single-cell RNA-sequencing for biomedical research and clinical applications , 2017, Genome Medicine.

[40]  R. Gebhardt,et al.  Heterogeneous distribution of glutamine synthetase among rat liver parenchymal cells in situ and in primary culture. , 1983, The EMBO journal.

[41]  T. Kietzmann Metabolic zonation of the liver: The oxygen gradient revisited , 2017, Redox biology.