Epigenetic modification and inheritance in sexual reversal of fish

Expression of genes related to gonadal development and construction of gonadal DNA methylation maps of Trachinotus blochii under hypoxia

Chronic hypoxia can affect the growth and metabolism of fish and potentially impact gonadal development through epigenetic regulation. Trachinotus blochii (Golden pompano) is widely cultured near the coast and is sensitive to low oxygen conditions. We found that hypoxia and reoxygenation processes acted on multiple targets on the HPG axis, leading to endocrine disorders. Changes in the expression of key genes in the brain (gnrh), pituitary (fsh and lh), ovaries (cyp19a1a, foxl2, and er), and testes (dmrt1, ar, sox9, and gsdf) were associated with significant decreases in estrogen and testosterone levels. Hypoxia and reoxygenation lead to changes in DNA methylation levels in the gonads. Hypoxia upregulated the expression of dnmt1, dnmt3a, dnmt3b, tet1, and tet2 in females and dnmt3a and dnmt3b in males, while reoxygenation down-regulated the expression of dnmt1, dnmt3a, dnmt3b, tet1, and tet2 in males. Whole genome methylation sequencing showed that the number of differentially methylated regions was highest on chromosome 10 (5192) and lowest on chromosome 24 (275). Differentially methylated genes in females and males, as well as between males and females, were enriched in the oxytocin signaling pathway, fatty acid metabolism pathway, and HIF-1a pathway. In summary, hypoxia and reoxygenation can induce endocrine disorders, affect the expression of HPG axis genes, change the methylation pattern and modification pattern of gonad DNA, and then have potential effects on gonad development.

Single generation epigenetic change in captivity and reinforcement in subsequent generations in a delta smelt (Hypomesus transpacificus) conservation hatchery

A refugial population of the endangered delta smelt (Hypomesus transpacificus) has been maintained at the Fish Conservation and Culture Laboratory (FCCL) at UC Davis since 2008. Despite intense genetic management, fitness differences between wild and cultured fish have been observed at the FCCL. To investigate the molecular underpinnings of hatchery domestication, we used whole-genome bisulfite sequencing to quantify epigenetic differences between wild and hatchery-origin delta smelt. Differentially methylated regions (DMRs) were identified from 104 individuals by comparing the methylation patterns in different generations of hatchery fish (G1, G2, G3) with their wild parents (G0). We discovered a total of 132 significant DMRs (p < .05) between G0 and G1, 132 significant DMRs between G0 and G2, and 201 significant DMRs between G0 and G3. Our results demonstrate substantial differences in methylation patterns emerged between the wild and hatchery-reared fish in the early generations in the hatchery, with a higher proportion of hypermethylated DMRs in hatchery-reared fish. The rearing environment was found to be a stronger predictor of individual clustering based on methylation patterns than family, sex or generation. Our study indicates a reinforcement of the epigenetic status with successive generations in the hatchery environment, as evidenced by an increase in methylation in hypermethylated DMRs and a decrease in methylation in hypomethylated DMRs over time. Lastly, our results demonstrated heterogeneity in inherited methylation pattern in families across generations. These insights highlight the long-term consequences of hatchery practices on the epigenetic landscape, potentially impacting wild fish populations.

The Role of Z Chromosome Localization Gene psmd9 in Spermatogenesis of Cynoglossus semilaevis

Proteasome 26S Subunit, Non-ATPase 9 (psmd9) plays an important role in the balance of protamine and the stability of the nucleolar structure during spermatogenesis. In this study, we cloned the psmd9 of Cynoglossus semilaevis and analyzed its expression pattern. psmd9 was identified on the Z chromosome of C. semilaevis, which is considered an interesting candidate gene for spermatogenesis. qRT-PCR and FISH experiments showed that the psmd9 gene was significantly highly expressed in the testes. It is worth noting that the expression level of psmd9 in male fish testes is significantly higher than that in pseudomales. In order to further explore the role of psmd9 in spermatogenesis, a male testicular cell line was used as the experimental material. The results of the psmd9-RNAi and overexpression experiments showed that psmd9 had a synergistic effect with spermatogenesis-related genes dnd1, cfap69, dnah3 and dnajb13, but had an antagonistic effect with ccne2. Our findings offer a scientific foundation for comprehending the role of psmd9 in the spermatogenesis regulatory network of C. semilaevis.

5 mC modification of steroid hormone biosynthesis-related genes orchestrates feminization of channel catfish induced by high-temperature

To elucidate the mechanism behind channel catfish feminization induced by high temperature, gonad samples were collected from XY pseudo-females and wild-type females and subjected to high-throughput sequencing for Whole-Genome-Bisulfite-Seq (WGBS) and transcriptome sequencing (RNA-Seq). The analysis revealed 50 differentially methylated genes between wild-type females and XY pseudo-females, identified through the analysis of KEGG pathways and GO enrichment in the promoter of the genome and differentially methylated regions (DMRs). Among these genes, multiple differential methylation sites observed within the srd5a2 gene. Repeatability tests confirmed 7 differential methylation sites in the srd5a2 gene in XY pseudo-females compared to normal males, with 1 specific differential methylation site (16608174) distinguishing XY pseudo-females from normal females. Interestingly, the expression of these genes in the transcriptome showed no difference between wild-type females and XY pseudo-females. Our study concluded that methylation of the srd5a2 gene sequence leads to decreased expression, which inhibits testosterone synthesis while promoting the synthesis of 17β-estradiol from testosterone. This underscores the significance of the srd5a2 gene in the sexual differentiation of channel catfish, as indicated by the ipu00140 KEGG pathway analysis.

Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective

Commencing with sperm-egg fusion, the early stages of metazoan development include the cleavage and formation of blastula and gastrula. These early embryonic events play a crucial role in ontogeny and are accompanied by a dramatic remodeling of the gene network, particularly encompassing the maternal-to-zygotic transition. Nonetheless, the gene expression dynamics governing early embryogenesis remain unclear in most metazoan lineages. We conducted transcriptomic profiling on two types of gametes (oocytes and sperms) and early embryos (ranging from the four-cell to the gastrula stage) of an economically valuable flatfish-the Chinese tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae). Comparative transcriptome analysis revealed that large-scale zygotic genome activation (ZGA) occurs in the blastula stage, aligning with previous findings in zebrafish. Through the comparison of the most abundant transcripts identified in each sample and the functional analysis of co-expression modules, we unveiled distinct functional enrichments across different gametes/developmental stages: actin- and immune-related functions in sperms; mitosis, transcription inhibition, and mitochondrial function in oocytes and in pre-ZGA embryos (four- to 1000-cell stage); and organ development in post-ZGA embryos (blastula and gastrula). These results provide insights into the intricate transcriptional regulation of early embryonic development in Cynoglossidae fish and expand our knowledge of developmental constraints in vertebrates.

Molecular functional characterization of the setdb1 and its potential target gene sox5 illuminate the histone modification-mediated orchestration of gonadal development in Chinese tongue sole (Cynoglossus semilaevis)

SET (SuVar3-9, Enhancer of Zeste, Trithorax) domain bifurcated histone lysine methyltransferase 1, setdb1, is the predominant histone lysine methyltransferase catalyzing H3K9me3. Prior studies have illustrated that setdb1 and H3K9me3 critically regulate sex differentiation and gametogenesis. However, the molecular details by which setdb1 is involved in these processes in fish have been poorly reported. Here, we cloned and characterized the setdb1 ORF (open reading frame) sequence from Chinese tongue sole (Cynoglossus semilaevis). The setdb1 ORF sequence was 3,669 bp, encoding a 1,222-amino-acid protein. Phylogenetic analysis showed that setdb1 was structurally conserved. qRT-PCR revealed that setdb1 had a high expression level in the testes at 12 mpf (months post fertilization). Single-cell RNA-seq data at 24 mpf indicated that setdb1 was generally expressed in spermatogenic cells at each stage except for sperm and was centrally expressed in oogonia. H3K9me3 modification was observed in gonads with the immunofluorescence technique. Furthermore, the overexpression experiment suggested that sox5 was a candidate target of setdb1. sox5 was abundantly expressed in male and pseudomale gonads at 24 mpf. Single-cell RNA-seq data showed that sox5 was mainly expressed in spermatogonia and its expression gradually declined with differentiation. Taken together, our findings imply that setdb1 regulates sox5 transcription in gonads, which provides molecular clues into histone modification-mediated orchestration of sex differentiation and gametogenesis.

