A single-cell transcriptome atlas for zebrafish development

Probing PAC1 receptor activation across species with an engineered sensor

Class-B1 G-protein-coupled receptors (GPCRs) are an important family of clinically relevant drug targets that remain difficult to investigate via high-throughput screening and in animal models. Here, we engineered PAClight1P78A, a novel genetically encoded sensor based on a class-B1 GPCR (the human PAC1 receptor, hmPAC1R) endowed with high dynamic range (ΔF/F0 = 1100%), excellent ligand selectivity, and rapid activation kinetics (τON = 1.15 s). To showcase the utility of this tool for in vitro applications, we thoroughly characterized and compared its expression, brightness and performance between PAClight1P78A-transfected and stably expressing cells. Demonstrating its use in animal models, we show robust expression and fluorescence responses upon exogenous ligand application ex vivo and in vivo in mice, as well as in living zebrafish larvae. Thus, the new GPCR-based sensor can be used for a wide range of applications across the life sciences empowering both basic research and drug development efforts.

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.

Keratinocytes drive the epithelial hyperplasia key to sea lice resistance in coho salmon

BACKGROUND

Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown.

RESULTS

We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite.

CONCLUSIONS

Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

Gap-junction-mediated bioelectric signaling required for slow muscle development and function in zebrafish

Bioelectric signaling, intercellular communication facilitated by membrane potential and electrochemical coupling, is emerging as a key regulator of animal development. Gap junction (GJ) channels can mediate bioelectric signaling by creating a fast, direct pathway between cells for the movement of ions and other small molecules. In vertebrates, GJ channels are formed by a highly conserved transmembrane protein family called the connexins. The connexin gene family is large and complex, creating challenges in identifying specific connexins that create channels within developing and mature tissues. Using the embryonic zebrafish neuromuscular system as a model, we identify a connexin conserved across vertebrate lineages, gjd4, which encodes the Cx46.8 protein, that mediates bioelectric signaling required for slow muscle development and function. Through mutant analysis and in vivo imaging, we show that gjd4/Cx46.8 creates GJ channels specifically in developing slow muscle cells. Using genetics, pharmacology, and calcium imaging, we find that spinal-cord-generated neural activity is transmitted to developing slow muscle cells, and synchronized activity spreads via gjd4/Cx46.8 GJ channels. Finally, we show that bioelectrical signal propagation within the developing neuromuscular system is required for appropriate myofiber organization and that disruption leads to defects in behavior. Our work reveals a molecular basis for GJ communication among developing muscle cells and reveals how perturbations to bioelectric signaling in the neuromuscular system may contribute to developmental myopathies. Moreover, this work underscores a critical motif of signal propagation between organ systems and highlights the pivotal role of GJ communication in coordinating bioelectric signaling during development.

Disruption of the creb3l1 gene causes defects in caudal fin regeneration and patterning in zebrafish Danio rerio

BACKGROUND

The gene cAMP-Responsive Element Binding protein 3-like-1 (CREB3L1) has been implicated in bone development in mice, with CREB3L1 knock-out mice exhibiting fragile bones, and in humans, with CREB3L1 mutations linked to osteogenesis imperfecta. However, the mechanism through which Creb3l1 regulates bone development is not fully understood.

RESULTS

To probe the role of Creb3l1 in organismal physiology, we used CRISPR-Cas9 genome editing to generate a Danio rerio (zebrafish) model of Creb3l1 deficiency. In contrast to mammalian phenotypes, the Creb3l1 deficient fish do not display abnormalities in osteogenesis, except for a decrease in the bifurcation pattern of caudal fin. Both, skeletal morphology and overall bone density appear normal in the mutant fish. However, the regeneration of caudal fin postamputation is significantly affected, with decreased overall regenerate and mineralized bone area. Moreover, the mutant fish exhibit a severe patterning defect during regeneration, with a significant decrease in bifurcation complexity of the fin rays and distalization of the bifurcation sites. Analysis of genes implicated in bone development showed aberrant patterning of shha and ptch2 in Creb3l1 deficient fish, linking Creb3l1 with Sonic Hedgehog signaling during fin regeneration.

CONCLUSIONS

Our results uncover a novel role for Creb3l1 in regulating tissue growth and patterning during regeneration.

scMultiome analysis identifies a single caudal hindbrain compartment in the developing zebrafish nervous system

BACKGROUND

A key step in nervous system development involves the coordinated control of neural progenitor specification and positioning. A long-standing model for the vertebrate CNS postulates that transient anatomical compartments - known as neuromeres - function to position neural progenitors along the embryonic anteroposterior neuraxis. Such neuromeres are apparent in the embryonic hindbrain - that contains six rhombomeres with morphologically apparent boundaries - but other neuromeres lack clear morphological boundaries and have instead been defined by different criteria, such as differences in gene expression patterns and the outcomes of transplantation experiments. Accordingly, the caudal hindbrain (CHB) posterior to rhombomere (r) 6 has been variably proposed to contain from two to five 'pseudo-rhombomeres', but the lack of comprehensive molecular data has precluded a detailed definition of such structures.

METHODS

We used single-cell Multiome analysis, which allows simultaneous characterization of gene expression and chromatin state of individual cell nuclei, to identify and characterize CHB progenitors in the developing zebrafish CNS.

RESULTS

We identified CHB progenitors as a transcriptionally distinct population, that also possesses a unique profile of accessible transcription factor binding motifs, relative to both r6 and the spinal cord. This CHB population can be subdivided along its dorsoventral axis based on molecular characteristics, but we do not find any molecular evidence that it contains multiple pseudo-rhombomeres. We further observe that the CHB is closely related to r6 at the earliest embryonic stages, but becomes more divergent over time, and that it is defined by a unique gene regulatory network.

CONCLUSIONS

We conclude that the early CHB represents a single neuromere compartment that cannot be molecularly subdivided into pseudo-rhombomeres and that it may share an embryonic origin with r6.

Meningeal lymphatic supporting cells govern the formation and maintenance of zebrafish mural lymphatic endothelial cells

The meninges are critical for the brain functions, but the diversity of meningeal cell types and intercellular interactions have yet to be thoroughly examined. Here we identify a population of meningeal lymphatic supporting cells (mLSCs) in the zebrafish leptomeninges, which are specifically labeled by ependymin. Morphologically, mLSCs form membranous structures that enwrap the majority of leptomeningeal blood vessels and all the mural lymphatic endothelial cells (muLECs). Based on its unique cellular morphologies and transcriptional profile, mLSC is characterized as a unique cell type different from all the currently known meningeal cell types. Because of the formation of supportive structures and production of pro-lymphangiogenic factors, mLSCs not only promote muLEC development and maintain the dispersed distributions of muLECs in the leptomeninges, but also are required for muLEC regeneration after ablation. This study characterizes a newly identified cell type in leptomeninges, mLSC, which is required for muLEC development, maintenance, and regeneration.

Mapping the chromatin accessibility landscape of zebrafish embryogenesis at single-cell resolution by SPATAC-seq

