Epigenetic modification and inheritance in sexual reversal of fish

Dynamic epimarks in sex-related genes predict gonad phenotype in the European sea bass, a fish with mixed genetic and environmental sex determination

The integration of genomic and environmental influences into methylation patterns to bring about a phenotype is of central interest in developmental epigenetics, but many details are still unclear. The sex ratios of the species used here, the European sea bass, are determined by genetic and temperature influences. We created four families from parents known to produce offspring with different sex ratios, exposed larvae to masculinizing temperatures and examined, in juvenile gonads, the DNA methylation of seven genes related to sexual development by a targeted sequencing approach. The genes most affected by both genetics and environment were cyp19a1a and dmrt1, with contrasting sex-specific methylation and temperature responses. The relationship between cyp19a1a methylation and expression is relevant to the epigenetic regulation of vertebrate sex, and we report the evidence of such relationship only below a methylation threshold, ~ 80%, and that it was sex-specific: negatively correlated in females but positively correlated in males. From parents to offspring, the methylation in gonads was midway between oocytes and sperm, with bias towards oocytes for amh-r2, er-β2, fsh-r and cyp19a1a. In contrast, dmrt1 levels resembled those of sperm. The methylation of individual CpGs from foxl2, er-β2 and nr3c1 were conserved from parents to offspring, whereas those of cyp19a1a, dmrt1 and amh-r2 were affected by temperature. Utilizing a machine-learning procedure based on the methylation levels of a selected set of CpGs, we present the first, to our knowledge, system based on epigenetic marks capable of predicting sex in an animal with ~ 90% accuracy and discuss possible applications.

Expression analysis and characterization of zglp1 in the Chinese tongue sole (Cynoglossus semilaevis)

Zinc finger GATA like protein-1 (ZGLP1) is a nuclear zinc finger protein that regulates the interaction between somatic cells and germ cells during gonad developmental process in mammals. In this study, the zglp1 of Chinese tongue sole, Cynoglossus semilaevis (cysezglp1), was cloned and characterized for the first time in fish. Cysezglp1 had an open reading frame with five exons and was located to chromosome 9. The open reading frame of cysezglp1 consisted of 1692 nucleotides and encoded a 583 amino acid polypeptide. The predicted protein contained two zinc finger structures (Znf1 and Znf2), one of which was highly homologous to the GATA-type zinc finger domain. Multiple sequence alignment showed that Znf1 was conserved across different species while Znf2 was more divergent. Through quantitative Real-time PCR (qRT-PCR), we found that cysezglp1 was predominantly expressed in gonads, and the expression level of the ovary was significantly higher than that of the testis. We compared expression level in different embryonic stages and found that cysezglp1 mRNAs were mainly expressed in the fertilized egg to the cleavage stage, subsequently declining in the blastula stage. Cysezglp1 expression was not detected from the gastrulation stage onward. In the ovary, cysezglp1 expression was detected at 120 days after hatching and expression gradually increased with the maturation of the ovary. In situ hybridization showed that the cysezglp1 was mainly expressed in oocytes. Taken together, our results suggest that cysezglp1 may play an important role in the process of oogenesis in Chinese tongue sole.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

The DNA Methylation Landscape of Stickleback Reveals Patterns of Sex Chromosome Evolution and Effects of Environmental Salinity

Epigenetic mechanisms such as DNA methylation are a key component of dosage compensation on sex chromosomes and have been proposed as an important source of phenotypic variation influencing plasticity and adaptive evolutionary processes, yet little is known about the role of DNA methylation in an ecological or evolutionary context in vertebrates. The threespine stickleback (Gasterosteus aculeatus) is an ecological and evolutionary model system that has been used to study mechanisms involved in the evolution of adaptive phenotypes in novel environments as well as the evolution heteromorphic sex chromosomes and dosage compensation in vertebrates. Using whole genome bisulfite sequencing, we compared genome-wide DNA methylation patterns between threespine stickleback males and females and between stickleback reared at different environmental salinities. Apparent hypermethylation of the younger evolutionary stratum of the stickleback X chromosome in females relative to males suggests a potential role of DNA methylation in the evolution of heteromorphic sex chromosomes. We also demonstrate that rearing salinity has genome-wide effects on DNA methylation levels, which has the potential to lead to the accumulation of epigenetic variation between natural populations in different environments.

Sex-specific dmrt1 and cyp19a1 methylation and alternative splicing in gonads of the protandrous hermaphrodite barramundi

