Abstract:Abstract Motivation Eukaryotic cells contain organelles called mitochondria that have their own genome. Most cells contain thousands of mitochondria which replicate, even in nondividing cells, by means of a relatively error-prone process resulting in somatic mutations in their genome. Because of the higher mutation rate compared to the nuclear genome, mitochondrial mutations have been used to track cellular lineage, particularly using single-cell sequencing that measures mitochondrial mutations in individual cells. However, existing methods to infer the cell lineage tree from mitochondrial mutations do not model “heteroplasmy,” which is the presence of multiple mitochondrial clones with distinct sets of mutations in an individual cell. Single-cell sequencing data thus provide a mixture of the mitochondrial clones in individual cells, with the ancestral relationships between these clones described by a mitochondrial clone tree. While deconvolution of somatic mutations from a mixture of evolutionarily related genomes has been extensively studied in the context of bulk sequencing of cancer tumor samples, the problem of mitochondrial deconvolution has the additional constraint that the mitochondrial clone tree must be concordant with the cell lineage tree. Results We formalize the problem of inferring a concordant pair of a mitochondrial clone tree and a cell lineage tree from single-cell sequencing data as the Nested Perfect Phylogeny Mixture (NPPM) problem. We derive a combinatorial characterization of the solutions to the NPPM problem, and formulate an algorithm, MERLIN, to solve this problem exactly using a mixed integer linear program. We show on simulated data that MERLIN outperforms existing methods that do not model mitochondrial heteroplasmy nor the concordance between the mitochondrial clone tree and the cell lineage tree. We use MERLIN to analyze single-cell whole-genome sequencing data of 5220 cells of a gastric cancer cell line and show that MERLIN infers a more biologically plausible cell lineage tree and mitochondrial clone tree compared to existing methods. Availability and implementation https://github.com/raphael-group/MERLIN.

Robust Lineage Reconstruction from High-Dimensional Single-Cell Data

Scart: Recognizing Cell Clusters and Constructing Trajectory from Single-Cell Epigenomic Data

LinRace: cell division history reconstruction of single cells using paired lineage barcode and gene expression data

Gene-expression memory-based prediction of cell lineages from scRNA-seq datasets

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

Deep Lineage: Single-Cell Lineage Tracing and Fate Inference Using Deep Learning

CIARA: a Cluster-Independent Algorithm for Identifying Markers of Rare Cell Types from Single-Cell Sequencing Data

scTrace+: enhance the cell fate inference by integrating the lineage-tracing and multi-faceted transcriptomic similarity information

CellRank for Directed Single-Cell Fate Mapping

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

Mapping Lineage-Traced Cells Across Time Points with Moslin

Bifurcation Analysis of Single-Cell Gene Expression Data Reveals Epigenetic Landscape.

LineageVAE: Reconstructing Historical Cell States and Transcriptomes toward Unobserved Progenitors

2015 – A MOUSE MODEL WITH HIGH CLONAL BARCODE DIVERSITY FOR JOINT LINEAGE, TRANSCRIPTOMIC, AND EPIGENOMIC PROFILING IN SINGLE CELLS

Reinforcement learning guides single-cell sequencing in decoding lineage and cell fate decisions

Fast, accurate reconstruction of cell lineages from large-scale fluorescence microscopy data

Joint inference of cell lineage and mitochondrial evolution from single-cell sequencing data

Inferring cell differentiation maps from lineage tracing data

Tree reconstruction guarantees from CRISPR-Cas9 lineage tracing data using Neighbor-Joining

ScisTree2: An Improved Method for Large-scale Inference of Cell Lineage Trees and Genotype Calling from Noisy Single Cell Data

Single-cell lineage tracing by integrating CRISPR-Cas9 mutations with transcriptomic data