Search

Mammalian embryogenesis is characterized by rapid cellular proliferation and diversification. Within a few weeks, a single-cell zygote gives rise to millions of cells expressing a panoply of molecular programs. Although intensively studied, a comprehensive delineation of the major cellular trajectories that comprise mammalian development in vivo remains elusive. Here, we set out to integrate several single-cell RNA-sequencing (scRNA-seq) datasets that collectively span mouse gastrulation and organogenesis, supplemented with new profiling of ~150,000 nuclei from approximately embryonic day 8.5 (E8.5) embryos staged in one-somite increments. Overall, we define cell states at each of 19 successive stages spanning E3.5 to E13.5 and heuristically connect them to their pseudoancestors and pseudodescendants. Although constructed through automated procedures, the resulting directed acyclic graph (TOME (trajectories of mammalian embryogenesis)) is largely consistent with our contemporary understanding of mammalian development. We leverage TOME to systematically nominate transcription factors (TFs) as candidate regulators of each cell type’s specification, as well as ‘cell-type homologs’ across vertebrate evolution.

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 thein vivodifferentiation 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.

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

Regeneration after peripheral nerve damage requires that axons re-grow to the correct target tissues in a process called target-specific regeneration. Although much is known about the mechanisms that promote axon re-growth, re-growing axons often fail to reach the correct targets, resulting in impaired nerve function. We know very little about how axons achieve target-specific regeneration, particularly in branched nerves that require distinct targeting decisions at branch points. The zebrafish vagus motor nerve is a branched nerve with a well-defined topographic organization. Here, we track regeneration of individual vagus axons after whole-nerve laser severing and find a robust capacity for target-specific, functional re-growth. We then develop a new single-cell chimera injury model for precise manipulation of axon-environment interactions and find that (1) the guidance mechanism used during regeneration is distinct from the nerve's developmental guidance mechanism, (2) target selection is specified by neurons' intrinsic memory of their position within the brain, and (3) targeting to a branch requires its pre-existing innervation. This work establishes the zebrafish vagus nerve as a tractable regeneration model and reveals the mechanistic basis of target-specific regeneration.

Planar cell polarity (PCP) plays crucial roles in developmental processes such as gastrulation, neural tube closure and hearing. Wnt pathway mutants are often classified as PCP mutants due to similarities between their phenotypes. Here, we show that in the zebrafish lateral line, disruptions of the PCP and Wnt pathways have differential effects on hair cell orientations. While mutations in the PCP genes vangl2 and scrib cause random orientations of hair cells, mutations in wnt11f1, gpc4 and fzd7a/b induce hair cells to adopt a concentric pattern. This concentric pattern is not caused by defects in PCP but is due to misaligned support cells. The molecular basis of the support cell defect is unknown but we demonstrate that the PCP and Wnt pathways work in parallel to establish proper hair cell orientation. Consequently, hair cell orientation defects are not solely explained by defects in PCP signaling, and some hair cell phenotypes warrant re-evaluation., Planar cell polarity (PCP) regulates hair cell orientation in the zebrafish lateral line. Here, the authors show that mutating Wnt pathway genes (wnt11f1, fzd7a/b, and gpc4) causes concentric hair cell patterns not regulated by PCP, thus showing PCP/Wnt pathway genes have different consequences on hair cell orientation.

The asymmetric localization of planar cell polarity (PCP) proteins is essential for the establishment of many planar polarized cellular processes, but the mechanisms that maintain these asymmetric distributions remain poorly understood. A body of evidence has tied oriented subapical microtubules (MTs) to the establishment of PCP protein polarity, yet recent studies have suggested that the MT cytoskeleton is later dispensable for the maintenance of this asymmetry. As MTs underlie the vesicular trafficking of membrane-bound proteins within cells, the requirement for MTs in the maintenance of PCP merited further investigation. We investigated the complex interactions between PCP proteins and the MT cytoskeleton in the polarized context of the floorplate of the zebrafish neural tube. We demonstrated that the progressive posterior polarization of the primary cilia of floorplate cells requires not only Vangl2 but also Fzd3a. We determined that GFP-Vangl2 asymmetrically localizes to anterior membranes whereas Fzd3a-GFP does not polarize on anterior or posterior membranes but maintains a cytosolic enrichment at the base of the primary cilium. Vesicular Fzd3a-GFP is rapidly trafficked along MTs primarily toward the apical membrane during a period of PCP maintenance, whereas vesicular GFP-Vangl2 is less frequently observed. Nocodazole-induced loss of MT polymerization disrupts basal body positioning as well as GFP-Vangl2 localization and reduces cytosolic Fzd3a-GFP movements. Removal of nocodazole after MT disruption restores MT polymerization but does not restore basal body polarity. Interestingly, GFP-Vangl2 repolarizes to anterior membranes and vesicular Fzd3a-GFP dynamics recover after multiple hours of recovery, even in the context of unpolarized basal bodies. Together our findings challenge previous work by revealing an ongoing role for MT-dependent transport of PCP proteins in maintaining both cellular and PCP protein asymmetry during development.

The cerebellum coordinates vestibular input into the hindbrain to control balance and movement, and its anatomical complexity is increasingly viewed as a high-throughput processing center for sensory and cognitive functions. Cerebellum development however is relatively simple, and arises from a specialized structure in the anterior hindbrain called the rhombic lip, which along with the ventricular zone of the rostral-most dorsal hindbrain region, give rise to the distinct cell types that constitute the cerebellum. Granule cells, being the most numerous cell types, arise from the rhombic lip and form a dense and distinct layer of the cerebellar cortex. In this short review, we describe the various strategies used by amniotes and anamniotes to generate and diversify granule cell types during cerebellar development.