Epigenetics is involved in sex differentiation of gonochoristic and hermaphroditic fish species, whereby two genes dmrt1 (pro-male) and cyp19a1 (pro-female) are known to play major roles. Barramundi, Lates calcarifer, is an important tropical aquaculture species that undergo natural and permanent male to female sex change, a process for which the exact underlying molecular mechanisms are still unknown. To elucidate whether DNA methylation is involved in sex control of barramundi, a next-generation bisulfite amplicon sequencing approach was used to target 146 CpG sites within proximal promoters and first exons of seven sex-related genes (dmrt1, cyp19a1, amh, foxl2, nr5a2, sox8 and sox9) of 24 testis and 18 ovaries of captive and wild adult barramundi. Moreover, comparative expression profiles of the key dmrt1 and cyp19a1 genes were further investigated using RT-qPCR and Sanger sequencing approaches, whereas expression levels of remaining targeted genes were based on available literature for the species. Results showed that cyp19a1 and amh were more methylated in males, whereas dmrt1 and nr5a2 were more methylated in females (P < 0.001), with no gender differences found for foxl2, sox8 or sox9 genes (P > 0.05). Sex-biased promoter DNA methylation was inversely related to gene expression only for dmrt1 and nr5a2, and directly related to amh expression, whereas no differences in cyp19a1 expression were found between testes and ovaries. Notably, unique sex-specific alternative splicing of dmrt1 and cyp19a1 were discovered, whereby males lacked the full-length aromatase coding cyp19a1 mRNA due to partial or total exon splicing, and females lacked the dmrt1 exon containing the DM-domain sequence. This study advances the current knowledge aiming to elucidate the genetic mechanisms within male and female gonads of this large protandrous hermaphrodite by providing the first evidence of epigenetics and alternative splicing simultaneously affecting key genes (cyp19a1 and dmrt1) central to sex differentiation pathways.

Genome-wide identification, evolution of DNA methyltransferases and their expression during gonadal development in Nile tilapia

DNA methyltransferases (dnmts) are responsible for DNA methylation and play important roles in organism development. In this study, seven dnmts genes (dnmt1, dnmt2, dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1, dnmt3bb.2) were identified in Nile tilapia. Comprehensive analyses of dnmts were performed using available genome databases from representative animal species. Phylogenetic analysis revealed that the dnmts family were highly conserved in teleosts. Based on transcriptome data from eight adult tilapia tissues, the dnmts were found to be dominantly expressed in the head kidney, testis and ovary. Analyses of the gonadal transcriptome data in different developmental stages revealed that all dnmts were expressed in both ovary and testis, and four de novo dnmts (dnmt3aa, dnmt3ab, dnmt3bb.1, dnmt3bb.2) showed higher expression in the testis than in the ovary. Furthermore, during sex reversal induced by Fadrozole, the expression of these four de novo dnmts increased significantly in treated group compared to female control group. By in situ hybridization, the seven dnmts were found to be expressed mainly in phase I and II oocytes of the ovary and spermatocytes of the testis. When gonads were incubated with a methyltransferase inhibitor (5-AzaCdR) in vitro, the expression of dnmts genes were down-regulated significantly, while the expression of cyp19a1a (a key gene in female pathway) and dmrt1 (a key gene in male pathway) increased significantly. Our results revealed the conservation of dnmts during evolution and indicated a potential role of dnmts in epigenetic regulation of gonadal development.

New locus reveals the genetic architecture of sex reversal in the Chinese tongue sole (Cynoglossus semilaevis)

Sex reversal in insects, amphibians, reptiles, and fishes is a complicated and interesting biological phenomenon. Sex reversal changes the sex ratio of populations and may complicate breeding schemes. In the Chinese tongue sole (Cynoglossus semilaevis), genetic females may change into pseudomales, thereby increasing aquaculture costs because of the lower growth rate of the males than that of the females. Here we identify a new locus associated with sex reversal; this single nucleotide polymorphism (SNP) is located in the third intron of the doublesex and mab-3 related transcription factor 1 (Dmrt1) gene on the Z chromosome (named Cyn_Z_8564889) and has two alleles, A and G. Cyn_Z_8564889 regulates sex reversal interactively with our previously detected SNP (Cyn_Z_6676874), with the genetic females simultaneously carrying the T allele of Cyn_Z_6676874 and the A allele of Cyn_Z_8564889 changing into pseudomales. Other Dmrt1 polymorphisms were detected, which formed two haplotypes. Two SNPs in the second exon of Dmrt1 result in amino acid changes, suggesting that Dmrt1 is essential in sex reversal. We also verified that pseudomales produce no or little W sperm. The interaction and linkage between Cyn_Z_6676874 and Cyn_Z_8564889 and the absence of W sperm from pseudomales unravel the genetic architecture of sex reversal in C. semilaevis.

Origin and transition of sex determination mechanisms in a gynogenetic hexaploid fish

Most vertebrates reproduce sexually, and plastic sex determination mechanisms including genotypic sex determination (GSD) and environmental sex determination (ESD) have been extensively revealed. However, why sex determination mechanisms evolve diversely and how they correlate with diverse reproduction strategies remain largely unclear. Here, we utilize the superiority of a hexaploid gibel carp (Carassius gibelio) that is able to reproduce by unisexual gynogenesis and contains a rare but diverse proportion of males to investigate these puzzles. A total of 2248 hexaploid specimens were collected from 34 geographic wild populations throughout mainland China, in which 24 populations were revealed to contain 186 males with various incidences ranging from 1.2 to 26.5%. Subsequently, the proportion of temperature-dependent sex determination (TSD) was revealed to be positively correlated to average annual temperature in wild populations, and male incidence in lab gynogenetic progenies was demonstrated to increase with the increasing of larval rearing temperature. Meanwhile, extra microchromosomes were confirmed to play genotypic male determination role as previously reported. Thereby, GSD and TSD were found to coexist in gibel carp, and the proportions of GSD were observed to be much higher than that of TSD in sympatric wild populations. Our findings uncover a potential new mechanism in the evolution of sex determination system in polyploid vertebrates with unisexual gynogenesis ability, and also reveal a possible association of sex determination mechanism transition between TSD and GSD and reproduction mode transition between unisexual gynogenesis and bisexual reproduction.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

