The reference genome of Macropodus opercularis (the paradise fish).

Amongst fishes, zebrafish (Danio rerio) has gained popularity as a model system over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species has a highly complex behavioral repertoire and has been the subject of many ethological investigations but lacks genomic resources. Here we report the reference genome assembly of M. opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of 483,077,705 base pairs (~483 Mb) on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ~90% of them to orthogroups.

Stat3 Has a Different Role in Axon Growth During Development Than It Does in Axon Regeneration After Injury.

Signal transducer and activator of transcription 3 (STAT3) is essential for neural development and regeneration as a key transcription factor and mitochondrial activator. However, the mechanism of Stat3 in axon development and regeneration has not been fully understood. In this study, using zebrafish posterior lateral line (PLL) axons, we demonstrate that Stat3 plays distinct roles in PLL axon embryonic growth and regeneration. Our experiments indicate that stat3 is required for PLL axon extension. In stat3 mutant zebrafish, the PLL axon ends were stalled at the level of the cloaca, and expression of stat3 rescues the PLL axon growth in a cell-autonomous manner. Jak/Stat signaling inhibition did not affect PLL axon growth indicating Jak/Stat was dispensable for PLL axon growth. In addition, we found that Stat3 was co-localized with mitochondria in PLL axons and important for the mitochondrial membrane potential and ATPase activity. The PLL axon growth defect of stat3 mutants was mimicked and rescued by rotenone and DCHC treatment, respectively, which suggests that Stat3 regulates PLL axon growth through mitochondrial Stat3. By contrast, mutation of stat3 or Jak/Stat signaling inhibition retarded PLL axon regeneration. Meanwhile, we also found Schwann cell migration was also inhibited in stat3 mutants. Taken together, Stat3 is required for embryonic PLL axon growth by regulating the ATP synthesis efficiency of mitochondria, whereas Stat3 stimulates PLL axon regeneration by regulating Schwann cell migration via Jak/Stat signaling. Our findings show a new mechanism of Stat3 in axon growth and regeneration.

Low mutation rate in epaulette sharks is consistent with a slow rate of evolution in sharks.

Sharks occupy diverse ecological niches and play critical roles in marine ecosystems, often acting as apex predators. They are considered a slow-evolving lineage and have been suggested to exhibit exceptionally low cancer rates. These two features could be explained by a low nuclear mutation rate. Here, we provide a direct estimate of the nuclear mutation rate in the epaulette shark (Hemiscyllium ocellatum). We generate a high-quality reference genome, and resequence the whole genomes of parents and nine offspring to detect de novo mutations. Using stringent criteria, we estimate a mutation rate of 7×10-10 per base pair, per generation. This represents one of the lowest directly estimated mutation rates for any vertebrate clade, indicating that this basal vertebrate group is indeed a slowly evolving lineage whose ability to restore genetic diversity following a sustained population bottleneck may be hampered by a low mutation rate.

The reference genome of the paradise fish (Macropodus opercularis).

Over the decades, a small number of model species, each representative of a larger taxa, have dominated the field of biological research. Amongst fishes, zebrafish (Danio rerio) has gained popularity over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species from Southeast Asia, has a highly complex behavioral repertoire and has been the subject of many ethological investigations, but lacks genomic resources. Here we report the reference genome assembly of Macropodus opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of ≈483 Mb on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ≈90% of them to orthogroups. Completeness analysis showed that 98.5% of the Actinopterygii core gene set (ODB10) was present as a complete ortholog in our reference genome with a further 1.2 % being present in a fragmented form. Additionally, we cloned multiple genes important during early development and using newly developed in situ hybridization protocols, we showed that they have conserved expression patterns.

Ontogenetically distinct neutrophils differ in function and transcriptional profile in zebrafish.

