Reconstructing Ordered Ontogenic Trajectories

A set of tools to reconstruct ordered ontogenic trajectories from single cell RNAseq data.

Wajid Jawaid 2017-07-09

Project Status: Active - The project has reached a stable, usable state and is being actively developed.

Reconstructing ordered ontogenic trajectories provides methods for:

  1. Identifying informative genes (crude),
  2. Finding a sparse graph structure between closely related cells by removing spurious edges and
  3. Traversing the graph in a biologiclal informed way using a Directed Non-markovian Monte-Carlo method.

The main goal of roots is to infer plausible developmental journeys guided by the user.




Here I take the mouse adult haematopoietic data from Nestorowa et al.. Data is downloaded and processed using the goggles() function as below.

blood <- read.table("",
                    sep = " ")
cellNames <- read.table("",
                        sep = " ", stringsAsFactors = FALSE)[,1]
rownames(blood) <- gsub("LT\\.", "LT-", cellNames)
geneNames <- read.table("",
                        sep = " ", stringsAsFactors = FALSE)[,1]
colnames(blood) <- geneNames
blood <- as.matrix(blood)
rm(cellNames, geneNames)
## Load metadata
meta <- read.csv("")
colnames(meta) <- c("cellType", "index", "name")
rownames(meta) <- meta$name
meta$col <- bglab::ggCol(meta$cellType)
nmeta <- data.frame(col=rep("#00000011", nrow(blood)), stringsAsFactors = FALSE,
                    row.names = rownames(blood))
nmeta[rownames(meta),"col"] <- meta$col
leg <- data.frame(cell=as.character(unique(meta$cellType)),
                  col=as.character(unique(meta$col)), stringsAsFactors = FALSE)
legOrd <- c(5, 8, 6, 7, 1, 4, 2, 3)
## Analyse
xx <- goggles(blood)
## Plot
plot(xx$l, pch=16, col = nmeta[rownames(xx$l), "col"])
legend("topright", legend = leg$cell[legOrd], fill=leg$col[legOrd], inset=0.02)

Output from goggle() function


Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.


1.0 by Wajid Jawaid, 2 years ago

Browse source code at

Authors: Wajid Jawaid [aut, cre]

Documentation:   PDF Manual  

GPL (>= 2) license

Imports animation, rARPACK, igraph

Suggests knitr, rmarkdown

See at CRAN