Epigenetic regulation of male fate commitment from an initially bipotential system

A fundamental goal in biology is to understand how distinct cell types containing the same genetic information arise from a single stem cell throughout development. Sex determination is a key developmental process that requires a unidirectional commitment of an initially bipotential gonad towards either the male or female fate. This makes sex determination a unique model to study cell fate commitment and differentiation in vivo. We have focused this review on the accumulating evidence that epigenetic mechanisms contribute to the bipotential state of the fetal gonad and to the regulation of chromatin accessibility during and immediately downstream of the primary sex-determining switch that establishes the male fate.

Distinct sperm nucleus behaviors between genotypic and temperature-dependent sex determination males are associated with replication and expression-related pathways in a gynogenetic fish

BACKGROUND

Coexistence and transition of diverse sex determination strategies have been revealed in some ectothermic species, but the variation between males caused by different sex determination strategies and the underlying mechanism remain unclear. Here, we used the gynogenetic gibel carp (Carassius gibelio) with both genotypic sex determination (GSD) and temperature-dependent sex determination (TSD) strategies to illustrate this issue.

RESULTS

We found out that males of GSD and TSD in gibel carp had similar morphology, testicular histology, sperm structure and sperm vitality. However, when maternal individuals were mated with males of GSD, sperm nucleus swelling and fusing with the female pronucleus were observed in the fertilized eggs. On the contrary, when maternal individuals were mated with males of TSD, sperm nucleus remained in the condensed status throughout the whole process. Subsequently, semen proteomics analysis unveiled that DNA replication and gene expression-related pathways were inhibited in the sperm from males of TSD compared to males of GSD, and most differentially expressed proteins associated with DNA replication, transcription and translation were down-regulated. Moreover, via BrdU incorporation and immunofluorescence detection, male nucleus replication was revealed to be present in the fertilized eggs by the sperm from males of GSD, but absent in the fertilized eggs by the sperm from males of TSD.

CONCLUSIONS

These findings indicate that DNA replication and gene expression-related pathways are associated with the distinct sperm nucleus development behaviors in fertilized eggs in response to the sperm from males of GSD and TSD. And this study is the first attempt to screen the differences between males determined via GSD and TSD in gynogenetic species, which might give a hint for understanding evolutionary adaption of diverse sex determination mechanisms in unisexual vertebrates.

Climate-driven shifts in adult sex ratios via sex reversals: the type of sex determination matters

Sex reversals whereby individuals of one genetic sex develop the phenotype of the opposite sex occur in ectothermic vertebrates with genetic sex-determination systems that are sensitive to extreme temperatures during sexual differentiation. Recent rises in global temperatures have led researchers to predict that sex reversals will become more common, resulting in the distortion of many populations' sex ratios. However, it is unclear whether susceptibility to climate-driven sex-ratio shifts depends on the type of sex determination that varies across species. First, we show here using individual-based theoretical models that XX/XY (male-heterogametic) and ZZ/ZW (female-heterogametic) sex-determination systems can respond differentially to temperature-induced sex reversals. Interestingly, the impacts of climate warming on adult sex ratio (ASR) depend on the effects of both genotypic and phenotypic sex on survival and reproduction. Second, we analyse the temporal changes of ASR in natural amphibian populations using data from the literature, and find that ASR shifted towards males in ZZ/ZW species over the past 60 years, but did not change significantly in XX/XY species. Our results highlight the fact that we need a better understanding of the interactions between genetic and environmental sex-determining mechanisms to predict the responses of ectotherms to climate change and the associated extinction risks.This article is part of the themed issue 'Adult sex ratios and reproductive decisions: a critical re-examination of sex differences in human and animal societies'.

Gonad Transcriptome Analysis of High-Temperature-Treated Females and High-Temperature-Induced Sex-Reversed Neomales in Nile Tilapia

Background: Nowadays, the molecular mechanisms governing TSD (temperature-dependent sex determination) or GSD + TE (genotypic sex determination + temperature effects) remain a mystery in fish. Methods: We developed three all-female families of Nile tilapia (Oreochromis niloticus), and the family with the highest male ratio after high-temperature treatment was used for transcriptome analysis. Results: First, gonadal histology analysis indicated that the histological morphology of control females (CF) was not significantly different from that of high-temperature-treated females (TF) at various development stages. However, the high-temperature treatment caused a lag of spermatogenesis in high-temperature-induced neomales (IM). Next, we sequenced the transcriptome of CF, TF, and IM Nile tilapia. 79, 11,117, and 11,000 differentially expressed genes (DEGs) were detected in the CF–TF, CF–IM, and TF–IM comparisons, respectively, and 44 DEGs showed identical expression changes in the CF–TF and CF–IM comparisons. Principal component analysis (PCA) indicated that three individuals in CF and three individuals in TF formed a cluster, and three individuals in IM formed a distinct cluster, which confirmed that the gonad transcriptome profile of TF was similar to that of CF and different from that of IM. Finally, six sex-related genes were validated by qRT-PCR. Conclusions: This study identifies a number of genes that may be involved in GSD + TE, which will be useful for investigating the molecular mechanisms of TSD or GSD + TE in fish.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.