Necrosis of macrophages in the tuberculous granuloma represents a major pathogenic event in tuberculosis. Through a zebrafish forward genetic screen, we identified the mTOR kinase, a master regulator of metabolism, as an early host resistance factor in tuberculosis. We found that mTOR complex 1 protects macrophages from mycobacterium-induced death by enabling infection-induced increases in mitochondrial energy metabolism fueled by glycolysis. These metabolic adaptations are required to prevent mitochondrial damage and death caused by the secreted mycobacterial virulence determinant ESAT-6. Thus, the host can effectively counter this early critical mycobacterial virulence mechanism simply by regulating energy metabolism, thereby allowing pathogen-specific immune mechanisms time to develop. Our findings may explain why Mycobacterium tuberculosis, albeit humanity’s most lethal pathogen, is successful in only a minority of infected individuals.

During vertebrate central nervous system development, most oligodendrocyte progenitor cells (OPCs) are specified in the ventral spinal cord and must migrate throughout the neural tube until they become evenly distributed, occupying non-overlapping domains. While this process of developmental OPC migration is well characterized, the nature of the molecular mediators that govern it remain largely unknown. Here, using zebrafish as a model, we demonstrate that Met signaling is required for initial developmental migration of OPCs, and, using cell-specific knock-down of Met signaling, show that Met acts cell-autonomously in OPCs. Taken together, these findings demonstrate in vivo, the role of Met signaling in OPC migration and provide new insight into how OPC migration is regulated during development.

Mammalian embryogenesis is characterized by rapid cellular proliferation and diversification. Within a few weeks, a single cell zygote gives rise to millions of cells expressing a panoply of molecular programs, including much of the diversity that will subsequently be present in adult tissues. Although intensively studied, a comprehensive delineation of the major cellular trajectories that comprise mammalian development in vivo remains elusive. Here we set out to integrate several single cell RNA-seq datasets (scRNA-seq) that collectively span mouse gastrulation and organogenesis. We define cell states at each of 19 successive stages spanning E3.5 to E13.5, heuristically connect them with their pseudo-ancestors and pseudo-descendants, and for a subset of stages, deconvolve their approximate spatial distributions. Despite being constructed through automated procedures, the resulting trajectories of mammalian embryogenesis (TOME) are largely consistent with our contemporary understanding of mammalian development. We leverage TOME to nominate transcription factors (TF) and TF motifs as key regulators of each branch point at which a new cell type emerges. Finally, to facilitate comparisons across vertebrates, we apply the same procedures to single cell datasets of zebrafish and frog embryogenesis, and nominate “cell type homologs” based on shared regulators and transcriptional states.

Necrosis of macrophages in the tuberculous granuloma represents a major pathogenic event in tuberculosis. Through a zebrafish forward genetic screen, we identified the mTOR kinase, a master regulator of metabolism, as an early host resistance factor in tuberculosis. We found that mTOR protects macrophages from mycobacterium-induced death by enabling infection-induced increases in mitochondrial energy metabolism fueled by glycolysis. These metabolic adaptations were required to prevent mitochondrial damage and death caused specifically by the mycobacterial ESX-1 secretion system. Our finding that the host can effectively counter this early critical mycobacterial virulence mechanism simply by regulating energy metabolism may help explain why Mycobacterium tuberculosis, albeit humanity’s most lethal pathogen, only causes disease in a minority of infected individuals. This work was supported by a Wellcome Trust Principal Research Fellowship (103950/Z/14) and an NIH MERIT award (R37 AI054503) (L.R.).

SUMMARYThe topographic map, in which the positions of neuron cell bodies correspond with the positions of their synaptic targets, is a major organizational motif in the nervous system. To understand how topographic axon targeting is controlled during development, we examine the mechanism underlying topographic innervation of the pharyngeal arches by the vagus motor nerve in zebrafish. We reveal that Retinoic Acid organizes the topographic map by specifying anterior-posterior identity in post-mitotic vagus motor neurons. We then show that chemoattractant signaling between hepatocyte growth factor (Hgf) and the Met receptor is required for pharyngeal arch innervation by the vagus motor nerve. Finally, we find that Retinoic Acid controls the spatiotemporal dynamics of Hgf/Met signaling to coordinate axon targeting with the developmental progression of the pharyngeal arches and show that experimentally altering the timing of Hgf/Met signaling is sufficient to redirect axon targeting and disrupt the topographic map. These findings establish a new mechanism of topographic map development in which regulation of chemoattractant signaling in both space and time guides axon targeting.

Information flow through neural circuits often requires their organization into topographic maps in which the positions of cell bodies and synaptic targets correspond. To understand how topographic map development is controlled, we examine the mechanism underlying targeting of vagus motor axons to the pharyngeal arches in zebrafish. We reveal that retinoic acid organizes topography by specifying anterior-posterior identity in vagus motor neurons. We then show that chemoattractant signaling between Hgf and Met is required for vagus innervation of the pharyngeal arches. Finally, we find that retinoic acid controls the spatiotemporal dynamics of Hgf/Met signaling to coordinate axon targeting with the developmental progression of the pharyngeal arches and show that experimentally altering the timing of Hgf/Met signaling is sufficient to redirect axon targeting and disrupt the topographic map. These findings establish a mechanism of topographic map development in which the regulation of chemoattractant signaling in space and time guides axon targeting.