Currently, the dynamic accessible elements that determine regulatory programs responsible for the unique identity and function of each cell type during zebrafish embryogenesis lack detailed study. Here we present SPATAC-seq: a split-pool ligation-based assay for transposase-accessible chromatin using sequencing. Using SPATAC-seq, we profiled chromatin accessibility in more than 800,000 individual nuclei across 20 developmental stages spanning the sphere stage to the early larval protruding mouth stage. Using this chromatin accessibility map, we identified 604 cell states and inferred their developmental relationships. We also identified 959,040 candidate cis-regulatory elements (cCREs) and delineated development-specific cCREs, as well as transcription factors defining diverse cell identities. Importantly, enhancer reporter assays confirmed that the majority of tested cCREs exhibited robust enhanced green fluorescent protein expression in restricted cell types or tissues. Finally, we explored gene regulatory programs that drive pigment and notochord cell differentiation. Our work provides a valuable open resource for exploring driver regulators of cell fate decisions in zebrafish embryogenesis.

A new atlas to study embryonic cell types in Xenopus

This paper introduces a single-cell atlas for pivotal developmental stages in Xenopus, encompassing gastrulation, neurulation, and early tailbud. Notably surpassing its predecessors, the new atlas enhances gene mapping, read counts, and gene/cell type nomenclature. Leveraging the latest Xenopus tropicalis genome version, alongside advanced alignment pipelines and machine learning for cell type assignment, this release maintains consistency with previous cell type annotations while rectifying nomenclature issues. Employing an unbiased approach for cell type assignment proves especially apt for embryonic contexts, given the considerable number of non-terminally differentiated cell types. An alternative cell type attribution here adopts a fuzzy, non-deterministic stance, capturing the transient nature of early embryo progenitor cells by presenting an ensemble of types in superposition. The value of the new resource is emphasized through numerous examples, with a focus on previously unexplored germ cell populations where we uncover novel transcription onset features. Offering interactive exploration via a user-friendly web portal and facilitating complete data downloads, this atlas serves as a comprehensive and accessible reference.

Loss of αBa-crystallin, but not αA-crystallin, increases age-related cataract in the zebrafish lens

The vertebrate eye lens is an unusual organ in that most of its cells lack nuclei and the ability to replace aging protein. The small heat shock protein α-crystallins evolved to become key components of this lens, possibly because of their ability to prevent aggregation of aging protein that would otherwise lead to lens opacity. Most vertebrates express two α-crystallins, αA- and αB-crystallin, and mutations in each are linked to human cataract. In a mouse knockout model only the loss of αA-crystallin led to early-stage lens cataract. We have used the zebrafish as a model system to investigate the role of α-crystallins during lens development. Interestingly, while zebrafish express one lens-specific αA-crystallin gene (cryaa), they express two αB-crystallin genes, with one evolving lens specificity (cryaba) and the other retaining the broad expression of its mammalian ortholog (cryabb). In this study we used individual mutant zebrafish lines for all three α-crystallin genes to determine the impact of their loss on age-related cataract. Surprisingly, unlike mouse knockout models, we found that the loss of the αBa-crystallin gene cryaba led to an increase in lens opacity compared to cryaa null fish at 24 months of age. Loss of αA-crystallin did not increase the prevalence of cataract. We also used single cell RNA-Seq and RT-qPCR data to show a shift in the lens expression of zebrafish α-crystallins between 5 and 10 days post fertilization (dpf), with 5 and 6 dpf lenses expressing cryaa almost exclusively, and expression of cryaba and cryabb becoming more prominent after 10 dpf. These data show that cryaa is the primary α-crystallin during early lens development, while the protective role for cryaba becomes more important during lens aging. This study is the first to quantify cataract prevalence in wild-type aging zebrafish, showing that lens opacities develop in approximately 25% of fish by 18 months of age. None of the three α-crystallin mutants showed a compensatory increase in the expression of the remaining two crystallins, or in the abundant βB1-crystallin. Overall, these findings indicate an ontogenetic shift in the functional importance of individual α-crystallins during zebrafish lens development. Our finding that the lens-specific zebrafish αBa-crystallin plays the leading role in preventing age-related cataract adds a new twist to our understanding of vertebrate lens evolution.

Zebrafish do not have calprotectin

The protein heterodimer calprotectin and its component proteins, S100A8 and S100A9, play important antibacterial and proinflammatory roles in the mammalian innate immune response. Gaining mechanistic insights into the regulation and biological function of calprotectin will help facilitate patient diagnostics and therapy and further our understanding of the host-microbe interface. Recent literature has identified zebrafish s100a10b as zebrafish calprotectin based on sequence similarity, genomic context, and transcriptional upregulation during the immune response to bacterial infections. The field would benefit from expanding the breadth of calprotectin studies into a zebrafish innate immunity model. Here, we carefully evaluated the possibility that zebrafish possess a calprotectin. We found that zebrafish do not possess an ortholog of mammalian S100A8 or S100A9. We then identified four zebrafish s100 proteins- including s100a10b-that are expressed in immune cells and upregulated during the immune response. We recombinantly expressed and purified these proteins and measured the antimicrobial and proinflammatory characteristics. We found that none of the zebrafish proteins exhibited activity comparable to mammalian calprotectin. Our work demonstrates conclusively that zebrafish have no ortholog of calprotectin, and the most plausible candidate proteins have not convergently evolved similar functions.

scRNA-seq reveals the diversity of the developing cardiac cell lineage and molecular players in heart rhythm regulation

We utilized scRNA-seq to delineate the diversity of cell types in the zebrafish heart. Transcriptome profiling of over 50,000 cells at 48 and 72 hpf defined at least 18 discrete cell lineages of the developing heart. Utilizing well-established gene signatures, we identified a population of cells likely to be the primary pacemaker and characterized the transcriptome profile defining this critical cell type. Two previously uncharacterized genes, atp1b3b and colec10, were found to be enriched in the sinoatrial cardiomyocytes. CRISPR/Cas9-mediated knockout of these two genes significantly reduced heart rate, implicating their role in cardiac development and conduction. Additionally, we describe other cardiac cell lineages, including the endothelial and neural cells, providing their expression profiles as a resource. Our results established a detailed atlas of the developing heart, providing valuable insights into cellular and molecular mechanisms, and pinpointed potential new players in heart rhythm regulation.

Lactate secreted by glycolytic conjunctival melanoma cells attracts and polarizes macrophages to drive angiogenesis in zebrafish xenografts

Conjunctival melanoma (CoM) is a rare but potentially lethal cancer of the eye, with limited therapeutic option for metastases. A better understanding how primary CoM disseminate to form metastases is urgently needed in order to develop novel therapies. Previous studies indicated that primary CoM tumors express Vascular Endothelial Growth Factor (VEGF) and may recruit pro-tumorigenic M2-like macrophages. However, due to a lack of proper models, the expected role of angiogenesis in the metastatic dissemination of CoM is still unknown. We show that cells derived from two CoM cell lines induce a strong angiogenic response when xenografted in zebrafish larvae. CoM cells are highly glycolytic and secrete lactate, which recruits and polarizes human and zebrafish macrophages towards a M2-like phenotype. These macrophages elevate the levels of proangiogenic factors such as VEGF, TGF-β, and IL-10 in the tumor microenvironment to induce an angiogenic response towards the engrafted CoM cells in vivo. Chemical ablation of zebrafish macrophages or inhibition of glycolysis in CoM cells terminates this response, suggesting that attraction of lactate-dependent macrophages into engrafted CoM cells drives angiogenesis and serves as a possible dissemination mechanism for glycolytic CoM cells.

