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

Single-cell RNA sequencing reveals tissue-specific transcriptomic changes induced by perfluorooctanesulfonic acid (PFOS) in larval zebrafish (Danio rerio)

Perfluorooctanesulfonic acid (PFOS) elicits adverse effects on numerous organs and developmental processes but the mechanisms underlying these effects are not well understood. Here, we use single-cell RNA-sequencing to assess tissue-specific transcriptomic changes in zebrafish (Danio rerio) larvae exposed to 16 µM PFOS or dimethylsulfoxide (0.01 %) from 3-72 h post fertilization (hpf). Data analysis was multi-pronged and included pseudo-bulk, untargeted clustering, informed pathway queries, and a cluster curated for hepatocyte biomarkers (fabp10a, and apoa2). Overall, 8.63 % (2390/27698) genes were significantly differentially expressed. Results from untargeted analysis revealed 22 distinct clusters that were manually annotated to specific tissues using a weight-of-evidence approach. The clusters with the highest number of significant differentially expressed genes (DEGs) were digestive organs, muscle, and otolith. Additionally, we assessed the distribution of pathway-specific genes known to be involved in PFOS toxicity: the PPAR pathway, β-oxidation of fatty acids, the Nfe2l2 pathway, and epigenetic modifications by DNA methylation, across clusters and identified the blood-related tissue to be the most sensitive. The curated hepatocyte cluster showed 220 significant DEGs and was enriched for the Notch signaling pathway. These findings provide insights into both established and novel sensitive target tissues and molecular mechanisms of developmental toxicity of PFOS.

A multi-kingdom genetic barcoding system for precise clone isolation

Cell-tagging strategies with DNA barcodes have enabled the analysis of clone size dynamics and clone-restricted transcriptomic landscapes in heterogeneous populations. However, isolating a target clone that displays a specific phenotype from a complex population remains challenging. Here we present a multi-kingdom genetic barcoding system, CloneSelect, which enables a target cell clone to be triggered to express a reporter gene for isolation through barcode-specific CRISPR base editing. In CloneSelect, cells are first stably tagged with DNA barcodes and propagated so that their subpopulation can be subjected to a given experiment. A clone that shows a phenotype or genotype of interest at a given time can then be isolated from the initial or subsequent cell pools stored during the experiment using CRISPR base editing. CloneSelect is scalable and compatible with single-cell RNA sequencing. We demonstrate the versatility of CloneSelect in human embryonic kidney 293T cells, mouse embryonic stem cells, human pluripotent stem cells, yeast cells and bacterial cells.

Parent-of-origin regulation by maternal auts2 shapes neurodevelopment and behavior in fish

BACKGROUND

Parental experience can influence progeny behavior through gamete-mediated non-genetic inheritance, that is, mechanisms that do not involve changes in inherited DNA sequence. However, underlying mechanisms remain poorly understood in vertebrates, especially for maternal effects. Here, we use the medaka, a model fish species, to investigate the role of auts2a, the ortholog of human AUTS2, a gene repressed in the fish oocyte following maternal stress and associated with neurodevelopmental disorders.

RESULTS

We show that auts2a expression in the oocyte influences long-term progeny behavior, including anxiety-like behavior and environment recognition capabilities. Using single-nuclei RNA-sequencing, we reveal that maternal auts2a influences gene expression in neural cell populations during neurodevelopment. We also show that maternal auts2a knock-out triggers differences in maternally inherited factors, including early embryonic transcriptional and post-transcriptional regulators.

CONCLUSIONS

Together, our results reveal the unsuspected role of an autism-related gene expressed in the mother's oocyte in shaping progeny neurodevelopment and behavior. Finally, we report that auts2a/AUTS2 is part of a group of evolutionarily conserved genes associated with human neurodevelopmental disorders and expressed in oocytes across species, from fish to mammals. These findings raise important questions about their potential role in the non-genetic regulation of progeny neurodevelopment and behavior in vertebrates.

Rbm24a dictates mRNA recruitment for germ granule assembly in zebrafish

The germ granules are ribonucleoprotein (RNP) biomolecular condensates that determine the fate of primordial germ cells (PGCs) and serve as a model for studying RNP granule assembly. Here, we show that the maternal RNA-binding protein Rbm24a is a key factor governing the specific sorting of mRNAs into germ granules. Mechanistically, Rbm24a interacts with the germ plasm component Buc to dictate the specific recruitment of germ plasm mRNAs into phase-separated condensates. Germ plasm particles lacking Rbm24a and mRNAs fail to undergo kinesin-dependent transport toward cleavage furrows where small granules fuse into large aggregates. Therefore, the loss of maternal Rbm24a causes a complete degradation of the germ plasm and the disappearance of PGCs. These findings demonstrate that the Rbm24a/Buc complex functions as a nucleating organizer of germ granules, highlighting an emerging mechanism for RNA-binding proteins in reading and recruiting RNA components into a phase-separated protein scaffold.

Ectoderm barcoding reveals neural and cochlear compartmentalization

Placodes and the neural crest are defining features of vertebrates. In this study, we investigate their lineages in mice using in utero approaches. We demonstrated that nanoinjection at embryonic day 7.5 targeted the ectoderm, including the future nervous system, placodes, and neural crest, allowing highly efficient manipulation of the future nervous system and inner ear. By using heritable DNA barcodes and high-throughput next-generation single-cell lineage tracing, we elucidated convergent differentiation pathways and identified distinct nervous system-, neural crest-, and otic placode-derived lineages. Clonal analyses identified early neural and cochlear compartmentalization, linking differentiated cell types to their progenitors or cellular siblings. This provides foundational insights for neuroscience and developmental biology.

Deletion of hepcidin disrupts iron homeostasis and hematopoiesis in zebrafish embryogenesis

Iron is essential for cell growth and hematopoiesis, which is regulated by hepcidin (hamp). However, the role of hamp in zebrafish hematopoiesis remains unclear. Here, we have created a stable hamp knockout zebrafish model using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease 9 system (CRISPR/Cas9 system). Our study revealed that hamp deletion led to maternal iron overload in embryos, significantly downregulating hemoglobin genes and reducing hemoglobin content. Single-cell RNA sequencing identified abnormal expression patterns in blood progenitor cells, with a specific progenitor subtype showing increased ferroptosis and delayed development. By crossing hamp knockout zebrafish with a gata1+ line (blood cells labeled fish line), we confirmed ferroptosis in blood progenitor cells. These findings underscore the crucial role of hamp in iron regulation and hematopoiesis, offering novel insights into developmental biology and potential therapeutic targets for blood disorders.

Regulatory and disruptive variants in the CLCN2 gene are associated with modified skin color pattern phenotypes in the corn snake

BACKGROUND