DNA Methylation Analysis Reveals Potential Mechanism in Takifugu rubripes Against Cryptocaryon irritans Infection

Takifugu rubripes (T. rubripes) is a valuable commercial fish, and Cryptocaryon irritans (C. irritans) has a significant impact on its aquaculture productivity. DNA methylation is one of the earliest discovered ways of gene epigenetic modification and also an important form of modification, as well as an essential type of alteration that regulates gene expression, including immune response. To further explore the anti-infection mechanism of T. rubripes in inhibiting this disease, we determined genome-wide DNA methylation profiles in the gill of T. rubripes using whole-genome bisulfite sequencing (WGBS) and combined with RNA sequence (RNA-seq). A total of 4659 differentially methylated genes (DMGs) in the gene body and 1546 DMGs in the promoter between the infection and control group were identified. And we identified 2501 differentially expressed genes (DEGs), including 1100 upregulated and 1401 downregulated genes. After enrichment analysis, we identified DMGs and DEGs of immune-related pathways including MAPK, Wnt, ErbB, and VEGF signaling pathways, as well as node genes prkcb, myca, tp53, and map2k2a. Based on the RNA-Seq results, we plotted a network graph to demonstrate the relationship between immune pathways and functional related genes, in addition to gene methylation and expression levels. At the same time, we predicted the CpG island and transcription factor of four immune-related key genes prkcb and mapped the gene structure. These unique discoveries could be helpful in the understanding of C. irritans pathogenesis, and the candidate genes screened may serve as optimum methylation-based biomarkers that can be utilized for the correct diagnosis and therapy T. rubripes in the development of the ability to resist C. irritans infection.

Fish reproduction in a warming world: vulnerable points in hormone regulation from sex determination to spawning

Reproduction in fishes is sensitive to temperature. Elevated temperatures and anomalous 'heat waves' associated with climate change have the potential to impact fish reproductive performance and, in some cases, even induce sex reversals. Here we examine how thermal sensitivity in the hormone pathways regulating reproduction provides a framework for understanding impacts of warmer conditions on fish reproduction. Such effects will differ depending on evolved variation in temperature sensitivity of endocrine pathways regulating reproductive processes of sex determination/differentiation, gametogenesis and spawning, as well as how developmental timing of those processes varies with reproductive ecology. For fish populations unable to shift geographical range, persistence under future climates may require changes in temperature responsiveness of the hormone pathways regulating reproductive processes. How thermal sensitivity in those hormone pathways varies among populations and species, how those pathways generate temperature maxima for reproduction, and how rapidly reproductive thermal tolerances can change via adaptation or transgenerational plasticity will shape which fishes are most at risk for impaired reproduction under rising temperatures. This article is part of the theme issue 'Endocrine responses to environmental variation: conceptual approaches and recent developments'.

Does stress make males? An experiment on the role of glucocorticoids in anuran sex reversal

Environmentally sensitive sex determination may help organisms adapt to environmental change but also makes them vulnerable to anthropogenic stressors, with diverse consequences for population dynamics and evolution. The mechanisms translating environmental stimuli to sex are controversial: although several fish experiments supported the mediator role of glucocorticoid hormones, results on some reptiles challenged it. We tested this hypothesis in amphibians by investigating the effect of corticosterone on sex determination in agile frogs (Rana dalmatina). This species is liable to environmental sex reversal whereby genetic females develop into phenotypic males. After exposing tadpoles during sex determination to waterborne corticosterone, the proportion of genetic females with testes or ovotestes increased from 11% to up to 32% at 3 out of 4 concentrations. These differences were not statistically significant except for the group treated with 10 nM corticosterone, and there was no monotonous dose-effect relationship. These findings suggest that corticosterone is unlikely to mediate sex reversal in frogs. Unexpectedly, animals originating from urban habitats had higher sex-reversal and corticosterone-release rates, reduced body mass and development speed, and lower survival compared to individuals collected from woodland habitats. Thus, anthropogenic environments may affect both sex and fitness, and the underlying mechanisms may vary across ectothermic vertebrates.

Three Copies of zbed1 Specific in Chromosome W Are Essential for Female-Biased Sexual Size Dimorphism in Cynoglossus semilaevis

The sex chromosome, especially specific in one sex, generally determines sexual size dimorphism (SSD), a phenomenon with dimorphic sexual difference in the body size. For Cynoglossus semilaevis, a flatfish in China, although the importance of chromosome W and its specific gene zbed1 in female-biased SSD have been suggested, its family members and regulation information are still unknown. At present, three zbed1 copies gene were identified on chromosome W, with no gametologs. Phylogenetic analysis for the ZBED family revealed an existence of ZBED9 in the fish. Nine members were uncovered from C. semilaevis, clustering into three kinds, ZBED1, ZBED4 and ZBEDX, which is less than the eleven kinds of ZBED members in mammals. The predominant expression of zbed1 in the female brain and pituitary tissues was further verified by qPCR. Transcription factor c/ebpα could significantly enhance the transcriptional activity of zbed1 promoter, which is opposite to its effect on the male determinant factor-dmrt1. When zbed1 was interfered with, piwil1, esr2 and wnt7b were up-regulated, while cell-cycle-related genes, including cdk4 and ccng1, were down-regulated. Thus, zbed1 is involved in cell proliferation by regulating esr2, piwil1, cell cycle and the Wnt pathway. Further research on their interactions would be helpful to understand fish SSD.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

The role of DNA methylation reprogramming during sex determination and sex reversal in the Pacific oyster Crassostrea gigas

DNA methylation is instrumental in vertebrate sex reversal. However, the mechanism of DNA methylation regulation regarding sex reversal in invertebrates is unclear. In this study, we used whole genome bisulfite sequencing (WGBS) to map single-base resolution methylation profiles of the Pacific oyster, including female-to-male (FMa-to-FMb) and male-to-female (MFa-to-MFb) sex reversal, as well as sex non-reversed males and females (MMa-to-MMb and FFa-to-FFb). The results showed that global DNA methylation levels increase during female-to-male sex reversals, with a particular increase in the proportion of high methylation levels (mCGs >0.75) and a decrease in the proportion of intermediate methylation levels (0.25 < mCGs <0.75). This increase in DNA methylation was mainly associated with the elevated expression of DNA methylase genes. Genome-wide methylation patterns of females were accurately remodeled to those of males after sex reversal, while the opposite was true for the male-to-female reversal. Those findings directly indicate that alterations in DNA methylation play a significant role in sex reversal in Pacific oysters. Comparative analysis of the DNA methylomes of pre- and post- sex reversal gonadal tissues (FMb-vs-FMa or MFb-vs-MFa) revealed that differentially methylated genes were mainly involved in the biological processes of sex determination or gonadal development. However critical genes such as Dmrt1, Foxl2 and Sox-like, which are involved in the putative sex determination pathway in Pacific oysters, showed almost an absence of methylation modifications, varying greatly from vertebrates. Additionally, comparative analysis of the DNA methylomes of sexual reversal and sex non-reversal (FMa-vs-FFa or MFa-vs-MMa) revealed that heat shock protein genes, such as Hsp68-like and Hsp70B, were important for the occurrence of sex reversal. These findings shed light on the epigenetic mechanisms underlying the maintenance of gonadal plasticity and the reversal of organ architecture in oysters.

Molecular Characterization, Expression Pattern, DNA Methylation and Gene Disruption of Figla in Blotched Snakehead (Channa maculata)