Genetic requirement of dact1/2 to regulate noncanonical Wnt signaling and calpain 8 during embryonic convergent extension and craniofacial morphogenesis

Wnt signaling plays a crucial role in the early embryonic patterning and development, to regulate convergent extension during gastrulation and the establishment of the dorsal axis. Further, Wnt signaling is a crucial regulator of craniofacial morphogenesis. The adapter proteins Dact1 and Dact2 modulate the Wnt signaling pathway through binding to Disheveled. However, the distinct relative functions of Dact1 and Dact2 during embryogenesis remain unclear. We found that dact1 and dact2 genes have dynamic spatiotemporal expression domains that are reciprocal to one another and to wnt11f2, that suggest distinct functions during zebrafish embryogenesis. We found that both dact1 and dact2 contribute to axis extension, with compound mutants exhibiting a similar convergent extension defect and craniofacial phenotype to the wnt11f2 mutant. Utilizing single-cell RNAseq and gpc4 mutant that disrupts noncanonical Wnt signaling, we identified dact1/2 specific roles during early development. Comparative whole transcriptome analysis between wildtype, gpc4 and dact1/2 mutants revealed a novel role for dact1/2 in regulating the mRNA expression of the classical calpain capn8. Over-expression of capn8 phenocopies dact1/2 craniofacial dysmorphology. These results identify a previously unappreciated role of capn8 and calcium-dependent proteolysis during embryogenesis. Taken together, our findings highlight the distinct and overlapping roles of dact1 and dact2 in embryonic craniofacial development, providing new insights into the multifaceted regulation of Wnt signaling.

Status of single-cell RNA sequencing for reproductive toxicology in zebrafish and the transcriptomic trade-off

The utilization of transcriptomic studies identifying profiles of gene expression, especially in toxicogenomics, has catapulted next-generation sequencing to the forefront of reproductive toxicology. An innovative yet underutilized RNA sequencing technique emerging into this field is single-cell RNA sequencing (scRNA-seq), which provides sequencing at the individual cellular level of gonad tissue. ScRNA-seq provides a novel and unique perspective for identifying distinct cellular profiles, including identification of rare cell subtypes. The specificity of scRNA-seq is a powerful tool for reproductive toxicity research, especially for translational animal models including zebrafish. Studies to date not only have focused on 'tissue atlassing' or characterizing what cell types make up different tissues but have also begun to include toxicant exposure as a factor that this review aims to explore. Future scRNA-seq studies will contribute to understanding exposure-induced outcomes; however, the trade-offs with traditional methods need to be considered.

Zebrafish as Versatile Model for Assessing Animal Venoms and Toxins: Current Applications and Future Prospects

Animal venoms and toxins hold promise as sources of novel drug candidates, therapeutic agents, and biomolecules. To fully harness their potential, it is crucial to develop reliable testing methods that provide a comprehensive understanding of their effects and mechanisms of action. However, traditional rodent assays encounter difficulties in mimicking venom-induced effects in human due to the impractical venom dosage levels. The search for reliable testing methods has led to the emergence of zebrafish (Danio rerio) as a versatile model organism for evaluating animal venoms and toxins. Zebrafish possess genetic similarities to humans, rapid development, transparency, and amenability to high-throughput assays, making it ideal for assessing the effects of animal venoms and toxins. This review highlights unique attributes of zebrafish and explores their applications in studying venom- and toxin-induced effects from various species, including snakes, jellyfish, cuttlefish, anemones, spiders, and cone snails. Through zebrafish-based research, intricate physiological responses, developmental alterations, and potential therapeutic interventions induced by venoms are revealed. Novel techniques such as CRISPR/Cas9 gene editing, optogenetics, and high-throughput screening hold great promise for advancing venom research. As zebrafish-based insights converge with findings from other models, the comprehensive understanding of venom-induced effects continues to expand, guiding the development of targeted interventions and promoting both scientific knowledge and practical applications.

Transcriptional background effects on a tumor driver gene in different pigment cell types of medaka

The Xiphophorus melanoma receptor kinase gene, xmrk, is a bona fide oncogene driving melanocyte tumorigenesis of Xiphophorus fish. When ectopically expressed in medaka, it not only induces development of several pigment cell tumor types in different strains of medaka but also induces different tumor types within the same animal, suggesting its oncogenic activity has a transcriptomic background effect. Although the central pathways that xmrk utilizes to lead to melanomagenesis are well documented, genes and genetic pathways that modulate the oncogenic effect and alter the course of disease have not been studied so far. To understand how the genetic networks between different histocytes of xmrk-driven tumors are composed, we isolated two types of tumors, melanoma and xanthoerythrophoroma, from the same xmrk transgenic medaka individuals, established the transcriptional profiles of both xmrk-driven tumors, and compared (1) genes that are co-expressed with xmrk in both tumor types, and (2) differentially expressed genes and their associated molecular functions, between the two tumor types. Transcriptomic comparisons between the two tumor types show melanoma and xanthoerythrophoroma are characterized by transcriptional features representing varied functions, indicating distinct molecular interactions between the driving oncogene and the cell-type-specific transcriptomes. Melanoma tumors exhibit gene signatures that are relevant to proliferation and invasion, while xanthoerythrophoroma tumors are characterized by expression profiles related to metabolism and DNA repair. We conclude the transcriptomic backgrounds, exemplified by cell-type-specific genes that are downstream of xmrk effected signaling pathways, contribute the potential to change the course of tumor development and may affect overall tumor outcomes.

An in vivo drug screen in zebrafish reveals that cyclooxygenase 2-derived prostaglandin D2 promotes spinal cord neurogenesis

The study of neurogenesis is essential to understanding fundamental developmental processes and for the development of cell replacement therapies for central nervous system disorders. Here, we designed an in vivo drug screening protocol in developing zebrafish to find new molecules and signalling pathways regulating neurogenesis in the ventral spinal cord. This unbiased drug screen revealed that 4 cyclooxygenase (COX) inhibitors reduced the generation of serotonergic interneurons in the developing spinal cord. These results fitted very nicely with available single-cell RNAseq data revealing that floor plate cells show differential expression of 1 of the 2 COX2 zebrafish genes (ptgs2a). Indeed, several selective COX2 inhibitors and two different morpholinos against ptgs2a reduced the number of serotonergic neurons in the ventral spinal cord and led to locomotor deficits. Single-cell RNAseq data and different pharmacological manipulations further revealed that COX2-floor plate-derived prostaglandin D2 promotes neurogenesis in the developing spinal cord by promoting mitotic activity in progenitor cells. Rescue experiments using a phosphodiesterase-4 inhibitor suggest that intracellular changes in cAMP levels underlie the effects of COX inhibitors on neurogenesis and locomotion. Our study provides compelling in vivo evidence showing that prostaglandin signalling promotes neurogenesis in the ventral spinal cord.