Snakes exhibit a broad variety of adaptive colors and color patterns, generated by the spatial arrangement of chromatophores, but little is known of the mechanisms responsible for these spectacular traits. Here, we investigate a mono-locus trait with two recessive alleles, motley and stripe, that both cause pattern aberrations in the corn snake.

RESULTS

We use mapping-by-sequencing to identify the genomic interval where the causal mutations reside. With our differential gene expression analyses, we find that CLCN2 (Chloride Voltage-Gated Channel 2), a gene within the genomic interval, is significantly downregulated in Motley embryonic skin. Furthermore, we identify the stripe allele as the insertion of an LTR-retrotransposon in CLCN2, resulting in a disruptive mutation of the protein. We confirm the involvement of CLCN2 in color pattern formation by producing knock-out snakes that present a phenotype similar to Stripe. In humans and mice, disruption of CLCN2 results in leukoencephalopathy, as well as retinal and testes degeneration. Our single-cell transcriptomic analyses in snakes reveal that CLCN2 is indeed expressed in chromatophores during embryogenesis and in the adult brain, but the behavior and fertility of Motley and Stripe corn snakes are not impacted.

CONCLUSIONS

Our genomic, transcriptomic, and functional analyses identify a plasma membrane anion channel to be involved in color pattern development in snakes and show that an active LTR-retrotransposon might be a key driver of trait diversification in corn snakes.

Unsupervised spatiotemporal classification of deformation patterns of embryonic tissues matches their fate map

During morphogenesis, embryonic tissues display fluid-like behavior with fluctuating strain rates. Digital cell lineages reconstructed from 4D images of developing zebrafish embryos are used to infer representative tissue deformation patterns and their association with developmental events. Finite deformation analysis along cell trajectories and unsupervised machine learning are applied to obtain reduced-order models condensing the collective cell motions, delineating tissue domains with distinct 4D biomechanical behavior. This reduced-order kinematic description is reproducible across specimens and matches fate maps of the zebrafish brain in wild-type and nodal pathway mutants (zoeptz57/tz57 ), shedding light into the morphogenetic defects causing these mutants' cyclopia. Furthermore, the inferred kinematic maps also match expression maps of the gene transcription factor goosecoid (gsc). In summary, this work introduces an objective analytical framework to systematically unravel the complex spatiotemporal patterns of embryonic tissue deformations and couple them with cell fate and gene expression maps.

Embryonic exposure of estrogen and BPA in zebrafish leads to ADHD-like and ASD-like phenotypes, respectively

Exposure to the estrogenic pollutant bisphenol A (BPA) during pregnancy and early childhood is a risk factor for numerous neurodevelopmental and psychiatric disorders in humans. To understand why early BPA exposure is associated with attention deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD), we have analyzed a series of zebrafish behaviors, neurodevelopmental process, and gene expression profiles, after a moderate level of estrogen (17β-estradiol, E2, as a positive control) and BPA treatments during embryogenesis (2-48 h post fertilization). E2 exposure-caused hyperactivity was likely due to elevated expression of cyp19a1b since blocking aromatase activity rescued the defect. Furthermore, E2 exposure resulted in impulsive behaviors, perhaps due to a reduced expression of brain th (crucial for dopamine synthesis), resembling the ADHD phenotypes. However, the hyperactivity upon BPA exposure was due to a reduction of GABAergic neurons, particularly in the midbrain. BPA-exposed fish were less-social, with increased repetitive behaviors and escape rate (during strobe light stimulation), like the ASD phenotypes. Taking advantage of published single-cell and bulk RNA-sequencing data related to zebrafish BPA exposure, we uncovered that embryonic midbrain GABAergic neurons express less stmn1a upon BPA exposure. When stmn1a function was partially lost, 14-day post-fertilization larvae became less social, further stressing the ASD phenotype after BPA exposure. Upon embryonic E2 and BPA exposure, we have unexpectedly unveiled zebrafish ADHD-like and ASD-like phenotypes, respectively, suggesting that women of childbearing age should be cautious to use BPA and estrogen related products.

Empirical methods that provide physical descriptions of dynamic cellular processes

We review empirical methods that can be used to provide physical descriptions of dynamic cellular processes during development and disease. Our focus will be nonspatial descriptions and the inference of underlying interaction networks including cell-state lineages, gene regulatory networks, and molecular interactions in living cells. Our overarching questions are: How much can we learn from just observing? To what degree is it possible to infer causal and/or precise mathematical relationships from observations? We restrict ourselves to data sets arising from only observations, or experiments in which minimal perturbations have taken place to facilitate observation of the systems as they naturally occur. We discuss analysis perspectives in order from those offering the least descriptive power but requiring the least assumptions such as statistical associations. We end with those that are most descriptive, but require stricter assumptions and more previous knowledge of the systems such as causal inference and dynamical systems approaches. We hope to provide and encourage the use of a wide array of options for quantitative cell biologists to learn as much as possible from their observations at all stages of understanding of their system of interest. Finally, we provide our own recipe of how to empirically determine quantitative relationships and growth laws from live-cell microscopy data, the resultant predictions of which can then be verified with perturbation experiments. We also include an extended supplement that describes further inference algorithms and theory for the interested reader.

Zebrafish and cellular models of SELENON-Congenital myopathy exhibit novel embryonic and metabolic phenotypes

BACKGROUND

SELENON-Congenital Myopathy (SELENON-CM) is a rare congenital myopathy caused by mutations of the SELENON gene characterized by axial muscle weakness and progressive respiratory insufficiency. Muscle histopathology may be non-specific, but commonly includes multiminicores or a dystrophic pattern. The SELENON gene encodes selenoprotein N (SelN), a selenocysteine-containing redox enzyme located in the endo/sarcoplasmic reticulum membrane where it colocalizes with mitochondria-associated membranes. However, the molecular mechanism(s) by which SelN deficiency cause SELENON-CM remain poorly understood. A hurdle is the lack of cellular and animal models that show easily assayable phenotypes.

METHODS

Using CRISPR-Cas9 we generated three zebrafish models of SELENON-CM, which were then studied by spontaneous coiling, hatching, and activity assays. We also performed selenon coexpression analysis using a single cell RNAseq zebrafish embryo-atlas. SelN-deficient myoblasts were generated and assayed for glutathione, reactive oxygen species, carbonylation, and nytrosylation levels. Finally, we tested Selenon-deficient myoblasts' metabolism using a Seahorse cell respirometer.

RESULTS

We report deep-phenotyping of SelN-deficient zebrafish and muscle cells. SelN-deficient zebrafish exhibit changes in embryonic muscle function and swimming activity in larvae. Analysis of single cell RNAseq data in a zebrafish embryo-atlas revealed coexpression of selenon and genes involved in the glutathione redox pathway. SelN-deficient zebrafish and mouse myoblasts exhibit altered glutathione and redox homeostasis, as well as abnormal patterns of energy metabolism, suggesting roles for SelN in these functions.

