Whole-organism clone tracing using single-cell sequencing

A mouse model with high clonal barcode diversity for joint lineage, transcriptomic, and epigenomic profiling in single cells

Cellular lineage histories and their molecular states encode fundamental principles of tissue development and homeostasis. Current lineage-recording mouse models have insufficient barcode diversity and single-cell lineage coverage for profiling tissues composed of millions of cells. Here, we developed DARLIN, an inducible Cas9 barcoding mouse line that utilizes terminal deoxynucleotidyl transferase (TdT) and 30 CRISPR target sites. DARLIN is inducible, generates massive lineage barcodes across tissues, and enables the detection of edited barcodes in ∼70% of profiled single cells. Using DARLIN, we examined fate bias within developing hematopoietic stem cells (HSCs) and revealed unique features of HSC migration. Additionally, we established a protocol for joint transcriptomic and epigenomic single-cell measurements with DARLIN and found that cellular clonal memory is associated with genome-wide DNA methylation rather than gene expression or chromatin accessibility. DARLIN will enable the high-resolution study of lineage relationships and their molecular signatures in diverse tissues and physiological contexts.

Retrospective analysis of enhancer activity and transcriptome history

Cell state changes in development and disease are controlled by gene regulatory networks, the dynamics of which are difficult to track in real time. In this study, we used an inducible DCM-RNA polymerase subunit b fusion protein which labels active genes and enhancers with a bacterial methylation mark that does not affect gene transcription and is propagated in S-phase. This DCM-RNA polymerase fusion protein enables transcribed genes and active enhancers to be tagged and then examined at later stages of development or differentiation. We apply this DCM-time machine (DCM-TM) technology to study intestinal homeostasis, revealing rapid and coordinated activation of enhancers and nearby genes during enterocyte differentiation. We provide new insights in absorptive-secretory lineage decision-making in intestinal stem cell (ISC) differentiation and show that ISCs retain a unique chromatin landscape required to maintain ISC identity and delineate future expression of differentiation-associated genes. DCM-TM has wide applicability in tracking cell states, providing new insights in the regulatory networks underlying cell state changes.

Comprehensive spatiotemporal mapping of single-cell lineages in developing mouse brain by CRISPR-based barcoding

A fundamental interest in developmental neuroscience lies in the ability to map the complete single-cell lineages within the brain. To this end, we developed a CRISPR editing-based lineage-specific tracing (CREST) method for clonal tracing in Cre mice. We then used two complementary strategies based on CREST to map single-cell lineages in developing mouse ventral midbrain (vMB). By applying snapshotting CREST (snapCREST), we constructed a spatiotemporal lineage landscape of developing vMB and identified six progenitor archetypes that could represent the principal clonal fates of individual vMB progenitors and three distinct clonal lineages in the floor plate that specified glutamatergic, dopaminergic or both neurons. We further created pandaCREST (progenitor and derivative associating CREST) to associate the transcriptomes of progenitor cells in vivo with their differentiation potentials. We identified multiple origins of dopaminergic neurons and demonstrated that a transcriptome-defined progenitor type comprises heterogeneous progenitors, each with distinct clonal fates and molecular signatures. Therefore, the CREST method and strategies allow comprehensive single-cell lineage analysis that could offer new insights into the molecular programs underlying neural specification.

Respiratory epithelial cell types, states and fates in the era of single-cell RNA-sequencing

Standalone and consortia-led single-cell atlases of healthy and diseased human airways generated with single-cell RNA-sequencing (scRNA-seq) have ushered in a new era in respiratory research. Numerous discoveries, including the pulmonary ionocyte, potentially novel cell fates, and a diversity of cell states among common and rare epithelial cell types have highlighted the extent of cellular heterogeneity and plasticity in the respiratory tract. scRNA-seq has also played a pivotal role in our understanding of host-virus interactions in coronavirus disease 2019 (COVID-19). However, as our ability to generate large quantities of scRNA-seq data increases, along with a growing number of scRNA-seq protocols and data analysis methods, new challenges related to the contextualisation and downstream applications of insights are arising. Here, we review the fundamental concept of cellular identity from the perspective of single-cell transcriptomics in the respiratory context, drawing attention to the need to generate reference annotations and to standardise the terminology used in literature. Findings about airway epithelial cell types, states and fates obtained from scRNA-seq experiments are compared and contrasted with information accumulated through the use of conventional methods. This review attempts to discuss major opportunities and to outline some of the key limitations of the modern-day scRNA-seq that need to be addressed to enable efficient and meaningful integration of scRNA-seq data from different platforms and studies, with each other as well as with data from other high-throughput sequencing-based genomic, transcriptomic and epigenetic analyses.

Cellular barcoding: From developmental tracing to anti-tumor drug discovery

