Is it possible to reconstruct an accurate cell lineage using CRISPR recorders?

Simulation of CRISPR-Cas9 editing on evolving barcode and accuracy of lineage tracing

We designed a simulation program that mimics the CRISPR-Cas9 editing on evolving barcode and double strand break repair procedure along with cell divisions. Emerging barcode mutations tend to build upon previously existing mutations, occurring sequentially with each generation. This process results in a unique mutation profile in each cell. We sample the barcodes in leaf cells and reconstruct the lineage, comparing it to the original lineage tree to test algorithm accuracy under different parameter settings. Our computational simulations validate the reasonable assumptions deduced from experimental observations, emphasizing that factors such as sampling size, barcode length, multiple barcodes, indel probabilities, and Cas9 activity are critical for accurate and successful lineage tracing. Among the many factors we found that sampling size and indel probabilities are two major ones that affect lineage tracing accuracy. Large segment deletions in early generations could greatly impact lineage accuracy. These simulation results offer insightful recommendations for enhancing the design and analysis of Cas9-mediated molecular barcodes in actual experiments.

Barcoding Notch signaling in the developing brain

Developmental signaling inputs are fundamental for shaping cell fates and behavior. However, traditional fluorescent-based signaling reporters have limitations in scalability and molecular resolution of cell types. We present SABER-seq, a CRISPR-Cas molecular recorder that stores transient developmental signaling cues as permanent mutations in cellular genomes for deconstruction at later stages via single-cell transcriptomics. We applied SABER-seq to record Notch signaling in developing zebrafish brains. SABER-seq has two components: a signaling sensor and a barcode recorder. The sensor activates Cas9 in a Notch-dependent manner with inducible control while the recorder accumulates mutations that represent Notch activity in founder cells. We combine SABER-seq with an expanded juvenile brain atlas to define cell types whose fates are determined downstream of Notch signaling. We identified examples wherein Notch signaling may have differential impact on terminal cell fates. SABER-seq is a novel platform for rapid, scalable and high-resolution mapping of signaling activity during development.

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for studying cellular differentiation, but accurately tracking cell fate transitions can be challenging, especially in disease conditions. Here we introduce PhyloVelo, a computational framework that estimates the velocity of transcriptomic dynamics by using monotonically expressed genes (MEGs) or genes with expression patterns that either increase or decrease, but do not cycle, through phylogenetic time. Through integration of scRNA-seq data with lineage information, PhyloVelo identifies MEGs and reconstructs a transcriptomic velocity field. We validate PhyloVelo using simulated data and Caenorhabditis elegans ground truth data, successfully recovering linear, bifurcated and convergent differentiations. Applying PhyloVelo to seven lineage-traced scRNA-seq datasets, generated using CRISPR-Cas9 editing, lentiviral barcoding or immune repertoire profiling, demonstrates its high accuracy and robustness in inferring complex lineage trajectories while outperforming RNA velocity. Additionally, we discovered that MEGs across tissues and organisms share similar functions in translation and ribosome biogenesis.

Modulating cell proliferation by asymmetric division: A conserved pattern in the early embryogenesis of nematode species

In the early stage of the nematode Caenorhabditis elegans embryogenesis, the zygote divides asymmetrically into a symmetric fast lineage and an asymmetric slow lineage, producing 16 and 8 cells respectively almost at the same time, followed by the onset of gastrulation. It was recently reported that this cell division pattern is optimal for rapid cell proliferation. In this work, we compare the cell lineages of 9 nematode species, revealing that this pattern is conserved for >60 million years. It further suggests that such lineage design has an important functional role and it might speed up embryonic development in the nematode kingdom, not limited to C. elegans , and independent of the maternal-zygotic transition dynamics.

Studying temporal dynamics of single cells: expression, lineage and regulatory networks

Learning how multicellular organs are developed from single cells to different cell types is a fundamental problem in biology. With the high-throughput scRNA-seq technology, computational methods have been developed to reveal the temporal dynamics of single cells from transcriptomic data, from phenomena on cell trajectories to the underlying mechanism that formed the trajectory. There are several distinct families of computational methods including Trajectory Inference (TI), Lineage Tracing (LT), and Gene Regulatory Network (GRN) Inference which are involved in such studies. This review summarizes these computational approaches which use scRNA-seq data to study cell differentiation and cell fate specification as well as the advantages and limitations of different methods. We further discuss how GRNs can potentially affect cell fate decisions and trajectory structures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-023-01090-5.

Charting Single Cell Lineage Dynamics and Mutation Networks via Homing CRISPR

