Dmrt1 genes at the crossroads: a widespread and central class of sexual development factors in fish

Transcriptome analysis of four types of gonadal tissues in largemouth bass (Micropterus salmoides) to reveal its sex-related genes

The sex determination system of largemouth bass (Micropterus salmoides, LMB) is XX/XY; however, the underlying molecular mechanisms involved in early sex differentiation, gonadal development, and exogenous hormone-induced sex reversal remain unknown. In this study, LMB at 15 days post-hatching (dph) were fed diets containing 20 mg/kg of 17α-methyltestosterone (17α-MT) or 30 mg/kg of 17β-estradiol (17β-E2) for 60 days, respectively. Serum steroid levels, histological observations of the gonads, and identification of sex-specific markers were employed to screen the gonads of 60-day-old normal female fish (XX-F), normal male fish (XY-M), 17β-E2 induced pseudo-female fish (XY-F), and 17α-MT-induced pseudo-male fish (XX-M) for transcriptome sequencing in order to uncover genes and pathway involved in the process of sexual reversal. The results from histology and serum sex steroid hormone analysis showed that both 17α-MT and 17β-E2 were capable of inducing sex reversal of LMB at 15 dph. Transcriptome results revealed a total of 2,753 genes exhibiting differential expression, and the expression pattern of these genes in the gonads of XX-M or XY-F resembled that of normal females or males. The male sex-biased genes that are upregulated in XX-M and downregulated in XY-F are referred to as key genes for male reversal, while the female sex-biased genes that are upregulated in XY-F and downregulated in XX-M are referred to as key genes for female reversal. Finally, 12 differentially expressed genes (DEGs) related to male sex reversal were screened, including star2, cyp17a, cyp11b1, dmrt1, amh, sox9a, katnal1, spata4, spata6l, spata7, spata18 and foxl3. 2 DEGs (foxl2a and cyp19a1b) were found to be associated with female sex reversal. The changes in these genes collectively influence the direction of sex differentiation of LMB. Among them, star2, dmrt1 and cyp19a1b with significantly altered expression levels may play potentially crucial role in the process of gender reversal. The expression patterns of 21 randomly selected genes were verified using qRT-PCR which confirmed the reliability and accuracy of the RNA-seq results. These findings not only enhance our understanding of the molecular basis underlying sex reversal but also provide crucial data support for future breeding research on unisexual LMB.

Medaka (Oryzias latipes) Dmrt3a Is Involved in Male Fertility

Research across various species has demonstrated that the doublesex and mab-3-related transcription factor 3 (dmrt3) plays pivotal roles in testis development. However, the precise molecular mechanisms of dmrt3 remain unclear. In this study, we investigated the role of dmrt3 (dmrt3a) in testis development using the model organism medaka (Oryzias latipes). SqRT-PCR and ISH analyses revealed that dmrt3a is predominantly expressed in the testis, especially in the spermatid and spermatozoon. Using CRISPR/Cas9, we generated two dmrt3a homozygous mutants (-8 bp and -11 bp), which exhibited significantly reduced fertilization rates and embryo production. Additionally, the number of germ cells and sperm motility were markedly decreased in the dmrt3a mutants, manifesting as the symptoms of asthenozoospermia and oligozoospermia. Interestingly, RNA-Seq analysis showed that the deficiency of dmrt3a could lead to a significant downregulation of numerous genes related to gonadal development and severe disruptions in mitochondrial function. These results suggested that dmrt3a is essential for spermatogenesis and spermatozoa energy production. This paper provides new insights and perspectives for further exploring the molecular mechanisms underlying spermatogenesis and addressing male reproductive issues.

A computational analysis of the molecular mechanisms underlying the effects of ibuprofen and dibutyl phthalate on gene expression in fish

The impact of emerging pollutants such as ibuprofen and dibutyl phthalate on aquatic species is a growing concern and the need for proper assessment and evaluation of these toxicants is imperative. The objective of this study was to examine the toxicogenomic impacts of ibuprofen and dibutyl phthalate on Clarias gariepinus, a widely distributed African catfish species. Results showed that exposure to the test compounds caused significant changes in gene expression, including upregulation of growth hormone, interleukin, melatonin receptors, 17β-Hydroxysteroid Dehydrogenase, heat shock protein, doublesex, and mab-3 related transcription factor. On the other hand, expression of forkhead Box Protein L2 and cytochrome P450 was downregulated, revealing a potential to induce female to male sex reversal. The binding affinities and hydrophobic interactions of the test compounds with the reference genes were also studied, showing that ibuprofen had the lowest binding energy and the highest affinity for the docked genes. Both compounds revealed a mutual molecular interaction with amino acids residues within the catalytic cavity of the docked genes. These results provide new insights into the toxic effects of ibuprofen and dibutyl phthalate on Clarias gariepinus, contributing to a better understanding of the environmental impact of these pollutants.

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.

Genome-wide investigation of the DMRT gene family sheds new insight into the regulation of sex differentiation in spotted knifejaw (Oplegnathus punctatus) with fusion chromosomes (Y)

The role of the DMRT family in male sex determination and differentiation is significant, but its regulatory role in spotted knifejaw with Y fusion chromosomes remains unclear. Through genome-wide scanning, transcriptome analysis, qPCR, FISH, and RNA interference (RNAi), we investigated the DMRT family and the dmrt1-based sex regulation network. Seven DMRTs were identified (DMRT1/2 (2a,2b)/6, DMRT4/5, DMRT3), and dmrt gene dispersion among chromosomes is possibly driven by three whole-genome duplications. Transcriptome analysis enriched genes were associated with sex regulation and constructed a network associated with dmrt1. qPCR and FISH results showed the expression dimorphism of sex-related genes in dmrt-related regulatory networks. RNAi experiments indicated a distinct sex regulation mode in spotted knifejaw. Dmrt1 knockdown upregulated male-related genes (sox9a, sox9b, dmrt1, amh, amhr2) and hsd11b2 expression, which is critical for androgen synthesis. Amhr2 is located on the heterozygous chromosome (Y) and is specifically localized in primary spermatocytes, and is extremely upregulated after dmrt1 knockdown which suggested besides the important role of dmrt1 in male differentiation, the amhr2 along with amhr2/amh system, also play important regulatory roles in maintaining high expression of the hsd11b2 and male differentiation. This study aims to further investigate sex regulatory mechanisms in species with fusion chromosomes.

hsd17b1 is a key gene for ovarian differentiation of the Siberian sturgeon

This is the first work using gonads from undifferentiated, genetically-sexed Siberian sturgeon describing expression changes in genes related to steroid synthesis and female and male sex differentiation. One factor identified as relevant for ovarian differentiation was the gene coding for the enzyme Hsd17b1, which converts estrone into estradiol-17β. hsd17b1 was highly activated in female gonads at 2.5 months of age, around the onset of sex differentiation, preceding activation of two other genes involved in estrogen production (cyp19a1 and foxl2). hsd17b1 was also strongly repressed in males. Two known foxl2 paralogs are found in Siberian sturgeon-foxl2 and foxl2l-but only foxl2 appeared to be associated with ovarian differentiation. With regard to the male pathway, neither 11-oxygenated androgens nor classic male genes (amh, dmrt1, sox9, and dhh) were found to be involved in male sex differentiation, leaving open the question of which genes participate in early male gonad development in this ancient fish. Taken together, these results indicate an estrogen-dependence of female sex differentiation and 11-oxygenated androgen-independence of male sex differentiation.

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

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