CONCLUSIONS

These data demonstrate a role for SelN in zebrafish early development and myoblast metabolism and provide a basis for cellular and animal model assays for SELENON-CM.

NSAID-mediated cyclooxygenase inhibition disrupts ectodermal derivative formation in axolotl embryos

Embryonic exposures to non-steroidal anti-inflammatory drugs (NSAIDs) have been linked to preterm birth, neural tube closure defects, abnormal enteric innervation, and craniofacial malformations, potentially due to disrupted neural tube or neural crest (NC) cell development. Naproxen (NPX), a common non-steroidal anti-inflammatory drug (NSAID) used to relieve pain and inflammation, exerts its effects through non-selective cyclooxygenase (COX) inhibition. Our lab has identified that the cyclooxygenase (COX-1 and COX-2) isoenzymes are expressed during the early stages of vertebrate embryonic development, and that global inhibition of COX-1 and COX-2 function disrupts NC cell migration and differentiation in Ambystoma mexicanum (axolotl) embryos. NC cells differentiate into various adult tissues including craniofacial cartilage, bone, and neurons in the peripheral and enteric nervous systems. To investigate the specific phenotypic and molecular effects of NPX exposure on NC development and differentiation, and to identify molecular links between COX inhibition and NC derivative anomalies, we exposed late neurula and early tailbud stage axolotl embryos to various concentrations of NPX and performed immunohistochemistry (IHC) for markers of migratory and differentiating NC cells. Our results reveal that NPX exposure impairs the migration of SOX9+ NC cells, leading to abnormal development of craniofacial cartilage structures, including Meckel's cartilage in the jaw. NPX exposure also alters the expression of markers associated with peripheral and central nervous system (PNS and CNS) development, suggesting concurrent neurodevelopmental changes.

The Mechanics of Building Functional Organs

Organ morphogenesis is multifaceted, multiscale, and fundamentally a robust process. Despite the complex and dynamic nature of embryonic development, organs are built with reproducible size, shape, and function, allowing them to support organismal growth and life. This striking reproducibility of tissue form exists because morphogenesis is not entirely hardwired. Instead, it is an emergent product of mechanochemical information flow, operating across spatial and temporal scales-from local cellular deformations to organ-scale form and function, and back. In this review, we address the mechanical basis of organ morphogenesis, as understood by observations and experiments in living embryos. To this end, we discuss how mechanical information controls the emergence of a highly conserved set of structural motifs that shape organ architectures across the animal kingdom: folds and loops, tubes and lumens, buds, branches, and networks. Moving forward, we advocate for a holistic conceptual framework for the study of organ morphogenesis, which rests on an interdisciplinary toolkit and brings the embryo center stage.

Decoding the blueprints of embryo development with single-cell and spatial omics

Embryonic development is a complex and intricately regulated process that encompasses precise control over cell differentiation, morphogenesis, and the underlying gene expression changes. Recent years have witnessed a remarkable acceleration in the development of single-cell and spatial omic technologies, enabling high-throughput profiling of transcriptomic and other multi-omic information at the individual cell level. These innovations offer fresh and multifaceted perspectives for investigating the intricate cellular and molecular mechanisms that govern embryonic development. In this review, we provide an in-depth exploration of the latest technical advancements in single-cell and spatial multi-omic methodologies and compile a systematic catalog of their applications in the field of embryonic development. We deconstruct the research strategies employed by recent studies that leverage single-cell sequencing techniques and underscore the unique advantages of spatial transcriptomics. Furthermore, we delve into both the current applications, data analysis algorithms and the untapped potential of these technologies in advancing our understanding of embryonic development. With the continuous evolution of multi-omic technologies, we anticipate their widespread adoption and profound contributions to unraveling the intricate molecular foundations underpinning embryo development in the foreseeable future.

Comparative analysis of neural crest development in the chick and mouse

A core framework of the gene regulatory network (GRN) governing neural crest (NC) cell development has been generated by integrating separate inputs from diverse model organisms rather than direct comparison. This has limited insights into the diversity of genes in the NC cell GRN and extent of conservation of newly identified transcriptional signatures in cell differentiation and invasion. Here, we address this by leveraging the strengths and accessibility of the avian embryo to precise developmental staging by egg incubation and use an integrated analysis of chick (HH13) and mouse (E9.5) embryo tissue samples collected during NC cell migration into pharyngeal arches 1-2 (PA1 and PA2). We successfully identify a cluster of NC cells containing both mouse and chick cells that share expression of Lmo4, Tfap2B, Sox10, and Twist1, and distinct genes that lack known conserved roles in NC. Importantly, we discovered a cluster of cells exhibiting a conserved transcriptional signature associated with the NC cell migratory wavefront in both mouse and chick, including KAZALD1, BAMBI, DES, and GPC3. We confirm their expression is restricted to leader mouse NCs by multiplexed FISH. Together, these data offer novel insights into the transcriptional programs that underlie NC cell migration and establish the foundation for future comparative functional analyses.

Single-cell morphological tracking of cell states to identify small-molecule modulators of liver differentiation

We have developed a single-cell assay that combines Cell Painting-a morphological profiling assay-with trajectory inference analysis. We have applied this morphological trajectory inference to the bi-potent HepaRG liver progenitor cell line allowing us to track liver cell fate and map small-molecule-induced changes using a morphological atlas of liver cell differentiation. Our overarching goal is to demonstrate the potential of Cell Painting to study biological processes as continuous trajectories at the single-cell level, enhancing resolution and biological understanding. This work has identified small-molecule Src family kinase inhibitors that promote the differentiation of HepaRG cells toward a hepatocyte-like lineage as well as primary human hepatic progenitor cells toward a hepatocyte-like phenotype in vitro. These findings could significantly advance research on liver cell regeneration mechanisms and facilitate the development of cell-based and small-molecule therapies.

NSAID-mediated cyclooxygenase inhibition disrupts ectodermal derivative formation in axolotl embryos

Embryonic exposures to non-stseroidal anti-inflammatory drugs (NSAIDs) have been linked to preterm birth, neural tube closure defects, abnormal enteric innervation, and craniofacial malformations, potentially due to disrupted neural tube or neural crest (NC) cell development. Naproxen (NPX), a common non-steroidal anti-inflammatory drug (NSAID) used to relieve pain and inflammation, exerts its effects through non-selective cyclooxygenase (COX) inhibition. Our lab has identified that the cyclooxygenase (COX-1 and COX-2) isoenzymes are expressed during the early stages of vertebrate embryonic development, and that global inhibition of COX-1 and COX-2 function disrupts NC cell migration and differentiation in Ambystoma mexicanum (axolotl) embryos. NC cells differentiate into various adult tissues including craniofacial cartilage, bone, and neurons in the peripheral and enteric nervous systems. To investigate the specific phenotypic and molecular effects of NPX exposure on NC development and differentiation, and to identify molecular links between COX inhibition and NC derivative anomalies, we exposed late neurula and early tailbud stage axolotl embryos to various concentrations of NPX and performed immunohistochemistry (IHC) for markers of migratory and differentiating NC cells. Our results reveal that NPX exposure impairs the migration of SOX9+ NC cells, leading to abnormal development of craniofacial cartilage structures, including Meckel's cartilage in the jaw. NPX exposure also alters the expression of markers associated with peripheral and central nervous system (PNS and CNS) development, suggesting concurrent neurodevelopmental changes.