Single cell lineage tracing, essential for unraveling cellular dynamics in disease evolution is critical for developing targeted therapies. CRISPR-Cas9, known for inducing permanent and cumulative mutations, is a cornerstone in lineage tracing. The novel homing guide RNA (hgRNA) technology enhances this by enabling dynamic retargeting and facilitating ongoing genetic modifications. Charting these mutations, especially through successive hgRNA edits, poses a significant challenge. Our solution, LINEMAP, is a computational framework designed to trace and map these mutations with precision. LINEMAP meticulously discerns mutation alleles at single-cell resolution and maps their complex interrelationships through a mutation evolution network. By utilizing a Markov Process model, we can predict mutation transition probabilities, revealing potential mutational routes and pathways. Our reconstruction algorithm, anchored in the Markov model's attributes, reconstructs cellular lineage pathways, shedding light on the cell's evolutionary journey to the minutiae of single-cell division. Our findings reveal an intricate network of mutation evolution paired with a predictive Markov model, advancing our capability to reconstruct single-cell lineage via hgRNA. This has substantial implications for advancing our understanding of biological mechanisms and propelling medical research forward.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

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

Lineage tracing technology using CRISPR/Cas9 genome editing has enabled simultaneous readouts of gene expressions and lineage barcodes in single cells, which allows for inference of cell lineage and cell types at the whole organism level. While most state-of-the-art methods for lineage reconstruction utilize only the lineage barcode data, methods that incorporate gene expressions are emerging. Effectively incorporating the gene expression data requires a reasonable model of how gene expression data changes along generations of divisions. Here, we present LinRace (Lineage Reconstruction with asymmetric cell division model), which integrates lineage barcode and gene expression data using asymmetric cell division model and infers cell lineages and ancestral cell states using Neighbor-Joining and maximum-likelihood heuristics. On both simulated and real data, LinRace outputs more accurate cell division trees than existing methods. With inferred ancestral states, LinRace can also show how a progenitor cell generates a large population of cells with various functionalities.

Machine learning based lineage tree reconstruction improved with knowledge of higher level relationships between cells and genomic barcodes

Tracking cells as they divide and progress through differentiation is a fundamental step in understanding many biological processes, such as the development of organisms and progression of diseases. In this study, we investigate a machine learning approach to reconstruct lineage trees in experimental systems based on mutating synthetic genomic barcodes. We refine previously proposed methodology by embedding information of higher level relationships between cells and single-cell barcode values into a feature space. We test performance of the algorithm on shallow trees (up to 100 cells) and deep trees (up to 10 000 cells). Our proposed algorithm can improve tree reconstruction accuracy in comparison to reconstructions based on a maximum parsimony method, but this comes at a higher computational time requirement.

LinRace: single cell lineage reconstruction using paired lineage barcode and gene expression data

Understanding how single cells divide and differentiate into different cell types in developed organs is one of the major tasks of developmental and stem cell biology. Recently, lineage tracing technology using CRISPR/Cas9 genome editing has enabled simultaneous readouts of gene expressions and lineage barcodes in single cells, which allows for the reconstruction of the cell division tree, and even the detection of cell types and differentiation trajectories at the whole organism level. While most state-of-the-art methods for lineage reconstruction utilize only the lineage barcode data, methods that incorporate gene expression data are emerging, aiming to improve the accuracy of lineage reconstruction. However, effectively incorporating the gene expression data requires a reasonable model on how gene expression data changes along generations of divisions. Here, we present LinRace (Lineage Reconstruction with asymmetric cell division model), a method that integrates the lineage barcode and gene expression data using the asymmetric cell division model and infers cell lineage under a framework combining Neighbor Joining and maximum-likelihood heuristics. On both simulated and real data, LinRace outputs more accurate cell division trees than existing methods. Moreover, LinRace can output the cell states (cell types) of ancestral cells, which is rarely performed with existing lineage reconstruction methods. The information on ancestral cells can be used to analyze how a progenitor cell generates a large population of cells with various functionalities. LinRace is available at: https://github.com/ZhangLabGT/LinRace.

The shaky foundations of simulating single-cell RNA sequencing data

BACKGROUND

With the emergence of hundreds of single-cell RNA-sequencing (scRNA-seq) datasets, the number of computational tools to analyze aspects of the generated data has grown rapidly. As a result, there is a recurring need to demonstrate whether newly developed methods are truly performant-on their own as well as in comparison to existing tools. Benchmark studies aim to consolidate the space of available methods for a given task and often use simulated data that provide a ground truth for evaluations, thus demanding a high quality standard results credible and transferable to real data.

RESULTS

Here, we evaluated methods for synthetic scRNA-seq data generation in their ability to mimic experimental data. Besides comparing gene- and cell-level quality control summaries in both one- and two-dimensional settings, we further quantified these at the batch- and cluster-level. Secondly, we investigate the effect of simulators on clustering and batch correction method comparisons, and, thirdly, which and to what extent quality control summaries can capture reference-simulation similarity.