Clonal evolution has gained immense attention in explaining cancer cell status, history, and fate during cancer progression. Current single-cell or spatial transcriptome technologies have broadened our understanding of various mechanisms underlying cancer initiation, relapse, and drug resistance. However, technical challenges still hinder a better understanding of the dynamics of distinctive phenotypic states and abnormal trajectories from normal physiological transition to malignant stages. Cellular barcoding enabled lineage tracing on parallelly massive cells at single-cell resolution through different mechanisms lately, enabling new insights into exploring developmental trajectories, cancer progression, and targeted therapies. This review summarizes the latest noteworthy and robust strategies for different types of cellular barcodes. To introduce the major characteristics, advantages and limitations of these different strategies, this review will further guide in choosing or improving cellular barcoding technologies and their applications in cancer research.

Formation of the Heart: Defining Cardiomyocyte Progenitors at Single-Cell Resolution

PURPOSE OF REVIEW

Formation of the heart requires the coordinated addition of multiple progenitor sources which have undergone different pathways of specification and differentiation. In this review, I aim to put into context how recent studies defining cardiac progenitor heterogeneity build on our understanding of early heart development and also discuss the questions raised by this new insight.

RECENT FINDINGS

With the development of sequencing technologies and imaging approaches, it has been possible to define, at high temporal resolution, the molecular profile and anatomical location of cardiac progenitors at the single-cell level, during the formation of the mammalian heart. Given the recent progress in our understanding of early heart development and technical advances in high-resolution time-lapse imaging and lineage analysis, we are now in a position of great potential, allowing us to resolve heart formation at previously impossible levels of detail. Understanding how this essential organ forms not only addresses questions of fundamental biological significance but also provides a blueprint for strategies to both treat and model heart disease.

Best Practices in Designing, Sequencing, and Identifying Random DNA Barcodes

Random DNA barcodes are a versatile tool for tracking cell lineages, with applications ranging from development to cancer to evolution. Here, we review and critically evaluate barcode designs as well as methods of barcode sequencing and initial processing of barcode data. We first demonstrate how various barcode design decisions affect data quality and propose a new design that balances all considerations that we are currently aware of. We then discuss various options for the preparation of barcode sequencing libraries, including inline indices and Unique Molecular Identifiers (UMIs). Finally, we test the performance of several established and new bioinformatic pipelines for the extraction of barcodes from raw sequencing reads and for error correction. We find that both alignment and regular expression-based approaches work well for barcode extraction, and that error-correction pipelines designed specifically for barcode data are superior to generic ones. Overall, this review will help researchers to approach their barcoding experiments in a deliberate and systematic way.

Toward a systems-level probing of tumor clonality

Cancer has been described as a genetic disease that clonally evolves in the face of selective pressures imposed by cell-intrinsic and extrinsic factors. Although classical models based on genetic data predominantly propose Darwinian mechanisms of cancer evolution, recent single-cell profiling of cancers has described unprecedented heterogeneity in tumors providing support for alternative models of branched and neutral evolution through both genetic and non-genetic mechanisms. Emerging evidence points to a complex interplay between genetic, non-genetic, and extrinsic environmental factors in shaping the evolution of tumors. In this perspective, we briefly discuss the role of cell-intrinsic and extrinsic factors that shape clonal behaviors during tumor progression, metastasis, and drug resistance. Taking examples of pre-malignant states associated with hematological malignancies and esophageal cancer, we discuss recent paradigms of tumor evolution and prospective approaches to further enhance our understanding of this spatiotemporally regulated process.

A new miniMOS tool kit capable of visualizing single copy insertion in C. elegans

The miniMOS technique has been widely used in the C. elegans community to generate single copy insertions. A worm is considered as a potential insertion candidate if it is resistant to G418 antibiotics and does not express a co-injected fluorescence marker. If the expression of the extrachromosomal array is very low, it is possible for a worm to be mistakenly identified as a miniMOS candidate, as this low expression level can still confer resistance to G418 without producing a detectable fluorescence signal from the co-injection marker. This may increase the workload for identifying the insertion locus in the subsequent steps. In the present study, we modified the plasmid platform for miniMOS insertion by incorporating a myo-2 promoter-driven TagRFP or a ubiquitous H2B::GFP expression cassette into the targeting vector and introducing two loxP sites flanking the selection cassettes. Based on this new miniMOS tool kit, the removable fluorescence reporters can be used to visualize the single copy insertions, greatly reducing insertion locus identification efforts. In our experience, this new platform greatly facilitates the isolation of the miniMOS mutants.

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.

FSF-GA: A Feature Selection Framework for Phenotype Prediction Using Genetic Algorithms

