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

Deciphering the genetic basis of sex differentiation in silver-lipped pearl oyster ( Pinctada maxima) based on integrative transcriptomic analysis

The silver-lipped pearl oyster ( Pinctada maxima) is the largest and most commercially valuable pearl-producing oyster, renowned for its ability to generate large, lustrous pearls. This species is a sequential hermaphrodite, with pearl production displaying notable sexual dimorphism. Consequently, understanding the molecular mechanisms governing sex determination and differentiation is crucial for advancing breeding strategies in the pearl oyster industry. To elucidate these mechanisms, this study conducted integrative transcriptomic analyses of P. maxima gonadal tissues using isoform sequencing (Iso-seq) and RNA sequencing (RNA-seq). Comparative analysis of ovarian and testicular tissues identified 2 768 differentially expressed genes (DEGs). Gene co-expression network analysis delineated four key modules, including three sex-specific modules and one shared module. Key genes implicated in sex determination and maintenance were identified, including FOXL2, NANOS1, and β-catenin, important for ovarian maintenance, and DMRT, SOX30, FEM1, and FOXJ1, crucial for testicular maintenance. These genes, widely studied in other taxa, were confirmed as hub genes in the sex-related modules of P. maxima. Interestingly, genes within the shared module were significantly enriched in the spliceosome pathway. Alternative splicing analysis highlighted its extensive role in gonadal tissues, with more pronounced activity observed in the testis compared to the ovary. Nearly half (47.83%, 375) of the identified genes undergoing differential alternative splicing (DASGs) also exhibited differential transcript usage (DTUGs), while only 17% of DTUGs overlapped with DEGs. Genes associated with sex differentiation, such as DMRT, β-catenin, and U2AF2, displayed sex-specific and/or sex-biased isoforms. These findings offer novel insights into the molecular basis of sex differentiation in P. maxima, which could inform the development of targeted breeding strategies aimed at sex control, thereby enhancing pearl quality and yield in aquaculture. This study offers a robust molecular foundation for advancing breeding programs and optimizing production in the pearl oyster industry.

Understanding Genetic Regulation of Sex Differentiation in Hermaphroditic Fish

As a fundamental taxonomic group within vertebrates, fish represent an invaluable resource for investigating the mechanisms underlying sex determination and differentiation owing to their extensive geographical distribution and rich biodiversity. Within this biological cohort, the processes of sex determination and differentiation are intricately governed by both genetic factors and the complex interplay of environmental cues. While variations in external environmental factors, particularly temperature, can exert a modulatory influence on sex differentiation in fish to a limited degree, genetic factors remain the primary determinants of sexual traits. Hermaphroditic fish display three distinct types of sexual transitions: protandry (male to female), protogyny (female-to-male), bidirectional sex change (both directions serially). These fish, characterized by their unique reproductive strategies and sexual plasticity, serve as exemplary natural models for elucidating the mechanisms of sex differentiation and sexual transitions in fish. The present review delves into the histological dynamics during gonadal development across three types of sequential hermaphroditic fish, meticulously delineating the pivotal characteristics at each stage, from the inception of primordial gonads to sexual specialization. Furthermore, it examines the regulatory genes and associated signaling pathways that orchestrate sex determination and differentiation. By systematically synthesizing these research advancements, this paper endeavors to offer a comprehensive and profound insight into the intricate mechanisms governing sex differentiation in sequential hermaphroditic fish.

Genome-wide identification, molecular evolution and gene expression of P450 gene family in Cyrtotrachelus buqueti

BACKGROUND

Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking.

RESULTS

In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds.

CONCLUSIONS

We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Characterization of Two Gonadal Genes, zar1 and wt1b, in Hermaphroditic Fish Asian Seabass (Lates calcarifer)

Zygote arrest-1 (Zar1) and Wilms' tumor 1 (Wt1) play an important role in oogenesis, with the latter also involved in testicular development and gender differentiation. Here, Lczar1 and Lcwt1b were identified in Asian seabass (Lates calcarifer), a hermaphrodite fish, as the valuable model for studying sex differentiation. The cloned cDNA fragments of Lczar1 were 1192 bp, encoding 336 amino acids, and contained a zinc-binding domain, while those of Lcwt1b cDNA were 1521 bp, encoding a peptide of 423 amino acids with a Zn finger domain belonging to Wt1b family. RT-qPCR analysis showed that Lczar1 mRNA was exclusively expressed in the ovary, while Lcwt1b mRNA was majorly expressed in the gonads in a higher amount in the testis than in the ovary. In situ hybridization results showed that Lczar1 mRNA was mainly concentrated in oogonia and oocytes at early stages in the ovary, but were undetectable in the testis. Lcwt1b mRNA was localized not only in gonadal somatic cells (the testis and ovary), but also in female and male germ cells in the early developmental stages, such as those of previtellogenic oocytes, spermatogonia, spermatocytes and spermatids. These results indicated that Lczar1 and Lcwt1b possibly play roles in gonadal development. Therefore, the findings of this study will provide a basis for clarifying the mechanism of Lczar1 and Lcwt1b in regulating germ cell development and the sex reversal of Asian seabass and even other hermaphroditic species.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

