A single-cell transcriptome atlas for zebrafish development

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

Zebrafish (Danio rerio) as a translational model for neuro-immune interactions in the enteric nervous system in autism spectrum disorders

Autism spectrum disorders (ASD) affect about 1% of the population and are strongly associated with gastrointestinal diseases creating shortcomings in quality of life. Multiple factors contribute to the development of ASD and although neurodevelopmental deficits are central, the pathogenesis of the condition is complex and the high prevalence of intestinal disorders is poorly understood. In agreement with the prominent research establishing clear bidirectional interactions between the gut and the brain, several studies have made it evident that such a relation also exists in ASD. Thus, dysregulation of the gut microbiota and gut barrier integrity may play an important role in ASD. However, only limited research has investigated how the enteric nervous system (ENS) and intestinal mucosal immune factors may impact on the development of ASD-related intestinal disorders. This review focuses on the mechanistic studies that elucidate the regulation and interactions between enteric immune cells, residing gut microbiota and the ENS in models of ASD. Especially the multifaceted properties and applicability of zebrafish (Danio rerio) for the study of ASD pathogenesis are assessed in comparison to studies conducted in rodent models and humans. Advances in molecular techniques and in vivo imaging, combined with genetic manipulation and generation of germ-free animals in a controlled environment, appear to make zebrafish an underestimated model of choice for the study of ASD. Finally, we establish the research gaps that remain to be explored to further our understanding of the complexity of ASD pathogenesis and associated mechanisms that may lead to intestinal disorders.

Regulative development as a model for origin of life and artificial life studies

Using the formal framework of the Free Energy Principle, we show how generic thermodynamic requirements on bidirectional information exchange between a system and its environment can generate complexity. This leads to the emergence of hierarchical computational architectures in systems that operate sufficiently far from thermal equilibrium. In this setting, the environment of any system increases its ability to predict system behavior by "engineering" the system towards increased morphological complexity and hence larger-scale, more macroscopic behaviors. When seen in this light, regulative development becomes an environmentally-driven process in which "parts" are assembled to produce a system with predictable behavior. We suggest on this basis that life is thermodynamically favorable and that, when designing artificial living systems, human engineers are acting like a generic "environment".

"Lighting up" fluoride: cellular imaging and zebrafish model interrogations using a simple ESIPT-based mycophenolic acid precursor-based probe

The discovery and implementation of media that derive from bioinspired designs and bear optical readouts featuring large Stokes shifts are of continued interest to a wide variety of researchers and clinicians. Myco-F, a novel mycophenolic acid precursor-based probe features a cleavable tert-butyldimethylsiloxy group to allow for fluoride detection. Myco-F exhibits high selectivity and specificity towards F^- (Stokes shift = 120 nm). All measurements were performed in complete aqueous media (LOD=0.38 μM). Myco-F enables detection of fluoride ions in living HEK293 cells and localizes in the eye region (among other regions) of the zebrafish. DFT calculations support the proposed ESIPT working photomechanism.

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 to study vertebrate aging. Here we describe the first multi-tissue, single-cell gene expression atlas of female and male turquoise killifish tissues comprising immune and metabolic cells from the blood, kidney, liver, and spleen. We were able to annotate 22 distinct cell types, define associated marker genes, and infer differentiation trajectories. Using this dataset, we found pervasive sex-dimorphic gene expression across cell types, especially in the liver. Sex-dimorphic genes tended to be involved in processes related to lipid metabolism, and indeed, we observed clear differences in lipid storage in female vs. male turquoise killifish livers. Importantly, we use machine-learning to predict sex using single-cell gene expression in our atlas and identify potential transcriptional markers for molecular sex identity in this species. As proof-of-principle, we show that our atlas can be used to deconvolute existing liver bulk RNA-seq data in this species to obtain accurate estimates of cell type proportions across biological conditions. We believe that this single-cell atlas can be a resource to the community that could notably be leveraged to identify cell type-specific genes for cell type-specific expression in transgenic animals.

RhoA GEF Mcf2lb regulates rosette integrity during collective cell migration

During development, multicellular rosettes serve as important cellular intermediates in the formation of diverse organ systems. Multicellular rosettes are transient epithelial structures that are defined by the apical constriction of cells towards the rosette center. Due to the important role these structures play during development, understanding the molecular mechanisms by which rosettes are formed and maintained is of high interest. Utilizing the zebrafish posterior lateral line primordium (pLLP) as a model system, we identify the RhoA GEF Mcf2lb as a regulator of rosette integrity. 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. Given the known role of RhoA in rosette formation, we asked whether Mcf2lb plays a role in regulating apical constriction of cells within rosettes. Live imaging and subsequent 3D analysis of mcf2lb mutant pLLP cells showed disrupted apical constriction and subsequent rosette organization. This in turn resulted in a unique posterior Lateral Line phenotype: 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 normally polarized. In contrast, signaling components that mediate apical constriction downstream of RhoA, Rock-2a and non-muscle Myosin II were diminished apically. Altogether our results suggest a model whereby Mcf2lb activates RhoA, which in turn activates downstream signaling machinery to induce and maintain apical constriction in cells incorporated into rosettes.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

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 hours 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 suggest new long-term cycling populations. Focused analyses of non-skeletal muscle and the endoderm identified transcriptional profiles of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and homologs of recently discovered human best4+ enterocytes. The transcriptional regulators of these populations remain unknown, so we reconstructed gene expression trajectories to suggest candidates. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Loss of RREB1 in pancreatic beta cells reduces cellular insulin content and affects endocrine cell gene expression

AIMS/HYPOTHESIS

Genome-wide studies have uncovered multiple independent signals at the RREB1 locus associated with altered type 2 diabetes risk and related glycaemic traits. However, little is known about the function of the zinc finger transcription factor Ras-responsive element binding protein 1 (RREB1) in glucose homeostasis or how changes in its expression and/or function influence diabetes risk.

METHODS

A zebrafish model lacking rreb1a and rreb1b was used to study the effect of RREB1 loss in vivo. Using transcriptomic and cellular phenotyping of a human beta cell model (EndoC-βH1) and human induced pluripotent stem cell (hiPSC)-derived beta-like cells, we investigated how loss of RREB1 expression and activity affects pancreatic endocrine cell development and function. Ex vivo measurements of human islet function were performed in donor islets from carriers of RREB1 type 2 diabetes risk alleles.

RESULTS

CRISPR/Cas9-mediated loss of rreb1a and rreb1b function in zebrafish supports an in vivo role for the transcription factor in beta cell mass, beta cell insulin expression and glucose levels. Loss of RREB1 also reduced insulin gene expression and cellular insulin content in EndoC-βH1 cells and impaired insulin secretion under prolonged stimulation. Transcriptomic analysis of RREB1 knockdown and knockout EndoC-βH1 cells supports RREB1 as a novel regulator of genes involved in insulin secretion. In vitro differentiation of RREB1KO/KO hiPSCs revealed dysregulation of pro-endocrine cell genes, including RFX family members, suggesting that RREB1 also regulates genes involved in endocrine cell development. Human donor islets from carriers of type 2 diabetes risk alleles in RREB1 have altered glucose-stimulated insulin secretion ex vivo, consistent with a role for RREB1 in regulating islet cell function.

CONCLUSIONS/INTERPRETATION