The zebrafish (Danio rerio) has emerged as a powerful model to study vertebrate development and disease. Its short generation time makes it amenable to genetic manipulation and analysis, and its small size and high fecundity make it especially well suited for large-scale forward genetic and chemical screens. Fast-developing zebrafish embryos are transparent, facilitating live imaging of a variety of developmental processes in wild-type and mutant animals. Published version

The authors would like to thank the Royal Society of London, the National Eye Research Centre, the Visual Research Trust, Fight for Sight, the W.H. Ross Foundation, the Rosetrees Trust, and the Glasgow Children’s Hospital Charity for supporting this work. This work was also supported by the Deanship of Scientific Research at King Saud University for funding this research (Research Project) grant number ‘RGP – VPP – 219’. Mutations in the RPGR-interacting protein 1 (RPGRIP1) gene cause recessive Leber congenital amaurosis (LCA), juvenile retinitis pigmentosa (RP) and cone-rod dystrophy. RPGRIP1 interacts with other retinal disease-causing proteins and has been proposed to have a role in ciliary protein transport; however, its function remains elusive. Here, we describe a new zebrafish model carrying a nonsense mutation in the rpgrip1 gene. Rpgrip1homozygous mutants do not form rod outer segments and display mislocalization of rhodopsin, suggesting a role for RPGRIP1 in rhodopsin-bearing vesicle trafficking. Furthermore, Rab8, the key regulator of rhodopsin ciliary trafficking, was mislocalized in photoreceptor cells of rpgrip1 mutants. The degeneration of rod cells is early onset, followed by the death of cone cells. These phenotypes are similar to that observed in LCA and juvenile RP patients. Our data indicate RPGRIP1 is necessary for rod outer segment development through regulating ciliary protein trafficking. The rpgrip1 mutant zebrafish may provide a platform for developing therapeutic treatments for RP patients. Publisher PDF

Normal cells acquire aggressive behavior by modifying signaling pathways. For instance, alteration of endocytosis profoundly impacts both proliferation and migration during tumorigenesis. Here we investigate the mechanisms that enable the endocytic machinery to coordinate these processes. We show that a membrane curvature-sensing protein, endophilin A3, promotes growth and migration of colon cancer cells through two competing mechanisms — an endocytosis pathway that is required for proliferation, and a GTPase regulatory pathway that controls cell motility. EndoA3 stimulates cell migration by binding the Rac GEF TIAM1 leading to activation of small GTPases. Competing interactions of EndoA3 with membrane versus TIAM1 modulate hyperproliferative and metastatic phenotypes. Disruption of EndoA3-membrane interactions stimulates TIAM1 and small GTPases in vitro, and further promotes prometastatic phenotypes in vivo. Together, these results uncover a coupling mechanism, by which EndoA3 promotes growth and migration of colon cancers, by linking membrane dynamics to GTPase regulation.

Neural network function is based upon the patterns and types of connections made between neurons. Neuronal synapses are adhesions specialized for communication and they come in two types, chemical and electrical. Communication at chemical synapses occurs via neurotransmitter release whereas electrical synapses utilize gap junctions for direct ionic and metabolic coupling. Electrical synapses are often viewed as symmetrical structures, with the same components making both sides of the gap junction. By contrast, we show that a broad set of electrical synapses in zebrafish, Danio rerio, require two gap-junction-forming Connexins for formation and function. We find that one Connexin functions presynaptically while the other functions postsynaptically in forming the channels. We also show that these synapses are required for the speed and coordination of escape responses. Our data identify a genetic basis for molecular asymmetry at vertebrate electrical synapses and show they are required for appropriate behavioral performance.

Many networks throughout the nervous system are organized into topographic maps, where the positions of neuron cell bodies in the projecting field correspond with the positions of their axons in the target field. Previous studies of topographic map development show evidence for spatial patterning mechanisms, in which molecular determinants expressed across the projecting and target fields are matched directly in a point-to-point mapping process. Here, we describe a novel temporal mechanism of topographic map formation that depends on spatially regulated differences in the timing of axon outgrowth and functions in parallel with spatial point-to-point mapping mechanisms. We focus on the vagus motor neurons, which are topographically arranged in both mammals and fish. We show that cell position along the anterior-posterior axis of hindbrain rhombomere 8 determines expression of hox5 genes, which are expressed in posterior, but not anterior, vagus motor neurons. Using live imaging and transplantation in zebrafish embryos, we additionally reveal that axon initiation is delayed in posterior vagus motor neurons independent of neuron birth time. We show that hox5 expression directs topographic mapping without affecting time of axon outgrowth and that time of axon outgrowth directs topographic mapping without affecting hox5 expression. The vagus motor neuron topographic map is therefore determined by two mechanisms that act in parallel: a hox5-dependent spatial mechanism akin to classic mechanisms of topographic map formation and a novel axon outgrowth-dependent temporal mechanism in which time of axon formation is spatially regulated to direct axon targeting.

Background Schwann cells, which arise from the neural crest, are the myelinating glia of the peripheral nervous system. During development neural crest and their Schwann cell derivatives engage in a sequence of events that comprise delamination from the neuroepithelium, directed migration, axon ensheathment and myelin membrane synthesis. At each step neural crest and Schwann cells are polarized, implying important roles for molecules that create cellular asymmetries. In this work we investigated the possibility that one polarity protein, Pard3, contributes to the polarized features of neural crest and Schwann cells that are associated with directed migration and myelination.

