Recording development with single cell dynamic lineage tracing

Dual-nuclease single-cell lineage tracing by Cas9 and Cas12a

Single-cell lineage tracing based on CRISPR-Cas9 gene editing enables the simultaneous linkage of cell states and lineage history at a high resolution. Despite its immense potential in resolving the cell fate determination and genealogy within an organism, existing implementations of this technology suffer from limitations in recording capabilities and considerable barcode dropout. Here, we introduce DuTracer, a versatile tool that utilizes two orthogonal gene editing systems to record cell lineage history at single-cell resolution in an inducible manner. DuTracer shows the ability to enhance lineage recording with minimized target dropouts and potentially deeper tree depths. Applying DuTracer in mouse embryoid bodies and neuromesodermal organoids illustrates the lineage relationship of different cell types and proposes potential lineage-biased molecular drivers, showcased by identifying transcription factor Foxb1 as a modulator in the cell fate determination of neuromesodermal progenitors. Collectively, DuTracer facilitates the precise and regulatory interrogation of cellular lineages of complex biological processes.

Advancements in prospective single-cell lineage barcoding and their applications in research

Single-cell lineage tracing (scLT) has emerged as a powerful tool, providing unparalleled resolution to investigate cellular dynamics, fate determination, and the underlying molecular mechanisms. This review thoroughly examines the latest prospective lineage DNA barcode tracing technologies. It further highlights pivotal studies that leverage single-cell lentiviral integration barcoding technology to unravel the dynamic nature of cell lineages in both developmental biology and cancer research. Additionally, the review navigates through critical considerations for successful experimental design in lineage tracing and addresses challenges inherent in this field, including technical limitations, complexities in data analysis, and the imperative for standardization. It also outlines current gaps in knowledge and suggests future research directions, contributing to the ongoing advancement of scLT studies.

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 obtains mutations in ancestral cells where Notch is active. We combine SABER-seq with an expanded juvenile brain atlas to identify cell types derived from Notch-active founders. Our data reveal rare examples where differential Notch activities in ancestral progenitors are detected in terminally differentiated neuronal subtypes. SABER-seq is a novel platform for rapid, scalable and high-resolution mapping of signaling activity during development.

Modelling human brain development and disease with organoids

Organoids are systems derived from pluripotent stem cells at the interface between traditional monolayer cultures and in vivo animal models. The structural and functional characteristics of organoids enable the modelling of early stages of brain development in a physiologically relevant 3D environment. Moreover, organoids constitute a tool with which to analyse how individual genetic variation contributes to the susceptibility and progression of neurodevelopmental disorders. This Roadmap article describes the features of brain organoids, focusing on the neocortex, and their advantages and limitations - in comparison with other model systems - for the study of brain development, evolution and disease. We highlight avenues for enhancing the physiological relevance of brain organoids by integrating bioengineering techniques and unbiased high-throughput analyses, and discuss future applications. As organoids advance in mimicking human brain functions, we address the ethical and societal implications of this technology.

The lives of cells, recorded

A paradigm for biology is emerging in which cells can be genetically programmed to write their histories into their own genomes. These records can subsequently be read, and the cellular histories reconstructed, which for each cell could include a record of its lineage relationships, extrinsic influences, internal states and physical locations, over time. DNA recording has the potential to transform the way that we study developmental and disease processes. Recent advances in genome engineering are driving the development of systems for DNA recording, and meanwhile single-cell and spatial omics technologies increasingly enable the recovery of the recorded information. Combined with advances in computational and phylogenetic inference algorithms, the DNA recording paradigm is beginning to bear fruit. In this Perspective, we explore the rationale and technical basis of DNA recording, what aspects of cellular biology might be recorded and how, and the types of discovery that we anticipate this paradigm will enable.

A multimodal zebrafish developmental atlas reveals the state-transition dynamics of late-vertebrate pluripotent axial progenitors

Elucidating organismal developmental processes requires a comprehensive understanding of cellular lineages in the spatial, temporal, and molecular domains. In this study, we introduce Zebrahub, a dynamic atlas of zebrafish embryonic development that integrates single-cell sequencing time course data with lineage reconstructions facilitated by light-sheet microscopy. This atlas offers high-resolution and in-depth molecular insights into zebrafish development, achieved through the sequencing of individual embryos across ten developmental stages, complemented by reconstructions of cellular trajectories. Zebrahub also incorporates an interactive tool to navigate the complex cellular flows and lineages derived from light-sheet microscopy data, enabling in silico fate-mapping experiments. To demonstrate the versatility of our multimodal resource, we utilize Zebrahub to provide fresh insights into the pluripotency of neuro-mesodermal progenitors (NMPs) and the origins of a joint kidney-hemangioblast progenitor population.