CONCLUSIONS

Our results suggest that most simulators are unable to accommodate complex designs without introducing artificial effects, they yield over-optimistic performance of integration and potentially unreliable ranking of clustering methods, and it is generally unknown which summaries are important to ensure effective simulation-based method comparisons.

TEMPO enables sequential genetic labeling and manipulation of vertebrate cell lineages

During development, regulatory factors appear in a precise order to determine cell fates over time. Consequently, to investigate complex tissue development, it is necessary to visualize and manipulate cell lineages with temporal control. Current strategies for tracing vertebrate cell lineages lack genetic access to sequentially produced cells. Here, we present TEMPO (Temporal Encoding and Manipulation in a Predefined Order), an imaging-readable genetic tool allowing differential labeling and manipulation of consecutive cell generations in vertebrates. TEMPO is based on CRISPR and powered by a cascade of gRNAs that drive orderly activation and inactivation of reporters and/or effectors. Using TEMPO to visualize zebrafish and mouse neurogenesis, we recapitulated birth-order-dependent neuronal fates. Temporally manipulating cell-cycle regulators in mouse cortex progenitors altered the proportion and distribution of neurons and glia, revealing the effects of temporal gene perturbation on serial cell fates. Thus, TEMPO enables sequential manipulation of molecular factors, crucial to study cell-type specification.

Quantitative fate mapping: A general framework for analyzing progenitor state dynamics via retrospective lineage barcoding

Natural and induced somatic mutations that accumulate in the genome during development record the phylogenetic relationships of cells; whether these lineage barcodes capture the complex dynamics of progenitor states remains unclear. We introduce quantitative fate mapping, an approach to reconstruct the hierarchy, commitment times, population sizes, and commitment biases of intermediate progenitor states during development based on a time-scaled phylogeny of their descendants. To reconstruct time-scaled phylogenies from lineage barcodes, we introduce Phylotime, a scalable maximum likelihood clustering approach based on a general barcoding mutagenesis model. We validate these approaches using realistic in silico and in vitro barcoding experiments. We further establish criteria for the number of cells that must be analyzed for robust quantitative fate mapping and a progenitor state coverage statistic to assess the robustness. This work demonstrates how lineage barcodes, natural or synthetic, enable analyzing progenitor fate and dynamics long after embryonic development in any organism.

Connecting past and present: single-cell lineage tracing

Central to the core principle of cell theory, depicting cells' history, state and fate is a fundamental goal in modern biology. By leveraging clonal analysis and single-cell RNA-seq technologies, single-cell lineage tracing provides new opportunities to interrogate both cell states and lineage histories. During the past few years, many strategies to achieve lineage tracing at single-cell resolution have been developed, and three of them (integration barcodes, polylox barcodes, and CRISPR barcodes) are noteworthy as they are amenable in experimentally tractable systems. Although the above strategies have been demonstrated in animal development and stem cell research, much care and effort are still required to implement these methods. Here we review the development of single-cell lineage tracing, major characteristics of the cell barcoding strategies, applications, as well as technical considerations and limitations, providing a guide to choose or improve the single-cell barcoding lineage tracing.

Cellular barcoding to decipher clonal dynamics in disease

Cellular barcodes are distinct DNA sequences that enable one to track specific cells across time or space. Recent advances in our ability to detect natural or synthetic cellular barcodes, paired with single-cell readouts of cell state, have markedly increased our knowledge of clonal dynamics and genealogies of the cells that compose a variety of tissues and organs. These advances hold promise to redefine our view of human disease. Here, we provide an overview of cellular barcoding approaches, discuss applications to gain new insights into disease mechanisms, and provide an outlook on future applications. We discuss unanticipated insights gained through barcoding in studies of cancer and blood cell production and describe how barcoding can be applied to a growing array of medical fields, particularly with the increasing recognition of clonal contributions in human diseases.

State-of-the-art CRISPR for in vivo and cell-based studies in Drosophila

For more than 100 years, the fruit fly, Drosophila melanogaster, has served as a powerful model organism for biological and biomedical research due to its many genetic and physiological similarities to humans and the availability of sophisticated technologies used to manipulate its genome and genes. The Drosophila research community quickly adopted CRISPR technologies and, in the 8 years since the first clustered regularly interspaced short palindromic repeats (CRISPR) publications in flies, has explored and innovated methods for mutagenesis, precise genome engineering, and beyond. Moreover, the short lifespan and ease of genetics have made Drosophila an ideal testing ground for in vivo applications and refinements of the rapidly evolving set of CRISPR-associated (CRISPR-Cas) tools. Here, we review innovations in delivery of CRISPR reagents, increased efficiency of cutting and homology-directed repair (HDR), and alternatives to standard Cas9-based approaches. While the focus is primarily on in vivo systems, we also describe the role of Drosophila cultured cells as both an indispensable first step in the process of assessing new CRISPR technologies and a platform for genome-wide CRISPR pooled screens.