Single-Cell Analysis of Rohon-Beard Neurons Implicates Fgf Signaling in Axon Maintenance and Cell Survival

Peripheral sensory neurons are a critical part of the nervous system that transmit a multitude of sensory stimuli to the central nervous system. During larval and juvenile stages in zebrafish, this function is mediated by Rohon-Beard somatosensory neurons (RBs). RBs are optically accessible and amenable to experimental manipulation, making them a powerful system for mechanistic investigation of sensory neurons. Previous studies provided evidence that RBs fall into multiple subclasses; however, the number and molecular makeup of these potential RB subtypes have not been well defined. Using a single-cell RNA sequencing (scRNA-seq) approach, we demonstrate that larval RBs in zebrafish fall into three, largely nonoverlapping classes of neurons. We also show that RBs are molecularly distinct from trigeminal neurons in zebrafish. Cross-species transcriptional analysis indicates that one RB subclass is similar to a mammalian group of A-fiber sensory neurons. Another RB subclass is predicted to sense multiple modalities, including mechanical stimulation and chemical irritants. We leveraged our scRNA-seq data to determine that the fibroblast growth factor (Fgf) pathway is active in RBs. Pharmacological and genetic inhibition of this pathway led to defects in axon maintenance and RB cell death. Moreover, this can be phenocopied by treatment with dovitinib, an FDA-approved Fgf inhibitor with a common side effect of peripheral neuropathy. Importantly, dovitinib-mediated axon loss can be suppressed by loss of Sarm1, a positive regulator of neuronal cell death and axonal injury. This offers a molecular target for future clinical intervention to fight neurotoxic effects of this drug.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Single Cell RNA Sequencing Provides Clues for the Developmental Genetic Basis of Syngnathidae's Evolutionary Adaptations

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 genomes 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.

Spatially resolved single-cell atlas of ascidian endostyle provides insight into the origin of vertebrate pharyngeal organs

The pharyngeal endoderm, an innovation of deuterostome ancestors, contributes to pharyngeal development by influencing the patterning and differentiation of pharyngeal structures in vertebrates; however, the evolutionary origin of the pharyngeal organs in vertebrates is largely unknown. The endostyle, a distinct pharyngeal organ exclusively present in basal chordates, represents a good model for understanding pharyngeal organ origins. Using Stereo-seq and single-cell RNA sequencing, we constructed aspatially resolved single-cell atlas for the endostyle of the ascidian Styela clava. We determined the cell composition of the hemolymphoid region, which illuminates a mixed ancestral structure for the blood and lymphoid system. In addition, we discovered a cluster of hair cell-like cells in zone 3, which has transcriptomic similarity with the hair cells of the vertebrate acoustico-lateralis system. These findings reshape our understanding of the pharynx of the basal chordate and provide insights into the evolutionary origin of multiplexed pharyngeal organs.

A molecular atlas of adult C. elegans motor neurons reveals ancient diversity delineated by conserved transcription factor codes

Motor neurons (MNs) constitute an ancient cell type targeted by multiple adult-onset diseases. It is therefore important to define the molecular makeup of adult MNs in animal models and extract organizing principles. Here, we generate a comprehensive molecular atlas of adult Caenorhabditis elegans MNs and a searchable database. Single-cell RNA sequencing of 13,200 cells reveals that ventral nerve cord MNs cluster into 29 molecularly distinct subclasses. Extending C. elegans Neuronal Gene Expression Map and Network (CeNGEN) findings, all MN subclasses are delineated by distinct expression codes of either neuropeptide or transcription factor gene families. Strikingly, combinatorial codes of homeodomain transcription factor genes succinctly delineate adult MN diversity in both C. elegans and mice. Further, molecularly defined MN subclasses in C. elegans display distinct patterns of connectivity. Hence, our study couples the connectivity map of the C. elegans motor circuit with a molecular atlas of its constituent MNs and uncovers organizing principles and conserved molecular codes of adult MN diversity.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Paternal starvation affects metabolic gene expression during zebrafish offspring development and lifelong fitness

Dietary restriction in the form of fasting is a putative key to a healthier and longer life, but these benefits may come at a trade-off with reproductive fitness and may affect the following generation(s). The potential inter- and transgenerational effects of long-term fasting and starvation are particularly poorly understood in vertebrates when they originate from the paternal line. We utilised the externally fertilising zebrafish amenable to a split-egg clutch design to explore the male-specific effects of fasting/starvation on fertility and fitness of offspring independently of maternal contribution. Eighteen days of fasting resulted in reduced fertility in exposed males. While average offspring survival was not affected, we detected increased larval growth rate in F1 offspring from starved males and more malformed embryos at 24 h post-fertilisation in F2 offspring produced by F1 offspring from starved males. Comparing the transcriptomes of F1 embryos sired by starved and fed fathers revealed robust and reproducible increased expression of muscle composition genes but lower expression of lipid metabolism and lysosome genes in embryos from starved fathers. A large proportion of these genes showed enrichment in the yolk syncytial layer suggesting gene regulatory responses associated with metabolism of nutrients through paternal effects on extra-embryonic tissues which are loaded with maternal factors. We compared the embryo transcriptomes to published adult transcriptome datasets and found comparable repressive effects of starvation on metabolism-associated genes. These similarities suggest a physiologically relevant, directed and potentially adaptive response transmitted by the father, independently from the offspring's nutritional state, which was defined by the mother.

Mapping the cellular landscape of Atlantic salmon head kidney by single cell and single nucleus transcriptomics

Single-cell transcriptomics is the current gold standard for global gene expression profiling, not only in mammals and model species, but also in non-model fish species. This is a rapidly expanding field, creating a deeper understanding of tissue heterogeneity and the distinct functions of individual cells, making it possible to explore the complexities of immunology and gene expression on a highly resolved level. In this study, we compared two single cell transcriptomic approaches to investigate cellular heterogeneity within the head kidney of healthy farmed Atlantic salmon (Salmo salar). We compared 14,149 cell transcriptomes assayed by single cell RNA-seq (scRNA-seq) with 18,067 nuclei transcriptomes captured by single nucleus RNA-Seq (snRNA-seq). Both approaches detected eight major cell populations in common: granulocytes, heamatopoietic stem cells, erythrocytes, mononuclear phagocytes, thrombocytes, B cells, NK-like cells, and T cells. Four additional cell types, endothelial, epithelial, interrenal, and mesenchymal cells, were detected in the snRNA-seq dataset, but appeared to be lost during preparation of the single cell suspension submitted for scRNA-seq library generation. We identified additional heterogeneity and subpopulations within the B cells, T cells, and endothelial cells, and revealed developmental trajectories of heamatopoietic stem cells into differentiated granulocyte and mononuclear phagocyte populations. Gene expression profiles of B cell subtypes revealed distinct IgM and IgT-skewed resting B cell lineages and provided insights into the regulation of B cell lymphopoiesis. The analysis revealed eleven T cell sub-populations, displaying a level of T cell heterogeneity in salmon head kidney comparable to that observed in mammals, including distinct subsets of cd4/cd8-negative T cells, such as tcrγ positive, progenitor-like, and cytotoxic cells. Although snRNA-seq and scRNA-seq were both useful to resolve cell type-specific expression in the Atlantic salmon head kidney, the snRNA-seq pipeline was overall more robust in identifying several cell types and subpopulations. While scRNA-seq displayed higher levels of ribosomal and mitochondrial genes, snRNA-seq captured more transcription factor genes. However, only scRNA-seq-generated data was useful for cell trajectory inference within the myeloid lineage. In conclusion, this study systematically outlines the relative merits of scRNA-seq and snRNA-seq in Atlantic salmon, enhances understanding of teleost immune cell lineages, and provides a comprehensive list of markers for identifying major cell populations in the head kidney with significant immune relevance.