Hox genes are classically ascribed to function in patterning the anterior-posterior axis of bilaterian animals; however, their role in directing molecular mechanisms underlying morphogenesis at the cellular level remains largely unstudied. We unveil a non-classical role for the zebrafish hoxb1b gene, which shares ancestral functions with mammalian Hoxa1, in controlling progenitor cell shape and oriented cell division during zebrafish anterior hindbrain neural tube morphogenesis. This is likely distinct from its role in cell fate acquisition and segment boundary formation. We show that, without affecting major components of apico-basal or planar cell polarity, Hoxb1b regulates mitotic spindle rotation during the oriented neural keel symmetric mitoses that are required for normal neural tube lumen formation in the zebrafish. This function correlates with a non-cell-autonomous requirement for Hoxb1b in regulating microtubule plus-end dynamics in progenitor cells in interphase. We propose that Hox genes can influence global tissue morphogenesis by control of microtubule dynamics in individual cells in vivo.

Brain pericytes are important regulators of brain vascular integrity, permeability and blood flow. Deficiencies of brain pericytes are associated with neonatal intracranial hemorrhage in human fetuses, as well as stroke and neurodegeneration in adults. Despite the important functions of brain pericytes, the mechanisms underlying their development are not well understood and little is known about how pericyte density is regulated across the brain. The Notch signaling pathway has been implicated in pericyte development, but its exact roles remain ill defined. Here, we report an investigation of the Notch3 receptor using zebrafish as a model system. We show that zebrafish brain pericytes express notch3 and that notch3 mutant zebrafish have a deficit of brain pericytes and impaired blood-brain barrier function. Conditional loss- and gain-of-function experiments provide evidence that Notch3 signaling positively regulates brain pericyte proliferation. These findings establish a new role for Notch signaling in brain vascular development whereby Notch3 signaling promotes expansion of the brain pericyte population.

Phenotypic robustness requires a process of developmental buffering that is largely not understood, but which can be disrupted by mutations. Here we show that in mef2cab1086 loss of function mutant embryos and early larvae, development of craniofacial hyoid bones, the opercle (Op) and branchiostegal ray (BR), becomes remarkably unstable; the large magnitude of the instability serves as a positive attribute to learn about features of this developmental buffering. The OpBR mutant phenotype variably includes bone expansion and fusion, Op duplication, and BR homeosis. Formation of a novel bone strut, or a bone bridge connecting the Op and BR together occurs frequently. We find no evidence that the phenotypic stability in the wild type is provided by redundancy between mef2ca and its co-ortholog mef2cb, or that it is related to the selector (homeotic) gene function of mef2ca. Changes in dorsal–ventral patterning of the hyoid arch also might not contribute to phenotypic instability in mutants. However, subsequent development of the bone lineage itself, including osteoblast differentiation and morphogenetic outgrowth, shows marked variation. Hence, steps along the developmental trajectory appear differentially sensitive to the loss of buffering, providing focus for the future study.

Muller glia are the resident radial glia in the vertebrate retina. The response of mammalian Muller glia to retinal damage often results in a glial scar and no functional replacement of lost neurons. Adult zebrafish Muller glia, in contrast, are considered tissue-specific stem cells that can self-renew and generate neurogenic progenitors to regenerate all retinal neurons after damage. Here, we demonstrate that regulation of TGFβ signaling by the corepressors Tgif1 and Six3b is critical for the proliferative response to photoreceptor destruction in the adult zebrafish retina. When function of these corepressors is disrupted, Muller glia and their progeny proliferate less, leading to a significant reduction in photoreceptor regeneration. Tgif1 expression and regulation of TGFβ signaling are implicated in the function of several types of stem cells, but this is the first demonstration that this regulatory network is necessary for regeneration of neurons.

The planar cell polarity (PCP) pathway is best known for its role in polarizing epithelial cells within the plane of a tissue but it also plays a role in a range of cell migration events during development. The mechanism by which the PCP pathway polarizes stationary epithelial cells is well characterized, but how PCP signaling functions to regulate more dynamic cell behaviors during directed cell migration is much less understood. Here, we review recent discoveries regarding the localization of PCP proteins in migrating cells and their impact on the cell biology of collective and individual cell migratory behaviors.

The basic Helix-Loop-Helix transcription factor Atoh1 is required for the specification of multiple neuron types in the mammalian hindbrain including tegmental, precerebellar output neurons and cerebellar granule neurons. How a single proneural gene specifies so many neuron types from a single progenitor zone, the upper rhombic lip (URL), is not known. Here we use the zebrafish to explore the role ofatoh1in cerebellar neurogenesis. Using transgenic reporters we show that zebrafishatoh1c-expressing cells give rise to tegmental and granule cell populations that, together with previously describedatoh1a-derived neuron populations, resemble the diversity ofatoh1derivatives observed in mammals. Using genetic mutants we find that of the three zebrafishatoh1paralogs,atoh1candatoh1aare required for the full complement of granule neurons in the zebrafish cerebellum. Interestingly,atoh1aandatoh1cspecify non-overlapping granule populations, indicating that fish use multipleatoh1genes to generate granule neuron diversity that is not detected in mammals. By live imaging of neurogenesis at the URL we show thatatoh1cis not required for the specification of granule neuron progenitors but promotes their delamination from the URL epithelium and this process is critical for neuronal maturation. This study thus provides a better understanding of how proneural transcription factors regulate neurogenesis in the vertebrate cerebellum.Summary statementatoh1genes specify distinct populations of tegmental and granule neurons in the zebrafish hindbrain and promote their delamination from the neuroepithelium, a process critical for neuronal maturation.