A New Generation of Lineage Tracing Dynamically Records Cell Fate Choices

Reconstructing the development of lineage relationships and cell fate mapping has been a fundamental problem in biology. Using advanced molecular biology and single-cell RNA sequencing, we have profiled transcriptomes at the single-cell level and mapped cell fates during development. Recently, CRISPR/Cas9 barcode editing for large-scale lineage tracing has been used to reconstruct the pseudotime trajectory of cells and improve lineage tracing accuracy. This review presents the progress of the latest CbLT (CRISPR-based Lineage Tracing) and discusses the current limitations and potential technical pitfalls in their application and other emerging concepts.

TedSim: temporal dynamics simulation of single-cell RNA sequencing data and cell division history

Recently, lineage tracing technology using CRISPR/Cas9 genome editing has enabled simultaneous readouts of gene expressions and lineage barcodes, which allows for the reconstruction of the cell division tree and makes it possible to reconstruct ancestral cell types and trace the origin of each cell type. Meanwhile, trajectory inference methods are widely used to infer cell trajectories and pseudotime in a dynamic process using gene expression data of present-day cells. Here, we present TedSim (single-cell temporal dynamics simulator), which simulates the cell division events from the root cell to present-day cells, simultaneously generating two data modalities for each single cell: the lineage barcode and gene expression data. TedSim is a framework that connects the two problems: lineage tracing and trajectory inference. Using TedSim, we conducted analysis to show that (i) TedSim generates realistic gene expression and barcode data, as well as realistic relationships between these two data modalities; (ii) trajectory inference methods can recover the underlying cell state transition mechanism with balanced cell type compositions; and (iii) integrating gene expression and barcode data can provide more insights into the temporal dynamics in cell differentiation compared to using only one type of data, but better integration methods need to be developed.

Deep distributed computing to reconstruct extremely large lineage trees

Phylogeny estimation (the reconstruction of evolutionary trees) has recently been applied to CRISPR-based cell lineage tracing, allowing the developmental history of an individual tissue or organism to be inferred from a large number of mutated sequences in somatic cells. However, current computational methods are not able to construct phylogenetic trees from extremely large numbers of input sequences. Here, we present a deep distributed computing framework to comprehensively trace accurate large lineages (FRACTAL) that substantially enhances the scalability of current lineage estimation software tools. FRACTAL first reconstructs only an upstream lineage of the input sequences and recursively iterates the same produce for its downstream lineages using independent computing nodes. We demonstrate the utility of FRACTAL by reconstructing lineages from >235 million simulated sequences and from >16 million cells from a simulated experiment with a CRISPR system that accumulates mutations during cell proliferation. We also successfully applied FRACTAL to evolutionary tree reconstructions and to an experiment using error-prone PCR (EP-PCR) for large-scale sequence diversification.

Orthogonal CRISPR-Cas tools for genome editing, inhibition, and CRISPR recording in zebrafish embryos

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas universe continues to expand. The type II CRISPR-Cas system from Streptococcus pyogenes (SpyCas9) is the most widely used for genome editing due to its high efficiency in cells and organisms. However, concentrating on a single CRISPR-Cas system imposes limits on target selection and multiplexed genome engineering. We hypothesized that CRISPR-Cas systems originating from different bacterial species could operate simultaneously and independently due to their distinct single-guide RNAs (sgRNAs) or CRISPR-RNAs (crRNAs), and protospacer adjacent motifs (PAMs). Additionally, we hypothesized that CRISPR-Cas activity in zebrafish could be regulated through the expression of inhibitory anti-CRISPR (Acr) proteins. Here, we use a simple mutagenesis approach to demonstrate that CRISPR-Cas systems from S. pyogenes (SpyCas9), Streptococcus aureus (SauCas9), Lachnospiraceae bacterium (LbaCas12a, previously known as LbCpf1) are orthogonal systems capable of operating simultaneously in zebrafish. CRISPR systems from Acidaminococcus sp. (AspCas12a, previously known as AsCpf1) and Neisseria meningitidis (Nme2Cas9) were also active in embryos. We implemented multichannel CRISPR recording using three CRISPR systems and show that LbaCas12a may provide superior information density compared with previous methods. We also demonstrate that type II Acrs (anti-CRISPRs) are effective inhibitors of SpyCas9 in zebrafish. Our results indicate that at least five CRISPR-Cas systems and two anti-CRISPR proteins are functional in zebrafish embryos. These orthogonal CRISPR-Cas systems and Acr proteins will enable combinatorial and intersectional strategies for spatiotemporal control of genome editing and genetic recording in animals.

