Abstract:Background Single-cell RNA sequencing (scRNA-seq) technologies have rapidly developed in recent years. The droplet-based single cell platforms enable the profiling of gene expression in tens of thousands of cells per sample. The goal of a typical scRNA-seq analysis is to identify different cell subpopulations and their respective marker genes. Additionally, trajectory analysis can be used to infer the developmental or differentiation trajectories of cells. Methods This article demonstrates a comprehensive workflow for performing trajectory inference and time course analysis on a multi-sample single-cell RNA-seq experiment of the mouse mammary gland. The workflow uses open-source R software packages and covers all steps of the analysis pipeline, including quality control, doublet prediction, normalization, integration, dimension reduction, cell clustering, trajectory inference, and pseudo-bulk time course analysis. Sample integration and cell clustering follows the Seurat pipeline while the trajectory inference is conducted using the monocle3 package. The pseudo-bulk time course analysis uses the quasi-likelihood framework of edgeR. Results Cells are ordered and positioned along a pseudotime trajectory that represented a biological process of cell differentiation and development. The study successfully identified genes that were significantly associated with pseudotime in the mouse mammary gland. Conclusions The demonstrated workflow provides a valuable resource for researchers conducting scRNA-seq analysis using open-source software packages. The study successfully demonstrated the usefulness of trajectory analysis for understanding the developmental or differentiation trajectories of cells. This analysis can be applied to various biological processes such as cell development or disease progression, and can help identify potential biomarkers or therapeutic targets.

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

PseudotimeDE: inference of differential gene expression along cell pseudotime with well-calibrated p-values from single-cell RNA sequencing data

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

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

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

Accounting for technical noise in differential expression analysis of single-cell RNA sequencing data

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

Identify Differential Genes and Cell Subclusters from Time-Series Scrna-Seq Data Using Sctitans.

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

A descriptive marker gene approach to single-cell pseudotime inference

Pseudotime analysis for time-series single-cell sequencing and imaging data

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

DTWscore: Differential Expression and Cell Clustering Analysis for Time-Series Single-Cell RNA-seq Data

Temporal modelling using single-cell transcriptomics

The Lomb-Scargle periodogram-based differentially expressed gene detection along pseudotime

CASi: A Framework for Cross-Timepoint Analysis of Single-Cell RNA Sequencing Data

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

Two-phase differential expression analysis for single cell RNA-seq.

Multi-scale Poisson process approaches for differential expression analysis of high-throughput sequencing data

CASi: A Multi-Timepoint Scrnaseq Data Analysis Framework

CCPE: cell cycle pseudotime estimation for single cell RNA-seq data