Single-cell sequencing provides clues about the developmental genetic basis of evolutionary adaptations in syngnathid fishes

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genome assemblies revealed suggestive gene content differences and provided the opportunity for detailed genetic analyses. We created a single-cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined the spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting that derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how the novelties of these fish evolved.

An ancient apical patterning system sets the position of the forebrain in chordates

The evolutionary origin of the vertebrate brain remains a major subject of debate, as its development from a dorsal tubular neuroepithelium is unique to chordates. To shed light on the evolutionary emergence of the vertebrate brain, we compared anterior neuroectoderm development across deuterostome species, using available single-cell datasets from sea urchin, amphioxus, and zebrafish embryos. We identified a conserved gene co-expression module, comparable to the anterior gene regulatory network (aGRN) controlling apical organ development in ambulacrarians, and spatially mapped it by multiplexed in situ hybridization to the developing retina and hypothalamus of chordates. Using functional approaches, we show Wnt signaling regulating this co-expression module in amphioxus, like the aGRN in echinoderms, and that its overactivation suppresses forebrain identity. This suggests a previously undescribed role for Wnt signaling in amphioxus in determining the position of the forebrain. We propose this Wnt-regulated gene co-expression module as a possible mechanism by which the brain set antero-dorsally early in chordate evolution.

Neural plate pre-patterning enables specification of intermediate neural progenitors in the spinal cord

Dorsal-ventral patterning of neural progenitors in the posterior neural tube, which gives rise to the spinal cord, has served as a model system to understand how extracellular signals organize developing tissues. While previous work has shown that signaling gradients diversify progenitor fates at the dorsal and ventral ends of the tissue, the basis of fate specification in intermediate regions has remained unclear. Here we use zebrafish to investigate the neural plate, which precedes neural tube formation, and show that its pre-patterning by a distinct signaling environment enables intermediate fate specification. Systematic spatial analysis of transcription factor (TF) expression and signaling activity using a reference-based mapping approach shows that the neural plate is partitioned into a striking complexity of TF co-expression states that, in part, correspond to the activity of gastrulation signals such as FGF and Wnt that persist through axis extension. Using in toto analysis of cellular movement combined with fate mapping, we find that dbx1b-expressing intermediate progenitors (p0) originate from a neural-plate specific state characterized by transient co-expression of the TFs pax3a, olig4 and her3. Finally, we show that this state is defined by Wnt signaling in the posterior neural plate and that ectopic Wnt activation within pax3a/olig4+ cells is sufficient to promote dbx1b expression. Our data broadly support a model in which neural progenitor specification occurs through the sequential use of multiple signals to progressively diversify the neural tissue as it develops. This has implications for in vitro differentiation of spinal cord cell types and for understanding signal-based patterning in other developmental contexts.

Benchmarking cross-species single-cell RNA-seq data integration methods: towards a cell type tree of life

Cross-species single-cell RNA-seq data hold immense potential for unraveling cell type evolution and transferring knowledge between well-explored and less-studied species. However, challenges arise from interspecific genetic variation, batch effects stemming from experimental discrepancies and inherent individual biological differences. Here, we benchmarked nine data-integration methods across 20 species, encompassing 4.7 million cells, spanning eight phyla and the entire animal taxonomic hierarchy. Our evaluation reveals notable differences between the methods in removing batch effects and preserving biological variance across taxonomic distances. Methods that effectively leverage gene sequence information capture underlying biological variances, while generative model-based approaches excel in batch effect removal. SATURN demonstrates robust performance across diverse taxonomic levels, from cross-genus to cross-phylum, emphasizing its versatility. SAMap excels in integrating species beyond the cross-family level, especially for atlas-level cross-species integration, while scGen shines within or below the cross-class hierarchy. As a result, our analysis offers recommendations and guidelines for selecting suitable integration methods, enhancing cross-species single-cell RNA-seq analyses and advancing algorithm development.

Drivers of vessel progenitor fate define intermediate mesoderm dimensions by inhibiting kidney progenitor specification

Proper organ formation depends on the precise delineation of organ territories containing defined numbers of progenitor cells. Kidney progenitors reside in bilateral stripes of posterior mesoderm that are referred to as the intermediate mesoderm (IM). Previously, we showed that the transcription factors Hand2 and Osr1 act to strike a balance between the specification of the kidney progenitors in the IM and the vessel progenitors in the laterally adjacent territory. Recently, the transcription factor Npas4l - an early and essential driver of vessel and blood progenitor formation - was shown to inhibit kidney development. Here we demonstrate how kidney progenitor specification is coordinated by hand2, osr1, and npas4l. We find that npas4l and the IM marker pax2a are transiently co-expressed in the posterior lateral mesoderm, and npas4l is necessary to inhibit IM formation. Consistent with the expression of npas4l flanking the medial and lateral sides of the IM, our findings suggest roles for npas4l in defining the IM boundaries at each of these borders. At the lateral IM border, hand2 promotes and osr1 inhibits the formation of npas4l-expressing lateral vessel progenitors, and hand2 requires npas4l to inhibit IM formation and to promote vessel formation. Meanwhile, npas4l appears to have an additional role in suppressing IM fate at the medial border: npas4l loss-of-function enhances hand2 mutant IM defects and results in excess IM generated outside of the lateral hand2-expressing territory. Together, our findings reveal that establishment of the medial and lateral boundaries of the IM requires inhibition of kidney progenitor specification by the neighboring drivers of vessel progenitor fate.

Recording Lineage History with Cellular Barcodes in the Mammary Epithelium and in Breast Cancer

Lineage tracing methods have extensively advanced our understanding of physiological cell behaviour in vivo and in situ and have vastly contributed to decipher the phylogeny and cellular hierarchies during normal and tumour development. In recent years, increasingly complex systems have been developed to track thousands of cells within a given tissue or even entire organisms. Cellular barcoding comprises all techniques designed to genetically label single cells with unique DNA sequences or with a combination of fluorescent proteins, in order to trace their history and lineage production in space and time. We distinguish these two types of cellular barcoding as genetic or optical barcodes. Furthermore, transcribed cellular barcodes can integrate the lineage information with single-cell profiling of each barcoded cell. This enables the potential identification of specific markers or signalling pathways defining distinct stem cell states during development, but also signals promoting tumour growth and metastasis or conferring therapy resistance.In this chapter, we describe recent advances in cellular barcoding technologies and outline experimental and computational challenges. We discuss the biological questions that can be addressed using single-cell dynamic lineage tracing, with a focus on the study of cellular hierarchies in the mammary epithelium and in breast cancer.