Figla is one of the earliest expressed genes in the oocyte during ovarian development. In this study, Figla was characterized in C. maculata, one of the main aquaculture species in China, and designated as CmFigla. The length of CmFigla cDNA was 1303 bp, encoding 197 amino acids that contained a conserved bHLH domain. CmFigla revealed a female-biased expression patterns in the gonads of adult fish, and CmFigla expression was far higher in ovaries than that in testes at all gonadal development stages, especially at 60~180 days post-fertilization (dpf). Furthermore, a noteworthy inverse relationship was observed between CmFigla expression and the methylation of its promoter in the adult gonads. Gonads at 90 dpf were used for in situ hybridization (ISH), and CmFigla transcripts were mainly concentrated in oogonia and the primary oocytes in ovaries, but undetectable in the testes. These results indicated that Figla would play vital roles in the ovarian development in C. maculata. Additionally, the frame-shift mutations of CmFigla were successfully constructed through the CRISPR/Cas9 system, which established a positive foundation for further investigation on the role of Figla in the ovarian development of C. maculata. Our study provides valuable clues for exploring the regulatory mechanism of Figla in the fish ovarian development and maintenance, which would be useful for the sex control and reproduction of fish in aquaculture.

Identification and Functional Analysis of Cynoglossus semilaevis Z-Linked E3 Ubiquitin Ligase rnf34

The high proportion of males in C. semilaevis hinders their industrial development. The genetic ZW individual can become a pseudomale by sex reversal. And the pseudomale can produce Z-sperm (with epigenetic information to cause sex reversal) while W-sperm is absent, which leads to an even higher male proportion in offspring. Recently, with the development of transcriptomic technologies, research on spermatogenesis in C. semilaevis has been focused on the ubiquitination pathway. In this study, we analyzed the function of the ubiquitin ligase rnf34 gene on the Z chromosome. A qPCR experiment showed that its expression level in the gonad was the highest among different tissues. In the ovary, the expression gradually increased with development from 40 days post-hatching (dph) to 1.5 years post-hatching (yph). In the testis, rnf34 showed increased expression from 40 dph to 6 months post-hatching (mpf) and stabilized up until 1.5 ypf. In situ hybridization showed that the mRNA of rnf34 was mainly distributed in the germ cells of the testis and the ovary. In vivo siRNA-mediated knockdown of the rnf34 gene in male fish affected the expression of a series of genes related to sex differentiation and spermatogenesis. These results provide genetic data on the molecular mechanisms of gonadal development and spermatogenesis in C. semilaevis.

Effect of high temperatures on sex ratio and differential expression analysis (RNA-seq) of sex-determining genes in Oreochromis niloticus from different river basins in Benin

The high temperature sex reversal process leading to functional phenotypic masculinization during development has been widely described in Nile tilapia (Oreochromis n iloticus) under laboratory or aquaculture conditions and in the wild. In this study, we selected five wild populations of O. niloticus from different river basins in Benin and produced twenty full-sib families of mixed-sex (XY and XX) by natural reproduction. Progenies were exposed to room temperature or high (36.5°C) temperatures between 10 and 30 days post-fertilization (dpf). In control groups, we observed sex ratios from 40% to 60% males as expected, except for 3 families from the Gobé region which showed a bias towards males. High temperature treatment significantly increased male rates in each family up to 88%. Transcriptome analysis was performed by RNA-sequencing (RNA-seq) on brains and gonads from control and treated batches of six families at 15 dpf and 40 dpf. Analysis of differentially expressed genes, differentially spliced genes, and correlations with sex reversal was performed. In 40 dpf gonads, genes involved in sex determination such as dmrt1, cyp11c1, amh, cyp19a1b, ara, and dax1 were upregulated. In 15 dpf brains, a negative correlation was found between the expression of cyp19a1b and the reversal rate, while at 40 dpf a negative correlation was found between the expression of foxl2, cyp11c1, and sf1 and positive correlation was found between dmrt1 expression and reversal rate. Ontology analysis of the genes affected by high temperatures revealed that male sex differentiation processes, primary male sexual characteristics, autophagy, and cilium organization were affected. Based on these results, we conclude that sex reversal by high temperature treatment leads to similar modifications of the transcriptomes in the gonads and brains in offspring of different natural populations of Nile tilapia, which thus may activate a common cascade of reactions inducing sex reversal in progenies.

Gamma irradiation-induced offspring masculinization is associated with epigenetic changes in female zebrafish

Sex ratio variation is a key topic in ecology, because of its direct effects on population dynamics and thus, on animal conservation strategies. Among factors affecting sex ratio, types of sex determination systems have a central role, since some species could have a sex determined by genetic factors, environmental factors or a mix of those two. Yet, most studies on the factors affecting sex determination have focused on temperature or endocrine-disrupting chemicals (EDCs), and much less is known regarding other factors. Exposure to gamma irradiation was found to trigger offspring masculinization in zebrafish. Here we aimed at deciphering the potential mechanisms involved, by focusing on stress (i.e. cortisol) and epigenetic regulation of key genes involved in sex differentiation in fish. Cortisol levels in exposed and control (F0) zebrafish females' gonads were similar. However, irradiation increased the DNA methylation level of foxl2a and cyp19a1a in females of the F0 and F1 generation, respectively, while no effects were detected in testis. Overall, our results suggest that parental exposure could alter offspring sex ratio, at least in part by inducing methylation changes in ovaries.

Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads

Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.

Alternative splicing of histone demethylase Kdm6bb mediates temperature-induced sex reversal in the Nile tilapia

Sex determination in many fish species is remarkably plastic and temperature sensitive. Nile tilapia display a genetic sex-determination system (XX/XY). However, high-temperature treatment during critical thermosensitive periods can induce XX females into XXm pseudo-males, and this phenomenon is termed temperature-induced sex reversal (TISR). To investigate the molecular mechanism of TISR in Nile tilapia, we performed Iso-seq analysis and found a dramatic effect of high temperature on gene alternative splicing (AS). Kdm6bb histone demethylase showed a novel AS at intron 5 that generates Kdm6bb_tv1 transcripts without intron 5 and Kdm6bb_tv2 with intron 5. Kdm6bb_tv1 encodes a full-length protein while Kdm6bb_tv2 encodes a truncated protein. Expression analysis revealed that intron 5 splicing of Kdm6bb is male and gonad biased at larval stage, and only gonad biased at adult stage. High-temperature treatment induced intron 5 splicing in the gonads of XX and XY fish, resulting in increased Kdm6bb_tv1 expression. To directly test the role of Kdm6bb_tv1 in Nile tilapia TISR, we knocked out expression of Kdm6bb_tv1. However, Kdm6bb_tv1-/- homozygous mutants showed embryonic lethality. Overexpression of Kdm6bb_tv1, but not Kdm6bb_tv2, induced sex reversal of XX females into pseudo-males. Overexpression of Kdm6bb_tv1, as with high-temperature treatment, modified the promotor region of Gsdf and Dmrt1 by demethylating the trimethylated lysine 27 of histone 3 (H3K27me3), thereby increasing expression. Collectively, these studies demonstrate that AS of Kdm6bb intron 5 increases the expression of Kdm6bb_tv1, which acts as a direct link between high temperature and activation of Gsdf and Dmrt1 expression, leading to male sex determination.

Sex-specific DNA methylation and transcription of zbtb38 and effects of gene-environment interactions on its natural antisense transcript in zebrafish