The current view of hematopoiesis considers leukocytes on a continuum with distinct developmental origins, and which exert non-overlapping functions. However, there is less known about the function and phenotype of ontogenetically distinct neutrophil populations. In this work, using a photoconvertible transgenic zebrafish line; Tg(mpx:Dendra2), we selectively label rostral blood island-derived and caudal hematopoietic tissue-derived neutrophils in vivo during steady state or upon injury. By comparing the migratory properties and single-cell expression profiles of both neutrophil populations at steady state we show that rostral neutrophils show higher csf3b expression and migration capacity than caudal neutrophils. Upon injury, both populations share a core transcriptional profile as well as subset-specific transcriptional signatures. Accordingly, both rostral and caudal neutrophils are recruited to the wound independently of their distance to the injury. While rostral neutrophils respond uniformly, caudal neutrophils respond heterogeneously. Collectively, our results reveal that co-existing neutrophils populations with ontogenically distinct origin display functional differences.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization.

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

A regulatory network of Sox and Six transcription factors initiate a cell fate transformation during hearing regeneration in adult zebrafish.

Using adult zebrafish inner ears as a model for sensorineural regeneration, we ablated the mechanosensory receptors and characterized the single-cell epigenome and transcriptome at consecutive time points during hair cell regeneration. We utilized deep learning on the regeneration-induced open chromatin sequences and identified cell-specific transcription factor (TF) motif patterns. Enhancer activity correlated with gene expression and identified potential gene regulatory networks. A pattern of overlapping Sox- and Six-family TF gene expression and binding motifs was detected, suggesting a combinatorial program of TFs driving regeneration and cell identity. Pseudotime analysis of single-cell transcriptomic data suggested that support cells within the sensory epithelium changed cell identity to a "progenitor" cell population that could differentiate into hair cells. We identified a 2.6 kb DNA enhancer upstream of the sox2 promoter that, when deleted, showed a dominant phenotype that resulted in a hair-cell-regeneration-specific deficit in both the lateral line and adult inner ear.

Zrsr2 Is Essential for the Embryonic Development and Splicing of Minor Introns in RNA and Protein Processing Genes in Zebrafish.

ZRSR2 (zinc finger CCCH-type, RNA binding motif and serine/arginine rich 2) is an essential splicing factor involved in 3' splice-site recognition as a component of both the major and minor spliceosomes that mediate the splicing of U2-type (major) and U12-type (minor) introns, respectively. Studies of ZRSR2-depleted cell lines and ZRSR2-mutated patient samples revealed its essential role in the U12-dependent minor spliceosome. However, the role of ZRSR2 during embryonic development is not clear, as its function is compensated for by Zrsr1 in mice. Here, we utilized the zebrafish model to investigate the role of zrsr2 during embryonic development. Using CRISPR/Cas9 technology, we generated a zrsr2-knockout zebrafish line, termed zrsr2hg129/hg129 (p.Trp167Argfs*9) and examined embryo development in the homozygous mutant embryos. zrsr2hg129/hg129 embryos displayed multiple developmental defects starting at 4 days post fertilization (dpf) and died after 8 dpf, suggesting that proper Zrsr2 function is required during embryonic development. The global transcriptome analysis of 3 dpf zrsr2hg129/hg129 embryos revealed that the loss of Zrsr2 results in the downregulation of essential metabolic pathways and the aberrant retention of minor introns in about one-third of all minor intron-containing genes in zebrafish. Overall, our study has demonstrated that the role of Zrsr2 as a component of the minor spliceosome is conserved and critical for proper embryonic development in zebrafish.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements.

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Chondroitin/dermatan sulfate glycosyltransferase genes are essential for craniofacial development.

Chondroitin/dermatan sulfate (CS/DS) proteoglycans are indispensable for animal development and homeostasis but the large number of enzymes involved in their biosynthesis have made CS/DS function a challenging problem to study genetically. In our study, we generated loss-of-function alleles in zebrafish genes encoding CS/DS biosynthetic enzymes and characterized the effect on development in single and double mutants. Homozygous mutants in chsy1, csgalnact1a, csgalnat2, chpfa, ust and chst7, respectively, develop to adults. However, csgalnact1a-/- fish develop distinct craniofacial defects while the chsy1-/- skeletal phenotype is milder and the remaining mutants display no gross morphological abnormalities. These results suggest a high redundancy for the CS/DS biosynthetic enzymes and to further reduce CS/DS biosynthesis we combined mutant alleles. The craniofacial phenotype is further enhanced in csgalnact1a-/-;chsy1-/- adults and csgalnact1a-/-;csgalnact2-/- larvae. While csgalnact1a-/-;csgalnact2-/- was the most affected allele combination in our study, CS/DS is still not completely abolished. Transcriptome analysis of chsy1-/-, csgalnact1a-/- and csgalnact1a-/-;csgalnact2-/- larvae revealed that the expression had changed in a similar way in the three mutant lines but no differential expression was found in any of fifty GAG biosynthesis enzymes identified. Thus, zebrafish larvae do not increase transcription of GAG biosynthesis genes as a consequence of decreased CS/DS biosynthesis. The new zebrafish lines develop phenotypes similar to clinical characteristics of several human congenital disorders making the mutants potentially useful to study disease mechanisms and treatment.

