Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo

Small Noncoding RNA Signatures for Determining the Developmental Potential of an Embryo at the Morula Stage

As part of the optimization of assisted reproductive technology programs, the aim of the study was to identify key small noncoding RNA (sncRNA) molecules that participate in maternal-to-zygotic transition and determine development potential and competence to form a healthy fetus. Small RNA deep sequencing followed by quantitative real-time RT-PCR was used to profile sncRNAs in 50 samples of spent culture medium from morula with different development potentials (no potential (degradation/developmental arrest), low potential (poor-quality blastocyst), and high potential (good/excellent quality blastocyst capable of implanting and leading to live birth)) obtained from 27 subfertile couples who underwent in vitro fertilization. We have shown that the quality of embryos at the morula stage is determined by secretion/uptake rates of certain sets of piRNAs and miRNAs, namely hsa_piR_011291, hsa_piR_019122, hsa_piR_001311, hsa_piR_015026, hsa_piR_015462, hsa_piR_016735, hsa_piR_019675, hsa_piR_020381, hsa_piR_020485, hsa_piR_004880, hsa_piR_000807, hsa-let-7b-5p, and hsa-let-7i-5p. Predicted gene targets of these sncRNAs included those globally decreased at the 8-cell–morula–blastocyst stage and critical to early embryo development. We show new original data on sncRNA profiling in spent culture medium from morula with different development potential. Our findings provide a view of a more complex network that controls human embryogenesis at the pre-implantation stage. Further research is required using reporter analysis to experimentally confirm interactions between identified sncRNA/gene target pairs.

A human cell atlas of fetal gene expression

The gene expression program underlying the specification of human cell types is of fundamental interest. We generated human cell atlases of gene expression and chromatin accessibility in fetal tissues. For gene expression, we applied three-level combinatorial indexing to >110 samples representing 15 organs, ultimately profiling ~4 million single cells. We leveraged the literature and other atlases to identify and annotate hundreds of cell types and subtypes, both within and across tissues. Our analyses focused on organ-specific specializations of broadly distributed cell types (such as blood, endothelial, and epithelial), sites of fetal erythropoiesis (which notably included the adrenal gland), and integration with mouse developmental atlases (such as conserved specification of blood cells). These data represent a rich resource for the exploration of in vivo human gene expression in diverse tissues and cell types.

Single-cell transcriptome profiling of an adult human cell atlas of 15 major organs

BACKGROUND

As core units of organ tissues, cells of various types play their harmonious rhythms to maintain the homeostasis of the human body. It is essential to identify the characteristics of cells in human organs and their regulatory networks for understanding the biological mechanisms related to health and disease. However, a systematic and comprehensive single-cell transcriptional profile across multiple organs of a normal human adult is missing.

RESULTS

We perform single-cell transcriptomes of 84,363 cells derived from 15 tissue organs of one adult donor and generate an adult human cell atlas. The adult human cell atlas depicts 252 subtypes of cells, including major cell types such as T, B, myeloid, epithelial, and stromal cells, as well as novel COCH+ fibroblasts and FibSmo cells, each of which is distinguished by multiple marker genes and transcriptional profiles. These collectively contribute to the heterogeneity of major human organs. Moreover, T cell and B cell receptor repertoire comparisons and trajectory analyses reveal direct clonal sharing of T and B cells with various developmental states among different tissues. Furthermore, novel cell markers, transcription factors, and ligand-receptor pairs are identified with potential functional regulations in maintaining the homeostasis of human cells among tissues.

CONCLUSIONS

The adult human cell atlas reveals the inter- and intra-organ heterogeneity of cell characteristics and provides a useful resource in uncovering key events during the development of human diseases in the context of the heterogeneity of cells and organs.

Exploring Additional Valuable Information From Single-Cell RNA-Seq Data

Single-cell RNA-seq (scRNA-seq) technologies are broadly applied to dissect the cellular heterogeneity and expression dynamics, providing unprecedented insights into single-cell biology. Most of the scRNA-seq studies mainly focused on the dissection of cell types/states, developmental trajectory, gene regulatory network, and alternative splicing. However, besides these routine analyses, many other valuable scRNA-seq investigations can be conducted. Here, we first review cell-to-cell communication exploration, RNA velocity inference, identification of large-scale copy number variations and single nucleotide changes, and chromatin accessibility prediction based on single-cell transcriptomics data. Next, we discuss the identification of novel genes/transcripts through transcriptome reconstruction approaches, as well as the profiling of long non-coding RNAs and circular RNAs. Additionally, we survey the integration of single-cell and bulk RNA-seq datasets for deconvoluting the cell composition of large-scale bulk samples and linking single-cell signatures to patient outcomes. These additional analyses could largely facilitate corresponding basic science and clinical applications.

Pax9 is essential for granulopoiesis but dispensable for erythropoiesis in zebrafish

Paired Box (Pax) gene family, a group of transcription regulators have been implicated in diverse physiological processes. However, their role during hematopoiesis which generate a plethora of blood cells remains largely unknown. Using a previously reported single cell transcriptomics data, we analyzed the expression of individual Pax family members in hematopoietic cells in zebrafish. We have identified that Pax9, which is an essential regulator for odontogenesis and palatogenesis, is selectively localized within a single cluster of the hematopoietic lineage. To further analyze the function of Pax9 in hematopoiesis, we generated two independent pax9 knock-out mutants using the CRISPR-Cas9 technique. We found that Pax9 appears to be an essential regulator for granulopoiesis but dispensable for erythropoiesis during development, as lack of pax9 selectively decreased the number of neutrophils with a concomitant decrease in the expression level of neutrophil markers. In addition, embryos, where pax9 was functionally disrupted by injecting morpholinos, failed to increase the number of neutrophils in response to pathogenic bacteria, suggesting that Pax9 is not only essential for developmental granulopoiesis but also emergency granulopoiesis. Due to the inability to initiate emergency granulopoiesis, innate immune responses were severely compromised in pax9 morpholino-mediated embryos, increasing their susceptibility and mortality. Taken together, our data indicate that Pax9 is essential for granulopoiesis and promotes innate immunity in zebrafish larvae.