Insights into Early Ontogenesis of Salmo salar: RNA Extraction, Housekeeping Gene Validation and Transcriptional Expression of Important Primordial Germ Cell and Sex-Determination Genes

The challenge in extracting high-quality RNA impedes the investigation of the transcriptome of developing salmonid embryos. Furthermore, the mRNA expression pattern of important PGC and SD genes during the initial embryonic development of Salmo salar is yet to be studied. So, in the present study, we aimed to isolate high-quality RNA from eggs and developing embryos to check vasa, dnd1, nanos3a, sdf1, gsdf, amh, cyp19a, dmrt1 and foxl2 expression by qPCR. Additionally, four HKGs (GAPDH, UB2L3, eEf1a and β-actin) were validated to select the best internal control for qPCR. High-quality RNA was extracted, which was confirmed by spectrophotometer, agarose gel electrophoresis and Agilent TapeStation analysis. UB2L3 was chosen as a reference gene because it exhibited lower intra- and inter-sample variation. vasa transcripts were expressed in all the developmental stages, while dnd1 was expressed only up to 40 d°C. Nanos3a was expressed in later stages and remained at its peak for a shorter period, while sdf1 showed an irregular pattern of mRNA expression. The mRNA expression levels of SD genes were observed to be upregulated during the later stages of development, prior to hatching. This study presents a straightforward methodology for isolating high-quality RNA from salmon eggs, and the resulting transcript profiles of significant PGC and SD genes in S. salar could aid in improving our comprehension of reproductive development in this commercially important species.

dmrtb1 is involved in the testicular development in Larimichthys crocea

In brief

dmrtb1 performs critical functions in sex determination/differentiation and gonadal development in many organisms, but its role in teleost is rarely studied. Through gene cloning, in situ hybridization, and RNA interference technology, the function of dmrtb1 in testicular development of large yellow croaker (Larimichthys crocea) was studied; our study will be helpful in understanding further the molecular regulation mechanism of Lcdmrtb1/Lcdmrt6 in testicular development in L. crocea, and our results enrich the theory of fish dmrts involved in reproductive regulation and provide a new idea for sex control breeding of L. crocea by manipulating reproductive pathway.

Abstract

Doublesex- and mab-3-related transcription factor B1 (dmrtb1/dmrt6) belongs to one of the members of DMRT family, which performs critical functions in sex determination and differentiation, gonadal development, and functional maintenance. However, knowledge of its exact mechanism remains unclear in teleost. Very little is known about the role of dmrtb1 in the gonad development of Larimichthys crocea. In this study, a dmrtb1 homolog in L. crocea named as Lcdmrtb1 with the full-length cDNA was isolated and characterized. Except for the conserved DM domain, the other regions had low homology. Of the tissues sampled, Lcdmrtb1 was only found to be highly expressed in the testis. In situ hybridization of testis revealed Lcdmrtb1 in both spermatogonia and spermatocytes. After Lcdmrtb1 interference in the testis cells (LYCT) of L. crocea, the expression levels of Lcdmrtb1 and Lcdmrt1 were significantly decreased; subsequently, testicular cell morphology changed from fibrous to round and their growth rate slowed. Similarly, the expression levels of Lcdmrtb1, Lcdmrt1, sox9a/b, and amh were significantly decreased after RNAi in the testis. Furthermore, it was discovered that the spermatogonia had disappeared, and the Sertoli cells had been reduced. The results of immunohistochemistry showed that the expression of Sox9 protein in the testis was not detected after dmrtb1 was knocked down. These results indicated that the absence of Lcdmrtb1 not only greatly inhibited cell growth and destroyed the morphology of testis cells but also down-regulated Lcdmrt1 expression in the testis. This study will be helpful in understanding further the molecular regulation mechanism of Lcdmrtb1/Lcdmrt6 in testicular development in L. crocea.

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.

Comprehensive Transcriptome Analysis of Gonadal and Somatic Tissues for Identification of Sex-Related Genes in the Largemouth Bass Micropterus salmoides

Largemouth bass (Micropterus salmoides) is an economically important fish. It can spawn many times during a breeding season, and there are no obvious morphological characteristics to distinguish male and female juvenile fish. So far, little is known about the genes regulating their sexual development in this species. Here, we performed RNA sequencing (RNA-Seq) analysis of the testis, ovary, and somatic tissue to identify sex-related genes in the largemouth bass. A total of 51,672 unigenes were obtained via the transcriptome analysis, and 5900 differential expression genes (DEGs), including 3028 up-regulated and 2872 down-regulated DEGs, were obtained in the somatic tissue, testis, and ovary. DEGs were retrieved by making comparisons: somatic tissue vs testis (1733-up and 1382-down), testis vs ovary (841-up and 807-down), and ovary vs somatic tissue (454-up and 683-down). Finally, functional annotation identified 22 key sex-related DEGs, including 13 testis-biased DEGs (dmrt1, cyp11b1, sox9, spata4, spata22, spata17, fshr, fem-1a, wt1, daz1, amh, vasa, and piwi1) and 9 ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, bmp15, fem-1b, fig. la, and piwi2). This result was further confirmed by the tissue expression detection via RT-PCR and RT-qPCR. Protein-protein interacting (PPI) network analysis revealed that the testis-specific dmrt1 interacts directly with the testis-biased DEGs (cyp11b1 and spata4) and the ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, and bmp15), suggesting that the dmrt1 as a sex-determining gene can play a dual role through inducing the testis-biased DEGs and inhibiting the ovary-biased DEGs during the testicular development. Our present results provide useful molecular data for a better understanding of sexual development in the largemouth bass.

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.

17β-Trenbolone binds to androgen receptor, decreases number of primordial germ cells, modulates expression of genes related to sexual differentiation, and affects sexual differentiation in zebrafish (Danio rerio)

