Microbe transmission from pet shop to lab-reared zebrafish reveals a pathogenic birnavirus.

Zebrafish are popular research organisms selected for laboratory use due in part to widespread availability from the pet trade. Many contemporary colonies of laboratory zebrafish are maintained in aquaculture facilities that monitor and aim to curb infections that can negatively affect colony health and confound experiments. The impact of laboratory control on the microbial constituents associated with zebrafish in research environments compared to the pet trade are unclear. Diseases of unknown causes are common in both environments. We conducted a metatranscriptomic survey to broadly compare the zebrafish-associated microbes in pet trade and laboratory environments. We detected many microbes in animals from the pet trade that were not found in laboratory animals. Cohousing experiments revealed several transmissible microbes including a newly described non-enveloped, double-stranded RNA virus in the Birnaviridae family we name Rocky Mountain birnavirus (RMBV). Infections were detected in asymptomatic animals from the pet trade, but when transmitted to laboratory animals RMBV was associated with pronounced antiviral responses and hemorrhagic disease. These experiments highlight the pet trade as a distinct source of diverse microbes that associate with zebrafish and establish a paradigm for the discovery of newly described pathogenic viruses and other infectious microbes that can be developed for study in the laboratory.

Distinct features of the regenerating heart uncovered through comparative single-cell profiling.

Adult humans respond to heart injury by forming a permanent scar, yet other vertebrates are capable of robust and complete cardiac regeneration. Despite progress towards characterizing the mechanisms of cardiac regeneration in fish and amphibians, the large evolutionary gulf between mammals and regenerating vertebrates complicates deciphering which cellular and molecular features truly enable regeneration. To better define these features, we compared cardiac injury responses in zebrafish and medaka, two fish species that share similar heart anatomy and common teleost ancestry but differ in regenerative capability. We used single-cell transcriptional profiling to create a time-resolved comparative cell atlas of injury responses in all major cardiac cell types across both species. With this approach, we identified several key features that distinguish cardiac injury response in the non-regenerating medaka heart. By comparing immune responses to injury, we found altered cell recruitment and a distinct pro-inflammatory gene program in medaka leukocytes, and an absence of the injury-induced interferon response seen in zebrafish. In addition, we found a lack of pro-regenerative signals, including nrg1 and retinoic acid, from medaka endothelial and epicardial cells. Finally, we identified alterations in the myocardial structure in medaka, where they lack primordial layer cardiomyocytes and fail to employ a cardioprotective gene program shared by regenerating vertebrates. Our findings reveal notable variation in injury response across nearly all major cardiac cell types in zebrafish and medaka, demonstrating how evolutionary divergence influences the hidden cellular features underpinning regenerative potential in these seemingly similar vertebrates.

Germ cells do not progress through spermatogenesis in the infertile zebrafish testis.

Vertebrate spermatogonial stem cells maintain sperm production over the lifetime of an animal but fertility declines with age. While morphological studies have greatly informed our understanding of typical spermatogenesis, the molecular and cellular mechanisms underlying spermatogenesis are not yet understood, particularly with respect to the onset of fertility. We used single-cell RNA sequencing to generate a developmental atlas of the zebrafish testis. Using 5 timepoints across the adult life of a zebrafish, we described cellular profiles in the testis during and after fertility. While all germ cell stages of spermatogenesis are detected in testes from fertile adult zebrafish, testes from older infertile males only contained spermatogonia and a reduced population of spermatocytes. These remaining germ cells are transcriptionally distinct from fertile spermatogonia. Immune cells including macrophages and lymphocytes drastically increase in abundance in infertile testes. Our developmental atlas reveals the cellular changes as the testis ages and defines a molecular roadmap for the regulation of male fertility.

Microbe transmission from pet shop to lab-reared zebrafish reveals a pathogenic birnavirus.

Zebrafish are popular research organisms selected for laboratory use due in part to widespread availability from the pet trade. Many contemporary colonies of laboratory zebrafish are maintained in aquaculture facilities that monitor and aim to curb infections that can negatively affect colony health and confound experiments. The impact of laboratory control on the microbial constituents associated with zebrafish in research environments compared to the pet trade are unclear. Diseases of unknown causes are common in both environments. We conducted a metagenomic survey to broadly compare the zebrafish-associated microbes in pet trade and laboratory environments. We detected many microbes in animals from the pet trade that were not found in laboratory animals. Co-housing experiments revealed several transmissible microbes including a newly described non-enveloped, double-stranded RNA virus in the Birnaviridae family we name Rocky Mountain birnavirus (RMBV). Infections were detected in asymptomatic animals from the pet trade, but when transmitted to laboratory animals RMBV was associated with pronounced antiviral responses and hemorrhagic disease. These experiments highlight the pet trade as a distinct source of diverse microbes that associate with zebrafish and establish a paradigm for the discovery of newly described pathogenic viruses and other infectious microbes that can be developed for study in the laboratory.