Defining the Teratoma as a Model for Multi-lineage Human Development

We propose that the teratoma, a recognized standard for validating pluripotency in stem cells, could be a promising platform for studying human developmental processes. Performing single-cell RNA sequencing (RNA-seq) of 179,632 cells across 23 teratomas from 4 cell lines, we found that teratomas reproducibly contain approximately 20 cell types across all 3 germ layers, that inter-teratoma cell type heterogeneity is comparable with organoid systems, and teratoma gut and brain cell types correspond well to similar fetal cell types. Furthermore, cellular barcoding confirmed that injected stem cells robustly engraft and contribute to all lineages. Using pooled CRISPR-Cas9 knockout screens, we showed that teratomas can enable simultaneous assaying of the effects of genetic perturbations across all germ layers. Additionally, we demonstrated that teratomas can be sculpted molecularly via microRNA (miRNA)-regulated suicide gene expression to enrich for specific tissues. Taken together, teratomas are a promising platform for modeling multi-lineage development, pan-tissue functional genetic screening, and tissue engineering.

Analyses of Avascular Mutants Reveal Unique Transcriptomic Signature of Non-conventional Endothelial Cells

Endothelial cells appear to emerge from diverse progenitors. However, to which extent their developmental origin contributes to define their cellular and molecular characteristics remains largely unknown. Here, we report that a subset of endothelial cells that emerge from the tailbud possess unique molecular characteristics that set them apart from stereotypical lateral plate mesoderm (LPM)-derived endothelial cells. Lineage tracing shows that these tailbud-derived endothelial cells arise at mid-somitogenesis stages, and surprisingly do not require Npas4l or Etsrp function, indicating that they have distinct spatiotemporal origins and are regulated by distinct molecular mechanisms. Microarray and single cell RNA-seq analyses reveal that somitogenesis- and neurogenesis-associated transcripts are over-represented in these tailbud-derived endothelial cells, suggesting that they possess a unique transcriptomic signature. Taken together, our results further reveal the diversity of endothelial cells with respect to their developmental origin and molecular properties, and provide compelling evidence that the molecular characteristics of endothelial cells may reflect their distinct developmental history.

Heterogeneity and 'memory' in stem cell populations

Modern single cell experiments have revealed unexpected heterogeneity in apparently functionally 'pure' cell populations. However, we are still lacking a conceptual framework to understand this heterogeneity. Here, we propose that cellular memories-changes in the molecular status of a cell in response to a stimulus, that modify the ability of the cell to respond to future stimuli-are an essential ingredient in any such theory. We illustrate this idea by considering a simple age-structured model of stem cell proliferation that takes account of mitotic memories. Using this model we argue that asynchronous mitosis generates heterogeneity that is central to stem cell population function. This model naturally explains why stem cell numbers increase through life, yet regenerative potency simultaneously declines.

Next-Generation Lineage Tracing and Fate Mapping to Interrogate Development

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

Hierarchical analysis of ontogenetic time to describe heterochrony and taxonomy of developmental stages

Even though an accurate description of early life stages is available for some teleostean species in form of embryonic and post-embryonic developmental tables, there is poor overlap between species-specific staging vocabularies beyond the taxonomic family level. What is called "embryonic period", "larval period", "metamorphosis", or "juvenile" is anatomically different across teleostean families. This problem, already pointed out 50 years ago, challenges the consistency of developmental biology, embryology, systematics, and hampers an efficient aquaculture diversification. We propose a general solution by producing a proof-of-concept hierarchical analysis of ontogenetic time using a set of four freshwater species displaying strongly divergent reproductive traits. With a parsimony analysis of a matrix where "operational taxonomic units" are species at a given ontogenetic time segment and characters are organs or structures which are coded present or absent at this time, we show that the hierarchies obtained have both very high consistency and retention index, indicating that the ontogenetic time is correctly grasped through a hierarchical graph. This allows to formally detect developmental heterochronies and might provide a baseline to name early life stages for any set of species. The present method performs a phylogenetic segmentation of ontogenetic time, which can be correctly seen as depicting ontophylogenesis.

Uncovering cellular networks in branching morphogenesis using single-cell transcriptomics

Single-cell RNA-sequencing (scRNA-seq) and related technologies to identify cell types and measure gene expression in space, in time, and within lineages have multiplied rapidly in recent years. As these techniques proliferate, we are seeing an increase in their application to the study of developing tissues. Here, we focus on single-cell investigations of branching morphogenesis. Branched organs are highly complex but typically develop recursively, such that a given developmental stage theoretically contains the entire spectrum of cell identities from progenitor to terminally differentiated. Therefore, branched organs are a highly attractive system for study by scRNA-seq. First, we provide an update on advances in the field of scRNA-seq analysis, focusing on spatial transcriptomics, computational reconstruction of differentiation trajectories, and integration of scRNA-seq with lineage tracing. In addition, we discuss the possibilities and limitations for applying these techniques to studying branched organs. We then discuss exciting advances made using scRNA-seq in the study of branching morphogenesis and differentiation in mammalian organs, with emphasis on the lung, kidney, and mammary gland. We propose ways that scRNA-seq could be used to address outstanding questions in each organ. Finally, we highlight the importance of physical and mechanical signals in branching morphogenesis and speculate about how scRNA-seq and related techniques could be applied to study tissue morphogenesis beyond just differentiation.

