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

Increased Netrin downstream of overactive Hedgehog signaling disrupts optic fissure formation

BACKGROUND

Uveal coloboma, a developmental eye defect, is caused by failed development of the optic fissure, a ventral structure in the optic stalk and cup where axons exit the eye and vasculature enters. The Hedgehog (Hh) signaling pathway regulates optic fissure development: loss-of-function mutations in the Hh receptor ptch2 produce overactive Hh signaling and can result in coloboma. We previously proposed a model where overactive Hh signaling disrupts optic fissure formation by upregulating transcriptional targets acting both cell- and non-cell-autonomously. Here, we examine the Netrin family of secreted ligands as candidate Hh target genes.

RESULTS

We find multiple Netrin ligands upregulated in the zebrafish ptch2 mutant during optic fissure development. Using a gain-of-function approach to overexpress Netrin in a spatiotemporally specific manner, we find that netrin1a or netrin1b overexpression is sufficient to cause coloboma and disrupt wild-type optic fissure formation. We used loss-of-function alleles, CRISPR/Cas9 mutagenesis, and morpholino knockdown to test if loss of Netrin can rescue coloboma in the ptch2 mutant: loss of netrin genes does not rescue the ptch2 mutant phenotype.

CONCLUSION

These results suggest that Netrin is sufficient but not required to disrupt optic fissure formation downstream of overactive Hh signaling in the ptch2 mutant.

Cspg4 sculpts oligodendrocyte precursor cell morphology

The extracellular matrix (ECM) provides critical biochemical and structural cues that regulate neural development. Chondroitin sulfate proteoglycans (CSPGs), a major ECM component, have been implicated in modulating oligodendrocyte precursor cell (OPC) proliferation, migration, and maturation, but their specific roles in oligodendrocyte lineage cell (OLC) development and myelination in vivo remain poorly understood. Here, we use zebrafish as a model system to investigate the spatiotemporal dynamics of ECM deposition and CSPG localization during central nervous system (CNS) development, with a focus on their relationship to OLCs. We demonstrate that ECM components, including CSPGs, are dynamically expressed in distinct spatiotemporal patterns coinciding with OLC development and myelination. We found that zebrafish lacking cspg4 function produced normal numbers of OLCs, which appeared to undergo proper differentiation. However, OPC morphology in mutant larvae was aberrant. Nevertheless, the number and length of myelin sheaths produced by mature oligodendrocytes were unaffected. These data indicate that Cspg4 regulates OPC morphogenesis in vivo , supporting the role of the ECM in neural development.

Requirement of a novel gene, drish, in the zebrafish retinal ganglion cell and primary motor axon development

BACKGROUND

During neurogenesis, growing axons must navigate through the complex extracellular environment and make correct synaptic connections for the proper functioning of neural circuits. The mechanisms underlying the formation of functional neural networks are still only partially understood.

RESULTS

Here we analyzed the role of a novel gene si:ch73-364h19.1/drish in the neural and vascular development of zebrafish embryos. We show that drish mRNA is expressed broadly and dynamically in multiple cell types including neural, glial, retinal progenitor and vascular endothelial cells throughout the early stages of embryonic development. To study Drish function during embryogenesis, we generated drish genetic mutant using CRISPR/Cas9 genome editing. drish loss-of-function mutant larvae displayed defects in early retinal ganglion cell, optic nerve and the retinal inner nuclear layer formation, as well as ectopic motor axon branching. In addition, drish mutant adults exhibited deficient retinal outer nuclear layer and showed defective light response and locomotory behavior. However, vascular patterning and blood circulation were not significantly affected.

CONCLUSIONS

Together, these data demonstrate important roles of zebrafish drish in the retinal ganglion cell, optic nerve and interneuron development and in spinal motor axon branching.

Creation of a novel zebrafish model with low DHA status to study the role of maternal nutrition during neurodevelopment