Optimization of a non-activating medium for short-term chilled storage of barramundi (Lates calcarifer) testicular spermatozoa

Reliable short-term chilled sperm storage is a critical prerequisite to using advanced reproductive techniques for captive breeding of barramundi (Asian sea bass; Lates calcarifer). Marine Ringer's solution (MRS) is a common non-activating medium (NAM) and has previously been used to store sperm from wild-caught barramundi. However, MRS-stored spermatozoa from captive-bred barramundi were observed to lyse within 30 min incubation. Therefore, this study aimed to optimize the composition of NAM for short-term chilled storage by characterizing and mimicking the biochemical profile of seminal and blood plasma of captive-bred barramundi. To further understand the effect of each component, osmolality was first examined to determine its effect on sperm viability. Thereafter, the effects of NaHCO3, pH, and Na+ and K+ concentrations on sperm motility were investigated. Optimization of the NAM formula was achieved through iterative adaptions. The increase in NAM osmolality from 260 to 400 mOsm/kg led to a significant improvement in sperm viability. Moreover, using HEPES instead of NaHCO3 as buffering agent significantly enhanced sperm motility and velocity. As a result, sperm samples diluted with optimized NAM (185 mM NaCl, 5.1 mM KCl, 1.6 mM CaCl2·2H2O, 1.1 mM MgSO4·7H2O, 10.0 mM HEPES, 5.6 mM D+ glucose, 400 mOsm/kg, pH 7.4) and stored at 4 °C showed no significant loss in total motility for up to 48 h and retained progressive motility for up to 72 h. The optimized NAM developed in this study significantly extended the functional lifespan of spermatozoa during chilled storage, permitting the ongoing development of advanced reproductive technologies for barramundi.

The Price of Human Evolution: Cancer-Testis Antigens, the Decline in Male Fertility and the Increase in Cancer

The increasing frequency of general and particularly male cancer coupled with the reduction in male fertility seen worldwide motivated us to seek a potential evolutionary link between these two phenomena, concerning the reproductive transcriptional modules observed in cancer and the expression of cancer-testis antigens (CTA). The phylostratigraphy analysis of the human genome allowed us to link the early evolutionary origin of cancer via the reproductive life cycles of the unicellulars and early multicellulars, potentially driving soma-germ transition, female meiosis, and the parthenogenesis of polyploid giant cancer cells (PGCCs), with the expansion of the CTA multi-families, very late during their evolution. CTA adaptation was aided by retrovirus domestication in the unstable genomes of mammals, for protecting male fertility in stress conditions, particularly that of humans, as compensation for the energy consumption of a large complex brain which also exploited retrotransposition. We found that the early and late evolutionary branches of human cancer are united by the immunity-proto-placental network, which evolved in the Cambrian and shares stress regulators with the finely-tuned sex determination system. We further propose that social stress and endocrine disruption caused by environmental pollution with organic materials, which alter sex determination in male foetuses and further spermatogenesis in adults, bias the development of PGCC-parthenogenetic cancer by default.

Development of Ovaries and Sex Change in Fish: Bringing Potential into Action

BACKGROUND

Encompassing about half of the 60,000 species of vertebrates, fish display the greatest diversity of sex determination mechanisms among metazoans. As such that phylum offers a unique playground to study the impressive variety of gonadal morphogenetic strategies, ranging from gonochorism, with either genetic or environmental sex determination, to unisexuality, with either simultaneous or consecutive hermaphroditism.

SUMMARY

From the two main types of gonads, the ovaries embrace the important role to produce the larger and non-motile gametes, which is the basis for the development of a future organism. The production of the egg cells is complex and involves the formation of follicular cells, which are necessary for the maturation of the oocytes and the production of feminine hormones. In this vein, our review focuses on the development of ovaries in fish with special emphasis on the germ cells, including those that transition from one sex to the other as part of their life cycle and those that are capable of transitioning to the opposite sex depending on environmental cues.