Distinct features of the regenerating heart uncovered through comparative single-cell profiling.

Adult humans respond to heart injury by forming a permanent scar, yet other vertebrates are capable of robust and complete cardiac regeneration. Despite progress towards characterizing the mechanisms of cardiac regeneration in fish and amphibians, the large evolutionary gulf between mammals and regenerating vertebrates complicates deciphering which cellular and molecular features truly enable regeneration. To better define these features, we compared cardiac injury responses in zebrafish and medaka, two fish species that share similar heart anatomy and common teleost ancestry but differ in regenerative capability. We used single-cell transcriptional profiling to create a time-resolved comparative cell atlas of injury responses in all major cardiac cell types across both species. With this approach, we identified several key features that distinguish cardiac injury response in the non-regenerating medaka heart. By comparing immune responses to injury, we found altered cell recruitment and a distinct pro-inflammatory gene program in medaka leukocytes, and an absence of the injury-induced interferon response seen in zebrafish. In addition, we found a lack of pro-regenerative signals, including nrg1 and retinoic acid, from medaka endothelial and epicardial cells. Finally, we identified alterations in the myocardial structure in medaka, where they lack embryonic-like primordial layer cardiomyocytes, and fail to employ a cardioprotective gene program shared by regenerating vertebrates. Our findings reveal notable variation in injury response across nearly all major cardiac cell types in zebrafish and medaka, demonstrating how evolutionary divergence influences the hidden cellular features underpinning regenerative potential in these seemingly similar vertebrates.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage.

Exposure to ultraviolet radiation (UVR) is harmful to living cells, leading organisms to evolve protective mechanisms against UVR-induced cellular damage and stress.1,2 UVR, particularly UVB (280-320 nm), can damage proteins and DNA, leading to errors during DNA repair and replication. Excessive UVR can induce cellular death. Aquatic organisms face risk of UV exposure as biologically harmful levels of UVB can penetrate >10 m in clear water.3 While melanin is the only known sunscreen in vertebrates, it often emerges late in embryonic development, rendering embryos of many species vulnerable during the earlier stages. Algae and microbes produce a class of sunscreening compounds known as mycosporine-like amino acids (MAAs).4 Fish eggs contain a similar compound called gadusol, whose role as a sunscreen has yet to be tested despite its discovery over 40 years ago.5 The recent finding that many vertebrate genomes contain a biosynthetic pathway for gadusol suggests that many fish may produce and use this molecule as a sunscreen.6 We generated a gadusol-deficient mutant zebrafish to investigate the role of gadusol in protecting fish embryos and larvae from UVR. Our results demonstrate that maternally provided gadusol is the primary sunscreen in embryonic and larval development, while melanin provides modest secondary protection. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage.

Ultraviolet radiation (UVR) and its deleterious effects on living cells selects for UVR-protective mechanisms. Organisms across the tree of life evolved a variety of natural sunscreens to prevent UVR-induced cellular damage and stress. However, in vertebrates, only melanin is known to act as a sunscreen. Here we demonstrate that gadusol, a transparent compound discovered over 40 years ago in fish eggs, is a maternally provided sunscreen required for survival of embryonic and larval zebrafish exposed to UVR. Mutating an enzyme involved in gadusol biosynthesis increases the formation of cyclobutane pyrimidine dimers, a hallmark of UVB-induced DNA damage. Compared to the contributions of melanin and the chorion, gadusol is the primary sunscreening mechanism in embryonic and larval fish. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Rational engineering of a synthetic insect-bacterial mutualism.