Docosahexaenoic acid (DHA), a dietary omega-3 fatty acid, is a major building block of brain cell membranes. Offspring rely on maternal DHA transfer to meet their neurodevelopmental needs, but DHA sources are lacking in the American diet. Low DHA status is linked to altered immune responses, white matter defects, impaired vision, and an increased risk of psychiatric disorders during development. However, the underlying cellular mechanisms involved are largely unknown, and advancements in the field have been limited by the existing tools and animal models. Zebrafish are an excellent model for studying neurodevelopmental mechanisms. Embryos undergo rapid external development and are optically transparent, enabling direct observation of individual cells and dynamic cell-cell interactions in a way that is not possible in rodents. Here, we create a novel DHA-deficient zebrafish model by 1) disrupting elovl2, a key gene in the DHA biosynthesis pathway, via CRISPR-Cas9 genome editing, and 2) feeding mothers a DHA-deficient diet. We show that low DHA status during development is associated with a small eye morphological phenotype and demonstrate that even the morphologically normal siblings exhibit dysregulated gene pathways related to vision and stress response. Future work using our zebrafish model could reveal the cellular and molecular mechanisms by which low DHA status leads to neurodevelopmental abnormalities and provide insight into maternal nutritional strategies that optimize infant brain health.

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.

Tuning apicobasal polarity and junctional recycling in the hemogenic endothelium orchestrates the morphodynamic complexity of emerging pre-hematopoietic stem cells

Hematopoietic stem cells emerge in the embryo from an aortic-derived tissue called the hemogenic endothelium (HE). The HE appears to give birth to cells of different nature and fate but the molecular principles underlying this complexity are largely unknown. Here we show, in the zebrafish embryo, that two cell types emerge from the aortic floor with radically different morphodynamics. With the support of live imaging, we bring evidence suggesting that the mechanics underlying the two emergence types rely, or not, on apicobasal polarity establishment. While the first type is characterized by reinforcement of apicobasal polarity and maintenance of the apical/luminal membrane until release, the second type emerges via a dynamic process reminiscent of trans-endothelial migration. Interfering with Runx1 function suggests that the balance between the two emergence types depends on tuning apicobasal polarity at the level of the HE. In support of this and unexpectedly, we show that Pard3ba - one of the four Pard3 proteins expressed in the zebrafish - is sensitive to interference with Runx1 activity, in aortic endothelial cells. This supports the idea of a signaling cross talk controlling cell polarity and its associated features, between aortic and hemogenic cells. In addition, using new transgenic fish lines that express Junctional Adhesion Molecules and functional interference, we bring evidence for the essential role of ArhGEF11/PDZ-RhoGEF in controlling the HE-endothelial cell dynamic interface, including cell-cell intercalation, which is ultimately required for emergence completion. Overall, we highlight critical cellular and dynamic events of the endothelial-to-hematopoietic transition that support emergence complexity, with a potential impact on cell fate.

Kif1a and intact microtubules maintain synaptic-vesicle populations at ribbon synapses in zebrafish hair cells

Sensory hair cells of the inner ear utilize specialized ribbon synapses to transmit sensory stimuli to the central nervous system. This sensory transmission necessitates rapid and sustained neurotransmitter release, which relies on a large pool of synaptic vesicles at the hair-cell presynapse. Work in neurons has shown that kinesin motor proteins traffic synaptic material along microtubules to the presynapse, but how new synaptic material reaches the presynapse in hair cells is not known. We show that the kinesin motor protein Kif1a and an intact microtubule network are necessary to enrich synaptic vesicles at the presynapse in hair cells. We use genetics and pharmacology to disrupt Kif1a function and impair microtubule networks in hair cells of the zebrafish lateral-line system. We find that these manipulations decrease synaptic-vesicle populations at the presynapse in hair cells. Using electron microscopy, along with in vivo calcium imaging and electrophysiology, we show that a diminished supply of synaptic vesicles adversely affects ribbon-synapse function. Kif1a mutants exhibit dramatic reductions in spontaneous vesicle release and evoked postsynaptic calcium responses. Additionally, we find that kif1a mutants exhibit impaired rheotaxis, a behavior reliant on the ability of hair cells in the lateral line to respond to sustained flow stimuli. Overall, our results demonstrate that Kif1a-based microtubule transport is critical to enrich synaptic vesicles at the active zone in hair cells, a process that is vital for proper ribbon-synapse function.

Novel cell states arise in embryonic cells devoid of key reprogramming factors

