Abstract:Among existing computational algorithms for single-cell RNA-seq analysis, clustering and trajectory inference are two major types of analysis that are routinely applied. For a given dataset, clustering and trajectory inference can generate vastly different visualizations that lead to very different interpretations of the data. To address this issue, we propose multiple scores to quantify the "clusterness" and "trajectoriness" of single-cell RNA-seq data, in other words, whether the data looks like a collection of distinct clusters or a continuum of progression trajectory. The scores we introduce are based on pairwise distance distribution, persistent homology, vector magnitude, Ripley's K, and degrees of connectivity. Using simulated datasets, we demonstrate that the proposed scores are able to effectively differentiate between cluster-like data and trajectory-like data. Using real single-cell RNA-seq datasets, we demonstrate the scores can serve as indicators of whether clustering analysis or trajectory inference is a more appropriate choice for biological interpretation of the data. Single-cell sequencing technologies have motivated development of numerous computational algorithms. Two main types of these algorithms are clustering and trajectory inference. When scientists have a scRNA-seq dataset, they usually pick one of these approaches based on what they think the data shows. If they think the data has distinct clusters of cells, they will analyze the data using clustering algorithms. If they think the data shows a continuous progression, they will use trajectory inference algorithms. However, sometimes using clustering and trajectory inference on the same data can lead to very different interpretations, where clustering algorithms produce distinct cell clusters while trajectory inference on the same data show continuous trajectories. This makes us wonder: which way of looking at the data is more appropriate? In this paper, we developed a pipeline for quantifying the "clusterness" and "trajectoriness" of scRNA-seq data, in other words, whether the data looks like a collection of distinct clusters or a continuum of progression trajectory. We think such geometric quantification is an important question that should be broadly discussed in the single-cell research community.

Unraveling the timeline of gene expression: A pseudotemporal trajectory analysis of single-cell RNA sequencing data

Gene trajectory inference for single-cell data by optimal transport metrics

Scart: Recognizing Cell Clusters and Constructing Trajectory from Single-Cell Epigenomic Data

Inference of Differentiation Time for Single Cell Transcriptomes Using Cell Population Reference Data

Beyond pseudotime: Following T-cell maturation in single-cell RNAseq time series

Concepts and Limitations for Learning Developmental Trajectories from Single Cell Genomics.

A statistical framework for differential pseudotime analysis with multiple single-cell RNA-seq samples

A descriptive marker gene approach to single-cell pseudotime inference

SCOUT: A New Algorithm for the Inference of Pseudo-Time Trajectory Using Single-Cell Data

scInTime: A Computational Method Leveraging Single-Cell Trajectory and Gene Regulatory Networks to Identify Master Regulators of Cellular Differentiation

Trajectory-based differential expression analysis for single-cell sequencing data

Modeling gene expression cascades during cell state transitions

Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics

Temporal modelling using single-cell transcriptomics

Time-resolved Single-cell Transcriptomic Sequencing

Inferring Population Dynamics from Single-Cell RNA-sequencing Time Series Data.

BLTSA: Pseudotime Prediction for Single Cells by Branched Local Tangent Space Alignment

Quantifying the clusterness and trajectoriness of single-cell RNA-seq data

Tempora: cell trajectory inference using time-series single-cell RNA sequencing data

Trajectory inference from single-cell genomics data with a process time model