There is increasing evidence for the involvement of epigenetics in sex determination, maintenance, and plasticity, from plants to humans. In our previous work, we reported a transgenerational feminization of a zebrafish population for which the first generation was exposed to cadmium, a metal with endocrine disrupting effects. In this study, starting from the previously performed whole methylome analysis, we focused on the zbtb38 gene and hypothesized that it could be involved in sex differentiation and Cd-induced offspring feminization. We observed sex-specific patterns of both DNA methylation and RNA transcription levels of zbtb38. We also discovered that the non-coding exon 3 of zbtb38 encodes for a natural antisense transcript (NAT). The activity of this NAT was found to be influenced by both genetic and environmental factors. Furthermore, increasing transcription levels of this NAT in parental gametes was highly correlated with offspring sex ratios. Since zbtb38 itself encodes for a transcription factor that binds methylated DNA, our results support a non-negligible role of zbtb38 not only in orchestrating the sex-specific transcriptome (i.e., sex differentiation) but also, via its NAT, offspring sex ratios.

The fusion gene hsf5-rnf43 in Nile tilapia: A potential regulator in the maintenance of testis function and sexual differentiation

We identified that the genes heat shock transcription factor 5 (hsf5) and ring finger protein 43 (rnf43) happened fusion in Nile tilapia (Oreochromis niloticus), called hsf5-rnf43, and provided the characteristic and functional analysis of hsf5-rnf43 gene in fish for the first time. Analysis of spatiotemporal expression showed that hsf5-rnf43 was specifically expressed in the testis and located in primary spermatocytes of adult Nile tilapia and gradually increased during testis development from 5 to 180 days after hatching. We also found DNA methylation regulated sex-biased expression of hsf5-rnf43 in the early development of Nile tilapia, and was affected by high temperature during the thermosensitive period of Nile tilapia sex differentiation. Therefore, we first reported that the fusion gene hsf5-rnf43 was sex-biased expressed in the testis regulated by DNA methylation and affected by high temperature, which may be involved in the maintenance of testis function and sex differentiation of Nile tilapia.

Molecular biological, physiological, cytological, and epigenetic mechanisms of environmental sex differentiation in teleosts: A systematic review

Human activities have been exerting widespread stress and environmental risks in aquatic ecosystems. Environmental stress, including temperature rise, acidification, hypoxia, light pollution, and crowding, had a considerable negative impact on the life histology of aquatic animals, especially on sex differentiation (SDi) and the resulting sex ratios. Understanding how the sex of fish responds to stressful environments is of great importance for understanding the origin and maintenance of sex, the dynamics of the natural population in the changing world, and the precise application of sex control in aquaculture. This review conducted an exhaustive search of the available literature on the influence of environmental stress (ES) on SDi. Evidence has shown that all types of ES can affect SDi and universally result in an increase in males or masculinization, which has been reported in 100 fish species and 121 cases. Then, this comprehensive review aimed to summarize the molecular biology, physiology, cytology, and epigenetic mechanisms through which ES contributes to male development or masculinization. The relationship between ES and fish SDi from multiple aspects was analyzed, and it was found that environmental sex differentiation (ESDi) is the result of the combined effects of genetic and epigenetic factors, self-physiological regulation, and response to environmental signals, which involves a sophisticated network of various hormones and numerous genes at multiple levels and multiple gradations in bipotential gonads. In both normal male differentiation and ES-induced masculinization, the stress pathway and epigenetic regulation play important roles; however, how they co-regulate SDi is unclear. Evidence suggests that the universal emergence or increase in males in aquatic animals is an adaptation to moderate ES. ES-induced sex reversal should be fully investigated in more fish species and extensively in the wild. The potential aquaculture applications and difficulties associated with ESDi have also been addressed. Finally, the knowledge gaps in the ESDi are presented, which will guide the priorities of future research.

Exploring the Effects of Rearing Densities on Epigenetic Modifications in the Zebrafish Gonads

Rearing density directly impacts fish welfare, which, in turn, affects productivity in aquaculture. Previous studies have indicated that high-density rearing during sexual development in fish can induce stress, resulting in a tendency towards male-biased sex ratios in the populations. In recent years, research has defined the relevance of the interactions between the environment and epigenetics playing a key role in the final phenotype. However, the underlying epigenetic mechanisms of individuals exposed to confinement remain elucidated. By using zebrafish (Danio rerio), the DNA methylation promotor region and the gene expression patterns of six genes, namely dnmt1, cyp19a1a, dmrt1, cyp11c1, hsd17b1, and hsd11b2, involved in the DNA maintenance methylation, reproduction, and stress were assessed. Zebrafish larvae were subjected to two high-density conditions (9 and 66 fish/L) during two periods of overlapping sex differentiation of this species (7 to 18 and 18 to 45 days post-fertilization, dpf). Results showed a significant masculinization in the populations of fish subjected to high densities from 18 to 45 dpf. In adulthood, the dnmt1 gene was differentially hypomethylated in ovaries and its expression was significantly downregulated in the testes of fish exposed to high-density. Further, the cyp19a1a gene showed downregulation of gene expression in the ovaries of fish subjected to elevated density, as previously observed in other studies. We proposed dnmt1 as a potential testicular epimarker and the expression of ovarian cyp19a1a as a potential biomarker for predicting stress originated from high densities during the early stages of development. These findings highlight the importance of rearing densities by long-lasting effects in adulthood conveying cautions for stocking protocols in fish hatcheries.

Exosome delivery to the testes for dmrt1 suppression: A powerful tool for sex-determining gene studies

Exosomes are endosome-derived extracellular vesicles about 100 nm in diameter. They are emerging as promising delivery platforms due to their advantages in biocompatibility and engineerability. However, research into and applications for engineered exosomes are still limited to a few areas of medicine in mammals. Here, we expanded the scope of their applications to sex-determining gene studies in early vertebrates. An integrated strategy for constructing the exosome-based delivery system was developed for efficient regulation of dmrt1, which is one of the most widely used sex-determining genes in metazoans. By combining classical methods in molecular biology and the latest technology in bioinformatics, isomiR-124a was identified as a dmrt1 inhibitor and was loaded into exosomes and a testis-targeting peptide was used to modify exosomal surface for efficient delivery. Results showed that isomiR-124a was efficiently delivered to the testes by engineered exosomes and revealed that dmrt1 played important roles in maintaining the regular structure and function of testis in juvenile fish. This is the first de novo development of an exosome-based delivery system applied in the study of sex-determining gene, which indicates an attractive prospect for the future applications of engineered exosomes in exploring more extensive biological conundrums.

The admixed brushtail possum genome reveals invasion history in New Zealand and novel imprinted genes

Combining genome assembly with population and functional genomics can provide valuable insights to development and evolution, as well as tools for species management. Here, we present a chromosome-level genome assembly of the common brushtail possum (Trichosurus vulpecula), a model marsupial threatened in parts of their native range in Australia, but also a major introduced pest in New Zealand. Functional genomics reveals post-natal activation of chemosensory and metabolic genes, reflecting unique adaptations to altricial birth and delayed weaning, a hallmark of marsupial development. Nuclear and mitochondrial analyses trace New Zealand possums to distinct Australian subspecies, which have subsequently hybridised. This admixture allowed phasing of parental alleles genome-wide, ultimately revealing at least four genes with imprinted, parent-specific expression not yet detected in other species (MLH1, EPM2AIP1, UBP1 and GPX7). We find that reprogramming of possum germline imprints, and the wider epigenome, is similar to eutherian mammals except onset occurs after birth. Together, this work is useful for genetic-based control and conservation of possums, and contributes to understanding of the evolution of novel mammalian epigenetic traits.

Estrodiol-17β and aromatase inhibitor treatment induced alternations of genome-wide DNA methylation pattern in Takifugu rubripes gonads