Human Brain Ancestral Barcodes

Dynamic CpG methylation "barcodes" were read from 15,000 to 21,000 single cells from three human male brains. To overcome sparse sequencing coverage, the barcode had ~31,000 rapidly fluctuating X-chromosome CpG sites (fCpGs), with at least 500 covered sites per cell and at least 30 common sites between cell pairs (average of ~48). Barcodes appear to start methylated and record mitotic ages because excitatory neurons and glial cells that emerge later in development were less methylated. Barcodes are different between most cells, with average pairwise differences (PWDs) of ~0.5 between cells. About 10 cell pairs per million were more closely related with PWDs < 0.05. Barcodes appear to record ancestry and reconstruct trees where more related cells had similar phenotypes, albeit some pairs had phenotypic differences. Inhibitory neurons showed more evidence of tangential migration than excitatory neurons, with related cells in different cortical regions. fCpG barcodes become polymorphic during development and can distinguish between thousands of human cells.

The molecular chronology of mammary epithelial cell fate switching

The adult mammary gland is maintained by lineage-restricted progenitor cells through pregnancy, lactation, involution, and menopause. Injury resolution, transplantation-associated mammary gland reconstitution, and tumorigenesis are unique exceptions, wherein mammary basal cells gain the ability to reprogram to a luminal state. Here, we leverage newly developed cell-identity reporter mouse strains, and time-resolved single-cell epigenetic and transcriptomic analyses to decipher the molecular programs underlying basal-to-luminal fate switching in vivo. We demonstrate that basal cells rapidly reprogram toward plastic cycling intermediates that appear to hijack molecular programs we find in bipotent fetal mammary stem cells and puberty-associatiated cap cells. Loss of basal-cell specifiers early in dedifferentiation coincides with activation of Notch and BMP, among others. Pharmacologic blockade of each pathway disrupts basal-to-luminal transdifferentiation. Our studies provide a comprehensive map and resource for understanding the coordinated molecular changes enabling terminally differentiated epithelial cells to transition between cell lineages and highlights the stunning rapidity by which epigenetic reprogramming can occur in response to disruption of tissue structure.

Beyond the bulk: overview and novel insights into the dynamics of muscle satellite cells during muscle regeneration

Skeletal muscle possesses remarkable regenerative capabilities, fully recovering within a month following severe acute damage. Central to this process are muscle satellite cells (MuSCs), a resident population of somatic stem cells capable of self-renewal and differentiation. Despite the highly predictable course of muscle regeneration, evaluating this process has been challenging due to the heterogeneous nature of myogenic precursors and the limited insight provided by traditional markers with overlapping expression patterns. Notably, recent advancements in single-cell technologies, such as single-cell (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq), have revolutionized muscle research. These approaches allow for comprehensive profiling of individual cells, unveiling dynamic heterogeneity among myogenic precursors and their contributions to regeneration. Through single-cell transcriptome analyses, researchers gain valuable insights into cellular diversity and functional dynamics of MuSCs post-injury. This review aims to consolidate classical and new insights into the heterogeneity of myogenic precursors, including the latest discoveries from novel single-cell technologies.

Cellular proliferation biases clonal lineage tracing and trajectory inference

MOTIVATION

Lineage tracing and trajectory inference from single-cell RNA-sequencing data hold tremendous potential for uncovering the genetic programs driving development and disease. Single cell datasets are thought to provide an unbiased view on the diverse cellular architecture of tissues. Sampling bias, however, can skew single cell datasets away from the cellular composition they are meant to represent.

RESULTS

We demonstrate a novel form of sampling bias, caused by a statistical phenomenon related to repeated sampling from a growing, heterogeneous population. Relative growth rates of cells influence the probability that they will be sampled in clones observed across multiple time points. We support our probabilistic derivations with a simulation study and an analysis of a real time-course of T-cell development. We find that this bias can impact fate probability predictions, and we explore how to develop trajectory inference methods which are robust to this bias.

AVAILABILITY AND IMPLEMENTATION

Source code for the simulated datasets and to create the figures in this manuscript is freely available in python at https://github.com/rbonhamcarter/simulate-clones. A python implementation of the extension of the LineageOT method is freely available at https://github.com/rbonhamcarter/LineageOT/tree/multi-time-clones.

Multiscale engineering of brain organoids for disease modeling

Brain organoids hold great potential for modeling human brain development and pathogenesis. They recapitulate certain aspects of the transcriptional trajectory, cellular diversity, tissue architecture and functions of the developing brain. In this review, we explore the engineering strategies to control the molecular-, cellular- and tissue-level inputs to achieve high-fidelity brain organoids. We review the application of brain organoids in neural disorder modeling and emerging bioengineering methods to improve data collection and feature extraction at multiscale. The integration of multiscale engineering strategies and analytical methods has significant potential to advance insight into neurological disorders and accelerate drug development.