(1) Background: Phenotype prediction is a pivotal task in genetics in order to identify how genetic factors contribute to phenotypic differences. This field has seen extensive research, with numerous methods proposed for predicting phenotypes. Nevertheless, the intricate relationship between genotypes and complex phenotypes, including common diseases, has resulted in an ongoing challenge to accurately decipher the genetic contribution. (2) Results: In this study, we propose a novel feature selection framework for phenotype prediction utilizing a genetic algorithm (FSF-GA) that effectively reduces the feature space to identify genotypes contributing to phenotype prediction. We provide a comprehensive vignette of our method and conduct extensive experiments using a widely used yeast dataset. (3) Conclusions: Our experimental results show that our proposed FSF-GA method delivers comparable phenotype prediction performance as compared to baseline methods, while providing features selected for predicting phenotypes. These selected feature sets can be used to interpret the underlying genetic architecture that contributes to phenotypic variation.

Zebrafish cutaneous injury models reveal that Langerhans cells engulf axonal debris in adult epidermis

Somatosensory neurons extend enormous peripheral axons to the skin, where they detect diverse environmental stimuli. Somatosensory peripheral axons are easily damaged due to their small caliber and superficial location. Axonal damage results in Wallerian degeneration, creating vast quantities of cellular debris that phagocytes must remove to maintain organ homeostasis. The cellular mechanisms that ensure efficient clearance of axon debris from stratified adult skin are unknown. Here, we established zebrafish scales as a tractable model to study axon degeneration in the adult epidermis. Using this system, we demonstrated that skin-resident immune cells known as Langerhans cells engulf the majority of axon debris. In contrast to immature skin, adult keratinocytes did not significantly contribute to debris removal, even in animals lacking Langerhans cells. Our study establishes a powerful new model for studying Wallerian degeneration and identifies a new function for Langerhans cells in maintenance of adult skin homeostasis following injury. These findings have important implications for pathologies that trigger somatosensory axon degeneration.

Targeting the metastatic niche: Single-cell lineage tracing in prime time

Identification of actionable drug targets remains a rate-limiting step of, and one of the most prominent barriers to successful drug development for metastatic cancers. CRISPR-Cas9, a tool for making targeted genomic edits, has given rise to various novel applications that have greatly accelerated discovery in developmental biology. Recent work has coupled a CRISPR-Cas9-based lineage tracing platform with single-cell transcriptomics in the unexplored context of cancer metastasis. In this perspective, we briefly reflect on the development of these distinct technological advances and the process by which they have become integrated. We also highlight the importance of single-cell lineage tracing in oncology drug development and suggest the profound capacity of a high-resolution, computational approach to reshape cancer drug discovery by enabling identification of novel metastasis-specific drug targets and mechanisms of resistance.

Geographically weighted linear combination test for gene-set analysis of a continuous spatial phenotype as applied to intratumor heterogeneity

Background: The impact of gene-sets on a spatial phenotype is not necessarily uniform across different locations of cancer tissue. This study introduces a computational platform, GWLCT, for combining gene set analysis with spatial data modeling to provide a new statistical test for location-specific association of phenotypes and molecular pathways in spatial single-cell RNA-seq data collected from an input tumor sample.

Methods: The main advantage of GWLCT consists of an analysis beyond global significance, allowing the association between the gene-set and the phenotype to vary across the tumor space. At each location, the most significant linear combination is found using a geographically weighted shrunken covariance matrix and kernel function. Whether a fixed or adaptive bandwidth is determined based on a cross-validation cross procedure. Our proposed method is compared to the global version of linear combination test (LCT), bulk and random-forest based gene-set enrichment analyses using data created by the Visium Spatial Gene Expression technique on an invasive breast cancer tissue sample, as well as 144 different simulation scenarios.

Results: In an illustrative example, the new geographically weighted linear combination test, GWLCT, identifies the cancer hallmark gene-sets that are significantly associated at each location with the five spatially continuous phenotypic contexts in the tumors defined by different well-known markers of cancer-associated fibroblasts. Scan statistics revealed clustering in the number of significant gene-sets. A spatial heatmap of combined significance over all selected gene-sets is also produced. Extensive simulation studies demonstrate that our proposed approach outperforms other methods in the considered scenarios, especially when the spatial association increases.

Conclusion: Our proposed approach considers the spatial covariance of gene expression to detect the most significant gene-sets affecting a continuous phenotype. It reveals spatially detailed information in tissue space and can thus play a key role in understanding the contextual heterogeneity of cancer cells.

BdLT-Seq as a barcode decay-based method to unravel lineage-linked transcriptome plasticity

