The Genetics and Epigenetics of Sex Change in Fish

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

Protein expression patterns and metal metabolites in a protogynous hermaphrodite fish, the ricefield eel (Monopterus albus)

BACKGROUND

The ricefield eel Monopterus albus undergoes a natural sex change from female to male during its life cycle, and previous studies have shown the potential mechanisms of this transition at the transcriptional and protein levels. However, the changes in protein levels have not been fully explored, especially in the intersexual stage.

RESULTS

In the present study, the protein expression patterns in the gonadal tissues from five different periods, the ovary (OV), early intersexual stage gonad (IE), middle intersexual stage gonad (IM), late intersexual stage gonad (IL), and testis (TE), were determined by untargeted proteomics sequencing. A total of 5125 proteins and 394 differentially expressed proteins (DEPs) were detected in the gonadal tissues. Of the 394 DEPs, there were 136 between the OV and IE groups, 20 between the IM and IE groups, 179 between the IL and IM groups, and 59 between the TE and IL groups. Three candidate proteins, insulin-like growth factor 2 mRNA-binding protein 3 isoform X1 (Igf2bp3), triosephosphate isomerase (Tpi), and Cu-Zn superoxide dismutase isoform X1 [(Cu-Zn) Sod1], were validated by western blotting to verify the reliability of the data. Furthermore, metal metabolite-related proteins were enriched in the IL vs. IM groups and TE vs. IL groups, which had close relationships with sex change, including Cu2+-, Ca2+-, Zn2+- and Fe2+/Fe3+-related proteins. Analysis of the combined transcriptome data revealed consistent protein/mRNA expression trends for two metal metabolite-related proteins/genes [LOC109953912 and calcium Binding Protein 39 Like (cab39l)]. Notably, we detected significantly higher levels of Cu2+ during the sex change process, suggesting that Cu2+ is a male-related metal metabolite that may have an important function in male reproductive development.

CONCLUSIONS

In summary, we analyzed the protein profiles of ricefield eel gonadal tissues in five sexual stages (OV, IE, IM, IL, and TE) and verified the plausibility of the data. After preforming the functional enrichment of metal metabolite-related DEPs, we detected the contents of the metal metabolites Zn2+, Cu2+, Ca2+, and Fe2+/Fe3+ at these five stages and screened for (Cu-Zn) Sod1 and Mmp-9 as possible key proteins in the sex reversal process.

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.

Reference transcriptome assembly of a protogynous sex change fish, harlequin sandsmelt (Parapercis pulchella)

The harlequin sandsmelt (Parapercis pulchella) is a female-to-male sex change fish in which functional females possess ovotestes that consist of both ovarian and testicular tissues. These features indicate that this species could be an excellent model for studying the flexibility of sex differentiation in vertebrates. However, genetic resources in this species have so far been limited. Therefore, in this study, the reference transcriptome of this fish was constructed through RNA-sequencing, de novo transcriptome assembly, superTranscripts construction, and functional annotations. To obtain as many genes as possible, RNA was extracted from various tissues (brains, gills, hearts, livers, guts, and gonads) and various sexual stages (females, individuals during sex change, and males) and then subjected to sequencing and downstream analyses. As a result, 91,884 representative transcripts with 32,627 protein-coding sequences were generated. 72.2% of protein-coding sequences (23,566 sequences) were functionally annotated. Also, our analysis shows that the superTranscripts method effectively removes redundant sequences from raw-assembled data compared with other strategies. The resultant dataset is a valuable resource for future molecular developmental studies on sex change in P. pulchella.

Chromosome-level genome assembly of humpback grouper using PacBio HiFi reads and Hi-C technologies

The humpback grouper (Cromileptes altivelis), a medium-sized coral reef teleost, is a naturally rare species distributed in the tropical waters of the Indian and Pacific Oceans. It has high market value, but artificial reproduction and breeding remain limited and need to be improved. Here, we assembled the genome with 1.08 Gb, with a contig N50 of 43.78 Mb. A total of 96.59% of the assembly anchored to 24 pseudochromosomes using Hi-C technology. It contained 24,442 protein-coding sequences, of which 99.3% were functionally annotated. The completeness of the assembly was estimated to be 97.3% using BUSCO. The phylogenomic analysis suggested that humpback grouper should be classified into the genus Epinephelus rather than Cromileptes. The comparative genomic analysis revealed that the gene families related to circadian entrainment were significantly expanded. The high-quality reference genome provides useful genomic tools for exploiting the genomic resource of humpback grouper and supports the functional genomic study of this species in the future.

Bidirectional sex-change behavior and physiological aspects in the Gorgeous goby Lythrypnus pulchellus (Gobiidae)