Clonal tracking in cancer and metastasis

The eradication of many cancers has proven challenging due to the presence of functionally and genetically heterogeneous clones maintained by rare cancer stem cells (CSCs), which contribute to disease progression, treatment refractoriness, and late relapse. The characterization of functional CSC activity has necessitated the development of modern clonal tracking strategies. This review describes viral-based and CRISPR-Cas9-based cellular barcoding, lineage tracing, and imaging-based approaches. DNA-based cellular barcoding technology is emerging as a powerful and robust strategy that has been widely applied to in vitro and in vivo model systems, including patient-derived xenograft models. This review also highlights the potential of these methods for use in the clinical and drug discovery contexts and discusses the important insights gained from such approaches.

Resolving tumor evolution: a phylogenetic approach

The evolutionary dynamics of cancer, characterized by its profound heterogeneity, demand sophisticated tools for a holistic understanding. This review delves into tumor phylogenetics, an essential approach bridging evolutionary biology with oncology, offering unparalleled insights into cancer's evolutionary trajectory. We provide an overview of the workflow, encompassing study design, data acquisition, and phylogeny reconstruction. Notably, the integration of diverse data sets emerges as a transformative step, enhancing the depth and breadth of evolutionary insights. With this integrated perspective, tumor phylogenetics stands poised to redefine our understanding of cancer evolution and influence therapeutic strategies.

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.

Engineering Plant Cell Fates and Functions for Agriculture and Industry

Many plant species are grown to enable access to specific organs or tissues, such as seeds, fruits, or stems. In some cases, a value is associated with a molecule that accumulates in a single type of cell. Domestication and subsequent breeding have often increased the yields of these target products by increasing the size, number, and quality of harvested organs and tissues but also via changes to overall plant growth architecture to suit large-scale cultivation. Many of the mutations that underlie these changes have been identified in key regulators of cellular identity and function. As key determinants of yield, these regulators are key targets for synthetic biology approaches to engineer new forms and functions. However, our understanding of many plant developmental programs and cell-type specific functions is still incomplete. In this Perspective, we discuss how advances in cellular genomics together with synthetic biology tools such as biosensors and DNA-recording devices are advancing our understanding of cell-specific programs and cell fates. We then discuss advances and emerging opportunities for cell-type-specific engineering to optimize plant morphology, responses to the environment, and the production of valuable compounds.

From Hematopoietic Stem Cells to Platelets: Unifying Differentiation Pathways Identified by Lineage Tracing Mouse Models

Platelets are the terminal progeny of megakaryocytes, primarily produced in the bone marrow, and play critical roles in blood homeostasis, clotting, and wound healing. Traditionally, megakaryocytes and platelets are thought to arise from multipotent hematopoietic stem cells (HSCs) via multiple discrete progenitor populations with successive, lineage-restricting differentiation steps. However, this view has recently been challenged by studies suggesting that (1) some HSC clones are biased and/or restricted to the platelet lineage, (2) not all platelet generation follows the "canonical" megakaryocytic differentiation path of hematopoiesis, and (3) platelet output is the default program of steady-state hematopoiesis. Here, we specifically investigate the evidence that in vivo lineage tracing studies provide for the route(s) of platelet generation and investigate the involvement of various intermediate progenitor cell populations. We further identify the challenges that need to be overcome that are required to determine the presence, role, and kinetics of these possible alternate pathways.

Punctuational evolution is pervasive in distal site metastatic colonization

The evolution of metastasis represents a lethal stage of cancer progression. Yet, the evolutionary kinetics of metastatic disease remain unresolved. Here, using single cell CRISPR-Cas9 lineage tracing data, we show that in metastatic disease, gradual molecular evolution is punctuated by episodes of rapid evolutionary change associated with lineage divergence. By measuring punctuational effects across the metastatic cascade, we show that punctuational effects contribute more to the molecular diversity at distal site metastases compared to the paired primary tumor, suggesting qualitatively different modes of evolution may drive primary and metastatic tumor progression. This is the first empirical evidence for distinct patterns of molecular evolution at early and late stages of metastasis and demonstrates the complex interplay of cell intrinsic and extrinsic factors that shape lethal cancer.

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

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 non-modifiable 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 time-scaled 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.

Bringing CAM photosynthesis to the table: Paving the way for resilient and productive agricultural systems in a changing climate