Cell plasticity is a core biological process underlying a myriad of molecular and cellular events taking place throughout organismal development and evolution. It has been postulated that cellular systems thrive to balance the organization of meta-stable states underlying this phenomenon, thereby maintaining a degree of populational homeostasis compatible with an ever-changing environment and, thus, life. Notably, albeit circumstantial evidence has been gathered in favour of the latter conceptual framework, a direct observation of meta-state dynamics and the biological consequences of such a process in generating non-genetic clonal diversity and divergent phenotypic output remains largely unexplored. To fill this void, here we develop a lineage-tracing technology termed Barcode decay Lineage Tracing-Seq. BdLT-Seq is based on episome-encoded molecular identifiers that, supported by the dynamic decay of the tracing information upon cell division, ascribe directionality to a cell lineage tree whilst directly coupling non-genetic molecular features to phenotypes in comparable genomic landscapes. We show that cell transcriptome states are both inherited, and dynamically reshaped following constrained rules encoded within the cell lineage in basal growth conditions, upon oncogene activation and throughout the process of reversible resistance to therapeutic cues thus adjusting phenotypic output leading to intra-clonal non-genetic diversity.

Integration of quantitative methods and mathematical approaches for the modeling of cancer cell proliferation dynamics

Physiological processes rely on the control of cell proliferation, and the dysregulation of these processes underlies various pathological conditions, including cancer. Mathematical modeling can provide new insights into the complex regulation of cell proliferation dynamics. In this review, we first examine quantitative experimental approaches for measuring cell proliferation dynamics in vitro and compare the various types of data that can be obtained in these settings. We then explore the toolbox of common mathematical modeling frameworks that can describe cell behavior, dynamics, and interactions of proliferation. We discuss how these wet-laboratory studies may be integrated with different mathematical modeling approaches to aid the interpretation of the results and to enable the prediction of cell behaviors, specifically in the context of cancer.

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.

Genetically engineered zebrafish as models of skeletal development and regeneration

Zebrafish (Danio rerio) are aquatic vertebrates with significant homology to their terrestrial counterparts. While zebrafish have a centuries-long track record in developmental and regenerative biology, their utility has grown exponentially with the onset of modern genetics. This is exemplified in studies focused on skeletal development and repair. Herein, the numerous contributions of zebrafish to our understanding of the basic science of cartilage, bone, tendon/ligament, and other skeletal tissues are described, with a particular focus on applications to development and regeneration. We summarize the genetic strengths that have made the zebrafish a powerful model to understand skeletal biology. We also highlight the large body of existing tools and techniques available to understand skeletal development and repair in the zebrafish and introduce emerging methods that will aid in novel discoveries in skeletal biology. Finally, we review the unique contributions of zebrafish to our understanding of regeneration and highlight diverse routes of repair in different contexts of injury. We conclude that zebrafish will continue to fill a niche of increasing breadth and depth in the study of basic cellular mechanisms of skeletal biology.

Updated Perspectives on Direct Vascular Cellular Reprogramming and Their Potential Applications in Tissue Engineered Vascular Grafts

Cardiovascular disease is a globally prevalent disease with far-reaching medical and socio-economic consequences. Although improvements in treatment pathways and revascularisation therapies have slowed disease progression, contemporary management fails to modulate the underlying atherosclerotic process and sustainably replace damaged arterial tissue. Direct cellular reprogramming is a rapidly evolving and innovative tissue regenerative approach that holds promise to restore functional vasculature and restore blood perfusion. The approach utilises cell plasticity to directly convert somatic cells to another cell fate without a pluripotent stage. In this narrative literature review, we comprehensively analyse and compare direct reprogramming protocols to generate endothelial cells, vascular smooth muscle cells and vascular progenitors. Specifically, we carefully examine the reprogramming factors, their molecular mechanisms, conversion efficacies and therapeutic benefits for each induced vascular cell. Attention is given to the application of these novel approaches with tissue engineered vascular grafts as a therapeutic and disease-modelling platform for cardiovascular diseases. We conclude with a discussion on the ethics of direct reprogramming, its current challenges, and future perspectives.

A Device for Isolation of Selected Single Adherent Cells

In recent years, extensive research has been focused on the field of single cell analysis. The isolation of single cells is the first step in this type of research. However, the techniques used for direct isolation and acquisition of single adherent cells are limited. Here, we present a method of obtaining selected single adherent cells using a separation device. Compared with other single cell isolation methods, this method has the advantages of simple operation, low cost, minimal cell damage, and preservation of cell morphology. Our methodology is, therefore, suitable for the collection of selected single adherent cells.

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.

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.

Recent advances in understanding cell types during human gastrulation