Zebrafish and cellular models of SELENON-Related Myopathy exhibit novel embryonic and metabolic phenotypes

SELENON-Related Myopathy (SELENON-RM) is a rare congenital myopathy caused by mutations of the SELENON gene characterized by axial muscle weakness and progressive respiratory insufficiency. Muscle histopathology commonly includes multiminicores or a dystrophic pattern but is often non-specific. 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 causes SELENON-RM are undetermined. A hurdle is the lack of cellular and animal models that show assayable phenotypes. Here 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 between selenon and genes involved in glutathione redox pathway. SelN-deficient zebrafish and mouse myoblasts exhibit changes in glutathione and redox homeostasis, suggesting a direct relationship with SelN function. We report changes in metabolic function abnormalities in SelN-null myotubes when compared to WT. These results suggest that SelN has functional roles during zebrafish early development and myoblast metabolism.

Transcriptomic and behavioral analyses reveal unique target tissues and molecular pathways associated with embryonic exposure to low level glyphosate and metal mixtures in zebrafish

Investigation of developmental molecular events following exposure to environmentally relevant agrochemical mixtures is critical to predicting their potential long-term ecological and human health risks. Here, we sought to uncover transcriptomic changes during zebrafish (Danio rerio) embryonic development following exposure to glyphosate and co-exposure to metals. Glyphosate is widely used globally with an allowable drinking water limit of 700 ppb. We examined effects of glyphosate (10 ppb) alone and when co-exposed to a metal mixture containing low levels of arsenic (4 ppb), lead (5 ppb), cadmium (2 ppb), and vanadium (15 ppb). This mixture was derived based on behavioral and morphological toxicity findings and environmentally relevant concentrations found in agricultural regions where glyphosate and metals are ubiquitously present. Gene expression patterns coupled to a single-cell transcriptomic dataset revealed that developmental exposure (28-72 h post fertilization) to glyphosate dysregulates expression of developmental genes specific to the central nervous system. Subsequent studies indicated significant suppression of larval zebrafish movement with 10 ppb glyphosate exposure. Studies with glyphosate + metals mixture and metals mixture alone showed unique developmental transcriptomic patterns and behavioral changes compared to glyphosate exposure alone. However, some outcomes (e.g., changes in expression of genes involved in epigenetic regulation and extracellular matrix patterning) were common across all three exposures compared to the control. Notably, glyphosate + metals co-exposure distinctly suppresses lysosomal transcripts and targets renal developmental genes. While further studies are required to uncover the precise nature of the interactions between glyphosate and metals, our study shows that glyphosate at very low levels is a behavioral and neurotoxicant that changes when metals are present. Given this herbicide affects distinctive physiological processes, including renal development and lysosomal dysregulation when co-exposed with metals, we conclude that environmental cation levels should be considered in glyphosate toxicity and risk assessment.

Bacterial phenotypic heterogeneity through the lens of single-cell RNA sequencing

Bacterial transcription is not monolithic. Microbes exist in a wide variety of cell states that help them adapt to their environment, acquire and produce essential nutrients, and engage in both competition and cooperation with their neighbors. While we typically think of bacterial adaptation as a group behavior, where all cells respond in unison, there is often a mixture of phenotypic responses within a bacterial population, where distinct cell types arise. A primary phenomenon driving these distinct cell states is transcriptional heterogeneity. Given that bacterial mRNA transcripts are extremely short-lived compared to eukaryotes, their transcriptional state is closely associated with their physiology, and thus the transcriptome of a bacterial cell acts as a snapshot of the behavior of that bacterium. Therefore, the application of single-cell transcriptomics to microbial populations will provide novel insight into cellular differentiation and bacterial ecology. In this review, we provide an overview of transcriptional heterogeneity in microbial systems, discuss the findings already provided by single-cell approaches, and plot new avenues of inquiry in transcriptional regulation, cellular biology, and mechanisms of heterogeneity that are made possible when microbial communities are analyzed at single-cell resolution.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Motor neurons are dispensable for the assembly of a sensorimotor circuit for gaze stabilization

Sensorimotor reflex circuits engage distinct neuronal subtypes, defined by precise connectivity, to transform sensation into compensatory behavior. Whether and how motor neuron populations specify the subtype fate and/or sensory connectivity of their pre-motor partners remains controversial. Here, we discovered that motor neurons are dispensable for proper connectivity in the vestibular reflex circuit that stabilizes gaze. We first measured activity following vestibular sensation in pre-motor projection neurons after constitutive loss of their extraocular motor neuron partners. We observed normal responses and topography indicative of unchanged functional connectivity between sensory neurons and projection neurons. Next, we show that projection neurons remain anatomically and molecularly poised to connect appropriately with their downstream partners. Lastly, we show that the transcriptional signatures that typify projection neurons develop independently of motor partners. Our findings comprehensively overturn a long-standing model: that connectivity in the circuit for gaze stabilization is retrogradely determined by motor partner-derived signals. By defining the contribution of motor neurons to specification of an archetypal sensorimotor circuit, our work speaks to comparable processes in the spinal cord and advances our understanding of general principles of neural development.

Bisphenol A Analogues Induce Neuroendocrine Disruption via Gut-Brain Regulation in Zebrafish

There is epidemiological evidence in humans that exposure to endocrine-disrupting chemicals such as bisphenol A (BPA) is tied to abnormal neuroendocrine function with both behavioral and intestinal symptoms. However, the underlying mechanism of this effect, particularly the role of gut-brain regulation, is poorly understood. We exposed zebrafish embryos to a concentration series (including environmentally relevant levels) of BPA and its analogues. The analogue bisphenol G (BPG) yielded the strongest behavioral impact on zebrafish larvae and inhibited the largest number of neurotransmitters, with an effective concentration of 0.5 μg/L, followed by bisphenol AF (BPAF) and BPA. In neurod1:EGFP transgenic zebrafish, BPG and BPAF inhibited the distribution of enteroendocrine cells (EECs), which is associated with decreased neurotransmitters level and behavioral activity. Immune staining of ace-α-tubulin suggested that BPAF inhibited vagal neural development at 50 and 500 μg/L. Single-cell RNA-Seq demonstrated that BPG disrupted the neuroendocrine system by inducing inflammatory responses in intestinal epithelial cells via TNFα-trypsin-EEC signaling. BPAF exposure activated apoptosis and inhibited neural developmental pathways in vagal neurons, consistent with immunofluorescence imaging studies. These findings show that both BPG and BPAF affect the neuroendocrine system through the gut-brain axis but by different mechanisms, revealing new insights into the modes of bisphenol-mediated neuroendocrine disruption.