Darwinian genomics and diversity in the tree of life.

Genomics encompasses the entire tree of life, both extinct and extant, and the evolutionary processes that shape this diversity. To date, genomic research has focused on humans, a small number of agricultural species, and established laboratory models. Fewer than 18,000 of ∼2,000,000 eukaryotic species (<1%) have a representative genome sequence in GenBank, and only a fraction of these have ancillary information on genome structure, genetic variation, gene expression, epigenetic modifications, and population diversity. This imbalance reflects a perception that human studies are paramount in disease research. Yet understanding how genomes work, and how genetic variation shapes phenotypes, requires a broad view that embraces the vast diversity of life. We have the technology to collect massive and exquisitely detailed datasets about the world, but expertise is siloed into distinct fields. A new approach, integrating comparative genomics with cell and evolutionary biology, ecology, archaeology, anthropology, and conservation biology, is essential for understanding and protecting ourselves and our world. Here, we describe potential for scientific discovery when comparative genomics works in close collaboration with a broad range of fields as well as the technical, scientific, and social constraints that must be addressed.

Singling out how genes are regulated during development.

In this issue of Cell Genomics, McGarvey et al. characterize chromatin accessibility and gene regulation at single-cell resolution in the zebrafish embryo 24 h after fertilization, providing a valuable resource. Their findings add another important piece to understanding the dynamic landscape of gene expression and regulation during early vertebrate development and hint at the dramatic changes single-cell genomics is bringing to the field of developmental biology.

Vestibular and Auditory Hair Cell Regeneration Following Targeted Ablation of Hair Cells With Diphtheria Toxin in Zebrafish.

Millions of Americans experience hearing or balance disorders due to loss of hair cells in the inner ear. The hair cells are mechanosensory receptors used in the auditory and vestibular organs of all vertebrates as well as the lateral line systems of aquatic vertebrates. In zebrafish and other non-mammalian vertebrates, hair cells turnover during homeostasis and regenerate completely after being destroyed or damaged by acoustic or chemical exposure. However, in mammals, destroying or damaging hair cells results in permanent impairments to hearing or balance. We sought an improved method for studying hair cell damage and regeneration in adult aquatic vertebrates by generating a transgenic zebrafish with the capacity for targeted and inducible hair cell ablation in vivo. This model expresses the human diphtheria toxin receptor (hDTR) gene under the control of the myo6b promoter, resulting in hDTR expressed only in hair cells. Cell ablation is achieved by an intraperitoneal injection of diphtheria toxin (DT) in adult zebrafish or DT dissolved in the water for larvae. In the lateral line of 5 days post fertilization (dpf) zebrafish, ablation of hair cells by DT treatment occurred within 2 days in a dose-dependent manner. Similarly, in adult utricles and saccules, a single intraperitoneal injection of 0.05 ng DT caused complete loss of hair cells in the utricle and saccule by 5 days post-injection. Full hair cell regeneration was observed for the lateral line and the inner ear tissues. This study introduces a new method for efficient conditional hair cell ablation in adult zebrafish inner ear sensory epithelia (utricles and saccules) and demonstrates that zebrafish hair cells will regenerate in vivo after this treatment.

Suppressing STAT3 activity protects the endothelial barrier from VEGF-mediated vascular permeability.