KEY MESSAGES

Clearly, establishing an individual as either a female or a male is not accomplished by the sole development of two types of gonads. In most cases, that dichotomy, be it final or transient, is accompanied by coordinated transformations across the entire organism, leading to changes in the physiological sex as a whole. These coordinated transformations require both molecular and neuroendocrine networks, but also anatomical and behavioural adjustments. Remarkably, fish managed to tame the ins and outs of sex reversal mechanisms to take the most advantages of changing sex as adaptive strategies in some situations.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

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

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

Comparative gonad transcriptome analysis in cobia (Rachycentron canadum)

Background: Cobia (Rachycentron canadum) is a species of fish with high commercial potential particularly due to fast growth rates. The evidence of sexual size dimorphism favoring females indicate potential benefits in having a monosex culture. However, the involvement of genetic factors responsible for sexual development and gonadal maintenance that produces phenotypic sex in cobia is largely unknown. Methods: In the present study, we performed transcriptome sequencing of cobia to identify sex-biased significantly differentially expressed genes (DEGs) in testes and ovaries. The reliability of the gonad transcriptome data was validated by qPCR analysis of eight selected significantly differential expressed sex-related candidate genes. Results: This comparative gonad transcriptomic analysis revealed that 7,120 and 4,628 DEGs are up-regulated in testes or ovaries, respectively. Further functional annotation analyses identified 76 important candidate genes involved in sex determination cascades or sex differentiation, including 42 known testis-biased DEGs (dmrt1, amh and sox9 etc.), and 34 known ovary-biased DEGs (foxl2, sox3 and cyp19a etc.). Moreover, eleven significantly enriched pathways functionally related to sex determination and sex differentiation were identified, including Wnt signaling pathway, oocyte meiosis, the TGF-beta signaling pathway and MAPK signaling pathway. Conclusion: This work represents the first comparative gonad transcriptome study in cobia. The putative sex-associated DEGs and pathways provide an important molecular basis for further investigation of cobia's sex determination, gonadal development as well as potential control breeding of monosex female populations for a possible aquaculture setting.

Sexually Dimorphic Gene Expression in X and Y Sperms Instructs Sexual Dimorphism of Embryonic Genome Activation in Yellow Catfish (Pelteobagrus fulvidraco)

Paternal factors play an important role in embryonic morphogenesis and contribute to sexual dimorphism in development. To assess the effect of paternal DNA on sexual dimorphism of embryonic genome activation, we compared X and Y sperm and different sexes of embryos before sex determination. Through transcriptome sequencing (RNA-seq) and whole-genome bisulfite sequencing (WGBS) of X and Y sperm, we found a big proportion of upregulated genes in Y sperm, supported by the observation that genome-wide DNA methylation level is slightly lower than in X sperm. Cytokine-cytokine receptor interaction, TGF-beta, and toll-like receptor pathways play important roles in spermatogenesis. Through whole-genome re-sequencing (WGRS) of parental fish and RNA-seq of five early embryonic stages, we found the low-blastocyst time point is a key to maternal transcriptome degradation and zygotic genome activation. Generally, sexual differences emerged from the bud stage. Moreover, through integrated analysis of paternal SNPs and gene expression, we evaluated the influence of paternal inheritance on sexual dimorphism of genome activation. Besides, we screened out gata6 and ddx5 as potential instructors for early sex determination and gonad development in yellow catfish. This work is meaningful for revealing the molecular mechanisms of sex determination and sexual dimorphism of fish species.

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

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

Expression and regulation of 42Sp50 in spotted scat (Scatophagus argus)

42Sp50 is an isoform of the eukaryotic translation elongation factor 1 A (eEF1A) and is vital for fish ovarian development. Spotted scat (Scatophagus argus) is a popular marine cultured fish species in Southern Asia and China, and its artificial reproduction is complicated, with a relatively low success ratio in practice. In this study, the 42Sp50 gene was cloned from spotted scat. Tissue distribution analysis showed that 42Sp50 was mainly expressed in the ovary. qRT-PCR showed that 42Sp50 expression levels gradually decreased insignificantly in the ovaries from phase II to IV. Western blot analysis showed that 42Sp50 was highly expressed in the ovary, while it was almost undetectable in the testis. Immunohistochemistry analysis stained 42Sp50 mainly in the cytoplasm of the previtellogenic oocytes in ovaries of normal XX-female and sex-reversed XY-female. Aside from fish and amphibians, 42Sp50 was also identified in some reptile species using genomic database searching. Analyses of the transcriptome data from four different fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), and Hong Kong catfish (Clarias fuscus)) revealed ovaries biased expression of 42Sp50 in all, similar to spotted scat. While the neighbor genes of 42Sp50 did not show ovary biased expression in the fish species analyzed. Bisulfite Sequencing PCR (BSP) results showed that the DNA methylation level of 42Sp50 promoter was low in ovaries, testes, and muscles. The luciferase reporter assay demonstrated that Dmrt4 activated 42Sp50 expression in the presence of Sf1 or Foxh1. These results suggest that 42Sp50 may be involved in regulating the early phase oocytes development of spotted scat.

