Water pollution remains a major environmental concern, with increased toxic by-products being released into water bodies. Many of these chemical contaminants persist in the environment and bio-accumulate in aquatic organisms. At present, toxicological tests are mostly based on laboratory tests, and effective methods for monitoring wild aquatic environments remain lacking. In the present study, we used a well-characterized toxic chemical, 3,3\',4,4\',5-polychlorinated biphenyl (PCB126), as an example to try to identify common biomarker genes to be used for predictive toxicity of this toxic substance. First, we used two laboratory fish models, the zebrafish (Danio rerio) and medaka (Oryzias latipes), to expose PCB126 to obtain liver transcriptomic data by RNA-seq. Comparative transcriptomic analyses indicated generally conserved and concerted changes from the two species, thus validating the transcriptomic data for biomarker gene selection. Based on the common up- and downregulated genes in the two species, we selected nine biomarker genes to further test in other fish species. The first validation experiment was carried out using the third fish species, Mozambique tilapia (Oreochromis mossambicus), and essentially, all these biomarker genes were validated for consistent responses with the two laboratory fish models. Finally, to develop universal PCR primers suitable for potentially all teleost fish species, we designed degenerate primers and tested them in the three fish species as well as in another fish species without a genomic sequence available: guppy (Poecilia reticulata). We found all the biomarker genes showed consistent response to PCB126 exposure in at least 50% of the species. Thus, our study provides a promising strategy to identify common biomarker genes to be used for teleost fish analyses. By using degenerate PCR primers and analyzing multiple biomarker genes, it is possible to develop diagnostic PCR arrays to predict water contamination from any wild fish species sampled in different water bodies.

The marine medaka Oryzias melastigma is a useful fish model for marine and estuarine ecotoxicology studies and can be applied to field-based population genomics because of its distribution in Asian estuaries and other coastal areas. We identified 769 full-length G protein-coupled receptor genes in the O. melastigma genome and classified them into five distinct classes. A phylogenetic comparison of GPCR genes in O. melastigma to humans and two other small fish species revealed a high-level orthological relationship. Purinergic and chemokine receptors were highly differentiated in humans whereas significant differentiation of chemosensory receptors was evident in fish species. Our results suggest that the GPCR gene families among the species used in this study exhibit evidence of sporadic evolutionary processes. These results may help improve our understanding of the advanced repertoires of GPCR and expand our knowledge of physiological mechanisms of fish in response to various environmental stimuli.

The progression of fatty liver disease to non-alcoholic steatohepatitis (NASH) is a leading cause of death in humans. Lifestyles and environmental chemical exposures can increase the susceptibility of humans to NASH. In humans, the presence of bisphenol A (BPA) in urine is associated with fatty liver disease, but whether ancestral BPA exposure leads to the activation of human NAFLD-NASH-associated genes in the unexposed descendants is unclear. In this study, using medaka fish as an animal model for human NAFLD, we investigated the transcriptional signatures of human NAFLD-NASH and their associated roles in the pathogenesis of the liver of fish that were not directly exposed, but their ancestors were exposed to BPA during embryonic and perinatal development three generations prior. Comparison of bulk RNA-Seq data of the liver in BPA lineage male and female medaka with publicly available human NAFLD-NASH patient data revealed transgenerational alterations in the transcriptional signature of human NAFLD-NASH in medaka liver. Twenty percent of differentially expressed genes (DEGs) were upregulated in both human NAFLD patients and medaka. Specifically in females, among the total shared DEGs in the liver of BPA lineage fish and NAFLD patient groups, 27.69% were downregulated, and 20% were upregulated. Of all DEGs, 52.31% of DEGs were found in ancestral BPA-lineage females, suggesting that NAFLD in females shared the majority of human NAFLD gene networks. Pathway analysis revealed beta-oxidation, lipoprotein metabolism, and HDL/LDL-mediated transport processes linked to downregulated DEGs in BPA lineage males and females. In contrast, the expression of genes encoding lipogenesis-related proteins was significantly elevated in the liver of BPA lineage females only. BPA lineage females exhibiting activation of myc, atf4, xbp1, stat4, and cancerous pathways, as well as inactivation of igf1, suggest their possible association with an advanced NAFLD phenotype. The present results suggest that gene networks involved in the progression of human NAFLD and the transgenerational NAFLD in medaka are conserved and that medaka can be an excellent animal model to understand the development and progression of liver disease and environmental influences in the liver.