Together, our results indicate that RREB1 regulates beta cell function by transcriptionally regulating the expression of genes involved in beta cell development and function.

Zebrafish as outgroup model to study evolution of scavenger receptor class B type I functions

BACKGROUND AND AIMS

Scavenger receptor class B1 (SCARB1) - also known as the high-density lipoprotein (HDL) receptor - is a multi-ligand scavenger receptor that is primarily expressed in liver and steroidogenic organs. This receptor is known for its function in reverse cholesterol transport (RCT) in mammals and hence disruption leads to a massive increase in HDL cholesterol in these species. The extracellular domain of SCARB1 - which is important for cholesterol handling - is highly conserved across multiple vertebrates, except in zebrafish.

METHODS

To examine the functional conservation of SCARB1 among vertebrates, two stable scarb1 knockout zebrafish lines, scarb1 715delA (scarb1 -1 nt) and scarb1 715_716insGG (scarb1 +2 nt), were created using CRISPR-Cas9 technology.

RESULTS

We demonstrate that, in zebrafish, SCARB1 deficiency leads to disruption of carotenoid-based pigmentation, reduced fertility, and a decreased larvae survival rate, whereas steroidogenesis was unaltered. The observed reduced fertility is driven by defects in female fertility (-50 %, p < 0.001). Importantly, these alterations were independent of changes in free (wild-type 2.4 ± 0.2 μg/μl versus scarb1^-/- 2.0 ± 0.1 μg/μl) as well as total (wild-type 4.2 ± 0.4 μg/μl versus scarb1^-/- 4.0 ± 0.3 μg/μl) plasma cholesterol levels. Uptake of HDL in the liver of scarb1^-/- zebrafish larvae was reduced (-86.7 %, p < 0.001), but this coincided with reduced perfusion of the liver. No effect was observed on lipoprotein uptake in the caudal vein. SCARB1 deficient canaries, which also lack carotenoids in their plumage, similarly as scarb1^-/- zebrafish, failed to show an increase in plasma free- and total cholesterol levels.

CONCLUSION

Our findings suggest that the specific function of SCARB1 in maintaining plasma cholesterol could be an evolutionary novelty that became prominent in mammals, while other known functions were already present earlier during vertebrate evolution.

Neuronal birthdate reveals topography in a vestibular brainstem circuit for gaze stabilization

Across the nervous system, neurons with similar attributes are topographically organized. This topography reflects developmental pressures. Oddly, vestibular (balance) nuclei are thought to be disorganized. By measuring activity in birthdated neurons, we revealed a functional map within the central vestibular projection nucleus that stabilizes gaze in the larval zebrafish. We first discovered that both somatic position and stimulus selectivity follow projection neuron birthdate. Next, with electron microscopy and loss-of-function assays, we found that patterns of peripheral innervation to projection neurons were similarly organized by birthdate. Finally, birthdate revealed spatial patterns of axonal arborization and synapse formation to projection neuron outputs. Collectively, we find that development reveals previously hidden organization to the input, processing, and output layers of a highly conserved vertebrate sensorimotor circuit. The spatial and temporal attributes we uncover constrain the developmental mechanisms that may specify the fate, function, and organization of vestibulo-ocular reflex neurons. More broadly, our data suggest that, like invertebrates, temporal mechanisms may assemble vertebrate sensorimotor architecture.

Zebrafish Galectin 3 binding protein is the target antigen of the microglial 4C4 monoclonal antibody

BACKGROUND

Two decades ago, the fish-specific monoclonal antibody 4C4 was found to be highly reactive to zebrafish microglia, the macrophages of the central nervous system. This has resulted in 4C4 being widely used, in combination with available fluorescent transgenic reporters to identify and isolate microglia. However, the target protein of 4C4 remains unidentified, which represents a major caveat. In addition, whether the 4C4 expression pattern is strictly restricted to microglial cells in zebrafish has never been investigated.

RESULTS

Having demonstrated that 4C4 is able to capture its native antigen from adult brain lysates, we used immunoprecipitation/mass-spectrometry, coupled to recombinant expression analyses, to identify its target. The cognate antigen was found to be a paralog of Galectin 3 binding protein (Lgals3bpb), known as MAC2-binding protein in mammals. Notably, 4C4 did not recognize other paralogs, demonstrating specificity. Moreover, our data show that Lgals3bpb expression, while ubiquitous in microglia, also identifies leukocytes in the periphery, including populations of gut and liver macrophages.

CONCLUSIONS

The 4C4 monoclonal antibody recognizes Lgals3bpb, a predicted highly glycosylated protein whose function in the microglial lineage is currently unknown. Identification of Lgals3bpb as a new pan-microglia marker will be fundamental in forthcoming studies using the zebrafish model.

Activated PI3K delta syndrome 1 mutations cause neutrophilia in zebrafish larvae

People with activated PI3 kinase delta syndrome 1 (APDS1) suffer from immune deficiency and severe bronchiectasis. APDS1 is caused by dominant activating mutations of the PIK3CD gene that encodes the PI3 kinase delta (PI3Kδ) catalytic subunit. Despite the importance of innate immunity defects in bronchiectasis, there has been limited investigation of neutrophils or macrophages in APDS1 patients or mouse models. Zebrafish embryos provide an ideal system to study neutrophils and macrophages. We used CRISPR-Cas9 and CRISPR-Cpf1, with oligonucleotide-directed homologous repair, to engineer zebrafish equivalents of the two most prevalent human APDS1 disease mutations. These zebrafish pik3cd alleles dominantly caused excessive neutrophilic inflammation in a tail-fin injury model. They also resulted in total body neutrophilia in the absence of any inflammatory stimulus but normal numbers of macrophages. Exposure of zebrafish to the PI3Kδ inhibitor CAL-101 reversed the total body neutrophilia. There was no apparent defect in neutrophil maturation or migration, and tail-fin regeneration was unimpaired. Overall, the finding is of enhanced granulopoeisis, in the absence of notable phenotypic change in neutrophils and macrophages.

Cell-type-specific responses to the microbiota across all tissues of the larval zebrafish

Animal development proceeds in the presence of intimate microbial associations, but the extent to which different host cells across the body respond to resident microbes remains to be fully explored. Using the vertebrate model organism, the larval zebrafish, we assessed transcriptional responses to the microbiota across the entire body at single-cell resolution. We find that cell types across the body, not limited to tissues at host-microbe interfaces, respond to the microbiota. Responses are cell-type-specific, but across many tissues the microbiota enhances cell proliferation, increases metabolism, and stimulates a diversity of cellular activities, revealing roles for the microbiota in promoting developmental plasticity. This work provides a resource for exploring transcriptional responses to the microbiota across all cell types of the vertebrate body and generating new hypotheses about the interactions between vertebrate hosts and their microbiota.

Post-GWAS screening of candidate genes for refractive error in mutant zebrafish models

