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 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 allow 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 necessary for efficient use of resources. Our analysis emphasizes the importance of selecting a protocol based on the biological questions and features of interest.

[1]  R. Reiter,et al.  Analysis of the daily changes of melatonin receptors in the rat liver , 2013, Journal of pineal research.

[2]  K. Bae,et al.  Identification of transcriptional regulatory elements required for the Mup2 expression in circadian clock mutant mice. , 2011, Biochemical and biophysical research communications.

[3]  Evan Z. Macosko,et al.  Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution , 2019, Science.

[4]  Davis J. McCarthy,et al.  A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor , 2016, F1000Research.

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

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

[7]  Céline Jouffe,et al.  Perturbed rhythmic activation of signaling pathways in mice deficient for Sterol Carrier Protein 2-dependent diurnal lipid transport and metabolism , 2016, Scientific Reports.

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

[9]  T. Baumert,et al.  Single-cell genomics and spatial transcriptomics: Discovery of novel cell states and cellular interactions in liver physiology and disease biology , 2020, Journal of hepatology.

[10]  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.

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

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

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

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

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

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

[17]  Richard Bonneau,et al.  High-definition spatial transcriptomics for in situ tissue profiling , 2019, Nature Methods.

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

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

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

[21]  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.

[22]  F S Fay,et al.  Visualization of single RNA transcripts in situ. , 1998, Science.

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

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

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

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

[27]  Albert K Groen,et al.  Complex interaction between circadian rhythm and diet on bile acid homeostasis in male rats , 2017, Chronobiology international.

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

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

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

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

[32]  Patrik L. Ståhl,et al.  Visualization and analysis of gene expression in tissue sections by spatial transcriptomics , 2016, Science.

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

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

[35]  Hannah A. Pliner,et al.  Reversed graph embedding resolves complex single-cell trajectories , 2017, Nature Methods.

[36]  D. McSwiggen,et al.  Single Molecule Fluorescence In Situ Hybridization (smFISH) Analysis in Budding Yeast Vegetative Growth and Meiosis , 2018, Journal of visualized experiments : JoVE.

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

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

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

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

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

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

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

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

[45]  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.

[46]  M. Lazar,et al.  The hepatic circadian clock fine-tunes the lipogenic response to feeding through RORα/γ , 2017, Genes & development.

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

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

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

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

[51]  R. Gebhardt,et al.  Liver zonation: Novel aspects of its regulation and its impact on homeostasis. , 2014, World journal of gastroenterology.

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

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