Targeted Single-Cell RNA and DNA Sequencing With Fluorescence-Activated Droplet Merger

Analyzing every cell in a diverse sample provides insight into population-level heterogeneity, but abundant cell types dominate the analysis and rarer populations are scarcely represented in the data. To focus on specific cell types, the current paradigm is to physically isolate subsets of interest prior to analysis; however, it remains difficult to isolate and then single-cell sequence such populations because of compounding losses. Here, we describe an alternative approach that selectively merges cells with reagents to achieve enzymatic reactions without having to physically isolate cells. We apply this technique to perform single-cell transcriptome and genome sequencing of specific cell subsets. Our method for analyzing heterogeneous populations obviates the need for pre- or post-enrichment and simplifies single-cell workflows, making it useful for other applications in single-cell biology, combinatorial chemical synthesis, and drug screening.

A single-cell analysis of the molecular lineage of chordate embryogenesis

Progressive unfolding of gene expression cascades underlies diverse embryonic lineage development. Here, we report a single-cell RNA sequencing analysis of the complete and invariant embryonic cell lineage of the tunicate Ciona savignyi from fertilization to the onset of gastrulation. We reconstructed a developmental landscape of 47 cell types over eight cell cycles in the wild-type embryo and identified eight fate transformations upon fibroblast growth factor (FGF) inhibition. For most FGF-dependent asymmetric cell divisions, the bipotent mother cell displays the gene signature of the default daughter fate. In convergent differentiation of the two notochord lineages, we identified additional gene pathways parallel to the master regulator T/Brachyury Last, we showed that the defined Ciona cell types can be matched to E6.5-E8.5 stage mouse cell types and display conserved expression of limited number of transcription factors. This study provides a high-resolution single-cell dataset to understand chordate early embryogenesis and cell lineage differentiation.

Of Cell Shapes and Motion: The Physical Basis of Animal Cell Migration

Motile cells have developed a variety of migration modes relying on diverse traction-force-generation mechanisms. Before the behavior of intracellular components could be easily imaged, cell movements were mostly classified by different types of cellular shape dynamics. Indeed, even though some types of cells move without any significant change in shape, most cell propulsion mechanisms rely on global or local deformations of the cell surface. In this review, focusing mostly on metazoan cells, we discuss how different types of local and global shape changes underlie distinct migration modes. We then discuss mechanical differences between force-generation mechanisms and finish by speculating on how they may have evolved.

Genetically encoded thiol redox-sensors in the zebrafish model: lessons for embryonic development and regeneration

Important roles for reactive oxygen species (ROS) and redox signaling in embryonic development and regenerative processes are increasingly recognized. However, it is difficult to obtain information on spatiotemporal dynamics of ROS production and signaling in vivo. The zebrafish is an excellent model for in vivo bioimaging and possesses a remarkable regenerative capacity upon tissue injury. Here, we review data obtained in this model system with genetically encoded redox-sensors targeting H₂O₂ and glutathione redox potential. We describe how such observations have prompted insight into regulation and downstream effects of redox alterations during tissue differentiation, morphogenesis and regeneration. We also discuss the properties of the different sensors and their consequences for the interpretation of in vivo imaging results. Finally, we highlight open questions and additional research fields that may benefit from further application of such sensor systems in zebrafish models of development, regeneration and disease.

Emerging Roles of Single-Cell Multi-Omics in Studying Developmental Temporal Patterning

The complexity of brain structure and function is rooted in the precise spatial and temporal regulation of selective developmental events. During neurogenesis, both vertebrates and invertebrates generate a wide variety of specialized cell types through the expansion and specification of a restricted set of neuronal progenitors. Temporal patterning of neural progenitors rests on fine regulation between cell-intrinsic and cell-extrinsic mechanisms. The rapid emergence of high-throughput single-cell technologies combined with elaborate computational analysis has started to provide us with unprecedented biological insights related to temporal patterning in the developing central nervous system (CNS). Here, we present an overview of recent advances in Drosophila and vertebrates, focusing both on cell-intrinsic mechanisms and environmental influences. We then describe the various multi-omics approaches that have strongly contributed to our current understanding and discuss perspectives on the various -omics approaches that hold great potential for the future of temporal patterning research.

Impact of data preprocessing on cell-type clustering based on single-cell RNA-seq data

Background

Advances in single-cell RNA-seq technology have led to great opportunities for the quantitative characterization of cell types, and many clustering algorithms have been developed based on single-cell gene expression. However, we found that different data preprocessing methods show quite different effects on clustering algorithms. Moreover, there is no specific preprocessing method that is applicable to all clustering algorithms, and even for the same clustering algorithm, the best preprocessing method depends on the input data.

Results

We designed a graph-based algorithm, SC3-e, specifically for discriminating the best data preprocessing method for SC3, which is currently the most widely used clustering algorithm for single cell clustering. When tested on eight frequently used single-cell RNA-seq data sets, SC3-e always accurately selects the best data preprocessing method for SC3 and therefore greatly enhances the clustering performance of SC3.

Conclusion

The SC3-e algorithm is practically powerful for discriminating the best data preprocessing method, and therefore largely enhances the performance of cell-type clustering of SC3. It is expected to play a crucial role in the related studies of single-cell clustering, such as the studies of human complex diseases and discoveries of new cell types.

Single-Cell Genomics and Epigenomics: Technologies and Applications in Plants