Gills and gut are the two primary osmoregulatory organs in fish. Recently, studies have expanded beyond the osmoregulatory mechanisms of these organs to explore the microbiota communities inhabiting them. It is now known that microbial communities in both organs shift in response to osmotic stress. However, there are limited studies identifying the major contributors and co-occurrence among these microbiota in both organs under seawater and freshwater transfer conditions. The current data mining report performed a bioinformatics analysis on two previous published datasets from our group, aiming to provide insights into host-bacteria relationships under osmotic stress. We divided the samples into four groups: control seawater gills (LSW); control seawater gut (TSW); freshwater transfer gills (LFW); and freshwater transfer gut (TFW). Our results showed that LSW had higher diversities, richness, and evenness compared to TSW. However, both the LFW and LSW did not show any significant differences after the freshwater transfer experiment. We further applied co-occurrence network analysis and, for the first time, reported on the interactions of taxa shaping the community structure in these two organs. Moreover, we identified enriched ectoine biosynthesis in seawater samples, suggesting its potential role in seawater environments. Increased mRNA expression levels of Na+/K+-atpase, and cftr, were observed in gills after 6 h of ectoine treatment. These findings provide a foundation for future studies on host-bacteria interactions under osmotic stress.

Bisphenol-A (BPA), a known endocrine-disrupting chemical (EDC) in plastics and resins, has been found to induce heritable health effects in fish and mammals, affecting directly exposed individuals and indirectly their progenies in subsequent generations. It is not clearly understood if subsequent generations of the BPA-exposed ancestors have increased sensitivity to the second hit by the chemicals of emerging concern. To understand this, the present study examined the effects of developmental exposure to perfluorooctanesulfonic acid (PFOS), which has been a global contaminant recently, in embryos whose ancestors were exposed to BPA. Two lineages of medaka (Oryzias latipes) were established: 1) the BPA lineage in which the F0 generation was exposed to 10 μg/L BPA during early development and 2) the control lineage with no BPA exposure in the F0 generation. These lineages were raised up to the F4 generation without further exposure. The embryos of the F4 generation were exposed to PFOS at 0, 0.002, 0.02, 0.2, 2, and 20 mg/L concentrations. Early developmental defects resulting in mortality, delayed hatching, teratogenic phenotypes, and altered gene expression were examined in both lineages. The expression level of genes encoding DNA methyltransferases and genes responsible for oxidative stress defense were determined. Following environmentally relevant PFOS exposure, organisms with a history of BPA exposure displayed significant changes in all categories of developmental defects mentioned above, including increased expression of genes related to oxidative stress, compared to individuals without BPA exposure. The present study provides initial evidence that a history of ancestral BPA exposure can alter sensitivity to developmental disorders following the second hit by PFOS exposure. The variable of ancestral BPA exposure could be considered in mechanistic, medical, and regulatory toxicology, and can also be applied to holistic environmental equity research.

Sex determination mechanisms differ widely among vertebrates, particularly in fish species, where diverse sex chromosomes and sex-determining genes have evolved. However, the sex-differentiation pathways activated by these sex-determining genes appear to be conserved. Gonadal soma-derived growth factor (Gsdf) is one of the genes conserved across teleost fish, especially in medaka fishes of the genus Oryzias, and is implicated in testis differentiation and germ cell proliferation. However, its role in sex differentiation remains unclear. In this study, we investigated Gsdf function in Oryzias hubbsi, a species with a ZW sex-determination system. We confirmed its male-dominant expression, as in other species. However, histological analyses revealed no male-to-female sex reversal in Gsdf-knockout fish, contrary to findings in other medaka species. Genetic sex determination remained intact without Gsdf function, indicating a Gsdf-independent sex-differentiation pathway in O. hubbsi. Instead, Gsdf loss led to germ cell overproliferation in both sexes and accelerated onset of meiosis in testes, suggesting a role in germ cell proliferation. Notably, the feminizing effect of germ cells observed in O. latipes was absent, suggesting diverse germ cell-somatic cell relationships in Oryzias gonad development. Our study highlights species-specific variations in the molecular pathways governing sex determination and differentiation, emphasizing the need for further exploration to elucidate the complexities of sexual development.

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.

