Trajectory inference across multiple conditions with condiments: differential topology, progression, differentiation, and expression.

In single-cell RNA-sequencing (scRNA-seq), gene expression is assessed individually for each cell, allowing the investigation of developmental processes, such as embryogenesis and cellular differentiation and regeneration, at unprecedented resolutions. In such dynamic biological systems, grouping cells into discrete groups is not reflective of the biology. Cellular states rather form a continuum, e.g., for the differentiation of stem cells into mature cell types. This process is often represented via a trajectory in a reduced-dimensional representation of the scRNA-seq dataset. While many methods have been suggested for trajectory inference, it is often unclear how to handle multiple biological groups or conditions, e.g., inferring and comparing the differentiation trajectories of wild-type and knock-out stem cell populations. In this manuscript, we present a method for the estimation and downstream interpretation of cell trajectories across multiple conditions. Our framework allows the interpretation of differences between conditions at the trajectory, cell population, and gene expression levels. We start by integrating datasets from multiple conditions into a single trajectory. By comparing the conditions along the trajectory's path, we can detect large-scale changes, indicative of differential progression. We also demonstrate how to detect subtler changes by finding genes that exhibit different behaviors between these conditions along a differentiation path.

[1]  David van Dijk,et al.  Enhancing experimental signals in single-cell RNA-sequencing data using graph signal processing , 2019 .

[2]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[3]  Bernhard Schölkopf,et al.  A Kernel Two-Sample Test , 2012, J. Mach. Learn. Res..

[4]  Paul J. Hoffman,et al.  Comprehensive Integration of Single-Cell Data , 2018, Cell.

[5]  Lior Pachter,et al.  A curated database reveals trends in single-cell transcriptomics , 2020, Database : the journal of biological databases and curation.

[6]  R. Satija,et al.  Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression , 2019, Genome Biology.

[7]  Yvan Saeys,et al.  Trajectory-based differential expression analysis for single-cell sequencing data , 2020, Nature Communications.

[8]  G. Nkiet,et al.  $k$-Sample problem based on generalized maximum mean discrepancy , 2018, 1811.09103.

[9]  K. Bui,et al.  The inner junction complex of the cilia is an interaction hub that involves tubulin post-translational modifications , 2020, eLife.

[10]  Leland McInnes,et al.  UMAP: Uniform Manifold Approximation and Projection for Dimension Reduction , 2018, ArXiv.

[11]  Xiuyuan Cheng,et al.  Detection of differentially abundant cell subpopulations discriminates biological states in scRNA-seq data , 2019, bioRxiv.

[12]  J. Kiefer,et al.  K-Sample Analogues of the Kolmogorov-Smirnov and Cramer-V. Mises Tests , 1959 .

[13]  T. Zacharewski,et al.  Single-Nuclei RNA Sequencing Assessment of the Hepatic Effects of 2,3,7,8-Tetrachlorodibenzo-p-dioxin , 2020, Cellular and molecular gastroenterology and hepatology.

[14]  Gennady Korotkevich,et al.  Fast gene set enrichment analysis , 2021 .

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

[16]  Leland McInnes,et al.  UMAP: Uniform Manifold Approximation and Projection , 2018, J. Open Source Softw..

[17]  Russell B. Fletcher,et al.  Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics , 2017, BMC Genomics.

[18]  SMARCB1 loss interacts with neuronal differentiation state to block maturation and impact cell stability. , 2020, Genes & development.

[19]  S. Stouffer Adjustment during army life , 1977 .

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

[21]  David Lopez-Paz,et al.  Revisiting Classifier Two-Sample Tests , 2016, ICLR.

[22]  Davide Risso,et al.  Rapid non-uniform adaptation to conformation-specific KRAS(G12C) inhibition , 2020, Nature.

[23]  Alexey Sergushichev,et al.  An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation , 2016 .

[24]  John C. Marioni,et al.  Milo: differential abundance testing on single-cell data using k-NN graphs , 2020, bioRxiv.

[25]  John C. Marioni,et al.  Testing for differential abundance in mass cytometry data , 2017, Nature Methods.

[26]  C. Ponting,et al.  Resolving the fibrotic niche of human liver cirrhosis at single cell level , 2019, Nature.

[27]  Mark S. Anderson,et al.  Combined transient ablation and single cell RNA sequencing reveals the development of medullary thymic epithelial cells , 2020, bioRxiv.

[28]  Yao-fei Pei,et al.  Silencing of LAMC2 Reverses Epithelial-Mesenchymal Transition and Inhibits Angiogenesis in Cholangiocarcinoma via Inactivation of the Epidermal Growth Factor Receptor Signaling Pathway. , 2019, The American journal of pathology.

[29]  Cole Trapnell,et al.  A pooled single-cell genetic screen identifies regulatory checkpoints in the continuum of the epithelial-to-mesenchymal transition , 2019, Nature Genetics.

[30]  Xiuyuan Cheng,et al.  Detection of differentially abundant cell subpopulations discriminates biological states in scRNA-seq data , 2019, bioRxiv.

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

[32]  R. Dobrushin Prescribing a System of Random Variables by Conditional Distributions , 1970 .

[33]  Yvan Saeys,et al.  dyngen: a multi-modal simulator for spearheading new single-cell omics analyses , 2020, bioRxiv.

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

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

[36]  Lai Guan Ng,et al.  Dimensionality reduction for visualizing single-cell data using UMAP , 2018, Nature Biotechnology.