First account of a transient intersex in spotted scat, Scatophagus argus: a marine gonochoristic fish

This study presents the first incidence of intersex associated with testis-ova in spotted scat (Scatophagus argus) reared in a controlled environment. The testis-ova is associated with the abnormal occurrence of primary oocytes (POs) in some male testis and is referred to as ectopic primary oocytes (Ecto-PO), whiles individuals with Ecto-PO are called "Ecto-PO gonad/individuals." We investigated gonads of 129 male spotted scat aged 4-12 and 18 months after hatch (mah) by histological studies for the presence of female sexual characteristics. A total of 20 out of 88 gonads representing 22.7% of male fish aged 6-12, or 15.5% of all male fish sampled, were found to have visible Ecto-PO. At least, the Ecto-PO had an average of 7 oocytes per gonadal section, indicating high severity. The Ecto-PO appears after sex differentiation and degenerates during sexual maturation. The Ecto-PO did not significantly influence the expression pattern of male and female sex-related genes performed using qPCR. Immunofluorescence of 42sp50 specifically stained the Ecto-PO without influence from the surrounding testicular tissues. In addition, temperature did not correlate with the proliferation of the Ecto-PO, but rather gonad developmental strategy. The results show that the naturally occurring Ecto-PO might not be detrimental to testis development and could be considered a frequent-high-level incidence of natural aberration. This study highlights the intricacy of fish sex differentiation and provides a new research chapter to ascertain the mystery behind the occurrence of Ecto-PO.

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.

Epigenetics underpins phenotypic plasticity of protandrous sex change in fish

Phenotypic plasticity is an important driver of species resilience. Often mediated by epigenetic changes, phenotypic plasticity enables individual genotypes to express variable phenotypes in response to environmental change. Barramundi (Lates calcarifer) are a protandrous (male-first) sequential hermaphrodite that exhibits plasticity in length-at-sex change between geographic regions. This plasticity is likely to be mediated by changes in DNA methylation (DNAm), a well-studied epigenetic modification. To investigate the relationships between length, sex, and DNAm in a sequential hermaphrodite, here, we compare DNAm in four conserved vertebrate sex-determining genes in male and female barramundi of differing lengths from three geographic regions of northern Australia. Barramundi first mature as male and later sex change to female upon the attainment of a larger body size; however, a general pattern of increasing female-specific DNAm markers with increasing length was not observed. Significant differences in DNAm between males and females of similar lengths suggest that female-specific DNAm arises rapidly during sex change, rather than gradually with fish growth. The findings also reveal that region-specific differences in length-at-sex change are accompanied by differences in DNAm and are consistent with variability in remotely sensed sea temperature and salinity. Together, these findings provide the first in situ evidence for epigenetically and environmentally mediated sex change in a protandrous hermaphrodite and offer significant insight into the molecular and ecological processes governing the marked and unique plasticity of sex in fish.

CHARGE syndrome-associated proteins FAM172A and CHD7 influence male sex determination and differentiation through transcriptional and alternative splicing mechanisms

To gain further insight into chromatin‐mediated regulation of mammalian sex determination, we analyzed the role of the CHARGE syndrome‐associated proteins FAM172A and CHD7. This study is based on our prior discoveries that a subset of corresponding mutant mice display complete male‐to‐female sex reversal, and that both of these proteins regulate co‐transcriptional alternative splicing in neural crest cells. Here, we report that FAM172A and CHD7 are present in the developing gonads when sex determination normally occurs in mice. The interactome of FAM172A in pre‐Sertoli cells again suggests a role at the chromatin‐spliceosome interface, like in neural crest cells. Accordingly, analysis of Fam172a‐mutant pre‐Sertoli cells revealed transcriptional and splicing dysregulation of hundreds of genes. Many of these genes are similarly affected in Chd7‐mutant pre‐Sertoli cells, including several known key regulators of sex determination and subsequent formation of testis cords. Among them, we notably identified Sry as a direct transcriptional target and WNT pathway‐associated Lef1 and Tcf7l2 as direct splicing targets. The identified molecular defects are also associated with the abnormal morphology of seminiferous tubules in mutant postnatal testes. Altogether, our results thus identify FAM172A and CHD7 as new players in the regulation of male sex determination and differentiation in mice, and further highlight the importance of chromatin‐mediated regulatory mechanisms in these processes.

