Abstract:Recent dynamic lineage tracing technologies combine CRISPR-based genome editing with single-cell sequencing to track cell divisions during development. A key computational problem in dynamic lineage tracing is to infer a cell lineage tree from the measured CRISPR-induced mutations. Three features of dynamic lineage tracing data distinguish this problem from standard phylogenetic tree inference. First, the CRISPR-editing process modifies a genomic location exactly once. This property is not well described by the time-reversible models commonly used in phylogenetics. Second, as a consequence of non-modifiability, the number of mutations per time unit decreases over time. Third, CRISPR-based genome-editing and single-cell sequencing results in high rates of both heritable and non-heritable (dropout) missing data. To model these features, we introduce the Probabilistic Mixed-type Missing (PMM) model. We describe an algorithm, LAML (Lineage Analysis via Maximum Likelihood), to search for the maximum likelihood (ML) tree under the PMM model. LAML combines an Expectation Maximization (EM) algorithm with a heuristic tree search to jointly estimate tree topology, branch lengths and missing data parameters. We derive a closed-form solution for the M-step in the case of no heritable missing data, and a block coordinate ascent approach in the general case which is more efficient than the standard General Time Reversible (GTR) phylogenetic model. On simulated data, LAML infers more accurate tree topologies and branch lengths than existing methods, with greater advantages on datasets with higher ratios of heritable to non-heritable missing data. We show that LAML provides unbiased estimates of branch lengths. In contrast, we demonstrate that maximum parsimony methods for lineage tracing data not only underestimate branch lengths, but also yield branch lengths which are not proportional to time, due to the nonlinear decay in the number of mutations on branches further from the root. On lineage tracing data from a mouse model of lung adenocarcinoma, we show that LAML infers phylogenetic distances that are more concordant with gene expression data compared to distances derived from maximum parsimony. The LAML tree topology is more plausible than existing published trees, with fewer total cell migrations between distant metastases and fewer reseeding events where cells migrate back to the primary tumor. Crucially, we identify three distinct time epochs of metastasis progression, which includes a burst of metastasis events to various anatomical sites during a single month.

Supplementary Material for: Diffusion pseudotime robustly reconstructs lineage branching

Diffusion pseudotime robustly reconstructs lineage branching

Model-Based Branching Point Detection In Single-Cell Data By K-Branches Clustering

Supplementary Figures and Text for “destiny – diffusion maps for single-cell time-course data”

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

Diffusion Maps for High-Dimensional Single-Cell Analysis of Differentiation Data

Maximum Likelihood Inference of Time-scaled Cell Lineage Trees with Mixed-type Missing Data

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

Detecting branching rate heterogeneity in multifurcating trees with applications in lineage tracing data

Cell cycle expression heterogeneity predicts degree of differentiation

diverse-seq: an application for alignment-free selecting and clustering biological sequences

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

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

Reconstructing probabilistic trees of cellular differentiation from single-cell RNA-seq data

Unveiling Clonal Cell Fate and Differentiation Dynamics: A Hybrid NeuralODE-Gillespie Approach

Reconstructing tumor clonal heterogeneity and evolutionary relationships based on tumor DNA sequencing data

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

Sciphy: A Bayesian phylogenetic framework using sequential genetic lineage tracing data.

Overcoming Expressional Drop-outs in Lineage Reconstruction from Single-Cell RNA-Sequencing Data

Branching topology of the human embryo transcriptome revealed by Entropy Sort Feature Weighting