Modern agricultural systems are directly threatened by global climate change and the resulting freshwater crisis. A considerable challenge in the coming years will be to develop crops that can cope with the consequences of declining freshwater resources and changing temperatures. One approach to meeting this challenge may lie in our understanding of plant photosynthetic adaptations and water use efficiency. Plants from various taxa have evolved crassulacean acid metabolism (CAM), a water-conserving adaptation of photosynthetic carbon dioxide fixation that enables plants to thrive under semi-arid or seasonally drought-prone conditions. Although past research on CAM has led to a better understanding of the inner workings of plant resilience and adaptation to stress, successful introduction of this pathway into C3 or C4 plants has not been reported. The recent revolution in molecular, systems, and synthetic biology, as well as innovations in high-throughput data generation and mining, creates new opportunities to uncover the minimum genetic tool kit required to introduce CAM traits into drought-sensitive crops. Here, we propose four complementary research avenues to uncover this tool kit. First, genomes and computational methods should be used to improve understanding of the nature of variations that drive CAM evolution. Second, single-cell 'omics technologies offer the possibility for in-depth characterization of the mechanisms that trigger environmentally controlled CAM induction. Third, the rapid increase in new 'omics data enables a comprehensive, multimodal exploration of CAM. Finally, the expansion of functional genomics methods is paving the way for integration of CAM into farming systems.

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.

Harnessing CRISPR-Cas adaptation for RNA recording and beyond

Prokaryotes encode clustered regularly interspaced short palindromic repeat (CRISPR) arrays and CRISPR-associated (Cas) genes as an adaptive immune machinery. CRISPR-Cas systems effectively protect hosts from the invasion of foreign enemies, such as bacteriophages and plasmids. During a process called 'adaptation', non-self-nucleic acid fragments are acquired as spacers between repeats in the host CRISPR array, to establish immunological memory. The highly conserved Cas1-Cas2 complexes function as molecular recorders to integrate spacers in a time course manner, which can subsequently be expressed as crRNAs complexed with Cas effector proteins for the RNAguided interference pathways. In some of the RNA-targeting type III systems, Cas1 proteins are fused with reverse transcriptase (RT), indicating that RT-Cas1-Cas2 complexes can acquire RNA transcripts for spacer acquisition. In this review, we summarize current studies that focus on the molecular structure and function of the RT-fused Cas1-Cas2 integrase, and its potential applications as a directional RNA-recording tool in cells. Furthermore, we highlight outstanding questions for RT-Cas1-Cas2 studies and future directions for RNA-recording CRISPR technologies. [BMB Reports 2024; 57(1): 40-49].

DNA Barcoding Technology for Lineage Recording and Tracing to Resolve Cell Fate Determination

In various biological contexts, cells receive signals and stimuli that prompt them to change their current state, leading to transitions into a future state. This change underlies the processes of development, tissue maintenance, immune response, and the pathogenesis of various diseases. Following the path of cells from their initial identity to their current state reveals how cells adapt to their surroundings and undergo transformations to attain adjusted cellular states. DNA-based molecular barcoding technology enables the documentation of a phylogenetic tree and the deterministic events of cell lineages, providing the mechanisms and timing of cell lineage commitment that can either promote homeostasis or lead to cellular dysregulation. This review comprehensively presents recently emerging molecular recording technologies that utilize CRISPR/Cas systems, base editing, recombination, and innate variable sequences in the genome. Detailing their underlying principles, applications, and constraints paves the way for the lineage tracing of every cell within complex biological systems, encompassing the hidden steps and intermediate states of organism development and disease progression.

Startle: A star homoplasy approach for CRISPR-Cas9 lineage tracing

CRISPR-Cas9-based genome editing combined with single-cell sequencing enables the tracing of the history of cell divisions, or cellular lineage, in tissues and whole organisms. Although standard phylogenetic approaches may be applied to reconstruct cellular lineage trees from this data, the unique features of the CRISPR-Cas9 editing process motivate the development of specialized models that describe the evolution of CRISPR-Cas9-induced mutations. Here, we introduce the "star homoplasy" evolutionary model that constrains a phylogenetic character to mutate at most once along a lineage, capturing the "non-modifiability" property of CRISPR-Cas9 mutations. We derive a combinatorial characterization of star homoplasy phylogenies and use this characterization to develop an algorithm, "Startle", that computes a maximum parsimony star homoplasy phylogeny. We demonstrate that Startle infers more accurate phylogenies on simulated lineage tracing data compared with existing methods and finds parsimonious phylogenies with fewer metastatic migrations on lineage tracing data from mouse metastatic lung adenocarcinoma.

Eyes open on stem cells