Reconstructing Progenitor State Hierarchy and Dynamics Using Lineage Barcoding Data

Measurements of cell phylogeny based on natural or induced mutations, known as lineage barcodes, in conjunction with molecular phenotype have become increasingly feasible for a large number of single cells. In this chapter, we delve into Quantitative Fate Mapping (QFM) and its computational pipeline, which enables the interrogation of the dynamics of progenitor cells and their fate restriction during development. The methods described here include inferring cell phylogeny with the Phylotime model, and reconstructing progenitor state hierarchy, commitment time, population size, and commitment bias with the ICE-FASE algorithm. Evaluation of adequate sampling based on progenitor state coverage statistics is emphasized for interpreting the QFM results. Overall, this chapter describes a general framework for characterizing the dynamics of cell fate changes using lineage barcoding data.

Fate specification triggers a positive feedback loop of TEAD-YAP and NANOG to promote epiblast formation in preimplantation embryos

In preimplantation embryos, epiblast (EPI) fate specification from the inner cell mass is controlled by the segregation of NANOG and GATA6 expression. TEAD-YAP interaction is activated during EPI formation and is required for pluripotency factor expression. These events occur asynchronously with similar timing during EPI formation, and their relationship remains elusive. Here, we examined the relationship between NANOG-GATA6 and TEAD-YAP. The nuclear accumulation of YAP takes place only in EPI-specified cells, and a positive feedback loop operates between NANOG and TEAD-YAP. The effects of TEAD-YAP on SOX2 upregulation in EPI-specified cells are likely indirect. EPI fate specification also alters the response of Nanog, Sox2 and Cdx2 to TEAD-YAP. These results suggest that EPI-fate specification alters the transcriptional network from the morula-like to the EPI-specified state and activates TEAD-YAP to trigger a positive feedback loop with NANOG, which stabilizes the EPI fate. The coordinated occurrence of these processes in individual cells likely supports proper EPI formation under the condition of asynchronous EPI-fate specification.

Activation-derepression synergy enables a bHLH network to coordinate a signal-specific fate response

Stem cells integrate multiple environmental signals to activate appropriate fate programs. To ensure coherent responses, alternative fates must be concomitantly inactivated. However, mechanisms that coordinate fates in a signal-specific manner are not fully understood. Here, we investigate the role of a network of basic-helix-loop-helix (bHLH) transcription factors in neural stem cells, which integrate leukemia inhibitory factor (LIF) and bone morphogenetic protein (BMP) signaling to synergistically induce glial fibrillary acidic protein (GFAP), a key astrocyte-fate determinant. Using quantitative RNA-fluorescence in situ hybridization (FISH) and ectopic expression, we find that multiple bHLHs that promote alternative fates also repress GFAP but are all suppressed by BMP and, to a lesser extent, LIF. Mathematical modeling shows that synergy arises from this coordinated derepression of GFAP combined with its activation by LIF signaling. Finally, we determine how coordinated and tunable derepression results from extensive cross-regulation among bHLHs. Activation-derepression synergy could be broadly utilized to couple signaling and fate, particularly across the numerous developmental systems regulated by bHLH factors.

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

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

tet2 and tet3 regulate cell fate specification and differentiation events during retinal development

Tet enzymes are epigenetic modifiers that impact gene expression via 5mC to 5hmC oxidation. Previous work demonstrated the requirement for Tet and 5hmC during zebrafish retinogenesis. tet2 -/- ;tet3 -/- mutants possessed defects in the formation of differentiated retinal neurons, but the mechanisms underlying these defects are unknown. Here, we leveraged scRNAseq technologies to better understand cell type-specific deficits and molecular signatures underlying the tet2 -/- ;tet3 -/- retinal phenotype. Our results identified defects in the tet2 -/- ;tet3 -/- retinae that included delayed specification of several retinal cell types, reduced maturity across late-stage cones, expansions of immature subpopulations of horizontal and bipolar cells, and altered biases of bipolar cell subtype fates at late differentiation stages. Together, these data highlight the critical role that tet2 and tet3 play as regulators of cell fate specification and terminal differentiation events during retinal development.

Maximizing single cell dissociation protocol for individual zebrafish embryo

Single-cell sequencing has revolutionized our understanding of cellular heterogeneity and cell state, enabling investigations across diverse fields such as developmental biology, immunology, and cancer biology. However, obtaining a high-quality single-cell suspension is still challenging, particularly when starting with limited materials like Zebrafish embryos, a powerful animal model for studying developmental processes and human diseases. Here, we present an optimized protocol for isolating single cells from individual zebrafish embryos, offering a valuable resource for researchers interested in working with limited starting material. The protocol facilitates unique investigations utilizing individual embryos, such as inter-individual genetic differences and embryo-specific lineage tracing analysis. Using a refined single-cell isolation protocol alongside zebrafish as a model organism, researchers can access a resource for exploring the emergence of all types and states of cells, advancing our understanding of cellular processes and disease mechanisms.

Applying single-cell and single-nucleus genomics to studies of cellular heterogeneity and cell fate transitions in the nervous system

Single-cell and single-nucleus genomic approaches can provide unbiased and multimodal insights. Here, we discuss what constitutes a molecular cell atlas and how to leverage single-cell omics data to generate hypotheses and gain insights into cell transitions in development and disease of the nervous system. We share points of reflection on what to consider during study design and implementation as well as limitations and pitfalls.

Untangling Zebrafish Genetic Annotation: Addressing Complexities and Nomenclature Issues in Orthologous Evaluation of TCOF1 and NOLC1