Genome-wide association studies (GWAS) have dissected numerous genetic factors underlying refractive errors (RE) such as myopia. Despite significant insights into understanding the genetic architecture of RE, few studies have validated and explored the functional role of candidate genes within these loci. To functionally follow-up on GWAS and characterize the potential role of candidate genes on the development of RE, we prioritized nine genes (TJP2, PDE11A, SHISA6, LAMA2, LRRC4C, KCNQ5, GNB3, RBFOX1, and GRIA4) based on biological and statistical evidence; and used CRISPR/cas9 to generate knock-out zebrafish mutants. These mutant fish were screened for abnormalities in axial length by spectral-domain optical coherence tomography and refractive status by eccentric photorefraction at the juvenile (2 months) and adult (4 months) developmental stage. We found a significantly increased axial length and myopic shift in refractive status in three of our studied mutants, indicating a potential involvement of the human orthologs (LAMA2, LRRC4C, and KCNQ5) in myopia development. Further, in-situ hybridization studies showed that all three genes are expressed throughout the zebrafish retina. Our zebrafish models provide evidence of a functional role of these three genes in refractive error development and offer opportunities to elucidate pathways driving the retina-to-sclera signaling cascade that leads to myopia.

Differentiation and Maturation of Muscle and Fat Cells in Cultivated Seafood: Lessons from Developmental Biology

Cultivated meat, also known as cultured or cell-based meat, is meat produced directly from cultured animal cells rather than from a whole animal. Cultivated meat and seafood have been proposed as a means of mitigating the substantial harms associated with current production methods, including damage to the environment, antibiotic resistance, food security challenges, poor animal welfare, and-in the case of seafood-overfishing and ecological damage associated with fishing and aquaculture. Because biomedical tissue engineering research, from which cultivated meat draws a great deal of inspiration, has thus far been conducted almost exclusively in mammals, cultivated seafood suffers from a lack of established protocols for producing complex tissues in vitro. At the same time, fish such as the zebrafish Danio rerio have been widely used as model organisms in developmental biology. Therefore, many of the mechanisms and signaling pathways involved in the formation of muscle, fat, and other relevant tissue are relatively well understood for this species. The same processes are understood to a lesser degree in aquatic invertebrates. This review discusses the differentiation and maturation of meat-relevant cell types in aquatic species and makes recommendations for future research aimed at recapitulating these processes to produce cultivated fish and shellfish.

Oxoeicosanoid signaling mediates early antimicrobial defense in zebrafish

5-oxoETE is a bioactive lipid derived from arachidonic acid generated when phospholipase A₂ activation coincides with oxidative stress. Through its G protein-coupled receptor OXER1, pure 5-oxoETE is a potent leukocyte chemoattractant. Yet, its physiological function has remained elusive owing to the unusual OXER1 conservation pattern. OXER1 is conserved from fish to primates but not in rodents, precluding genetic loss-of-function studies in mouse. To determine its physiological role, we combine transcriptomic, lipidomic, and intravital imaging assays with genetic perturbations of the OXER1 ortholog hcar1-4 in zebrafish. Pseudomonas aeruginosa infection induces the synthesis of 5-oxoETE and its receptor, along with other inflammatory pathways. Hcar1-4 deletion attenuates neutrophil recruitment and decreases post-infection survival, which could be rescued by ectopic expression of hcar1-4 or human OXER1. By revealing 5-oxoETE as dominant lipid regulator of the early antimicrobial response in a non-rodent vertebrate, our work expands the current, rodent-centric view of early inflammation.

Stress resilience is established during development and is regulated by complement factors

Individuals in a population respond differently to stressful situations. While resilient individuals recover efficiently, others are susceptible to the same stressors. However, it remains challenging to determine if resilience is established as a trait during development or acquired later in life. Using a behavioral paradigm in zebrafish larvae, we show that resilience is a stable and heritable trait, which is determined and exhibited early in life. Resilient larvae show unique stress-induced transcriptional response, and larvae with mutations in resilience-associated genes, such as neuropeptide Y and miR218, are less resilient. Transcriptome analysis shows that resilient larvae downregulate multiple factors of the innate immune complement cascade in response to stress. Perturbation of critical complement factors leads to an increase in resilience. We conclude that resilience is established as a stable trait early during development and that neuropeptides and the complement pathway play positive and negative roles in determining resilience, respectively.

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 have 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 mature hindbrain consists of alternating odd- versus even-type rhombomeres, our scMultiome analyses do not detect extensive odd versus even 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 and 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 mature rhombomeres requires resolution of mixed transcription and chromatin states.

Construction of a cross-species cell landscape at single-cell level

Individual cells are basic units of life. Despite extensive efforts to characterize the cellular heterogeneity of different organisms, cross-species comparisons of landscape dynamics have not been achieved. Here, we applied single-cell RNA sequencing (scRNA-seq) to map organism-level cell landscapes at multiple life stages for mice, zebrafish and Drosophila. By integrating the comprehensive dataset of > 2.6 million single cells, we constructed a cross-species cell landscape and identified signatures and common pathways that changed throughout the life span. We identified structural inflammation and mitochondrial dysfunction as the most common hallmarks of organism aging, and found that pharmacological activation of mitochondrial metabolism alleviated aging phenotypes in mice. The cross-species cell landscape with other published datasets were stored in an integrated online portal-Cell Landscape. Our work provides a valuable resource for studying lineage development, maturation and aging.

Analyses of binding partners and functional domains for the developmentally essential protein Hmx3a/HMX3

HMX3 is a homeodomain protein with essential roles in CNS and ear development. Homeodomains are DNA-binding domains and hence homeodomain-containing proteins are usually assumed to be transcription factors. However, intriguingly, our recent data suggest that zebrafish Hmx3a may not require its homeodomain to function, raising the important question of what molecular interactions mediate its effects. To investigate this, we performed a yeast two-hybrid screen and identified 539 potential binding partners of mouse HMX3. Using co-immunoprecipitation, we tested whether a prioritized subset of these interactions are conserved in zebrafish and found that Tle3b, Azin1b, Prmt2, Hmgb1a, and Hmgn3 bind Hmx3a. Next, we tested whether these proteins bind the products of four distinct hmx3a mutant alleles that all lack the homeodomain. Embryos homozygous for two of these alleles develop abnormally and die, whereas zebrafish homozygous for the other two alleles are viable. We found that all four mutations abrogate binding to Prmt2 and Tle3b, whereas Azin1b binding was preserved in all cases. Interestingly, Hmgb1a and Hmgn3 had more affinity for products of the viable mutant alleles. These data shed light on how HMX3/Hmx3a might function at a molecular level and identify new targets for future study in these vital developmental processes.

Zebrafish Erc1b mediates motor innervation and organization of craniofacial muscles in control of jaw movement

BACKGROUND

Movement of the lower jaw, a common behavior observed among vertebrates, is required for eating and processing food. This movement is controlled by signals sent from the trigeminal motor nerve through neuromuscular junctions (NMJs) to the masticatory muscles. Dysfunctional jaw movements contribute to craniomandibular disorders, yet the pathophysiology of these disorders is not well understood, as limited studies have been conducted on the molecular mechanisms of jaw movement.

RESULTS

Using erc1b/kimm533 genetic loss of function mutant, we evaluated lower jaw muscle organization and innervation by the cranial motor nerves in developing zebrafish. Using time-lapse confocal imaging of the erc1b mutant in a transgenic fluorescent reporter line, we found delayed trigeminal nerve growth and disrupted nerve branching architecture during muscle innervation. By automated 3D image analysis of NMJ distribution, we identified an increased number of small, disorganized NMJ clusters in erc1b mutant larvae compared to WT siblings. Using genetic replacement experiments, we determined the Rab GTPase binding domain of Erc1b is required for cranial motor nerve branching, but not NMJ organization or muscle attachment.