Treatment with a DNA methyltransferase inhibitor feminizes zebrafish and induces long-term expression changes in the gonads

BACKGROUND

The role of epigenetic modifications such as DNA methylation during vertebrate sexual development is far from being clear. Using the zebrafish model, we tested the effects of one of the most common DNA methyltransferase (dnmt) inhibitor, 5-aza-2'-deoxycytidine (5-aza-dC), which is approved for the treatment of acute myeloid leukaemia and is under active investigation for the treatment of solid tumours. Several dose-response experiments were carried out during two periods, including not only the very first days of development (0-6 days post-fertilization, dpf), as done in previous studies, but also, and as a novelty, the period of gonadal development (10-30 dpf).

RESULTS

Early treatment with 5-aza-dC altered embryonic development, delayed hatching and increased teratology and mortality, as expected. The most striking result, however, was an increase in the number of females, suggesting that alterations induced by 5-aza-dC treatment can affect sexual development as well. Results were confirmed when treatment coincided with gonadal development. In addition, we also found that the adult gonadal transcriptome of 5-aza-dC-exposed females included significant changes in the expression of key reproduction-related genes (e.g. cyp11a1, esr2b and figla), and that several pro-female-related pathways such as the Fanconi anaemia or the Wnt signalling pathways were downregulated. Furthermore, an overall inhibition of genes implicated in epigenetic regulatory mechanisms (e.g. dnmt1, dicer, cbx4) was also observed.

CONCLUSIONS

Taken together, our results indicate that treatment with a DNA methylation inhibitor can also alter the sexual development in zebrafish, with permanent alterations of the adult gonadal transcriptome, at least in females. Our results show the importance of DNA methylation for proper control of sexual development, open new avenues for the potential control of sex ratios in fish (aquaculture, population control) and call attention to possibly hidden long-term effects of dnmt therapy when used, for example, in the treatment of prepuberal children affected by some types of cancer.

Epigenetic considerations in aquaculture

Epigenetics has attracted considerable attention with respect to its potential value in many areas of agricultural production, particularly under conditions where the environment can be manipulated or natural variation exists. Here we introduce key concepts and definitions of epigenetic mechanisms, including DNA methylation, histone modifications and non-coding RNA, review the current understanding of epigenetics in both fish and shellfish, and propose key areas of aquaculture where epigenetics could be applied. The first key area is environmental manipulation, where the intention is to induce an ‘epigenetic memory’ either within or between generations to produce a desired phenotype. The second key area is epigenetic selection, which, alone or combined with genetic selection, may increase the reliability of producing animals with desired phenotypes. Based on aspects of life history and husbandry practices in aquaculture species, the application of epigenetic knowledge could significantly affect the productivity and sustainability of aquaculture practices. Conversely, clarifying the role of epigenetic mechanisms in aquaculture species may upend traditional assumptions about selection practices. Ultimately, there are still many unanswered questions regarding how epigenetic mechanisms might be leveraged in aquaculture.

Rainbow trout exposed to benzo[a]pyrene yields conserved microRNA binding sites in DNA methyltransferases across 500 million years of evolution

The objective of this study was to examine the regulation of DNA methylation following acute (24 h) and prolonged (14 d) exposure to low (1 ng/L) and high (10 ng/L) benzo[a]pyrene. However, with the recent release of the rainbow trout genome, we were able to conduct a more detailed analysis regarding the regulation of the enzymes involved in DNA methylation; DNA methyltransferases (DNMTs). Bioinformatic approaches were used to identify candidate microRNA (miRNA) that potentially bind to the DNMT1 and DNMT3a 3′UTR. Results indicated a significant decrease in global methylation in both liver and muscle, with an associated decrease in DNA methyltransferase activity and DNMT3a transcript abundance. There was a significant increase in one specific candidate miRNA (miR29a) that was predicted to bind to DNMT3a. Taking a comparative genomics approach, the binding sites of miR29a to the DNMT3a 3′UTR was compared across species, spanning fish to mammals, and revealed a highly conserved binding motif that has been maintained since the vertebrate ancestor, approximately 500 million years ago. This research establishes that miRNA act as an essential mediator between the environment and DNA methylation patterns via DNMTs, which is further confirmed by a genomic regulatory mechanism that has been deeply conserved throughout evolution.

Molecular characterization of Pod1 during sex development in Chinese tongue sole (Cynoglossus semilaevis)