The development of genomics and epigenomics has allowed rapid advances in our understanding of plant biology. However, conventional bulk analysis dilutes cell-specific information by providing only average information, thereby limiting the resolution of genomic and functional genomic studies. Recent advances in single-cell sequencing technology concerning genomics and epigenomics open new avenues to dissect cell heterogeneity in multiple biological processes. Recent applications of these approaches to plants have provided exciting insights into diverse biological questions. We highlight the methodologies underlying the current techniques of single-cell genomics and epigenomics before covering their recent applications, potential significance, and future perspectives in plant biology.

Combinatorial Action of Temporally Segregated Transcription Factors

Combinatorial action of transcription factors (TFs) with partially overlapping expression is a widespread strategy to generate novel gene-expression patterns and, thus, cellular diversity. Known mechanisms underlying combinatorial activity require co-expression of TFs within the same cell. Here, we describe the mechanism by which two TFs that are never co-expressed generate a new, intersectional expression pattern in C. elegans embryos: lineage-specific priming of a gene by a transiently expressed TF generates a unique intersection with a second TF acting on the same gene four cell divisions later; the second TF is expressed in multiple cells but only activates transcription in those where priming occurred. Early induction of active transcription is necessary and sufficient to establish a competent state, maintained by broadly expressed regulators in the absence of the initial trigger. We uncover additional cells diversified through this mechanism. Our findings define a mechanism for combinatorial TF activity with important implications for generation of cell-type diversity.

•

Lineage-specific priming enables asymmetric gene expression in L/R neuron pairs

•

Transient, lineage-specific TFs prime a locus for later activation by a bilateral TF

•

An early active transcriptional state is necessary and sufficient for priming

•

Maintenance of asymmetric primed state occurs in a symmetric regulatory environment

Methodologies for Following EMT In Vivo at Single Cell Resolution

An epithelial-mesenchymal transition (EMT) occurs in almost every metazoan embryo at the time mesoderm begins to differentiate. Several embryos have a long record as models for studying an EMT given that a known population of cells enters the EMT at a known time thereby enabling a detailed study of the process. Often, however, it is difficult to learn the molecular details of these model EMT systems because the transitioning cells are a minority of the population of cells in the embryo and in most cases there is an inability to isolate that population. Here we provide a method that enables an examination of genes expressed before, during, and after the EMT with a focus on just the cells that undergo the transition. Single cell RNA-seq (scRNA-seq) has advanced as a technology making it feasible to study the trajectory of gene expression specifically in the cells of interest, in vivo, and without the background noise of other cell populations. The sea urchin skeletogenic cells constitute only 5% of the total number of cells in the embryo yet with scRNA-seq it is possible to study the genes expressed by these cells without background noise. This approach, though not perfect, adds a new tool for uncovering the mechanism of EMT in this cell type.

Adult sox10+ Cardiomyocytes Contribute to Myocardial Regeneration in the Zebrafish

During heart regeneration in the zebrafish, fibrotic tissue is replaced by newly formed cardiomyocytes derived from preexisting ones. It is unclear whether the heart is composed of several cardiomyocyte populations bearing different capacity to replace lost myocardium. Here, using sox10 genetic fate mapping, we identify a subset of preexistent cardiomyocytes in the adult zebrafish heart with a distinct gene expression profile that expanded after cryoinjury. Genetic ablation of sox10⁺ cardiomyocytes impairs cardiac regeneration, revealing that these cells play a role in heart regeneration.

Label-free microendoscopy using a micro-needle imaging probe for in vivo deep tissue imaging

We report a label-free imaging method for microendoscopy that uses a needle-type imaging probe. We inserted a thin GRIN lens that had been attached to a fiber bundle into a medical-grade needle that was used as an imaging probe. The introduction of the needle probe into biological tissue allows for direct access to deep regions that we otherwise could not achieve because of the multiple light scattering. To minimize invasiveness, we introduced the illuminating probe on the tissue surface, using an oblique back-illumination configuration. We achieved three-dimensional depth imaging by changing the depth of penetration. Since only the imaging probe goes deep into the tissue while leaving the illumination channels outside, the achievable signal depends on the location of the illumination channels. We explored this point and investigated the optimal condition for the illumination distance in a systematic way. We also applied this method to ex vivo, as well as in vivo, imaging of a mouse brain, and confirmed that we had visualized the microvasculature embedded deep within the brain.

An improved zebrafish transcriptome annotation for sensitive and comprehensive detection of cell type-specific genes

The zebrafish is ideal for studying embryogenesis and is increasingly applied to model human disease. In these contexts, RNA-sequencing (RNA-seq) provides mechanistic insights by identifying transcriptome changes between experimental conditions. Application of RNA-seq relies on accurate transcript annotation for a genome of interest. Here, we find discrepancies in analysis from RNA-seq datasets quantified using Ensembl and RefSeq zebrafish annotations. These issues were due, in part, to variably annotated 3' untranslated regions and thousands of gene models missing from each annotation. Since these discrepancies could compromise downstream analyses and biological reproducibility, we built a more comprehensive zebrafish transcriptome annotation that addresses these deficiencies. Our annotation improves detection of cell type-specific genes in both bulk and single cell RNA-seq datasets, where it also improves resolution of cell clustering. Thus, we demonstrate that our new transcriptome annotation can outperform existing annotations, providing an important resource for zebrafish researchers.

Single-Cell RNA-Seq of the Developing Cardiac Outflow Tract Reveals Convergent Development of the Vascular Smooth Muscle Cells