Many insects maintain mutualistic associations with bacterial endosymbionts, but little is known about how they originate in nature. In this study, we describe the establishment and manipulation of a synthetic insect-bacterial symbiosis in a weevil host. Following egg injection, the nascent symbiont colonized many tissues, including prototypical somatic and germinal bacteriomes, yielding maternal transmission over many generations. We then engineered the nascent symbiont to overproduce the aromatic amino acids tyrosine and phenylalanine, which facilitate weevil cuticle strengthening and accelerated larval development, replicating the function of mutualistic symbionts that are widely distributed among weevils and other beetles in nature. Our work provides empirical support for the notion that mutualistic symbioses can be initiated in insects by the acquisition of environmental bacteria. It also shows that certain bacterial genera, including the Sodalis spp. used in our study, are predisposed to develop these associations due to their ability to maintain benign infections and undergo vertical transmission in diverse insect hosts, facilitating the partner-fidelity feedback that is critical for the evolution of obligate mutualism. These experimental advances provide a new platform for laboratory studies focusing on the molecular mechanisms and evolutionary processes underlying insect-bacterial symbiosis.

Orthogonal CRISPR-Cas tools for genome editing, inhibition, and CRISPR recording in zebrafish embryos.

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas universe continues to expand. The type II CRISPR-Cas system from Streptococcus pyogenes (SpyCas9) is the most widely used for genome editing due to its high efficiency in cells and organisms. However, concentrating on a single CRISPR-Cas system imposes limits on target selection and multiplexed genome engineering. We hypothesized that CRISPR-Cas systems originating from different bacterial species could operate simultaneously and independently due to their distinct single-guide RNAs (sgRNAs) or CRISPR-RNAs (crRNAs), and protospacer adjacent motifs (PAMs). Additionally, we hypothesized that CRISPR-Cas activity in zebrafish could be regulated through the expression of inhibitory anti-CRISPR (Acr) proteins. Here, we use a simple mutagenesis approach to demonstrate that CRISPR-Cas systems from S. pyogenes (SpyCas9), Streptococcus aureus (SauCas9), Lachnospiraceae bacterium (LbaCas12a, previously known as LbCpf1) are orthogonal systems capable of operating simultaneously in zebrafish. CRISPR systems from Acidaminococcus sp. (AspCas12a, previously known as AsCpf1) and Neisseria meningitidis (Nme2Cas9) were also active in embryos. We implemented multichannel CRISPR recording using three CRISPR systems and show that LbaCas12a may provide superior information density compared with previous methods. We also demonstrate that type II Acrs (anti-CRISPRs) are effective inhibitors of SpyCas9 in zebrafish. Our results indicate that at least five CRISPR-Cas systems and two anti-CRISPR proteins are functional in zebrafish embryos. These orthogonal CRISPR-Cas systems and Acr proteins will enable combinatorial and intersectional strategies for spatiotemporal control of genome editing and genetic recording in animals.

Functional and ultrastructural analysis of reafferent mechanosensation in larval zebrafish.

All animals need to differentiate between exafferent stimuli, which are caused by the environment, and reafferent stimuli, which are caused by their own movement. In the case of mechanosensation in aquatic animals, the exafferent inputs are water vibrations in the animal's proximity, which need to be distinguishable from the reafferent inputs arising from fluid drag due to locomotion. Both of these inputs are detected by the lateral line, a collection of mechanosensory organs distributed along the surface of the body. In this study, we characterize in detail how hair cells-the receptor cells of the lateral line-in zebrafish larvae discriminate between such reafferent and exafferent signals. Using dye labeling of the lateral line nerve, we visualize two parallel descending inputs that can influence lateral line sensitivity. We combine functional imaging with ultra-structural EM circuit reconstruction to show that cholinergic signals originating from the hindbrain transmit efference copies (copies of the motor command that cancel out self-generated reafferent stimulation during locomotion) and that dopaminergic signals from the hypothalamus may have a role in threshold modulation, both in response to locomotion and salient stimuli. We further gain direct mechanistic insight into the core components of this circuit by loss-of-function perturbations using targeted ablations and gene knockouts. We propose that this simple circuit is the core implementation of mechanosensory reafferent suppression in these young animals and that it might form the first instantiation of state-dependent modulation found at later stages in development.

Linking Virus Discovery to Immune Responses Visualized during Zebrafish Infections.