Pod1 encodes a Class II bHLH transcription factor involved in the development of a number of tissues such as gonad, spleen, lungs and heart. However, to date, little is known about its function in teleosts. In this study, we cloned and characterized Pod1 gene from Cynoglossus semilaevis. This gene contains three exons and two introns, with the full-length cDNA of 918 nucleotides that encodes a 183 amino acid protein with a conserved bHLH domain. Realtime quantitative PCR revealed that Pod1 was predominantly expressed in the testes of C. semilaevis. In different stages of testes development, Pod1 expression was undetectable up to 120 days after hatching (dah), and then increased at 210 dah and 1 year after hatching (yah). Furthermore, in situ hybridization (ISH) analysis revealed that Pod1 was mainly localized in the germ cells of testes, but was not detected in ovarian cells; which suggested its possible functions in spermatogenesis of C. semilaevis. The methylation profile analysis of Pod1 genomic sequence in the gonads showed that the differences in their putative promoter regions of Pod1 among ovary, male and pseudo-male testes were not obvious. Thus, further research might be needed to evaluate whether Pod1 expression is regulated by epigenetic level.

Environmental epigenetics in zebrafish

It is widely accepted that the epigenome can act as the link between environmental cues, both external and internal, to the organism and phenotype by converting the environmental stimuli to phenotypic responses through changes in gene transcription outcomes. Environmental stress endured by individual organisms can also enforce epigenetic variations in offspring that had never experienced it directly, which is termed transgenerational inheritance. To date, research in the environmental epigenetics discipline has used a wide range of both model and non-model organisms to elucidate the various epigenetic mechanisms underlying the adaptive response to environmental stimuli. In this review, we discuss the advantages of the zebrafish model for studying how environmental toxicant exposures affect the regulation of epigenetic processes, especially DNA methylation, which is the best-studied epigenetic mechanism. We include several very recent studies describing the state-of-the-art knowledge on this topic in zebrafish, together with key concepts in the function of DNA methylation during vertebrate embryogenesis.

Small ocean temperature increases elicit stage-dependent changes in DNA methylation and gene expression in a fish, the European sea bass

In natural fish populations, temperature increases can result in shifts in important phenotypic traits. DNA methylation is an epigenetic mechanism mediating phenotypic changes. However, whether temperature increases of the magnitude predicted by the latest global warming models can affect DNA methylation is unknown. Here, we exposed European sea bass to moderate temperature increases in different periods within the first two months of age. We show that increases of even 2 °C in larvae significantly changed global DNA methylation and the expression of ecologically-relevant genes related to DNA methylation, stress response, muscle and organ formation, while 4 °C had no effect on juveniles. Furthermore, DNA methylation changes were more marked in larvae previously acclimated to a different temperature. The expression of most genes was also affected by temperature in the larvae but not in juveniles. In conclusion, this work constitutes the first study of DNA methylation in fish showing that temperature increases of the magnitude predicted by the latest global warming models result in stage-dependent alterations in global DNA methylation and gene expression levels. This study, therefore, provides insights on the possible consequences of climate change in fish mediated by genome-wide epigenetic modifications.

Stress and sex: does cortisol mediate sex change in fish?

Cortisol is the main glucocorticoid (GC) in fish and the hormone most directly associated with stress. Recent research suggests that this hormone may act as a key factor linking social environmental stimuli and the onset of sex change by initiating a shift in steroidogenesis from estrogens to androgens. For many teleost fish, sex change occurs as a usual part of the life cycle. Changing sex is known to enhance the lifetime reproductive success of these fish and the modifications involved (behavioral, gonadal and morphological) are well studied. However, the exact mechanism behind the transduction of the environmental signals into the molecular cascade that underlies this singular process remains largely unknown. We here synthesize current knowledge regarding the role of cortisol in teleost sex change with a focus on two well-described transformations: temperature-induced masculinization and socially regulated sex change. Three non-mutually exclusive pathways are considered when describing the potential role of cortisol in mediating teleost sex change: cross-talk between GC and androgen pathways, inhibition of aromatase expression and upregulation of amh (the gene encoding anti-Müllerian hormone). We anticipate that understanding the role of cortisol in the initial stages of sex change will further improve our understanding of sex determination and differentiation across vertebrates, and may lead to new tools to control fish sex ratios in aquaculture.

Effects of different dietary DHA:EPA ratios on gonadal steroidogenesis in the marine teleost, tongue sole (Cynoglossus semilaevis)

The present study was conducted to investigate the effects of dietary DHA and EPA on gonadal steroidogenesis in mature females and males, with a feeding trial on tongue sole, a typical marine teleost with sexual dimorphism. Three experimental diets differing basically in DHA:EPA ratio, that is, 0·68 (diet D:E-0·68), 1·09 (D:E-1·09) and 2·05 (D:E-2·05), were randomly assigned to nine tanks of 3-year-old tongue sole (ten females and fifteen males in each tank). The feeding trail lasted for 90 d before and during the spawning season. Fish were reared in a flowing seawater system and fed to apparent satiation twice daily. Compared with diet D:E-0·68, diet D:E-1·09 significantly enhanced the oestradiol production in females, whereas diet D:E-2·05 significantly enhanced the testosterone production in males. In ovaries, diet D:E-1·09 induced highest mRNA expression of follicle-stimulating hormone receptor (FSHR), steroidogenic acute regulatory protein, 17α-hydroxylase (P450c17) and 3β-hydroxysteroid dehydrogenase (3β-HSD). In testes, diet 2·05 resulted in highest mRNA expression of FSHR, cholesterol side-chain cleavage enzyme, P450c17 and 3β-HSD. Fatty acid profiles in fish tissues reflected closely those of diets. Female fish had more gonadal EPA content but less DHA content than male fish, whereas there was a reverse observation in liver. In conclusion, the dietary DHA:EPA ratio, possibly combined with the dietary EPA:arachidonic acid ratio, differentially regulated sex steroid hormone synthesis in mature female and male tongue soles. Females seemed to require more EPA but less DHA for the gonadal steroidogenesis than males. The results are beneficial to sex-specific nutritive strategies in domestic teleost.