Forward genetic screens have illuminated how genes encode the information necessary for life (Crick et al. 1961; Brenner 1974; Nusslein-Volhard and Wieschaus 1980; Meyerowitz and Pruitt 1985; Haffter et al. 1996; Nolan et al. 2000). However, the subsequent identification of the causal mutation has been the bottleneck in the forward genetics process. Indeed, only one-third of the mutants identified in the first large-scale forward screens undertaken in a vertebrate model (Haffter et al. 1996) have been cloned. This problem has been solved in invertebrate systems through the use of whole-genome sequencing (WGS) of mutant animals to identify candidate genes; the main advantage of both Caenorhabditis elegans and Drosophila is that the genomes are small, animals are isogenic, and chemically induced mutations are rare enough that the changes can be identified and quickly confirmed to be causative (Sarin et al. 2008; Blumenstiel et al. 2009). In the zebrafish, the genome is large (greater than 10 times larger than worm or fly), making it relatively expensive to use WGS. Additionally, the zebrafish genome is highly polymorphic, with each strain, and even each individual, carrying numerous polymorphisms. Thus sequencing a single animal is not sufficient to distinguish potential causative mutations from other polymorphisms. Here we describe an RNA-seq–based bulk segregant analysis (BSA) approach that allows for the inexpensive mapping and identification of candidate mutations from forward genetic screens. While we have used zebrafish as a model, this methodology is applicable to any model system with a sequenced genome. BSA identifies regions of the genome that are linked to a causative mutation in a group of phenotypically mutant animals. This is accomplished by identifying regions of homozygosity within mutants at genetic markers found throughout the genome (Supplemental Fig. 1). BSA using PCR-based testing of infrequent microsatellite markers has been the standard for the initial mapping of zebrafish mutations (Geisler et al. 2007; Zhou and Zon 2011). The approach is both laborious and low-resolution and requires subsequent fine mapping using several hundreds to thousands of individual animals, each requiring multiple rounds of testing (Zhou and Zon 2011). This process is costly in terms of reagents and time. Next-generation sequencing (NGS) provides a means to identify the most abundant class of marker in the genome, single nucleotide polymorphisms (SNPs), in order to map mutations; additionally, in the same experiment, the data identify candidate mutations within the region of linkage that may be causal for the phenotype. Several WGS approaches have been successfully applied to the mapping and identification of candidate zebrafish mutations (Bowen et al. 2012; Leshchiner et al. 2012; Obholzer et al. 2012; Voz et al. 2012); however, because of the size of the zebrafish genome and the cost of WGS, these approaches have relied on relatively low coverage sequencing (two- to sevenfold). RNA-seq offers an alternative method to perform NGS mapping and presents several advantages over WGS: (1) It effectively reduces the representation of the genome, thereby decreasing the amount of sequencing needed to obtain high coverage and thus high-quality information; (2) the effect of candidate mutations on transcript splicing can be directly assessed in mutants; and (3) the effect of mutations on the expression levels of genes can be directly identified. Overall, the RNA-seq approach offers a number of advantages at reduced cost. Here we have developed and validated in vivo methodology and an in silico bioinformatics pipeline using RNA-seq to map and identify mutations. The in vivo preparation is simple and straightforward, and the bioinformatics pipeline is constructed from existing open-source programs and custom scripts. We have validated our approach on several independent mutations, demonstrating the ability to map and identify mutations that are deleterious due to amino acid changes, altered splicing, or altered expression levels. Importantly, we developed a simple computational platform that allows data processing within one integrated application. Our methodology greatly increases the power of forward genetics approaches in model systems with large, polymorphic genomes.

MYC-nick is a cytoplasmic, transcriptionally inactive member of the MYC oncoprotein family, generated by a proteolytic cleavage of full-length MYC. MYC-nick promotes migration and survival of cells in response to chemotherapeutic agents or withdrawal of glucose. Here we report that MYC-nick is abundant in colonic and intestinal tumors derived from mouse models with mutations in the Wnt, TGF-β, and PI3K pathways. Moreover, MYC-nick is elevated in colon cancer cells deleted for FBWX7, which encodes the major E3 ligase of full-length MYC frequently mutated in colorectal cancers. MYC-nick promotes the migration of colon cancer cells assayed in 3D cultures or grown as xenografts in a zebrafish metastasis model. MYC-nick accelerates migration by activating the Rho GTPase Cdc42 and inducing fascin expression. MYC-nick, fascin, and Cdc42 are frequently up-regulated in cells present at the invasive front of human colorectal tumors, suggesting a coordinated role for these proteins in tumor migration.

Hair cells possess a single primary cilium, called the kinocilium, early in development. While the kinocilium is lost in auditory hair cells of most species it is maintained in vestibular hair cells. It has generally been believed that the primary role of the kinocilium and cilia-associated genes in hair cells is in the establishment of the polarity of actin-based stereocilia, the hair cell mechanotransduction apparatus. Through genetic screening and testing of candidate genes in zebrafish (Danio rerio) we have found that mutations in multiple cilia genes implicated in intraflagellar transport (dync2h1, wdr35, ift88, and traf3ip), and the ciliary transition zone (cc2d2a, mks1, and cep290) lead to resistance to aminoglycoside-induced hair cell death. These genes appear to have differing roles in hair cells, as mutations in intraflagellar transport genes, but not transition zone genes, lead to defects in kinocilia formation and processes dependent upon hair cell mechanotransduction activity. These mutants highlight a novel role of cilia-associated genes in hair cells, and provide powerful tools for further study.