Estradiol-17β (E2) and aromatase inhibitor (AI) exposure can change the phenotypic sex of fish gonads. To investigated whether alterations in DNA methylation is involved in this process, the level of genome-wide DNA methylation in Takifugu rubripes gonads was quantitatively analyzed during the E2-induced feminization and AI-induced masculinization processes in this study. The methylation levels of the total cytosine (C) in control-XX(C-XX), control-XY (C-XY), E2-treated-XY (E-XY) and AI-treated-XX (AI-XX) were 9.11%, 9.19%, 8.63% and 9.23%, respectively. In the C-XX vs C-XY comparison, 4,196 differentially methylated regions (DMRs) overlapped with the gene body of 2,497 genes and 608 DMRs overlapped with the promoter of 575 genes. In the E-XY vs C-XY comparison, 6,539 DMRs overlapped with the gene body of 3,416 genes and 856 DMRs overlapped with the promoter of 776 genes. In the AI-XX vs C-XX comparison, 2,843 DMRs overlapped with the gene body of 1,831 genes and 461 DMRs overlapped with the promoter of 421 genes. Gonadal genomic methylation mainly occurred at CG sites and the genes that overlapped with DMRs on CG context were most enriched in the signaling pathways related to gonad differentiation, such as the Wnt, TGF-β, MAPK, CAM and GnRH pathways. The DNA methylation levels of steroid synthesis genes and estrogen receptor genes promoter or gene body were negative correlated with their expression. After bisulfite sequencing verification, the DNA methylation level of the amhr2 promoter in XY was increased after E2 treatment, which consistent with the data from the genome-wide DNA methylation sequencing. In C-XY group, the expression of amhr2 was significantly higher than that in E-XY (p < 0.05). Additionally, dnmt1, which is responsible for methylation maintenance, expressed at similar level in four groups (p > 0.05). dnmt3, tet2, and setd1b, which were responsible for methylation modification, expressed at significantly higher levels in E-XY compared to the C-XY (p < 0.05). Dnmt3 and tet2 were expressed at significantly higher levels in AI-XX than that in C-XX (p < 0.05). These results indicated that E2 and AI treatment lead to the aberrant genome-wide DNA methylation level and expression level of dnmt3, tet2, and setd1b in T. rubripes gonad.

Phosphoproteomics Reveal New Candidates in Abnormal Spermatogenesis of Pseudomales in Cynoglossus semilaevis

Phosphorylation is a post-translational modification that contributes to versatile protein functions in spermatogenesis, and the variations they generate usually results in abnormal spermatogenesis or sperm dysfunction. The sex-reversal phenomenon exists in Chinese tongue sole under certain conditions such that individuals with a ZW genotype can acquire a male phenotype and are thus called pseudomales. Pseudomale tongue sole can reach sexual maturity but produce only Z-type sperm, and the Z sperm carries paternal epigenetic information. Whether phosphorylation plays a role in the sperm abnormality of pseudomales is unknown. In this study, a phosphoproteomic analysis was performed to compare protein phosphorylation profiles between pseudomale and male testes. Altogether, we identified 14,253 phosphopeptides matching with 4843 proteins, with 1329 differentially phosphorylated peptides corresponding to 1045 differentially phosphorylated proteins (DPPs). Phosphorylation at 781 sites was upregulated and at 548 sites was downregulated. Four motifs were identified among differentially phosphorylated peptides, which were "SP", "SD", "RxxS", and "TP". Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that the cell cycle and DNA/RNA processing were significantly enriched with the genes encoding DPPs. To analyze DPP function in depth, a protein-protein interaction network was constructed, and Ran-binding protein 2 was found to play a central role in spermatogenesis by regulating several processes such as the cell cycle, eukaryotic translation, ubiquitination, and minichromosome maintenance. In kinase-associated network analyses, two "mitogen-activated protein kinase (Mapk)-centered" clusters were identified that may account for abnormal spermatogenesis in pseudomales. One cluster was centered on Mapk6, which predominantly regulated the cell cycle by interacting with several cyclin-dependent kinases, and the other was centered on the "testis-expressed kinase 1-like (Tesk1l)/Pim1l-Mapk4l- testis-expressed 14 (Tex14)" kinase cascade, which might contribute to spermatogenesis by regulating β-catenin. Taken together, these data suggested the new candidates involved in pseudomale sperm abnormalities and provided clues to discover the phosphorylated regulatory mechanism underlying tongue sole spermatogenesis.

Starvation-induced changes in sex ratio involve alterations in sex-related gene expression and methylation in Pacific oyster Crassostrea gigas

Aquatic animals are subject to varying degrees of starvation stress in their natural habitats due to food limitations. Consequently, starvation is a crucial environmental factor for sex determination in many species; however, limited research has been conducted on the effects of starvation on sex determination in shellfish. Here, four full sibling families of Pacific oysters were established and subjected to starvation stress. The results demonstrated that starvation caused the sex ratio (female to male) to change from 1:0.78 to 1:1.44 and resulted in a delay in gonadal development. Further studies revealed that the expression levels of DNA methylation-related genes Dnmt1 (DNA methyltransferase 1), Dnmt3a/b (DNA methyltransferase 3a/b) and Tet3 (tet methylcytosine dioxygenase 3) decreased under starvation stress. Conversely, the upregulation of Dmrt1 (doublesex and mab-3 related transcription factor 1), a gene typically associated with males, in females suggests that the increased proportion of males may be linked to starvation-induced high expression of this particular gene. In addition, the gene Dgkd (diacylglycerol kinase delta), which is involved in the regulation of second messenger protein kinase C, was differentially methylated between males and females, with the methylation level of this gene gradually increasing with male development, while the methylation level of this gene decreased under starvation stress. Correlation analysis of Dgkd methylation levels with expression levels showed a negative correlation between DNA methylation and gene expression. Finally dual fluorescence reporter experiments confirmed that DNA methylation suppressed Dgkd expression in vitro. In summary, the results suggest that starvation may regulate Dgkd gene expression through DNA methylation and thus affect Dmrt1 expression, thereby causing sex reversal in the Pacific oyster. The outcomes resolved how environmental factors are involved in sex reversal from an epigenetic perspective and provided a theoretical basis for sex control breeding in the Pacific oyster.

Study of DNA methylation of hsd17β, er and reproductive endocrine disrupting effects in female Chlamys farreri under benzo[a]pyrene stress

Benzo[a]pyrene (B[a]P) is one kind of persistent organic pollutants (POPs) in the marine environment which has multiple toxic effects. However, epigenetic studies correlated with reproductive endocrine disruption in invertebrates have not been explored. In our study, Chlamys farreri in the mature stage were exposed to B[a]P (0, 0.4, 2 and 10 μg/L) for 5 and 10 d to explore the effects on reproductive endocrine and DNA methylation. The results proved that B[a]P stress significantly restrained the growth of mature oocytes, reduced the content of sex hormones, and affected the expression of genes related to ovarian development. Histological observation showed that the ovarian microstructure was damaged. The detection of SAM/SAH, dnmts, GNMT in the ovary showed that the level of global DNA methylation fluctuated. Significant hypermethylation of the hsd17β promoter region in the ovary was associated with a significant downregulation of its gene expression. In summary, our results suggested that exposure to B[a]P might affect DNA methylation to regulate key reproductive genes, interfere with the synthesis of sex hormones, and inhibit ovarian development. These findings provide a basis for a better understanding of how epigenetic mechanisms are involved in the response of marine invertebrates to POPs stress, opening up new avenues for incorporating environmental epigenetic approaches into marine invertebrate management and conservation plans.

The hunger strikes back: an epigenetic memory for autophagy

Historical and demographical human cohorts of populations exposed to famine, as well as animal studies, revealed that exposure to food deprivation is associated to lasting health-related effects for the exposed individuals, as well as transgenerational effects in their offspring that affect their diseases' risk and overall longevity. Autophagy, an evolutionary conserved catabolic process, serves as cellular response to cope with nutrient starvation, allowing the mobilization of an internal source of stored nutrients and the production of energy. We review the evidence obtained in multiple model organisms that support the idea that autophagy induction, including through dietary regimes based on reduced food intake, is in fact associated to improved health span and extended lifespan. Thereafter, we expose autophagy-induced chromatin remodeling, such as DNA methylation and histone posttranslational modifications that are known heritable epigenetic marks, as a plausible mechanism for transgenerational epigenetic inheritance of hunger.

Prenatal exposure to acrylamide differently affected the sex ratio, aromatase and apoptosis in female adult offspring of two subsequent generations

