Embryo-scale reverse genetics at single-cell resolution.

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

The role of retrograde intraflagellar transport genes in aminoglycoside-induced hair cell death.

Sensory hair cells are susceptible to numerous insults, including certain therapeutic medications like aminoglycoside antibiotics, and hearing and balance disorders are often a dose-limiting side effect of these medications. We show that mutations in multiple genes in both the retrograde intraflagellar transport (IFT) motor and adaptor complexes lead to resistance to aminoglycoside-induced hair cell death. These mutations also lead to defects in the entry of both aminoglycosides and the vital dye FM1-43 into hair cells, both processes that depend on hair cell mechanotransduction activity. However, the trafficking of proteins important for mechanotransduction activity is not altered by these mutations. Our data suggest that both retrograde IFT motor and adaptor complex genes are playing a role in aminoglycoside toxicity through affecting aminoglycoside uptake into hair cells.

Cumulative mitochondrial activity correlates with ototoxin susceptibility in zebrafish mechanosensory hair cells.

Mitochondria play a prominent role in mechanosensory hair cell damage and death. Although hair cells are thought to be energetically demanding cells, how mitochondria respond to these demands and how this might relate to cell death is largely unexplored. Using genetically encoded indicators, we found that mitochondrial calcium flux and oxidation are regulated by mechanotransduction and demonstrate that hair cell activity has both acute and long-term consequences on mitochondrial function. We tested whether variation in mitochondrial activity reflected differences in the vulnerability of hair cells to the toxic drug neomycin. We observed that susceptibility did not correspond to the acute level of mitochondrial activity but rather to the cumulative history of that activity.

The occhiolino (occ) mutant Zebrafish, a model for development of the optical function in the biological lens.

BACKGROUND: Zebrafish visual function depends on quality optics. An F3 screen for developmental mutations in the Zebrafish nervous system was conducted in wild-type (wt) AB Zebrafish exposed to 3 mM of N-ethyl-N-nitrosourea (ENU). RESULTS: Mutant offspring, identified in an F3 screen, were characterized by a small pupil, resulting from retinal hypertrophy or hyperplasia and a small lens. Deficits in visual function made feeding difficult after hatching at approximately 5-6 days postfertilization (dpf). Special feeding conditions were necessary for survival of the occhiolino (occ) mutants after 6 dpf. Optokinetic response (OKR) tests measured defects in visual function in the occ mutant, although electroretinograms (ERGs) were normal in the mutant and wt. Consistent with the ERGs, histology found normal retinal structure in the occ mutant and wt Zebrafish. However, lens development was abnormal. Multiphoton imaging of the developmental stages of live embryos confirmed the formation of a secondary mass of lens cells in the developing eye of the mutant Zebrafish at 3-4 dpf, and laminin immunohistochemistry indicated the lens capsule was thin and disorganized in the mutant Zebrafish. CONCLUSIONS: The occ Zebrafish is a novel disease model for visual defects associated with abnormal lens development. Developmental Dynamics 246:915-924, 2017. © 2017 Wiley Periodicals, Inc.

Fluorescent aminoglycosides reveal intracellular trafficking routes in mechanosensory hair cells.

Aminoglycosides (AGs) are broad-spectrum antibiotics that are associated with kidney damage, balance disorders, and permanent hearing loss. This damage occurs primarily by killing of proximal tubule kidney cells and mechanosensory hair cells, though the mechanisms underlying cell death are not clear. Imaging molecules of interest in living cells can elucidate how molecules enter cells, traverse intracellular compartments, and interact with sites of activity. Here, we have imaged fluorescently labeled AGs in live zebrafish mechanosensory hair cells. We determined that AGs enter hair cells via both nonendocytic and endocytic pathways. Both routes deliver AGs from the extracellular space to lysosomes, and structural differences between AGs alter the efficiency of this delivery. AGs with slower delivery to lysosomes were immediately toxic to hair cells, and impeding lysosome delivery increased AG-induced death. Therefore, pro-death cascades induced at early time points of AG exposure do not appear to derive from the lysosome. Our findings help clarify how AGs induce hair cell death and reveal properties that predict toxicity. Establishing signatures for AG toxicity may enable more efficient evaluation of AG treatment paradigms and structural modifications to reduce hair cell damage. Further, this work demonstrates how following fluorescently labeled drugs at high resolution in living cells can reveal important details about how drugs of interest behave.

Mitochondrial calcium uptake underlies ROS generation during aminoglycoside-induced hair cell death.