The mechanisms that allow cells to bypass anti–vascular endothelial growth factor A (VEGFA) therapy remain poorly understood. Here we use zebrafish to investigate this question and first show that vegfaa mutants display a severe vascular phenotype that can surprisingly be rescued to viability by vegfaa messenger RNA injections at the 1-cell stage. Using vegfaa mutants as an in vivo test tube, we found that zebrafish Vegfbb, Vegfd, and Pgfb can also rescue these animals to viability. Taking advantage of a new vegfr1 tyrosine kinase–deficient mutant, we determined that Pgfb rescues vegfaa mutants via Vegfr1. Altogether, these data reveal potential resistance routes against current anti-VEGFA therapies. In order to circumvent this resistance, we engineered and validated new dominant negative Vegfa molecules that by trapping Vegf family members can block vascular development. Thus, our results show that Vegfbb, Vegfd, and Pgfb can sustain vascular development in the absence of VegfA, and our newly engineered Vegf molecules expand the toolbox for basic research and antiangiogenic therapy.

Somites are formed from the presomitic mesoderm (PSM) and give rise to the axial skeleton and skeletal muscles. The PSM is dynamic; somites are generated at the anterior end, while the posterior end is continually renewed with new cells entering from the tailbud progenitor region. Which genes control the conversion of tailbud progenitors into PSM and how is this process coordinated with cell movement? Using loss- and gain-of-function experiments and heat-shock transgenics we show in zebrafish that the transcription factor Mesogenin 1 (Msgn1), acting with Spadetail (Spt), has a central role. Msgn1 allows progression of the PSM differentiation program by switching off the progenitor maintenance genes ntl, wnt3a, wnt8 and fgf8 in the future PSM cells as they exit from the tailbud, and subsequently induces expression of PSM markers such as tbx24. msgn1 is itself positively regulated by Ntl/Wnt/Fgf, creating a negative-feedback loop that might be crucial to regulate homeostasis of the progenitor population until somitogenesis ends. Msgn1 drives not only the changes in gene expression in the nascent PSM cells but also the movements by which they stream out of the tailbud into the PSM. Loss of Msgn1 reduces the flux of cells out of the tailbud, producing smaller somites and an enlarged tailbud, and, by delaying exhaustion of the progenitor population, results in supernumerary tail somites. Through its combined effects on gene expression and cell movement, Msgn1 (with Spt) plays a key role both in genesis of the paraxial mesoderm and in maintenance of the progenitor population from which it derives.

Myogenic regulatory factors of the myod family (MRFs) are transcription factors essential for mammalian skeletal myogenesis. Here we show that a mutation in the zebrafish myod gene delays and reduces early somitic and pectoral fin myogenesis, reduces miR-206 expression, and leads to a persistent reduction in somite size until at least the independent feeding stage. A mutation in myog, encoding a second MRF, has little obvious phenotype at early stages, but exacerbates the loss of somitic muscle caused by lack of Myod. Mutation of both myod and myf5 ablates all skeletal muscle. Haploinsufficiency of myod leads to reduced embryonic somite muscle bulk. Lack of Myod causes a severe reduction in cranial musculature, ablating most muscles including the protractor pectoralis, a putative cucullaris homologue. This phenotype is accompanied by a severe dysmorphology of the cartilaginous skeleton and failure of maturation of several cranial bones, including the opercle. As myod expression is restricted to myogenic cells, the data show that myogenesis is essential for proper skeletogenesis in the head.

Disc1 is a schizophrenia risk gene that engages multiple signaling pathways during neurogenesis and brain development. Using the zebrafish as a tool, we analyze the function of zebrafish Disc1 (zDisc1) at the earliest stages of brain and body development. We define a “tool” as a biological system that gives insight into mechanisms underlying a human disorder, although the system does not phenocopy the disorder. A zDisc1 peptide binds to GSK3β, and zDisc1 directs early brain development and neurogenesis, by promoting β-catenin-mediated Wnt signaling and inhibiting GSK3β activity. zDisc1 loss-of-function embryos additionally display a convergence and extension phenotype, demonstrated by abnormal movement of dorsolateral cells during gastrulation, through changes in gene expression, and later through formation of abnormal, U-shaped muscle segments, and a truncated tail. These phenotypes are caused by alterations in the noncanonical Wnt pathway, via Daam and Rho signaling. The convergence and extension phenotype can be rescued by a dominant negative GSK3β construct, suggesting that zDisc1 inhibits GSK3β activity during noncanonical Wnt signaling. This is the first demonstration that Disc1 modulates the noncanonical Wnt pathway and suggests a previously unconsidered mechanism by which Disc1 may contribute to the etiology of neuropsychiatric disorders.—De Rienzo, G., Bishop, J. A., Mao, Y., Pan, L., Ma, T. P., Moens, C. B., Tsai, L. H., Sive, H. Disc1 regulates both β-catenin-mediated and noncanonical Wnt signaling during vertebrate embryogenesis.

Piwi proteins function in an RNAi-like pathway that silences transposons. Piwi-associated RNAs, also known as piRNAs, act as a guide to identify Piwi targets. The tudor domain-containing protein Tdrd1 has been linked to this pathway but its function has thus far remained unclear. We show that zebrafish Tdrd1 is required for efficient Piwi-pathway activity and proper nuage formation. Furthermore, we find that Tdrd1 binds both zebrafish Piwi proteins, Ziwi and Zili, and reveals sequence specificity in the interaction between Tdrd1 tudor domains and symmetrically dimethylated arginines (sDMAs) in Zili. Finally, we show that Tdrd1 complexes contain piRNAs and RNA molecules that are longer than piRNAs. We name these longer transcripts Tdrd1-associated transcripts (TATs). TATs likely represent cleaved Piwi pathway targets and may serve as piRNA biogenesis intermediates. Altogether, our data suggest that Tdrd1 acts as a molecular scaffold for Piwi proteins, bound through specific tudor domain-sDMA interactions, piRNAs and piRNA targets.