Gastrulation is a fundamental process during embryonic development, conserved across all multicellular animals [1]. In the majority of metazoans, gastrulation is characterised by large scale morphogenetic remodeling, leading to the conversion of an early pluripotent embryonic cell layer into the three primary 'germ layers': an outer ectoderm, inner endoderm and intervening mesoderm layer. The morphogenesis of these three layers of cells is closely coordinated with cellular diversification, laying the foundation for the generation of the hundreds of distinct specialized cell types in the animal body. The process of gastrulation has for a long time attracted tremendous attention in a broad range of experimental systems ranging from sponges to mice. In humans the process of gastrulation starts approximately 14 days after fertilization and continues for slightly over a week. However our understanding of this important process, as it pertains to human, is limited. Donations of human fetal material at these early stages are exceptionally rare, making it nearly impossible to study human gastrulation directly. Therefore, our understanding of human gastrulation is predominantly derived from animal models such as the mouse [2,3] and from studies of limited collections of fixed whole samples and histological sections of human gastrulae [4-7], some of which date back to over a century ago. More recently we have been gaining valuable molecular insights into human gastrulation using in vitro models of hESCs [8-12] and increasingly, in vitro cultured human and non-human primate embryos [13-16]. However, while methods have been developed to culture human embryos into this stage (and probably beyond), current ethical standards prohibit the culture of human embryos past 14 days again limiting our ability to experimentally probe human gastrulation. This review discusses recent molecular insights from the study of a rare CS 7 human gastrula obtained as a live sample and raises several questions arising from this recent study that it will be interesting to address in the future using emerging models of human gastrulation.

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.

Machine-learning-optimized Cas12a barcoding enables the recovery of single-cell lineages and transcriptional profiles

The development of CRISPR-based barcoding methods creates an exciting opportunity to understand cellular phylogenies. We present a compact, tunable, high-capacity Cas12a barcoding system called dual acting inverted site array (DAISY). We combined high-throughput screening and machine learning to predict and optimize the 60-bp DAISY barcode sequences. After optimization, top-performing barcodes had ∼10-fold increased capacity relative to the best random-screened designs and performed reliably across diverse cell types. DAISY barcode arrays generated ∼12 bits of entropy and ∼66,000 unique barcodes. Thus, DAISY barcodes-at a fraction of the size of Cas9 barcodes-achieved high-capacity barcoding. We coupled DAISY barcoding with single-cell RNA-seq to recover lineages and gene expression profiles from ∼47,000 human melanoma cells. A single DAISY barcode recovered up to ∼700 lineages from one parental cell. This analysis revealed heritable single-cell gene expression and potential epigenetic modulation of memory gene transcription. Overall, Cas12a DAISY barcoding is an efficient tool for investigating cell-state dynamics.

Origin and function of activated fibroblast states during zebrafish heart regeneration

The adult zebrafish heart has a high capacity for regeneration following injury. However, the composition of the regenerative niche has remained largely elusive. Here, we dissected the diversity of activated cell states in the regenerating zebrafish heart based on single-cell transcriptomics and spatiotemporal analysis. We observed the emergence of several transient cell states with fibroblast characteristics following injury, and we outlined the proregenerative function of collagen-12-expressing fibroblasts. To understand the cascade of events leading to heart regeneration, we determined the origin of these cell states by high-throughput lineage tracing. We found that activated fibroblasts were derived from two separate sources: the epicardium and the endocardium. Mechanistically, we determined Wnt signalling as a regulator of the endocardial fibroblast response. In summary, our work identifies specialized activated fibroblast cell states that contribute to heart regeneration, thereby opening up possible approaches to modulating the regenerative capacity of the vertebrate heart.

Single-cell RNA sequencing and spatiotemporal analysis of the regenerating zebrafish heart identify transient proregenerative fibroblast-like cells that are derived from the epicardium and the endocardium. Wnt signalling regulates the endocardial fibroblast response.

CoSpar identifies early cell fate biases from single-cell transcriptomic and lineage information

A goal of single-cell genome-wide profiling is to reconstruct dynamic transitions during cell differentiation, disease onset and drug response. Single-cell assays have recently been integrated with lineage tracing, a set of methods that identify cells of common ancestry to establish bona fide dynamic relationships between cell states. These integrated methods have revealed unappreciated cell dynamics, but their analysis faces recurrent challenges arising from noisy, dispersed lineage data. In this study, we developed coherent, sparse optimization (CoSpar) as a robust computational approach to infer cell dynamics from single-cell transcriptomics integrated with lineage tracing. Built on assumptions of coherence and sparsity of transition maps, CoSpar is robust to severe downsampling and dispersion of lineage data, which enables simpler experimental designs and requires less calibration. In datasets representing hematopoiesis, reprogramming and directed differentiation, CoSpar identifies early fate biases not previously detected, predicting transcription factors and receptors implicated in fate choice. Documentation and detailed examples for common experimental designs are available at https://cospar.readthedocs.io/ .

Hematopoietic Stem and Progenitor Cells (HSPCs) and Hematopoietic Microenvironment: Molecular and Bioinformatic Studies of the Zebrafish Models