Recently, the murine cornea has reemerged as a robust stem cell (SC) model, allowing individual SC tracing in living animals. The cornea has pioneered seminal discoveries in SC biology and regenerative medicine, from the first corneal transplantation in 1905 to the identification of limbal SCs and their transplantation to successfully restore vision in the early 1990s. Recent experiments have exposed unexpected properties attributed to SCs and progenitors and revealed flexibility in the differentiation program and a key role for the SC niche. Here, we discuss the limbal SC model and its broader relevance to other tissues, disease, and therapy.

Minimizing cell number fluctuations in self-renewing tissues with a stem-cell niche

Self-renewing tissues require that a constant number of proliferating cells is maintained over time. This maintenance can be ensured at the single-cell level or the population level. Maintenance at the population level leads to fluctuations in the number of proliferating cells over time. Often, it is assumed that those fluctuations can be reduced by increasing the number of asymmetric divisions, i.e., divisions where only one of the daughter cells remains proliferative. Here, we study a model of cell proliferation that incorporates a stem-cell niche of fixed size, and explicitly model the cells inside and outside the niche. We find that in this model, fluctuations are minimized when the difference in growth rate between the niche and the rest of the tissue is maximized and all divisions are symmetric divisions, producing either two proliferating or two nonproliferating daughters. We show that this optimal state leaves visible signatures in clone size distributions and could thus be detected experimentally.

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.

Best practices for single-cell analysis across modalities

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

Microenvironmental control of hematopoietic stem cell fate via CXCL8 and protein kinase C

Altered hematopoietic stem cell (HSC) fate underlies primary blood disorders but microenvironmental factors controlling this are poorly understood. Genetically barcoded genome editing of synthetic target arrays for lineage tracing (GESTALT) zebrafish were used to screen for factors expressed by the sinusoidal vascular niche that alter the phylogenetic distribution of the HSC pool under native conditions. Dysregulated expression of protein kinase C delta (PKC-δ, encoded by prkcda) increases the number of HSC clones by up to 80% and expands polyclonal populations of immature neutrophil and erythroid precursors. PKC agonists such as cxcl8 augment HSC competition for residency within the niche and expand defined niche populations. CXCL8 induces association of PKC-δ with the focal adhesion complex, activating extracellular signal-regulated kinase (ERK) signaling and expression of niche factors in human endothelial cells. Our findings demonstrate the existence of reserve capacity within the niche that is controlled by CXCL8 and PKC and has significant impact on HSC phylogenetic and phenotypic fate.

A reference cell tree will serve science better than a reference cell atlas

Single-cell biology is facing a crisis of sorts. Vast numbers of single-cell molecular profiles are being generated, clustered and annotated. However, this is overwhelmingly ad hoc, and we continue to lack a principled, unified, and well-moored system for defining, naming, and organizing cell types. In this perspective, we argue against an atlas or periodic table-like discretization as the right metaphor for a reference taxonomy of cell types. In its place, we advocate for a data-driven, tree-based nomenclature that is rooted in a "consensus ontogeny" spanning the life cycle of a given species. We explore how such a reference cell tree, inclusive of both lineage histories and molecular states, could be constructed, represented, and segmented in practice. Analogous to the taxonomic classification of species, a consensus ontogeny would provide a universal, stable, and extendable framework for precise scientific communication, both contemporaneously and across the ages.

The heart field transcriptional landscape at single-cell resolution

Organogenesis requires the orchestrated development of multiple cell lineages that converge, interact, and specialize to generate coherent functional structures, exemplified by transformation of the cardiac crescent into a four-chambered heart. Cardiomyocytes originate from the first and second heart fields, which make different regional contributions to the definitive heart. In this review, a series of recent single-cell transcriptomic analyses, together with genetic tracing experiments, are discussed, providing a detailed panorama of the cardiac progenitor cell landscape. These studies reveal that first heart field cells originate in a juxtacardiac field adjacent to extraembryonic mesoderm and contribute to the ventrolateral side of the cardiac primordium. In contrast, second heart field cells are deployed dorsomedially from a multilineage-primed progenitor population via arterial and venous pole pathways. Refining our knowledge of the origin and developmental trajectories of cells that build the heart is essential to address outstanding challenges in cardiac biology and disease.

Molecular recording: transcriptional data collection into the genome

Advances in regenerative medicine depend upon understanding the complex transcriptional choreography that guides cellular development. Transcriptional molecular recorders, tools that record different transcriptional events into the genome of cells, hold promise to elucidate both the intensity and timing of transcriptional activity at single-cell resolution without requiring destructive multitime point assays. These technologies are dependent on DNA writers, which translate transcriptional signals into stable genomic mutations that encode the duration, intensity, and order of transcriptional events. In this review, we highlight recent progress toward more informative and multiplexable transcriptional recording through the use of three different types of DNA writing - recombineering, Cas1-Cas2 acquisition, and prime editing - and the architecture of the genomic data generated.