Cardiac outflow tract (OFT) is a major hotspot for congenital heart diseases. A thorough understanding of the cellular diversity, transitions, and regulatory networks of normal OFT development is essential to decipher the etiology of OFT malformations. We performed single-cell transcriptomic sequencing of 55,611 mouse OFT cells from three developmental stages that generally correspond to the early, middle, and late stages of OFT remodeling and septation. Known cellular transitions, such as endothelial-to-mesenchymal transition, have been recapitulated. In particular, we identified convergent development of the vascular smooth muscle cell (VSMC) lineage where intermediate cell subpopulations were found to be involved in either myocardial-to-VSMC trans-differentiation or mesenchymal-to-VSMC transition. Finally, we uncovered transcriptional regulators potentially governing cellular transitions. Our study provides a single-cell reference map of cell states for normal OFT development and paves the way for further studies of the etiology of OFT malformations at the single-cell level.

Morphological landscape of endothelial cell networks reveals a functional role of glutamate receptors in angiogenesis

Angiogenesis plays a key role in several diseases including cancer, ischemic vascular disease, and Alzheimer's disease. Chemical genetic screening of endothelial tube formation provides a robust approach for identifying signalling components that impact microvascular network morphology as well as endothelial cell biology. However, the analysis of the resulting imaging datasets has been limited to a few phenotypic features such as the total tube length or the number of branching points. Here we developed a high content analysis framework for detailed quantification of various aspects of network morphology including network complexity, symmetry and topology. By applying our approach to a high content screen of 1,280 characterised drugs, we found that drugs that result in a similar phenotype share the same mechanism of action or common downstream signalling pathways. Our multiparametric analysis revealed that a group of glutamate receptor antagonists enhances branching and network connectivity. Using an integrative meta-analysis approach, we validated the link between these receptors and angiogenesis. We further found that the expression of these genes is associated with the prognosis of Alzheimer's patients. In conclusion, our work shows that detailed image analysis of complex endothelial phenotypes can reveal new insights into biological mechanisms modulating the morphogenesis of endothelial networks and identify potential therapeutics for angiogenesis-related diseases.

An Experiment on Ab Initio Discovery of Biological Knowledge from scRNA-Seq Data Using Machine Learning

Expectations of machine learning (ML) are high for discovering new patterns in high-throughput biological data, but most such practices are accustomed to relying on existing knowledge conditions to design experiments. Investigations of the power and limitation of ML in revealing complex patterns from data without the guide of existing knowledge have been lacking. In this study, we conducted systematic experiments on such ab initio knowledge discovery with ML methods on single-cell RNA-sequencing data of early embryonic development. Results showed that a strategy combining unsupervised and supervised ML can reveal major cell lineages with minimum involvement of prior knowledge or manual intervention, and the ab initio mining enabled a new discovery of human early embryonic cell differentiation. The study illustrated the feasibility, significance, and limitation of ab initio ML knowledge discovery on complex biological problems.

Application of single-cell RNA sequencing in embryonic development

Embryonic development is a complex process that is regulated by a series of precise cellular behaviours. The limited number of cells in the early stages of embryonic development represents a challenge for studying early gene regulation and maintaining cell sternness. Single-cell sequencing is a new technology for high-throughput sequencing analysis at the single-cell level that not only reflects the heterogeneity between cells but also reveals gene expression characteristics in different cells from limited samples. Currently, the widespread application of single-cell RNA sequencing technology is gradually changing our understanding of disease pathogenesis. This article reviews the application of single-cell RNA sequencing in embryonic development in recent years and provides innovative ideas for research on embryonic development and the treatment of diseases related to embryonic development.

PRL3-DDX21 Transcriptional Control of Endolysosomal Genes Restricts Melanocyte Stem Cell Differentiation

Melanocytes, replenished throughout life by melanocyte stem cells (MSCs), play a critical role in pigmentation and melanoma. Here, we reveal a function for the metastasis-associated phosphatase of regenerating liver 3 (PRL3) in MSC regeneration. We show that PRL3 binds to the RNA helicase DDX21, thereby restricting productive transcription by RNAPII at master transcription factor (MITF)-regulated endolysosomal vesicle genes. In zebrafish, this mechanism controls premature melanoblast expansion and differentiation from MSCs. In melanoma patients, restricted transcription of this endolysosomal vesicle pathway is a hallmark of PRL3-high melanomas. Our work presents the conceptual advance that PRL3-mediated control of transcriptional elongation is a differentiation checkpoint mechanism for activated MSCs and has clinical relevance for the activity of PRL3 in regenerating tissue and cancer.

Mutations in MYLPF Cause a Novel Segmental Amyoplasia that Manifests as Distal Arthrogryposis

We identified ten persons in six consanguineous families with distal arthrogryposis (DA) who had congenital contractures, scoliosis, and short stature. Exome sequencing revealed that each affected person was homozygous for one of two different rare variants (c.470G>T [p.Cys157Phe] or c.469T>C [p.Cys157Arg]) affecting the same residue of myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF). In a seventh family, a c.487G>A (p.Gly163Ser) variant in MYLPF arose de novo in a father, who transmitted it to his son. In an eighth family comprised of seven individuals with dominantly inherited DA, a c.98C>T (p.Ala33Val) variant segregated in all four persons tested. Variants in MYLPF underlie both dominant and recessively inherited DA. Mylpf protein models suggest that the residues associated with dominant DA interact with myosin whereas the residues altered in families with recessive DA only indirectly impair this interaction. Pathological and histological exam of a foot amputated from an affected child revealed complete absence of skeletal muscle (i.e., segmental amyoplasia). To investigate the mechanism for this finding, we generated an animal model for partial MYLPF impairment by knocking out zebrafish mylpfa. The mylpfa mutant had reduced trunk contractile force and complete pectoral fin paralysis, demonstrating that mylpf impairment most severely affects limb movement. mylpfa mutant muscle weakness was most pronounced in an appendicular muscle and was explained by reduced myosin activity and fiber degeneration. Collectively, our findings demonstrate that partial loss of MYLPF function can lead to congenital contractures, likely as a result of degeneration of skeletal muscle in the distal limb.