Exposure to 17β-trenbolone caused a skewed sex ratio in fish. However, the molecular initiating event and key molecular event(s) remain unknown. In this study, zebrafish were exposed to 17β-trenbolone at nominal concentrations of 2 ng/L, 20 ng/L, 200 ng/L, and 2000 ng/L from fertilization to 60 days post fertilization (dpf). First, the sex ratio at 60 dpf was calculated to evaluate adverse outcomes on sexual differentiation. 17β-Trenbolone caused a skewed sex ratio toward males, with intersex individuals observed in the 20 ng/L group and all-male populations found in the 200 ng/L and 2000 ng/L groups. Then, the distribution and number of primordial germ cells, the expression of sex differentiation-related genes, and plasma vitellogenin concentrations were detected in wild-type zebrafish and the EGFP-nanos-3'UTR transgenic line using whole-mount in situ hybridization, real-time PCR, EGFP fluorescence quantification, and enzyme-linked immunosorbent assay. The results indicated that 17β-trenbolone exposure decreased the number of primordial germ cells at 1 dpf and 3 dpf, decreased expression of ovarian differentiation-related genes foxl2 and cyp19a1a at 60 dpf, increased expression of testis differentiation-related genes dmrt1, sox9a, and amh at 60 dpf, and decreased plasma vitellogenin levels at 60 dpf, revealing the key molecular events at different time points involved in affected sexual differentiation by 17β-trenbolone. Finally, molecular docking showed that 17β-trenbolone docked into ligand-binding domain of zebrafish androgen receptor with high binding energy (-3.72 kcal/mol), suggesting that binding to androgen receptor is the molecular initiating event affecting sexual differentiation by 17β-trenbolone. We found that 17β-trenbolone can bind to the zebrafish androgen receptor, decrease the number of primordial germ cells during the early embryonic stage, modulate the expression of genes related to sexual differentiation during gonadal differentiation, and eventually cause a skewed sex ratio toward males in adults.

Involvement of Transforming Growth Factor Beta Family Genes in Gonadal Differentiation in Japanese Eel, Anguilla japonica, According to Sex-Related Gene Expressions

The gonochoristic feature with environmental sex determination that occurs during the yellow stage in the eel provides an interesting model to investigate the mechanisms of gonadal development. We previously studied various sex-related genes during gonadal sex differentiation in Japanese eels. In the present study, the members of transforming growth factor beta (TGF-β) superfamily were investigated. Transcript levels of anti-Müllerian hormone, its receptor, gonadal soma-derived factor (amh, amhr2, and gsdf, respectively) measured by real-time polymerase chain reaction (qPCR) showed a strong sexual dimorphism. Transcripts were dominantly expressed in the testis, and their levels significantly increased with testicular differentiation. In contrast, the expressions of amh, amhr2, and gsdf transcripts were low in the ovary of E2-feminized female eels. In situ hybridization detected gsdf (but not amh) transcript signals in undifferentiated gonads. amh and gsdf signals were localized to Sertoli cells and had increased significantly with testicular differentiation. Weak gsdf and no amh signals were detected in early ovaries of E2-feminized female eels. Transcript levels of amh and gsdf (not amhr2) decreased during human chorionic gonadotropin (HCG)-induced spermatogenesis in males. This study suggests that amh, amhr2, and especially gsdf might be involved in the gene pathway regulating testicular differentiation of Japanese eels.

Characterization, expression and CpG methylation analysis of Dmrt1 and its response to steroid hormone in blotched snakehead (Channa maculata)

Dmrt1 is an important transcriptional regulator that plays critical role in male gonadogenesis, testicular differentiation and development. In this study, Dmrt1 was cloned from blotched snakehead (Channa maculata), which is designated as CmDmrt1. CmDmrt1 encoded a putative protein with 293 amino acids and presented an extremely conserved DM domain. It was nearly expressed in the gonads, and the expression was more than 15 times higher in the testis than in the ovary. 1851 bp promoter sequence of CmDmrt1 was characterized and the methylation levels of the CpG sites were analyzed to detect sex-related differences. A significant negative correlation between CmDmrt1 expression and CpG methylation level of its promoter was found in the testis and ovary. During gonadal development, CmDmrt1 transcription displayed strong male-biased expression patterns, increased with the maturation of testis and reached the peak at 195 days after hatching (dah), which indicates a significant role of Dmrt1 in spermatogenesis. Steroid treatment could influence CmDmrt1 expression, and long-term 17β-estradiol (E2) treatment could induce the male-to-female secondary sex reversal (SSR), which resulted in the differentiated testis transformed to ovary or ovotestis. Meanwhile, CmDmrt1 expression was down-regulated to fairly low level in the ovary of the SSR XY fish, which was similar to that in normal XX females ovary. Our research illustrates that Dmrt1 is linked to testis differentiation and spermatogenesis in blotched snakehead, providing information for functional studies on sex differentiation and gonadal development of C. maculata, and scientific basis for the production practice of all-male snakehead breeding.

Evolution of master sex determiners: TGF-β signalling pathways at regulatory crossroads

To date, more than 20 different vertebrate master sex-determining genes have been identified on different sex chromosomes of mammals, birds, frogs and fish. Interestingly, six of these genes are transcription factors (Dmrt1- or Sox3- related) and 13 others belong to the TGF-β signalling pathway (Amh, Amhr2, Bmpr1b, Gsdf and Gdf6). This pattern suggests that only a limited group of factors/signalling pathways are prone to become top regulators again and again. Although being clearly a subordinate member of the sex-regulatory network in mammals, the TGF-β signalling pathway made it to the top recurrently and independently. Facing this rolling wave of TGF-β signalling pathways, this review will decipher how the TGF-β signalling pathways cope with the canonical sex gene regulatory network and challenge the current evolutionary concepts accounting for the diversity of sex-determining mechanisms. 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)'.

Allelic diversification after transposable element exaptation promoted gsdf as the master sex determining gene of sablefish

Concepts of evolutionary biology suggest that morphological change may occur by rare punctual but rather large changes, or by more steady and gradual transformations. It can therefore be asked whether genetic changes underlying morphological, physiological, and/or behavioral innovations during evolution occur in a punctual manner, whereby a single mutational event has prominent phenotypic consequences, or if many consecutive alterations in the DNA over longer time periods lead to phenotypic divergence. In the marine teleost, sablefish (Anoplopoma fimbria), complementary genomic and genetic studies led to the identification of a sex locus on the Y Chromosome. Further characterization of this locus resulted in identification of the transforming growth factor, beta receptor 1a (tgfbr1a) gene, gonadal somatic cell derived factor (gsdf), as the main candidate for fulfilling the master sex determining (MSD) function. The presence of different X and Y Chromosome copies of this gene indicated that the male heterogametic (XY) system of sex determination in sablefish arose by allelic diversification. The gsdfY gene has a spatio-temporal expression profile characteristic of a male MSD gene. We provide experimental evidence demonstrating a pivotal role of a transposable element (TE) for the divergent function of gsdfY. By insertion within the gsdfY promoter region, this TE generated allelic diversification by bringing cis-regulatory modules that led to transcriptional rewiring and thus creation of a new MSD gene. This points out, for the first time in the scenario of MSD gene evolution by allelic diversification, a single, punctual molecular event in the appearance of a new trigger for male development.

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.

Brief exposure to warm temperatures reduces intron retention in Kdm6b in a species with temperature-dependent sex determination