The capacity for embryonic cells to differentiate relies on a large-scale reprogramming of the oocyte and sperm nucleus into a transient totipotent state. In zebrafish, this reprogramming step is achieved by the pioneer factors Nanog, Pou5f3, and Sox19b (NPS). Yet, it remains unclear whether cells lacking this reprogramming step are directed towards wild type states or towards novel developmental canals in the Waddington landscape of embryonic development. Here we investigate the developmental fate of embryonic cells mutant for NPS by analyzing their single-cell gene expression profiles. We find that cells lacking the first developmental reprogramming steps can acquire distinct cell states. These states are manifested by gene expression modules that result from a failure of nuclear reprogramming, the persistence of the maternal program, and the activation of somatic compensatory programs. As a result, most mutant cells follow new developmental canals and acquire new mixed cell states in development. In contrast, a group of mutant cells acquire primordial germ cell-like states, suggesting that NPS-dependent reprogramming is dispensable for these cell states. Together, these results demonstrate that developmental reprogramming after fertilization is required to differentiate most canonical developmental programs, and loss of the transient totipotent state canalizes embryonic cells into new developmental states in vivo.

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.

Perturbed development of calb2b expressing dI6 interneurons and motor neurons underlies locomotor defects observed in calretinin knock-down zebrafish larvae

Calcium binding proteins are essential for neural development and cellular activity. Calretinin, encoded by calb2a and calb2b, plays a role during early zebrafish development and has been proposed as a marker for distinct neuronal populations within the locomotor network. We generated a calb2b:hs:eGFP transgenic reporter line to characterize calretinin expressing cells in the developing spinal cord and describe morphological and behavioral defects in calretinin knock-down larvae. eGFP was detected in primary and secondary motor neurons, as well as in dI6 and V0v interneurons. Knock-down of calretinin lead to disturbed development of motor neurons and dI6 interneurons, revealing a crucial role during early development of the locomotor network. Primary motor neurons showed delayed axon outgrowth and the distinct inhibitory CoLo neurons, originating from the dI6 lineage, were absent. These observations explain the locomotor defects we observed in calretinin knock-down animals where the velocity, acceleration and coordination were affected during escapes. Altogether, our analysis suggests an essential role for calretinin during the development of the circuits regulating escape responses and fast movements within the locomotor network.

The development of brain pericytes requires expression of the transcription factor nkx3.1 in intermediate precursors

Brain pericytes are one of the critical cell types that regulate endothelial barrier function and activity, thus ensuring adequate blood flow to the brain. The genetic pathways guiding undifferentiated cells into mature pericytes are not well understood. We show here that pericyte precursor populations from both neural crest and head mesoderm of zebrafish express the transcription factor nkx3.1 develop into brain pericytes. We identify the gene signature of these precursors and show that an nkx3.1-, foxf2a-, and cxcl12b-expressing pericyte precursor population is present around the basilar artery prior to artery formation and pericyte recruitment. The precursors later spread throughout the brain and differentiate to express canonical pericyte markers. Cxcl12b-Cxcr4 signaling is required for pericyte attachment and differentiation. Further, both nkx3.1 and cxcl12b are necessary and sufficient in regulating pericyte number as loss inhibits and gain increases pericyte number. Through genetic experiments, we have defined a precursor population for brain pericytes and identified genes critical for their differentiation.

EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

Abcg2a is the functional homolog of human ABCG2 expressed at the zebrafish blood-brain barrier

BACKGROUND

A principal protective component of the mammalian blood-brain barrier (BBB) is the high expression of the multidrug efflux transporters P-glycoprotein (P-gp, encoded by ABCB1) and ABCG2 (encoded by ABCG2) on the lumenal surface of endothelial cells. The zebrafish P-gp homolog Abcb4 is expressed at the BBB and phenocopies human P-gp. Comparatively little is known about the four zebrafish homologs of the human ABCG2 gene: abcg2a, abcg2b, abcg2c, and abcg2d. Here we report the functional characterization and brain tissue distribution of zebrafish ABCG2 homologs.

METHODS