Treacher Collins syndrome (TCS) is a genetic disorder affecting facial development, primarily caused by mutations in the TCOF1 gene. TCOF1, along with NOLC1, play important roles in ribosomal RNA transcription and processing. Previously, a zebrafish model of TCS successfully recapitulated the main characteristics of the syndrome by knocking down the expression of a gene on chromosome 13 (coding for Uniprot ID B8JIY2), which was identified as the TCOF1 orthologue. However, database updates renamed this gene as nolc1 and the zebrafish database (ZFIN) identified a different gene on chromosome 14 as the TCOF1 orthologue (coding for Uniprot ID E7F9D9). NOLC1 and TCOF1 are large proteins with unstructured regions and repetitive sequences that complicate alignments and comparisons. Also, the additional whole genome duplication of teleosts sets further difficulty. In this study, we present evidence that endorses that NOLC1 and TCOF1 are paralogs, and that the zebrafish gene on chromosome 14 is a low-complexity LisH domain-containing factor that displays homology to NOLC1 but lacks essential sequence features to accomplish TCOF1 nucleolar functions. Our analysis also supports the idea that zebrafish, as has been suggested for other non-tetrapod vertebrates, lack the TCOF1 gene that is associated with tripartite nucleolus. Using BLAST searches in a group of teleost genomes, we identified fish-specific sequences similar to E7F9D9 zebrafish protein. We propose naming them "LisH-containing Low Complexity Proteins" (LLCP). Interestingly, the gene on chromosome 13 (nolc1) displays the sequence features, developmental expression patterns, and phenotypic impact of depletion that are characteristic of TCOF1 functions. These findings suggest that in teleost fish, the nucleolar functions described for both NOLC1 and TCOF1 mediated by their repeated motifs, are carried out by a single gene, nolc1. Our study, which is mainly based on computational tools available as free web-based algorithms, could help to solve similar conflicts regarding gene orthology in zebrafish.

A single-nucleus transcriptomic atlas characterizes cell types and their molecular features in the ovary of adult Nile tilapia

In fish species, there is limited analysis of signature transcriptome profiles at the single-cell level in gonadal cells. Here, the molecular signatures of distinct ovarian cell categories in adult Nile tilapia (Oreochromis niloticus) were analysed using single-nucleus RNA sequencing (snRNA-seq). We identified four cell types (oogonia, oocytes, granulosa cell, and thecal cell) based on their specifically expressed genes and biological functions. Similarly, we found some key pathways involved in ovarian development that may affect germline-somatic interactions. A cell-to-cell communication network between the distinct cell types was constructed. We found that the bidirectional communication is mandatory for the development of germ cells and somatic cells in fish ovaries, and the granulosa cells and thecal cells play a central regulating role in the cell network in fish ovary. Additionally, we identified some novel candidate marker genes for various types of ovarian cells and also validated them using in situ hybridization. Our work reveals an ovarian atlas at the cellular and molecular levels and contributes to providing insights into oogenesis and gonad development in fish.

Ionizing radiation inhibits zebrafish embryo hatching through induction of tissue inhibitors of metalloproteinases (TIMPs) expression

Ionizing radiation (IR) has garnered growing attention because of its biological effects on aquatic organisms and humans. Here, we identify the most impacted organs and uncover the molecular mechanisms causing the changes in the context of vertebrate development using single-cell RNA sequencing. Alterations in cellular composition and biological functions were explored using transcriptomic profiling of zebrafish embryos exposed to 5 Gy. Single-cell RNA sequencing analyses unveiled notable shifts in the proportions of brain/central nervous system and hatching gland clusters. Although IR exposure led to increased expression of hatching enzymes, a significant but mild delay in hatching was observed following 5 Gy IR exposure. Gene Ontology analysis showed an increased expression of tissue inhibitors of metalloproteinases (TIMPs), known as matrix metalloproteinase inhibitors, which was confirmed via whole-mount in situ hybridization. Correlation analysis linked TIMPs to transcription factors cebpb and cebpd, which were significantly correlated post-IR exposure. Although no morphological changes were observed in some organs, including the brain, the study reveals substantial alterations in developing vertebrates. Notably, despite increased hatching enzymes, elevated TIMPs in the hatching gland suggest a regulatory mechanism impacting hatching activity. This research contributes to comprehending the ecological repercussions of IR exposure, emphasizing the importance of safety measures for aquatic ecosystems and overall environmental health.

Cell state transitions are decoupled from cell division during early embryo development

As tissues develop, cells divide and differentiate concurrently. Conflicting evidence shows that cell division is either dispensable or required for formation of cell types. Here, to determine the role of cell division in differentiation, we arrested the cell cycle in zebrafish embryos using two independent approaches and profiled them at single-cell resolution. We show that cell division is dispensable for differentiation of all embryonic tissues from early gastrulation to the end of segmentation. However, arresting cell division does slow down differentiation in some cell types, and it induces global stress responses. While differentiation is robust to blocking cell division, the proportions of cells across cell states are not, but show evidence of partial compensation. This work clarifies our understanding of the role of cell division in development and showcases the utility of combining embryo-wide perturbations with single-cell RNA sequencing to uncover the role of common biological processes across multiple tissues.

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

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

The molecular basis for functional divergence of duplicated SOX factors controlling endoderm formation and left-right patterning in zebrafish