Hematopoietic stem cells (HSCs) reside in a specialized microenvironment in a peculiar anatomic location which regulates the maintenance of stem cells and controls its functions. Recent scientific progress in experimental technologies have enabled the specific detection of epigenetic factors responsible for the maintenance and quiescence of the hematopoietic niche, which has improved our knowledge of regulatory mechanisms. The aberrant role of RNA-binding proteins and their impact on the disruption of stem cell biology have been reported by a number of recent studies. Despite recent modernization in hematopoietic microenvironment research avenues, our comprehension of the signaling mechanisms and interactive pathways responsible for integration of the hematopoietic niche is still limited. In the past few decades, zebrafish usage with regards to exploratory studies of the hematopoietic niche has expanded our knowledge for deeper understanding of novel cellular interactions. This review provides an update on the functional roles of different genetic and epigenetic factors and molecular signaling events at different sections of the hematopoietic microenvironment. The explorations of different molecular approaches and interventions of latest web-based tools being used are also outlined. This will help us to get more mechanistic insights and develop therapeutic options for the malignancies.

Toward the dissection of hematopoietic stem cell fates and their determinants

Hematopoietic stem cell (HSC) functions have long been difficult to study under physiological conditions. Recently, genetic in vivo approaches have been developed for lineage tracing of differentiating progeny emerging from HSC over time (output), and for high-resolution, endogenous barcoding to uncover the lineages that HSC contribute to (fate). Such fate measurements have in principle led to the recognition of three major fate groups of HSC: multilineage, myelo-erythroid-restricted, and inactive, that is, no or no known progeny, in addition to a minor group of megakaryocyte-restricted HSC. The most recent RNA-barcoding experiments have begun to directly link fate measurements with transcriptome reading in HSC clones and single HSC, which yielded insights into transcriptional signatures associated with fate patterns. Here, we discuss these findings in light of the structure of the hematopoietic differentiation hierarchy, and we provide an outlook on strategies to dissect molecular determinants of HSC fates.

Advances and applications of single-cell omics technologies in plant research

Single-cell sequencing approaches reveal the intracellular dynamics of individual cells and answer biological questions with high-dimensional catalogs of millions of cells, including genomics, transcriptomics, chromatin accessibility, epigenomics, and proteomics data across species. These emerging yet thriving technologies have been fully embraced by the field of plant biology, with a constantly expanding portfolio of applications. Here, we introduce the current technical advances used for single-cell omics, especially single-cell genome and transcriptome sequencing. Firstly, we overview methods for protoplast and nucleus isolation and genome and transcriptome amplification. Subsequently, we use well-executed benchmarking studies to highlight advances made through the application of single-cell omics techniques. Looking forward, we offer a glimpse of additional hurdles and future opportunities that will introduce broad adoption of single-cell sequencing with revolutionary perspectives in plant biology.

Developmental dynamics of two bipotent thymic epithelial progenitor types

T cell development in the thymus is essential for cellular immunity and depends on the organotypic thymic epithelial microenvironment. In comparison with other organs, the size and cellular composition of the thymus are unusually dynamic, as exemplified by rapid growth and high T cell output during early stages of development, followed by a gradual loss of functional thymic epithelial cells and diminished naive T cell production with age1-10. Single-cell RNA sequencing (scRNA-seq) has uncovered an unexpected heterogeneity of cell types in the thymic epithelium of young and aged adult mice11-18; however, the identities and developmental dynamics of putative pre- and postnatal epithelial progenitors have remained unresolved1,12,16,17,19-27. Here we combine scRNA-seq and a new CRISPR-Cas9-based cellular barcoding system in mice to determine qualitative and quantitative changes in the thymic epithelium over time. This dual approach enabled us to identify two principal progenitor populations: an early bipotent progenitor type biased towards cortical epithelium and a postnatal bipotent progenitor population biased towards medullary epithelium. We further demonstrate that continuous autocrine provision of Fgf7 leads to sustained expansion of thymic microenvironments without exhausting the epithelial progenitor pools, suggesting a strategy to modulate the extent of thymopoietic activity.

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.

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.

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.

Single cell lineage reconstruction using distance-based algorithms and the R package, DCLEAR

Background

DCLEAR is an R package used for single cell lineage reconstruction. The advances of CRISPR-based gene editing technologies have enabled the prediction of cell lineage trees based on observed edited barcodes from each cell. However, the performance of existing reconstruction methods of cell lineage trees was not accessed until recently. In response to this problem, the Allen Institute hosted the Cell Lineage Reconstruction Dream Challenge in 2020 to crowdsource relevant knowledge from across the world. Our team won sub-challenges 2 and 3 in the challenge competition.

Results