The Gorgeous goby Lythrypnus pulchellus shows extreme sexual plasticity with the bidirectional sex-change ability socially controlled in adults. Therefore, this study describes how the hierarchical status affects hormone synthesis through newborn hormone waste products in water and tests the influence of body size and social dominance establishment in sex reversal duration and direction. The associated changes in behavior and hormone levels are described under laboratory conditions in male-male and female-female pairs of similar and different body sizes, recording the changes until spawning. The status establishment occurred in a relatively shorter time period in male and female pairs of different sizes (1-3 days) compared to those of similar size (3-5 days), but the earlier one did not significantly affect the overall time of sex change (verified by pair spawning). The changes in gonads, hormones, and papilla occurred in sex-changer individuals, but the first one was observed in behavior. Courtship started at 3-5 days in male pairs and from 2 h to 1 day in female pairs of both groups of different and similar sizes. Hormones did not gradually move in the new sexual phenotype direction during the sex-change time course. Nonetheless, estradiol regulated sex change and 11-ketotestosterone enabled bidirectional sex change and was modulated by agonistic interactions. Cortisol is associated with status and gonadal sex change. In general, similar mechanisms underlie sex change in both directions with a temporal change sequence in phases. These results shed new light on sex-change mechanisms. Further studies should be performed to determine whether these localized changes exist in the steroid hormone synthesis along the brain-pituitary gonad axis during social and bidirectional sex changes in L. pulchellus.

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.

Transcriptomic signatures associated with underlying rapid changes in the early phase brain of bi-directional sex change in Trimma okinawae

Teleost fish exhibit remarkable sexual plasticity and divergent developmental systems, including sequential hermaphroditism. One of the more fascinating models of sexual plasticity is socially controlled sex change, which is often observed in coral reef fish. The Okinawa rubble goby, Trimma okinawae, is a bi-directional sex-changing fish. It can rapidly change sex in either direction based on social circumstances. Although behavioural and neuroendocrine sex change occurs immediately and is believed to trigger gonadal changes, the underlying mechanisms remain poorly understood. In this study, we conducted a de novo transcriptome analysis of the T. okinawae brain and identified genes that are differentially expressed between the sexes and genes that were immediately controlled by social stimulation implicating sex change. Several genes showed concordant expression shifts regardless of the sex change direction and were associated with histone modification in nerve cells. These genes are known to function in the neuroendocrine control of reproduction in nerve cells. Overall, we identified genes associated with the initiation of sex change, which provides insight into the regulation of sex change and sexual plasticity.

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.

A Strategy to Characterize the Global Landscape of Histone Post-Translational Modifications Within Tissues of Nonmodel Organisms

Histone post-translational modifications (PTMs) are epigenetic marks that play a critical role in the expression and maintenance of DNA, but they remain largely uninvestigated in nonmodel organisms due to technical challenges. To begin alleviating this issue, we developed a workflow for histone PTM analysis in Mozambique tilapia (Oreochromis mossambicus), being a widespread and environmentally hardy fish, using mass spectrometry methods. By incorporating multiple protein digestion methods into the preparation of each sample, we reliably quantified 214 biologically relevant histone PTMs. All of these histone PTMs, collectively referred to as the global histone PTM landscape, were characterized in the gills, kidney, and testes of this fish. By comparing the global histone PTM landscape between the three tissues, we found that 91.59% of histone PTMs were tissue-dependent. The workflow and tools for histone PTM analysis described in this study are now publicly available and enable comprehensive investigation into the influence of environmental stress on histone PTMs in nonmodel organisms. Given the functionality and flexibility of histone PTMs, we anticipate that the study of histone PTMs in ecologically relevant contexts will provide ground-breaking insights into comparative physiology and evolution.

Screening and characterization of sex-specific sequences through 2b-RAD sequencing in American shad (Alosa sapidissima)

American shad (Alosa sapidissima), introduced from the United States, has become one of the most expensive farmed fish in the aquatic product market of China. The shad reveals significant sexual dimorphism in growth and behaviors. For the study, five male-specific tags were identified in two-generation breeding populations of Alosa sapidissima and were verified by PCR amplification. Averages of 10,245,091 and 8,685,704 raw and enzyme reads were obtained by high-throughput sequencing of the 2b-RAD library, respectively. 301,022 unique tags were obtained from the sequences of twenty samples with sequencing depths of 0 to 500. Finally, 274,324 special tags and 29,327 SNPs were selected with a sequencing depth of 3 to 500. Eleven preliminary screening male-specific tags and three male heterogametic SNP loci were isolated. After verification by PCR amplification, five male-specific sequences of 27 bp located on chromosome 3 were screened out. Chromosome 3 could be assumed to be the sex chromosome of Alosa sapidissima. Sex-specific markers will provide invaluable and systematic animal germplasm resources to allow for the precise identification of neo-males for the all-female breeding of Alosa sapidissima in commercial aquaculture.

Natural cortisol production is not linked to the sexual fate of European sea bass