Epigenetic Measurement of Key Vertebrate Population Biology Parameters

The age, sex, and sexual maturity of individual animals are key parameters in assessing wild populations and informing conservation management strategies. These parameters represent the reproductive potential of a population and can indicate recovery rates or vulnerabilities. Natural populations of wild animals are difficult to study; logistically, economically, and due to the impacts of invasive biomonitoring. Genetic and epigenetic analyses offer a low impact, low cost, and information-rich alternative. As epigenetic mechanisms are intrinsically linked with both biological aging and reproductive processes, DNA methylation can be used as a suitable biomarker for population biology study. This review assesses published research utilizing DNA methylation analysis in relation to three key population parameters: age, sex, and sexual maturity. We review studies on wild vertebrates that investigate epigenetic age relationships, with successful age estimation assays designed for mammals, birds, and fish. For both determination of sex and identification of sexual maturity, very little has been explored regarding DNA methylation-based assays. Related research, however, confirms the links between DNA methylation and these processes. Future development of age estimation assays for underrepresented and key conservation taxa is suggested, as is the experimental development and design of DNA methylation-based assays for both sex and sexual maturity identification, further expanding the genomics toolkit for population biology studies.

Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits

DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.

Role of Alternative Splicing in Sex Determination in Vertebrates

During the process of sex determination, a germ-cell-containing undifferentiated gonad is converted into either a male or a female reproductive organ. Both the composition of sex chromosomes and the environment determine sex in vertebrates. It is assumed that transcription level regulation drives this cascade of mechanisms; however, transcription factors can alter gene expression beyond transcription initiation by controlling pre-mRNA splicing and thereby mRNA isoform production. Using the key time window in sex determination and gonad development in mice, it has been reported that new non-transcriptional events, such as alternative splicing, could play a key role in sex determination in mammals. We know the role of key regulatory factors, like WT1(+/-KTS) or FGFR2(b/c) in pre-mRNA splicing and sex determination, indicating that important steps in the vertebrate sex determination process probably operate at a post-transcriptional level. Here, we discuss the role of pre-mRNA splicing regulators in sex determination in vertebrates, focusing on the new RNA-seq data reported from mice fetal gonadal transcriptome.

Potential role of DNA methylation of cyp19a1a promoter during sex change in protogynous orange-spotted grouper, Epinephelus coioides

