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.

Time-Series Multi-Level Probabilistic Graphical Model for Representing Lineages over Uncertain Data

Representing and Processing Lineages over Uncertain Data Based on the Bayesian Network

Uncertain Knowledge Representation And Inference For Lineage Processing Over Uncertain Data

A Working Model for Uncertain Data with Lineage

Probabilistic Object Deputy Model for Uncertain Data and Lineage Management

Attribute Level Lineage in Uncertain Data with Dependencies

Bayesian multiple Gaussian graphical models for multilevel variables from unknown classes

Learning uncertain knowledge from uncertain data

Bayesian graphical models for computational network biology

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

Recording How-Provenance on Probabilistic Databases.

Bayesian Structure Learning in Multi-layered Genomic Networks

A Probabilistic Integrated Gene Network Model of Multiple Hematopoietic Lineage Specification

Qualitative Representation and Fusion of Probabilistic Causalities in Time-Series Environments

FamLink2 - A comprehensive tool for likelihood computations in pedigrees analyses involving linked DNA markers accounting for genotype uncertainties

Probabilistic logic methods and some applications to biology and medicine

Semantically-linked Bayesian Network: a Framework for Probabilistic Inference over Multiple Bayesian Networks

Leveraging graphical model techniques to study evolution on phylogenetic networks

Probabilistic Graphical Model Representation in Phylogenetics

Probabilistic Boolean networks: a rule-based uncertainty model for gene regulatory networks