Mapping single-cell-resolution cell phylogeny reveals cell population dynamics during organ development

Mapping the cell phylogeny of a complex multicellular organism relies on somatic mutations accumulated from zygote to adult. Available cell barcoding methods can record about three mutations per barcode, enabling only low-resolution mapping of the cell phylogeny of complex organisms. Here we developed SMALT, a substitution mutation-aided lineage-tracing system that outperforms the available cell barcoding methods in mapping cell phylogeny. We applied SMALT to Drosophila melanogaster and obtained on average more than 20 mutations on a three-kilobase-pair barcoding sequence in early-adult cells. Using the barcoding mutations, we obtained high-quality cell phylogenetic trees, each comprising several thousand internal nodes with 84-93% median bootstrap support. The obtained cell phylogenies enabled a population genetic analysis that estimates the longitudinal dynamics of the number of actively dividing parental cells (Np) in each organ through development. The Np dynamics revealed the trajectory of cell births and provided insight into the balance of symmetric and asymmetric cell division.

Spatial Organization of Chromatin: Emergence of Chromatin Structure During Development

Nuclei are central hubs for information processing in eukaryotic cells. The need to fit large genomes into small nuclei imposes severe restrictions on genome organization and the mechanisms that drive genome-wide regulatory processes. How a disordered polymer such as chromatin, which has vast heterogeneity in its DNA and histone modification profiles, folds into discernibly consistent patterns is a fundamental question in biology. Outstanding questions include how genomes are spatially and temporally organized to regulate cellular processes with high precision and whether genome organization is causally linked to transcription regulation. The advent of next-generation sequencing, super-resolution imaging, multiplexed fluorescent in situ hybridization, and single-molecule imaging in individual living cells has caused a resurgence in efforts to understand the spatiotemporal organization of the genome. In this review, we discuss structural and mechanistic properties of genome organization at different length scales and examine changes in higher-order chromatin organization during important developmental transitions.

Current approaches to fate mapping and lineage tracing using image data

Visualizing, tracking and reconstructing cell lineages in developing embryos has been an ongoing effort for well over a century. Recent advances in light microscopy, labelling strategies and computational methods to analyse complex image datasets have enabled detailed investigations into the fates of cells. Combined with powerful new advances in genomics and single-cell transcriptomics, the field of developmental biology is able to describe the formation of the embryo like never before. In this Review, we discuss some of the different strategies and applications to lineage tracing in live-imaging data and outline software methodologies that can be applied to various cell-tracking challenges.

Benchmarked approaches for reconstruction of in vitro cell lineages and in silico models of C. elegans and M. musculus developmental trees

The recent advent of CRISPR and other molecular tools enabled the reconstruction of cell lineages based on induced DNA mutations and promises to solve the ones of more complex organisms. To date, no lineage reconstruction algorithms have been rigorously examined for their performance and robustness across dataset types and number of cells. To benchmark such methods, we decided to organize a DREAM challenge using in vitro experimental intMEMOIR recordings and in silico data for a C. elegans lineage tree of about 1,000 cells and a Mus musculus tree of 10,000 cells. Some of the 22 approaches submitted had excellent performance, but structural features of the trees prevented optimal reconstructions. Using smaller sub-trees as training sets proved to be a good approach for tuning algorithms to reconstruct larger trees. The simulation and reconstruction methods here generated delineate a potential way forward for solving larger cell lineage trees such as in mouse.

LineageOT is a unified framework for lineage tracing and trajectory inference

Understanding the genetic and epigenetic programs that control differentiation during development is a fundamental challenge, with broad impacts across biology and medicine. Measurement technologies like single-cell RNA-sequencing and CRISPR-based lineage tracing have opened new windows on these processes, through computational trajectory inference and lineage reconstruction. While these two mathematical problems are deeply related, methods for trajectory inference are not typically designed to leverage information from lineage tracing and vice versa. Here, we present LineageOT, a unified framework for lineage tracing and trajectory inference. Specifically, we leverage mathematical tools from graphical models and optimal transport to reconstruct developmental trajectories from time courses with snapshots of both cell states and lineages. We find that lineage data helps disentangle complex state transitions with increased accuracy using fewer measured time points. Moreover, integrating lineage tracing with trajectory inference in this way could enable accurate reconstruction of developmental pathways that are impossible to recover with state-based methods alone.

CeLaVi: an interactive cell lineage visualization tool