Exposure to aminoglycoside antibiotics can lead to the generation of toxic levels of reactive oxygen species (ROS) within mechanosensory hair cells of the inner ear that have been implicated in hearing and balance disorders. Better understanding of the origin of aminoglycoside-induced ROS could focus the development of therapies aimed at preventing this event. In this work, we used the zebrafish lateral line system to monitor the dynamic behavior of mitochondrial and cytoplasmic oxidation occurring within the same dying hair cell following exposure to aminoglycosides. The increased oxidation observed in both mitochondria and cytoplasm of dying hair cells was highly correlated with mitochondrial calcium uptake. Application of the mitochondrial uniporter inhibitor Ru360 reduced mitochondrial and cytoplasmic oxidation, suggesting that mitochondrial calcium drives ROS generation during aminoglycoside-induced hair cell death. Furthermore, targeting mitochondria with free radical scavengers conferred superior protection against aminoglycoside exposure compared with identical, untargeted scavengers. Our findings suggest that targeted therapies aimed at preventing mitochondrial oxidation have therapeutic potential to ameliorate the toxic effects of aminoglycoside exposure.

Cilia-Associated Genes Play Differing Roles in Aminoglycoside-Induced Hair Cell Death in Zebrafish.

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 zebrafish merovingian mutant reveals a role for pH regulation in hair cell toxicity and function.

Control of the extracellular environment of inner ear hair cells by ionic transporters is crucial for hair cell function. In addition to inner ear hair cells, aquatic vertebrates have hair cells on the surface of their body in the lateral line system. The ionic environment of these cells also appears to be regulated, although the mechanisms of this regulation are less understood than those of the mammalian inner ear. We identified the merovingian mutant through genetic screening in zebrafish for genes involved in drug-induced hair cell death. Mutants show complete resistance to neomycin-induced hair cell death and partial resistance to cisplatin-induced hair cell death. This resistance is probably due to impaired drug uptake as a result of reduced mechanotransduction ability, suggesting that the mutants have defects in hair cell function independent of drug treatment. Through genetic mapping we found that merovingian mutants contain a mutation in the transcription factor gcm2. This gene is important for the production of ionocytes, which are cells crucial for whole body pH regulation in fish. We found that merovingian mutants showed an acidified extracellular environment in the vicinity of both inner ear and lateral line hair cells. We believe that this acidified extracellular environment is responsible for the defects seen in hair cells of merovingian mutants, and that these mutants would serve as a valuable model for further study of the role of pH in hair cell function.

oca2 Regulation of chromatophore differentiation and number is cell type specific in zebrafish.

We characterized a zebrafish mutant that displays defects in melanin synthesis and in the differentiation of melanophores and iridophores of the skin and retinal pigment epithelium. Positional cloning and candidate gene sequencing link this mutation to a 410-kb region on chromosome 6, containing the oculocutaneous albinism 2 (oca2) gene. Quantification of oca2 mutant melanophores shows a reduction in the number of differentiated melanophores compared with wildtype siblings. Consistent with the analysis of mouse Oca2-deficient melanocytes, zebrafish mutant melanophores have immature melanosomes which are partially rescued following treatment with vacuolar-type ATPase inhibitor/cytoplasmic pH modifier, bafilomycin A1. Melanophore-specific gene expression is detected at the correct time and in anticipated locations. While oca2 zebrafish display unpigmented gaps on the head region of mutants 3 days post-fertilization, melanoblast quantification indicates that oca2 mutants have the correct number of melanoblasts, suggesting a differentiation defect explains the reduced melanophore number. Unlike melanophores, which are reduced in number in oca2 mutants, differentiated iridophores are present at significantly higher numbers. These data suggest distinct mechanisms for oca2 in establishing differentiated chromatophore number in developing zebrafish.

Modulation of dorsal root ganglion development by ErbB signaling and the scaffold protein Sorbs3.

