Optogenetic control of Nodal signaling patterns.

A crucial step in early embryogenesis is the establishment of spatial patterns of signaling activity. Tools to perturb morphogen signals with high resolution in space and time can help reveal how embryonic cells decode these signals to make appropriate fate decisions. Here, we present new optogenetic reagents and an experimental pipeline for creaHng designer Nodal signaling patterns in live zebrafish embryos. Nodal receptors were fused to the light-sensitive heterodimerizing pair Cry2/CIB1N, and the Type II receptor was sequestered to the cytosol. The improved optoNodal2 reagents eliminate dark activity and improve response kinetics, without sacrificing dynamic range. We adapted an ultra-widefield microscopy platform for parallel light patterning in up to 36 embryos and demonstrated precise spatial control over Nodal signaling activity and downstream gene expression. Patterned Nodal activation drove precisely controlled internalization of endodermal precursors. Further, we used patterned illumination to generate synthetic signaling patterns in Nodal signaling mutants, rescuing several characteristic developmental defects. This study establishes an experimental toolkit for systematic exploration of Nodal signaling patterns in live embryos.

A vertebrate family without a functional Hypocretin/Orexin arousal system.

The Hypocretin/Orexin signaling pathway suppresses sleep and promotes arousal, whereas the loss of Hypocretin/Orexin results in narcolepsy, including the involuntary loss of muscle tone (cataplexy).1 Here, we show that the South Asian fish species Chromobotia macracanthus exhibits a sleep-like state during which individuals stop swimming and rest on their side. Strikingly, we discovered that the Hypocretin/Orexin system is pseudogenized in C. macracanthus, but in contrast to Hypocretin-deficient mammals, C. macracanthus does not suffer from sudden behavioral arrests. Similarly, zebrafish mutations in hypocretin/orexin show no evident signs of cataplectic-like episodes. Notably, four additional species in the Botiidae family also lack a functional Hypocretin/Orexin system. These findings identify the first vertebrate family that does not rely on a functional Hypocretin/Orexin system for the regulation of sleep and arousal.

A single-cell time-lapse of mouse prenatal development from gastrula to birth.

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Molecular mechanisms controlling the biogenesis of the TGF-ß signal Vg1.

The TGF-beta signals Vg1 (Dvr1/Gdf3) and Nodal form heterodimers to induce vertebrate mesendoderm. The Vg1 proprotein is a monomer retained in the endoplasmic reticulum (ER) and is processed and secreted upon heterodimerization with Nodal, but the mechanisms underlying Vg1 biogenesis are largely elusive. Here, we clarify the mechanisms underlying Vg1 retention, processing, secretion, and signaling and introduce a Synthetic Processing (SynPro) system that enables the programmed cleavage of ER-resident and extracellular proteins. First, we find that Vg1 can be processed by intra- or extracellular proteases. Second, Vg1 can be processed without Nodal but requires Nodal for secretion and signaling. Third, Vg1-Nodal signaling activity requires Vg1 processing, whereas Nodal can remain unprocessed. Fourth, Vg1 employs exposed cysteines, glycosylated asparagines, and BiP chaperone-binding motifs for monomer retention in the ER. These observations suggest two mechanisms for rapid mesendoderm induction: Chaperone-binding motifs help store Vg1 as an inactive but ready-to-heterodimerize monomer in the ER, and the flexibility of Vg1 processing location allows efficient generation of active heterodimers both intra- and extracellularly. These results establish SynPro as an in vivo processing system and define molecular mechanisms and motifs that facilitate the generation of active TGF-beta heterodimers.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup.

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

Gene module reconstruction elucidates cellular differentiation processes and the regulatory logic of specialized secretion.