CONCLUSIONS

We identified Erc1b/ERC1 as a novel component of a genetic pathway contributing to muscle organization, trigeminal nerve outgrowth, and NMJ spatial distribution during development that is required for jaw movement.

Ancient multiplicity in cyclic nucleotide-gated (CNG) cation channel repertoire was reduced in the ancestor of Olfactores before re-expansion by whole genome duplications in vertebrates

Cyclic nucleotide-gated (CNG) cation channels are important heterotetrameric proteins in the retina, with different subunit composition in cone and rod photoreceptor cells: three CNGA3 and one CNGB3 in cones and three CNGA1 and one CNGB1 in rods. CNGA and CNGB subunits form separate subfamilies. We have analyzed the evolution of the CNG gene family in metazoans, with special focus on vertebrates by using sequence-based phylogeny and conservation of chromosomal synteny to deduce paralogons resulting from the early vertebrate whole genome duplications (WGDs). Our analyses show, unexpectedly, that the CNGA subfamily had four sister subfamilies in the ancestor of bilaterians and cnidarians that we named CNGC, CNGD, CNGE and CNGF. Of these, CNGC, CNGE and CNGF were lost in the ancestor of Olfactores while CNGD was lost in the vertebrate ancestor. The remaining CNGA and CNGB genes were expanded by a local duplication of CNGA and the subsequent chromosome duplications in the basal vertebrate WGD events. Upon some losses, this resulted in the gnathostome ancestor having three members in the visual CNGA subfamily (CNGA1-3), a single CNGA4 gene, and two members in the CNGB subfamily (CNGB1 and CNGB3). The nature of chromosomal rearrangements in the vertebrate CNGA paralogon was resolved by including the genomes of a non-teleost actinopterygian and an elasmobranch. After the teleost-specific WGD, additional duplicates were generated and retained for CNGA1, CNGA2, CNGA3 and CNGB1. Furthermore, teleosts retain a local duplicate of CNGB3. The retention of duplicated CNG genes is explained by their subfunctionalisation and photoreceptor-specific expression. In conclusion, this study provides evidence for four previously unknown CNG subfamilies in metazoans and further evidence that the early vertebrate WGD events were instrumental in the evolution of the vertebrate visual and central nervous systems.

Urotensin II-related peptides, Urp1 and Urp2, control zebrafish spine morphology

The spine provides structure and support to the body, yet how it develops its characteristic morphology as the organism grows is little understood. This is underscored by the commonality of conditions in which the spine curves abnormally such as scoliosis, kyphosis, and lordosis. Understanding the origin of these spinal curves has been challenging in part due to the lack of appropriate animal models. Recently, zebrafish have emerged as promising tools with which to understand the origin of spinal curves. Using zebrafish, we demonstrate that the urotensin II-related peptides (URPs), Urp1 and Urp2, are essential for maintaining spine morphology. Urp1 and Urp2 are 10-amino acid cyclic peptides expressed by neurons lining the central canal of the spinal cord. Upon combined genetic loss of Urp1 and Urp2, adolescent-onset planar curves manifested in the caudal region of the spine. Highly similar curves were caused by mutation of Uts2r3, an URP receptor. Quantitative comparisons revealed that urotensin-associated curves were distinct from other zebrafish spinal curve mutants in curve position and direction. Last, we found that the Reissner fiber, a proteinaceous thread that sits in the central canal and has been implicated in the control of spine morphology, breaks down prior to curve formation in mutants with perturbed cilia motility but was unaffected by loss of Uts2r3. This suggests a Reissner fiber-independent mechanism of curvature in urotensin-deficient mutants. Overall, our results show that Urp1 and Urp2 control zebrafish spine morphology and establish new animal models of spine deformity.

Knockdown of myorg leads to brain calcification in zebrafish

Primary familial brain calcification (PFBC) is a neurogenetic disorder characterized by bilateral calcified deposits in the brain. We previously identified that MYORG as the first pathogenic gene for autosomal recessive PFBC, and established a Myorg-KO mouse model. However, Myorg-KO mice developed brain calcifications until nine months of age, which limits their utility as a facile PFBC model system. Hence, whether there is another typical animal model for mimicking PFBC phenotypes in an early stage still remained unknown. In this study, we profiled the mRNA expression pattern of myorg in zebrafish, and used a morpholino-mediated blocking strategy to knockdown myorg mRNA at splicing and translation initiation levels. We observed multiple calcifications throughout the brain by calcein staining at 2–4 days post-fertilization in myorg-deficient zebrafish, and rescued the calcification phenotype by replenishing myorg cDNA. Overall, we built a novel model for PFBC via knockdown of myorg by antisense oligonucleotides in zebrafish, which could shorten the observation period and replenish the Myorg-KO mouse model phenotype in mechanistic and therapeutic studies.

The CXCR4-CXCL12 axis promotes T cell reconstitution via efficient hematopoietic immigration

T cells play a critical role in immunity to protect against pathogens and malignant cells. T cell immunodeficiency is detrimental, especially when T cell perturbation occurs during severe infection, irradiation, chemotherapy, and age-related thymic atrophy. Therefore, strategies that enhance T cell reconstitution provide considerable benefit and warrant intensive investigation. Here, we report the construction of a T cell ablation model in Tg(coro1a:DenNTR) zebrafish via metronidazole administration. The nascent T cells are mainly derived from the hematopoietic cells migrated from the kidney, the functional homolog of bone marrow and the complete recovery time is 6.5 days post-treatment. The cxcr4b gene is upregulated in the responsive hematopoietic cells. Functional interference of CXCR4 via both genetic and chemical manipulations does not greatly affect T lymphopoiesis, but delays T cell regeneration by disrupting hematopoietic migration. In contrast, cxcr4b accelerates the replenishment of hematopoietic cells in the thymus. Consistently, Cxcl12b, a ligand of Cxcr4, is increased in the thymic epithelial cells of the injured animals. Decreased or increased expression of Cxcl12b results in compromised or accelerated T cell recovery, respectively, similar to those observed with Cxcr4b. Taken together, our study reveals a role of CXCR4-CXCL12 signaling in promoting T cell recovery and provides a promising target for the treatment of immunodeficiency due to T cell injury.

Genetic duplication of tissue factor reveals subfunctionalization in venous and arterial hemostasis

Tissue factor (TF) is an evolutionarily conserved protein necessary for initiation of hemostasis. Zebrafish have two copies of the tissue factor gene (f3a and f3b) as the result of an ancestral teleost fish duplication event (so called ohnologs). In vivo physiologic studies of TF function have been difficult given early lethality of TF knockout in the mouse. We used genome editing to produce knockouts of both f3a and f3b in zebrafish. Since ohnologs arose through sub- or neofunctionalization, they can unmask unknown functions of non-teleost genes and could reveal whether mammalian TF has developmental functions distinct from coagulation. Here we show that a single copy of either f3a or f3b is necessary and sufficient for normal lifespan. Complete loss of TF results in lethal hemorrhage by 2-4 months despite normal embryonic and vascular development. Larval vascular endothelial injury reveals predominant roles for TFa in venous circulation and TFb in arterial circulation. Finally, we demonstrate that loss of TF predisposes to a stress-induced cardiac tamponade independent of its role in fibrin formation. Overall, our data suggest partial subfunctionalization of TFa and TFb. This multigenic zebrafish model has the potential to facilitate study of the role of TF in different vascular beds.