Dynamics of macrophage polarization support Salmonella persistence in a whole living organism

Numerous intracellular bacterial pathogens interfere with macrophage function, including macrophage polarization, to establish a niche and persist. However, the spatiotemporal dynamics of macrophage polarization during infection within host remain to be investigated. Here, we implement a model of persistent Salmonella Typhimurium infection in zebrafish, which allows visualization of polarized macrophages and bacteria in real time at high resolution. While macrophages polarize toward M1-like phenotype to control early infection, during later stages, Salmonella persists inside non-inflammatory clustered macrophages. Transcriptomic profiling of macrophages showed a highly dynamic signature during infection characterized by a switch from pro-inflammatory to anti-inflammatory/pro-regenerative status and revealed a shift in adhesion program. In agreement with this specific adhesion signature, macrophage trajectory tracking identifies motionless macrophages as a permissive niche for persistent Salmonella. Our results demonstrate that zebrafish model provides a unique platform to explore, in a whole organism, the versatile nature of macrophage functional programs during bacterial acute and persistent infections.

Astrocyte growth is driven by the Tre1/S1pr1 phospholipid-binding G protein-coupled receptor

Astrocytes play crucial roles in regulating neural circuit function by forming a dense network of synapse-associated membrane specializations, but signaling pathways regulating astrocyte morphogenesis remain poorly defined. Here, we show the Drosophila lipid-binding G protein-coupled receptor (GPCR) Tre1 is required for astrocytes to establish their intricate morphology in vivo. The lipid phosphate phosphatases Wunen/Wunen2 also regulate astrocyte morphology and, via Tre1, mediate astrocyte-astrocyte competition for growth-promoting lipids. Loss of s1pr1, the functional analog of Tre1 in zebrafish, disrupts astrocyte process elaboration, and live imaging and pharmacology demonstrate that S1pr1 balances proper astrocyte process extension/retraction dynamics during growth. Loss of Tre1 in flies or S1pr1 in zebrafish results in defects in simple assays of motor behavior. Tre1 and S1pr1 are thus potent evolutionarily conserved regulators of the elaboration of astrocyte morphological complexity and, ultimately, astrocyte control of behavior.

Loss of αBa-crystallin, but not αA-crystallin, increases age-related cataract in the zebrafish lens

The vertebrate eye lens is an unusual organ in that most of its cells lack nuclei and the ability to replace aging protein. The small heat shock protein α-crystallins evolved to become key components of this lens, possibly because of their ability to prevent aggregation of aging protein that would otherwise lead to lens opacity. Most vertebrates express two α-crystallins, αA- and αB-crystallin, and mutations in each are linked to human cataract. In a mouse knockout model only the loss of αA-crystallin led to early-stage lens cataract. We have used the zebrafish as a model system to investigate the role of α-crystallins during lens development. Interestingly, while zebrafish express one lens-specific αA-crystallin gene (cryaa), they express two αB-crystallin genes, with one evolving lens specificity (cryaba) and the other retaining the broad expression of its mammalian ortholog (cryabb). In this study we used individual mutant zebrafish lines for all three α-crystallin genes to determine the impact of their loss on age-related cataract. Surprisingly, unlike mouse knockout models, we found that the loss of the αBa-crystallin gene cryaba led to an increase in lens opacity compared to cryaa null fish at 24 months of age. Loss of αA-crystallin did not increase the prevalence of cataract. We also used single cell RNA-Seq and RT-qPCR data to show a shift in the lens expression of zebrafish α-crystallins between 5 and 10 days post fertilization (dpf), with 5 and 6 dpf lenses expressing cryaa almost exclusively, and expression of cryaba and cryabb becoming more prominent after 10 dpf. These data show that cryaa is the primary α-crystallin during early lens development, while the protective role for cryaba becomes more important during lens aging. This study is the first to quantify cataract prevalence in wild-type zebrafish, showing that lens opacities develop in approximately 25% of fish by 18 months of age. None of the three α-crystallin mutants showed a compensatory increase in the expression of the remaining two crystallins, or in the abundant βB1-crystallin. Overall, these findings indicate an ontogenetic shift in the functional importance of individual α-crystallins during zebrafish lens development. Our finding that the lens-specific zebrafish αBa-crystallin plays the leading role in preventing age-related cataract adds a new twist to our understanding of vertebrate lens evolution.

RhoA GEF Mcf2lb regulates rosette integrity during collective cell migration

Multicellular rosettes are transient epithelial structures that serve as important cellular intermediates in the formation of diverse organs. Using the zebrafish posterior lateral line primordium (pLLP) as a model system, we investigated the role of the RhoA GEF Mcf2lb in rosette morphogenesis. The pLLP is a group of ∼150 cells that migrates along the zebrafish trunk and is organized into epithelial rosettes; these are deposited along the trunk and will differentiate into sensory organs called neuromasts (NMs). Using single-cell RNA-sequencing and whole-mount in situ hybridization, we showed that mcf2lb is expressed in the pLLP during migration. Live imaging and subsequent 3D analysis of mcf2lb mutant pLLP cells showed disrupted apical constriction and subsequent rosette organization. This resulted in an excess number of deposited NMs along the trunk of the zebrafish. Cell polarity markers ZO-1 and Par-3 were apically localized, indicating that pLLP cells are properly polarized. In contrast, RhoA activity, as well as signaling components downstream of RhoA, Rock2a and non-muscle Myosin II, were diminished apically. Thus, Mcf2lb-dependent RhoA activation maintains the integrity of epithelial rosettes.

The zebrafish as a potential model for vaccine and adjuvant development

INTRODUCTION

Zebrafishes represent a proven model for human diseases and systems biology, exhibiting physiological and genetic similarities and having innate and adaptive immune systems. However, they are underexplored for human vaccinology, vaccine development, and testing. Here we summarize gaps and challenges.

AREAS COVERED

Zebrafish models have four potential applications: 1) Vaccine safety: The past successes in using zebrafishes to test xenobiotics could extend to vaccine and adjuvant formulations for general safety or target organs due to the zebrafish embryos' optical transparency. 2) Innate immunity: The zebrafish offers refined ways to examine vaccine effects through signaling via Toll-like or NOD-like receptors in zebrafish myeloid cells. 3) Adaptive immunity: Zebrafishes produce IgM, IgD,and two IgZ immunoglobulins, but these are understudied, due to a lack of immunological reagents for challenge studies. 4) Systems vaccinology: Due to the availability of a well-referenced zebrafish genome, transcriptome, proteome, and epigenome, this model offers potential here.

EXPERT OPINION

It remains unproven whether zebrafishes can be employed for testing and developing human vaccines. We are still at the hypothesis-generating stage, although it is possible to begin outlining experiments for this purpose. Through transgenic manipulation, zebrafish models could offer new paths for shaping animal models and systems vaccinology.