During differentiation, cells become structurally and functionally specialized, but comprehensive views of the underlying remodeling processes are elusive. Here, we leverage scRNA-seq developmental trajectories to reconstruct differentiation using two secretory tissues as a model system - the zebrafish notochord and hatching gland. First, we present an approach to integrate expression and functional similarities for gene module identification, revealing dozens of gene modules representing known and newly associated differentiation processes and their temporal ordering. Second, we focused on the unfolded protein response (UPR) transducer module to study how general versus cell-type specific secretory functions are regulated. By profiling loss- and gain-of-function embryos, we found that the UPR transcription factors creb3l1, creb3l2, and xbp1 are master regulators of a general secretion program. creb3l1/creb3l2 additionally activate an extracellular matrix secretion program, while xbp1 partners with bhlha15 to activate a gland-specific secretion program. Our study offers a multi-source integrated approach for functional gene module identification and illustrates how transcription factors confer general and specialized cellular functions.

Basic science under threat: Lessons from the Skirball Institute.

Support for basic science has been eclipsed by initiatives aimed at specific medical problems. The latest example is the dismantling of the Skirball Institute at NYU School of Medicine. Here, we reflect on the achievements and mission underlying the Skirball to gain insight into the dividends of maintaining a basic science vision within the academic enterprises.

Gene family evolution underlies cell-type diversification in the hypothalamus of teleosts.

Hundreds of cell types form the vertebrate brain but it is largely unknown how similar cellular repertoires are between or within species or how cell-type diversity evolves. To examine cell-type diversity across and within species, we performed single-cell RNA sequencing of ~130,000 hypothalamic cells from zebrafish (Danio rerio) and surface and cave morphs of Mexican tetra (Astyanax mexicanus). We found that over 75% of cell types were shared between zebrafish and Mexican tetra, which diverged from a common ancestor over 150 million years ago. Shared cell types displayed shifts in paralogue expression that were generated by subfunctionalization after genome duplication. Expression of terminal effector genes, such as neuropeptides, was more conserved than the expression of their associated transcriptional regulators. Species-specific cell types were enriched for the expression of species-specific genes and characterized by the neofunctionalization of expression patterns of members of recently expanded or contracted gene families. Comparisons between surface and cave morphs revealed differences in immune repertoires and transcriptional changes in neuropeptidergic cell types associated with genomic differences. The single-cell atlases presented here are a powerful resource to explore hypothalamic cell types and reveal how gene family evolution and shifts in paralogue expression contribute to cellular diversity.

The pattern of nodal morphogen signaling is shaped by co-receptor expression.

Embryos must communicate instructions to their constituent cells over long distances. These instructions are often encoded in the concentration of signals called morphogens. In the textbook view, morphogen molecules diffuse from a localized source to form a concentration gradient, and target cells adopt fates by measuring the local morphogen concentration. However, natural patterning systems often incorporate numerous co-factors and extensive signaling feedback, suggesting that embryos require additional mechanisms to generate signaling patterns. Here, we examine the mechanisms of signaling pattern formation for the mesendoderm inducer Nodal during zebrafish embryogenesis. We find that Nodal signaling activity spans a normal range in the absence of signaling feedback and relay, suggesting that diffusion is sufficient for Nodal gradient formation. We further show that the range of endogenous Nodal ligands is set by the EGF-CFC co-receptor Oep: in the absence of Oep, Nodal activity spreads to form a nearly uniform distribution throughout the embryo. In turn, increasing Oep levels sensitizes cells to Nodal ligands. We recapitulate these experimental results with a computational model in which Oep regulates the diffusive spread of Nodal ligands by setting the rate of capture by target cells. This model predicts, and we confirm in vivo, the surprising observation that a failure to replenish Oep transforms the Nodal signaling gradient into a travelling wave. These results reveal that patterns of Nodal morphogen signaling are shaped by co-receptor-mediated restriction of ligand spread and sensitization of responding cells.

Emergence of Neuronal Diversity during Vertebrate Brain Development.