The destiny of the resistance/susceptibility against GCRV is controlled by epigenetic mechanisms in CIK cells

Hemorrhagic disease caused by grass carp reovirus (GCRV) has severely threatened the grass carp (Ctenopharyngodon idella) cultivation industry. It is noteworthy that the resistance against GCRV infection was reported to be inheritable, and identified at both individual and cellular levels. Therefore, this work was inspired and dedicated to unravel the molecular mechanisms of fate decision post GCRV infection in related immune cells. Foremost, the resistant and susceptible CIK (C. idella kidney) monoclonal cells were established by single cell sorting, subculturing and infection screening successively. RNA-Seq, MeDIP-Seq and small RNA-Seq were carried out with C1 (CIK cells), R2 (resistant cells) and S3 (susceptible cells) groups. It was demonstrated that genome-wide DNA methylation, mRNA and microRNA expression levels in S3 were the highest among three groups. Transcriptome analysis elucidated that pathways associated with antioxidant activity, cell proliferation regulation, apoptosis activity and energy consuming might contribute to the decision of cell fates post infection. And a series of immune-related genes were identified differentially expressed across resistant and susceptible groups, which were negatively modulated by DNA methylation or microRNAs. To conclude, this study systematically uncovered the regulatory mechanism on the resistance from epigenetic perspective and provided potential biomarkers for future studies on resistance breeding.

Hatching enzymes disrupt aberrant gonadal degeneration by the autophagy/apoptosis cell fate decision

Environmental stressors, gonadal degenerative diseases and tumour development can significantly alter the oocyte physiology, and species fertility and fitness. To expand the molecular understanding about oocyte degradation, we isolated several spliced variants of Japanese anchovy hatching enzymes (AcHEs; ovastacin homologue) 1 and 2, and analysed their potential in oocyte sustenance. Particularly, AcHE1b, an ovary-specific, steroid-regulated, methylation-dependent, stress-responsive isoform, was neofunctionalized to regulate autophagic oocyte degeneration. AcHE1a and 2 triggered apoptotic degeneration in vitellogenic and mature oocytes, respectively. Progesterone, starvation, and high temperature elevated the total degenerating oocyte population and AcHE1b transcription by hyper-demethylation. Overexpression, knockdown and intracellular zinc ion chelation study confirmed the functional significance of AcHE1b in autophagy induction, possibly to mitigate the stress effects in fish, via ion-homeostasis. Our finding chronicles the importance of AcHEs in stress-influenced apoptosis/autophagy cell fate decision and may prove significant in reproductive failure assessments, gonadal health maintenance and ovarian degenerative disease therapy.

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination

Background

DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.

Results

Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.

Conclusions

Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0136-2) contains supplementary material, which is available to authorized users.

Dmrt1 induces the male pathway in a turtle species with temperature-dependent sex determination

The molecular mechanism underlying temperature-dependent sex determination (TSD) has been a long-standing mystery; in particular, the thermosensitive genetic triggers for gonadal sex differentiation are largely unknown. Here, we have characterized a conserved DM domain gene, Dmrt1, in the red-eared slider turtle Trachemys scripta (T. scripta), which exhibits TSD. We found that Dmrt1 has a temperature-dependent, sexually dimorphic expression pattern, preceding gonadal sex differentiation, and is capable of responding rapidly to temperature shifts and aromatase inhibitor treatment. Most importantly, loss- and gain-of-function analyses provide solid evidence that Dmrt1 is both necessary and sufficient to initiate male development in T. scripta Furthermore, the DNA methylation dynamics of the Dmrt1 promoter are tightly correlated with temperature and could mediate the impact of temperature on sex determination. Collectively, our findings demonstrate that Dmrt1 is a candidate master male sex-determining gene in this TSD species, consistent with the idea that DM domain genes are conserved during the evolution of sex determination mechanisms.

Identification and analysis of the β-catenin1 gene in half-smooth tongue sole (Cynoglossus semilaevis)

β-catenin is a key signalling molecule in the canonical Wnt pathway, which plays a role in cell adhesion, embryogenesis and sex determination. However, little is known about its function in teleosts. We cloned and characterized the full-length β-catenin1 gene from half-smooth tongue sole (Cynoglossus semilaevis), which was designated CS-β-catenin1. The CS-β-catenin1 cDNA consists of 2,346 nucleotides and encodes a protein with 782 amino acids. Although CS-β-catenin1 was transcribed in the gonads of both sexes, the level was significantly higher in ovaries compared to testes. Furthermore, the mRNA level of CS-β-catenin1 was significantly upregulated at 160 days and constantly increased until 2 years of age. In situ hybridization revealed that CS-β-catenin1 mRNA was mainly localized in oocyte cells, especially in stage I, II and III oocytes. When CS-β-catenin1 expression was inhibited by injection of quercetin in the ovaries, levels of CS-Figla and CS-foxl2 mRNA were significantly down-regulated, and CS-dmrt1 was up-regulated, which suggested that CS-β-catenin1 is a potential upstream gene of CS-Figla and is involved in the development of the ovaries, i.e., folliculogenesis.

