Zebrafish Information Network, the knowledgebase for Danio rerio research

Complete persistence of the primary somatosensory system in zebrafish

The somatosensory system detects peripheral stimuli that are translated into behaviors necessary for survival. Fishes and amphibians possess two somatosensory systems in the trunk: the primary somatosensory system, formed by the Rohon-Beard neurons, and the secondary somatosensory system, formed by the neural crest cell-derived neurons of the Dorsal Root Ganglia. Rohon-Beard neurons have been characterized as a transient population that mostly disappears during the first days of life and is functionally replaced by the Dorsal Root Ganglia. Here, I follow Rohon-Beard neurons in vivo and show that the entire repertoire remains present in zebrafish from 1-day post-fertilization until the juvenile stage, 15-days post-fertilization. These data indicate that zebrafish retain two complete somatosensory systems until at least a developmental stage when the animals display complex behavioral repertoires.

Examining the liver-pancreas crosstalk reveals a role for the molybdenum cofactor in β-cell regeneration

Regeneration of insulin-producing β-cells is an alternative avenue to manage diabetes, and it is crucial to unravel this process in vivo during physiological responses to the lack of β-cells. Here, we aimed to characterize how hepatocytes can contribute to β-cell regeneration, either directly or indirectly via secreted proteins or metabolites, in a zebrafish model of β-cell loss. Using lineage tracing, we show that hepatocytes do not directly convert into β-cells even under extreme β-cell ablation conditions. A transcriptomic analysis of isolated hepatocytes after β-cell ablation displayed altered lipid- and glucose-related processes. Based on the transcriptomics, we performed a genetic screen that uncovers a potential role of the molybdenum cofactor (Moco) biosynthetic pathway in β-cell regeneration and glucose metabolism in zebrafish. Consistently, molybdenum cofactor synthesis 2 (Mocs2) haploinsufficiency in mice indicated dysregulated glucose metabolism and liver function. Together, our study sheds light on the liver-pancreas crosstalk and suggests that the molybdenum cofactor biosynthesis pathway should be further studied in relation to glucose metabolism and diabetes.

Dexamethasone reduces cisplatin-induced hair cell damage by inducing cisplatin resistance through metallothionein-2

PURPOSE

Hair cell damage is a common side effect caused by the anticancer drug cisplatin (CDDP), which reduces patient quality of life. One CDDP resistance mechanism that occurs in recurrent cancers is heavy metal detoxification by metallothionein-2 (mt2). Here, we show that in zebrafish larvae, dexamethasone (DEX) reduces CDDP-induced hair cell damage by enhancing mt2 expression.

METHODS

Transgenic zebrafish (cldn: gfp; atoh1: rfp) that express green and red fluorescent proteins in neuromasts and hair cells, respectively, were used. The zebrafish were pretreated with DEX at 52 h post-fertilization (hpf) for 8 h, followed by CDDP treatment for 12 h. The lateral line hair cells of CDDP-treated zebrafish at 72 hpf were observed by fluorescence microscopy.

RESULTS

Reporting odds ratio (ROR) analysis using an adverse event database indicated an association between a decrease in CDDP-induced ototoxicity and DEX as an antiemetic treatment for cancer chemotherapy. Pretreatment with DEX protected 72 hpf zebrafish hair cells from CDDP-induced damage. The expression of mt2 mRNA was significantly increased by the combination of 10 µM DEX with CDDP. Gene editing of mt2 reversed the protective effect of DEX against CDDP-induced damage in hair cells.

CONCLUSION

DEX protects hair cells from CDDP-induced damage through increased mt2 expression, which is a resistance mechanism for platinum-based anticancer drugs.

Fibroblast Growth Factor (FGF) 13

Fibroblast Growth Factor (FGF) 13, also referred to as FGF homologous factor (FHF) 2, is a member of the FGF11 subfamily that is characterized as having sequence similarities to classical FGF receptor (FGFR)-binding FGFs, but functionally do not bind FGFRs. In this primer mini-review, we summarize current knowledge regarding FGF13 expression, mutant analyses, and gene and protein structure. Similar to other FHFs, FGF13 has been considered a non-secreted protein that lacks an amino signal and is prominently expressed in developing and mature neurons of the central and peripheral nervous systems, as well as the heart. The expression of FGF13 is not limited to early embryonic stages and has been shown to persist in adult tissues. As well, FGF13 is known to localize subcellularly, both within the cytoplasm and the nucleus. FGF13 is extremely adaptable, as it interacts with MAPK scaffolding protein islet brain 2 (IB2), stabilizes microtubules, or binds to voltage-gated sodium channels. Fgf13 mutant mouse lines display various neurological pathologies. Through sequence mapping, FGF13 is considered a candidate causative gene that is mutated in multiple human X-linked neurological diseases.

Zebrafish Models for Skeletal and Extraskeletal Osteogenesis Imperfecta Features: Unveiling Pathophysiology and Paving the Way for Drug Discovery

In the last decades, the easy genetic manipulation, the external fertilization, the high percentage of homology with human genes and the reduced husbandry costs compared to rodents, made zebrafish a valid model for studying human diseases and for developing new therapeutical strategies. Since zebrafish shares with mammals the same bone cells and ossification types, it became widely used to dissect mechanisms and possible new therapeutic approaches in the field of common and rare bone diseases, such as osteoporosis and osteogenesis imperfecta (OI), respectively. OI is a heritable skeletal disorder caused by defects in gene encoding collagen I or proteins/enzymes necessary for collagen I synthesis and secretion. Nevertheless, OI patients can be also characterized by extraskeletal manifestations such as dentinogenesis imperfecta, muscle weakness, cardiac valve and pulmonary abnormalities and skin laxity. In this review, we provide an overview of the available zebrafish models for both dominant and recessive forms of OI. An updated description of all the main similarities and differences between zebrafish and mammal skeleton, muscle, heart and skin, will be also discussed. Finally, a list of high- and low-throughput techniques available to exploit both larvae and adult OI zebrafish models as unique tools for the discovery of new therapeutic approaches will be presented.

Zebrafish models of human-duplicated gene SRGAP2 reveal novel functions in microglia and visual system development