In this study, we aimed to investigate the relationship between cortisol and the determination of sexual fate in the commercially important European sea bass (Dicentrarchus labrax). To test our hypothesis, we designed two temperature-based experiments (19 ℃, 21 ℃ and 23 ℃, experiment 1; 16 ℃ and 21 ℃, experiment 2) to assess the effects of these thermal treatments on European sea bass sex determination and differentiation. In the fish from the first experiment, we evaluated whether blood cortisol levels and expression of stress key regulatory genes were different between differentiating (149 to 183 dph) males and females. In the second experiment, we assessed whether cortisol accumulated in scales over time during the labile period for sex determination as well as the neuroanatomical localisation of brain cells expressing brain aromatase (cyp19a1b) and corticotropin-releasing factor (crf) differed between males and females undergoing molecular sex differentiation (117 to 124 dph). None of the gathered results allowed to detect differences between males and females regarding cortisol production and regulatory mechanisms. Altogether, our data provide strong physiological, molecular and histochemical evidence, indicating that in vivo cortisol regulation has no major effects on the sex of European sea bass.

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.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

Polymorphism of Sex Determination Amongst Wild Populations Suggests its Rapid Turnover Within the Nile Tilapia Species

Sex-determining regions have been identified in the Nile tilapia on linkage groups (LG) 1, 20 and 23, depending on the domesticated strains used. Sex determining studies on wild populations of this species are scarce. Previous work on two wild populations, from Lake Volta (Ghana) and from Lake Koka (Ethiopia), found the sex-determining region on LG23. These populations have a Y-specific tandem duplication containing two copies of the Anti-Müllerian Hormone amh gene (named amhY and amhΔY). Here, we performed a whole-genome short-reads analysis using male and female pools on a third wild population from Lake Hora (Ethiopia). We found no association of sex with LG23, and no duplication of the amh gene. Furthermore, we found no evidence of sex linkage on LG1 or on any other LGs. Long read whole genome sequencing of a male from each population confirmed the absence of a duplicated region on LG23 in the Lake Hora male. In contrast, long reads established the structure of the Y haplotype in Koka and Kpandu males and the order of the genes in the duplicated region. Phylogenies constructed on the nuclear and mitochondrial genomes, showed a closer relationship between the two Ethiopian populations compared to the Ghanaian population, implying an absence of the LG23Y sex-determination region in Lake Hora males. Our study supports the hypothesis that the amh region is not the sex-determining region in Hora males. The absence of the Y amh duplication in the Lake Hora population reflects a rapid change in sex determination within Nile tilapia populations. The genetic basis of sex determination in the Lake Hora population remains unknown.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Genome-Wide Identification, Phylogeny, and Expression Profile of the Dmrt (Doublesex and Mab-3 Related Transcription Factor) Gene Family in Channel Catfish (Ictalurus punctatus)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM domain related to sex determination and differentiation, which has been systematically described in various teleost fish, but less in channel catfish (Ictalurus punctatus), an important global aquaculture species in the US and China. In this study, seven Dmrt genes from channel catfish genome were identified and analyzed using bioinformatics methods. Seven IpDmrt genes were distributed unevenly across five chromosomes. Synteny analysis revealed that Dmrt1, Dmrt2a, Dmrt2b, Dmrt3, Dmrt4, and Dmrt5 were relatively conserved in teleost fish. Tissue distribution analysis showed that IpDmrt1, IpDmrt2b, IpDmrt5, and IpDmrt6 exhibited sexually dimorphic expression patterns and, among them, IpDmrt1 and IpDmrt6 had high expression levels in the testes, while IpDmrt2b and IpDmrt5 had more significant expression levels in the ovaries than in other tissues. After 17β-estradiol treatment, IpDmrt2b and IpDmrt5 were significantly up regulated, while the expression of IpDmrt1 and IpDmrt6 was significantly repressed in XY channel catfish ovaries compared with XX channel catfish ovaries. The present study provides a comprehensive insight into the Dmrt gene family of channel catfish. The results suggest that IpDmrt1 and IpDmrt6 may play an important role in testis differentiation/development, while IpDmrt2b and IpDmrt5 are critical in ovary development in this species.

Direct Full-Length RNA Sequencing Reveals an Important Role of Epigenetics During Sexual Reversal in Chinese Soft-Shelled Turtle

Sex dimorphism is a key feature of Chinese soft-shelled turtle (Pelodiscus sinensis). The males (M) have higher econosmic value than females (F) due to wider calipash and faster growth. Exogenous hormones like estradiol and methyltestosterone can induce sexual reversal to form new phenotypes (pseudo-female, PF; pseudo-male, PM) without changing the genotype. The possibility of inducing sexual reversal is particularly important in aquaculture breeding, but the underlying biological mechanisms remain unclear. Here we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in P. sinensis. Nanopore sequencing of the four gender types (M, F, PF, and PM) showed that the distribution of read length and gene expression was more similar between same-sex phenotypes than same-sex genotypes. Compared to turtles with an M phenotype, alternative splicing was more pronounced in F turtles, especially at alternative 3′ splice sites, alternative 5′ splice sites, and alternative first exons. Furthermore, the two RNA methylation modifications m5C and m6A were differentially distributed across gender phenotypes, with the M type having more modification sites in coding sequence regions, but fewer modification sites in 3′UTR regions. Quantitative analysis of enriched m6A RNAs revealed that the N6-methylated levels of Odf2, Pacs2, and Ak1 were significantly higher in M phenotype individuals, while the N6-methylated levels of Ube2o were reduced after sexual reversal from both M and F phenotypes. Taken together, these findings reveal an important role of epigenetics during sexual reversal in Chinese soft-shelled turtles.