To determine substrates of the transporters, we stably expressed each in HEK-293 cells and performed cytotoxicity and fluorescent efflux assays with known ABCG2 substrates. To assess the expression of transporter homologs, we used a combination of RNAscope in situ hybridization probes and immunohistochemistry to stain paraffin-embedded sections of adult and larval zebrafish.

RESULTS

We found Abcg2a had the greatest substrate overlap with ABCG2, and Abcg2d appeared to be the least functionally similar. We identified abcg2a as the only homolog expressed at the adult and larval zebrafish BBB, based on its localization to claudin-5 positive brain vasculature.

CONCLUSIONS

These results demonstrate the conserved function of zebrafish Abcg2a and suggest that zebrafish may be an appropriate model organism for studying the role of ABCG2 at the BBB.

The cytoskeleton adaptor protein Sorbs1 controls the development of lymphatic and venous vessels in zebrafish

BACKGROUND

Lymphangiogenesis, the formation of lymphatic vessels, is tightly linked to the development of the venous vasculature, both at the cellular and molecular levels. Here, we identify a novel role for Sorbs1, the founding member of the SoHo family of cytoskeleton adaptor proteins, in vascular and lymphatic development in the zebrafish.

RESULTS

We show that Sorbs1 is required for secondary sprouting and emergence of several vascular structures specifically derived from the axial vein. Most notably, formation of the precursor parachordal lymphatic structures is affected in sorbs1 mutant embryos, severely impacting the establishment of the trunk lymphatic vessel network. Interestingly, we show that Sorbs1 interacts with the BMP pathway and could function outside of Vegfc signaling. Mechanistically, Sorbs1 controls FAK/Src signaling and subsequently impacts on the cytoskeleton processes regulated by Rac1 and RhoA GTPases. Inactivation of Sorbs1 altered cell-extracellular matrix (ECM) contacts rearrangement and cytoskeleton dynamics, leading to specific defects in endothelial cell migratory and adhesive properties.

CONCLUSIONS

Overall, using in vitro and in vivo assays, we identify Sorbs1 as an important regulator of venous and lymphatic angiogenesis independently of the Vegfc signaling axis. These results provide a better understanding of the complexity found within context-specific vascular and lymphatic development.

Plxnd1-mediated mechanosensing of blood flow controls the caliber of the Dorsal Aorta via the transcription factor Klf2

The cardiovascular system generates and responds to mechanical forces. The heartbeat pumps blood through a network of vascular tubes, which adjust their caliber in response to the hemodynamic environment. However, how endothelial cells in the developing vascular system integrate inputs from circulatory forces into signaling pathways to define vessel caliber is poorly understood. Using vertebrate embryos and in vitro-assembled microvascular networks of human endothelial cells as models, flow and genetic manipulations, and custom software, we reveal that Plexin-D1, an endothelial Semaphorin receptor critical for angiogenic guidance, employs its mechanosensing activity to serve as a crucial positive regulator of the Dorsal Aorta's (DA) caliber. We also uncover that the flow-responsive transcription factor KLF2 acts as a paramount mechanosensitive effector of Plexin-D1 that enlarges endothelial cells to widen the vessel. These findings illuminate the molecular and cellular mechanisms orchestrating the interplay between cardiovascular development and hemodynamic forces.

Progress in characterizing ABC multidrug transporters in zebrafish

Zebrafish have proved to be invaluable for modeling complex physiological processes shared by all vertebrate animals. Resistance of cancers and other diseases to drug treatment can occur owing to expression of the ATP-dependent multidrug transporters ABCB1, ABCG2, and ABCC1, either because of expression of these transporters by the target cells to reduce intracellular concentrations of cytotoxic drugs at barrier sites such as the blood-brain barrier (BBB) to limit penetration of drugs into privileged compartments, or by affecting the absorption, distribution, and excretion of drugs administered orally, through the skin, or directly into the bloodstream. We describe the drug specificity, cellular localization, and function of zebrafish orthologs of multidrug resistance ABC transporters with the goal of developing zebrafish models to explore the physiological and pathophysiological functions of these transporters. Finally, we provide context demonstrating the utility of zebrafish in studying cancer drug resistance. Our ultimate goal is to improve treatment of cancer and other diseases which are affected by ABC multidrug resistance transporters.