Recent expansion of duplicated genes unique in the Homo lineage likely contributed to brain evolution and other human-specific traits. One hallmark example is the expansion of the human SRGAP2 family, resulting in a human-specific paralog SRGAP2C . Introduction of SRGAP2C in mouse models is associated with altering cortical neuronal migration, axon guidance, synaptogenesis, and sensory-task performance. Truncated, human-specific SRGAP2C heterodimerizes with the full-length ancestral gene product SRGAP2A and antagonizes its functions. However, the significance of SRGAP2 duplication beyond neocortex development has not been elucidated due to the embryonic lethality of complete Srgap2 knockout in mice. Using zebrafish, we showed that srgap2 knockout results in viable offspring that phenocopy "humanized" SRGAP2C larvae. Specifically, human SRGAP2C protein interacts with zebrafish Srgap2, demonstrating similar Srgap2 functional antagonism observed in mice. Shared traits between knockout and humanized zebrafish larvae include altered morphometric features (i.e., reduced body length and inter-eye distance) and differential expression of synapse-, axogenesis-, vision-related genes. Through single-cell transcriptome analysis, we further observed a skewed balance of excitatory and inhibitory neurons that likely contributes to increased susceptibility to seizures displayed by Srgap2 mutant larvae, a phenotype resembling SRGAP2 loss-of-function in a child with early infantile epileptic encephalopathy. Single-cell data also pointed to strong microglia expression of srgap2 with mutants exhibiting altered membrane dynamics and likely delayed maturation of microglial cells. srgap2 -expressing microglia cells were also detected in the developing eye together with altered expression of genes related to axogenesis and synaptogenesis in mutant retinal cells. Consistent with the perturbed gene expression in the retina, we found that SRGAP2 mutant larvae exhibited increased sensitivity to broad and fine visual cues. Finally, comparing the transcriptomes of relevant cell types between human (+ SRGAP2C ) and non-human primates (- SRGAP2C ) revealed significant overlaps of gene alterations with mutant cells in our zebrafish models; this suggests that SRGAP2C plays similar roles altering microglia and the visual system in modern humans. Together, our functional characterization of zebrafish Srgap2 and human SRGAP2C in zebrafish uncovered novel gene functions and highlights the strength of cross-species analysis in understanding the development of human-specific features.

Toward universal cell embeddings: integrating single-cell RNA-seq datasets across species with SATURN

Analysis of single-cell datasets generated from diverse organisms offers unprecedented opportunities to unravel fundamental evolutionary processes of conservation and diversification of cell types. However, interspecies genomic differences limit the joint analysis of cross-species datasets to homologous genes. Here we present SATURN, a deep learning method for learning universal cell embeddings that encodes genes' biological properties using protein language models. By coupling protein embeddings from language models with RNA expression, SATURN integrates datasets profiled from different species regardless of their genomic similarity. SATURN can detect functionally related genes coexpressed across species, redefining differential expression for cross-species analysis. Applying SATURN to three species whole-organism atlases and frog and zebrafish embryogenesis datasets, we show that SATURN can effectively transfer annotations across species, even when they are evolutionarily remote. We also demonstrate that SATURN can be used to find potentially divergent gene functions between glaucoma-associated genes in humans and four other species.

Optimization of Zebrafish Larvae 6-OHDA Exposure for Neurotoxin Induced Dopaminergic Marker Reduction

Parkinson's disease (PD) is a neurodegenerative disorder that is clinically assessed by motor symptoms associated with the loss of midbrain dopaminergic neurons affecting the quality of life for over 8.5 million people worldwide. The neurotoxin 6-hydroxydopamine (6-OHDA) has been used to chemically induce a PD-like state in zebrafish larvae by several laboratories; however, highly variable concentration, methodology, and reagents have resulted in conflicting results suggesting a need to investigate these issues of reproducibility. We propose a protocol that addresses the differences in methodology and induces changes in 6 days postfertilization (dpf) larvae utilizing a 24-h exposure at 3 dpf with 30 μM 6-OHDA. Despite ∼50% lethality, no morphological or development differences in surviving fish are observed. Definition of our model is defined by downregulation of the expression of th1 by reverse transcriptase-quantitative polymerase chain reaction, a marker for dopaminergic neurons and a reduction in movement. Additionally, we observed a downregulation of pink1 and an upregulation of sod1 and sod2, indicators of mitochondrial dysfunction and response to reactive oxygen species, respectively.

Mapping the chromatin accessibility landscape of zebrafish embryogenesis at single-cell resolution by SPATAC-seq

Currently, the dynamic accessible elements that determine regulatory programs responsible for the unique identity and function of each cell type during zebrafish embryogenesis lack detailed study. Here we present SPATAC-seq: a split-pool ligation-based assay for transposase-accessible chromatin using sequencing. Using SPATAC-seq, we profiled chromatin accessibility in more than 800,000 individual nuclei across 20 developmental stages spanning the sphere stage to the early larval protruding mouth stage. Using this chromatin accessibility map, we identified 604 cell states and inferred their developmental relationships. We also identified 959,040 candidate cis-regulatory elements (cCREs) and delineated development-specific cCREs, as well as transcription factors defining diverse cell identities. Importantly, enhancer reporter assays confirmed that the majority of tested cCREs exhibited robust enhanced green fluorescent protein expression in restricted cell types or tissues. Finally, we explored gene regulatory programs that drive pigment and notochord cell differentiation. Our work provides a valuable open resource for exploring driver regulators of cell fate decisions in zebrafish embryogenesis.

Deficiency in Prader-Willi syndrome gene necdin leads to attenuated cardiac contractility

Prader-Willi syndrome (PWS) is a genetic disorder characterized by behavioral disturbances, hyperphagia, and intellectual disability. Several surveys indicate that PWS is also associated with cardiac abnormalities, possibly contributing to a high incidence of sudden death. However, the pathological mechanisms underlying cardiac dysfunction in PWS remain unclear. In this study, we found that deficiency in necdin, an intronless gene within PWS region, led to heart systolic and diastolic dysfunction in mice. Through yeast two-hybrid screening, we identified an interaction between necdin and non-muscle myosin regulatory light chain 12a/b (MYL12 A/B). We further showed that necdin stabilized MYL12 A/B via SGT1-heat shock protein 90 (HSP90) chaperone machinery. The zebrafish lacking the MYL12 A/B analog, MYL12.1, exhibited impaired heart function, while cardiac-specific overexpression of MYL12A normalized the heart dysfunction in necdin-deficient mice. Our findings revealed necdin dysfunction as a contributing factor to cardiomyopathy in PWS patients and emphasized the importance of HSP90 chaperone machinery and non-muscle myosin in heart fitness.

Rbpms2 promotes female fate upstream of the nutrient sensing Gator2 complex component Mios