Neurogenesis comprises many highly regulated processes including proliferation, differentiation, and maturation. However, the transcriptional landscapes underlying brain development are poorly characterized. We describe a developmental single-cell catalog of ∼220,000 zebrafish brain cells encompassing 12 stages from embryo to larva. We characterize known and novel gene markers for ∼800 clusters and provide an overview of the diversification of neurons and progenitors across these time points. We also introduce an optimized GESTALT lineage recorder that enables higher expression and recovery of Cas9-edited barcodes to query lineage segregation. Cell type characterization indicates that most embryonic neural progenitor states are transitory and transcriptionally distinct from neural progenitors of post-embryonic stages. Reconstruction of cell specification trajectories reveals that late-stage retinal neural progenitors transcriptionally overlap cell states observed in the embryo. The zebrafish brain development atlas provides a resource to define and manipulate specific subsets of neurons and to uncover the molecular mechanisms underlying vertebrate neurogenesis.

Zebrafish dscaml1 Deficiency Impairs Retinal Patterning and Oculomotor Function.

Down syndrome cell adhesion molecules (dscam and dscaml1) are essential regulators of neural circuit assembly, but their roles in vertebrate neural circuit function are still mostly unexplored. We investigated the functional consequences of dscaml1 deficiency in the larval zebrafish (sexually undifferentiated) oculomotor system, where behavior, circuit function, and neuronal activity can be precisely quantified. Genetic perturbation of dscaml1 resulted in deficits in retinal patterning and light adaptation, consistent with its known roles in mammals. Oculomotor analyses revealed specific deficits related to the dscaml1 mutation, including severe fatigue during gaze stabilization, reduced saccade amplitude and velocity in the light, greater disconjugacy, and impaired fixation. Two-photon calcium imaging of abducens neurons in control and dscaml1 mutant animals confirmed deficits in saccade-command signals (indicative of an impairment in the saccadic premotor pathway), whereas abducens activation by the pretectum-vestibular pathway was not affected. Together, we show that loss of dscaml1 resulted in impairments in specific oculomotor circuits, providing a new animal model to investigate the development of oculomotor premotor pathways and their associated human ocular disorders.SIGNIFICANCE STATEMENTDscaml1 is a neural developmental gene with unknown behavioral significance. Using the zebrafish model, this study shows that dscaml1 mutants have a host of oculomotor (eye movement) deficits. Notably, the oculomotor phenotypes in dscaml1 mutants are reminiscent of human ocular motor apraxia, a neurodevelopmental disorder characterized by reduced saccade amplitude and gaze stabilization deficits. Population-level recording of neuronal activity further revealed potential subcircuit-specific requirements for dscaml1 during oculomotor behavior. These findings underscore the importance of dscaml1 in the development of visuomotor function and characterize a new model to investigate potential circuit deficits underlying human oculomotor disorders.

Convergent Temperature Representations in Artificial and Biological Neural Networks.

Discoveries in biological neural networks (BNNs) shaped artificial neural networks (ANNs) and computational parallels between ANNs and BNNs have recently been discovered. However, it is unclear to what extent discoveries in ANNs can give insight into BNN function. Here, we designed and trained an ANN to perform heat gradient navigation and found striking similarities in computation and heat representation to a known zebrafish BNN. This included shared ON- and OFF-type representations of absolute temperature and rates of change. Importantly, ANN function critically relied on zebrafish-like units. We furthermore used the accessibility of the ANN to discover a new temperature-responsive cell type in the zebrafish cerebellum. Finally, constraining the ANN by the C. elegans motor repertoire retuned sensory representations indicating that our approach generalizes. Together, these results emphasize convergence of ANNs and BNNs on stereotypical representations and that ANNs form a powerful tool to understand their biological counterparts.

Zebrafish oxytocin neurons drive nocifensive behavior via brainstem premotor targets.