In the present study, we investigated the effect of acrylamide (ACR) exposure during pregnancy on the ovary of female adult offspring of two subsequent generations. Sixty-day-old Wistar albino female rats were given different doses of ACR (2.5 and 10 mg/kg/day) from day 6 of pregnancy until giving birth. Females from the first generation (AF1) were fed ad libitum, and thereafter, a subgroup was euthanized at 8 weeks of age and ovary samples were obtained. The remaining females were maintained until they reached sexual maturity (50 days old) and then treated in the same way as the previous generation to obtain the second generation of females (AF2). The histopathological examination indicated a high frequency of corpora lutea along with an increased number of antral follicles that reached the selectable stage mainly at a dose of 2.5 mg/kg/day. Interestingly, ACR exposure significantly increased the mRNA levels of CYP19 gene and its corresponding CYP19 protein expression in AF1 females. The TUNEL assay showed a significantly high rate of apoptosis in stromal cells except for dose of 2.5 mg/kg/day. However, in AF2 females, ACR exposure significantly increased the number of degenerating follicles and cysts while the number of growing follicles was reduced. Moreover, in both ACR-treated groups, estradiol-producing enzyme CYP19A gene and its corresponding protein were significantly reduced, and an excessive apoptosis was produced. We concluded that the ovarian condition of AF1 females had considerable similarity to the typical early perimenopausal stage, whereas that of AF2 females was similar to the late perimenopausal stage in women.

The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis)

Fish sex determination can be affected by environmental temperature. This process relies on temperature-sensitive proteins such as heat shock proteins (HSPs). Our previous work found that heat shock cognate proteins (HSCs) may participate in high-temperature associated sex reversal of Chinese tongue sole (Cynoglossus semilaevis). However, the role of hsc genes in responding to high temperature and affecting sex determination/differentiation remains unclear. Here, by using C. semilaevis as model, we identified hsc70 and hsc70-like. hsc70 was abundant in the gonads with a testicular-higher expression at all gonadal development stages except for 6 months post fertilization (mpf). Intriguingly, hsc70-like showed higher expression in testes from 6 mpf on. Both long-term heat treatment during the temperature-sensitive sex-determining period and short-term heat stress at the end of this period caused different expression of hsc70/hsc70-like between sexes. The dual-luciferase assay results also suggested that these genes can respond to high temperature rapidly in vitro. Heat treatment of C. semilaevis testis cells overexpressed with hsc70/hsc70-like could affect the expression of sex-related genes sox9a and cyp19a1a. Our results indicated that hsc70 and hsc70-like were key regulators linking external high-temperature signals with sex differentiation in vivo and provide a new idea for understanding the mechanism by which high temperature affects sex determination/differentiation in teleosts.

Development of epigenetic biomarkers for the identification of sex and thermal stress in fish using DNA methylation analysis and machine learning procedures

The sex ratio is a key ecological demographic parameter crucial for population viability. However, the epigenetic mechanisms operating during gonadal development regulating gene expression and the sex ratio remain poorly understood. Moreover, there is interest in the development of epigenetic markers associated with a particular phenotype or as sentinels of environmental effects. Here, we profiled DNA methylation and gene expression of 10 key genes related to sex development and stress, including steroidogenic enzymes, and growth and transcription factors. We provide novel information on the sex-related differences and on the influence of elevated temperature on these genes in zebrafish, a species with mixed genetic and environmental influences on sex ratios. We identified both positive (e.g., amh, cyp11c and hsd11b2) and negative (e.g., cyp11a1 and dmrt1) correlations in unexposed males, and negative correlation (amh) in exposed females between DNA methylation and gene expression levels. Further, we combined DNA methylation analysis with machine learning procedures and found a series of informative CpGs capable not only of correctly identifying sex (based on cyp19a1a DNA methylation levels) but also of identifying whether males and females had been exposed to abnormally elevated temperature when young (based on amh and foxl2a DNA methylation levels, respectively). This was achieved in the absence of conspicuous morphological alterations of the gonads. These DNA methylation-based epigenetic biomarkers represent molecular resources that can correctly recapitulate past thermal history and pave the way for similar findings in other species to assess potential ecological effects of environmental disturbances in the context of climate change.

Three dimensions of thermolabile sex determination

The molecular mechanism of temperature-dependent sex determination (TSD) is a long-standing mystery. How is the thermal signal sensed, captured and transduced to regulate key sex genes? Although there is compelling evidence for pathways via which cells capture the temperature signal, there is no known mechanism by which cells transduce those thermal signals to affect gene expression. Here we propose a novel hypothesis we call 3D-TSD (the three dimensions of thermolabile sex determination). We postulate that the genome has capacity to remodel in response to temperature by changing 3D chromatin conformation, perhaps via temperature-sensitive transcriptional condensates. This could rewire enhancer-promoter interactions to alter the expression of key sex-determining genes. This hypothesis can accommodate monogenic or multigenic thermolabile sex-determining systems, and could be combined with upstream thermal sensing and transduction to the epigenome to commit gonadal fate.

Whole-genome resequencing and bisulfite sequencing provide new insights into the feeding habit domestication in mandarin fish (Siniperca chuatsi)

Mandarin fish (Siniperca chuatsi) is one of the most economically important fish in China. However, it has the peculiar feeding habit that it feeds solely on live prey fish since first-feeding, while refuses dead prey fish or artificial diets. After the specific training procedure, partial individuals could accept dead prey fish and artificial diets. The genetic basis of individual difference in artificial diet feeding habit is still unknown. In the present study, the resequencing was performed between 10 individuals which could be domesticated to accept artificial diets and 10 individuals which could not. Through the selective sweep analysis based on heterozygosity (Hp) and population differentiation coefficient (Fst), 57 candidate windows were identified as the putative selected regions for feeding habit domestication of mandarin fish, involved in 149 genes. These genes were related to memory, vision and olfaction function, which could be potential targets of molecular marker assistant breeding of artificial diet feeding trait. Beside of the DNA sequence, we also explored the potential role of DNA methylation in feeding habit domestication in mandarin fish. Whole-genome bisulfite sequencing was performed between the individuals which could be domesticated to accept artificial diets and those could not. 5,976 differentially methylated regions were identified, referring to 3,522 genes, such as the genes involved in cAMP signaling pathway. The DNA methylation changes of these genes might contribute to the adaption of artificial diets in mandarin fish. In conclusion, the putative selected regions and the differentially methylated regions were identified in the whole genome, providing new insights into the feeding habit domestication from live prey fish to artificial diets in mandarin fish. And the involved genes were identified as the candidate genes for molecular breeding of artificial diet utilization in mandarin fish.

Genomic imprinting-like monoallelic paternal expression determines sex of channel catfish

The X and Y chromosomes of channel catfish have the same gene contents. Here, we report allelic hypermethylation of the X chromosome within the sex determination region (SDR). Accordingly, the X-borne hydin-1 gene was silenced, whereas the Y-borne hydin-1 gene was expressed, making monoallelic expression of hydin-1 responsible for sex determination, much like genomic imprinting. Treatment with a methylation inhibitor, 5-aza-dC, erased the epigenetic marks within the SDR and caused sex reversal of genetic females into phenotypic males. After the treatment, hydin-1 and six other genes related to cell cycle control and proliferative growth were up-regulated, while three genes related to female sex differentiation were down-regulated in genetic females, providing additional support for epigenetic sex determination in catfish. This mechanism of sex determination provides insights into the plasticity of genetic sex determination in lower vertebrates and its connection with temperature sex determination where DNA methylation is broadly involved.

Transcriptomic Analysis Revealed Candidate Genes Involved in Pseudomale Sperm Abnormalities in Chinese Tongue Sole (Cynoglossus semilaevis)