Reproductive success relies on proper establishment and maintenance of biological sex. In many animals, including mammals, the primary gonad is initially ovary biased. We previously showed the RNA binding protein (RNAbp), Rbpms2, is required for ovary fate in zebrafish. Here, we identified Rbpms2 targets in oocytes (Rbpms2-bound oocyte RNAs; rboRNAs). We identify Rbpms2 as a translational regulator of rboRNAs, which include testis factors and ribosome biogenesis factors. Further, genetic analyses indicate that Rbpms2 promotes nucleolar amplification via the mTorc1 signaling pathway, specifically through the mTorc1-activating Gap activity towards Rags 2 (Gator2) component, Missing oocyte (Mios). Cumulatively, our findings indicate that early gonocytes are in a dual poised, bipotential state in which Rbpms2 acts as a binary fate-switch. Specifically, Rbpms2 represses testis factors and promotes oocyte factors to promote oocyte progression through an essential Gator2-mediated checkpoint, thereby integrating regulation of sexual differentiation factors and nutritional availability pathways in zebrafish oogenesis.

Tissue Accumulation and Biotransformation of 6PPD-Quinone in Adult Zebrafish and Its Effects on the Intestinal Microbial Community

The pronounced lethality of N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine quinone (6PPD-quinone or 6PPDQ) toward specific salmonids, while sparing other fish species, has received considerable attention. However, the underlying cause of this species-specific toxicity remains unresolved. This study explored 6PPDQ toxicokinetics and intestinal microbiota composition in adult zebrafish during a 14-day exposure to environmentally realistic concentrations, followed by a 7-day recovery phase. Predominant accumulation occurred in the brain, intestine, and eyes, with the lowest levels in the liver. Six metabolites were found to undergo hydroxylation, with two additionally undergoing O-sulfonation. Semiquantitative analyses revealed that the predominant metabolite featured a hydroxy group situated on the phenyl ring adjacent to the quinone. This was further validated by assessing enzyme activity and determining in silico binding interactions. Notably, the binding affinity between 6PPDQ and zebrafish phase I and II enzymes exceeded that with the corresponding coho salmon enzymes by 1.04-1.53 times, suggesting a higher potential for 6PPDQ detoxification in tolerant species. Whole-genome sequencing revealed significant increases in the genera Nocardioides and Rhodococcus after exposure to 6PPDQ. Functional annotation and pathway enrichment analyses predicted that these two genera would be responsible for the biodegradation and metabolism of xenobiotics. These findings offer crucial data for comprehending 6PPDQ-induced species-specific toxicity.

pIGLET: Safe harbor landing sites for reproducible and efficient transgenesis in zebrafish

Standard zebrafish transgenesis involves random transgene integration with resource-intensive screening. While phiC31 integrase-based attP/attB recombination has streamlined transgenesis in mice and Drosophila, validated attP-based landing sites for universal applications are lacking in zebrafish. Here, we developed phiC31 Integrase Genomic Loci Engineered for Transgenesis (pIGLET) as transgenesis approach, with two attP landing sites pIGLET14a and pIGLET24b from well-validated Tol2 transgenes. Both sites facilitate diverse transgenesis applications including reporters and Cre/loxP transgenes. The pIGLET14a and pIGLET24b landing sites consistently yield 25 to 50% germline transmission, substantially reducing the resources needed for transgenic line generation. Transgenesis into these sites enables reproducible expression patterns in F0 zebrafish embryos for enhancer discovery and testing of gene regulatory variants. Together, our new landing sites streamline targeted, reproducible zebrafish transgenesis as a robust platform for various applications while minimizing the workload for generating transgenic lines.

A nox2/cybb zebrafish mutant with defective myeloid cell reactive oxygen species production displays normal initial neutrophil recruitment to sterile tail injuries

Reactive oxygen species are important effectors and modifiers of the acute inflammatory response, recruiting phagocytes including neutrophils to sites of tissue injury. In turn, phagocytes such as neutrophils are both consumers and producers of reactive oxygen species. Phagocytes including neutrophils generate reactive oxygen species in an oxidative burst through the activity of a multimeric phagocytic nicotinamide adenine dinucleotide phosphate oxidase complex. Mutations in the NOX2/CYBB (previously gp91phox) nicotinamide adenine dinucleotide phosphate oxidase subunit are the commonest cause of chronic granulomatous disease, a disease characterized by infection susceptibility and an inflammatory phenotype. To model chronic granulomatous disease, we made a nox2/cybb zebrafish (Danio rerio) mutant and demonstrated it to have severely impaired myeloid cell reactive oxygen species production. Reduced early survival of nox2 mutant embryos indicated an essential requirement for nox2 during early development. In nox2/cybb zebrafish mutants, the dynamics of initial neutrophil recruitment to both mild and severe surgical tailfin wounds was normal, suggesting that excessive neutrophil recruitment at the initiation of inflammation is not the primary cause of the "sterile" inflammatory phenotype of chronic granulomatous disease patients. This nox2 zebrafish mutant adds to existing in vivo models for studying reactive oxygen species function in myeloid cells including neutrophils in development and disease.

Evaluating the toxicity of estetrol, 17α-ethinylestradiol, and their combination with drospirenone on zebrafish larvae: A behavioural and proteomic study

OBJECTIVE

To characterise and compare the toxicity of estetrol (E4) and 17α-ethinylestradiol (EE2), and their respective mixture with the progestin drospirenone (DRSP) in zebrafish (Danio rerio) embryos.

METHODS

Zebrafish embryos were exposed to E4, EE2, DRSP, E4+DRSP, and EE2+DRSP in a fish embryo acute toxicity (FET) test. A second test examined behavioural responses and, using label-free proteomics, identified changes in protein expression in response to hormonal treatments, across a range of concentrations, including those that are considered to be environmentally relevant.

RESULTS

In the FET test, no effects were found from E4 at concentrations ≤100 mg/L, while EE2 induced mortality and morphological abnormalities at concentrations of 1-2 mg/L. In the behavioural test, exposure to 30 ng/L EE2 (∼200 × predicted environmental concentration - PEC) resulted in hypoactivity in fish larvae and exposure to 0.3 ng/L EE2 (∼2 × PEC) led to quantitative changes in protein abundance, revealing potential impacts on RNA processing and protein synthesis machinery. Exposure to E4 did not alter behaviour, but several groups of proteins were modulated, mainly at 710 ng/L (∼200 × PEC), including proteins involved in oxidative phosphorylation. When combined with DRSP, EE2 induced reduced effects on behaviour and proteomic responses, suggesting an antagonistic effect of DRSP. E4+DRSP induced no significant effects on behaviour or proteomic profiles at tested concentrations.

CONCLUSIONS

These findings suggest that E4-based combined oral contraceptives present a more favourable environmental profile than EE2-based contraceptives, particularly during the early developmental stages of fish.

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

CRISPR-RfxCas13d screening uncovers Bckdk as a post-translational regulator of the maternal-to-zygotic transition in teleosts