Cell competition controls differentiation in mouse embryos and stem cells

Cell competition is a short-range intercellular communication, in which cells compare their fitness with that of their neighbors and eliminate the cells with relatively lower fitness. It is considered important for the formation and maintenance of healthy tissues; however, its exact role during development, especially in mammals, has been obscure. Recent studies in mouse embryonic epiblast and skin tissues revealed that cell differentiation in early embryos and stem cell proliferation tends to produce suboptimal cells, especially during early developmental stages. Cell competition occurs at multiple stages and via multiple mechanisms during development to ensure elimination of such low-quality cells. Thus, quality control via cell competition supports correct development by overcoming the heterogeneity produced during cell differentiation and stem cell proliferation.

Investigating higher-order interactions in single-cell data with scHOT

Single-cell genomics has transformed our ability to examine cell fate choice. Examining cells along a computationally ordered 'pseudotime' offers the potential to unpick subtle changes in variability and covariation among key genes. We describe an approach, scHOT-single-cell higher-order testing-which provides a flexible and statistically robust framework for identifying changes in higher-order interactions among genes. scHOT can be applied for cells along a continuous trajectory or across space and accommodates various higher-order measurements including variability or correlation. We demonstrate the use of scHOT by studying coordinated changes in higher-order interactions during embryonic development of the mouse liver. Additionally, scHOT identifies subtle changes in gene-gene correlations across space using spatially resolved transcriptomics data from the mouse olfactory bulb. scHOT meaningfully adds to first-order differential expression testing and provides a framework for interrogating higher-order interactions using single-cell data.

Cellular and developmental basis of orofacial clefts

During craniofacial development, defective growth and fusion of the upper lip and/or palate can cause orofacial clefts (OFCs), which are among the most common structural birth defects in humans. The developmental basis of OFCs includes morphogenesis of the upper lip, primary palate, secondary palate, and other orofacial structures, each consisting of diverse cell types originating from all three germ layers: the ectoderm, mesoderm, and endoderm. Cranial neural crest cells and orofacial epithelial cells are two major cell types that interact with various cell lineages and play key roles in orofacial development. The cellular basis of OFCs involves defective execution in any one or several of the following processes: neural crest induction, epithelial-mesenchymal transition, migration, proliferation, differentiation, apoptosis, primary cilia formation and its signaling transduction, epithelial seam formation and disappearance, periderm formation and peeling, convergence and extrusion of palatal epithelial seam cells, cell adhesion, cytoskeleton dynamics, and extracellular matrix function. The latest cellular and developmental findings may provide a basis for better understanding of the underlying genetic, epigenetic, environmental, and molecular mechanisms of OFCs.

Epigenetic regulator function through mouse gastrulation

During ontogeny, proliferating cells become restricted in their fate through the combined action of cell-type specific transcription factors and ubiquitous epigenetic machinery, which recognize universally available histone residues or nucleotides but are nonetheless deployed in a highly context-dependent manner^1,2. The molecular functions of these regulators are generally well understood, but assigning direct developmental roles is hampered by complex mutant phenotypes that often emerge following gastrulation^3,4. Recently, single-cell RNA sequencing (scRNA-seq) and analytical approaches have explored this highly conserved process across numerous model organisms^5–8, including mouse^9–18. To elaborate on these strategies, we investigated a panel of ten essential regulators using a combined zygotic perturbation, scRNA-seq platform where many mutant embryos can be assayed simultaneously to recover robust transcriptional and morphological information. Deeper analysis of central Polycomb Repressive Complex (PRC) 1 and 2 members indicate substantial cooperativity, but distinguishes a PRC2-dominant role in restricting the germline that emerges from gross molecular changes within the initial conceptus. We believe our experimental framework will eventually allow for a fully quantitative view of how cellular diversity emerges using an identical genetic template and from a single totipotent cell.

Biological and Medical Importance of Cellular Heterogeneity Deciphered by Single-Cell RNA Sequencing