The discovery of new viruses currently outpaces our capacity for experimental examination of infection biology. To better couple virus discovery with immunology, we genetically modified zebrafish to visually report on virus infections. After generating a strain that expresses green fluorescent protein (GFP) under an interferon-stimulated gene promoter, we repeatedly observed transgenic larvae spontaneously expressing GFP days after hatching. RNA sequencing comparisons of co-housed GFP-positive and GFP-negative zebrafish revealed a naturally occurring picornavirus that induced a canonical interferon-mediated response and hundreds of antiviral defense genes not observed following immunostimulatory treatments or experimental infections with other viruses. Among the many genes induced by picornavirus infection was a large set encoding guanosine triphosphatase (GTPase) of immunity-associated proteins (GIMAPs). The GIMAP gene family is massively expanded in fish genomes and may also play a crucial role in antiviral responses in mammals, including humans. We subsequently detected zebrafish picornavirus in publicly available sequencing data from seemingly asymptomatic zebrafish in many research institutes and found that it altered gene expression in a previous study of zebrafish development. Experiments revealed a horizontal mode of virus transmission, highlighting a system for studying the spread of picornavirus infections within and between individuals. Our study describes a naturally occurring picornavirus that elicits strong antiviral responses in zebrafish and provides new strategies for simultaneously discovering viruses and their impact on vertebrate hosts.

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.

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.

Recording development with single cell dynamic lineage tracing.

Every animal grows from a single fertilized egg into an intricate network of cell types and organ systems. This process is captured in a lineage tree: a diagram of every cell's ancestry back to the founding zygote. Biologists have long sought to trace this cell lineage tree in individual organisms and have developed a variety of technologies to map the progeny of specific cells. However, there are billions to trillions of cells in complex organisms, and conventional approaches can only map a limited number of clonal populations per experiment. A new generation of tools that use molecular recording methods integrated with single cell profiling technologies may provide a solution. Here, we summarize recent breakthroughs in these technologies, outline experimental and computational challenges, and discuss biological questions that can be addressed using single cell dynamic lineage tracing.

Accurate analysis of genuine CRISPR editing events with ampliCan.

We present ampliCan, an analysis tool for genome editing that unites highly precise quantification and visualization of genuine genome editing events. ampliCan features nuclease-optimized alignments, filtering of experimental artifacts, event-specific normalization, and off-target read detection and quantifies insertions, deletions, HDR repair, as well as targeted base editing. It is scalable to thousands of amplicon sequencing-based experiments from any genome editing experiment, including CRISPR. It enables automated integration of controls and accounts for biases at every step of the analysis. We benchmarked ampliCan on both real and simulated data sets against other leading tools, demonstrating that it outperformed all in the face of common confounding factors.

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.

Simultaneous single-cell profiling of lineages and cell types in the vertebrate brain.

The lineage relationships among the hundreds of cell types generated during development are difficult to reconstruct. A recent method, GESTALT, used CRISPR-Cas9 barcode editing for large-scale lineage tracing, but was restricted to early development and did not identify cell types. Here we present scGESTALT, which combines the lineage recording capabilities of GESTALT with cell-type identification by single-cell RNA sequencing. The method relies on an inducible system that enables barcodes to be edited at multiple time points, capturing lineage information from later stages of development. Sequencing of ∼60,000 transcriptomes from the juvenile zebrafish brain identified >100 cell types and marker genes. Using these data, we generate lineage trees with hundreds of branches that help uncover restrictions at the level of cell types, brain regions, and gene expression cascades during differentiation. scGESTALT can be applied to other multicellular organisms to simultaneously characterize molecular identities and lineage histories of thousands of cells during development and disease.

The primary role of zebrafish nanog is in extra-embryonic tissue.

The role of the zebrafish transcription factor Nanog has been controversial. It has been suggested that Nanog is primarily required for the proper formation of the extra-embryonic yolk syncytial layer (YSL) and only indirectly regulates gene expression in embryonic cells. In an alternative scenario, Nanog has been proposed to directly regulate transcription in embryonic cells during zygotic genome activation. To clarify the roles of Nanog, we performed a detailed analysis of zebrafish nanog mutants. Whereas zygotic nanog mutants survive to adulthood, maternal-zygotic (MZnanog) and maternal mutants exhibit developmental arrest at the blastula stage. In the absence of Nanog, YSL formation and epiboly are abnormal, embryonic tissue detaches from the yolk, and the expression of dozens of YSL and embryonic genes is reduced. Epiboly defects can be rescued by generating chimeric embryos of MZnanog embryonic tissue with wild-type vegetal tissue that includes the YSL and yolk cell. Notably, cells lacking Nanog readily respond to Nodal signals and when transplanted into wild-type hosts proliferate and contribute to embryonic tissues and adult organs from all germ layers. These results indicate that zebrafish Nanog is necessary for proper YSL development but is not directly required for embryonic cell differentiation.