wnt16 regulates spine and muscle morphogenesis through parallel signals from notochord and dermomyotome

Bone and muscle are coupled through developmental, mechanical, paracrine, and autocrine signals. Genetic variants at the CPED1-WNT16 locus are dually associated with bone- and muscle-related traits. While Wnt16 is necessary for bone mass and strength, this fails to explain pleiotropy at this locus. Here, we show wnt16 is required for spine and muscle morphogenesis in zebrafish. In embryos, wnt16 is expressed in dermomyotome and developing notochord, and contributes to larval myotome morphology and notochord elongation. Later, wnt16 is expressed at the ventral midline of the notochord sheath, and contributes to spine mineralization and osteoblast recruitment. Morphological changes in wnt16 mutant larvae are mirrored in adults, indicating that wnt16 impacts bone and muscle morphology throughout the lifespan. Finally, we show that wnt16 is a gene of major effect on lean mass at the CPED1-WNT16 locus. Our findings indicate that Wnt16 is secreted in structures adjacent to developing bone (notochord) and muscle (dermomyotome) where it affects the morphogenesis of each tissue, thereby rendering wnt16 expression into dual effects on bone and muscle morphology. This work expands our understanding of wnt16 in musculoskeletal development and supports the potential for variants to act through WNT16 to influence bone and muscle via parallel morphogenetic processes.

In humans, genetic variants (DNA sequences that vary amongst individuals) have been identified that appear to influence two tissues, bone and skeletal muscle. However, how single genes and genetic variants exert dual influence on both tissues is not well understood. In this study, we found that the wnt16 gene is necessary for specifying the size and shape of both muscle and bone during development in zebrafish. We also disentangled how wnt16 affects both tissues: distinct cellular populations adjacent to muscle and bone secrete Wnt16, where it acts as a signal guiding the size and shape of each tissue. This is important because in humans, genetic variants near the WNT16 gene have effects on both bone- and muscle-related traits. This study expands our understanding of the role of WNT16 in bone and muscle development, and helps to explain how genetic variants near WNT16 affect traits for both tissues. Moreover, WNT16 is actively being explored as a target for osteoporosis therapies, thus our study could have implications with regard to the potential of targeting WNT16 to treat bone and muscle simultaneously.

Mycn regulates intestinal development through ribosomal biogenesis in a zebrafish model of Feingold syndrome 1

Feingold syndrome type 1, caused by loss-of-function of MYCN, is characterized by varied phenotypes including esophageal and duodenal atresia. However, no adequate model exists for studying the syndrome's pathological or molecular mechanisms, nor is there a treatment strategy. Here, we developed a zebrafish Feingold syndrome type 1 model with nonfunctional mycn, which had severe intestinal atresia. Single-cell RNA-seq identified a subcluster of intestinal cells that were highly sensitive to Mycn, and impaired cell proliferation decreased the overall number of intestinal cells in the mycn mutant fish. Bulk RNA-seq and metabolomic analysis showed that expression of ribosomal genes was down-regulated and that amino acid metabolism was abnormal. Northern blot and ribosomal profiling analysis showed abnormal rRNA processing and decreases in free 40S, 60S, and 80S ribosome particles, which led to impaired translation in the mutant. Besides, both Ribo-seq and western blot analysis showed that mTOR pathway was impaired in mycn mutant, and blocking mTOR pathway by rapamycin treatment can mimic the intestinal defect, and both L-leucine and Rheb, which can elevate translation via activating TOR pathway, could rescue the intestinal phenotype of mycn mutant. In summary, by this zebrafish Feingold syndrome type 1 model, we found that disturbance of ribosomal biogenesis and blockage of protein synthesis during development are primary causes of the intestinal defect in Feingold syndrome type 1. Importantly, our work suggests that leucine supplementation may be a feasible and easy treatment option for this disease.

In vivo tracking of KCC2b expression during early brain development

The neuronal chloride (Cl-) exporter, KCC2, regulates neuron excitability and development and undergoes a stereotypical pattern of delayed upregulation as neurons mature. KCC2 upregulation favors neural inhibition by establishing a negative Cl- gradient, ensuring GABA-induced Cl- currents are inward and inhibitory. We developed a zebrafish fluorescent reporter line, KCC2b:mCitrine, to track KCC2 expression in vivo during early brain development. KCC2b:mCitrine was first detected at 16 h postfertilization and by day 6 labeled most central and peripheral neurons and processes. At 20 h, expression was greatest in the soma-dense basal neuroepithelium but largely absent in apical and mantle zones where differentiation and migration primarily occur, and time lapse imaging at this stage supports a postmigration upregulation of KCC2b. Central dopamine neurons showed low KCC2b expression as observed in other species. KCC2b:mCitrine fluorescence was stable over minutes in most neurons, but brightness transients observed in single cells fit our expectation for real-time tracking of KCC2b upregulation in new neurons. To further assess whether fluorescence brightness tracks KCC2b expression, zebrafish embryos were exposed to bisphenol-A (BPA), which is known to suppress KCC2 expression. Fluorescence decreased after 6 days of BPA exposure but not after 2 or 4 days, suggesting that it is an accurate but delayed indicator of KCC2b expression. KCC2b:mCitrine zebrafish present a new method for visualizing KCC2b's complex dynamics during brain development, and potentially screening compounds aimed at modulating KCC2 expression.

A conserved transcriptional fingerprint of multi-neurotransmitter neurons necessary for social behavior

Background

An essential determinant of a neuron’s functionality is its neurotransmitter phenotype. We previously identified a defined subpopulation of cholinergic neurons required for social orienting behavior in zebrafish.

Results

We transcriptionally profiled these neurons and discovered that they are capable of synthesizing both acetylcholine and GABA. We also established a constellation of transcription factors and neurotransmitter markers that can be used as a “transcriptomic fingerprint” to recognize a homologous neuronal population in another vertebrate.

Conclusion

Our results suggest that this transcriptomic fingerprint and the cholinergic-GABAergic neuronal subtype that it defines are evolutionarily conserved.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08879-w.

Transcriptional profiling from whole embryos to single neuroblast lineages in Drosophila

Embryonic development results in the production of distinct tissue types, and different cell types within each tissue. A major goal of developmental biology is to uncover the "parts list" of cell types that comprise each organ. Here we perform single cell RNA sequencing (scRNA-seq) of the Drosophila embryo to identify the genes that characterize different cell and tissue types during development. We assay three different timepoints, revealing a coordinated change in gene expression within each tissue. Interestingly, we find that the elav and Mhc genes, whose protein products are widely used as markers for neurons and muscles, respectively, show broad pan-embryonic expression, indicating the importance of post-transcriptional regulation. We next focus on the central nervous system (CNS), where we identify genes whose expression is enriched at each stage of neuronal differentiation: from neural progenitors, called neuroblasts, to their immediate progeny ganglion mother cells (GMCs), followed by new-born neurons, young neurons, and the most mature neurons. Finally, we ask whether the clonal progeny of a single neuroblast (NB7-1) share a similar transcriptional identity. Surprisingly, we find that clonal identity does not lead to transcriptional clustering, showing that neurons within a lineage are diverse, and that neurons with a similar transcriptional profile (e.g. motor neurons, glia) are distributed among multiple neuroblast lineages. Although each lineage consists of diverse progeny, we were able to identify a previously uncharacterized gene, Fer3, as an excellent marker for the NB7-1 lineage. Within the NB7-1 lineage, neurons which share a temporal identity (e.g. Hunchback, Kruppel, Pdm, and Castor temporal transcription factors in the NB7-1 lineage) have shared transcriptional features, allowing for the identification of candidate novel temporal factors or targets of the temporal transcription factors. In conclusion, we have characterized the embryonic transcriptome for all major tissue types and for three stages of development, as well as the first transcriptomic analysis of a single, identified neuroblast lineage, finding a lineage-enriched transcription factor.