Animals with temperature-dependent sex determination (TSD) respond to thermal cues during early embryonic development to trigger gonadal differentiation. TSD has primarily been studied using constant temperature incubations, where embryos are exposed to constant male- or female-producing temperatures, and these studies have identified genes that display sex-specific expression in response to incubation temperature. Kdm6b, a histone demethylase gene, has received specific attention as it is among the initial genes to respond to incubation temperature and is necessary for testis development. Interestingly, Kdm6b retains an intron when eggs are incubated at a constant male-producing temperature, but the role of thermal variability in this developmental process is relatively understudied. Species with TSD regularly experience thermal cues that fluctuate between male- and female-producing temperatures throughout development but it is unclear how Kdm6b responds to such variable temperatures. In this study, we investigate temperature-sensitive splicing in Kdm6b by exposing embryos to male- and female-producing thermal conditions. We show a rapid decrease in levels of the intron retaining transcript of Kdm6b upon exposure to female-producing conditions. These results demonstrate that, under ecologically relevant conditions, temperature-sensitive splicing can differentially regulate genes critical to TSD.

Comparison of Gonadal Transcriptomes Uncovers Reproduction-Related Genes with Sexually Dimorphic Expression Patterns in Diodon hystrix

Diodon hystrix is a new and emerging aquaculture species in south China. However, due to the lack of understanding of reproductive regulation, the management of breeding and reproduction under captivity remains a barrier for the commercial aquaculture of D. hystrix. More genetic information is needed to identify genes critical for gonadal development. Here, the first gonadal transcriptomes of D. hystrix were analyzed and 151.89 million clean reads were generated. All reads were assembled into 57,077 unigenes, and 24,574 could be annotated. By comparing the gonad transcriptomes, 11,487 differentially expressed genes were obtained, of which 4599 were upregulated and 6888 were downregulated in the ovaries. Using enrichment analyses, many functional pathways were found to be associated with reproduction regulation. A set of sex-biased genes putatively involved in gonad development and gametogenesis were identified and their sexually dimorphic expression patterns were characterized. The detailed transcriptomic data provide a useful resource for further research on D. hystrix reproductive manipulation.

Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus

The genetic regulatory network governing vertebrate gonadal differentiation appears less conserved than previously thought. Here, we investigated the gonadal development of Astyanax mexicanus Pachón cavefish by looking at primordial germ cells (PGCs) migration and proliferation, gonad histology, and gene expression patterns. We showed that PGCs are first detected at the 80% epiboly stage and then reach the gonadal primordium at 1 day post-fertilization (dpf). However, in contrast to the generally described absence of PGCs proliferation during their migration phase, PGCs number in cavefish doubles between early neurula and 8-9 somites stages. Combining both gonadal histology and vasa (germ cell marker) expression patterns, we observed that ovarian and testicular differentiation occurs around 65 dpf in females and 90 dpf in males, respectively, with an important inter-individual variability. The expression patterns of dmrt1, gsdf, and amh revealed a conserved predominant male expression during cavefish gonadal development, but none of the ovarian differentiation genes, i. e., foxl2a, cyp19a1a, and wnt4b displayed an early sexually dimorphic expression, and surprisingly all these genes exhibited predominant expression in adult testes. Altogether, our results lay the foundation for further research on sex determination and differentiation in A. mexicanus and contribute to the emerging picture that the vertebrate sex differentiation downstream regulatory network is less conserved than previously thought, at least in teleost fishes.

Dmrt1 regulates the immune response by repressing the TLR4 signaling pathway in goat male germline stem cells

Double sex and mab-3-related transcription factor 1 (Dmrt1), which is expressed in goat male germline stem cells (mGSCs) and Sertoli cells, is one of the most conserved transcription factors involved in sex determination. In this study, we highlighted the role of Dmrt1 in balancing the innate immune response in goat mGSCs. Dmrt1 recruited promyelocytic leukemia zinc finger (Plzf), also known as zinc finger and BTB domain-containing protein 16 (Zbtb16), to repress the Toll-like receptor 4 (TLR4)-dependent inflammatory signaling pathway and nuclear factor (NF)-κB. Knockdown of Dmrt1 in seminiferous tubules resulted in widespread degeneration of germ and somatic cells, while the expression of proinflammatory factors were significantly enhanced. We also demonstrated that Dmrt1 stimulated proliferation of mGSCs, but repressed apoptosis caused by the immune response. Thus, Dmrt1 is sufficient to reduce inflammation in the testes, thereby establishing the stability of spermatogenesis and the testicular microenvironment.

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.

Comparative Genomics Studies on the dmrt Gene Family in Fish

Doublesex and mab-3-related transcription factor (dmrt) genes are widely distributed across various biological groups and play critical roles in sex determination and neural development. Here, we applied bioinformatics methods to exam cross-species changes in the dmrt family members and evolutionary relationships of the dmrt genes based on genomes of 17 fish species. All the examined fish species have dmrt1-5 while only five species contained dmrt6. Most fish harbored two dmrt2 paralogs (dmrt2a and dmrt2b), with dmrt2b being unique to fish. In the phylogenetic tree, 147 DMRT are categorized into eight groups (DMRT1-DMRT8) and then clustered in three main groups. Selective evolutionary pressure analysis indicated purifying selections on dmrt1-3 genes and the dmrt1-3-2(2a) gene cluster. Similar genomic conservation patterns of the dmrt1-dmrt3-dmrt2(2a) gene cluster with 20-kb upstream/downstream regions in fish with various sex-determination systems were observed except for three regions with remarkable diversity. Synteny analysis revealed that dmrt1, dmrt2a, dmrt2b, and dmrt3-5 were relatively conserved in fish during the evolutionary process. While dmrt6 was lost in most species during evolution. The high conservation of the dmrt1-dmrt3-dmrt2(2a) gene cluster in various fish genomes suggests their crucial biological functions while various dmrt family members and sequences across fish species suggest different biological roles during evolution. This study provides a molecular basis for fish dmrt functional analysis and may serve as a reference for in-depth phylogenomics.

The effects of 11-ketotestosterone implants on transcript levels of gonadotropin receptors, and foxl2 and dmrt1 genes in the Previtellogenic ovary of cultured beluga (Huso huso)

The in vivo effect of 11-ketotestosterone (11KT) on transcript levels of the gonadotropin receptors (fshr and lhr) and sex differentiation-related genes (dmrt1 and foxl2) was examined in the ovaries of immature female beluga. For this purpose, six fish were treated with implants containing 2.5 mg 11KT and a placebo group of six females of the same age and gametogenic stage were given a blank implant. The implants were intraperitoneally inserted into 4-year-old females at the previtellogenic stage (mean body weight 5580 ± 165 g) and maintained under culture conditions for 8 weeks. Ovary samples for gene expression analysis of lhr, fshr, dmrt1 and foxl2 were collected by biopsy at 3 and 8 weeks post implantation. Diameters of oocytes increased in response to 11KT treatment, both at 3 and at 8 weeks post implantation, but no obvious changes were evident in cytology. Three weeks of 11KT treatment did not affect target gene expression, but a tendency for a time-dependent decrease of lhr and dmrt1 mRNA levels was observed in both treatment and placebo groups. By 8 weeks of treatment, however, 11KT implants provoked the upregulation of fshr and foxl2 transcript levels. Furthermore, lhr and dmrt1 transcript abundances recovered by 8 weeks of exposure in both blank- and 11KT-implanted beluga. These results suggest that 11KT, either directly or indirectly, may affect gametogenesis and regulate some key components of the reproductive axis in female beluga.