Vascular permeability triggered by inflammation or ischemia promotes edema, exacerbates disease progression and impairs tissue recovery. Vascular endothelial growth factor (VEGF) is a potent inducer of vascular permeability. VEGF plays an integral role in regulating vascular barrier function physiologically and in pathologies, including cancer, stroke, cardiovascular disease, retinal conditions and COVID-19-associated pulmonary edema, sepsis and acute lung injury. Understanding temporal molecular regulation of VEGF-induced vascular permeability will facilitate developing therapeutics to inhibit vascular permeability, while preserving tissue-restorative angiogenesis. Here, we demonstrate that VEGF signals through signal transducer and activator of transcription 3 (STAT3) to promote vascular permeability. We show that genetic STAT3 ablation reduces vascular permeability in STAT3-deficient endothelium of mice and VEGF-inducible zebrafish crossed with CRISPR/Cas9-generated Stat3 knockout zebrafish. Intercellular adhesion molecule 1 (ICAM-1) expression is transcriptionally regulated by STAT3, and VEGF-dependent STAT3 activation is regulated by JAK2. Pyrimethamine, an FDA-approved antimicrobial agent that inhibits STAT3-dependent transcription, substantially reduces VEGF-induced vascular permeability in zebrafish, mouse and human endothelium. Collectively, our findings suggest that VEGF/VEGFR-2/JAK2/STAT3 signaling regulates vascular barrier integrity, and inhibition of STAT3-dependent activity reduces VEGF-induced vascular permeability. This article has an associated First Person interview with the first author of the paper.

The Warburg effect is necessary to promote glycosylation in the blastema during zebrafish tail regeneration.

Throughout their lifetime, fish maintain a high capacity for regenerating complex tissues after injury. We utilized a larval tail regeneration assay in the zebrafish Danio rerio, which serves as an ideal model of appendage regeneration due to its easy manipulation, relatively simple mixture of cell types, and superior imaging properties. Regeneration of the embryonic zebrafish tail requires development of a blastema, a mass of dedifferentiated cells capable of replacing lost tissue, a crucial step in all known examples of appendage regeneration. Using this model, we show that tail amputation triggers an obligate metabolic shift to promote glucose metabolism during early regeneration similar to the Warburg effect observed in tumor forming cells. Inhibition of glucose metabolism did not affect the overall health of the embryo but completely blocked the tail from regenerating after amputation due to the failure to form a functional blastema. We performed a time series of single-cell RNA sequencing on regenerating tails with and without inhibition of glucose metabolism. We demonstrated that metabolic reprogramming is required for sustained TGF-ß signaling and blocking glucose metabolism largely mimicked inhibition of TGF-ß receptors, both resulting in an aberrant blastema. Finally, we showed using genetic ablation of three possible metabolic pathways for glucose, that metabolic reprogramming is required to provide glucose specifically to the hexosamine biosynthetic pathway while neither glycolysis nor the pentose phosphate pathway were necessary for regeneration.

The Formin Fmn2b Is Required for the Development of an Excitatory Interneuron Module in the Zebrafish Acoustic Startle Circuit.

The formin family member Fmn2 is a neuronally enriched cytoskeletal remodeling protein conserved across vertebrates. Recent studies have implicated Fmn2 in neurodevelopmental disorders, including sensory processing dysfunction and intellectual disability in humans. Cellular characterization of Fmn2 in primary neuronal cultures has identified its function in the regulation of cell-substrate adhesion and consequently growth cone translocation. However, the role of Fmn2 in the development of neural circuits in vivo, and its impact on associated behaviors have not been tested. Using automated analysis of behavior and systematic investigation of the associated circuitry, we uncover the role of Fmn2b in zebrafish neural circuit development. As reported in other vertebrates, the zebrafish ortholog of Fmn2 is also enriched in the developing zebrafish nervous system. We find that Fmn2b is required for the development of an excitatory interneuron pathway, the spiral fiber neuron, which is an essential circuit component in the regulation of the Mauthner cell (M-cell)-mediated acoustic startle response. Consistent with the loss of the spiral fiber neurons tracts, high-speed video recording revealed a reduction in the short latency escape events while responsiveness to the stimuli was unaffected. Taken together, this study provides evidence for a circuit-specific requirement of Fmn2b in eliciting an essential behavior in zebrafish. Our findings underscore the importance of Fmn2 in neural development across vertebrate lineages and highlight zebrafish models in understanding neurodevelopmental disorders.