Decoding genome recombination and sex reversal

Over the past 440 years since the discovery of the medicinal value of swamp eels, much progress has been made in the study of their biology. The fish is emerging as an important model animal in sexual development, in addition to economic and pharmaceutical implications. Tracing genomic history that shapes speciation of the fish has led to discovery of the whole genome-wide chromosome fission/fusion events. Natural intersex differentiation is a compelling feature for sexual development research. Notably, identification of progenitors of germline stem cells that have bipotential to differentiate into either male or female germline stem cells provides new insight into sex reversal. Here, we review these advances that have propelled the field forwards and present unsolved issues that will guide future investigations to finally elucidate vertebrate sexual development using the new model.

Effects of endocrine disrupting chemicals on gonad development: Mechanistic insights from fish and mammals

Over the past century, evidence has emerged that endocrine disrupting chemicals (EDCs) have an impact on reproductive health. An increased frequency of reproductive disorders has been observed worldwide in both wildlife and humans that is correlated with accidental exposures to EDCs and their increased production. Epidemiological and experimental studies have highlighted the consequences of early exposures and the existence of key windows of sensitivity during development. Such early in life exposures can have an immediate impact on gonadal and reproductive tract development, as well as on long-term reproductive health in both males and females. Traditionally, EDCs were thought to exert their effects by modifying the endocrine pathways controlling reproduction. Advances in knowledge of the mechanisms regulating sex determination, differentiation and gonadal development in fish and rodents have led to a better understanding of the molecular mechanisms underlying the effects of early exposure to EDCs on reproduction. In this manuscript, we review the key developmental stages sensitive to EDCs and the state of knowledge on the mechanisms by which model EDCs affect these processes, based on the roadmap of gonad development specific to fish and mammals.

OUP accepted manuscript

Epigenetic DNA methylation plays an important role in bacteria by influencing gene expression and allowing discrimination between self-DNA and intruders such as phages and plasmids. Restriction–modification (RM) systems use a methyltransferase (MTase) to modify a specific sequence motif, thus protecting host DNA from cleavage by a cognate restriction endonuclease (REase) while leaving invading DNA vulnerable. Other REases occur solitarily and cleave methylated DNA. REases and RM systems are frequently mobile, influencing horizontal gene transfer by altering the compatibility of the host for foreign DNA uptake. However, whether mobile defence systems affect pre-existing host defences remains obscure. Here, we reveal an epigenetic conflict between an RM system (PcaRCI) and a methylation-dependent REase (PcaRCII) in the plant pathogen Pectobacterium carotovorum RC5297. The PcaRCI RM system provides potent protection against unmethylated plasmids and phages, but its methylation motif is targeted by the methylation-dependent PcaRCII. This potentially lethal co-existence is enabled through epigenetic silencing of the PcaRCII-encoding gene via promoter methylation by the PcaRCI MTase. Comparative genome analyses suggest that the PcaRCII-encoding gene was already present and was silenced upon establishment of the PcaRCI system. These findings provide a striking example for selfishness of RM systems and intracellular competition between different defences.

Estrogenic signaling and sociosexual behavior in wild sex-changing bluehead wrasses, Thalassoma bifasciatum

Estrogenic signaling is an important focus in studies of gonadal and brain sexual differentiation in fishes and vertebrates generally. This study examined variation in estrogenic signaling (1) across three sexual phenotypes (female, female-mimic initial phase [IP] male, and terminal phase [TP] male), (2) during socially-controlled female-to-male sex change, and (3) during tidally-driven spawning cycles in the protogynous bluehead wrasse (Thalassoma bifasciatum). We analyzed relative abundances of messenger RNAs (mRNAs) for the brain form of aromatase (cyp19a1b) and the three nuclear estrogen receptors (ER) (ERα, ERβa, and ERβb) by qPCR. Consistent with previous reports, forebrain/midbrain cyp19a1b was highest in females, significantly lower in TP males, and lowest in IP males. By contrast, ERα and ERβb mRNA abundances were highest in TP males and increased during sex change. ERβa mRNA did not vary significantly. Across the tidally-driven spawning cycle, cyp19a1b abundances were higher in females than TP males. Interestingly, cyp19a1b levels were higher in TP males close (~1 h) to the daily spawning period when sexual and aggressive behaviors rise than males far from spawning (~10-12 h). Together with earlier findings, our results suggest alterations in neural estrogen signaling are key regulators of socially-controlled sex change and sexual phenotype differences. Additionally, these patterns suggest TP male-typical sociosexual behaviors may depend on intermediate rather than low estrogenic signaling. We discuss these results and the possibility that an inverted-U shaped relationship between neural estrogen and male-typical behaviors is more common than presently appreciated.