Single cell RNA-seq analysis reveals temporally-regulated and quiescence-regulated gene expression in Drosophila larval neuroblasts

The mechanisms that generate neural diversity during development remains largely unknown. Here, we use scRNA-seq methodology to discover new features of the Drosophila larval CNS across several key developmental timepoints. We identify multiple progenitor subtypes - both stem cell-like neuroblasts and intermediate progenitors - that change gene expression across larval development, and report on new candidate markers for each class of progenitors. We identify a pool of quiescent neuroblasts in newly hatched larvae and show that they are transcriptionally primed to respond to the insulin signaling pathway to exit from quiescence, including relevant pathway components in the adjacent glial signaling cell type. We identify candidate "temporal transcription factors" (TTFs) that are expressed at different times in progenitor lineages. Our work identifies many cell type specific genes that are candidates for functional roles, and generates new insight into the differentiation trajectory of larval neurons.

dock8 deficiency attenuates microglia colonization in early zebrafish larvae

Microglia are tissue-resident macrophages that carry out immune functions in the brain. The deficiency or dysfunction of microglia has been implicated in many neurodegenerative disorders. DOCK8, a member of the DOCK family, functions as a guanine nucleotide exchange factor and plays key roles in immune regulation and neurological diseases. The functions of DOCK8 in microglia development are not fully understood. Here, we generated zebrafish dock8 mutants by CRISPR/Cas9 genome editing and showed that dock8 mutations attenuate microglia colonization in the zebrafish midbrain at early larvae stages. In vivo time-lapse imaging revealed that the motility of macrophages was reduced in the dock8 mutant. We further found that cdc42/cdc42l, which encode the small GTPase activated by Dock8, also regulate microglia colonization in zebrafish. Collectively, our study suggests that the Dock8-Cdc42 pathway is required for microglia colonization in zebrafish larvae.

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Single-cell transcriptomic analysis of vascular endothelial cells in zebrafish embryos

Vascular endothelial cells exhibit substantial phenotypic and transcriptional heterogeneity which is established during early embryogenesis. However, the molecular mechanisms involved in establishing endothelial cell diversity are still not well understood. Zebrafish has emerged as an advantageous model to study vascular development. Despite its importance, the single-cell transcriptomic profile of vascular endothelial cells during zebrafish development is still missing. To address this, we applied single-cell RNA-sequencing (scRNA-seq) of vascular endothelial cells isolated from zebrafish embryos at the 24 hpf stage. Six distinct clusters or subclusters related to vascular endothelial cells were identified which include arterial, two venous, cranial, endocardial and endothelial progenitor cell subtypes. Furthermore, we validated our findings by characterizing novel markers for arterial, venous, and endocardial cells. We experimentally confirmed the presence of two transcriptionally different venous cell subtypes, demonstrating heterogeneity among venous endothelial cells at this early developmental stage. This dataset will be a valuable resource for future functional characterization of vascular endothelial cells and interrogation of molecular mechanisms involved in the establishment of their heterogeneity and cell-fate decisions.

FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma

Simple Summary

Despite successful treatment of primary uveal melanoma (UM), metastases still occur in approximately 50% of the patients. Unfortunately, little is known about the mechanism behind metastasized UM. By reanalyzing publicly available single-cell RNA sequencing data of embryonic zebrafish larvae and validating the results with UM data, we have identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. The most significant finding is FOXD1, which is nearly exclusively expressed in high-risk UM and is associated with poor survival. FOXD1 is a novel gene which could be involved in the metastatic capability of UM. Elucidating its function and role in metastatic UM could help to understand and develop treatment for UM.

Abstract

Uveal melanoma (UM) is a deadly ocular malignancy, originating from uveal melanocytes. Although much is known regarding prognostication in UM, the exact mechanism of metastasis is mostly unknown. Metastatic tumor cells are known to express a more stem-like RNA profile which is seen often in cell-specific embryonic development to induce tumor progression. Here, we identified novel transcription regulators by reanalyzing publicly available single cell RNA sequencing experiments. We identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. Our most significant finding is FOXD1, as this gene is nearly exclusively expressed in high-risk UM and its expression is associated with a poor prognosis. Even within the BAP1-mutated UM, the expression of FOXD1 is correlated with poor survival. FOXD1 is a novel factor which could potentially be involved in the metastatic capacity of high-risk UM. Elucidating the function of FOXD1 in UM could provide insight into the malignant transformation of uveal melanocytes, especially in high-risk UM.

Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology

The vertebrate sodium-iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters (slc5a7); these clades are the NIS clade of sodium/iodide and lactate cotransporters (slc5a5, slc5a6, slc5a8, slc5a8, and slc5a12) and the SGLT clade of sodium/glucose cotransporters (slc5a1, slc5a2, slc5a3, slc5a4, slc5a10, and slc5a11). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single-cell transcriptomics (scRNA-seq) on zebrafish sex-genotyped gonads revealed that NIS clade-expressing cells included germ cells (slc5a5, slc5a6a, and slc5a6b) and gonadal soma cells (slc5a8l). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Fibronectin and Integrin α5 play overlapping and independent roles in regulating the development of pharyngeal endoderm and cartilage

Craniofacial skeletal elements are derived from cranial neural crest cells (CNCCs), which migrate along discrete paths and populate distinct pharyngeal arches, structures that are separated by the neighboring endodermal pouches (EPs). Interactions between the CNCCs and the endoderm are critical for proper craniofacial development. In zebrafish, integrin α5 (Itga5) functions in the endoderm to regulate formation of specifically the first EP (EP1) and the development of the hyoid cartilage. Here we show that fibronectin (Fn), a major component of the extracellular matrix (ECM), is also required for these developmental processes, and that the penetrance of defects in mutants is temperature-dependent. fn1a^-/- embryos exhibited defects that are similar to, but much more severe than, those of itga5^-/- embryos, and a loss of integrin av (itgav) function enhanced both endoderm and cartilage defects in itga5^-/- embryos, suggesting that Itga5 and Itgav cooperate to transmit signals from Fn to regulate the development of endoderm and cartilage. Whereas the endodermal defects in itga5; itga5v^-/- double mutant embryos were comparable to those of fn1a^-/- mutants, the cartilage defects were much milder. Furthermore, Fn assembly was detected in migrating CNCCs, and the epithelial organization and differentiation of CNCC-derived arches were impaired in fn1a^-/- embryos, indicating that Fn1 exerts functions in arch development that are independent of Itga5 and Itgav. Additionally, reduction of itga5 function in fn1a^-/- embryos led to profound defects in body axis elongation, as well as in endoderm and cartilage formation, suggesting that other ECM proteins signal through Itga5 to regulate development of the endoderm and cartilage. Thus, our studies reveal that Fn1a and Itga5 have both overlapping and independent functions in regulating development of the pharyngeal endoderm and cartilage.