Recent innovations in genetics and imaging are providing the means to reconstruct cell lineages, either by tracking cell divisions using live microscopy, or by deducing the history of cells using molecular recorders. A cell lineage on its own, however, is simply a description of cell divisions as branching events. A major goal of current research is to integrate this description of cell relationships with information about the spatial distribution and identities of the cells those divisions produce. Visualizing, interpreting and exploring these complex data in an intuitive manner requires the development of new tools. Here we present CeLaVi, a web-based visualization tool that allows users to navigate and interact with a representation of cell lineages, whilst simultaneously visualizing the spatial distribution, identities and properties of cells. CeLaVi’s principal functions include the ability to explore and manipulate the cell lineage tree; to visualise the spatial distribution of cell clones at different depths of the tree; to colour cells in the 3D viewer based on lineage relationships; to visualise various cell qualities on the 3D viewer (e.g. gene expression, cell type) and to annotate selected cells/clones. All these capabilities are demonstrated with four different example data sets. CeLaVi is available at http://www.celavi.pro.

Imaging cell lineage with a synthetic digital recording system

During multicellular development, spatial position and lineage history play powerful roles in controlling cell fate decisions. Using a serine integrase-based recording system, we engineered cells to record lineage information in a format that can be read out in situ. The system, termed integrase-editable memory by engineered mutagenesis with optical in situ readout (intMEMOIR), allowed in situ reconstruction of lineage relationships in cultured mouse cells and flies. intMEMOIR uses an array of independent three-state genetic memory elements that can recombine stochastically and irreversibly, allowing up to 59,049 distinct digital states. It reconstructed lineage trees in stem cells and enabled simultaneous analysis of single-cell clonal history, spatial position, and gene expression in Drosophila brain sections. These results establish a foundation for microscopy-readable lineage recording and analysis in diverse systems.

ESTIMATION OF CELL LINEAGE TREES BY MAXIMUM-LIKELIHOOD PHYLOGENETICS

CRISPR technology has enabled cell lineage tracing for complex multicellular organisms through insertion-deletion mutations of synthetic genomic barcodes during organismal development. To reconstruct the cell lineage tree from the mutated barcodes, current approaches apply general-purpose computational tools that are agnostic to the mutation process and are unable to take full advantage of the data's structure. We propose a statistical model for the CRISPR mutation process and develop a procedure to estimate the resulting tree topology, branch lengths, and mutation parameters by iteratively applying penalized maximum likelihood estimation. By assuming the barcode evolves according to a molecular clock, our method infers relative ordering across parallel lineages, whereas existing techniques only infer ordering for nodes along the same lineage. When analyzing transgenic zebrafish data from McKenna, Findlay and Gagnon et al. (2016), we find that our method recapitulates known aspects of zebrafish development and the results are consistent across samples.

Single-cell lineages reveal the rates, routes, and drivers of metastasis in cancer xenografts

Detailed phylogenies of tumor populations can recount the history and chronology of critical events during cancer progression, such as metastatic dissemination. We applied a Cas9-based, single-cell lineage tracer to study the rates, routes, and drivers of metastasis in a lung cancer xenograft mouse model. We report deeply resolved phylogenies for tens of thousands of cancer cells traced over months of growth and dissemination. This revealed stark heterogeneity in metastatic capacity, arising from preexisting and heritable differences in gene expression. We demonstrate that these identified genes can drive invasiveness and uncovered an unanticipated suppressive role for KRT17 We also show that metastases disseminated via multidirectional tissue routes and complex seeding topologies. Overall, we demonstrate the power of tracing cancer progression at subclonal resolution and vast scale.

Phylodynamics for cell biologists

Multicellular organisms are composed of cells connected by ancestry and descent from progenitor cells. The dynamics of cell birth, death, and inheritance within an organism give rise to the fundamental processes of development, differentiation, and cancer. Technical advances in molecular biology now allow us to study cellular composition, ancestry, and evolution at the resolution of individual cells within an organism or tissue. Here, we take a phylogenetic and phylodynamic approach to single-cell biology. We explain how "tree thinking" is important to the interpretation of the growing body of cell-level data and how ecological null models can benefit statistical hypothesis testing. Experimental progress in cell biology should be accompanied by theoretical developments if we are to exploit fully the dynamical information in single-cell data.

Next-Generation Lineage Tracing and Fate Mapping to Interrogate Development

Lineage tracing and fate mapping, overlapping yet distinct disciplines to follow cells and their progeny, have evolved rapidly over the last century. Lineage tracing aims to identify all progeny arising from an individual cell, placing them within a lineage hierarchy. The recent emergence of genomic technologies, such as single-cell and spatial transcriptomics, has fostered sophisticated new methods to reconstruct lineage relationships at high resolution. In contrast, fate maps, schematics showing which parts of the embryo will develop into which tissue, have remained relatively static since the 1970s. However, fate maps provide spatial information, often lost in lineage reconstruction, that can offer fundamental mechanistic insight into development. Here, we broadly review the origins of fate mapping and lineage tracing approaches. We focus on the most recent developments in lineage tracing, permitted by advances in single-cell genomics. Finally, we explore the current potential to leverage these new technologies to synthesize high-resolution fate maps and discuss their potential for interrogating development at new depths.