The multipotent cells of the vertebrate neural crest (NC) arise at the dorsal aspect of the neural tube, then migrate throughout the developing embryo and differentiate into diverse cell types, including the sensory neurons and glia of the dorsal root ganglia (DRG). As multiple cell types are derived from this lineage, it is ideal for examining mechanisms of fate restriction during development. We have isolated a mutant, ouchless, that specifically fails to develop DRG neurons, although other NC derivatives develop normally. This mutation affects the expression of Sorbs3, a scaffold protein known to interact with proteins involved in focal adhesions and several signaling pathways. ouchless mutants share some phenotypic similarities with mutants in ErbB receptors, EGFR homologs that are implicated in diverse developmental processes and associated with several cancers; and ouchless interacts genetically with an allele of erbb3 in DRG neurogenesis. However, the defect in ouchless DRG neurogenesis is distinct from ErbB loss of function in that it is not associated with a loss of glia. Both ouchless and neurogenin1 heterozygous fish are sensitized to the effects of ErbB chemical inhibitors, which block the development of DRG in a dose-dependent manner. Inhibitors of MEK show similar effects on DRG neurogenesis. We propose a model in which Sorbs3 helps to integrate ErbB signals to promote DRG neurogenesis through the activation of MAPK and upregulation of neurogenin1.

Maintenance of melanophore morphology and survival is cathepsin and vps11 dependent in zebrafish.

Here, we characterize a Danio rerio zebrafish pigment cell mutant (melanophore integrity mutant), which displays a defect in maintenance of melanophore and iridophore number. Mapping and candidate gene analysis links the melanophore integrity mutant mutation to the vacuolar protein sorting 11 (vps11(w66)) gene. Quantification of vps11(w66) chromatophores during larval stages suggests a decrease in number as compared to wildtype siblings. TUNEL analysis and treatment with the caspase inhibitor, zVAD-fmk, indicate that vps11(w66) chromatophore death is caspase independent. Western blot analysis of PARP-1 cleavage patterns in mutant lysates suggests that increases in pH dependent cathepsin activity is involved in the premature chromatophore death observed in vps11(w66) mutants. Consistently, treatment with ALLM and Bafilomycin A1 (cathepsin/calpain and vacuolar-type H+-ATPase inhibitors, respectively), restore normal melanophore morphology and number in vps11(w66) mutants. Last, LC3B western blot analysis indicates an increase in autophagosome marker, LC3B II in vps11(w66) mutants as compared to wildtype control, but not in ALLM or Bafilomycin A1 treated mutants. Taken together, these data suggest that vps11 promotes normal melanophore morphology and survival by inhibiting cathepsin release and/or activity.

Loss of Slc4a1b chloride/bicarbonate exchanger function protects mechanosensory hair cells from aminoglycoside damage in the zebrafish mutant persephone.

Mechanosensory hair cell death is a leading cause of hearing and balance disorders in the human population. Hair cells are remarkably sensitive to environmental insults such as excessive noise and exposure to some otherwise therapeutic drugs. However, individual responses to damaging agents can vary, in part due to genetic differences. We previously carried out a forward genetic screen using the zebrafish lateral line system to identify mutations that alter the response of larval hair cells to the antibiotic neomycin, one of a class of aminoglycoside compounds that cause hair cell death in humans. The persephone mutation confers resistance to aminoglycosides. 5 dpf homozygous persephone mutants are indistinguishable from wild-type siblings, but differ in their retention of lateral line hair cells upon exposure to neomycin. The mutation in persephone maps to the chloride/bicarbonate exchanger slc4a1b and introduces a single Ser-to-Phe substitution in zSlc4a1b. This mutation prevents delivery of the exchanger to the cell surface and abolishes the ability of the protein to import chloride across the plasma membrane. Loss of function of zSlc4a1b reduces hair cell death caused by exposure to the aminoglycosides neomycin, kanamycin, and gentamicin, and the chemotherapeutic drug cisplatin. Pharmacological block of anion transport with the disulfonic stilbene derivatives DIDS and SITS, or exposure to exogenous bicarbonate, also protects hair cells against damage. Both persephone mutant and DIDS-treated wild-type larvae show reduced uptake of labeled aminoglycosides. persephone mutants also show reduced FM1-43 uptake, indicating a potential impact on mechanotransduction-coupled activity in the mutant. We propose that tight regulation of the ionic environment of sensory hair cells, mediated by zSlc4a1b activity, is critical for their sensitivity to aminoglycoside antibiotics.

The metalloproteinase inhibitor Reck is essential for zebrafish DRG development.

The neural crest is a migratory, multipotent cell lineage that contributes to myriad tissues, including sensory neurons and glia of the dorsal root ganglia (DRG). To identify genes affecting cell fate specification in neural crest, we performed a forward genetic screen for mutations causing DRG deficiencies in zebrafish. This screen yielded a mutant lacking all DRG, which we named sensory deprived (sdp). We identified a total of four alleles of sdp, all of which possess lesions in the gene coding for reversion-inducing cysteine-rich protein containing Kazal motifs (Reck). Reck is an inhibitor of metalloproteinases previously shown to regulate cell motility. We found reck function to be both necessary for DRG formation and sufficient to rescue the sdp phenotype. reck is expressed in neural crest cells and is required in a cell-autonomous fashion for appropriate sensory neuron formation. In the absence of reck function, sensory neuron precursors fail to migrate to the position of the DRG, suggesting that this molecule is crucial for proper migration and differentiation.