Gap junction mediated bioelectric coordination is required for slow muscle development, organization, and function

Bioelectrical signaling, intercellular communication facilitated by membrane potential and electrochemical coupling, is emerging as a key regulator of animal development. Gap junction (GJ) channels can mediate bioelectric signaling by creating a fast, direct pathway between cells for the movement of ions and other small molecules. In vertebrates, GJ channels are formed by a highly conserved transmembrane protein family called the Connexins. The connexin gene family is large and complex, presenting a challenge in identifying the specific Connexins that create channels within developing and mature tissues. Using the embryonic zebrafish neuromuscular system as a model, we identify a connexin conserved across vertebrate lineages, gjd4, which encodes the Cx46.8 protein, that mediates bioelectric signaling required for appropriate slow muscle development and function. Through a combination of mutant analysis and in vivo imaging we show that gjd4/Cx46.8 creates GJ channels specifically in developing slow muscle cells. Using genetics, pharmacology, and calcium imaging we find that spinal cord generated neural activity is transmitted to developing slow muscle cells and synchronized activity spreads via gjd4/Cx46.8 GJ channels. Finally, we show that bioelectrical signal propagation within the developing neuromuscular system is required for appropriate myofiber organization, and that disruption leads to defects in behavior. Our work reveals the molecular basis for GJ communication among developing muscle cells and reveals how perturbations to bioelectric signaling in the neuromuscular system_may contribute to developmental myopathies. Moreover, this work underscores a critical motif of signal propagation between organ systems and highlights the pivotal role played by GJ communication in coordinating bioelectric signaling during development.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei

The innate immune system is the first line of defense in multicellular organisms. Danio rerio is widely considered a promising model for IIS-related research, with the most amount of scRNAseq data available among Teleostei. We summarized the scRNAseq and spatial transcriptomics experiments related to the IIS for zebrafish and other Teleostei from the GEO NCBI and the Single-Cell Expression Atlas. We found a considerable number of scRNAseq experiments at different stages of zebrafish development in organs such as the kidney, liver, stomach, heart, and brain. These datasets could be further used to conduct large-scale meta-analyses and to compare the IIS of zebrafish with the mammalian one. However, only a small number of scRNAseq datasets are available for other fish (turbot, salmon, cavefish, and dark sleeper). Since fish biology is very diverse, it would be a major mistake to use zebrafish alone in fish immunology studies. In particular, there is a special need for new scRNAseq experiments involving nonmodel Teleostei, e.g., long-lived species, cancer-resistant fish, and various fish ecotypes.

SorCS2 binds progranulin to regulate motor neuron development

Motor neuron (MN) development and nerve regeneration requires orchestrated action of a vast number of molecules. Here, we identify SorCS2 as a progranulin (PGRN) receptor that is required for MN diversification and axon outgrowth in zebrafish and mice. In zebrafish, SorCS2 knockdown also affects neuromuscular junction morphology and fish motility. In mice, SorCS2 and PGRN are co-expressed by newborn MNs from embryonic day 9.5 until adulthood. Using cell-fate tracing and nerve segmentation, we find that SorCS2 deficiency perturbs cell-fate decisions of brachial MNs accompanied by innervation deficits of posterior nerves. Additionally, adult SorCS2 knockout mice display slower motor nerve regeneration. Interestingly, primitive macrophages express high levels of PGRN, and their interaction with SorCS2-positive motor axon is required during axon pathfinding. We further show that SorCS2 binds PGRN to control its secretion, signaling, and conversion into granulins. We propose that PGRN-SorCS2 signaling controls MN development and regeneration in vertebrates.

Single-cell analysis reveals distinct fibroblast plasticity during tenocyte regeneration in zebrafish

Despite their importance in tissue maintenance and repair, fibroblast diversity and plasticity remain poorly understood. Using single-cell RNA sequencing, we uncover distinct sclerotome-derived fibroblast populations in zebrafish, including progenitor-like perivascular/interstitial fibroblasts, and specialized fibroblasts such as tenocytes. To determine fibroblast plasticity in vivo, we develop a laser-induced tendon ablation and regeneration model. Lineage tracing reveals that laser-ablated tenocytes are quickly regenerated by preexisting fibroblasts. By combining single-cell clonal analysis and live imaging, we demonstrate that perivascular/interstitial fibroblasts actively migrate to the injury site, where they proliferate and give rise to new tenocytes. By contrast, perivascular fibroblast-derived pericytes or specialized fibroblasts, including tenocytes, exhibit no regenerative plasticity. Active Hedgehog (Hh) signaling is required for the proliferation of activated fibroblasts to ensure efficient tenocyte regeneration. Together, our work highlights the functional diversity of fibroblasts and establishes perivascular/interstitial fibroblasts as tenocyte progenitors that promote tendon regeneration in a Hh signaling-dependent manner.

scMultiome analysis identifies embryonic hindbrain progenitors with mixed rhombomere identities

Rhombomeres serve to position neural progenitors in the embryonic hindbrain, thereby ensuring appropriate neural circuit formation, but the molecular identities of individual rhombomeres and the mechanism whereby they form has not been fully established. Here, we apply scMultiome analysis in zebrafish to molecularly resolve all rhombomeres for the first time. We find that rhombomeres become molecularly distinct between 10hpf (end of gastrulation) and 13hpf (early segmentation). While the embryonic hindbrain transiently contains alternating odd- versus even-type rhombomeres, our scMultiome analyses do not detect extensive odd versus even molecular characteristics in the early hindbrain. Instead, we find that each rhombomere displays a unique gene expression and chromatin profile. Prior to the appearance of distinct rhombomeres, we detect three hindbrain progenitor clusters (PHPDs) that correlate with the earliest visually observed segments in the hindbrain primordium that represent prospective rhombomere r2/r3 (possibly including r1), r4, and r5/r6, respectively. We further find that the PHPDs form in response to Fgf and RA morphogens and that individual PHPD cells co-express markers of multiple mature rhombomeres. We propose that the PHPDs contain mixed-identity progenitors and that their subdivision into individual rhombomeres requires the resolution of mixed transcription and chromatin states.

Unique activities of two overlapping PAX6 retinal enhancers

Enhancers play a critical role in development by precisely modulating spatial, temporal, and cell type-specific gene expression. Sequence variants in enhancers have been implicated in diseases; however, establishing the functional consequences of these variants is challenging because of a lack of understanding of precise cell types and developmental stages where the enhancers are normally active. PAX6 is the master regulator of eye development, with a regulatory landscape containing multiple enhancers driving the expression in the eye. Whether these enhancers perform additive, redundant or distinct functions is unknown. Here, we describe the precise cell types and regulatory activity of two PAX6 retinal enhancers, HS5 and NRE. Using a unique combination of live imaging and single-cell RNA sequencing in dual enhancer-reporter zebrafish embryos, we uncover differences in the spatiotemporal activity of these enhancers. Our results show that although overlapping, these enhancers have distinct activities in different cell types and therefore likely nonredundant functions. This work demonstrates that unique cell type-specific activities can be uncovered for apparently similar enhancers when investigated at high resolution in vivo.