Animals have evolved specialized neural circuits to defend themselves from pain- and injury-causing stimuli. Using a combination of optical, behavioral and genetic approaches in the larval zebrafish, we describe a novel role for hypothalamic oxytocin (OXT) neurons in the processing of noxious stimuli. In vivo imaging revealed that a large and distributed fraction of zebrafish OXT neurons respond strongly to noxious inputs, including the activation of damage-sensing TRPA1 receptors. OXT population activity reflects the sensorimotor transformation of the noxious stimulus, with some neurons encoding sensory information and others correlating more strongly with large-angle swims. Notably, OXT neuron activation is sufficient to generate this defensive behavior via the recruitment of brainstem premotor targets, whereas ablation of OXT neurons or loss of the peptide attenuates behavioral responses to TRPA1 activation. These data highlight a crucial role for OXT neurons in the generation of appropriate defensive responses to noxious input.

Distributed Plasticity Drives Visual Habituation Learning in Larval Zebrafish.

Habituation is a simple form of learning where animals learn to reduce their responses to repeated innocuous stimuli [1]. Habituation is thought to occur via at least two temporally and molecularly distinct mechanisms, which lead to short-term memories that last for seconds to minutes and long-term memories that last for hours or longer [1, 2]. Here, we focus on long-term habituation, which, due to the extended time course, necessitates stable alterations to circuit properties [2-4]. In its simplest form, long-term habituation could result from a plasticity event at a single point in a circuit, and many studies have focused on identifying the site and underlying mechanism of plasticity [5-10]. However, it is possible that these individual sites are only one of many points in the circuit where plasticity is occurring. Indeed, studies of short-term habituation in C. elegans indicate that in this paradigm, multiple genetically separable mechanisms operate to adapt specific aspects of behavior [11-13]. Here, we use a visual assay in which larval zebrafish habituate their response to sudden reductions in illumination (dark flashes) [14, 15]. Through behavioral analyses, we find that multiple components of the dark-flash response habituate independently of one another using different molecular mechanisms. This is consistent with a modular model in which habituation originates from multiple independent processes, each adapting specific components of behavior. This may allow animals to more specifically or flexibly habituate based on stimulus context or internal states.

Phenotypic Landscape of Schizophrenia-Associated Genes Defines Candidates and Their Shared Functions.

Genomic studies have identified hundreds of candidate genes near loci associated with risk for schizophrenia. To define candidates and their functions, we mutated zebrafish orthologs of 132 human schizophrenia-associated genes. We created a phenotype atlas consisting of whole-brain activity maps, brain structural differences, and profiles of behavioral abnormalities. Phenotypes were diverse but specific, including altered forebrain development and decreased prepulse inhibition. Exploration of these datasets identified promising candidates in more than 10 gene-rich regions, including the magnesium transporter cnnm2 and the translational repressor gigyf2, and revealed shared anatomical sites of activity differences, including the pallium, hypothalamus, and tectum. Single-cell RNA sequencing uncovered an essential role for the understudied transcription factor znf536 in the development of forebrain neurons implicated in social behavior and stress. This phenotypic landscape of schizophrenia-associated genes prioritizes more than 30 candidates for further study and provides hypotheses to bridge the divide between genetic association and biological mechanism.

Individual long non-coding RNAs have no overt functions in zebrafish embryogenesis, viability and fertility.

Hundreds of long non-coding RNAs (lncRNAs) have been identified as potential regulators of gene expression, but their functions remain largely unknown. To study the role of lncRNAs during vertebrate development, we selected 25 zebrafish lncRNAs based on their conservation, expression profile or proximity to developmental regulators, and used CRISPR-Cas9 to generate 32 deletion alleles. We observed altered transcription of neighboring genes in some mutants, but none of the lncRNAs were required for embryogenesis, viability or fertility. Even RNAs with previously proposed non-coding functions (cyrano and squint) and other conserved lncRNAs (gas5 and lnc-setd1ba) were dispensable. In one case (lnc-phox2bb), absence of putative DNA regulatory-elements, but not of the lncRNA transcript itself, resulted in abnormal development. LncRNAs might have redundant, subtle, or context-dependent roles, but extrapolation from our results suggests that the majority of individual zebrafish lncRNAs have no overt roles in embryogenesis, viability and fertility.