The Effects of CpG Densities around Transcription Start Sites on Sex-Biased Gene Expression in Poecilia reticulata

As most genes are shared between females and males, DNA methylation is assumed to play a crucial role in sex-biased gene expression. DNA methylation exclusively occurs at CpG dinucleotides, and therefore, we would expect that CpG density around transcription start sites (TSSs) relate to sex-biased gene expression. Here we investigated the relationship between CpG densities around TSSs and the ratio of gene expression levels between sexes in the guppy (Poecilia reticulata), which displays remarkable sexual dimorphisms. We found that genes with sex-biased gene expression had different CpG densities downstream of TSSs compared with genes lacking sex-biased gene expression. Intriguingly, male-biased expression genes with intermediate CpG density downstream of TSSs exhibited greater differences in gene expression between sexes in the gonad and tail. Our findings suggested the possibility that CpGs around TSSs, especially in the downstream regions, play a crucial role in sex-biased gene expression through DNA methylation.

The elephant shark methylome reveals conservation of epigenetic regulation across jawed vertebrates

BACKGROUND

Methylation of CG dinucleotides constitutes a critical system of epigenetic memory in bony vertebrates, where it modulates gene expression and suppresses transposon activity. The genomes of studied vertebrates are pervasively hypermethylated, with the exception of regulatory elements such as transcription start sites (TSSs), where the presence of methylation is associated with gene silencing. This system is not found in the sparsely methylated genomes of invertebrates, and establishing how it arose during early vertebrate evolution is impeded by a paucity of epigenetic data from basal vertebrates.

METHODS

 We perform whole-genome bisulfite sequencing to generate the first genome-wide methylation profiles of a cartilaginous fish, the elephant shark Callorhinchus milii. Employing these to determine the elephant shark methylome structure and its relationship with expression, we compare this with higher vertebrates and an invertebrate chordate using published methylation and transcriptome data.  Results: Like higher vertebrates, the majority of elephant shark CG sites are highly methylated, and methylation is abundant across the genome rather than patterned in the mosaic configuration of invertebrates. This global hypermethylation includes transposable elements and the bodies of genes at all expression levels. Significantly, we document an inverse relationship between TSS methylation and expression in the elephant shark, supporting the presence of the repressive regulatory architecture shared by higher vertebrates.

CONCLUSIONS

 Our demonstration that methylation patterns in a cartilaginous fish are characteristic of higher vertebrates imply the conservation of this epigenetic modification system across jawed vertebrates separated by 465 million years of evolution. In addition, these findings position the elephant shark as a valuable model to explore the evolutionary history and function of vertebrate methylation.

Identification of exosomes and its signature miRNAs of male and female Cynoglossus semilaevis

Exosomes are small membrane particles which are widely found in various cell lines and physiological fluids in mammalian. MicroRNAs (miRNAs) enclosed in exosomes have been identified as proper signatures for many diseases and response to therapies. However, the composition of exosomes and enclosed miRNAs in fishes has not been investigated. Cynoglossus semilaevis is an important commercial flatfish with ambiguous distinction between males and females before sex maturation, which leads to screening difficulty in reproduction and cultivation. An effective detection method was required for sex differentiation of C. semilaevis. In this work, we successfully identified exosomes in C. semilaevis serum. The analysis of nucleotide composition showed that miRNA dominated in exosomes. Thereafter the miRNA profiles in exosomes from males and females were sequenced and compared to identify the signature miRNAs corresponding to sex differentiation. The functions of signature miRNAs were analyzed by target matching and annotation. Furthermore, 7 miRNAs with high expression in males were selected from signature miRNAs as the markers for sex identification with their expression profiles verified by real time quantitative PCR. Exosomes were first found in fish serum in this work. Investigation of marker miRNAs supplies an effective index for the filtration of male and female C. semilaevis in cultivation.

Genome editing reveals dmrt1 as an essential male sex-determining gene in Chinese tongue sole (Cynoglossus semilaevis)

Chinese tongue sole is a marine fish with ZW sex determination. Genome sequencing suggested that the Z-linked dmrt1 is a putative male determination gene, but direct genetic evidence is still lacking. Here we show that TALEN of dmrt1 efficiently induced mutations of this gene. The ZZ dmrt1 mutant fish developed ovary-like testis, and the spermatogenesis was disrupted. The female-related genes foxl2 and cyp19a1a were significantly increased in the gonad of the ZZ dmrt1 mutant. Conversely, the male-related genes Sox9a and Amh were significantly decreased. The dmrt1 deficient ZZ fish grew much faster than ZZ male control. Notably, we obtained an intersex ZW fish with a testis on one side and an ovary on the other side. This fish was chimeric for a dmrt1 mutation in the ovary, and wild-type dmrt1 in the testis. Our data provide the first functional evidence that dmrt1 is a male determining gene in tongue sole.