Expression Profiles of dmrt1 in Schizothorax kozlovi, and Their Relation to CpG Methylation of Its Promoter and Temperature

To elucidate the role of dmrt1 in sex differentiation of a teleost fish Schizothorax kozlovi, the full-length sequences of its cDNA and promoter were cloned by rapid amplification of cDNA ends (RACE) and genome walking. The relative mRNA expression levels were determined by quantitative real-time PCR (RT-PCR). The 1095-bp dmrt1 cDNA was predicted to encode a protein of 264 amino acids. It was expressed only in the gonads, and the expression was 17-times higher in the testis than in the ovary. The 1215-bp promoter sequence of dmrt1 was cloned and analyzed to detect sex-related differences in its methylation levels. A significant negative relationship between the dmrt1 expression and CpG methylation of its promoter were found in the testes and ovaries of S. kozlovi. Significant differences in dmrt1 expression levels were also found between the larval and juvenile stages. No significant differences in expression were found during the entire larval stage, and in the individuals among three different temperature groups (10°C, 14°C, and 18°C). Considering that the sex of sampled larval fish cannot be distinguished, correlations between dmrt1 expression and effects of temperature on sex differentiation in S. kozlovi need further study.

Spatial and temporal expression pattern of sex-related genes in ovo-testis of the self-fertilizing mangrove killifish (Kryptolebias marmoratus)

In vertebrates, sex determination and differentiation comprehend a fine balance between female and male factors, leading the bipotential anlage to develop towards ovary or testis, respectively. Nevertheless, the mangrove killifish, (Kryptolebias marmoratus) a simultaneous hermaphroditic species, could overcome those antagonistic pathways and evolved to develop and maintain reproductively active ovarian and testicular tissues in the same organ. Morphological and mRNA localization analyzes of developing and adult gonads demonstrate that genes related to testis (dmrt1 and amh) and ovary differentiation (foxl2 and sox9a) follow the same expression pattern observed in gonochoristic species, thus functioning as two independent organs. In addition, Amh expression patterns make it a strong candidate for initiation of the formation and maintenance of the testicular tissue in the hermaphroditic gonad. Differently from described so far, foxl3 seems to have an important role in oogenesis as well as spermatogenesis and gonadal structure.

Disruption of dmrt1 rescues the all-male phenotype of the cyp19a1a mutant in zebrafish - a novel insight into the roles of aromatase/estrogens in gonadal differentiation and early folliculogenesis

Sex determination and differentiation are complex processes controlled by many different factors; however, the relationships among these factors are poorly understood. Zebrafish gonadal differentiation exhibits high plasticity involving multiple factors and pathways, which provides an excellent model for investigating the interactions between them. Ovarian aromatase (cyp19a1a) and dmrt1 are key factors in directing vertebrate ovary and testis differentiation, respectively. Knockout of zebrafish cyp19a1a leads to all-male offspring, whereas the loss of dmrt1 results in a female-biased sex ratio. In the present study, we established dmrt1^-/- ;cyp19a1a^-/- double mutant zebrafish and discovered that the introduction of the dmrt1 mutation into the cyp19a1a mutant could rescue the all-male phenotype of the latter. Interestingly, despite the lack of aromatase/estrogens, the follicles in the ovary of the rescued cyp19a1a mutant could develop normally up to the previtellogenic stage. Further evidence suggested the ovarian aromatase directed ovarian differentiation by suppressing dmrt1 expression via nuclear estrogen receptors (nERs). Our results provide solid evidence for an interaction between cyp19a1a and dmrt1 in zebrafish gonadal differentiation, and for the dispensability of estrogens in controlling early folliculogenesis.

Molecular characterization, computational analysis and expression profiling of Dmrt1 gene in Indian major carp, Labeo rohita (Hamilton 1822)

Sexual dimorphism of fish morphology, physiology and behavior is diverse and complex in nature. Doublesex and mab-3 related transcription factor (Dmrt) is a large protein family whose function is sexual development and differentiation in vertebrates. Here, we report a full-length cDNA sequence of Labeo rohita (rohu) Dmrt1 of 907 bp length having 798 bp of open reading frame encoding 265 amino acids. The molecular weight of rohu DMRT1 protein was found to be 28.74 KDa and isoelectric point was 7.53. DMRT1 protein contains 23 positively and 24 negatively charged amino acids with a GRAVY score of -0.618. A characteristic DM domain was found in DMRT1 protein, which is a novel DNA-binding domain. Phylogenetic analysis showed maximum similarity with Cyprinus carpio when compared with DMRT1 of other vertebrates. Molecular docking study identified active sites to be targeted for drug designing. Rohu DMRT1 was observed to interact with other proteins such as FOXL2, CYP19a1a, AMH and SOX9a. Differential expression study revealed higher expression in testis tissue implying its role in male sex differentiation and testicular development. The information generated in the present work could facilitate further research to resolve the issues related to gonadal maturation and reproduction of commercially important aquaculture species.

Phylogeography, more than elevation, accounts for sex chromosome differentiation in Swiss populations of the common frog (Rana temporaria)

Sex chromosomes in vertebrates range from highly heteromorphic (as in most birds and mammals) to strictly homomorphic (as in many fishes, amphibians, and nonavian reptiles). Reasons for these contrasted evolutionary trajectories remain unclear, but species such as common frogs with polymorphism in the extent of sex chromosome differentiation may potentially deliver important clues. By investigating 92 common frog populations from a wide range of elevations throughout Switzerland, we show that sex chromosome differentiation strongly correlates with alleles at the candidate sex-determining gene Dmrt1. Y-specific Dmrt1 haplotypes cluster into two main haplogroups, Y_A and Y_B , with a phylogeographic signal that parallels mtDNA haplotypes: Y_A populations, with mostly well-differentiated sex chromosomes, occur primarily south of the main alpine ridge that bisects Switzerland, whereas Y_B populations, with mostly undifferentiated (proto-)sex chromosomes, occur north of this ridge. Elevation has only a marginal effect, opposing previous suggestions of a major role for climate on sex chromosome differentiation. The Y-haplotype effect might result from differences in the penetrance of alleles at the sex-determining locus (such that sex reversal and ensuing X-Y recombination are more frequent in Y_B populations), and/or fixation of an inversion on Y_A (as supported by the empirical observation that Y_A haplotypes might not recombine in XY_A females).

Comparative transcriptome analysis of male and female gonads reveals sex-biased genes in spotted scat (Scatophagus argus)