Accurate estimation of molecular counts from amplicon sequence data with unique molecular identifiers

MOTIVATION

Amplicon sequencing is widely applied to explore heterogeneity and rare variants in genetic populations. Resolving true biological variants and quantifying their abundance is crucial for downstream analyses, but measured abundances are distorted by stochasticity and bias in amplification, plus errors during polymerase chain reaction (PCR) and sequencing. One solution attaches unique molecular identifiers (UMIs) to sample sequences before amplification. Counting UMIs instead of sequences provides unbiased estimates of abundance. While modern methods improve over naïve counting by UMI identity, most do not account for UMI reuse or collision, and they do not adequately model PCR and sequencing errors in the UMIs and sample sequences.

RESULTS

We introduce Deduplication and Abundance estimation with UMIs (DAUMI), a probabilistic framework to detect true biological amplicon sequences and accurately estimate their deduplicated abundance. DAUMI recognizes UMI collision, even on highly similar sequences, and detects and corrects most PCR and sequencing errors in the UMI and sampled sequences. DAUMI performs better on simulated and real data compared to other UMI-aware clustering methods.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/DormanLab/AmpliCI.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Size-Based Sorting and In Situ Clonal Expansion of Single Cells Using Microfluidics

Separation and clonal culture and growth kinetics analysis of target cells in a mixed population is critical for pathological research, disease diagnosis, and cell therapy. However, long-term culture with time-lapse imaging of the isolated cells for clonal analysis is still challenging. This paper reports a microfluidic device with four-level filtration channels and a pneumatic microvalve for size sorting and in situ clonal culture of single cells. The valve was on top of the filtration channels and used to direct fluid flow by membrane deformation during separation and long-term culture to avoid shear-induced cell deformation. Numerical simulations were performed to evaluate the influence of device parameters affecting the pressure drop across the filtration channels. Then, a droplet model was employed to evaluate the impact of cell viscosity, cell size, and channel width on the pressure drop inducing cell deformation. Experiments showed that filtration channels with a width of 7, 10, 13, or 17 μm successfully sorted K562 cells into four different size ranges at low driving pressure. The maximum efficiency of separating K562 cells from media and whole blood was 98.6% and 89.7%, respectively. Finally, the trapped single cells were cultured in situ for 4-7 days with time-lapse imaging to obtain the lineage trees and growth curves. Then, the time to the first division, variation of cell size before and after division, and cell fusion were investigated. This proved that cells at the G1 and G2 phases were of significantly distinct sizes. The microfluidic device for size sorting and clonal expansion will be of tremendous application potential in single-cell studies.

TiDeTree: a Bayesian phylogenetic framework to estimate single-cell trees and population dynamic parameters from genetic lineage tracing data

The development of organisms and tissues is dictated by an elaborate balance between cell division, apoptosis and differentiation: the cell population dynamics. To quantify these dynamics, we propose a phylodynamic inference approach based on single-cell lineage recorder data. We developed a Bayesian phylogenetic framework-time-scaled developmental trees (TiDeTree)-that uses lineage recorder data to estimate time-scaled single-cell trees. By implementing TiDeTree within BEAST 2, we enable joint inference of the time-scaled trees and the cell population dynamics. We validated TiDeTree using simulations and showed that performance further improves when including multiple independent sources of information into the inference, such as frequencies of editing outcomes or experimental replicates. We benchmarked TiDeTree against state-of-the-art methods and show comparable performance in terms of tree topology, plus direct assessment of uncertainty and co-estimation of additional parameters. To demonstrate TiDeTree's use in practice, we analysed a public dataset containing lineage data from approximately 100 stem cell colonies. We estimated a time-scaled phylogeny for each colony; as well as the cell division and apoptosis rates underlying the growth dynamics of all colonies. We envision that TiDeTree will find broad application in the analysis of single-cell lineage tracing data, which will improve our understanding of cellular processes during development.

Mastering the use of cellular barcoding to explore cancer heterogeneity

Tumours are often composed of a multitude of malignant clones that are genomically unique, and only a few of them may have the ability to escape cancer therapy and grow as symptomatic lesions. As a result, tumours with a large degree of genomic diversity have a higher chance of leading to patient death. However, clonal fate can be driven by non-genomic features. In this context, new technologies are emerging not only to track the spatiotemporal fate of individual cells and their progeny but also to study their molecular features using various omics analysis. In particular, the recent development of cellular barcoding facilitates the labelling of tens to millions of cancer clones and enables the identification of the complex mechanisms associated with clonal fate in different microenvironments and in response to therapy. In this Review, we highlight the recent discoveries made using lentiviral-based cellular barcoding techniques, namely genetic and optical barcoding. We also emphasize the strengths and limitations of each of these technologies and discuss some of the key concepts that must be taken into consideration when one is designing barcoding experiments. Finally, we suggest new directions to further improve the use of these technologies in cancer research.

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.