ASER: Animal Sex Reversal Database

Sex reversal, representing extraordinary sexual plasticity during the life cycle, not only triggers reproduction in animals but also affects reproductive and endocrine system-related diseases and cancers in humans. Sex reversal has been broadly reported in animals; however, an integrated resource hub of sex reversal information is still lacking. Here, we constructed a comprehensive database named ASER (Animal Sex Reversal) by integrating sex reversal-related data of 18 species from teleostei to mammalia. We systematically collected 40,018 published papers and mined the sex reversal-associated genes (SRGs), including their regulatory networks, from 1611 core papers. We annotated homologous genes and computed conservation scores for whole genomes across the 18 species. Furthermore, we collected available RNA-seq datasets and investigated the expression dynamics of SRGs during sex reversal or sex determination processes. In addition, we manually annotated 550 in situ hybridization (ISH), fluorescence in situ hybridization (FISH), and immunohistochemistry (IHC) images of SRGs from the literature and described their spatial expression in the gonads. Collectively, ASER provides a unique and integrated resource for researchers to query and reuse organized data to explore the mechanisms and applications of SRGs in animal breeding and human health. The ASER database is publicly available at http://aser.ihb.ac.cn/.

A new experimental model for the investigation of sequential hermaphroditism

The stunning sexual transformation commonly triggered by age, size or social context in some fishes is one of the best examples of phenotypic plasticity thus far described. To date our understanding of this process is dominated by studies on a handful of subtropical and tropical teleosts, often in wild settings. Here we have established the protogynous New Zealand spotty wrasse, Notolabrus celidotus, as a temperate model for the experimental investigation of sex change. Captive fish were induced to change sex using aromatase inhibition or manipulation of social groups. Complete female-to-male transition occurred over 60 days in both cases and time-series sampling was used to quantify changes in hormone production, gene expression and gonadal cellular anatomy. Early-stage decreases in plasma 17β-estradiol (E2) concentrations or gonadal aromatase (cyp19a1a) expression were not detected in spotty wrasse, despite these being commonly associated with the onset of sex change in subtropical and tropical protogynous (female-to-male) hermaphrodites. In contrast, expression of the masculinising factor amh (anti-Müllerian hormone) increased during early sex change, implying a potential role as a proximate trigger for masculinisation. Collectively, these data provide a foundation for the spotty wrasse as a temperate teleost model to study sex change and cell fate in vertebrates.

MicroRNA ssa-mir-196a-4 deceases lgr8 expression in testis development of Chinese tongue sole (Cynoglossus semilaevis)

MicroRNAs (miRNAs) contribute to gonadal development in animals. However, there is little information about miRNA regulation function involved in gonadal development in fish. Our group previously identified sex-related miRNAs of Chinese tongue sole (Cynoglossus semilaevis) during sex determination and differentiation by small RNA sequencing. In the present study, we characterized ssa-mir-196a-4 and its expression in testis and verified its interaction with lgr8. miRNA ssa-mir-196a-4 precursor was predicted to have a typical hairpin structure and highly conserved among various fish species. Fluorescence in situ hybridization (FISH) of ssa-mir-196a-4 in the testis of Chinese tongue sole showed that it is mainly expressed in the cytoplasm of Sertoli cells. We determined that ssa-mir-196a-4 interacted with lgr8 by bioinformatics analysis using miRanda software. According to the dual-luciferase gene reporter assay, lgr8 is a direct target of ssa-mir-196a-4. Overexpression of ssa-mir-196a-4 in the cells of the testis cell line of Chinese tongue sole decreased the expression levels of lgr8 messenger RNA (mRNA) and protein by targeting its coding sequence (CDS) region. These results suggest that ssa-mir-196a-4 acts as a post-transcriptional regulator of lgr8 and plays an important role in developing testes of Chinese tongue sole.

Advances in Reproductive Endocrinology and Neuroendocrine Research Using Catfish Models

Catfishes, belonging to the order siluriformes, represent one of the largest groups of freshwater fishes with more than 4000 species and almost 12% of teleostean population. Due to their worldwide distribution and diversity, catfishes are interesting models for ecologists and evolutionary biologists. Incidentally, catfish emerged as an excellent animal model for aquaculture research because of economic importance, availability, disease resistance, adaptability to artificial spawning, handling, culture, high fecundity, hatchability, hypoxia tolerance and their ability to acclimate to laboratory conditions. Reproductive system in catfish is orchestrated by complex network of nervous, endocrine system and environmental factors during gonadal growth as well as recrudescence. Lot of new information on the molecular mechanism of gonadal development have been obtained over several decades which are evident from significant number of scientific publications pertaining to reproductive biology and neuroendocrine research in catfish. This review aims to synthesize key findings and compile highly relevant aspects on how catfish can offer insight into fundamental mechanisms of all the areas of reproduction and its neuroendocrine regulation, from gametogenesis to spawning including seasonal reproductive cycle. In addition, the state-of-knowledge surrounding gonadal development and neuroendocrine control of gonadal sex differentiation in catfish are comprehensively summarized in comparison with other fish models.