Chinese tongue sole (Cynoglossus semilaevis) has a ZZ/ZW sex determination system, but the genotypic female (ZW) can be sex-reversed into phenotypic males, namely, pseudomales. Pseudomale fish can produce only Z-type sperm but not W sperm. However, the molecular mechanism is unclear. To screen the key genes involved in pseudomale sperm abnormalities, we analysed the transcriptomic profiles of pseudomale and male sperm. In comparison to male sperm, 592 differentially expressed genes (DEGs) were identified in pseudomale sperm, including 499 upregulated and 93 downregulated genes. KEGG analysis indicated that the FoxO signalling pathway, especially the foxo3a and foxo6-like genes, may play an important role in spermatogenesis. The DEGs were mainly distributed on sex chromosomes, with 158 downregulated genes on the Z chromosome and 41 upregulated genes on the W chromosome. A specific area (14-15 M) on the Z chromosome was identified, which enriched eight DEGs inside the ~1 M region. In addition, there were five gene alleles on the sex chromosomes, which showed the opposite transcription pattern (upregulated for the W allele, downregulated for the Z allele). This study has provided valuable data for screening candidate genes involved in the pseudomale sperm abnormality.

Single-cell-resolution transcriptome map revealed novel genes involved in testicular germ cell progression and somatic cells specification in Chinese tongue sole with sex reversal

Female-to-male sex reversals (pseudomales) are common in lower vertebrates and have been found in natural populations, which is a concern under rapid changes in environmental conditions. Pseudomales can exhibit altered spermatogenesis. However, the regulatory mechanisms underlying pseudomale spermatogenesis remain unclear. Here, we characterized spermatogenesis in Chinese tongue sole (Cynoglossus semilaevis), a species with genetic and environmental sex determination, based on a high-resolution single-cell RNA-seq atlas of cells derived from the testes of genotypic males and pseudomales. We identified five germ cell types and six somatic cell types and obtained a single-cell atlas of dynamic changes in gene expression during spermatogenesis in Chinese tongue sole, including alterations in pseudomales. We detected decreased levels of Ca²⁺ signaling pathway-related genes in spermatogonia, insufficient meiotic initiation in spermatocytes, and a malfunction of somatic niche cells in pseudomales. However, a cluster of CaSR genes and MAPK signaling factors were upregulated in undifferentiated spermatogonia of pseudomales. Additionally, we revealed that Z chromosome-specific genes, such as piwil2, dhx37, and ehmt1, were important for spermatogenesis. These results improve our understanding of reproduction after female-to-male sex-reversal and provide new insights into the adaptability of reproductive strategies in lower vertebrates.

Chromosome-level Assembly, Dosage Compensation and Sex-biased Gene Expression in the Small Brown Planthopper, Laodelphax striatellus

In insects, sex chromosome differentiation often results in unequal gene dosages between sexes. Dosage compensation mechanisms evolve to balance gene expression, but the degree and mechanism of regulation often vary by insect species. In hemipteran species, the small brown planthopper (SBPH), Laodelphax striatellus, is an injurious crop pest, with a sex chromosome type XX in females and XO in males. This species offers the opportunity to study dosage compensation and sex-biased gene expression. In this study, we generated a chromosome-level genome of SBPH using Oxford Nanopore Technologies and high-throughput chromatin conformation capture (Hi-C) technology. We also sequenced RNA-seq data from 16 tissue samples to annotate the genome and analyze gene dosage compensation. We finally obtained a 510.2 megabases (Mb) genome with 99.12% of the scaffolds anchored on 15 chromosomes (14 autosomes and 1 X chromosome) and annotated 16,160 protein-coding genes based on full-length cDNA sequencing data. Furthermore, we found complete dosage compensation in all L. striatellus somatic tissues, but lack of dosage compensation in gonad tissue testis. We also found that female-biased genes were significantly enriched on the X chromosome in all tissues, whereas male-biased genes in gonad tissues were enriched on autosomes. This study not only provides a high-quality genome assembly but also lays a foundation for a better understanding of the sexual regulatory network in hemipteran insects.

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

CircRNA Identification and CircRNA-miRNA-mRNA Network in Cynoglossus semilaevis Sexual Size Dimorphism

Sexual size dimorphism (SSD), which is the sexual differences in body size, has been widely reported in various species including fishes. For Chinese tongue sole (Cynoglossus semilaevis), a flatfish exhibiting typically female-biased SSD, little is known for its epigenetic regulation mechanism, especially the role of circRNAs. Here, we identified the differently expressed abundances of circRNAs in females, males, and pseudo-males to explore the potential functions of circRNAs in Chinese tongue sole SSD. In total, 14,745 novel circRNAs were screened, among which 1461 DE circRNAs were identified from the brain, gonad, liver, and muscle in female, male, and pseudo-male individuals. The ceRNA network was subsequently constructed, including 10 circRNAs, 26 mRNAs, and 11 miRNAs. These DE mRNAs were mainly related to the mRNA surveillance pathway, metabolic pathways, and cellular senescence. Importantly, the ceRNA network has revealed that several circRNAs such as novel_circ_004374 and novel_circ_014597 may regulate homeodomain interacting protein kinase 2 (hipk2) expression by sponging miR-130-x. It is also worth exploring whether or how novel_circ_008696 regulates SET Domain Containing 2, histone lysine methyltransferase (setd2), which in turn affects the epigenetic patterns of different sexual individuals. The present study not only enriches the knowledge on the potential roles of circRNA in the physiological process, but also provides new clues for the explanation of fish SSD. In future studies, the precise function and involvement of circRNAs in female-biased SSD will require more efforts.

Identification and Expression Pattern of cyp26b1 Gene in Gonad of the Chinese Tongue Sole (Cynoglossus semilaevis)

Simple Summary

In fish, it is obvious that the asynchronous development of the gonads and sexual dimorphism limit the development of aquaculture, so the research into sex-differentiation and gonadal growth is very important. Due to the sexual reversal phenomenon (genetic females becoming phenotypic males), the Chinese tongue sole (Cynoglossus semilaevis) is a great model for investigating sex-differentiation. Herein, we report one gene involved in sex-differentiation and gonadal growth of the Chinese Tongue Sole. The gene cyp26b1 (cytochrome P450 family 26 subfamily b member 1) is a metabolizing Retinoic Acid (RA) enzyme. Since it regulates RA to control sex determination and differentiation, cyp26b1 is considered a critical part of mammals’ ovary-antagonizing and testis-determining downstream passageway of Sry (sex-determining region Y) and Sox9 (sry-box transcription factor 9). In fish, the related research is reported only on the Japanese flounder (Paralichthys olivaceus) and zebrafish (Danio rerio). In the current investigation, the identification and expression pattern of the cyp26b1 gene in the Chinese tongue sole suggested that cyp26b1 might impact sex-differentiation and gonadal development.

Abstract

As an RA-metabolizing enzyme, cyp26b1 has a substantial impact on RA-signaling pathways. The cyp26b1 gene from the Chinese tongue sole was cloned and identified in this investigation. The cyp26b1 ORF was 1536 bp in length and encoded a 512 amino acid protein. A quantitative real-time PCR (qPCR) indicated that the cyp26b1 expression is no significant sexual dimorphism in the gonads at the 80 days post-hatching (dph) stages. After 4 months post-hatching (mph), the expression of cyp26b1 showed sexual dimorphism and lower level of expression in the ovaries than in the testes. An in situ hybridization demonstrated that cyp26b1 mRNA was primarily located in the testis. Interestingly, the cyp26b1 mRNA probe was also detected in the ovaries. These results suggested that cyp26b1 participates in the sex-differentiation and gonadal development of the Chinese tongue sole.

Fishing for fish environmental DNA: Ecological applications, methodological considerations, surveying designs, and ways forward