The Maternal-to-Zygotic transition (MZT) is a reprograming process encompassing zygotic genome activation (ZGA) and the clearance of maternally-provided mRNAs. While some factors regulating MZT have been identified, there are thousands of maternal RNAs whose function has not been ascribed yet. Here, we have performed a proof-of-principle CRISPR-RfxCas13d maternal screening targeting mRNAs encoding protein kinases and phosphatases in zebrafish and identified Bckdk as a novel post-translational regulator of MZT. Bckdk mRNA knockdown caused epiboly defects, ZGA deregulation, H3K27ac reduction and a partial impairment of miR-430 processing. Phospho-proteomic analysis revealed that Phf10/Baf45a, a chromatin remodeling factor, is less phosphorylated upon Bckdk depletion. Further, phf10 mRNA knockdown also altered ZGA and Phf10 constitutively phosphorylated rescued the developmental defects observed after bckdk mRNA depletion. Altogether, our results demonstrate the competence of CRISPR-RfxCas13d screenings to uncover new regulators of early vertebrate development and shed light on the post-translational control of MZT mediated by protein phosphorylation.

Mouse Genome Informatics: an integrated knowledgebase system for the laboratory mouse

Mouse Genome Informatics (MGI) is a federation of expertly curated information resources designed to support experimental and computational investigations into genetic and genomic aspects of human biology and disease using the laboratory mouse as a model system. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are core MGI databases that share data and system architecture. MGI serves as the central community resource of integrated information about mouse genome features, variation, expression, gene function, phenotype, and human disease models acquired from peer-reviewed publications, author submissions, and major bioinformatics resources. To facilitate integration and standardization of data, biocuration scientists annotate using terms from controlled metadata vocabularies and biological ontologies (e.g. Mammalian Phenotype Ontology, Mouse Developmental Anatomy, Disease Ontology, Gene Ontology, etc.), and by applying international community standards for gene, allele, and mouse strain nomenclature. MGI serves basic scientists, translational researchers, and data scientists by providing access to FAIR-compliant data in both human-readable and compute-ready formats. The MGI resource is accessible at https://informatics.jax.org. Here, we present an overview of the core data types represented in MGI and highlight recent enhancements to the resource with a focus on new data and functionality for MGD and GXD.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

danRerLib: a Python package for zebrafish transcriptomics

Summary

Understanding the pathways and biological processes underlying differential gene expression is fundamental for characterizing gene expression changes in response to an experimental condition. Zebrafish, with a transcriptome closely mirroring that of humans, are frequently utilized as a model for human development and disease. However, a challenge arises due to the incomplete annotations of zebrafish pathways and biological processes, with more comprehensive annotations existing in humans. This incompleteness may result in biased functional enrichment findings and loss of knowledge. danRerLib, a versatile Python package for zebrafish transcriptomics researchers, overcomes this challenge and provides a suite of tools to be executed in Python including gene ID mapping, orthology mapping for the zebrafish and human taxonomy, and functional enrichment analysis utilizing the latest updated Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. danRerLib enables functional enrichment analysis for GO and KEGG pathways, even when they lack direct zebrafish annotations through the orthology of human-annotated functional annotations. This approach enables researchers to extend their analysis to a wider range of pathways, elucidating additional mechanisms of interest and greater insight into experimental results.

Availability and implementation

danRerLib, along with comprehensive documentation and tutorials, is freely available. The source code is available at https://github.com/sdsucomptox/danrerlib/ with associated documentation and tutorials at https://sdsucomptox.github.io/danrerlib/. The package has been developed with Python 3.9 and is available for installation on the package management systems PIP (https://pypi.org/project/danrerlib/) and Conda (https://anaconda.org/sdsu_comptox/danrerlib) with additional installation instructions on the documentation website.

Zebrafish nampt-a mutants are viable despite perturbed primitive hematopoiesis

BACKGROUND

Nicotinamide phosphoribosyltransferase (Nampt) is required for recycling NAD+ in numerous cellular contexts. Morpholino-based knockdown of zebrafish nampt-a has been shown to cause abnormal development and defective hematopoiesis concomitant with decreased NAD+ levels. However, surprisingly, nampt-a mutant zebrafish were recently found to be viable, suggesting a discrepancy between the phenotypes in knockdown and knockout conditions. Here, we address this discrepancy by directly comparing loss-of-function approaches that result in identical defective transcripts in morphants and mutants.

RESULTS

Using CRISPR/Cas9-mediated mutagenesis, we generated nampt-a mutant lines that carry the same mis-spliced mRNA as nampt-a morphants. Despite reduced NAD+ levels and perturbed expression of specific blood markers, nampt-a mutants did not display obvious developmental defects and were found to be viable. In contrast, injection of nampt-a morpholinos into wild-type or mutant nampt-a embryos caused aberrant phenotypes. Moreover, nampt-a morpholinos caused additional reduction of blood-related markers in nampt-a mutants, suggesting that the defects observed in nampt-a morphants can be partially attributed to off-target effects of the morpholinos.

CONCLUSIONS

Our findings show that zebrafish nampt-a mutants are viable despite reduced NAD+ levels and a perturbed hematopoietic gene expression program, indicating strong robustness of primitive hematopoiesis during early embryogenesis.

Shedding a Light on Dark Genes: A Comparative Expression Study of PRR12 Orthologues during Zebrafish Development

Haploinsufficiency of the PRR12 gene is implicated in a human neuro-ocular syndrome. Although identified as a nuclear protein highly expressed in the embryonic mouse brain, PRR12 molecular function remains elusive. This study explores the spatio-temporal expression of zebrafish PRR12 co-orthologs, prr12a and prr12b, as a first step to elucidate their function. In silico analysis reveals high evolutionary conservation in the DNA-interacting domains for both orthologs, with significant syntenic conservation observed for the prr12b locus. In situ hybridization and RT-qPCR analyses on zebrafish embryos and larvae reveal distinct expression patterns: prr12a is expressed early in zygotic development, mainly in the central nervous system, while prr12b expression initiates during gastrulation, localizing later to dopaminergic telencephalic and diencephalic cell clusters. Both transcripts are enriched in the ganglion cell and inner neural layers of the 72 hpf retina, with prr12b widely distributed in the ciliary marginal zone. In the adult brain, prr12a and prr12b are found in the cerebellum, amygdala and ventral telencephalon, which represent the main areas affected in autistic patients. Overall, this study suggests PRR12's potential involvement in eye and brain development, laying the groundwork for further investigations into PRR12-related neurobehavioral disorders.

Gluconeogenesis in the extraembryonic yolk syncytial layer of the zebrafish embryo