Transcriptomic analysis revealed gene expression profiles during the sex differentiation of Chinese tongue sole (Cynoglossus semilaevis)

Sex differentiation in aquatic fish is important both for theoretical study and practical production, as growth dimorphism frequently appears in different sexes, especially in marine fish. The deciphered genome, identification of the male-determining gene dmrt1 and established genotypic sex screening method make Chinese tongue sole (Cynoglossus semilaevis) an ideal model to study sex differentiation in fish. In this study, comparative gonadal transcriptomic analyses were conducted for genetic females and males at 48, 68, and 108 days post hatching (dph), representing pre-, during- and post-gonadal differentiation stages, although the gonad is not completely differentiated and isolable in 48 and 68 dph individuals, while it is in 108 dph individuals. Altogether, 28 libraries were constructed, and a mean of 46.64 M clean reads was obtained. Differentially expressed gene (DEG) analysis revealed that 179 genes had similar expression patterns in males and females in all three stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the enriched pathways included ubiquitin-mediated proteolysis, lysosomes, and RNA degradation. Moreover, weighted gene coexpression network analyses (WGCNA) identified 14 modules, one of which was closely correlated with female differentiation, exhibiting female-biased expression in all three stages (48, 68, 108 dph). An illustrated core gene interaction network of this module identified 50 genes, most of which are on W chromosomes. Six genes, including two ubiquitin conjugating enzymes, were selected for further investigation, and their female-biased expression was confirmed in even earlier stages, at 10 and 30 dph. These data facilitate our understanding of sex differentiation in fish and provide a genomic rationale for screening candidate genes (preferentially W-linked genes) that could be involved in the female differentiation process.

Protein expression reveals a molecular sexual identity of avian primordial germ cells at pre-gonadal stages

In poultry, in vitro propagated primordial germ cells (PGCs) represent an important tool for the cryopreservation of avian genetic resources. However, several studies have highlighted sexual differences exhibited by PGCs during in vitro propagation, which may compromise their reproductive capacities. To understand this phenomenon, we compared the proteome of pregonadal migratory male (ZZ) and female (ZW) chicken PGCs propagated in vitro by quantitative proteomic analysis using a GeLC-MS/MS strategy. Many proteins were found to be differentially abundant in chicken male and female PGCs indicating their early sexual identity. Many of the proteins more highly expressed in male PGCs were encoded by genes localised to the Z sex chromosome. This suggests that the known lack of dosage compensation of the transcription of Z-linked genes between sexes persists at the protein level in PGCs, and that this may be a key factor of their autonomous sex differentiation. We also found that globally, protein differences do not closely correlate with transcript differences indicating a selective translational mechanism in PGCs. Male and female PGC expressed protein sets were associated with differential biological processes and contained proteins known to be biologically relevant for male and female germ cell development, respectively. We also discovered that female PGCs have a higher capacity to uptake proteins from the cell culture medium than male PGCs. This study presents the first evidence of an early predetermined sex specific cell fate of chicken PGCs and their sexual molecular specificities which will enable the development of more precise sex-specific in vitro culture conditions for the preservation of avian genetic resources.

Updated Genome Assembly of Bighead Carp (Hypophthalmichthys nobilis) and Its Differences Between Male and Female on Genomic, Transcriptomic, and Methylation Level

Cyprinidae is one of the largest family in freshwater fishes, and it is most intensively cultured fish taxon of the world. However, studies about sex determination in this large family is still rear, and one of the reasons is lack of high quality and complete genome. Here, we used nanopore to sequence the genome of a male bighead carp, obtaining contig N50 = 24.25 Mb, which is one of the best assemblies in Cyprinidae. Five males and five females were re-sequenced, and a male-specific region on LG19 was confirmed. We find this region holds many male-specific markers in other Cyprinidae fishes, such as grass carp and silver carp. Transcriptome analyses of hypothalamus and pituitary tissues showed that several sex-specific differentially expressed genes were associated with steroid biosynthesis. The UCH64E gene, located in the male-specific region on LG19, showed higher expression levels in male than female tissues of bighead carp. The methyl-RAD of hypothalamus tissues between males and females indicated that the sexual methylation differences are significant in bighead carp. We also compared the methylation sites recognized using methyl-RAD and nanopore raw reads and found that approximately 73% of the methylation sites identified using methyl-RAD were within nanopore CpG sites.

Epigenetic Regulation and Environmental Sex Determination in Cichlid Fishes