Lef1 is required for progenitor cell identity in the zebrafish lateral line primordium.

The zebrafish posterior lateral line (pLL) is a sensory system that comprises clusters of mechanosensory organs called neuromasts (NMs) that are stereotypically positioned along the surface of the trunk. The NMs are deposited by a migrating pLL primordium, which is organized into polarized rosettes (proto-NMs). During migration, mature proto-NMs are deposited from the trailing part of the primordium, while progenitor cells in the leading part give rise to new proto-NMs. Wnt signaling is active in the leading zone of the primordium and global Wnt inactivation leads to dramatic disorganization of the primordium and a loss of proto-NM formation. However, the exact cellular events that are regulated by the Wnt pathway are not known. We identified a mutant strain, lef1(nl2), that contains a lesion in the Wnt effector gene lef1. lef1(nl2) mutants lack posterior NMs and live imaging reveals that rosette renewal fails during later stages of migration. Surprisingly, the overall primordium patterning, as assayed by the expression of various markers, appears unaltered in lef1(nl2) mutants. Lineage tracing and mosaic analyses revealed that the leading cells (presumptive progenitors) move out of the primordium and are incorporated into NMs; this results in a decrease in the number of proliferating progenitor cells and eventual primordium disorganization. We concluded that Lef1 function is not required for initial primordium organization or migration, but is necessary for proto-NM renewal during later stages of pLL formation. These findings revealed a novel role for the Wnt signaling pathway during mechanosensory organ formation in zebrafish.

UDP xylose synthase 1 is required for morphogenesis and histogenesis of the craniofacial skeleton.

UDP-xylose synthase (Uxs1) is strongly conserved from bacteria to humans, but because no mutation has been studied in any animal, we do not understand its roles in development. Furthermore, no crystal structure has been published. Uxs1 synthesizes UDP-xylose, which initiates glycosaminoglycan attachment to a protein core during proteoglycan formation. Crystal structure and biochemical analyses revealed that an R233H substitution mutation in zebrafish uxs1 alters an arginine buried in the dimer interface, thereby destabilizing and, as enzyme assays show, inactivating the enzyme. Homozygous uxs1 mutants lack Alcian blue-positive, proteoglycan-rich extracellular matrix in cartilages of the neurocranium, pharyngeal arches, and pectoral girdle. Transcripts for uxs1 localize to skeletal domains at hatching. GFP-labeled neural crest cells revealed defective organization and morphogenesis of chondrocytes, perichondrium, and bone in uxs1 mutants. Proteoglycans were dramatically reduced and defectively localized in uxs1 mutants. Although col2a1a transcripts over-accumulated in uxs1 mutants, diminished quantities of Col2a1 protein suggested a role for proteoglycans in collagen secretion or localization. Expression of col10a1, indian hedgehog, and patched was disrupted in mutants, reflecting improper chondrocyte/perichondrium signaling. Up-regulation of sox9a, sox9b, and runx2b in mutants suggested a molecular mechanism consistent with a role for proteoglycans in regulating skeletal cell fate. Together, our data reveal time-dependent changes to gene expression in uxs1 mutants that support a signaling role for proteoglycans during at least two distinct phases of skeletal development. These investigations are the first to examine the effect of mutation on the structure and function of Uxs1 protein in any vertebrate embryos, and reveal that Uxs1 activity is essential for the production and organization of skeletal extracellular matrix, with consequent effects on cartilage, perichondral, and bone morphogenesis.

Kit and foxd3 genetically interact to regulate melanophore survival in zebrafish.