Yolk-consuming (lecithotrophic) embryos of oviparous animals, such as those of fish, need to make do with the maternally derived yolk. However, in many cases, yolk possesses little carbohydrates and sugars, including glucose, the essential monosaccharide. Interestingly, increases in the glucose content were found in embryos of some teleost fishes; however, the origin of this glucose has been unknown. Unveiling new metabolic strategies in fish embryos has a potential for better aquaculture technologies. In the present study, using zebrafish, we assessed how these embryos obtain the glucose. We employed stable isotope (13C)-labeled substrates and injected them to the zebrafish embryos. Our liquid chromatography-mass spectrometry-based isotope tracking revealed that among all tested substrate, glutamate was most actively metabolized to produce glucose in the zebrafish embryos. Expression analysis for gluconeogenic genes found that many of these were expressed in the yolk syncytial layer (YSL), an extraembryonic tissue found in teleost fishes. Generation 0 (G0) knockout of pck2, a gene encoding the key enzyme for gluconeogenesis from Krebs cycle intermediates, reduced gluconeogenesis from glutamate, suggesting that this gene is responsible for gluconeogenesis from glutamate in the zebrafish embryos. These results showed that teleost YSL undergoes gluconeogenesis, likely contributing to the glucose supplementation to the embryos with limited glucose source. Since many other animal lineages lack YSL, further comparative analysis will be interesting.

A Foxf1-Wnt-Nr2f1 cascade promotes atrial cardiomyocyte differentiation in zebrafish