Lineage tracing meets single-cell omics: opportunities and challenges

A fundamental goal of developmental and stem cell biology is to map the developmental history (ontogeny) of differentiated cell types. Recent advances in high-throughput single-cell sequencing technologies have enabled the construction of comprehensive transcriptional atlases of adult tissues and of developing embryos from measurements of up to millions of individual cells. Parallel advances in sequencing-based lineage-tracing methods now facilitate the mapping of clonal relationships onto these landscapes and enable detailed comparisons between molecular and mitotic histories. Here we review recent progress and challenges, as well as the opportunities that emerge when these two complementary representations of cellular history are synthesized into integrated models of cell differentiation.

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

Recent studies combine two novel technologies, single-cell RNA-sequencing and CRISPR-Cas9 barcode editing for elucidating developmental lineages at the whole organism level. While these studies provided several insights, they face several computational challenges. First, lineages are reconstructed based on noisy and often saturated random mutation data. Additionally, due to the randomness of the mutations, lineages from multiple experiments cannot be combined to reconstruct a species-invariant lineage tree. To address these issues we developed a statistical method, LinTIMaT, which reconstructs cell lineages using a maximum-likelihood framework by integrating mutation and expression data. Our analysis shows that expression data helps resolve the ambiguities arising in when lineages are inferred based on mutations alone, while also enabling the integration of different individual lineages for the reconstruction of an invariant lineage tree. LinTIMaT lineages have better cell type coherence, improve the functional significance of gene sets and provide new insights on progenitors and differentiation pathways.

Base editors for simultaneous introduction of C-to-T and A-to-G mutations

We describe base editors that combine both cytosine and adenine base-editing functions. A codon-optimized fusion of the cytosine deaminase PmCDA1, the adenosine deaminase TadA and a Cas9 nickase (Target-ACEmax) showed a high median simultaneous C-to-T and A-to-G editing activity at 47 genomic targets. On-target as well as DNA and RNA off-target activities of Target-ACEmax were similar to those of existing single-function base editors.

Endogenous CRISPR/Cas9 arrays for scalable whole-organism lineage tracing

The past decade has seen a renewed appreciation of the central importance of cellular lineages to many questions in biology (especially organogenesis, stem cells and tumor biology). This has been driven in part by a renaissance in genetic clonal-labeling techniques. Recent approaches are based on accelerated mutation of DNA sequences, which can then be sequenced from individual cells to re-create a ‘phylogenetic’ tree of cell lineage. However, current approaches depend on making transgenic alterations to the genome in question, which limit their application. Here, we introduce a new method that completely avoids the need for prior genetic engineering, by identifying endogenous CRISPR/Cas9 target arrays suitable for lineage analysis. In both mouse and zebrafish, we identify the highest quality compact arrays as judged by equal base composition, 5′ G sequence, minimal likelihood of residing in the functional genome, minimal off targets and ease of amplification. We validate multiple high-quality endogenous CRISPR/Cas9 arrays, demonstrating their utility for lineage tracing. Our pragmatically scalable technique thus can produce deep and broad lineages in vivo, while removing the dependence on genetic engineering.

Inference of single-cell phylogenies from lineage tracing data using Cassiopeia

The pairing of CRISPR/Cas9-based gene editing with massively parallel single-cell readouts now enables large-scale lineage tracing. However, the rapid growth in complexity of data from these assays has outpaced our ability to accurately infer phylogenetic relationships. First, we introduce Cassiopeia-a suite of scalable maximum parsimony approaches for tree reconstruction. Second, we provide a simulation framework for evaluating algorithms and exploring lineage tracer design principles. Finally, we generate the most complex experimental lineage tracing dataset to date, 34,557 human cells continuously traced over 15 generations, and use it for benchmarking phylogenetic inference approaches. We show that Cassiopeia outperforms traditional methods by several metrics and under a wide variety of parameter regimes, and provide insight into the principles for the design of improved Cas9-enabled recorders. Together, these should broadly enable large-scale mammalian lineage tracing efforts. Cassiopeia and its benchmarking resources are publicly available at www.github.com/YosefLab/Cassiopeia.

Inclusion of temporal information in single cell transcriptomics