Estrogen has a pivotal role in early female differentiation and further ovarian development. Aromatase (Cyp19a) is responsible for the conversion of androgens to estrogens in vertebrates. In teleosts, cyp19a1a and it paralog cyp19a1b are mainly expressed in the ovary and hypothalamus, respectively. Decreased plasma estrogen levels and lower cyp19a1a expression are associated with the initiation of female-to-male sex change in protogynous grouper. However, an 17α-methyltestosterone (MT)-induced the sex change from a female to a precocious male is a transient phase, and a reversible sex change (induced male-to-female) occurs after chemical withdrawal. Thus, we used this characteristic to study the epigenetic modification of cyp19a1a promoter in orange-spotted grouper. CpG-rich region with a CpG island is located on the putative regulatory region of distal cyp19a1a promoter. Our results showed that cyp19a1a promoter exhibited tissue-specific methylation status. Low methylation levels of distal cyp19a1a promoter and hypomethylated (0-40%) clones of cyp19a1a promoter region were widely observed in the ovary but not shown in testis and other tissues. In femaleness, higher numbers of hypomethylated clones of cyp19a1a promoter region were observed in the vitellogenic oocyte stage compared to the primary oocyte stage. Furthermore, decreased numbers of hypomethylated clones of cyp19a1a promoter region were associated with the maleness during the female-to-male sex change. DNA methylation inhibitor (5-aza-2'-deoxycytidine) delayed the spermatogenesis process (according to germ cell stage and numbers: by decrease of sperm and increase of spermatocytes) but did not influence the reversed sex change in MT-induced bi-directional sex change. These results suggest that epigenetic modification of cyp19a1a promoter may play an important role during the sex change in orange-spotted grouper.

High temperature induced masculinization of zebrafish by down-regulation of sox9b and esr1 via DNA methylation

Elevated temperature could influence the sex differentiation by altering the expression of sex-related genes in fish. However, the underlying mechanisms by which the gene expression is altered remain poorly understood. Here, we aimed to explore the role of DNA methylation in sex differentiation of zebrafish (Danio rerio) in response to elevated temperature. The results showed that high temperature (33°C) exposure of fish from 20 to 30 days post fertilization (dpf), compared to normal temperature (28°C), resulted in male-biased sex ratio and decreased expression of female-related genes including cyp19a1a, sox9b and esr1. Meanwhile, the expressions of DNA methyltransferases dnmt3a1 and dnmt3a2, and the DNA methylation levels in sox9b and esr1 promoter were significantly increased by high temperature, strongly implying that DNA methylation is involved in high temperature-induced masculinization of zebrafish. Co-treatment with 5-aza-2'-deoxycytidine (a DNA methylation inhibitor) attenuated the high temperature-induced masculinizing effect, recovered the expression of esr1 and sox9b, suppressed the transcription of dnmt3a1 and dnmt3a2, and decreased the methylation of esr1 and sox9b promoter, further confirming that DNA methylation plays an important role in high temperature-induced masculinization of zebrafish. Furthermore, the methylation of sox9b promoter decreased the enrichment of transcription factor CREB (cAMP-responsive element binding proteins). Overall, these findings suggest that high temperature induce masculinization of zebrafish by down-regulation of female-related genes via DNA methylation, providing a new insight in understanding the epigenetic mechanism of thermal-mediated sex differentiation in fish.

Differences in gene expression and variable splicing events of ovaries between large and small litter size in Chinese Xiang pigs

BACKGROUND

Although lots of quantitative trait loci (QTLs) and genes present roles in litter size of some breeds, the information might not make it clear for the huge diversity of reproductive capability in pig breeds. To elucidate the inherent mechanisms of heterogeneity of reproductive capability in litter size of Xiang pig, we performed transcriptome analysis for the expression profile in ovaries using RNA-seq method.

RESULTS

We identified 1,419 up-regulated and 1,376 down-regulated genes in Xiang pigs with large litter size. Among them, 1,010 differentially expressed genes (DEGs) were differently spliced between two groups with large or small litter sizes. Based on GO and KEGG analysis, numerous members of genes were gathered in ovarian steroidogenesis, steroid biosynthesis, oocyte maturation and reproduction processes.

CONCLUSIONS

Combined with gene biological function, twelve genes were found out that might be related with the reproductive capability of Xiang pig, of which, eleven genes were recognized as hub genes. These genes may play a role in promoting litter size by elevating steroid and peptide hormones supply through the ovary and facilitating the processes of ovulation and in vivo fertilization.

Scale Drop Disease Virus (SDDV) and Lates calcarifer Herpes Virus (LCHV) Coinfection Downregulate Immune-Relevant Pathways and Cause Splenic and Kidney Necrosis in Barramundi Under Commercial Farming Conditions

Marine farming of barramundi (Lates calcarifer) in Southeast Asia is currently severely affected by viral diseases. To better understand the biological implications and gene expression response of barramundi in commercial farming conditions during a disease outbreak, the presence of pathogens, comparative RNAseq, and histopathology targeting multiple organs of clinically “sick” and “healthy” juveniles were investigated. Coinfection of scale drop disease virus (SDDV) and L. calcarifer herpes virus (LCHV) were detected in all sampled fish, with higher SDDV viral loads in sick than in healthy fish. Histopathology showed that livers in sick fish often had moderate to severe abnormal fat accumulation (hepatic lipidosis), whereas the predominant pathology in the kidneys shows moderate to severe inflammation and glomerular necrosis. The spleen was the most severely affected organ, with sick fish presenting severe multifocal and coalescing necrosis. Principal component analysis (PC1 and PC2) explained 70.3% of the observed variance and strongly associated the above histopathological findings with SDDV loads and with the sick phenotypes, supporting a primary diagnosis of the fish being impacted by scale drop disease (SDD). Extracted RNA from kidney and spleen of the sick fish were also severely degraded likely due to severe inflammation and tissue necrosis, indicating failure of these organs in advanced stages of SDD. RNAseq of sick vs. healthy barramundi identified 2,810 and 556 differentially expressed genes (DEGs) in the liver and muscle, respectively. Eleven significantly enriched pathways (e.g., phagosome, cytokine-cytokine-receptor interaction, ECM-receptor interaction, neuroactive ligand-receptor interaction, calcium signaling, MAPK, CAMs, etc.) and gene families (e.g., tool-like receptor, TNF, lectin, complement, interleukin, chemokine, MHC, B and T cells, CD molecules, etc.) relevant to homeostasis and innate and adaptive immunity were mostly downregulated in sick fish. These DEGs and pathways, also previously identified in L. calcarifer as general immune responses to other pathogens and environmental stressors, suggest a failure of the clinically sick fish to cope and overcome the systemic inflammatory responses and tissue degeneration caused by SDD.

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.

The divergence of alternative splicing between ohnologs in teleost fishes

BACKGROUND

Gene duplication and alternative splicing (AS) are two distinct mechanisms generating new materials for genetic innovations. The evolutionary link between gene duplication and AS is still controversial, due to utilizing duplicates from inconsistent ages of duplication events in earlier studies. With the aid of RNA-seq data, we explored evolutionary scenario of AS divergence between duplicates with ohnologs that resulted from the teleost genome duplication event in zebrafish, medaka, and stickleback.

RESULTS

Ohnologs in zebrafish have fewer AS forms compared to their singleton orthologs, supporting the function-sharing model of AS divergence between duplicates. Ohnologs in stickleback have more AS forms compared to their singleton orthologs, which supports the accelerated model of AS divergence between duplicates. The evolution of AS in ohnologs in medaka supports a combined scenario of the function-sharing and the accelerated model of AS divergence between duplicates. We also found a small number of ohnolog pairs in each of the three teleosts showed significantly asymmetric AS divergence. For example, the well-known ovary-factor gene cyp19a1a has no AS form but its ohnolog cyp19a1b has multiple AS forms in medaka, suggesting that functional divergence between duplicates might have result from AS divergence.

CONCLUSIONS

We found that a combined scenario of function-sharing and accelerated models for AS evolution in ohnologs in teleosts and rule out the independent model that assumes a lack of correlation between gene duplication and AS. Our study thus provided insights into the link between gene duplication and AS in general and ohnolog divergence in teleosts from AS perspective in particular.

Transgenerational epigenetic sex determination: Environment experienced by female fish affects offspring sex ratio

Sex determination is a complex process that can be influenced by environment in various taxa. Disturbed environments can affect population sex ratios and thus threaten their viability. Emerging evidences support a role of epigenetic mechanisms, notably DNA methylation, in environmental sex determination (ESD). In this work, using zebrafish as model and a transgenerational experiment comprising 4 successive generations, we report a strength link between the promotor methylation level of three genes in female gonads and population sex ratio. One generation of zebrafish was exposed throughout its lifetime to cadmium (Cd), a non-essential metal, at an environmentally relevant concentration. The subsequent generations were not exposed. At the first and the third generation a subset of individuals was exposed to an elevated temperature, a well-known masculinizing factor in zebrafish. While heat was associated to an increase in the methylation level of cyp19a1a gene and population masculinization, foxl2a/dmrt1 methylation levels appeared to be influenced by Cd and fish density leading to offspring feminization. Ancestral Cd exposure indeed led to a progressive feminization of the population over generations and affected the sex plastic response of zebrafish in response to heat. The effect of Cd on the methylation level of foxl2a was observed until the third generation, supporting potential transgenerational inheritance. Our results support (i) a key role of cyp19a1a methylation in SD in zebrafish in response to environmental cues and (ii) the fact that the environment experienced by parents, namely mothers in the present case, can affect their offspring sex ratio via environment-induced DNA methylation changes in gonads.

A Transient Hermaphroditic Stage in Early Male Gonadal Development in Little Yellow Croaker, Larimichthys polyactis

Animal taxa show remarkable variability in sexual reproduction, where separate sexes, or gonochorism, is thought to have evolved from hermaphroditism for most cases. Hermaphroditism accounts for 5% in animals, and sequential hermaphroditism has been found in teleost. In this study, we characterized a novel form of the transient hermaphroditic stage in little yellow croaker (Larimichthys polyactis) during early gonadal development. The ovary and testis were indistinguishable from 7 to 40 days post-hatching (dph). Morphological and histological examinations revealed an intersex stage of male gonads between 43 and 80 dph, which consist of germ cells, somatic cells, efferent duct, and early primary oocytes (EPOs). These EPOs in testis degenerate completely by 90 dph through apoptosis yet can be rescued by exogenous 17-β-estradiol. Male germ cells enter the mitotic flourishing stage before meiosis is initiated at 180 dph, and they undergo normal spermatogenesis to produce functional sperms. This transient hermaphroditic stage is male-specific, and the ovary development appears to be normal in females. This developmental pattern is not found in the sister species Larimichthys crocea or any other closely related species. Further examinations of serum hormone levels indicate that the absence of 11-ketotestosterone and elevated levels of 17-β-estradiol delineate the male intersex gonad stage, providing mechanistic insights on this unique phenomenon. Our research is the first report on male-specific transient hermaphroditism and will advance the current understanding of fish reproductive biology. This unique gonadal development pattern can serve as a useful model for studying the evolutionary relationship between hermaphroditism and gonochorism, as well as teleost sex determination and differentiation strategies.

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.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Immune genes, IL1β and Casp9, show sexual dimorphic methylation patterns in zebrafish gonads

There is crosstalk between the immune and reproductive systems in which sexual dimorphism is a common pattern in vertebrates. In recent years, epigenetics has emerged as a way to study the molecular mechanisms involved in gonadal development, those responsible for integrating environmental information that contribute to assigning a specific sexual phenotype (either an ovary or a testis). The knowledge of epigenetic mechanisms in certain molecular processes allows the development of epigenetic markers. In fish gonads, the existence of reproduction-immune system interactions is known, although the epigenetic mechanisms involved are far from clear. Here, we used the zebrafish (Danio rerio) as a model to study the DNA methylation patterns in gonads of two well-known innate immune genes: IL1β and Casp9. DNA methylation levels were studied by a candidate gene approach at single nucleotide resolution and gene expression analyses were also carried out. Results showed that there was clear sexual dimorphism in the DNA methylation levels of the two immune genes studied, being significantly higher in the testes when compared to the ovaries. In summary, and although further research is needed, this paper presents sexual dimorphic methylation patterns of two immune-related genes, thus sex-biased differences in methylation profiles should considered when analyzing immune responses in fish. Data showed here can help to develop epimarkers with forthcoming applications in livestock and fish farming production, for example, in immune fish diseases or sexual control programs as epigenetic molecular tools to predict environmental pressure in the gonads.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.)

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