Studying environmental sex determination (ESD) in cichlids provides a phylogenetic and comparative approach to understand the evolution of the underlying mechanisms, their impact on the evolution of the overlying systems, and the neuroethology of life history strategies. Natural selection normally favors parents who invest equally in the development of male and female offspring, but evolution may favor deviations from this 50:50 ratio when environmental conditions produce an advantage for doing so. Many species of cichlids demonstrate ESD in response to water chemistry (temperature, pH, and oxygen concentration). The relative strengths of and the exact interactions between these factors vary between congeners, demonstrating genetic variation in sensitivity. The presence of sizable proportions of the less common sex towards the environmental extremes in most species strongly suggests the presence of some genetic sex-determining loci acting in parallel with the ESD factors. Sex determination and differentiation in these species does not seem to result in the organization of a final and irreversible sexual fate, so much as a life-long ongoing battle between competing male- and female-determining genetic and hormonal networks governed by epigenetic factors. We discuss what is and is not known about the epigenetic mechanism behind the differentiation of both gonads and sex differences in the brain. Beyond the well-studied tilapia species, the 2 best-studied dwarf cichlid systems showing ESD are the South American genus Apistogramma and the West African genus Pelvicachromis. Both species demonstrate male morphs with alternative reproductive tactics. We discuss the further neuroethology opportunities such systems provide to the study of epigenetics of alternative life history strategies and other behavioral variation.

Epigenetic mechanisms in sex determination and in the evolutionary transitions between sexual systems

The hypothesis that epigenetic mechanisms of gene expression regulation have two main roles in vertebrate sex is presented. First, and within a given generation, by contributing to the acquisition and maintenance of (i) the male or female function once during the lifetime in individuals of gonochoristic species; and (ii) the male and female function in the same individual, either at the same time in simultaneous hermaphrodites, or first as one sex and then as the other in sequential hermaphrodites. Second, if environmental conditions change, epigenetic mechanisms may have also a role across generations, by providing the necessary phenotypic plasticity to facilitate the transition: (i) from one sexual system to another, or (ii) from one sex-determining mechanism to another. Furthermore, if the environmental change lasts enough time, epimutations could facilitate assimilation into genetic changes that stabilize the new sexual system or sex-determining mechanism. Examples supporting these assertions are presented, caveats or difficulties and knowledge gaps identified, and possible ways to test this hypothesis suggested.

This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)’.

Might Gene Duplication and Neofunctionalization Contribute to the Sexual Lability Observed in Fish?

Sex determination and differentiation varies widely across vertebrates, but is most dramatically diverse in fishes. Among fishes sex reversal and sex change are observed in 41 teleost families spanning 7 orders. These sex-changing fish perhaps highlight better than any other system that sex determination is not the narrow and fixed construct we once thought, but a plastic trait that is better viewed as a reaction norm. However, while this stunning transformation is increasingly understood, a fundamental question arises, which is why some fish species have retained this inherent plasticity in sexual fate, while others have not? Here, we explore our current understanding of sex change in fish, some of the factors that permit and constrain sex reversal, and posit that gene duplication and neofunctionalization contribute to the sexual lability observed in fish.

Environmental Cues and Mechanisms Underpinning Sex Change in Fish

Fishes are the only vertebrates that undergo sex change during their lifetime, but even within this group, a unique reproductive strategy is displayed by only 1.5% of the teleosts. This lability in alternating sexual fate is the result of the simultaneous suppression and activation of opposing male and female networks. Here, we provide a brief review summarizing recent advances in our understanding of the environmental cues that trigger sex change and their perception, integration, and translation into molecular cascades that convert the sex of an individual. We particularly focus on molecular events underpinning the complex behavioral and morphological transformation involved in sex change, dissecting the main molecular players and regulatory networks that shape the transformation of one sex into the opposite. We show that histological changes and molecular pathways governing gonadal reorganization are better described than the neuroendocrine basis of sex change and that, despite important advances, information is lacking for the majority of hermaphrodite species. We highlight significant gaps in our knowledge of how sex change takes place and suggest future research directions.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

Epigenetics drive the evolution of sex chromosomes in animals and plants

We review how epigenetics affect sex chromosome evolution in animals and plants. In a few species, sex is determined epigenetically through the action of Y-encoded small RNAs. Epigenetics is also responsible for changing the sex of individuals through time, even in species that carry sex chromosomes, and could favour species adaptation through breeding system plasticity. The Y chromosome accumulates repeats that become epigenetically silenced which leads to an epigenetic conflict with the expression of Y genes and could accelerate Y degeneration. Y heterochromatin can be lost through ageing, which activates transposable elements and lowers male longevity. Y chromosome degeneration has led to the evolution of meiotic sex chromosome inactivation in eutherians (placentals) and marsupials, and dosage compensation mechanisms in animals and plants. X-inactivation convergently evolved in eutherians and marsupials via two independently evolved non-coding RNAs. In Drosophila, male X upregulation by the male specific lethal (MSL) complex can spread to neo-X chromosomes through the transposition of transposable elements that carry an MSL-binding motif. We discuss similarities and possible differences between plants and animals and suggest future directions for this dynamic field of research. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Molecular and cellular regulation on sex change in hermaphroditic fish, with a special focus on protandrous black porgy, Acanthopagrus schlegelii

