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

**Authors:** Thinh N. Tran, Gary D. Bader, Thomas Lengauer, Qing Nie, Thomas Lengauer, Qing Nie, Thomas Lengauer, Qing Nie, Thomas Lengauer, Qing Nie

PMC · DOI: 10.1371/journal.pcbi.1008205 · PLoS Computational Biology · 2020-09-09

## TL;DR

Tempora is a new method that uses time-series data and biological pathways to better understand how cells change over time in processes like tissue development.

## Contribution

Tempora introduces a novel approach that combines time information and pathway data to improve cell trajectory inference from time-series scRNA-seq data.

## Key findings

- Tempora outperformed existing methods in accuracy and speed when reconstructing known developmental trajectories.
- Using cluster-level analysis and pathway information improves signal amplification and reduces computational demands.
- The method enhances interpretability by identifying time-dependent pathways and known cell types.

## Abstract

Single-cell RNA sequencing (scRNA-seq) can map cell types, states and transitions during dynamic biological processes such as tissue development and regeneration. Many trajectory inference methods have been developed to order cells by their progression through a dynamic process. However, when time series data is available, most of these methods do not consider the available time information when ordering cells and are instead designed to work only on a single scRNA-seq data snapshot. We present Tempora, a novel cell trajectory inference method that orders cells using time information from time-series scRNA-seq data. In performance comparison tests, Tempora inferred known developmental lineages from three diverse tissue development time series data sets, beating state of the art methods in accuracy and speed. Tempora works at the level of cell clusters (types) and uses biological pathway information to help identify cell type relationships. This approach increases gene expression signal from single cells, processing speed, and interpretability of the inferred trajectory. Our results demonstrate the utility of a combination of time and pathway information to supervise trajectory inference for scRNA-seq based analysis.

Single-cell RNA sequencing (scRNA-seq) enables an unparalleled ability to map the heterogeneity of dynamic multicellular processes, such as tissue development, tumor growth, wound response and repair, and inflammation. Multiple methods have been developed to order cells along a pseudotime axis that represents a trajectory through such processes using the concept that cells that are closely related in a lineage will have similar transcriptomes. However, time series experiments provide another useful information source to order cells, from earlier to later time point. By introducing a novel use of biological pathway prior information, our Tempora algorithm improves the accuracy and speed of cell trajectory inference from time-series scRNA-seq data as measured by reconstructing known developmental trajectories from three diverse data sets. By analyzing scRNA-seq data at the cluster (cell type) level instead of at the single-cell level and by using known pathway information, Tempora amplifies gene expression signals from one cell using similar cells in a cluster and similar genes within a pathway. This approach also reduces computational time and resources needed to analyze large data sets because it works with a relatively small number of clusters instead of a potentially large number of cells. Finally, it eases interpretation, via operating on a relatively small number of clusters which usually represent known cell types, as well as by identifying time-dependent pathways. Tempora is useful for finding novel insights in dynamic processes.

## Full-text entities

- **Genes:** Cldn5 (claudin 5) [NCBI Gene 12741] {aka MBEC1, Tmvcf}, Prdx1 (peroxiredoxin 1) [NCBI Gene 18477] {aka MSP23, NkefA, OSF-3, OSF3, PAG, Paga}, Dlx1 (distal-less homeobox 1) [NCBI Gene 13390] {aka DII B, Dlx, Dlx-1}, Myog (myogenin) [NCBI Gene 17928] {aka MYF4, bHLHc3, myo}, Sox2 (SRY (sex determining region Y)-box 2) [NCBI Gene 20674] {aka Sox-2, lcc, ysb}, Nes (nestin) [NCBI Gene 18008] {aka ESTM46, Ifaprc2, Marc2, RC2}, Olig3 (oligodendrocyte transcription factor 3) [NCBI Gene 94222] {aka Bhlhb7, bHLHe20, oligo3}, Cxcl12 (C-X-C motif chemokine ligand 12) [NCBI Gene 20315] {aka Pbsf, Scyb12, Sdf1, Tlsf, Tpar1}, Cdk1 (cyclin dependent kinase 1) [NCBI Gene 12534] {aka Cdc2, Cdc2a, p34<CDC2>}, Prrx2 (paired related homeobox 2) [NCBI Gene 20204] {aka Prx2, S8}, Sphk1 (sphingosine kinase 1) [NCBI Gene 20698] {aka 1110006G24Rik, Sk1, Spk1}, Mki67 (antigen identified by monoclonal antibody Ki 67) [NCBI Gene 17345] {aka D630048A14Rik, Ki-67, Ki67}, Rgs5 (regulator of G-protein signaling 5) [NCBI Gene 19737] {aka 1110070A02Rik}, MYOG (myogenin) [NCBI Gene 4656] {aka MYF4, bHLHc3, myf-4}, Aif1 (allograft inflammatory factor 1) [NCBI Gene 11629] {aka AIF-1, D17H6S50E, G1, Iba1}
- **Diseases:** PEA (MESH:D058606), cancer (MESH:D009369), inflammation (MESH:D007249), HSMM (MESH:D005207), AP/RP (MESH:D012174)
- **Chemicals:** Monocle (-), GABA (MESH:D005680)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090], Danio rerio (leopard danio, species) [taxon 7955]
- **Mutations:** S10A, S11A, Pro with 2, S11D

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC7505465/full.md

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC7505465/full.md

## References

55 references — full list in the complete paper: https://tomesphere.com/paper/PMC7505465/full.md

---
Source: https://tomesphere.com/paper/PMC7505465