The DCLEAR package contained the R codes, which was submitted in response to sub-challenges 2 and 3. Our method consists of two steps: (1) distance matrix estimation and (2) the tree reconstruction from the distance matrix. We proposed two novel methods for distance matrix estimation as outlined in the DCLEAR package. Using our method, we find that two of the more sophisticated distance methods display a substantially improved level of performance compared to the traditional Hamming distance method. DCLEAR is open source and freely available from R CRAN and from under the GNU General Public License, version 3.

Conclusions

DCLEAR is a powerful resource for single cell lineage reconstruction.

Genetic mosaicism in the human brain: from lineage tracing to neuropsychiatric disorders

Genetic mosaicism is the result of the accumulation of somatic mutations in the human genome starting from the first postzygotic cell generation and continuing throughout the whole life of an individual. The rapid development of next-generation and single-cell sequencing technologies is now allowing the study of genetic mosaicism in normal tissues, revealing unprecedented insights into their clonal architecture and physiology. The somatic variant repertoire of an adult human neuron is the result of somatic mutations that accumulate in the brain by different mechanisms and at different rates during development and ageing. Non-pathogenic developmental mutations function as natural barcodes that once identified in deep bulk or single-cell sequencing can be used to retrospectively reconstruct human lineages. This approach has revealed novel insights into the clonal structure of the human brain, which is a mosaic of clones traceable to the early embryo that contribute differentially to the brain and distinct areas of the cortex. Some of the mutations happening during development, however, have a pathogenic effect and can contribute to some epileptic malformations of cortical development and autism spectrum disorder. In this Review, we discuss recent findings in the context of genetic mosaicism and their implications for brain development and disease.

Characterization of Pipefish Immune Cell Populations Through Single-Cell Transcriptomics

Teleost adaptive immune systems have evolved with more flexibility than previously assumed. A particularly enigmatic system to address immune system modifications in the evolutionary past is represented by the Syngnathids, the family of pipefishes, seahorses and seadragons. These small fishes with their unique male pregnancy have lost the spleen as an important immune organ as well as a functional major histocompatibility class II (MHC II) pathway. How these evolutionary changes have impacted immune cell population dynamics have up to this point remained unexplored. Here, we present the first immune cell repertoire characterization of a syngnathid fish (Syngnathus typhle) using single-cell transcriptomics. Gene expression profiles of individual cells extracted from blood and head-kidney clustered in twelve putative cell populations with eight belonging to those with immune function. Upregulated cell marker genes identified in humans and teleosts were used to define cell clusters. While the suggested loss of CD4+ T-cells accompanied the loss of the MHC II pathway was supported, the upregulation of specific subtype markers within the T-cell cluster indicates subpopulations of regulatory T-cells (il2rb) and cytotoxic T-cells (gzma). Utilizing single-cell RNA sequencing this report is the first to characterize immune cell populations in syngnathids and provides a valuable foundation for future cellular classification and experimental work within the lineage.

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.

Clonal relations in the mouse brain revealed by single-cell and spatial transcriptomics

The mammalian brain contains many specialized cells that develop from a thin sheet of neuroepithelial progenitor cells. Single-cell transcriptomics revealed hundreds of molecularly diverse cell types in the nervous system, but the lineage relationships between mature cell types and progenitor cells are not well understood. Here we show in vivo barcoding of early progenitors to simultaneously profile cell phenotypes and clonal relations in the mouse brain using single-cell and spatial transcriptomics. By reconstructing thousands of clones, we discovered fate-restricted progenitor cells in the mouse hippocampal neuroepithelium and show that microglia are derived from few primitive myeloid precursors that massively expand to generate widely dispersed progeny. We combined spatial transcriptomics with clonal barcoding and disentangled migration patterns of clonally related cells in densely labeled tissue sections. Our approach enables high-throughput dense reconstruction of cell phenotypes and clonal relations at the single-cell and tissue level in individual animals and provides an integrated approach for understanding tissue architecture.

Ratz et al. present an easy-to-use method to barcode progenitor cells, enabling profiling of cell phenotypes and clonal relations using single-cell and spatial transcriptomics, providing an integrated approach for understanding brain architecture.

Temporal signaling, population control, and information processing through chromatin-mediated gene regulation

Chromatin regulation is a key pathway cells use to regulate gene expression in response to temporal stimuli, and is becoming widely used as a platform for synthetic biology applications. Here, we build a mathematical framework for analyzing the response of genetic circuits containing chromatin regulators to temporal signals in mammalian cell populations. Chromatin regulators can silence genes in an all-or-none fashion at the single-cell level, with individual cells stochastically transitioning between active, reversibly silent, and irreversibly silent gene states at constant rates over time. We integrate this mode of regulation with classical gene regulatory motifs, such as autoregulatory and incoherent feedforward loops, to determine the types of responses achievable with duration-dependent signaling. We demonstrate that repressive regulators without long-term epigenetic memory can filter out high frequency noise, and as part of an autoregulatory loop can precisely tune the fraction of cells in a population that expresses a gene of interest. Additionally, we find that repressive regulators with epigenetic memory can sum up and encode the total duration of their recruitment in the fraction of cells irreversibly silenced and, when included in a feed forward loop, enable perfect adaptation. Last, we use an information theoretic approach to show that all-or-none stochastic silencing can be used by populations to transmit information reliably and with high fidelity even in very simple genetic circuits. Altogether, we show that chromatin-mediated gene control enables a repertoire of complex cell population responses to temporal signals and can transmit higher information levels than previously measured in gene regulation.

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.