The Model of the Conserved Epigenetic Regulation of Sex

Epigenetics integrates genomic and environmental information to produce a given phenotype. Here, the model of Conserved Epigenetic Regulation of Sex (CERS) is discussed. This model is based on our knowledge on genes involved in sexual development and on epigenetic regulation of gene expression activation and silencing. This model was recently postulated to be applied to the sexual development of fish, and it states that epigenetic and gene expression patterns are more associated with the development of a particular gonadal phenotype, e.g., testis differentiation, rather than with the intrinsic or extrinsic causes that lead to the development of this phenotype. This requires the existence of genes with different epigenetic modifications, for example, changes in DNA methylation levels associated with the development of a particular sex. Focusing on DNA methylation, the identification of CpGs, the methylation of which is linked to sex, constitutes the basis for the identification of Essential Epigenetic Marks (EEM). EEMs are defined as the number and identity of informative epigenetic marks that are strictly necessary, albeit perhaps not sufficient, to bring about a specific, measurable, phenotype of interest. Here, we provide a summary of the genes where DNA methylation has been investigated so far, focusing on fish. We found that cyp19a1a and dmrt1, two key genes for ovary and testis development, respectively, consistently show an inverse relationship between their DNA methylation and expression levels, thus following CERS predictions. However, in foxl2a, a pro-female gene, and amh, a pro-male gene, such relationship is not clear. The available data of other genes related to sexual development such as sox9, gsdf, and amhr2 are also discussed. Next, we discuss the use of CERS to make testable predictions of how sex is epigenetically regulated and to better understand sexual development, as well as the use of EEMs as tools for the diagnosis and prognosis of sex. We argue that CERS can aid in focusing research on the epigenetic regulation of sexual development not only in fish but also in vertebrates in general, particularly in reptiles with temperature sex-determination, and can be the basis for possible practical applications including sex control in aquaculture and also in conservation biology.

The Genetics and Epigenetics of Sex Change in Fish

Fish show extraordinary sexual plasticity, changing sex naturally as part of their life cycle or reversing sex because of environmental stressors. This plasticity shows that sexual fate is not an irreversible process but the result of an ongoing tug-of-war for supremacy between male and female signaling networks. The behavioral, gonadal, and morphological changes involved in this process are well described, yet the molecular events that underpin those changes remain poorly understood. Epigenetic modifications emerge as a critical link between environmental stimuli, the onset of sex change, and subsequent maintenance of sexual phenotype. Here we synthesize current knowledge of sex change, focusing on the genetic and epigenetic processes that are likely involved in the initiation and regulation of sex change. We anticipate that better understanding of sex change in fish will shed new light on sex determination and development in vertebrates and on how environmental perturbations affect sexual fate.