The facial branchiomotor neurons (FBMNs) undergo a characteristic tangential migration in the vertebrate hindbrain. We previously used a morpholino knockdown approach to reveal that zebrafish prickle1b (pk1b) is required for this migration. Here we report that FBMN migration is also blocked in a pk1b mutant with a disruption in the consensus farnesylation motif. We confirmed that this lipid modification is required during FBMN migration by disrupting the function of farnesyl biosynthetic enzymes. Furthermore, farnesylation of a tagged Pk1b is required for its nuclear localization. Using a unique rescue approach, we have demonstrated that Pk1b nuclear localization and farnesylation are required during FBMN migration. Our data suggest that Pk1b acts at least partially independently of core planar cell polarity molecules at the plasma membrane, and might instead be acting at the nucleus. We also found that the neuronal transcriptional silencer REST is necessary for FBMN migration, and we provide evidence that interaction between Pk1b and REST is required during this process. Finally, we demonstrate that REST protein, which is normally localized in the nuclei of migrating FBMNs, is depleted from the nuclei of Pk1b-deficient neurons. We conclude that farnesylation-dependent nuclear localization of Pk1b is required to regulate REST localization and thus FBMN migration.

How control of subcellular events in single cells determines morphogenesis on the scale of the tissue is largely unresolved. The stereotyped cross-midline mitoses of progenitors in the zebrafish neural keel [1–4] provide a unique experimental paradigm for defining the role and control of single-cell orientation for tissue-level morphogenesis in vivo. We show here that the coordinated orientation of individual progenitor cell division in the neural keel is the cellular determinant required for morphogenesis into a neural tube epithelium with a single straight lumen. We find that Scribble is required for oriented cell division and that its function in this process is independent of canonical apicobasal and planar polarity pathways. We identify a role for Scribble in controlling clustering of α-catenin foci in dividing progenitors. Loss of either Scrib or N-cadherin results in abnormally oriented mitoses, reduced cross-midline cell divisions, and similar neural tube defects. We propose that Scribble-dependent nascent cell-cell adhesion clusters between neuroepithelial progenitors contribute to define orientation of their cell division. Finally, our data demonstrate that while oriented mitoses of individual cells determine neural tube architecture, the tissue can in turn feed back on its constituent cells to define their polarization and cell division orientation to ensure robust tissue morphogenesis.

Retinoic acid (RA) is an important developmental signaling molecule responsible for the patterning of multiple vertebrate tissues. RA is also a potent teratogen, causing multi-organ birth defects in humans. Endogenous RA levels must therefore be tightly controlled in the developing embryo. We used a microarray approach to identify genes that function as negative feedback regulators of retinoic acid signaling. We screened for genes expressed in early somite-stage embryos that respond oppositely to treatment with RA versus RA antagonists and validated them by RNA in situ hybridization. Focusing on genes known to be involved in RA metabolism, we determined that dhrs3a, which encodes a member of the short-chain dehydrogenase/reductase protein family, is both RA dependent and strongly RA inducible. Dhrs3a is known to catalyze the reduction of the RA precursor all-trans retinaldehyde to vitamin A; however, a developmental function has not been demonstrated. Using morpholino knockdown and mRNA over-expression, we demonstrate that Dhrs3a is required to limit RA levels in the embryo, primarily within the central nervous system. Dhrs3a is thus an RA-induced feedback inhibitor of RA biosynthesis. We conclude that retinaldehyde availability is an important level at which RA biosynthesis is regulated in vertebrate embryos.

Spinal muscular atrophy (SMA), a recessive genetic disease, affects lower motoneurons leading to denervation, atrophy, paralysis and in severe cases death. Reduced levels of survival motor neuron (SMN) protein cause SMA. As a first step towards generating a genetic model of SMA in zebrafish, we identified three smn mutations. Two of these alleles, smnY262stop and smnL265stop, were stop mutations that resulted in exon 7 truncation, whereas the third, smnG264D, was a missense mutation corresponding to an amino acid altered in human SMA patients. Smn protein levels were low/undetectable in homozygous mutants consistent with unstable protein products. Homozygous mutants from all three alleles were smaller and survived on the basis of maternal Smn dying during the second week of larval development. Analysis of the neuromuscular system in these mutants revealed a decrease in the synaptic vesicle protein, SV2. However, two other synaptic vesicle proteins, synaptotagmin and synaptophysin were unaffected. To address whether the SV2 decrease was due specifically to Smn in motoneurons, we tested whether expressing human SMN protein exclusively in motoneurons in smn mutants could rescue the phenotype. For this, we generated a transgenic zebrafish line with human SMN driven by the motoneuron-specific zebrafish hb9 promoter and then generated smn mutant lines carrying this transgene. We found that introducing human SMN specifically into motoneurons rescued the SV2 decrease observed in smn mutants. Our analysis indicates the requirement for Smn in motoneurons to maintain SV2 in presynaptic terminals indicating that Smn, either directly or indirectly, plays a role in presynaptic integrity.