Spatial components of molecular tissue biology

Methods for profiling RNA and protein expression in a spatially resolved manner are rapidly evolving, making it possible to comprehensively characterize cells and tissues in health and disease. To maximize the biological insights obtained using these techniques, it is critical to both clearly articulate the key biological questions in spatial analysis of tissues and develop the requisite computational tools to address them. Developers of analytical tools need to decide on the intrinsic molecular features of each cell that need to be considered, and how cell shape and morphological features are incorporated into the analysis. Also, optimal ways to compare different tissue samples at various length scales are still being sought. Grouping these biological problems and related computational algorithms into classes across length scales, thus characterizing common issues that need to be addressed, will facilitate further progress in spatial transcriptomics and proteomics.

Temporal modelling using single-cell transcriptomics

Methods for profiling genes at the single-cell level have revolutionized our ability to study several biological processes and systems including development, differentiation, response programmes and disease progression. In many of these studies, cells are profiled over time in order to infer dynamic changes in cell states and types, sets of expressed genes, active pathways and key regulators. However, time-series single-cell RNA sequencing (scRNA-seq) also raises several new analysis and modelling issues. These issues range from determining when and how deep to profile cells, linking cells within and between time points, learning continuous trajectories, and integrating bulk and single-cell data for reconstructing models of dynamic networks. In this Review, we discuss several approaches for the analysis and modelling of time-series scRNA-seq, highlighting their steps, key assumptions, and the types of data and biological questions they are most appropriate for.

RNA-Seq of Single Fish Cells - Seeking Out the Leukocytes Mediating Immunity in Teleost Fishes

The immune system is a complex and sophisticated biological system, spanning multiple levels of complexity, from the molecular level to that of tissue. Our current understanding of its function and complexity, of the heterogeneity of leukocytes, is a result of decades of concentrated efforts to delineate cellular markers using conventional methods of antibody screening and antigen identification. In mammalian models, this led to in-depth understanding of individual leukocyte subsets, their phenotypes, and their roles in health and disease. The field was further propelled forward by the development of single-cell (sc) RNA-seq technologies, offering an even broader and more integrated view of how cells work together to generate a particular response. Consequently, the adoption of scRNA-seq revealed the unexpected plasticity and heterogeneity of leukocyte populations and shifted several long-standing paradigms of immunology. This review article highlights the unprecedented opportunities offered by scRNA-seq technology to unveil the individual contributions of leukocyte subsets and their crosstalk in generating the overall immune responses in bony fishes. Single-cell transcriptomics allow identifying unseen relationships, and formulating novel hypotheses tailored for teleost species, without the need to rely on the limited number of fish-specific antibodies and pre-selected markers. Several recent studies on single-cell transcriptomes of fish have already identified previously unnoticed expression signatures and provided astonishing insights into the diversity of teleost leukocytes and the evolution of vertebrate immunity. Without a doubt, scRNA-seq in tandem with bioinformatics tools and state-of-the-art methods, will facilitate studying the teleost immune system by not only defining key markers, but also teaching us about lymphoid tissue organization, development/differentiation, cell-cell interactions, antigen receptor repertoires, states of health and disease, all across time and space in fishes. These advances will invite more researchers to develop the tools necessary to explore the immunology of fishes, which remain non-conventional animal models from which we have much to learn.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

Towards Tabula Gallus

The Tabula Gallus is a proposed project that aims to create a map of every cell type in the chicken body and chick embryos. Chickens (Gallus gallus) are one of the most recognized model animals that recapitulate the development and physiology of mammals. The Tabula Gallus will generate a compendium of single-cell transcriptome data from Gallus gallus, characterize each cell type, and provide tools for the study of the biology of this species, similar to other ongoing cell atlas projects (Tabula Muris and Tabula Sapiens/Human Cell Atlas for mice and humans, respectively). The Tabula Gallus will potentially become an international collaboration between many researchers. This project will be useful for the basic scientific study of Gallus gallus and other birds (e.g., cell biology, molecular biology, developmental biology, neuroscience, physiology, oncology, virology, behavior, ecology, and evolution). It will eventually be beneficial for a better understanding of human health and diseases.