Single cell transcriptomics has emerged as a powerful method for dissecting cell type diversity and for understanding mechanisms of cell fate decisions. However, inclusion of temporal information remains challenging, since each cell can be measured only once by sequencing analysis. Here, we discuss recent progress and current efforts towards inclusion of temporal information in single cell transcriptomics. Even from snapshot data, temporal dynamics can be computationally inferred via pseudo-temporal ordering of single cell transcriptomes. Temporal information can also come from analysis of intronic reads or from RNA metabolic labeling, which can provide additional evidence for pseudo-time trajectories and enable more fine-grained analysis of gene regulatory interactions. These approaches measure dynamics on short timescales of hours. Emerging methods for high-throughput lineage tracing now enable information storage over long timescales by using CRISPR/Cas9 to record information in the genome, which can later be read out by sequencing.

Defining endogenous barcoding sites for CRISPR/Cas9-based cell lineage tracing in zebrafish

There is a growing interest in developing experimental methods for tracking the developmental cell lineages of a complex organism. The recently developed CRISPR/Cas9-based barcoding method is, although highly promising, difficult to scale up because it relies on exogenous barcoding sequences that are engineered into the genome. In this study, we characterized 78 high-quality endogenous sites in the zebrafish genome that can be used as CRISPR/Cas9-based barcoding sites. The 78 sites are all highly expressed in most of the cell types according to single-cell RNA sequencing (scRNA-seq) data. Hence, the barcoding information of the 78 endogenous sites is recovered by the available scRNA-seq platforms, enabling simultaneous characterization of cell type and cell lineage information.

Quantifying Hematopoietic Stem Cell Clonal Diversity by Selecting Informative Amplicon Barcodes

Hematopoietic stem cells (HSCs) are functionally and genetically diverse and this diversity decreases with age and disease. Numerous systems have been developed to quantify HSC diversity by genetic barcoding, but no framework has been established to empirically validate barcode sequences. Here we have developed an analytical framework, Selection of informative Amplicon Barcodes from Experimental Replicates (SABER), that identifies barcodes that are unique among a large set of experimental replicates. Amplicon barcodes were sequenced from the blood of 56 adult zebrafish divided into training and validation sets. Informative barcodes were identified and samples with a high fraction of informative barcodes were chosen by bootstrapping. There were 4.2 ± 1.8 barcoded HSC clones per sample in the training set and 3.5 ± 2.1 in the validation set (p = 0.3). SABER reproducibly quantifies functional HSCs and can accommodate a wide range of experimental group sizes. Future large-scale studies aiming to understand the mechanisms of HSC clonal evolution will benefit from this new approach to identifying informative amplicon barcodes.

High-throughput dense reconstruction of cell lineages

The first meeting exclusively dedicated to the 'High-throughput dense reconstruction of cell lineages' took place at Janelia Research Campus (Howard Hughes Medical Institute) from 14 to 18 April 2019. Organized by Tzumin Lee, Connie Cepko, Jorge Garcia-Marques and Isabel Espinosa-Medina, this meeting echoed the recent eruption of new tools that allow the reconstruction of lineages based on the phylogenetic analysis of DNA mutations induced during development. Combined with single-cell RNA sequencing, these tools promise to solve the lineage of complex model organisms at single-cell resolution. Here, we compile the conference consensus on the technological and computational challenges emerging from the use of the new strategies, as well as potential solutions.

Recording development with single cell dynamic lineage tracing

Every animal grows from a single fertilized egg into an intricate network of cell types and organ systems. This process is captured in a lineage tree: a diagram of every cell's ancestry back to the founding zygote. Biologists have long sought to trace this cell lineage tree in individual organisms and have developed a variety of technologies to map the progeny of specific cells. However, there are billions to trillions of cells in complex organisms, and conventional approaches can only map a limited number of clonal populations per experiment. A new generation of tools that use molecular recording methods integrated with single cell profiling technologies may provide a solution. Here, we summarize recent breakthroughs in these technologies, outline experimental and computational challenges, and discuss biological questions that can be addressed using single cell dynamic lineage tracing.

Is it possible to reconstruct an accurate cell lineage using CRISPR recorders?

Cell lineages provide the framework for understanding how cell fates are decided during development. Describing cell lineages in most organisms is challenging; even a fruit fly larva has ~50,000 cells and a small mammal has >1 billion cells. Recently, the idea of applying CRISPR to induce mutations during development, to be used as heritable markers for lineage reconstruction, has been proposed by several groups. While an attractive idea, its practical value depends on the accuracy of the cell lineages that can be generated. Here, we use computer simulations to estimate the performance of these approaches under different conditions. We incorporate empirical data on CRISPR-induced mutation frequencies in Drosophila. We show significant impacts from multiple biological and technical parameters - variable cell division rates, skewed mutational outcomes, target dropouts and different sequencing strategies. Our approach reveals the limitations of published CRISPR recorders, and indicates how future implementations can be optimised.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).