Nr2f transcription factors (TFs) are conserved regulators of vertebrate atrial cardiomyocyte (AC) differentiation. However, little is known about the mechanisms directing Nr2f expression in ACs. Here, we identified a conserved enhancer 3' to the nr2f1a locus, which we call 3'reg1-nr2f1a (3'reg1), that can promote Nr2f1a expression in ACs. Sequence analysis of the enhancer identified putative Lef/Tcf and Foxf TF binding sites. Mutation of the Lef/Tcf sites within the 3'reg1 reporter, knockdown of Tcf7l1a, and manipulation of canonical Wnt signaling support that Tcf7l1a is derepressed via Wnt signaling to activate the transgenic enhancer and promote AC differentiation. Similarly, mutation of the Foxf binding sites in the 3'reg1 reporter, coupled with gain- and loss-of-function analysis supported that Foxf1 promotes expression of the enhancer and AC differentiation. Functionally, we find that Wnt signaling acts downstream of Foxf1 to promote expression of the 3'reg1 reporter within ACs and, importantly, both Foxf1 and Wnt signaling require Nr2f1a to promote a surplus of differentiated ACs. CRISPR-mediated deletion of the endogenous 3'reg1 abrogates the ability of Foxf1 and Wnt signaling to produce surplus ACs in zebrafish embryos. Together, our data support that downstream members of a conserved regulatory network involving Wnt signaling and Foxf1 function on a nr2f1a enhancer to promote AC differentiation in the zebrafish heart.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

fhl2b mediates extraocular muscle protection in zebrafish models of muscular dystrophies and its ectopic expression ameliorates affected body muscles

In muscular dystrophies, muscle fibers loose integrity and die, causing significant suffering and premature death. Strikingly, the extraocular muscles (EOMs) are spared, functioning well despite the disease progression. Although EOMs have been shown to differ from body musculature, the mechanisms underlying this inherent resistance to muscle dystrophies remain unknown. Here, we demonstrate important differences in gene expression as a response to muscle dystrophies between the EOMs and trunk muscles in zebrafish via transcriptomic profiling. We show that the LIM-protein Fhl2 is increased in response to the knockout of desmin, plectin and obscurin, cytoskeletal proteins whose knockout causes different muscle dystrophies, and contributes to disease protection of the EOMs. Moreover, we show that ectopic expression of fhl2b can partially rescue the muscle phenotype in the zebrafish Duchenne muscular dystrophy model sapje, significantly improving their survival. Therefore, Fhl2 is a protective agent and a candidate target gene for therapy of muscular dystrophies.

Paternal starvation affects metabolic gene expression during zebrafish offspring development and lifelong fitness

Dietary restriction in the form of fasting is a putative key to a healthier and longer life, but these benefits may come at a trade-off with reproductive fitness and may affect the following generation(s). The potential inter- and transgenerational effects of long-term fasting and starvation are particularly poorly understood in vertebrates when they originate from the paternal line. We utilised the externally fertilising zebrafish amenable to a split-egg clutch design to explore the male-specific effects of fasting/starvation on fertility and fitness of offspring independently of maternal contribution. Eighteen days of fasting resulted in reduced fertility in exposed males. While average offspring survival was not affected, we detected increased larval growth rate in F1 offspring from starved males and more malformed embryos at 24 h post-fertilisation in F2 offspring produced by F1 offspring from starved males. Comparing the transcriptomes of F1 embryos sired by starved and fed fathers revealed robust and reproducible increased expression of muscle composition genes but lower expression of lipid metabolism and lysosome genes in embryos from starved fathers. A large proportion of these genes showed enrichment in the yolk syncytial layer suggesting gene regulatory responses associated with metabolism of nutrients through paternal effects on extra-embryonic tissues which are loaded with maternal factors. We compared the embryo transcriptomes to published adult transcriptome datasets and found comparable repressive effects of starvation on metabolism-associated genes. These similarities suggest a physiologically relevant, directed and potentially adaptive response transmitted by the father, independently from the offspring's nutritional state, which was defined by the mother.

Molecular interactions of paraben family of pollutants with embryonic neuronal proteins of Danio rerio: A step ahead in computational toxicity towards adverse outcome pathway

The paraben family of endocrine disruptors exhibit persistent behaviours in aquatic matrices, having bio-accumulative effects and necessitating toxicity analysis and safe use, as well as prevention of food web penetration. In this study, the toxicity effects of 9 different parabens (Methyl, Ethyl, Propyl, Butyl, Heptyl, Isopropyl, Isobutyl, benzyl parabens and p-hydroxybenzoic acid) were studied against 17 neuronal proteins (Neurog1, Ascl1a, DLA, Syn2a, Ntn1a, Pitx2, and SoxB1, Her/Hes, Zic family) expressed during the early embryonic developmental stage of Danio rerio. The neuronal genes were selected as a biomarker to study the inhibitory effects on the cascade of genes expressed in the early developmental stage. The study uses trRossetta software to predict protein structures of neuronal genes, followed by structural refinement, energy minimisation, and active site prediction, evaluated using energy value, RC plot and ERRAT scores of PROCHECK and ERRAT programs. Compared to raw structures, highly confident predicted structures and quality scores were observed for refined protein with few exceptions. Based on the polarity and charge of the aminoacids, the probable pockets were identified using active site prediction, which were then used for molecular docking analysis. Further, the ADMET analysis, ligand likeliness and toxicological test revealed the paraben family of compounds as one of the most susceptible toxic and mutagenic compounds. The molecular docking results showed an interesting pattern of increasing binding affinity with increase in the carbon chains of paraben molecules. Benzyl Paraben showed higher binding affinities across all 17 neuronal proteins. Finally, gene co-occurrence/co-expression and protein-protein interaction studies using the STRING database depict that all proteins are functionally related and play essential roles in standard biological processes or pathways, conserved and expressed in diverse organisms. The interaction between paraben compounds and neuronal genes indicates high risks of inhibiting reactions in embryonic stages, emphasising the need for effective treatment measures and strict regulations.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

Genetic regulation of enteric nervous system development in zebrafish

The enteric nervous system (ENS) is a complex series of interconnected neurons and glia that reside within and along the entire length of the gastrointestinal tract. ENS functions are vital to gut homeostasis and digestion, including local control of peristalsis, water balance, and intestinal cell barrier function. How the ENS develops during embryological development is a topic of great concern, as defects in ENS development can result in various diseases, the most common being Hirschsprung disease, in which variable regions of the infant gut lack ENS, with the distal colon most affected. Deciphering how the ENS forms from its progenitor cells, enteric neural crest cells, is an active area of research across various animal models. The vertebrate animal model, zebrafish, has been increasingly leveraged to understand early ENS formation, and over the past 20 years has contributed to our knowledge of the genetic regulation that underlies enteric development. In this review, I summarize our knowledge regarding the genetic regulation of zebrafish enteric neuronal development, and based on the most current literature, present a gene regulatory network inferred to underlie its construction. I also provide perspectives on areas for future zebrafish ENS research.

Cis-regulatory control of mammalian Trps1 gene expression

TRPS1 serves as the causative gene for tricho-rhino phalangeal syndrome, known for its craniofacial and skeletal abnormalities. The Trps1 gene encodes a protein that represses Wnt signaling through strong interactions with Wnt signaling inhibitors. The identification of genomic cis-acting regulatory sequences governing Trps1 expression is crucial for understanding its role in embryogenesis. Nevertheless, to date, no investigations have been conducted concerning these aspects of Trps1. To identify deeply conserved noncoding elements (CNEs) within the Trps1 locus, we employed a comparative genomics approach, utilizing slowly evolving fish such as coelacanth and spotted gar. These analyses resulted in the identification of eight CNEs in the intronic region of the Trps1 gene. Functional characterization of these CNEs in zebrafish revealed their regulatory potential in various tissues, including pectoral fins, heart, and pharyngeal arches. RNA in-situ hybridization experiments revealed concordance between the reporter expression pattern induced by the identified set of CNEs and the spatial expression pattern of the trps1 gene in zebrafish. Comparative in vivo data from zebrafish and mice for CNE7/hs919 revealed conserved functions of these enhancers. Each of these eight CNEs was further investigated in cell line-based reporter assays, revealing their repressive potential. Taken together, in vivo and in vitro assays suggest a context-dependent dual functionality for the identified set of Trps1-associated CNE enhancers. This functionally characterized set of CNE-enhancers will contribute to a more comprehensive understanding of the developmental roles of Trps1 and can aid in the identification of noncoding DNA variants associated with human diseases.

Pluripotency of a founding field: rebranding developmental biology

The field of developmental biology has declined in prominence in recent decades, with off-shoots from the field becoming more fashionable and highly funded. This has created inequity in discovery and opportunity, partly due to the perception that the field is antiquated or not cutting edge. A 'think tank' of scientists from multiple developmental biology-related disciplines came together to define specific challenges in the field that may have inhibited innovation, and to provide tangible solutions to some of the issues facing developmental biology. The community suggestions include a call to the community to help 'rebrand' the field, alongside proposals for additional funding apparatuses, frameworks for interdisciplinary innovative collaborations, pedagogical access, improved science communication, increased diversity and inclusion, and equity of resources to provide maximal impact to the community.

Plxnd1-mediated mechanosensing of blood flow controls the caliber of the Dorsal Aorta via the transcription factor Klf2

The cardiovascular system generates and responds to mechanical forces. The heartbeat pumps blood through a network of vascular tubes, which adjust their caliber in response to the hemodynamic environment. However, how endothelial cells in the developing vascular system integrate inputs from circulatory forces into signaling pathways to define vessel caliber is poorly understood. Using vertebrate embryos and in vitro-assembled microvascular networks of human endothelial cells as models, flow and genetic manipulations, and custom software, we reveal that Plexin-D1, an endothelial Semaphorin receptor critical for angiogenic guidance, employs its mechanosensing activity to serve as a crucial positive regulator of the Dorsal Aorta's (DA) caliber. We also uncover that the flow-responsive transcription factor KLF2 acts as a paramount mechanosensitive effector of Plexin-D1 that enlarges endothelial cells to widen the vessel. These findings illuminate the molecular and cellular mechanisms orchestrating the interplay between cardiovascular development and hemodynamic forces.

Environmental and molecular control of tissue-specific ionocyte differentiation in zebrafish

Organisms adjust their physiology to cope with environmental fluctuations and maintain fitness. These adaptations occur via genetic changes over multiple generations or through acclimation, a set of reversible phenotypic changes that confer resilience to the individual. Aquatic organisms are subject to dramatic seasonal fluctuations in water salinity, which can affect the function of lateral line mechanosensory hair cells. To maintain hair cell function when salinity decreases, ion-regulating cells, Neuromast-associated ionocytes (Nm ionocytes), increase in number and invade lateral line neuromasts. How environmental changes trigger this adaptive differentiation of Nm ionocytes and how these cells are specified is still unknown. Here, we identify Nm ionocyte progenitors as foxi3a/foxi3b-expressing skin cells and show that their differentiation is associated with sequential activation of different Notch pathway components, which control ionocyte survival. We demonstrate that new Nm ionocytes are rapidly specified by absolute salinity levels, independently of stress response pathways. We further show that Nm ionocyte differentiation is selectively triggered by depletion of specific ions, such as Ca2+ and Na+/Cl-, but not by low K+ levels, and is independent of media osmolarity. Finally, we demonstrate that hair cell activity plays a role in Nm ionocyte recruitment and that systemic factors are not necessary for Nm ionocyte induction. In summary, we have identified how environmental changes activate a signaling cascade that triggers basal skin cell progenitors to differentiate into Nm ionocytes and invade lateral line organs. This adaptive behavior is an example of physiological plasticity that may prove essential for survival in changing climates.

The Monarch Initiative in 2024: an analytic platform integrating phenotypes, genes and diseases across species

Bridging the gap between genetic variations, environmental determinants, and phenotypic outcomes is critical for supporting clinical diagnosis and understanding mechanisms of diseases. It requires integrating open data at a global scale. The Monarch Initiative advances these goals by developing open ontologies, semantic data models, and knowledge graphs for translational research. The Monarch App is an integrated platform combining data about genes, phenotypes, and diseases across species. Monarch's APIs enable access to carefully curated datasets and advanced analysis tools that support the understanding and diagnosis of disease for diverse applications such as variant prioritization, deep phenotyping, and patient profile-matching. We have migrated our system into a scalable, cloud-based infrastructure; simplified Monarch's data ingestion and knowledge graph integration systems; enhanced data mapping and integration standards; and developed a new user interface with novel search and graph navigation features. Furthermore, we advanced Monarch's analytic tools by developing a customized plugin for OpenAI's ChatGPT to increase the reliability of its responses about phenotypic data, allowing us to interrogate the knowledge in the Monarch graph using state-of-the-art Large Language Models. The resources of the Monarch Initiative can be found at monarchinitiative.org and its corresponding code repository at github.com/monarch-initiative/monarch-app.

A3D Model Organism Database (A3D-MODB): a database for proteome aggregation predictions in model organisms

Protein aggregation has been associated with aging and different pathologies and represents a bottleneck in the industrial production of biotherapeutics. Numerous past studies performed in Escherichia coli and other model organisms have allowed to dissect the biophysical principles underlying this process. This knowledge fuelled the development of computational tools, such as Aggrescan 3D (A3D) to forecast and re-design protein aggregation. Here, we present the A3D Model Organism Database (A3D-MODB) http://biocomp.chem.uw.edu.pl/A3D2/MODB, a comprehensive resource for the study of structural protein aggregation in the proteomes of 12 key model species spanning distant biological clades. In addition to A3D predictions, this resource incorporates information useful for contextualizing protein aggregation, including membrane protein topology and structural model confidence, as an indirect reporter of protein disorder. The database is openly accessible without any need for registration. We foresee A3D-MOBD evolving into a central hub for conducting comprehensive, multi-species analyses of protein aggregation, fostering the development of protein-based solutions for medical, biotechnological, agricultural and industrial applications.

Metabolomic analysis of Drosophila melanogaster larvae lacking pyruvate kinase

Pyruvate kinase (Pyk) is a rate-limiting enzyme that catalyzes the final metabolic reaction in glycolysis. The importance of this enzyme, however, extends far beyond ATP production, as Pyk is also known to regulate tissue growth, cell proliferation, and development. Studies of this enzyme in Drosophila melanogaster are complicated by the fact that the fly genome encodes 6 Pyk paralogs whose functions remain poorly defined. To address this issue, we used sequence distance and phylogenetic approaches to demonstrate that the gene Pyk encodes the enzyme most similar to the mammalian Pyk orthologs, while the other 5 Drosophila Pyk paralogs have significantly diverged from the canonical enzyme. Consistent with this observation, metabolomic studies of 2 different Pyk mutant strains revealed that larvae lacking Pyk exhibit a severe block in glycolysis, with a buildup of glycolytic intermediates upstream of pyruvate. However, our analysis also unexpectedly reveals that pyruvate levels are unchanged in Pyk mutants, indicating that larval metabolism maintains pyruvate pool size despite severe metabolic limitations. Consistent with our metabolomic findings, a complementary RNA-seq analysis revealed that genes involved in lipid metabolism and protease activity are elevated in Pyk mutants, again indicating that loss of this glycolytic enzyme induces compensatory changes in other aspects of metabolism. Overall, our study provides both insight into how Drosophila larval metabolism adapts to disruption of glycolytic metabolism as well as immediate clinical relevance, considering that Pyk deficiency is the most common congenital enzymatic defect in humans.

Interpretable trajectory inference with single-cell Linear Adaptive Negative-binomial Expression (scLANE) testing

The rapid proliferation of trajectory inference methods for single-cell RNA-seq data has allowed researchers to investigate complex biological processes by examining underlying gene expression dynamics. After estimating a latent cell ordering, statistical models are used to determine which genes exhibit changes in expression that are significantly associated with progression through the biological trajectory. While a few techniques for performing trajectory differential expression exist, most rely on the flexibility of generalized additive models in order to account for the inherent nonlinearity of changes in gene expression. As such, the results can be difficult to interpret, and biological conclusions often rest on subjective visual inspections of the most dynamic genes. To address this challenge, we propose scLANE testing, which is built around an interpretable generalized linear model and handles nonlinearity with basis splines chosen empirically for each gene. In addition, extensions to estimating equations and mixed models allow for reliable trajectory testing under complex experimental designs. After validating the accuracy of scLANE under several different simulation scenarios, we apply it to a set of diverse biological datasets and display its ability to provide novel biological information when used downstream of both pseudotime and RNA velocity estimation methods.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Loss of Panx1 function in zebrafish alters motor behavior in a lab-on-chip model of Parkinson's disease

Pannexin 1 (Panx1) forms ATP-permeable membrane channels that play roles in purinergic signaling in the nervous system. A link between Panx1 activity and neurodegenerative disorders including Parkinson's disease (PD) has been suggested, but experimental evidence is limited. Here, a zebrafish model of PD was produced by exposing panx1a+/+ and panx1a-/- zebrafish larvae to 6-hydroxydopamine (6-OHDA). Electrical stimulation in a microfluidic chip and quantitative real-time-qPCR of zebrafish larvae tested the role of Panx1 in both pathological and normal conditions. After 72-h treatment with 6-OHDA, the electric-induced locomotor activity of 5 days post fertilization (5dpf) panx1a+/+ larvae were reduced, while the stimulus did not affect locomotor activity of age-matched panx1a-/- larvae. A RT-qPCR analysis showed an increase in the expression of genes that are functionally related to dopaminergic signaling, like the tyrosine hydroxylase (th2) and the leucine-rich repeat kinase 2 (lrrk2). Extending the 6-OHDA treatment duration to 120 h caused a significant reduction in the locomotor response of 7dpf panx1a-/- larvae compared to the untreated panx1a-/- group. The RT-qPCR data showed a reduced expression of dopaminergic signaling genes in both genotypes. It was concluded that the absence of Panx1a channels compromised dopaminergic signaling in 6-OHDA-treated zebrafish larvae and that the increase in the expression of dopaminergic genes was transient, most likely due to a compensatory upregulation. We propose that zebrafish Panx1a models offer opportunities to shed light on PD's physiological and molecular basis. Panx1a might play a role on the progression of PD, and therefore deserves further investigation.

Prechordal structures act cooperatively in early trabeculae development of gnathostome skull

The vertebrate skull is formed by mesoderm and neural crest (NC) cells. The mesoderm contributes to the skull chordal domain, with the notochord playing an essential role in this process. The NC contributes to the skull prechordal domain, prompting investigation into the embryonic structures involved in prechordal neurocranium cartilage formation. The trabeculae cartilage, a structure of the prechordal neurocranium, arises at the convergence of prechordal plate (PCP), ventral midline (VM) cells of the diencephalon, and dorsal oral ectoderm. This study examines the molecular participation of these embryonic structures in gnathostome trabeculae development. PCP-secreted SHH induces its expression in VM cells of the diencephalon, initiating a positive feedback loop involving SIX3 and GLI1. SHH secreted by the VM cells of the diencephalon acts on the dorsal oral ectoderm, stimulating condensation of NC cells to form trabeculae. SHH from the prechordal region affects the expression of SOX9 in NC cells. BMP7 and SHH secreted by PCP induce NKX2.1 expression in VM cells of the diencephalon, but this does not impact trabeculae formation. Molecular cooperation between PCP, VM cells of the diencephalon, and dorsal oral ectoderm is crucial for craniofacial development by NC cells in the prechordal domain.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Single-cell RNAseq analysis of spinal locomotor circuitry in larval zebrafish

Identification of the neuronal types that form the specialized circuits controlling distinct behaviors has benefited greatly from the simplicity offered by zebrafish. Electrophysiological studies have shown that in addition to connectivity, understanding of circuitry requires identification of functional specializations among individual circuit components, such as those that regulate levels of transmitter release and neuronal excitability. In this study, we use single-cell RNA sequencing (scRNAseq) to identify the molecular bases for functional distinctions between motoneuron types that are causal to their differential roles in swimming. The primary motoneuron, in particular, expresses high levels of a unique combination of voltage-dependent ion channel types and synaptic proteins termed functional 'cassettes.' The ion channel types are specialized for promoting high-frequency firing of action potentials and augmented transmitter release at the neuromuscular junction, both contributing to greater power generation. Our transcriptional profiling of spinal neurons further assigns expression of this cassette to specific interneuron types also involved in the central circuitry controlling high-speed swimming and escape behaviors. Our analysis highlights the utility of scRNAseq in functional characterization of neuronal circuitry, in addition to providing a gene expression resource for studying cell type diversity.

scMultiome analysis identifies embryonic hindbrain progenitors with mixed rhombomere identities

Rhombomeres serve to position neural progenitors in the embryonic hindbrain, thereby ensuring appropriate neural circuit formation, but the molecular identities of individual rhombomeres and the mechanism whereby they form has not been fully established. Here, we apply scMultiome analysis in zebrafish to molecularly resolve all rhombomeres for the first time. We find that rhombomeres become molecularly distinct between 10hpf (end of gastrulation) and 13hpf (early segmentation). While the embryonic hindbrain transiently contains alternating odd- versus even-type rhombomeres, our scMultiome analyses do not detect extensive odd versus even molecular characteristics in the early hindbrain. Instead, we find that each rhombomere displays a unique gene expression and chromatin profile. Prior to the appearance of distinct rhombomeres, we detect three hindbrain progenitor clusters (PHPDs) that correlate with the earliest visually observed segments in the hindbrain primordium that represent prospective rhombomere r2/r3 (possibly including r1), r4, and r5/r6, respectively. We further find that the PHPDs form in response to Fgf and RA morphogens and that individual PHPD cells co-express markers of multiple mature rhombomeres. We propose that the PHPDs contain mixed-identity progenitors and that their subdivision into individual rhombomeres requires the resolution of mixed transcription and chromatin states.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

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.

Aon: a service to augment Alliance Genome Resource data with additional species

OBJECTIVE

Cross-species comparative genomics requires access to accurate homology data across the entire range of annotated genes. The Alliance of Genome Resources (AGR) provides an open-source and comprehensive database of homology data calculated using a wide array of algorithms at differing stringencies to elucidate orthologous relationships. However, the current AGR application program interface (API) is limited to five homology endpoints for nine species. While AGR provides a robust resource for several canonical species, its utility can be greatly enhanced by increased filtering and data processing options and incorporating additional species.

RESULTS

Here, we describe a novel API tool, AON, that expands access to the AGR orthology resource by creating a data structure that supports 50 additional endpoints. More importantly, it provides users with a framework for adding bespoke endpoints, custom species, and additional orthology data. We demonstrate AON's functionality by incorporating the service into the GeneWeaver ecosystem for supporting cross-species data analysis.

SGMS2 in primary osteoporosis with facial nerve palsy

Pathogenic heterozygous variants in SGMS2 cause a rare monogenic form of osteoporosis known as calvarial doughnut lesions with bone fragility (CDL). The clinical presentations of SGMS2-related bone pathology range from childhood-onset osteoporosis with low bone mineral density and sclerotic doughnut-shaped lesions in the skull to a severe spondylometaphyseal dysplasia with neonatal fractures, long-bone deformities, and short stature. In addition, neurological manifestations occur in some patients. SGMS2 encodes sphingomyelin synthase 2 (SMS2), an enzyme involved in the production of sphingomyelin (SM). This review describes the biochemical structure of SM, SM metabolism, and their molecular actions in skeletal and neural tissue. We postulate how disrupted SM gradient can influence bone formation and how animal models may facilitate a better understanding of SGMS2-related osteoporosis.

The RELT Family of Proteins: An Increasing Awareness of Their Importance for Cancer, the Immune System, and Development

This review highlights Receptor Expressed in Lymphoid Tissues (RELT), a Tumor Necrosis Factor Superfamily member, and its two paralogs, RELL1 and RELL2. Collectively, these three proteins are referred to as RELTfms and have gained much interest in recent years due to their association with cancer and other human diseases. A thorough knowledge of their physiological functions, including the ligand for RELT, is lacking, yet emerging evidence implicates RELTfms in a variety of processes including cytokine signaling and pathways that either promote cell death or survival. T cells from mice lacking RELT exhibit increased responses against tumors and increased inflammatory cytokine production, and multiple lines of evidence indicate that RELT may promote an immunosuppressive environment for tumors. The relationship of individual RELTfms in different cancers is not universal however, as evidence indicates that individual RELTfms may be risk factors in certain cancers yet appear to be protective in other cancers. RELTfms are important for a variety of additional processes related to human health including microbial pathogenesis, inflammation, behavior, reproduction, and development. All three proteins have been strongly conserved in all vertebrates, and this review aims to provide a clearer understanding of the current knowledge regarding these interesting proteins.