Vast global declines of freshwater and marine fish diversity and population abundance pose serious threats to both ecosystem sustainability and human livelihoods. Environmental DNA (eDNA)-based biomonitoring provides robust, efficient, and cost-effective assessment of species occurrences and population trends in diverse aquatic environments. Thus, it holds great potential for improving conventional surveillance frameworks to facilitate fish conservation and fisheries management. However, the many technical considerations and rapid developments underway in the eDNA arena can overwhelm researchers and practitioners new to the field. Here, we systematically analysed 416 fish eDNA studies to summarize research trends in terms of investigated targets, research aims, and study systems, and reviewed the applications, rationales, methodological considerations, and limitations of eDNA methods with an emphasis on fish and fisheries research. We highlighted how eDNA technology may advance our knowledge of fish behaviour, species distributions, population genetics, community structures, and ecological interactions. We also synthesized the current knowledge of several important methodological concerns, including the qualitative and quantitative power eDNA has to recover fish biodiversity and abundance, and the spatial and temporal representations of eDNA with respect to its sources. To facilitate ecological applications implementing fish eDNA techniques, recent literature was summarized to generate guidelines for effective sampling in lentic, lotic, and marine habitats. Finally, we identified current gaps and limitations, and pointed out newly emerging research avenues for fish eDNA. As methodological optimization and standardization improve, eDNA technology should revolutionize fish monitoring and promote biodiversity conservation and fisheries management that transcends geographic and temporal boundaries.

Environmentally-induced sex reversal in fish with chromosomal vs. polygenic sex determination

Sex ratio depends on sex determination mechanisms and is a key demographic parameter determining population viability and resilience to natural and anthropogenic stressors. There is increasing evidence that the environment can alter sex ratio even in genetically sex-determined species (GSD), as elevated temperature can cause female-to-male sex reversal (neomales). Alarmingly, neomales are being discovered in natural populations of several fish, amphibian and reptile species worldwide. Understanding the basis of neomale development is important for conservation biology. Among GSD species, it is unknown whether those with chromosomal sex determination (CSD), the most common system, will better resist the influence of high temperature than those with polygenic sex determination (PSD). Here, we compared the effects of elevated temperature in two wild zebrafish strains, Nadia (NA) and Ekkwill (EKW), which have CSD with a ZZ/ZW system, against the AB laboratory strain, which has PSD. First, we uncovered novel sex genotypes and the results showed that, at control temperature, the masculinization rate roughly doubled with the addition of each Z chromosome, while some ZW and WW fish of the wild strains became neomales. Surprisingly, we found that at elevated temperatures WW fish were just as likely as ZW fish to become neomales and that all strains were equally susceptible to masculinization. These results demonstrate that the Z chromosome is not essential for male development and that the dose of W buffers masculinization at the control temperature but not at elevated temperature. Furthermore, at the elevated temperature the testes of neomales, but not of normal males, contained more spermatozoa than at the control temperature. Our results show in an unprecedented way that, in a global warming scenario, CSD species may not necessarily be better protected against the masculinizing effect of elevated temperature than PSD species, and reveal genotype-by-temperature interactions in male sex determination and spermatogenesis.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

Sex-specific DNA methylation: impact on human health and development

Human evolution has shaped gender differences between males and females. Over the years, scientific studies have proposed that epigenetic modifications significantly influence sex-specific differences. The evolution of sex chromosomes with epigenetics as the driving force may have led to one sex being more adaptable than the other when exposed to various factors over time. Identifying and understanding sex-specific differences, particularly in DNA methylation, will help determine how each gender responds to factors, such as disease susceptibility, environmental exposure, brain development and neurodegeneration. From a medicine and health standpoint, sex-specific methylation studies have shed light on human disease severity, progression, and response to therapeutic intervention. Interesting findings in gender incongruent individuals highlight the role of genetic makeup in influencing DNA methylation differences. Sex-specific DNA methylation studies will empower the biotechnology and pharmaceutical industry with more knowledge to identify biomarkers, design and develop sex bias drugs leading to better treatment in men and women based on their response to different diseases.

Tandem Mass Tag-Based Quantitative Proteomics Analysis of Gonads Reveals New Insight into Sexual Reversal Mechanism in Chinese Soft-Shelled Turtles

Chinese soft-shelled turtles display obvious sex dimorphism. The exogenous application of hormones (estradiol and methyltestosterone) can change the direction of gonadal differentiation of P. sinensis to produce sex reversed individuals. However, the molecular mechanism remains unclear. In this study, TMT-based quantitative proteomics analysis of four types of P. sinensis (female, male, pseudo-female, and pseudo-male) gonads were compared. Quantitative analysis of 6107 labeled proteins in the four types of P. sinensis gonads was performed. We identified 440 downregulated and 423 upregulated proteins between pseudo-females and males, as well as 394 downregulated and 959 upregulated proteins between pseudo-males and females. In the two comparisons, the differentially expressed proteins, including K7FKG1, K7GIQ2, COL4A6, K7F2U2, and K7FF80, were enriched in some important pathways, such as focal adhesion, endocytosis, apoptosis, extracellular matrix-receptor interaction, and the regulation of actin cytoskeleton, which were upregulated in pseudo-female vs. male and downregulated in pseudo-male vs. female. In pathways such as ribosome and spliceosome, the levels of RPL28, SRSF3, SNRNP40, and HNRNPK were increased from male to pseudo-female, while they decreased from female to pseudo-male. All differentially expressed proteins after sexual reversal were divided into six clusters, according to their altered levels in the four types of P. sinensis, and associated with cellular processes, such as embryonic development and catabolic process, that were closely related to sexual reversal. These data will provide clues for the sexual reversal mechanism in P. sinensis.

Epigenetic Regulation of Phenotypic Sexual Plasticity Inducing Skewed Sex Ratio in Zebrafish

The plasticity of sexual phenotype in response to environmental conditions results in biased sex ratios, and their variation has an effect on population dynamics. Epigenetic modifications can modulate sex ratio variation in species, where sex is determined by genetic and environmental factors. However, the role of epigenetic mechanisms underlying skewed sex ratios is far from being clear and is still an object of debate in evolutionary developmental biology. In this study, we used zebrafish as a model animal to investigate the effect of DNA methylation on sex ratio variation in sex-biased families in response to environmental temperature. Two sex-biased families with a significant difference in sex ratio were selected for genome-wide DNA methylation analysis using reduced representation bisulfite sequencing (RRBS). The results showed significant genome-wide methylation differences between male-biased and female-biased families, with a greater number of methylated CpG sites in testes than ovaries. Likewise, pronounced differences between testes and ovaries were identified within both families, where the male-biased family exhibited a higher number of methylated sites than the female-biased family. The effect of temperature showed more methylated positions in the high incubation temperature than the control temperature. We found differential methylation of many reproduction-related genes (e.g., sox9a, nr5a2, lhx8a, gata4) and genes involved in epigenetic mechanisms (e.g., dnmt3bb.1, dimt1l, hdac11, h1m) in both families. We conclude that epigenetic modifications can influence the sex ratio variation in zebrafish families and may generate skewed sex ratios, which could have a negative consequence for population fitness in species with genotype-environment interaction sex-determining system under rapid environmental changes.

Temporal variation in DNA methylation during gonadal development in a reptile with temperature-dependent sex determination

DNA methylation plays a significant role in transducing external environmental signals to a cellular response in reptiles; however, whether the methylation patterns are conserved across species remains unclear. Here, we examined the genome-wide DNA methylation differentiation between male and female hatchling gonads of the temperature-dependent sex determination (TSD) Mauremys mutica (M. mutica) using methylation-dependent restriction-site associated DNA sequencing (MethylRAD-seq) to test differentially methylated genes underlying sexual development. Several categories, including heat shock genes (HSP90A, HSP30C), histone- (KDM8) and ubiquitin-related genes (TRIM39), kinases (WNK3) and gonad differentiation or gonadal development related genes (HSD17B8, HSD17B12), were identified as candidates for future study. Additionally, we identified several regulatory pathways potentially mediating TSD thermosensitivity such as the GnRH signaling pathway and calcium signaling pathway. These findings provide evidence that sexually dimorphic DNA methylation may be associated with sex determination or sex differentiation in TSD M. mutica.