Single Locus Maintains Large Variation of Sex Reversal in Half-Smooth Tongue Sole (Cynoglossus semilaevis)

Sex determination is a fundamental biological process for individual sex development and population sex ratios. However, for some species, the primary sex might be altered during development, and individuals can develop into the opposite sex. Sex reversal may happen in insects, reptiles, amphibians, and fishes. In half-smooth tongue sole (Cynoglossus semilaevis), some genetically female fish irreversibly reverse to pseudomales, resulting in higher costs in aquaculture owing to a lower growth rate of male fish during a 2-yr growth period. Here, we identified a locus with large controlling effect on sex reversal in the half-smooth tongue sole through genome-wide association study with high-density single nucleotide polymorphisms (SNPs). This SNP is located at the third intron of the F-box and leucine rich repeat protein 17 (FBXL17) gene on the Z chromosome, and it has two alleles, A and T. Genetic females with Z^AW genotypes will never reverse into phenotypic males, but those with Z^TW genotypes can sometimes undergo sex reversal. This SNP explains 82.7% of the genetic variation, or 58.4% of the phenotypic variation. Based on our results, a reproductive management program could be developed to improve the phenotypic female ratio in aquaculture, and elucidate the mechanism of sex reversal in half-smooth tongue sole. We expect that these findings will have a substantial impact on the population management in many harvested species where sex reversal occurs.

Heat-induced masculinization in domesticated zebrafish is family-specific and yields a set of different gonadal transcriptomes

Understanding environmental influences on sex ratios is important for the study of the evolution of sex-determining mechanisms and for evaluating the effects of global warming and chemical pollution. Fishes exhibit sexual plasticity, but the underlying mechanisms of environmental effects on their reproduction are unclear even in the well-established teleost research model, the zebrafish. Here we established the conditions to study the effects of elevated temperature on zebrafish sex. We showed that sex ratio response to elevated temperature is family-specific and typically leads to masculinization (female-to-male sex reversal), resulting in neomales. These results uncovered genotype-by-environment interactions that support a polygenic sex determination system in domesticated (laboratory) zebrafish. We found that some heat-treated fish had gene expression profiles similar to untreated controls of the same sex, indicating that they were resistant to thermal effects. Further, most neomales had gonadal transcriptomes similar to that of regular males. Strikingly, we discovered heat-treated females that displayed a normal ovarian phenotype but with a "male-like" gonadal transcriptome. Such major transcriptomic reprogramming with preserved organ structure has never been reported. Juveniles were also found to have a male-like transcriptome shortly after exposure to heat. These findings were validated by analyzing the expression of genes and signaling pathways associated with sex differentiation. Our results revealed a lasting thermal effect on zebrafish gonads, suggesting new avenues for detection of functional consequences of elevated temperature in natural fish populations in a global warming scenario.

The Reversible Sex of Gonochoristic Fish: Insights and Consequences

Fish sex reversal is a means to understand sex determination and differentiation, but it is also used to control sex in aquaculture. This review discusses sex reversal in gonochoristic fish, with the coexistence of genetic and environmental influences. The different periods of fish sensitivity to sex reversal treatments are presented with the mechanisms implicated. The old players of sex differentiation are revisited with transcriptome data and loss of function studies following hormone- or temperature-induced sex reversal. We also discuss whether cortisol is the universal mediator of sex reversal in fish due to its implication in ovarian meiosis and 11KT increase. The large plasticity in fish for sex reversal is also evident in the brain, with a reversibility existing even in adulthood. Studies on epigenetics are presented, since it links the environment, gene expression, and sex reversal, notably the association of DNA methylation in sex reversal. Manipulations with exogenous factors reverse the primary sex in many fish species under controlled conditions, but several questions arise on whether this can occur under wild conditions and what is the ecological significance. Cases of sex reversal in wild fish populations are shown and their fitness and future perspectives are discussed.

Bending Genders: The Biology of Natural Sex Change in Fish

Sexual fate is no longer seen as an irreversible deterministic switch set during early embryonic development but as an ongoing battle for primacy between male and female developmental trajectories. That sexual fate is not final and must be actively maintained via continuous suppression of the opposing sexual network creates the potential for flexibility into adulthood. In many fishes, sexuality is not only extremely plastic, but sex change is a usual and adaptive part of the life cycle. Sequential hermaphrodites begin life as one sex, changing sometime later to the other, and include species capable of protandrous (male-to-female), protogynous (female-to-male), or serial (bidirectional) sex change. Natural sex change involves coordinated transformations across multiple biological systems, including behavioural, anatomical, neuroendocrine, and molecular axes. We here review the biological processes underlying this amazing transformation, focussing particularly on its molecular basis, which remains poorly understood, but where new genomic technologies are significantly advancing our understanding of how sex change is initiated and progressed at the molecular level. Knowledge of how a usually committed developmental process remains plastic in sequentially hermaphroditic fishes is relevant to understanding the evolution and functioning of sexual developmental systems in vertebrates generally, as well as pathologies of sexual development in humans.