Spotlighting adult stem cells: advances, pitfalls, and challenges

The existence of stem cells (SCs) at the tip of the cellular differentiation hierarchy has fascinated the scientific community ever since their discovery in the early 1950s to 1960s. Despite the remarkable success of the SC theory and the development of SC-based treatments, fundamental features of SCs remain enigmatic. Recent advances in single-cell lineage tracing, live imaging, and genomic technologies have allowed capture of life histories and transcriptional signatures of individual cells, leaving SCs much less space to 'hide'. Focusing on epithelial SCs and comparing them to other SCs, we discuss new paradigms of the SC niche, dynamics, and pathology, highlighting key open questions in SC biology that need to be resolved for harnessing SC potential in regenerative medicine.

Current understanding of osteoarthritis pathogenesis and relevant new approaches

Osteoarthritis (OA) is the most common degenerative joint disease that causes painful swelling and permanent damage to the joints in the body. The molecular mechanisms of OA are currently unknown. OA is a heterogeneous disease that affects the entire joint, and multiple tissues are altered during OA development. To better understand the pathological mechanisms of OA, new approaches, methods, and techniques need to be used to understand OA pathogenesis. In this review, we first focus on the epigenetic regulation of OA, with a particular focus on DNA methylation, histone modification, and microRNA regulation, followed by a summary of several key mediators in OA-associated pain. We then introduce several innovative techniques that have been and will continue to be used in the fields of OA and OA-associated pain, such as CRISPR, scRNA sequencing, and lineage tracing. Next, we discuss the timely updates concerning cell death regulation in OA pathology, including pyroptosis, ferroptosis, and autophagy, as well as their individual roles in OA and potential molecular targets in treating OA. Finally, our review highlights new directions on the role of the synovial lymphatic system in OA. An improved understanding of OA pathogenesis will aid in the development of more specific and effective therapeutic interventions for OA.

A time-resolved, multi-symbol molecular recorder via sequential genome editing

DNA is naturally well suited to serve as a digital medium for in vivo molecular recording. However, contemporary DNA-based memory devices are constrained in terms of the number of distinct 'symbols' that can be concurrently recorded and/or by a failure to capture the order in which events occur1. Here we describe DNA Typewriter, a general system for in vivo molecular recording that overcomes these and other limitations. For DNA Typewriter, the blank recording medium ('DNA Tape') consists of a tandem array of partial CRISPR-Cas9 target sites, with all but the first site truncated at their 5' ends and therefore inactive. Short insertional edits serve as symbols that record the identity of the prime editing guide RNA2 mediating the edit while also shifting the position of the 'type guide' by one unit along the DNA Tape, that is, sequential genome editing. In this proof of concept of DNA Typewriter, we demonstrate recording and decoding of thousands of symbols, complex event histories and short text messages; evaluate the performance of dozens of orthogonal tapes; and construct 'long tape' potentially capable of recording as many as 20 serial events. Finally, we leverage DNA Typewriter in conjunction with single-cell RNA-seq to reconstruct a monophyletic lineage of 3,257 cells and find that the Poisson-like accumulation of sequential edits to multicopy DNA tape can be maintained across at least 20 generations and 25 days of in vitro clonal expansion.

What is a cell type and how to define it?

Cell types are the basic functional units of an organism. Cell types exhibit diverse phenotypic properties at multiple levels, making them challenging to define, categorize, and understand. This review provides an overview of the basic principles of cell types rooted in evolution and development and discusses approaches to characterize and classify cell types and investigate how they contribute to the organism's function, using the mammalian brain as a primary example. I propose a roadmap toward a conceptual framework and knowledge base of cell types that will enable a deeper understanding of the dynamic changes of cellular function under healthy and diseased conditions.

Specification of the haematopoietic stem cell lineage: From blood-fated mesodermal angioblasts to haemogenic endothelium

Haematopoietic stem and progenitor cells emerge from specialized haemogenic endothelial cells in select vascular beds during embryonic development. Specification and commitment to the blood lineage, however, occur before endothelial cells are endowed with haemogenic competence, at the time of mesoderm patterning and production of endothelial cell progenitors (angioblasts). Whilst early blood cell fate specification has long been recognized, very little is known about the mechanisms that induce endothelial cell diversification and progressive acquisition of a blood identity by a subset of these cells. Here, we review the endothelial origin of the haematopoietic system and the complex developmental journey of blood-fated angioblasts. We discuss how recent technological advances will be instrumental to examine the diversity of the embryonic anatomical niches, signaling pathways and downstream epigenetic and transcriptional processes controlling endothelial cell heterogeneity and blood cell fate specification. Ultimately, this will give essential insights into the ontogeny of the cells giving rise to haematopoietic stem cells, that may aid in the development of novel strategies for their in vitro production for clinical purposes.