Endoderm, one of three primary germ layers of vertebrate embryos, makes major contributions to the respiratory and gastrointestinal tracts and associated organs, including liver and pancreas. In mammals, the transcription factor SOX17 is vital for endoderm organ formation and can induce endoderm progenitor identity. Duplication of ancestral sox17 in the teleost lineage produced the paralogues sox32 and sox17 in zebrafish. Sox32 is required for specification of endoderm and progenitors of the left-right organiser (Kupffer's Vesicle, KV), with Sox17 a downstream target of Sox32 that is implicated in further KV development. Phenotypic evidence therefore suggests functional similarities between zebrafish Sox32 and Sox17 and mammalian SOX17. Here, we directly compare these orthologues and paralogues, using the early zebrafish embryo as a biological platform for functional testing. Our results indicate that, unlike Sox32, human SOX17 cannot induce endoderm specification in zebrafish. Furthermore, using hybrid protein functional analyses, we show that Sox32 specificity for the endoderm gene regulatory network is linked to evolutionary divergence in its DNA-binding HMG domain from its paralogue Sox17. Additionally, changes in the C-terminal regions of Sox32 and Sox17 underpin their differing target specificities. Finally, we establish that specific conserved peptides in the C-terminal domain are essential for the role of Sox17 in establishing correct organ asymmetry. Overall, our results illuminate the molecular basis for functional divergence of Sox32 and Sox17 in vertebrate endoderm development and left-right patterning, and the evolution of SoxF transcription factor function.

Critical role of Spatio-Temporally Regulated Maternal RNAs in Zebrafish Embryogenesis

The maternal-to-zygotic transition shifts regulatory control from maternal to zygotic messenger RNAs (mRNA) through maternal mRNA degradation. While temporal aspects of maternal mRNA decay are known, spatial mechanisms remain underexplored. Using CRISPR-Cas9 and CRISPR-Cas13d systems, we functionally dissected the contribution of maternal versus zygotic fractions and overcame challenges of studying embryonic lethal genes. We identified differentially distributed maternal mRNAs in specific cells and evidenced the critical role of five maternal mRNAs, cth1, arl4d, abi1b, foxa and lhx1a, in embryogenesis. Further, we focused on the functionally uncharacterized cth1 gene, revealing its essential role in gametogenesis and embryogenesis. Cth1 acts as a spatio-temporal RNA decay factor regulating mRNA stability and accumulation of its targets in a spatio-temporal manner through 3'UTR recognition during early development. Furthermore, Cth1 3'UTR drives its spatio-temporal RNA localization. Our findings provide new insights into spatio-temporal RNA decay mechanisms and highlight dual CRISPR-Cas strategies in studying embryonic development.

Comparative Analysis of Neural Crest Development in the Chick and Mouse

A core framework of the gene regulatory network (GRN) governing neural crest (NC) cell development has been generated by integrating separate inputs from diverse model organisms rather than direct comparison. This has limited insights into the diversity of genes in the NC cell GRN and extent of conservation of newly identified transcriptional signatures in cell differentiation and invasion. Here, we address this by leveraging the strengths and accessibility of the avian embryo to precise developmental staging by egg incubation and use an integrated analysis of chick (HH13) and mouse (E9.5) embryo tissue samples collected during NC cell migration into pharyngeal arches 1-2 (PA1 and PA2). We successfully identify a cluster of NC cells containing both mouse and chick cells that share expression of Lmo4 , Tfap2B , Sox10 , and Twist1 , and distinct genes that lack known conserved roles in NC. Importantly, we discovered a cluster of cells exhibiting a conserved transcriptional signature associated with the NC cell migratory wavefront in both mouse and chick, including KAZALD1, BAMBI, DES, and GPC3. We confirm their expression is restricted to leader mouse NCs by multiplexed FISH. Together, these data offer novel insights into the transcriptional programs that underlie NC cell migration and establish the foundation for future comparative functional analyses.

Single-cell sequencing: Current applications in various tuberculosis specimen types

Tuberculosis (TB) is a chronic disease caused by Mycobacterium tuberculosis (M.tb) and responsible for millions of deaths worldwide each year. It has a complex pathogenesis that primarily affects the lungs but can also impact systemic organs. In recent years, single-cell sequencing technology has been utilized to characterize the composition and proportion of immune cell subpopulations associated with the pathogenesis of TB disease since it has a high resolution that surpasses conventional techniques. This paper reviews the current use of single-cell sequencing technologies in TB research and their application in analysing specimens from various sources of TB, primarily peripheral blood and lung specimens. The focus is on how these technologies can reveal dynamic changes in immune cell subpopulations, genes and proteins during disease progression after M.tb infection. Based on the current findings, single-cell sequencing has significant potential clinical value in the field of TB research. Next, we will focus on the real-world applications of the potential targets identified through single-cell sequencing for diagnostics, therapeutics and the development of effective vaccines.

The dynamics and functional impact of tRNA repertoires during early embryogenesis in zebrafish

Embryogenesis entails dramatic shifts in mRNA translation and turnover that reprogram gene expression during cellular proliferation and differentiation. Codon identity modulates mRNA stability during early vertebrate embryogenesis, but how the composition of tRNA pools is matched to translational demand is unknown. By quantitative profiling of tRNA repertoires in zebrafish embryos during the maternal-to-zygotic transition, we show that zygotic tRNA repertoires are established after the onset of gastrulation, succeeding the major wave of zygotic mRNA transcription. Maternal and zygotic tRNA pools are distinct, but their reprogramming does not result in a better match to the codon content of the zygotic transcriptome. Instead, we find that an increase in global translation at gastrulation sensitizes decoding rates to tRNA supply, thus destabilizing maternal mRNAs enriched in slowly translated codons. Translational activation and zygotic tRNA expression temporally coincide with an increase of TORC1 activity at gastrulation, which phosphorylates and inactivates the RNA polymerase III repressor Maf1a/b. Our data indicate that a switch in global translation, rather than tRNA reprogramming, determines the onset of codon-dependent maternal mRNA decay during zebrafish embryogenesis.

Detecting significant expression patterns in single-cell and spatial transcriptomics with a flexible computational approach

Gene expression data holds the potential to shed light on multiple biological processes at once. However, data analysis methods for single cell sequencing mostly focus on finding cell clusters or the principal progression line of the data. Data analysis for spatial transcriptomics mostly addresses clustering and finding spatially variable genes. Existing data analysis methods are effective in finding the main data features, but they might miss less pronounced, albeit significant, processes, possibly involving a subset of the samples. In this work we present SPIRAL: Significant Process InfeRence ALgorithm. SPIRAL is based on Gaussian statistics to detect all statistically significant biological processes in single cell, bulk and spatial transcriptomics data. The algorithm outputs a list of structures, each defined by a set of genes working simultaneously in a specific population of cells. SPIRAL is unique in its flexibility: the structures are constructed by selecting subsets of genes and cells based on statistically significant and consistent differential expression. Every gene and every cell may be part of one structure, more or none. SPIRAL also provides several visual representations of structures and pathway enrichment information. We validated the statistical soundness of SPIRAL on synthetic datasets and applied it to single cell, spatial and bulk RNA-sequencing datasets. SPIRAL is available at https://spiral.technion.ac.il/ .

Mapping lineage-traced cells across time points with moslin

Simultaneous profiling of single-cell gene expression and lineage history holds enormous potential for studying cellular decision-making. Recent computational approaches combine both modalities into cellular trajectories; however, they cannot make use of all available lineage information in destructive time-series experiments. Here, we present moslin, a Gromov-Wasserstein-based model to couple cellular profiles across time points based on lineage and gene expression information. We validate our approach in simulations and demonstrate on Caenorhabditis elegans embryonic development how moslin predicts fate probabilities and putative decision driver genes. Finally, we use moslin to delineate lineage relationships among transiently activated fibroblast states during zebrafish heart regeneration.

Single Cell Sequencing Provides Clues about the Developmental Genetic Basis of Evolutionary Adaptations in Syngnathid Fishes

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genome assemblies revealed suggestive gene content differences and provide the opportunity for detailed genetic analyses. We created a single cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how their novelties evolved.

The pericardium forms as a distinct structure during heart formation

The heart integrates diverse cell lineages into a functional unit, including the pericardium, a mesothelial sac that supports heart movement, homeostasis, and immune responses. However, despite its critical roles, the developmental origins of the pericardium remain uncertain due to disparate models. Here, using live imaging, lineage tracking, and single-cell transcriptomics in zebrafish, we find the pericardium forms within the lateral plate mesoderm from dedicated anterior mesothelial progenitors and distinct from the classic heart field. Imaging of transgenic reporters in zebrafish documents lateral plate mesoderm cells that emerge lateral of the classic heart field and among a continuous mesothelial progenitor field. Single-cell transcriptomics and trajectories of hand2-expressing lateral plate mesoderm reveal distinct populations of mesothelial and cardiac precursors, including pericardial precursors that are distinct from the cardiomyocyte lineage. The mesothelial gene expression signature is conserved in mammals and carries over to post-natal development. Light sheet-based live-imaging and machine learning-supported cell tracking documents that during heart tube formation, pericardial precursors that reside at the anterior edge of the heart field migrate anteriorly and medially before fusing, enclosing the embryonic heart to form a single pericardial cavity. Pericardium formation proceeds even upon genetic disruption of heart tube formation, uncoupling the two structures. Canonical Wnt/β-catenin signaling modulates pericardial cell number, resulting in a stretched pericardial epithelium with reduced cell number upon canonical Wnt inhibition. We connect the pathological expression of secreted Wnt antagonists of the SFRP family found in pediatric dilated cardiomyopathy to increased pericardial stiffness: sFRP1 in the presence of increased catecholamines causes cardiomyocyte stiffness in neonatal rats as measured by atomic force microscopy. Altogether, our data integrate pericardium formation as an independent process into heart morphogenesis and connect disrupted pericardial tissue properties such as pericardial stiffness to pediatric cardiomyopathies.

Robust Double Emulsions for Multicolor Fluorescence-Activated Cell Sorting

Cell-cell interactions are essential for the proper functioning of multicellular organisms. For example, T cells interact with antigen-presenting cells (APCs) through specific T-cell receptor (TCR)-antigen interactions during an immune response. Fluorescence-activated droplet sorting (FADS) is a high-throughput technique for efficiently screening cellular interaction events. Unfortunately, current droplet sorting instruments have significant limitations, most notably related to analytical throughput and complex operation. In contrast, commercial fluorescence-activated cell sorters offer superior speed, sensitivity, and multiplexing capabilities, although their use as droplet sorters is poorly defined and underutilized. Herein, we present a universally applicable and simple-to-implement workflow for generating double emulsions and performing multicolor cell sorting using a commercial FACS instrument. This workflow achieves a double emulsion detection rate exceeding 90%, enabling multicellular encapsulation and high-throughput immune cell activation sorting for the first time. We anticipate that the presented droplet sorting strategy will benefit cell biology laboratories by providing access to an advanced microfluidic toolbox with minimal effort and cost investment.

Inferring cell differentiation maps from lineage tracing data

During development, mulitpotent cells differentiate through a hierarchy of increasingly restricted progenitor cell types until they realize specialized cell types. A cell differentiation map describes this hierarchy, and inferring these maps is an active area of research spanning traditional single marker lineage studies to data-driven trajectory inference methods on single-cell RNA-seq data. Recent high-throughput lineage tracing technologies profile lineages and cell types at scale, but current methods to infer cell differentiation maps from these data rely on simple models with restrictive assumptions about the developmental process. We introduce a mathematical framework for cell differentiation maps based on the concept of potency, and develop an algorithm, Carta, that infers an optimal cell differentiation map from single-cell lineage tracing data. The key insight in Carta is to balance the trade-off between the complexity of the cell differentiation map and the number of unobserved cell type transitions on the lineage tree. We show that Carta more accurately infers cell differentiation maps on both simulated and real data compared to existing methods. In models of mammalian trunk development and mouse hematopoiesis, Carta identifies important features of development that are not revealed by other methods including convergent differentiation of specialized cell types, progenitor differentiation dynamics, and the refinement of routes of differentiation via new intermediate progenitors.

Single‑cell RNA sequencing analysis of human embryos from the late Carnegie to fetal development

BACKGROUND

The cell development atlas of transition stage from late Carnegie to fetal development (7-9 weeks) remain unclear. It can be seen that the early period of human embryos (7-9 weeks) is a critical research gap. Therefore, we employed single‑cell RNA sequencing to identify cell types and elucidate differentiation relationships.

RESULTS

The single‑cell RNA sequencing analysis determines eighteen cell clusters in human embryos during the 7-9 weeks period. We uncover two distinct pathways of cellular development and differentiation. Initially, mesenchymal progenitor cells differentiated into osteoblast progenitor cells and neural stem cells, respectively. Neural stem cells further differentiated into neurons. Alternatively, multipotential stem cells differentiated into adipocyte, hematopoietic stem cells and neutrophil, respectively. Additionally, COL1A2-(ITGA1 + ITGB1) mediated the cell communication between mesenchymal progenitor cells and osteoblast progenitor cells. NCAM1-FGFR1 facilitated the cell communication between mesenchymal progenitor cells and neural stem cells. Notably, NCAM1-NCAM1 as a major contributor mediated the cell communication between neural stem cells and neurons. Moreover, CGA-FSHR simultaneously mediated the communication between multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil. Distinct cell clusters activated specific transcription factors such as HIC1, LMX1B, TWIST1, and et al., which were responsible for their specific functions. These coregulators, such as HOXB13, VSX2, PAX5, and et al., may mediate cell development and differentiation in human embryos.

CONCLUSIONS

We provide the cell development atlas for human embryos (7-9 weeks). Two distinct cell development and differentiation pathways are revealed.

Achieving single-cell-resolution lineage tracing in zebrafish by continuous barcoding mutations during embryogenesis

Unraveling the lineage relationships of all descendants from a zygote is fundamental to advancing our understanding of developmental and stem cell biology. However, existing cell barcoding technologies in zebrafish lack the resolution to capture the majority of cell divisions during embryogenesis. A recently developed method, a substitution mutation-aided lineage-tracing system (SMALT), successfully reconstructed high-resolution cell phylogenetic trees for Drosophila melanogaster. Here, we implement the SMALT system in zebrafish, recording a median of 14 substitution mutations on a one-kilobase-pair barcoding sequence for one-day post-fertilization embryos. Leveraging this system, we reconstruct four cell lineage trees for zebrafish fin cells, encompassing both original and regenerated fin. Each tree consists of hundreds of internal nodes with a median bootstrap support of 99%. Analysis of the obtained cell lineage trees reveals that regenerated fin cells mainly originate from cells in the same part of the fins. Through multiple times sampling germ cells from the same individual, we show the stability of the germ cell pool and the early separation of germ cell and somatic cell progenitors. Our system offers the potential for reconstructing high-quality cell phylogenies across diverse tissues, providing valuable insights into development and disease in zebrafish.

Establishment and transcriptome analysis of single blastomere-derived cell lines from zebrafish

Maintaining chromosome euploidy in zebrafish embryonic cells is challenging because of the degradation of genomic integrity during cell passaging. In this study, we report the derivation of zebrafish cell lines from single blastomeres. These cell lines have a stable chromosome status attributed to BMP4 and exhibit continuous proliferation in vitro. Twenty zebrafish cell lines are successfully established from single blastomeres. Single-cell transcriptome sequencing analysis confirms the fidelity of gene expression profiles throughout long-term culturing of at least 45 passages. The long-term cultured cells are specialized into epithelial cells, exhibiting similar expression patterns validated by integrative transcriptomic analysis. Overall, this work provides a protocol for establishing zebrafish cell lines from single blastomeres, which can serve as valuable tools for in vitro investigations of epithelial cell dynamics in terms of life-death balance and cell fate determination during normal homeostasis.