In teleost fish, sex can be determined by genetic factors, environmental factors, or both. Unlike in gonochoristic fish, in which sex is fixed in adults, sex can change in adults of hermaphroditic fish species. Thus, sex is generated during the initial gonadal differentiation stage (primary sex differentiation) and later during sexual fate alternation (secondary sex differentiation) in hermaphroditic fish species. Depending on the species, sex phase alternation can be induced by endogenous cues (such as individual age and body size) or by social cues (such as sex ratio or relative body size within the population). In general, the fluctuation in plasma estradiol-17β (E2) levels is correlated with the sexual fate alternation in hermaphroditic fish. Hormonal treatments can artificially induce sexual phase alternation in sequential hermaphroditic fishes, but in a transient and reversible manner. This is the case for the E2-induced female phase in protandrous black porgy and the methyltestosterone (MT)- or aromatase inhibitor (AI)-induced male phase in protogynous grouper. Recent reviews have focused on the different forms of sex change in fish who undergo sequential sex change, especially in terms of gene expression and the role of hormones. In this review, we use the protandrous black porgy, a nonsocial cue-influenced hermaphroditic species, with digonic gonads (ovarian and testis separated by a connective tissue), as a model to describe our findings and discuss the molecular and cellular regulation of sexual fate determination in hermaphroditic fish.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

Anemonefish, a model for Eco-Evo-Devo

Anemonefish, are a group of about 30 species of damselfish (Pomacentridae) that have long aroused the interest of coral reef fish ecologists. Combining a series of original biological traits and practical features in their breeding that are described in this paper, anemonefish are now emerging as an experimental system of interest for developmental biology, ecology and evolutionary sciences. They are small sized and relatively easy to breed in specific husbandries, unlike the large-sized marine fish used for aquaculture. Because they live in highly structured social groups in sea anemones, anemonefish allow addressing a series of relevant scientific questions such as the social control of growth and sex change, the mechanisms controlling symbiosis, the establishment and variation of complex color patterns, and the regulation of aging. Combined with the use of behavioral experiments, that can be performed in the lab or directly in the wild, as well as functional genetics and genomics, anemonefish provide an attractive experimental system for Eco-Evo-Devo.

Effects of global warming on sex ratios in fishes

In fishes, sex is determined by genetics, the environment or an interaction of both. Temperature is among the most important environmental factors that can affect sex determination. As a consequence, changes in temperature at critical developmental stages can induce biases in primary sex ratios in some species. However, early sex ratios can also be biased by sex-specific tolerances to environmental stresses that may, in some cases, be amplified by changes in water temperature. Sex-specific reactions to environmental stress have been observed at early larval stages before gonad formation starts. It is therefore necessary to distinguish between temperature effects on sex determination, generally acting through the stress axis or epigenetic mechanisms, and temperature effects on sex-specific mortality. Both are likely to affect sex ratios and hence population dynamics. Moreover, in cases where temperature effects on sex determination lead to genotype-phenotype mismatches, long-term effects on population dynamics are possible, for example temperature-induced masculinization potentially leading to the loss of Y chromosomes or feminization to male-biased operational sex ratios in future generations. To date, most studies under controlled conditions conclude that if temperature affects sex ratios, elevated temperatures mostly lead to a male bias. The few studies that have been performed on wild populations seem to confirm this general trend. Recent findings suggest that transgenerational plasticity could mitigate the effects of warming on sex ratios in some populations.

Sex-Dependent RNA Editing and N6-adenosine RNA Methylation Profiling in the Gonads of a Fish, the Olive Flounder (Paralichthys olivaceus)

Adenosine-to-inosine (A-to-I) editing and N6-methyladenosine (m6A) are two of the most abundant RNA modifications. Here, we examined the characteristics of the RNA editing and transcriptome-wide m6A modification profile in the gonads of the olive flounder, Paralichthys olivaceus, an important maricultured fish in Asia. The gonadal differentiation and development of the flounder are controlled by genetic as well as environmental factors, and the epigenetic mechanism may play an important role. In total, 742 RNA editing events were identified, 459 of which caused A to I conversion. Most A-to-I sites were located in 3′UTRs, while 61 were detected in coding regions (CDs). The number of editing sites in the testis was higher than that in the ovary. Transcriptome-wide analyses showed that more than one-half of the transcribed genes presented an m6A modification in the flounder gonads, and approximately 60% of the differentially expressed genes (DEGs) between the testis and ovary appeared to be negatively correlated with m6A methylation enrichment. Further analyses revealed that the mRNA expression of some sex-related genes (e.g., dmrt1 and amh) in the gonads may be regulated by changes in mRNA m6A enrichment. Functional enrichment analysis indicated that the RNA editing and m6A modifications were enriched in several canonical pathways (e.g., Wnt and MAPK signaling pathways) in fish gonads and in some pathways whose roles have not been investigated in relation to fish sex differentiation and gonadal development (e.g., PPAR and RNA degradation pathways). There were 125 genes that were modified by both A-to-I editing and m6A, but the two types of modifications mostly occurred at different sites. Our results suggested that the presence of sex-specific RNA modifications may be involved in the regulation of gonadal development and gametogenesis.