The progression of fatty liver disease to non-alcoholic steatohepatitis (NASH) is a leading cause of death in humans. Lifestyles and environmental chemical exposures can increase the susceptibility of humans to NASH. In humans, the presence of bisphenol A (BPA) in urine is associated with fatty liver disease, but whether ancestral BPA exposure leads to the activation of human NAFLD-NASH-associated genes in the unexposed descendants is unclear. In this study, using medaka fish as an animal model for human NAFLD, we investigated the transcriptional signatures of human NAFLD-NASH and their associated roles in the pathogenesis of the liver of fish who were not directly exposed but their ancestors were exposed to BPA during embryonic and perinatal development three generations prior. Comparison of bulk RNA-Seq data of the liver in BPA lineage male and female medaka with publicly available human NAFLD-NASH patient data revealed transgenerational alterations in the transcriptional signature of human NAFLD-NASH in medaka liver. Twenty percent of differentially expressed genes (DEGs) were upregulated in both human NAFLD patients and medaka. Specifically in females, among the total shared DEGs in the liver of BPA lineage fish and NAFLD patient groups, 27.69% DEGs were downregulated and 20% DEGs were upregulated. Off all DEGs, 52.31% DEGs were found in ancestral BPA-lineage females, suggesting that NAFLD in females shared majority of human NAFLD gene networks. Pathway analysis revealed beta-oxidation, lipoprotein metabolism, and HDL/LDL-mediated transport processes linked to downregulated DEGs in BPA lineage males and females. In contrast, the expression of genes encoding lipogenesis-related proteins was significantly elevated in the liver of BPA lineage females only. BPA lineage females exhibiting activation of myc, atf4, xbp1, stat4, and cancerous pathways, as well as inactivation of igf1, suggest their possible association with an advanced NAFLD phenotype. The present results suggest that gene networks involved in the progression of human NAFLD and the transgenerational NAFLD in medaka are conserved and that medaka can be an excellent animal model to understand the development and progression of liver disease and environmental influences in the liver.

Formation of the synaptonemal complex (SC) is a prerequisite for proper recombination and chromosomal segregation during meiotic prophase I. One mechanism that ensures SC formation is chromosomal movement, which is driven by the force derived from cytoskeletal motors. Here, we report the phenotype of medaka mutants lacking the telomere repeat binding bouquet formation protein 1 (TERB1), which, in combination with the SUN/KASH protein, mediates chromosomal movement by connecting telomeres and cytoskeletal motors. Mutations in the terb1 gene exhibit defects in SC formation in medaka. Although SC formation was initiated, as seen by the punctate lateral elements and fragmented transverse filaments, it was not completed in the terb1 mutant meiocytes. The mutant phenotype further revealed that the introduction of double strand breaks was independent of synapsis completion. In association with these phenotypes, meiocytes in both the ovaries and testes exhibited an aberrant arrangement of homologous chromosomes. Interestingly, although oogenesis halted at the zygotene-like stage in terb1 mutant, testes continued to produce sperm-like cells with aberrant DNA content. This indicates that the mechanism of meiotic checkpoint is sexually different in medaka, similar to the mammalian checkpoint in which oogenesis proceeds while spermatogenesis is arrested. Moreover, our results suggest that spermatogenesis is mechanistically dissociable from meiosis.

Motilin (MLN) is a peptide hormone originally isolated from the mucosa of the porcine intestine. Its orthologs have been identified in various vertebrates. Although MLN regulates gastrointestinal motility in tetrapods from amphibians to mammals, recent studies indicate that MLN is not involved in the regulation of isolated intestinal motility in zebrafish, at least in vitro. To determine the unknown function of MLN in teleosts, we examined the expression of MLN and the MLN receptor (MLNR) at the cellular level in Japanese medaka (Oryzias latipes). Quantitative PCR revealed that mln mRNA was limitedly expressed in the gut, whereas mlnr mRNA was not detected in the gut but was expressed in the brain and kidney. By in situ hybridization and immunohistochemistry, mlnr mRNA was detected in the dopaminergic neurons of the area postrema in the brain and the noradrenaline-producing cells in the interrenal gland of the kidney. Furthermore, we observed efferent projections of mlnr-expressing dopaminergic neurons in the lobus vagi (XL) and nucleus motorius nervi vagi (NXm) of the medulla oblongata by establishing a transgenic medaka expressing the enhanced green fluorescence protein driven by the mlnr promoter. The expression of dopamine receptor mRNAs in the XL and cholinergic neurons in NXm was confirmed by in situ hybridization. These results indicate novel sites of MLN activity other than the gastrointestinal tract. MLN may exert central and peripheral actions through the regulation of catecholamine release in medaka.