Model-Based Trajectory Inference for Single-Cell RNA Sequencing Using Deep Learning with a Mixture Prior
Published:
Background
In the analysis of single-cell RNA sequencing (scRNA-seq) data, trajectory inference or pseudotime analysis methods are proposed to infer the cell developmental lineages and order cells on their pseudotime within these lineages. As scRNA-seq datasets continue to expand, there is a strong need to integrate different datasets and perform trajectory inference jointly to align cells experiencing similar dynamic changes but are from different labs or under various conditions. By a joint cell trajectory inference of our own data with other publicly available scRNA-seq datasets, we aim to obtain a complete picture of the neurogenesis in the mouse neocortex and expand our understanding of the progenitor heterogeneity.
Results
We develop a new quantitative framework to perform a joint cell lineage analysis, which uses a statistical hierarchical mixture model for the trajectory structure on the cells' low-dimensional embedding. Our method outperforms other state-of-the-art trajectory inference methods on both real and synthetic data under various trajectory topologies. We also apply it to jointly analyze two single-cell RNA sequencing datasets on the mouse neocortex, with a conditional variational auto-encoder and adjusting for batch and cell-cycle effects. The inferred developmental trajectory enables us to study the progenitor cell types and identify the subtle development of IPCs. After correcting post-estimation bias, we have detected around 1000 genes whose co-expression patterns with Eomes change significantly across embryonic days.
Conclusions
(1) Our method successfully merges cells from different datasets and learns a shared trajectory of neurogenesis, keeping biologically meaningful differences between the two datasets. (2) We identify progenitor cell types based on integrated data and can further divide the IPCs into post-mitotic and mitotic IPCs. (3) We are able to identify temporally expressed transcription regulators in cortical progenitors.