Conserved regulation of Nodal-mediated left-right patterning in zebrafish and mouse.

Nodal is the major effector of left-right axis development. In mice, Nodal forms heterodimers with Gdf1 and is inhibited by Cerl2/Dand5 at the node, and by Lefty1 in the lateral plate mesoderm (LPM). Studies in zebrafish have suggested some parallels, but also differences, between left-right patterning in mouse and zebrafish. To address these discrepancies, we generated single and double zebrafish mutants for southpaw (spaw, the Nodal ortholog), dand5 and lefty1, and performed biochemical and activity assays with Spaw and Vg1/Gdf3 (the Gdf1 ortholog). Contrary to previous findings, spaw mutants failed to initiate spaw expression in the LPM, and asymmetric heart looping was absent, similar to mouse Nodal mutants. In blastoderm assays, Vg1 and Spaw were interdependent for target gene induction, and contrary to previous results, formed heterodimers. Loss of Dand5 or Lefty1 caused bilateral spaw expression, similar to mouse mutants, and Lefty1 was replaceable with a uniform Nodal signaling inhibitor. Collectively, these results indicate that Dand5 activity biases Spaw-Vg1 heterodimer activity to the left, Spaw around Kupffer's vesicle induces the expression of spaw in the LPM and global Nodal inhibition maintains the left bias of Spaw activity, demonstrating conservation between zebrafish and mouse mechanisms of left-right patterning.

Large-scale reconstruction of cell lineages using single-cell readout of transcriptomes and CRISPR-Cas9 barcodes by scGESTALT.

Lineage relationships among the large number of heterogeneous cell types generated during development are difficult to reconstruct in a high-throughput manner. We recently established a method, scGESTALT, that combines cumulative editing of a lineage barcode array by CRISPR-Cas9 with large-scale transcriptional profiling using droplet-based single-cell RNA sequencing (scRNA-seq). The technique generates edits in the barcode array over multiple timepoints using Cas9 and pools of single-guide RNAs (sgRNAs) introduced during early and late zebrafish embryonic development, which distinguishes it from similar Cas9 lineage-tracing methods. The recorded lineages are captured, along with thousands of cellular transcriptomes, to build lineage trees with hundreds of branches representing relationships among profiled cell types. Here, we provide details for (i) generating transgenic zebrafish; (ii) performing multi-timepoint barcode editing; (iii) building scRNA-seq libraries from brain tissue; and (iv) concurrently amplifying lineage barcodes from captured single cells. Generating transgenic lines takes 6 months, and performing barcode editing and generating single-cell libraries involve 7 d of hands-on time. scGESTALT provides a scalable platform to map lineage relationships between cell types in any system that permits genome editing during development, regeneration, or disease.

Scale-invariant patterning by size-dependent inhibition of Nodal signalling.

Individuals can vary substantially in size, but the proportions of their body plans are often maintained. We generated smaller zebrafish by removing 30% of their cells at the blastula stages and found that these embryos developed into normally patterned individuals. Strikingly, the proportions of all germ layers adjusted to the new embryo size within 2 hours after cell removal. As Nodal-Lefty signalling controls germ-layer patterning, we performed a computational screen for scale-invariant models of this activator-inhibitor system. This analysis predicted that the concentration of the highly diffusive inhibitor Lefty increases in smaller embryos, leading to a decreased Nodal activity range and contracted germ-layer dimensions. In vivo studies confirmed that Lefty concentration increased in smaller embryos, and embryos with reduced Lefty levels or with diffusion-hindered Lefty failed to scale their tissue proportions. These results reveal that size-dependent inhibition of Nodal signalling allows scale-invariant patterning.