Daw1 regulates the timely onset of cilia motility during development

Motile cilia generate cell propulsion and extracellular fluid flows that are crucial for airway clearance, fertility and left-right patterning. Motility is powered by dynein arm complexes that are assembled in the cytoplasm then imported into the cilium. Studies in Chlamydomonas reinhardtii showed that ODA16 is a cofactor which promotes dynein arm import. Here, we demonstrate that the zebrafish homolog of ODA16, Daw1, facilitates the onset of robust cilia motility during development. Without Daw1, cilia showed markedly reduced motility during early development; however, motility subsequently increased to attain close to wild-type levels. Delayed motility onset led to differential effects on early and late cilia-dependent processes. Remarkably, abnormal body axis curves, which formed during the first day of development due to reduced cilia motility, self-corrected when motility later reached wild-type levels. Zebrafish larva therefore possess the ability to survey and correct body shape abnormalities. This work defines Daw1 as a factor which promotes the onset of timely cilia motility and can explain why human patients harboring DAW1 mutations exhibit significant laterality perturbations but mild airway and fertility complications.

Summary: Daw1 promotes the onset of timely cilia motility for robust axial straightening during zebrafish development.

Weighted gene co-expression network analysis of embryos and first instar larvae of the horseshoe crab Tachypleus tridentatus uncovers development gene networks

Horseshoe crabs are marine chelicerates that have existed on Earth for about 450 million years, and they are often used as an experimental model for studying marine invertebrate embryology. In this study, we performed transcriptome gene expression profiling of four continuous embryonic stages (Stages 18-21) and first instar larvae of Tachypleus tridentatus. A mean of 50,742,995 high-quality clean reads was obtained from each library. We then conducted weighted gene co-expression network analysis (WGCNA) for 13,698 genes with fragments per kilobase of exon per million mapped fragments values >5. We identified 17 modules, six of which likely play critical roles in development. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differentially expressed genes was performed on the biologically significant modules. We found that several pathways, such as hedgehog signaling pathway, VEGF signaling pathway, dorso-ventral axis formation, may be involved in the embryonic development process of T. tridentatus. We also identified hub genes that were highly connected in the six critical modules. This is the first study to apply WGCNA to horseshoe crabs to identify hub genes that may play critical roles in development, and our results provide new insight into the mechanisms underlying early development in horseshoe crabs.

Expression Pattern of nos1 in the Developing Nervous System of Ray-Finned Fish

Fish have colonized nearly all aquatic niches, making them an invaluable resource to understand vertebrate adaptation and gene family evolution, including the evolution of complex neural networks and modulatory neurotransmitter pathways. Among ancient regulatory molecules, the gaseous messenger nitric oxide (NO) is involved in a wide range of biological processes. Because of its short half-life, the modulatory capability of NO is strictly related to the local activity of nitric oxide synthases (Nos), enzymes that synthesize NO from L-arginine, making the localization of Nos mRNAs a reliable indirect proxy for the location of NO action domains, targets, and effectors. Within the diversified actinopterygian nos paralogs, nos1 (alias nnos) is ubiquitously present as a single copy gene across the gnathostome lineage, making it an ideal candidate for comparative studies. To investigate variations in the NO system across ray-finned fish phylogeny, we compared nos1 expression patterns during the development of two well-established experimental teleosts (zebrafish and medaka) with an early branching holostean (spotted gar), an important evolutionary bridge between teleosts and tetrapods. Data reported here highlight both conserved expression domains and species-specific nos1 territories, confirming the ancestry of this signaling system and expanding the number of biological processes implicated in NO activities.

Aldh2 is a lineage-specific metabolic gatekeeper in melanocyte stem cells

Melanocyte stem cells (McSCs) in zebrafish serve as an on-demand source of melanocytes during growth and regeneration, but metabolic programs associated with their activation and regenerative processes are not well known. Here, using live imaging coupled with scRNA-sequencing, we discovered that, during regeneration, quiescent McSCs activate a dormant embryonic neural crest transcriptional program followed by an aldehyde dehydrogenase (Aldh) 2 metabolic switch to generate progeny. Unexpectedly, although ALDH2 is well known for its aldehyde-clearing mechanisms, we find that, in regenerating McSCs, Aldh2 activity is required to generate formate – the one-carbon (1C) building block for nucleotide biosynthesis – through formaldehyde metabolism. Consequently, we find that disrupting the 1C cycle with low doses of methotrexate causes melanocyte regeneration defects. In the absence of Aldh2, we find that purines are the metabolic end product sufficient for activated McSCs to generate progeny. Together, our work reveals McSCs undergo a two-step cell state transition during regeneration, and that the reaction products of Aldh2 enzymes have tissue-specific stem cell functions that meet metabolic demands in regeneration.

Summary: In zebrafish melanocyte regeneration, quiescent McSCs respond by re-expressing a neural crest identity, followed by an Aldh2-dependent metabolic switch to generate progeny.

Balancing the transcriptome: leveraging sample similarity to improve measures of gene specificity

The spatial and temporal domain of a gene's expression can range from ubiquitous to highly specific. Quantifying the degree to which this expression is unique to a specific tissue or developmental timepoint can provide insight into the etiology of genetic diseases. However, quantifying specificity remains challenging as measures of specificity are sensitive to similarity between samples in the sample set. For example, in the Gene-Tissue Expression project (GTEx), brain subregions are overrepresented at 13 of 54 (24%) unique tissues sampled. In this dataset, existing specificity measures have a decreased ability to identify genes specific to the brain relative to other organs. To solve this problem, we leverage sample similarity information to weight samples such that overrepresented tissues do not have an outsized effect on specificity estimates. We test this reweighting procedure on 4 measures of specificity, Z-score, Tau, Tsi and Gini, in the GTEx data and in single cell datasets for zebrafish and mouse. For all of these measures, incorporating sample similarity information to weight samples results in greater stability of sets of genes called as specific and decreases the overall variance in the change of specificity estimates as sample sets become more unbalanced. Furthermore, the genes with the largest improvement in their specificity estimate's stability are those with functions related to the overrepresented sample types. Our results demonstrate that incorporating similarity information improves specificity estimates' stability to the choice of the sample set used to define the transcriptome, providing more robust and reproducible measures of specificity for downstream analyses.

Connexinplexity: the spatial and temporal expression of connexin genes during vertebrate organogenesis