Scatophagus argus is a new emerging aquaculture fish in East and Southeast Asia. To date, research on reproductive development and regulation in S. argus is lacking. Additionally, genetic and genomic information about reproduction, such as gonadal transcriptome data, is also lacking. Herein, we report the first gonadal transcriptomes of S. argus and identify genes potentially involved in reproduction and gonadal development. A total of 136,561 unigenes were obtained by sequencing of testes (n = 3) and ovaries (n = 3) at stage III. Genes upregulated in males and females known to be involved in gonadal development and gametogenesis were identified, including male-biased dmrt1, amh, gsdf, wt1a, sox9b, and nanos2, and female-biased foxl2, gdf9, bmp15, sox3, zar1, and figla. Serum estradiol-17β and 11-ketotestosterone levels were biased in female and male fish, respectively. Sexual dimorphism of serum steroid hormone levels were interpreted after expression analysis of 20 steroidogenesis-related genes, including cyp19a1a and cyp11b2. This gonadal transcript dataset will help investigate functional genes related to reproduction in S. argus.

Transcriptional effects in the estuarine guppy Poecilia vivipara exposed to sanitary sewage in laboratory and in situ

The urban growth has increased sanitary sewage discharges in coastal ecosystems, negatively affecting the aquatic biota. Mangroves, one of the most human-affected coastal biomes, are areas for reproduction and nursing of several species. In order to evaluate the effects of sanitary sewage effluents in mangrove species, this study assessed the hepatic transcriptional responses of guppy fish Poecilia vivipara exposed to sanitary sewage 33% (v:v), using suppressive subtraction hybridization (SSH), high throughput sequencing of RNA (Ion-proton) and quantification of transcript levels by qPCR of some identified genes in fish kept in a sewage-contaminated environment. Genes identified are related predominantly to xenobiotic biotransformation, immune system and sexual differentiation. The qPCR results confirmed the induction of cytochrome P450 1A (CYP1A), glutathione S transferase A-like (GST A-like) methyltransferase (MET) and UDP glycosyltransferase 1A (UDPGT1A), and repression of complement component C3 (C3), doublesex and mab-3 related transcription factor 1 (DMRT1), and transferrin (TF) in the laboratory experiment. In the field exposure, the transcript levels of CYP1A, DMRT1, MET, GST A-like and UDPGT1A were higher in fishes exposed at the contaminated sites compared to the reference site. Chemical analysis in fish from the laboratory and in situ experiments, and surface sediment from the sewage-contaminated sites revealed relevant levels of polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyl (PCBs) and linear alkylbenzenes (LABs). These data reinforce the use of P. vivipara as a sentinel for monitoring environmental contamination in coastal regions.

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.

Dmrt1 (doublesex and mab-3-related transcription factor 1) expression during gonadal development and spermatogenesis in the Japanese eel

Dmrt1, doublesex- and mab-3-related transcription factor-1, has been suggested to play critical roles in male gonadogenesis, testicular differentiation and development, including spermatogenesis, among different vertebrates. Vasa is a putative molecular marker of germ cells in vertebrates. In this study, we cloned the full-length dmrt1 cDNA from Japanese eel, and the protein comprised 290 amino acids and presented an extremely conserved Doublesex and Mab-3 (DM) domain. Vasa proteins were expressed in gonadal germ cells in a stage-specific manner, and were expressed at high levels in PGC and spermatogonia, low levels in spermatocytes, and were absent in spermatids and spermatozoa of Japanese eels. Dmrt1 proteins were abundantly expressed in spermatogonia B cells, spermatocytes, spermatids, but not in spermatozoa, spermatogonia A and Sertoli cells. To our knowledge, this study is the first to show a restricted expression pattern for the Dmrt1 protein in spermatogonia B cells, but not spermatogonia A cells, of teleosts. Therefore, Dmrt1 might play vital roles at the specific stages during spermatogenesis from spermatogonia B cells to spermatids in the Japanese eel. Moreover, the Dmrt1 protein exhibited a restricted localization in differentiating oogonia in the early differentiating gonad (ovary-like structure) of male Japanese eels and in E2-induced feminized Japanese eels. We proposed that dmrt1 may be not only required for spermatogenesis but might also play a role in oogenesis in the Japanese eel.

Conservation and diversity in expression of candidate genes regulating socially-induced female-male sex change in wrasses

Fishes exhibit remarkably diverse, and plastic, patterns of sexual development, most striking of which is sequential hermaphroditism, where individuals readily reverse sex in adulthood. How this stunning example of phenotypic plasticity is controlled at a genetic level remains poorly understood. Several genes have been implicated in regulating sex change, yet the degree to which a conserved genetic machinery orchestrates this process has not yet been addressed. Using captive and in-the-field social manipulations to initiate sex change, combined with a comparative qPCR approach, we compared expression patterns of four candidate regulatory genes among three species of wrasses (Labridae)-a large and diverse teleost family where female-to-male sex change is pervasive, socially-cued, and likely ancestral. Expression in brain and gonadal tissues were compared among the iconic tropical bluehead wrasse (Thalassoma bifasciatum) and the temperate spotty (Notolabrus celidotus) and kyusen (Parajulus poecilepterus) wrasses. In all three species, gonadal sex change was preceded by downregulation of cyp19a1a (encoding gonadal aromatase that converts androgens to oestrogens) and accompanied by upregulation of amh (encoding anti-müllerian hormone that primarily regulates male germ cell development), and these genes may act concurrently to orchestrate ovary-testis transformation. In the brain, our data argue against a role for brain aromatase (cyp19a1b) in initiating behavioural sex change, as its expression trailed behavioural changes. However, we find that isotocin (it, that regulates teleost socio-sexual behaviours) expression correlated with dominant male-specific behaviours in the bluehead wrasse, suggesting it upregulation mediates the rapid behavioural sex change characteristic of blueheads and other tropical wrasses. However, it expression was not sex-biased in temperate spotty and kyusen wrasses, where sex change is more protracted and social groups may be less tightly-structured. Together, these findings suggest that while key components of the molecular machinery controlling gonadal sex change are phylogenetically conserved among wrasses, neural pathways governing behavioural sex change may be more variable.

Characterization of the European Sea Bass (Dicentrarchus labrax) Gonadal Transcriptome During Sexual Development

The European sea bass is one of the most important cultured fish in Europe and has a marked sexual growth dimorphism in favor of females. It is a gonochoristic species with polygenic sex determination, where a combination between still undifferentiated genetic factors and environmental temperature determines sex ratios. The molecular mechanisms responsible for gonadal sex differentiation are still unknown. Here, we sampled fish during the gonadal developmental period (110 to 350 days post fertilization, dpf), and performed a comprehensive transcriptomic study by using a species-specific microarray. This analysis uncovered sex-specific gonadal transcriptomic profiles at each stage of development, identifying larger number of differentially expressed genes in ovaries when compared to testis. The expression patterns of 54 reproduction-related genes were analyzed. We found that hsd17β10 is a reliable marker of early ovarian differentiation. Further, three genes, pdgfb, snx1, and nfy, not previously related to fish sex differentiation, were tightly associated with testis development in the sea bass. Regarding signaling pathways, lysine degradation, bladder cancer, and NOD-like receptor signaling were enriched for ovarian development while eight pathways including basal transcription factors and steroid biosynthesis were enriched for testis development. Analysis of the transcription factor abundance showed an earlier increase in females than in males. Our results show that, although many players in the sex differentiation pathways are conserved among species, there are peculiarities in gene expression worth exploring. The genes identified in this study illustrate the diversity of players involved in fish sex differentiation and can become potential biomarkers for the management of sex ratios in the European sea bass and perhaps other cultured species.