Lineage tracing in human tissues

The dynamical process of cell division that underpins homeostasis in the human body cannot be directly observed in vivo, but instead is measurable from the pattern of somatic genetic or epigenetic mutations that accrue in tissues over an individual's lifetime. Because somatic mutations are heritable, they serve as natural lineage tracing markers that delineate clonal expansions. Mathematical analysis of the distribution of somatic clone sizes gives a quantitative readout of the rates of cell birth, death, and replacement. In this review we explore the broad range of somatic mutation types that have been used for lineage tracing in human tissues, introduce the mathematical concepts used to infer dynamical information from these clone size data, and discuss the insights of this lineage tracing approach for our understanding of homeostasis and cancer development. We use the human colon as a particularly instructive exemplar tissue. There is a rich history of human somatic cell dynamics surreptitiously written into the cell genomes that is being uncovered by advances in sequencing and careful mathematical analysis lineage of tracing data. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

As in Real Estate, Location Matters: Cellular Expression of Complement Varies Between Macular and Peripheral Regions of the Retina and Supporting Tissues

The cellular events that dictate the initiation of the complement pathway in ocular degeneration, such as age-related macular degeneration (AMD), is poorly understood. Using gene expression analysis (single cell and bulk), mass spectrometry, and immunohistochemistry, we dissected the role of multiple retinal and choroidal cell types in determining the complement homeostasis. Our scRNA-seq data show that the cellular response to early AMD is more robust in the choroid, particularly in fibroblasts, pericytes and endothelial cells. In late AMD, complement changes were more prominent in the retina especially with the expression of the classical pathway initiators. Notably, we found a spatial preference for these differences. Overall, this study provides insights into the heterogeneity of cellular responses for complement expression and the cooperation of neighboring cells to complete the pathway in healthy and AMD eyes. Further, our findings provide new cellular targets for therapies directed at complement.

Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution. We found that the loss of the initial, stable alveolar-type2-like state was accompanied by a transient increase in plasticity. This was followed by the adoption of distinct transcriptional programs that enable rapid expansion and, ultimately, clonal sweep of stable subclones capable of metastasizing. Finally, tumors develop through stereotypical evolutionary trajectories, and perturbing additional tumor suppressors accelerates progression by creating novel trajectories. Our study elucidates the hierarchical nature of tumor evolution and, more broadly, enables in-depth studies of tumor progression.

Interactive, integrated analysis of single-cell transcriptomic and phylogenetic data with PhyloVision

Recent advances in CRISPR-Cas9 engineering and single-cell assays have enabled the simultaneous measurement of single-cell transcriptomic and phylogenetic profiles. However, there are few computational tools enabling users to integrate and derive insight from a joint analysis of these two modalities. Here, we describe "PhyloVision": an open-source software for interactively exploring data from both modalities and for identifying and interpreting heritable gene modules whose concerted expression are associated with phylogenetic relationships. PhyloVision provides a feature-rich, interactive, and shareable web-based report for investigating these modules while also supporting several other data and meta-data exploration capabilities. We demonstrate the utility of PhyloVision using a published dataset of metastatic lung adenocarcinoma cells, whose phylogeny was resolved using a CRISPR-Cas9-based lineage-tracing system. Together, we anticipate that PhyloVision and the methods it implements will be a useful resource for scalable and intuitive data exploration for any assay that simultaneously measures cell state and lineage.

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.

Engineering Brain Organoids: Toward Mature Neural Circuitry with an Intact Cytoarchitecture

The emergence of brain organoids as a model system has been a tremendously exciting development in the field of neuroscience. Brain organoids are a gateway to exploring the intricacies of human-specific neurogenesis that have so far eluded the neuroscience community. Regardless, current culture methods have a long way to go in terms of accuracy and reproducibility. To perfectly mimic the human brain, we need to recapitulate the complex in vivo context of the human fetal brain and achieve mature neural circuitry with an intact cytoarchitecture. In this review, we explore the major challenges facing the current brain organoid systems, potential technical breakthroughs to advance brain organoid techniques up to levels similar to an in vivo human developing brain, and the future prospects of this technology.

Re-building limbs, one cell at a time

New techniques for visualizing and interrogating single cells hold the key to unlocking the underlying mechanisms of salamander limb regeneration.

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.