The present review discusses recent progress in single-cell RNA sequencing (scRNA-seq), which can describe cellular heterogeneity in various organs, bodily fluids, and pathologies (e.g., cancer and Alzheimer's disease). We outline scRNA-seq techniques that are suitable for investigating cellular heterogeneity that is present in cell populations with very high resolution of the transcriptomic landscape. We summarize scRNA-seq findings and applications of this technology to identify cell types, activity, and other features that are important for the function of different bodily organs. We discuss future directions for scRNA-seq techniques that can link gene expression, protein expression, cellular function, and their roles in pathology. We speculate on how the field could develop beyond its present limitations (e.g., performing scRNA-seq in situ and in vivo). Finally, we discuss the integration of machine learning and artificial intelligence with cutting-edge scRNA-seq technology, which could provide a strong basis for designing precision medicine and targeted therapy in the future.

Single-cell analyses identify distinct and intermediate states of zebrafish pancreatic islet development

Pancreatic endocrine islets are vital for glucose homeostasis. However, the islet developmental trajectory and its regulatory network are not well understood. To define the features of these specification and differentiation processes, we isolated individual islet cells from TgBAC(neurod1:EGFP) transgenic zebrafish and analyzed islet developmental dynamics across four different embryonic stages using a single-cell RNA-seq strategy. We identified proliferative endocrine progenitors, which could be further categorized by different cell cycle phases with the G1/S subpopulation displaying a distinct differentiation potential. We identified endocrine precursors, a heterogeneous intermediate-state population consisting of lineage-primed alpha, beta and delta cells that were characterized by the expression of lineage-specific transcription factors and relatively low expression of terminally differentiation markers. The terminally differentiated alpha, beta, and delta cells displayed stage-dependent differentiation states, which were related to their functional maturation. Our data unveiled distinct states, events and molecular features during the islet developmental transition, and provided resources to comprehensively understand the lineage hierarchy of islet development at the single-cell level.

The transcription factor Hypermethylated in Cancer 1 (Hic1) regulates neural crest migration via interaction with Wnt signaling

The transcription factor Hypermethylated in Cancer 1 (HIC1) is associated with both tumorigenesis and the complex human developmental disorder Miller-Dieker Syndrome. While many studies have characterized HIC1 as a tumor suppressor, HIC1 function in development is less understood. Loss-of-function mouse alleles show embryonic lethality accompanied with developmental defects, including craniofacial abnormalities that are reminiscent of human Miller-Dieker Syndrome patients. However, the tissue origin of the defects has not been reported. In this study, we use the power of the Xenopus laevis model system to explore Hic1 function in early development. We show that hic1 mRNA is expressed throughout early Xenopus development and has a spatial distribution within the neural plate border and in migrating neural crest cells in branchial arches. Targeted manipulation of hic1 levels in the dorsal ectoderm that gives rise to neural and neural crest tissues reveals that both overexpression and knockdown of hic1 result in craniofacial defects with malformations of the craniofacial cartilages. Neural crest specification is not affected by altered hic1 levels, but migration of the cranial neural crest is impaired both in vivo and in tissue explants. Mechanistically, we find that Hic1 regulates cadherin expression profiles and canonical Wnt signaling. Taken together, these results identify Hic1 as a novel regulator of the canonical Wnt pathway during neural crest migration.

Single-cell lineage analysis reveals genetic and epigenetic interplay in glioblastoma drug resistance

BACKGROUND

Tumors can evolve and adapt to therapeutic pressure by acquiring genetic and epigenetic alterations that may be transient or stable. A precise understanding of how such events contribute to intratumoral heterogeneity, dynamic subpopulations, and overall tumor fitness will require experimental approaches to prospectively label, track, and characterize resistant or otherwise adaptive populations at the single-cell level. In glioblastoma, poor efficacy of receptor tyrosine kinase (RTK) therapies has been alternatively ascribed to genetic heterogeneity or to epigenetic transitions that circumvent signaling blockade.

RESULTS

We combine cell lineage barcoding and single-cell transcriptomics to trace the emergence of drug resistance in stem-like glioblastoma cells treated with RTK inhibitors. Whereas a broad variety of barcoded lineages adopt a Notch-dependent persister phenotype that sustains them through early drug exposure, rare subclones acquire genetic changes that enable their rapid outgrowth over time. Single-cell analyses reveal that these genetic subclones gain copy number amplifications of the insulin receptor substrate-1 and substrate-2 (IRS1 or IRS2) loci, which activate insulin and AKT signaling programs. Persister-like cells and genomic amplifications of IRS2 and other loci are evident in primary glioblastomas and may underlie the inefficacy of targeted therapies in this disease.

CONCLUSIONS

A method for combined lineage tracing and scRNA-seq reveals the interplay between complementary genetic and epigenetic mechanisms of resistance in a heterogeneous glioblastoma tumor model.

Lineage reconstruction from clonal correlations

Animals begin life as a single cell that divides and differentiates to form a complex body. In doing so, cells make a sequence of fate decisions, often depicted as a tree. A goal in developmental biology is to chart the structure of this tree across tissues, typically by tagging cells and tracking their offspring. Recent advances in DNA sequencing enable tracking thousands of cells simultaneously using unique DNA barcodes, but one can construct false differentiation hierarchies from barcode data. Here, we apply the theory of branching processes to derive conditions under which barcode statistics correctly encode developmental hierarchy. We use this formal basis to develop a practical pipeline for analyzing lineage barcoding experiments. The pipeline is demonstrated in studying hematopoiesis.

A central task in developmental biology is to learn the sequence of fate decisions that leads to each mature cell type in a tissue or organism. Recently, clonal labeling of cells using DNA barcodes has emerged as a powerful approach for identifying cells that share a common ancestry of fate decisions. Here we explore the idea that stochasticity of cell fate choice during tissue development could be harnessed to read out lineage relationships after a single step of clonal barcoding. By considering a generalized multitype branching process, we determine the conditions under which the final distribution of barcodes over observed cell types encodes their bona fide lineage relationships. We then propose a method for inferring the order of fate decisions. Our theory predicts a set of symmetries of barcode covariance that serves as a consistency check for the validity of the method. We show that broken symmetries may be used to detect multiple paths of differentiation to the same cell types. We provide computational tools for general use. When applied to barcoding data in hematopoiesis, these tools reconstruct the classical hematopoietic hierarchy and detect couplings between monocytes and dendritic cells and between erythrocytes and basophils that suggest multiple pathways of differentiation for these lineages.

Single cell genomics and developmental biology: moving beyond the generation of cell type catalogues

Major developmental processes such as gastrulation and early embryogenesis rely on a complex network of cell-cell interactions, chromatin remodeling, and transcriptional regulators. This makes it challenging to study early development when using bulk populations of cells. Recent advances in single-cell technologies have allowed researchers to better understand the interactions between different molecular modalities and the heterogeneities within classically defined cell types. As new single-cell technologies mature, they have the potential of providing a step-change in our understanding of embryogenesis. In this review, we summarize recent advances in single-cell technologies with particular focus on those that lend insight into early organogenesis. We then discuss current pitfalls and implications for future research.

Advancing Stem Cell Research through Multimodal Single-Cell Analysis

Technological advances play a key role in furthering our understanding of stem cell biology, and advancing the prospects of regenerative therapies. Highly parallelized methods, developed in the last decade, can profile DNA, RNA, or proteins in thousands of cells and even capture data across two or more modalities (multiomics). This allows unbiased and precise definition of molecular cell states, thus allowing classification of cell types, tracking of differentiation trajectories, and discovery of underlying mechanisms. Despite being based on destructive techniques, novel experimental and bioinformatic approaches enable embedding and extraction of temporal information, which is essential for deconvolution of complex data and establishing cause and effect relationships. Here, we provide an overview of recent studies pertinent to stem cell biology, followed by an outlook on how further advances in single-cell molecular profiling and computational analysis have the potential to shape the future of both basic and translational research.

Lineage tracing meets single-cell omics: opportunities and challenges

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

Neural circuit reorganisation after spinal cord injury in zebrafish

Spinal cord injuries disrupt signalling from the brain leading to loss of limb, locomotion, sexual and bladder function, usually irreversible in humans. In zebrafish, recovery of function occurs in a few days for larvae or a few weeks for adults due to regrowth of axons and de novo neurogenesis. Together with its genetic amenability and optical clarity, this makes zebrafish a powerful animal model to study circuit reorganisation after spinal cord injuries. With the fast evolution of techniques, we can forecast significative improvements of our knowledge of the mechanisms leading to successful or failed recovery of spinal cord function. We review here the present knowledge on the subject, the new technological approaches and we propose future directions of research.

The maternal coordinate system: Molecular-genetics of embryonic axis formation and patterning in the zebrafish

Axis specification of the zebrafish embryo begins during oogenesis and relies on proper formation of well-defined cytoplasmic domains within the oocyte. Upon fertilization, maternally-regulated cytoplasmic flow and repositioning of dorsal determinants establish the coordinate system that will build the structure and developmental body plan of the embryo. Failure of specific genes that regulate the embryonic coordinate system leads to catastrophic loss of body structures. Here, we review the genetic principles of axis formation and discuss how maternal factors orchestrate axis patterning during zebrafish early embryogenesis. We focus on the molecular identity and functional contribution of genes controlling critical aspects of oogenesis, egg activation, blastula, and gastrula stages. We examine how polarized cytoplasmic domains form in the oocyte, which set off downstream events such as animal-vegetal polarity and germ line development. After gametes interact and form the zygote, cytoplasmic segregation drives the animal-directed reorganization of maternal determinants through calcium- and cell cycle-dependent signals. We also summarize how maternal genes control dorsoventral, anterior-posterior, mesendodermal, and left-right cell fate specification and how signaling pathways pattern these axes and tissues during early development to instruct the three-dimensional body plan. Advances in reverse genetics and phenotyping approaches in the zebrafish model are revealing positional patterning signatures at the single-cell level, thus enhancing our understanding of genotype-phenotype interactions in axis formation. Our emphasis is on the genetic interrogation of novel and specific maternal regulatory mechanisms of axis specification in the zebrafish.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

Unravelling the heterogeneity and dynamic relationships of tumor-infiltrating T cells by single-cell RNA sequencing analysis

T cells are crucial for the success of immune-based cancer therapy. Reinvigorating antitumor T cell activity by blocking checkpoint inhibitory receptors has provided clinical benefits for many cancer patients. However, the efficacy of these treatments varies in cancer patients and the mechanisms underlying these diverse responses remain elusive. The density and status of tumor-infiltrating T cells have been shown to positively correlate with patient response to checkpoint blockades. Therefore, further understanding of the heterogeneity, clonal expansion, migration, and effector functions of tumor-infiltrating T cells will provide fundamental insights into antitumor immune responses. To this end, recent advances in single-cell RNA sequencing technology have enabled profound and extensive characterization of intratumoral immune cells and have improved our understanding of their dynamic relationships. Here, we summarize recent progress in single-cell RNA sequencing technology and current strategies to uncover heterogeneous tumor-infiltrating T cell subsets. In particular, we discuss how the coupling of deep transcriptome information with T cell receptor (TCR)-based lineage tracing has furthered our understanding of intratumoral T cell populations. We also discuss the functional implications of various T cell subsets in tumors and highlight the identification of novel T cell markers with therapeutic or prognostic potential.

Single-Cell Analysis of Human Retina Identifies Evolutionarily Conserved and Species-Specific Mechanisms Controlling Development

The development of single-cell RNA sequencing (scRNA-seq) has allowed high-resolution analysis of cell-type diversity and transcriptional networks controlling cell-fate specification. To identify the transcriptional networks governing human retinal development, we performed scRNA-seq analysis on 16 time points from developing retina as well as four early stages of retinal organoid differentiation. We identified evolutionarily conserved patterns of gene expression during retinal progenitor maturation and specification of all seven major retinal cell types. Furthermore, we identified gene-expression differences between developing macula and periphery and between distinct populations of horizontal cells. We also identified species-specific patterns of gene expression during human and mouse retinal development. Finally, we identified an unexpected role for ATOH7 expression in regulation of photoreceptor specification during late retinogenesis. These results provide a roadmap to future studies of human retinal development and may help guide the design of cell-based therapies for treating retinal dystrophies.

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

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

Tools and Concepts for Interrogating and Defining Cellular Identity

Defining the mechanisms that generate specialized cell types and coordinate their functions is critical for understanding organ development and renewal. New tools and discoveries are challenging and refining our definitions of a cell type. A rapidly growing toolkit for single-cell analyses has expanded the number of markers that can be assigned to a cell simultaneously, revealing heterogeneity within cell types that were previously regarded as homogeneous populations. Additionally, cell types defined by specific molecular markers can exhibit distinct, context-dependent functions; for example, between tissues in homeostasis and those responding to damage. Here we review the current technologies used to identify and characterize cells, and we discuss how experimental and pathological perturbations are adding increasing complexity to our definitions of cell identity.