Embryo-scale reverse genetics at single-cell resolution

The maturation of single-cell transcriptomic technologies has facilitated the generation of comprehensive cellular atlases from whole embryos1-4. A majority of these data, however, has been collected from wild-type embryos without an appreciation for the latent variation that is present in development. Here we present the 'zebrafish single-cell atlas of perturbed embryos': single-cell transcriptomic data from 1,812 individually resolved developing zebrafish embryos, encompassing 19 timepoints, 23 genetic perturbations and a total of 3.2 million cells. The high degree of replication in our study (eight or more embryos per condition) enables us to estimate the variance in cell type abundance organism-wide and to detect perturbation-dependent deviance in cell type composition relative to wild-type embryos. Our approach is sensitive to rare cell types, resolving developmental trajectories and genetic dependencies in the cranial ganglia neurons, a cell population that comprises less than 1% of the embryo. Additionally, time-series profiling of individual mutants identified a group of brachyury-independent cells with strikingly similar transcriptomes to notochord sheath cells, leading to new hypotheses about early origins of the skull. We anticipate that standardized collection of high-resolution, organism-scale single-cell data from large numbers of individual embryos will enable mapping of the genetic dependencies of zebrafish cell types, while also addressing longstanding challenges in developmental genetics, including the cellular and transcriptional plasticity underlying phenotypic diversity across individuals.

Electrical synapse structure requires distinct isoforms of a postsynaptic scaffold

Electrical synapses are neuronal gap junction (GJ) channels associated with a macromolecular complex called the electrical synapse density (ESD), which regulates development and dynamically modifies electrical transmission. However, the proteomic makeup and molecular mechanisms utilized by the ESD that direct electrical synapse formation are not well understood. Using the Mauthner cell of zebrafish as a model, we previously found that the intracellular scaffolding protein ZO1b is a member of the ESD, localizing postsynaptically, where it is required for GJ channel localization, electrical communication, neural network function, and behavior. Here, we show that the complexity of the ESD is further diversified by the genomic structure of the ZO1b gene locus. The ZO1b gene is alternatively initiated at three transcriptional start sites resulting in isoforms with unique N-termini that we call ZO1b-Alpha, -Beta, and -Gamma. We demonstrate that ZO1b-Beta and ZO1b-Gamma are broadly expressed throughout the nervous system and localize to electrical synapses. By contrast, ZO1b-Alpha is expressed mainly non-neuronally and is not found at synapses. We generate mutants in all individual isoforms, as well as double mutant combinations in cis on individual chromosomes, and find that ZO1b-Beta is necessary and sufficient for robust GJ channel localization. ZO1b-Gamma, despite its localization to the synapse, plays an auxiliary role in channel localization. This study expands the notion of molecular complexity at the ESD, revealing that an individual genomic locus can contribute distinct isoforms to the macromolecular complex at electrical synapses. Further, independent scaffold isoforms have differential contributions to developmental assembly of the interneuronal GJ channels. We propose that ESD molecular complexity arises both from the diversity of unique genes and from distinct isoforms encoded by single genes. Overall, ESD proteomic diversity is expected to have critical impacts on the development, structure, function, and plasticity of electrical transmission.

Widespread sex dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism for aging research. Here, we describe a multitissue, single-cell gene expression atlas of female and male blood, kidney, liver, and spleen. We annotate 22 cell types, define marker genes, and infer differentiation trajectories. We find pervasive sex-dimorphic gene expression across cell types. Sex-dimorphic genes tend to be linked to lipid metabolism, consistent with clear differences in lipid storage in female vs. male turquoise killifish livers. We use machine learning to predict sex using single-cell gene expression and identify potential markers for molecular sex identity. As a proof of principle, we show that our atlas can be used to deconvolute existing bulk RNA sequencing (RNA-seq) data to obtain accurate estimates of cell type proportions. This atlas can be a resource to the community that could be leveraged to develop cell-type-specific expression in transgenic animals.

Insulin-like growth factor receptor / mTOR signaling elevates global translation to accelerate zebrafish fin regenerative outgrowth

Zebrafish robustly regenerate fins, including their characteristic bony ray skeleton. Amputation activates intra-ray fibroblasts and dedifferentiates osteoblasts that migrate under a wound epidermis to establish an organized blastema. Coordinated proliferation and re-differentiation across lineages then sustains progressive outgrowth. We generate a single cell transcriptome dataset to characterize regenerative outgrowth and explore coordinated cell behaviors. We computationally identify sub-clusters representing most regenerative fin cell lineages, and define markers of osteoblasts, intra- and inter-ray fibroblasts and growth-promoting distal blastema cells. A pseudotemporal trajectory and in vivo photoconvertible lineage tracing indicate distal blastemal mesenchyme restores both intra- and inter-ray fibroblasts. Gene expression profiles across this trajectory suggest elevated protein production in the blastemal mesenchyme state. O-propargyl-puromycin incorporation and small molecule inhibition identify insulin growth factor receptor (IGFR)/mechanistic target of rapamycin kinase (mTOR)-dependent elevated bulk translation in blastemal mesenchyme and differentiating osteoblasts. We test candidate cooperating differentiation factors identified from the osteoblast trajectory, finding IGFR/mTOR signaling expedites glucocorticoid-promoted osteoblast differentiation in vitro. Concordantly, mTOR inhibition slows but does not prevent fin regenerative outgrowth in vivo. IGFR/mTOR may elevate translation in both fibroblast- and osteoblast-lineage cells during the outgrowth phase as a tempo-coordinating rheostat.

Cellular heterogeneity of the developing worker honey bee (Apis mellifera) pupa: a single cell transcriptomics analysis

It is estimated that animals pollinate 87.5% of flowering plants worldwide and that managed honey bees (Apis mellifera) account for 30-50% of this ecosystem service to agriculture. In addition to their important role as pollinators, honey bees are well-established insect models for studying learning and memory, behavior, caste differentiation, epigenetic mechanisms, olfactory biology, sex determination, and eusociality. Despite their importance to agriculture, knowledge of honey bee biology lags behind many other livestock species. In this study, we have used scRNA-Seq to map cell types to different developmental stages of the worker honey bee (prepupa at day 11 and pupa at day 15) and sought to determine their gene expression signatures. To identify cell-type populations, we examined the cell-to-cell network based on the similarity of the single-cells transcriptomic profiles. Grouping similar cells together we identified 63 different cell clusters of which 17 clusters were identifiable at both stages. To determine genes associated with specific cell populations or with a particular biological process involved in honey bee development, we used gene coexpression analysis. We combined this analysis with literature mining, the honey bee protein atlas, and gene ontology analysis to determine cell cluster identity. Of the cell clusters identified, 17 were related to the nervous system and sensory organs, 7 to the fat body, 19 to the cuticle, 5 to muscle, 4 to compound eye, 2 to midgut, 2 to hemocytes, and 1 to malpighian tubule/pericardial nephrocyte. To our knowledge, this is the first whole single-cell atlas of honey bees at any stage of development and demonstrates the potential for further work to investigate their biology at the cellular level.