We have investigated the role of foxd3 activity in conjunction with signaling by the kit tyrosine kinase receptor in zebrafish black pigment cell (melanophore) development. As loss-of-function of these molecules individually has distinct effects on melanophore number, we have examined the phenotype of double mutants. Individuals with a null mutation in kit have fewer melanophores than wild-type, with cells lost through death. When kit mutants are injected with foxd3 antisense morpholino oligonucleotides or crossed with a foxd3 zebrafish mutant, they have more melanophores than their uninjected or foxd3+ counterparts. Examination of foxd3 loss-of-function in two additional kit mutants that differentially alter kit-dependent migration and survival indicates a change in melanophore number in survival mutants only. Consistently, TUNEL (terminal deoxynucleotidyl transferase-mediated deoxyuridinetriphosphate nick end-labeling) analysis confirms a partial rescue of melanophores from cell death. Ectopic expression of foxd3 indicates that foxd3 promotes early melanophore death only when kit is inactive. Taken together, these data suggest a kit-dependent role for foxd3 in the regulation of melanophore survival.

Identification of genetic and chemical modulators of zebrafish mechanosensory hair cell death.

Inner ear sensory hair cell death is observed in the majority of hearing and balance disorders, affecting the health of more than 600 million people worldwide. While normal aging is the single greatest contributor, exposure to environmental toxins and therapeutic drugs such as aminoglycoside antibiotics and antineoplastic agents are significant contributors. Genetic variation contributes markedly to differences in normal disease progression during aging and in susceptibility to ototoxic agents. Using the lateral line system of larval zebrafish, we developed an in vivo drug toxicity interaction screen to uncover genetic modulators of antibiotic-induced hair cell death and to identify compounds that confer protection. We have identified 5 mutations that modulate aminoglycoside susceptibility. Further characterization and identification of one protective mutant, sentinel (snl), revealed a novel conserved vertebrate gene. A similar screen identified a new class of drug-like small molecules, benzothiophene carboxamides, that prevent aminoglycoside-induced hair cell death in zebrafish and in mammals. Testing for interaction with the sentinel mutation suggests that the gene and compounds may operate in different pathways. The combination of chemical screening with traditional genetic approaches is a new strategy for identifying drugs and drug targets to attenuate hearing and balance disorders.

Specification of epibranchial placodes in zebrafish.

In all vertebrates, the neurogenic placodes are transient ectodermal thickenings that give rise to sensory neurons of the cranial ganglia. Epibranchial (EB) placodes generate neurons of the distal facial, glossopharyngeal and vagal ganglia, which convey sensation from the viscera, including pharyngeal endoderm structures, to the CNS. Recent studies have implicated signals from pharyngeal endoderm in the initiation of neurogenesis from EB placodes; however, the signals underlying the formation of placodes are unknown. Here, we show that zebrafish embryos mutant for fgf3 and fgf8 do not express early EB placode markers, including foxi1 and pax2a. Mosaic analysis demonstrates that placodal cells must directly receive Fgf signals during a specific crucial period of development. Transplantation experiments and mutant analysis reveal that cephalic mesoderm is the source of Fgf signals. Finally, both Fgf3 and Fgf8 are sufficient to induce foxi1-positive placodal precursors in wild-type as well as Fgf3-plus Fgf8-depleted embryos. We propose a model in which mesoderm-derived Fgf3 and Fgf8 signals establish both the EB placodes and the development of the pharyngeal endoderm, the subsequent interaction of which promotes neurogenesis. The coordinated interplay between craniofacial tissues would thus assure proper spatial and temporal interactions in the shaping of the vertebrate head.

Endoderm-derived Fgf3 is necessary and sufficient for inducing neurogenesis in the epibranchial placodes in zebrafish.

In vertebrates, epibranchial placodes are transient ectodermal thickenings that contribute sensory neurons to the epibranchial ganglia. These ganglia innervate internal organs and transmit information on heart rate, blood pressure and visceral distension from the periphery to the central nervous system. Despite their importance, the molecular mechanisms that govern the induction and neurogenesis of the epibranchial placodes are only now being elucidated. In this study, we demonstrate that endoderm is required for neurogenesis of the zebrafish epibranchial placodes. Mosaic analyses confirm that endoderm is the source of the neurogenic signal. Using a morpholino knockdown approach, we find that fgf3 is required for the majority of placode cells to undergo neurogenesis. Tissue transplants demonstrate that fgf3 activity is specifically required in the endodermal pouches. Furthermore, ectopic fgf3 expression is sufficient for inducing phox2a-positive neurons in wild-type and endoderm-deficient embryos. Surprisingly, ectodermal foxi1 expression, a marker for the epibranchial placode precursors, is present in both endoderm-deficient embryos and fgf3 morphants, indicating that neither endoderm nor Fgf3 is required for initial placode induction. Based on these findings, we propose a model for epibranchial placode development in which Fgf3 is a major endodermal determinant required for epibranchial placode neurogenesis.