Genomic organization and sexually dimorphic expression of the Dmrt1 gene in largemouth bass (Micropterus salmoides)

Doublesex and Mab-3 related transcription factor (Dmrt) genes play important roles in the process of sex determination and differentiation. In this study, a Dmrt1 gene open reading frame sequence was obtained from the gonadal transcriptome data of largemouth bass (Micropterus salmoides), and identified by cloning and sequencing. The ORF of Dmrt1 is 900 bp long, encodes 298 amino acids, and contains the DM region which is characteristic of Dmrt1. Full gDNA sequence of Dmrt1 was composed of five exons and four introns. RT-PCR and Q-PCR analysis of Dmrt1 were conducted in eight tissues and three developmental stages of mature male and female individuals. In situ hybridization was used to locate the expression of Dmrt1 in the testis and ovary of largemouth bass. The results showed that Dmrt1 was highly expressed in the testis of mature fish, but only weakly expressed in other tissues such as heart, liver, and brain, and exhibited gender dimorphism in the gonads of male and female fish at different stages. Furthermore, the expression level in female fish was very low and decreased gradually with ovary maturation. In situ hybridization indicated positive signals were found in early oocytes, but not in mature oocytes, while strong positive signals were found in all types of mature testis cells. The study showed that the sequence and structure of Dmrt1 were highly conserved and exhibited significant gender dimorphism in largemouth bass, as in other fish species. It is suggested that Dmrt1 plays an important role in sex determination and differentiation in largemouth bass.

A novel evolutionary conserved mechanism of RNA stability regulates synexpression of primordial germ cell-specific genes prior to the sex-determination stage in medaka

Dmrt1 is a highly conserved transcription factor, which is critically involved in regulation of gonad development of vertebrates. In medaka, a duplicate of dmrt1-acting as master sex-determining gene-has a tightly timely and spatially controlled gonadal expression pattern. In addition to transcriptional regulation, a sequence motif in the 3' UTR (D3U-box) mediates transcript stability of dmrt1 mRNAs from medaka and other vertebrates. We show here that in medaka, two RNA-binding proteins with antagonizing properties target this D3U-box, promoting either RNA stabilization in germ cells or degradation in the soma. The D3U-box is also conserved in other germ-cell transcripts, making them responsive to the same RNA binding proteins. The evolutionary conservation of the D3U-box motif within dmrt1 genes of metazoans-together with preserved expression patterns of the targeting RNA binding proteins in subsets of germ cells-suggest that this new mechanism for controlling RNA stability is not restricted to fishes but might also apply to other vertebrates.

Sex-Chromosome Recombination in Common Frogs Brings Water to the Fountain-of-Youth

According to the canonical model of sex-chromosome evolution, the degeneration of Y or W chromosomes (as observed in mammals and birds, respectively) results from an arrest of recombination in the heterogametic sex, driven by the fixation of sexually antagonistic mutations. However, sex chromosomes have remained homomorphic in many lineages of fishes, amphibians, and nonavian reptiles. According to the "fountain-of-youth" model, this homomorphy results from occasional events of sex reversal. If recombination arrest in males is controlled by maleness per se (and not by genotype), then Y chromosomes are expected to recombine in XY females, preventing their long-term degeneration. Here, we provide field support for the fountain-of-youth, by showing that sex-chromosome recombination in Rana temporaria only depends on phenotypic sex: naturally occurring XX males show the same restriction of recombination as XY males (average map length ∼2 cM), while XY females recombine as much as XX females (average map length ∼150 cM). Our results challenge several common assumptions regarding the evolution of sex chromosomes, including the role of sexually antagonistic genes as drivers of recombination arrest, and that of chromosomal inversions as underlying mechanisms.

Molecular identification of dmrt1 and its promoter CpG methylation in correlation with gene expression during gonad development in Culter alburnus

Dmrt1, a member of the Dmrt family, is an important transcription regulator of gender determination. To study the biological function of dmrt1 in sexual differentiation and its potential implication in breeding technology, we obtained the full-length cDNA and proximal promoter sequence of dmrt1 in Culter alburnus, and analyzed the impact of promoter CpG methylation on the gene expression pattern of dmrt1 during gonad development. Dmrt1 was 922 bp in length and consisted a 150 bp 5'-UTR, a 28 bp 3'-UTR, and a 744 bp open reading frame (ORF). Based on the coding sequence of the dmrt1 gene, the deduced amino acid sequence was detected, and the protein structure of this gene was predicted in C. alburnus. The results indicate that the structure and function of dmrt1 were highly conservative compared to other vertebrates. The expression level of dmrt1 mRNA in different tissues was explored by qRT-PCR, which was only highly expressed in the testes and almost undetectable in other tissues. The CpG methylation pattern of the dmrt1 promoter was studied using DNA sequencing of sodium bisulfite in adult testes and ovaries, and it was found that dmrt1 promoter CpGs were not methylated in the testes, whereas hypermethylated in the ovaries. These findings demonstrate that DNA methylation can regulate sexual dimorphic expression of dmrt1, and therefore epigenetic modifications may play a critical role in the gonad differentiation of C. alburnus.

Natural sex change in fish

Sexual fate can no longer be considered an irreversible deterministic process that once established during early embryonic development, plays out unchanged across an organism's life. Rather, it appears to be a dynamic process, with sexual phenotype determined through an ongoing battle for supremacy between antagonistic male and female developmental pathways. That sexual fate is not final and is actively regulated via the suppression or activation of opposing genetic networks creates the potential for flexibility in sexual phenotype in adulthood. Such flexibility is seen in many fish, where sex change is a usual and adaptive part of the life cycle. Many fish are sequential hermaphrodites, beginning life as one sex and changing sometime later to the other. Sequential hermaphrodites include species capable of female-to-male (protogynous), male-to-female (protandrous), or bidirectional (serial) sex change. These natural forms of sex change involve coordinated transformations across multiple biological systems, including behavioral, anatomical, neuroendocrine and molecular axes. Here we review the biological processes underlying this amazing transformation, focusing particularly on the molecular aspects, where new genomic technologies are beginning to help us understand how sex change is initiated and regulated at the molecular level.

The Gene Toolkit Implicated in Functional Sex in Sparidae Hermaphrodites: Inferences From Comparative Transcriptomics