During vertebrate development, the hindbrain is transiently segmented into 7 distinct rhombomeres (r). Hindbrain segmentation takes place within the context of the complex morphogenesis required for neurulation, which in zebrafish involves a characteristic cross-midline division that distributes progenitor cells bilaterally in the forming neural tube. The Eph receptor tyrosine kinase EphA4 and the membrane-bound Ephrin (Efn) ligand EfnB2a, which are expressed in complementary segments in the early hindbrain, are required for rhombomere boundary formation. We showed previously that EphA4 promotes cell–cell affinity within r3 and r5, and proposed that preferential adhesion within rhombomeres contributes to boundary formation. Here we show that EfnB2a is similarly required in r4 for normal cell affinity and that EphA4 and EfnB2a regulate cell affinity independently within their respective rhombomeres. Live imaging of cell sorting in mosaic embryos shows that both proteins function during cross-midline cell divisions in the hindbrain neural keel. Consistent with this, mosaic EfnB2a over-expression causes widespread cell sorting and disrupts hindbrain organization, but only if induced at or before neural keel stage. We propose a model in which Eph and Efn-dependent cell affinity within rhombomeres serve to maintain rhombomere organization during the potentially disruptive process of teleost neurulation.

TILLING, for Targeting Induced Local Lesions in Genomes, is a reverse genetics strategy that identifies mutations in specific genes of interest in chemically mutagenized populations. First described in 2000 for mutation detection in Arabidopsis, TILLING is now used in a wide range of plants including soybean, rice, barley and maize as well as for animal model systems, including Arabidopsis, Drosophila, Caenorhabditis elegans, rat, medaka and zebrafish and for the discovery of naturally occurring polymorphisms in humans. This review summarizes current TILLING methodologies as they have been applied to the zebrafish, ongoing TILLING projects and resources in the zebrafish community, and the future of zebrafish TILLING.

In the postgenomic era, the ability to quickly, efficiently, and inexpensively assign function to the zebrafish proteome is critical. Clustered regularly interspaced short palindromic repeats (CRISPRs) have revolutionized the ability to perform reverse genetics because of its simplicity and broad applicability. The CRISPR system is composed of an engineered, gene-specific single guide RNA (sgRNA) and a Cas9 enzyme that causes double-stranded breaks in DNA at the targeted site. This simple, two-part system, when injected into one-cell stage zebrafish embryos, efficiently mutates target loci at a frequency such that injected embryos phenocopy known mutant phenotypes. This property allows for CRISPR-based F0 screening in zebrafish, which provides a means to screen through a large number of candidate genes for their role in a phenotype of interest. While there are important considerations for any successful genetic screen, CRISPR screening has significant benefits over conventional methods and can be accomplished in any lab with modest molecular biology experience.

The basic helix-loop-helix (bHLH) transcription factor Myod directly regulates gene expression throughout the program of skeletal muscle differentiation. It is not known how a Myod-driven myogenic program is modulated to achieve muscle fiber-type-specific gene expression. Pbx homeodomain proteins mark promoters of a subset of Myod target genes,including myogenin (Myog); thus, Pbx proteins might modulate the program of myogenesis driven by Myod. By inhibiting Pbx function in zebrafish embryos, we show that Pbx proteins are required in order for Myod to induce the expression of a subset of muscle genes in the somites. In the absence of Pbx function, expression of myog and of fast-muscle genes is inhibited, whereas slow-muscle gene expression appears normal. By knocking down Pbx or Myod function in combination with another bHLH myogenic factor,Myf5, we show that Pbx is required for Myod to regulate fast-muscle, but not slow-muscle, development. Furthermore, we show that Sonic hedgehog requires Myod in order to induce both fast- and slow-muscle markers but requires Pbx only to induce fast-muscle markers. Our results reveal that Pbx proteins modulate Myod activity to drive fast-muscle gene expression, thus showing that homeodomain proteins can direct bHLH proteins to establish a specific cell-type identity.

Piwi proteins specify an animal-specific subclass of the Argonaute family that, in vertebrates, is specifically expressed in germ cells. We demonstrate that zebrafish Piwi (Ziwi) is expressed in both the male and the female gonad and is a component of a germline-specifying structure called nuage. Loss of Ziwi function results in a progressive loss of germ cells due to apoptosis during larval development. In animals that have reduced Ziwi function, germ cells are maintained but display abnormal levels of apoptosis in adults. In mammals, Piwi proteins associate with approximately 29-nucleotide-long, testis-specific RNA molecules called piRNAs. Here we show that zebrafish piRNAs are present in both ovary and testis. Many of these are derived from transposons, implicating a role for piRNAs in the silencing of repetitive elements in vertebrates. Furthermore, we show that piRNAs are Dicer independent and that their 3′ end likely carries a 2′O-Methyl modification.

The medial habenular nuclei of the zebrafish diencephalon, which lie bilateral to the pineal complex, exhibit left-right differences in their neuroanatomy, gene expression profiles and axonal projections to the unpaired midbrain target - the interpeduncular nucleus (IPN). Efferents from the left habenula terminate along the entire dorsoventral extent of the IPN, whereas axons from the right habenula project only to the ventral IPN. How this left-right difference in connectivity is established and the factors involved in differential target recognition are unknown. Prior to IPN innervation, we find that only the left habenula expresses the zebrafish homologue of Neuropilin1a (Nrp1a), a receptor for class III Semaphorins (Sema3s). Directional asymmetry of nrp1a expression relies on Nodal signaling and the presence of the left-sided parapineal organ. Loss of Nrp1a, through parapineal ablation or depletion by antisense morpholinos, prevents left habenular neurons from projecting to the dorsal IPN. Selective depletion of Sema3D, but not of other Sema family members, similarly disrupts innervation of the dorsal IPN. Conversely, Sema3D overexpression results in left habenular projections that extend to the dorsal IPN, as well as beyond the target. The results indicate that Sema3D acts in concert with Nrp1a to guide neurons on the left side of the brain to innervate the target nucleus differently than those on the right side.