Biallelic variants in KARS1 are associated with neurodevelopmental disorders and hearing loss recapitulated by the knockout zebrafish.

PURPOSE: Pathogenic variants in Lysyl-tRNA synthetase 1 (KARS1) have increasingly been recognized as a cause of early-onset complex neurological phenotypes. To advance the timely diagnosis of KARS1-related disorders, we sought to delineate its phenotype and generate a disease model to understand its function in vivo. METHODS: Through international collaboration, we identified 22 affected individuals from 16 unrelated families harboring biallelic likely pathogenic or pathogenic in KARS1 variants. Sequencing approaches ranged from disease-specific panels to genome sequencing. We generated loss-of-function alleles in zebrafish. RESULTS: We identify ten new and four known biallelic missense variants in KARS1 presenting with a moderate-to-severe developmental delay, progressive neurological and neurosensory abnormalities, and variable white matter involvement. We describe novel KARS1-associated signs such as autism, hyperactive behavior, pontine hypoplasia, and cerebellar atrophy with prevalent vermian involvement. Loss of kars1 leads to upregulation of p53, tissue-specific apoptosis, and downregulation of neurodevelopmental related genes, recapitulating key tissue-specific disease phenotypes of patients. Inhibition of p53 rescued several defects of kars1-/- knockouts. CONCLUSION: Our work delineates the clinical spectrum associated with KARS1 defects and provides a novel animal model for KARS1-related human diseases revealing p53 signaling components as potential therapeutic targets.

Cnr2 Is Important for Ribbon Synapse Maturation and Function in Hair Cells and Photoreceptors.

The role of the cannabinoid receptor 2 (CNR2) is still poorly described in sensory epithelia. We found strong cnr2 expression in hair cells (HCs) of the inner ear and the lateral line (LL), a superficial sensory structure in fish. Next, we demonstrated that sensory synapses in HCs were severely perturbed in larvae lacking cnr2. Appearance and distribution of presynaptic ribbons and calcium channels (Cav1.3) were profoundly altered in mutant animals. Clustering of membrane-associated guanylate kinase (MAGUK) in post-synaptic densities (PSDs) was also heavily affected, suggesting a role for cnr2 for maintaining the sensory synapse. Furthermore, vesicular trafficking in HCs was strongly perturbed suggesting a retrograde action of the endocannabinoid system (ECs) via cnr2 that was modulating HC mechanotransduction. We found similar perturbations in retinal ribbon synapses. Finally, we showed that larval swimming behaviors after sound and light stimulations were significantly different in mutant animals. Thus, we propose that cnr2 is critical for the processing of sensory information in the developing larva.

A New Zebrafish Model for Pseudoxanthoma Elasticum.

Calcification of various tissues is a significant health issue associated with aging, cancer and autoimmune diseases. There are both environmental and genetic factors behind this phenomenon and understanding them is essential for the development of efficient therapeutic approaches. Pseudoxanthoma elasticum (PXE) is a rare genetic disease, a prototype for calcification disorders, resulting from the dysfunction of ABCC6, a transport protein found in the membranes of cells. It is identified by excess calcification in a variety of tissues (e.g., eyes, skin, arteries) and currently it has no cure, known treatments target the symptoms only. Preclinical studies of PXE have been successful in mice, proving the usefulness of animal models for the study of the disease. Here, we present a new zebrafish (Danio rerio) model for PXE. By resolving some ambiguous assemblies in the zebrafish genome, we show that there are two functional and one non-functional paralogs for ABCC6 in zebrafish (abcc6a, abcc6b.1, and abcc6b.2, respectively). We created single and double mutants for the functional paralogs and characterized their calcification defects with a combination of techniques. Zebrafish deficient in abcc6a show defects in their vertebral calcification and also display ectopic calcification foci in their soft tissues. Our results also suggest that the impairment of abcc6b.1 does not affect this biological process.