Sex-biased gene expression is the mode through which sex dimorphism arises from a nearly identical genome, especially in organisms without genetic sex determination. Teleost fishes show great variations in the way the sex phenotype forms. Among them, Sparidae, that might be considered as a model family displays a remarkable diversity of reproductive modes. In this study, we sequenced and analyzed the sex-biased transcriptome in gonads and brain (the tissues with the most profound role in sexual development and reproduction) of two sparids with different reproductive modes: the gonochoristic common dentex, Dentex dentex, and the protandrous hermaphrodite gilthead seabream, Sparus aurata. Through comparative analysis with other protogynous and rudimentary protandrous sparid transcriptomes already available, we put forward common male and female-specific genes and pathways that are probably implicated in sex-maintenance in this fish family. Our results contribute to the understanding of the complex processes behind the establishment of the functional sex, especially in hermaphrodite species and set the groundwork for future experiments by providing a gene toolkit that can improve efforts to control phenotypic sex in finfish in the ever-increasingly important field of aquaculture.

Differential Expression of dmrt1 in Astyanax scabripinnis (Teleostei, Characidade) Is Correlated with B Chromosome Occurrence

Astyanax is an abundant fish genus in South America. Some species of this group are characterized by the presence of B chromosomes and absence of morphologically differentiated sex chromosomes. In this study, we used quantitative real-time polymerase chain reaction to characterize mRNA expression of dmrt1 in Astyanax scabripinnis gonads. Maturing gonads of males with the B chromosome overexpressed dmrt1. Our findings suggest that B chromosomes may have an adaptive role in A. scabripinnis sex determination and maintenance.

A comparative transcriptomic study on developmental gonads provides novel insights into sex change in the protandrous black porgy (Acanthopagrus schlegelii)

Protandrous black porgy (Acanthopagrus schlegelii) is a popular and valuable commercial marine fish in China and East Asian countries. Controlling and managing its breeding has been an imperative step towards obtaining a sustainable supply of this fish in aquaculture production systems. Therefore, study on the molecular mechanisms of sex change in black porgy has both scientific and commercial importance. Previously, we identified some candidate genes related to sex determination and differentiation from a high-quality genome assembly of the black porgy. In the present study, transcriptome sequencing of developmental gonads (including testis, ovotestis and ovary) of black porgy was performed to further investigate the sex-change mechanisms. Our results showed that the highly expressed male-related genes (dmrt1, piwi1, piwi2, sox9, sox30 and amh) at the male phase were significantly down-regulated to a substantial degree at the intersexual stage, and the female-related genes (jnk1, vasa, wnt4, figla and foxl2) were distinctly up-regulated when the fish grows into a female adult, suggesting the potential roles of these genes in sex change of the black porgy. These data also support a previous hypothesis that the femaleness will be switched on when the testis is entering the degenerated stage through the diminished dmrt1 expression. Our transcriptome data provide a very useful genomic resource for future studies on sex change and practical aquaculture in the black porgy.

Identification of genes involved in gonadal sex differentiation and the dimorphic expression pattern in Takifugu rubripes gonad at the early stage of sex differentiation

Quantifying the expression of mRNAs in the gonads at the critical stage of molecular sex differentiation stage might help to clarify the regulatory network during early sex differentiation and provide new information on the role of sex-related genes in gonadal function. In this study, transcriptomic analysis of sex-related genes expression profiles in fugu gonads at 60 and 90 days after hatching (dah) was conducted firstly, and a total of 112,504,991 clean reads, encompassing 28.35 Gb of sequences were retrieved. Twenty-three thousand eight hundred ten genes were found to be expressed in juvenile fugu gonads, and we mainly focused on the differentially expressed genes that have the potential to be involved in the gonadal sex differentiation. For 60-dah juveniles, we identified 1014 genes that were upregulated in the ovary and 1570 that were upregulated in the testis. For 90-dah juveniles, we identified 1287 genes that were upregulated in the ovary and 1500 that were upregulated in the testis. The dimorphic expression patterns of 15 genes in gonads at 30 and 40 dah were further investigate using qPCR. Cyp11b and star were expressed at higher levels in XY than in XX, while cyp11a1 and cyp19a1a were expressed at higher levels in XX than in XY at 30 dah. At 40 dah, the levels of gsdf, dmrt1, dmrt3, cyp11c1, star, and hsd3b expression were higher in XY, while the levels of foxl2, cyp19a1a, wnt9b, and foxD4 expression were higher in XX. Sox9, cyp11a1, cyp17a1, cyp17a2, and nr5a2 were expressed at similar levels in XX and XY at 40 dah. This is the first report of gonadal transcriptome of fugu at early sex differentiation stage, and our results provide an archive for further study on molecular mechanism underlying sex differentiation in this species.

Hypothalamus and pituitary transcriptome profiling of male and female Hong Kong grouper (Epinephelus akaara)

Hong Kong grouper (Epinephelus akaara) is an important commercially cultured marine fish in Asia, and a protogynous hermaphrodite with the "diandry" pattern. In order to explore the gene expression patterns of hypothalamus and pituitary between male and female Hong Kong grouper, we used RNA-seq technology to investigate transcriptomes of both tissues in immature and mature male and female adults. This produced 227,227,148 and 215,858,948 high quality reads from hypothalamus and pituitary, which were jointly assembled into 199,203 unigenes. Among them, 30,786 unigenes were mapped to known genes. Differential expression analysis revealed 275 unigenes that were differentially expressed between immature male and female adults and 561 between mature male and female adults. According to annotation and KEGG information, these differentially expressed genes (DEGs) were involved in development, metabolism, and regulation of transcription. One DEG, amino-terminal enhancer of split (AES), a member of the Groucho/transducin-like enhancer of split family of transcriptional regulators that played important roles in neurogenesis, segmentation, and sex determination, was significantly upregulated in male individuals in both immature and mature adult comparisons, indicating it may be involved in male reproductive function during development. Our report, for the first time, uses RNA-seq technology to investigate transcriptomes of both hypothalamus and pituitary in teleost fish, and provides a basis for further studies of molecular mechanism of sex determination and development in Hong Kong grouper.

Primary males guide the femaleness through the regulation of testicular Dmrt1 and ovarian Cyp19a1a in protandrous black porgy

Controlling the development of the sexes is critically important for the broodstock management in aquaculture. Sex steroids are widely used for sex control of fish. However, hermaphroditic fish have a plastic sex, and a stable sex is difficult to maintain with sex steroids. We used the black porgy (Acanthopagrus schlegelii) as a model to understand the possible mechanism of sexual fate decision. Low exogenous estradiol (E2) induced male development. In contrast, high exogenous E2 induced the regression of the testis and the development of the ovary and resulted in an unstable expression of femaleness (passive femaleness, with ovaries containing only the primary oocytes). The removal of testicular tissue by surgery resulted in the early development of vitellogenic oocytes and active femaleness. Our data also demonstrated that the male-to-female sex change is blocked by the maintenance of male function with gonadotropin-induced dmrt1 expression in the testis. Furthermore, our data also indicated that ovarian cyp19a1a expression is regulated by the testis through epigenetic modifications. Therefore, the primary male guides the femaleness in the protandrous black porgy and the transition of sexual fate from male to female is determined by the status of the testicular tissue.