Animal development requires coordinated communication between cells. The Connexin family of proteins is a major contributor to intercellular communication in vertebrates by forming gap junction channels that facilitate the movement of ions, small molecules, and metabolites between cells. Additionally, individual hemichannels can provide a conduit to the extracellular space for paracrine and autocrine signaling. Connexin-mediated communication is widely used in epithelial, neural, and vascular development and homeostasis, and most tissues likely use this form of communication. In fact, Connexin disruptions are of major clinical significance contributing to disorders developing from all major germ layers. Despite the fact that Connexins serve as an essential mode of cellular communication, the temporal and cell-type-specific expression patterns of connexin genes remain unknown in vertebrates. A major challenge is the large and complex connexin gene family. To overcome this barrier, we determined the expression of all connexins in zebrafish using single-cell RNA-sequencing of entire animals across several stages of organogenesis. Our analysis of expression patterns has revealed that few connexins are broadly expressed, but rather, most are expressed in tissue- or cell-type-specific patterns. Additionally, most tissues possess a unique combinatorial signature of connexin expression with dynamic temporal changes across the organism, tissue, and cell. Our analysis has identified new patterns for well-known connexins and assigned spatial and temporal expression to genes with no-existing information. We provide a field guide relating zebrafish and human connexin genes as a critical step toward understanding how Connexins contribute to cellular communication and development throughout vertebrate organogenesis.

TFAP2 paralogs facilitate chromatin access for MITF at pigmentation and cell proliferation genes

In developing melanocytes and in melanoma cells, multiple paralogs of the Activating-enhancer-binding Protein 2 family of transcription factors (TFAP2) contribute to expression of genes encoding pigmentation regulators, but their interaction with Microphthalmia transcription factor (MITF), a master regulator of these cells, is unclear. Supporting the model that TFAP2 facilitates MITF's ability to activate expression of pigmentation genes, single-cell seq analysis of zebrafish embryos revealed that pigmentation genes are only expressed in the subset of mitfa-expressing cells that also express tfap2 paralogs. To test this model in SK-MEL-28 melanoma cells we deleted the two TFAP2 paralogs with highest expression, TFAP2A and TFAP2C, creating TFAP2 knockout (TFAP2-KO) cells. We then assessed gene expression, chromatin accessibility, binding of TFAP2A and of MITF, and the chromatin marks H3K27Ac and H3K27Me3 which are characteristic of active enhancers and silenced chromatin, respectively. Integrated analyses of these datasets indicate TFAP2 paralogs directly activate enhancers near genes enriched for roles in pigmentation and proliferation, and directly repress enhancers near genes enriched for roles in cell adhesion. Consistently, compared to WT cells, TFAP2-KO cells proliferate less and adhere to one another more. TFAP2 paralogs and MITF co-operatively activate a subset of enhancers, with the former necessary for MITF binding and chromatin accessibility. By contrast, TFAP2 paralogs and MITF do not appear to co-operatively inhibit enhancers. These studies reveal a mechanism by which TFAP2 profoundly influences the set of genes activated by MITF, and thereby the phenotype of pigment cells and melanoma cells.

Heat Shock Procedure Affects Cell Division-Associated Genes in Gynogenetic Manipulation

Heat shock procedure is crucial for gynogenetic manipulation leading to diploidization of the maternal genomes; however, the underlying molecular mechanism especially the transcriptomic changes during this procedure has still not been unveiled yet. Here, the artificial gynogenesis of zebrafish (Danio rerio) using inactivated sperm from rare minnow (Gobiocypris rarus) was conducted. We found that artificial gynogenetic manipulation, including pseudo-fertilization and heat shock, decreased hatching rates, whereas heat shock treatment alone had medium hatching rates. The first cleavage changed the expression of genes associated with RNA transcription and protein synthesis. A co-expression network regulated by hub genes GIT1, Sepsecs, and FLNB was significantly correlated with heat shock procedure. The cyclin family and cyclin-dependent kinase-related genes were lowly expressed in embryos from gynogenetic zebrafish, and genes involved in controlling the cell cycle and genomic stability were significantly altered by the gynogenetic treatment. Our results show the effects of artificial gynogenesis on embryos and describe changes in gene expression that suggest drastic changes take place in cell division by heat shock procedure. These findings will contribute to an understanding of the molecular basis for germplasm improving, including the purifying effect and allogynogenetic biological effect by gynogenesis.

Gain of gene regulatory network interconnectivity at the origin of vertebrates

In this manuscript, we address an essential question in developmental and evolutionary biology: How have changes in gene regulatory networks contributed to the invertebrate-to-vertebrate transition? To address this issue, we perturbed four signaling pathways critical for body plan formation in the cephalochordate amphioxus and in zebrafish and compared the effects of such perturbations on gene expression and gene regulation in both species. Our data reveal that many developmental genes have gained response to these signaling pathways in the vertebrate lineage. Moreover, we show that the interconnectivity between these pathways is much higher in zebrafish than in amphioxus. We conclude that this increased signaling pathway complexity likely contributed to vertebrate morphological novelties during evolution.

Signaling pathways control a large number of gene regulatory networks (GRNs) during animal development, acting as major tools for body plan formation [A. Pires-daSilva, R. J. Sommer, Nat. Rev. Genet. 4, 39–49 (2003)], although only a few of these pathways operate during this period [J. J. Sanz-Ezquerro, A. E. Münsterberg, S. Stricker, Front. Cell Dev. Biol. 5, 76 (2017)]. Moreover, most of them have been largely conserved during metazoan evolution [L. S. Babonis, M. Q. Martindale, Philos. Trans. R. Soc. Lond. B Biol. Sci. 372, 20150477 (2017)]. How evolution has generated a vast diversity of animal morphologies with such a limited number of tools is still largely unknown. Here, we show that gain of interconnectivity between signaling pathways and the GRNs they control may have critically contributed to the origin of vertebrates. We perturbed the retinoic acid, Wnt, FGF, and Nodal signaling pathways during gastrulation in the invertebrate chordate amphioxus and zebrafish and compared the effects on gene expression and cis-regulatory elements (CREs). We found that multiple developmental genes gain response to these pathways through vertebrate-specific CREs. Moreover, in contrast to amphioxus, many of these CREs responded to multiple pathways in zebrafish, which reflects their high interconnectivity. Furthermore, we found that vertebrate-specific cell types are more enriched in highly interconnected genes than in tissues with more ancient origin. Thus, the increase of CREs in vertebrates integrating inputs from different signaling pathways probably contributed to gene expression complexity and to the formation of new cell types and morphological novelties in this lineage.

Exploring Zebrafish Larvae as a COVID-19 Model: Probable Abortive SARS-CoV-2 Replication in the Swim Bladder

Animal models are essential to understanding COVID-19 pathophysiology and for preclinical assessment of drugs and other therapeutic or prophylactic interventions. We explored the small, cheap, and transparent zebrafish larva as a potential host for SARS-CoV-2. Bath exposure, as well as microinjection in the coelom, pericardium, brain ventricle, or bloodstream, resulted in a rapid decrease of SARS-CoV-2 RNA in wild-type larvae. However, when the virus was inoculated in the swim bladder, viral RNA stabilized after 24 h. By immunohistochemistry, epithelial cells containing SARS-CoV-2 nucleoprotein were observed in the swim bladder wall. Our data suggest an abortive infection of the swim bladder. In some animals, several variants of concern were also tested with no evidence of increased infectivity in our model. Low infectivity of SARS-CoV-2 in zebrafish larvae was not due to the host type I interferon response, as comparable viral loads were detected in type I interferon-deficient animals. A mosaic overexpression of human ACE2 was not sufficient to increase SARS-CoV-2 infectivity in zebrafish embryos or in fish cells in vitro. In conclusion, wild-type zebrafish larvae appear mostly non-permissive to SARS-CoV-2, except in the swim bladder, an aerial organ sharing similarities with the mammalian lung.