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

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

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

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

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

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.

Transcriptomic Analysis on the Effects of Altered Water Temperature Regime on the Fish Ovarian Development of Coreius guichenoti under the Impact of River Damming

Field investigation indicated that the reduction in fish spawning was associated with the alteration in water temperatures, even a 2-3 °C monthly difference due to reservoir operations. However, the physiological mechanism that influences the development of fish ovary (DFO) remains unclear. Thus, experiments of Coreius guichenoti were conducted at three different temperatures, optimal temperature (~20 °C, N) for fish spawning, lower (~17 °C, L), and higher (~23 °C, H), to reveal the effects of altered water temperature on the DFO. Comparisons were made between the L and N (LvsN) conditions and H and N (HvsN) conditions. Transcriptomic analysis differentially expressed transcripts (DETs) related to heat stress were observed only in LvsN conditions, indicating that the DFO showed a stronger response to changes in LvsN than in HvsN conditions. Upregulation of DETs of vitellogenin receptors in N temperature showed that normal temperature was conducive to vitellogenin entry into the oocytes. Other temperature-sensitive DETs, including microtubule, kinesin, dynein, and actin, were closely associated with cell division and material transport. LvsN significantly impacted cell division and nutrient accumulation in the yolk, whereas HvsN only influenced cell division. Our results highlight the impact of altered water temperature on the DFO, thereby providing insights for future reservoir operations regarding river damming and climate change and establishing fish conservation measures.

Thermal Stress Induces Metabolic Responses in Juvenile Qingtian Paddy Field Carp Cyprinus carpio var qingtianensis

Extreme fluctuations in water temperature lead to significant economic losses for the aquaculture industry. Cyprinus carpio var qingtianensis (locally called Qingtian paddy field carp), is a local variety commonly found in Zhejiang province, China. Unlike traditional aquaculture environments, the water temperature range between day and night in the rice field environment is much larger, and the high temperature in summer may exceed the growth threshold of fish because there is no manual intervention; therefore, the study of how the Qingtian paddy field carp (PF carp) adapts to high-temperature conditions can shed light how the species adapt to the rice field environment. To investigate the molecular mechanisms of this fish under thermal stress, the liver metabolomics of Qiangtian paddy field carp (PF carp) were analyzed. In this study, metabolomics was used to examine the metabolic reaction of PF carp (102 days old, 104.69 ± 3.08 g in weight, 14.65 ± 0.46 cm in length) at water temperatures of 28 °C (control group, CG), 34 °C (experimental group (EG) 34), and 38 °C (EG38). The results show that 175 expression profile metabolites (DEMs), including 115 upregulated and 60 downregulated metabolites, were found in the CG vs. EG34. A total of 354 DEMs were inspected in CG vs. EG38, with 85 metabolites downregulated and 269 metabolites upregulated. According to the pathway enrichment study, various pathways were altered by thermal stress, including those of lipid, amino-acid, and carbohydrate metabolism. Our study presents a potential metabolic profile for PF carp under thermal stress. It also demonstrates how the host responds to thermal stress on a metabolic and molecular level.

Molecular Cloning and Expression Responses of Jarid2b to High-Temperature Treatment in Nile Tilapia (Oreochromis niloticus)

Nile tilapia is a GSD + TE (Genetic Sex Determination + Temperature Effect) fish, and high-temperature treatment during critical thermosensitive periods (TSP) can induce the sex reversal of Nile tilapia genetic females, and brain transcriptomes have revealed the upregulation of Jarid2 (Jumonji and AT-rich domain containing 2) expression after 36 °C high-temperature treatment for 12 days during TSP. It was shown that JARID2 forms a complex with polycomb repressive complex 2 (PRC2) that catalyzed H3K27me3, which was strongly associated with transcriptional repression. In this study, Jarid2b was cloned and characterized in Nile tilapia, which was highly conserved among the analyzed fish species. The expression of Jarid2b was upregulated in the gonad of 21 dpf XX genetic females after 12-day high-temperature treatment and reached a similar level to that of males. Similar responses to high-temperature treatment also appeared in the brain, heart, liver, muscle, eye, and skin tissues. Interestingly, Jarid2b expression was only in response to high-temperature treatment, and not to 17α-methyltestosterone (MT) or letrozole treatments; although, these treatments can also induce the sex reversal of genetic Nile tilapia females. Further studies revealed that Jarid2b responded rapidly at the 8th hour after high-temperature treatment. Considering that JARID2 can recruit PRC2 and establish H3K27me3, we speculated that it might be an upstream gene participating in the regulation of Nile tilapia GSD + TE through regulating the H3K27 methylation level at the locus of many sex differentiation-related genes.

Comparative transcriptome analysis of heat-induced domesticated zebrafish during gonadal differentiation

BACKGROUND

The influence of environmental factors, especially temperature, on sex ratio is of great significance to elucidate the mechanism of sex determination. However, the molecular mechanisms by which temperature affects sex determination remains unclear, although a few candidate genes have been found to play a role in the process. In this study, we conducted transcriptome analysis of the effects induced by high temperature on zebrafish during gonad differentiation period.

RESULTS

Totals of 1171, 1022 and 2921 differentially expressed genes (DEGs) between high temperature and normal temperature were identified at 35, 45 and 60 days post-fertilization (dpf) respectively, revealing that heat shock proteins (HSPs) and DNA methyltransferases (DNMTs) were involved in the heat-exposed sex reversal. The Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway that were enriched in individuals after heat treatment included Fanconi anemia (FA) pathway, cell cycle, oocyte meiosis and homologous recombination.

CONCLUSIONS

Our study provides the results of comparative transcriptome analyses between high temperature and normal temperature, and reveals that the molecular mechanism of heat-induced masculinization in zebrafish is strongly related to the expression of HSPs and DNMTs and FA pathway during gonad differentiation.

A Transcriptomic Analysis of Gonads from the Low-Temperature-Induced Masculinization of Takifugu rubripes

Simple Summary

Our study analyzed the differentiation of transcriptomes in normal and sex-reverse Takifugu rubripes, and screened out 13 differentially expressed genes related to sex differentiation. This is the first report on the gonadal transcriptome of pseudo-males in Takifugu rubripes. Our results provide an important contribution to the molecular mechanism of masculinization in a cultured fish subject to low-temperature treatment.

Abstract

The phenotypic sex of fish is usually plastic. Low-temperature treatment induces the masculinization of Takifugu rubripes, resulting in pseudo-males (PM) with the physiological sex of a male (M) and genetic sex of a female (F). For a comparison of gonadal transcriptomes, we collected gonads from three groups of T. rubripes (F, M, and PM) for high-throughput transcriptome sequencing. The results provided 467,640,218 raw reads (70.15 Gb) and a total of 436,151,088 clean reads (65.43 Gb), with an average length of 150 bp. Only 79 differentially expressed genes (DEGs) were identified between F and PM, whereas 12,041 and 11,528 DEGs were identified between F and M, and PM and M, respectively. According to the functional annotation of DEGs, 13 DEGs related to gonadal development were screened (LOC101066759, dgat1, limk1, fbxl3, col6a3, fgfr3, dusp22b, svil, abhd17b, srgap3, tmem88b, bud4, and mustn10) which might participate in formating PM. A quantitative PCR of the DEGs confirmed the reliability of the RNA-seq. Our results provide an important contribution to the genome sequence resources for T. rubripes and insight into the molecular mechanism of masculinization in a cultured fish subject to low-temperature treatment.

Tissue and Temperature-Specific RNA-Seq Analysis Reveals Genomic Versatility and Adaptive Potential in Wild Sea Turtle Hatchlings (Caretta caretta)

Simple Summary

Digital transcriptomics is rapidly emerging as a powerful new technology for modelling the environmental dynamics of the adaptive landscape in diverse lineages. This is particularly valuable in taxa such as turtles and tortoises (order Testudines) which contain a large fraction of endangered species at risk due to anthropogenic impacts on the environment, including pollution, overharvest, habitat degradation, and climate change. Sea turtles (family Cheloniidae) in particular invite a genomics-enabled approach to investigating their remarkable portfolio of adaptive evolution. Our de novo transcriptome assemblies and measurements of tissue- and temperature- specific global gene expression in the loggerhead sea turtle (Caretta caretta) reveal the genomic basis for potential resiliency in this endangered flagship species, and are crucial to future management and conservation strategies with attention to changing climates. We summarize the interactions among differentially expressed genes by producing network visualizations, and highlight the shared biological pathways related to development, migration, immunity, and longevity reported in the avian and reptilian literature. Our original results for loggerhead sea turtles provide a large, new comparative genomic resource for the investigation of genotype–phenotype relationships in amniotes.

Abstract

Background: Digital transcriptomics is rapidly emerging as a powerful new technology for modelling the environmental dynamics of the adaptive landscape in diverse lineages. This is particularly valuable in taxa such as turtles and tortoises (order Testudines) which contain a large fraction of endangered species at risk due to anthropogenic impacts on the environment, including pollution, overharvest, habitat degradation, and climate change. Sea turtles (family Cheloniidae) in particular invite a genomics-enabled approach to investigating their remarkable portfolio of adaptive evolution. The sex of the endangered loggerhead sea turtle (Caretta caretta) is subject to temperature-dependent sex determination (TSD), a mechanism by which exposure to temperatures during embryonic development irreversibly determines sex. Higher temperatures produce mainly female turtles and lower temperatures produce mainly male turtles. Incubation temperature can have long term effects on the immunity, migratory ability, and ultimately longevity of hatchlings. We perform RNA-seq differential expression analysis to investigate tissue- and temperature-specific gene expression within brain (n = 7) and gonadal (n = 4) tissue of male and female loggerhead hatchlings. Results: We assemble tissue- and temperature-specific transcriptomes and identify differentially expressed genes relevant to sexual development and life history traits of broad adaptive interest to turtles and other amniotic species. We summarize interactions among differentially expressed genes by producing network visualizations, and highlight shared biological pathways related to migration, immunity, and longevity reported in the avian and reptile literature. Conclusions: The measurement of tissue- and temperature-specific global gene expression of an endangered, flagship species such as the loggerhead sea turtle (Caretta caretta) reveals the genomic basis for potential resiliency and is crucial to future management and conservation strategies with attention to changing climates. Brain and gonadal tissue collected from experimentally reared loggerhead male and female hatchlings comprise an exceedingly rare dataset that permits the identification of genes enriched in functions related to sexual development, immunity, longevity, and migratory behavior and will serve as a large, new genomic resource for the investigation of genotype–phenotype relationships in amniotes.

Developmental aspects of the hypothalamic-pituitary network related to reproduction in teleost fish

The hypothalamic-pituitary-gonadal axis is the main system that regulates reproduction in vertebrates through a complex network that involves different neuropeptides, neurotransmitters, and pituitary hormones. Considering that this axis is established early on life, the main goal of the present work is to gather information on its development and the actions of its components during early life stages. This review focuses on fish because their neuroanatomical characteristics make them excellent models to study neuroendocrine systems. The following points are discussed: i) developmental functions of the neuroendocrine components of this network, and ii) developmental disruptions that may impact adult reproduction. The importance of the components of this network and their susceptibility to external/internal signals that can alter their specific early functions and/or even the establishment of the reproductive axis, indicate that more studies are necessary to understand this complex and dynamic network.

Where the Ends Meet: An Overview of Sex Determination in Atheriniform Fishes

Atheriniform fishes have recently emerged as attractive models for evolutionary, ecological, and molecular/physiological studies on sex determination. Many species in this group have marked temperature-dependent sex determination (TSD) and yet many species also have a sex determinant gene that provides a strong drive for male differentiation. Thus, in these species the 2 forms of sex determination that were once considered to be mutually exclusive, environmental (ESD) and genotypic (GSD) sex determination, can coexist at environmentally relevant conditions. Here, we review the current knowledge on sex determination in atheriniform fishes with emphasis on the molecular and physiological mechanisms of ESD and GSD, the coexistence and cross-talk between these 2 mechanisms, the possibility of extragonadal transduction of environmental information and/or extragonadal onset of sex determination, and the results of field studies applying novel tools such as otolith increment analysis and molecular markers of genetic sex developed for selected New World and Old World atheriniform species. We also discuss the existence of molecular and histological mechanisms to prevent the discrepant differentiation in parts of the gonads because of ambiguous or conflicting environmental and genetic signals and particularly the possibility that the female is the default state in these species.

Early response to heat stress in Chinese tongue sole (Cynoglossus semilaevis): performance of different sexes, candidate genes and networks

BACKGROUND

Temperature is known to affect living organisms and alter the expression of responsive genes, which affects a series of life processes, such as development, reproduction and metabolism. Several genes and gene families have been involved in high temperature responses, such as heat shock protein (hsp) family, Jumonji family and genes related to cortisol synthesis. Gonad is a vital organ related to the existence of a species. However, the comprehensive understanding of gonadal responses to environmental temperature is limited.

RESULTS

To explore the effects of environmental temperature on genes and gene networks in gonads, we performed acute heat treatment (48 h) on Chinese tongue sole (Cynoglossus semilaevis). Gonadal transcriptome analysis was conducted on females, pseudomales and males exposed to high (28 °C) and normal (22 °C) temperatures. A total of 1226.24 million clean reads were obtained from 18 libraries. Principal component analysis (PCA) and differentially expressed gene (DEG) analysis revealed different performance of sex responses to heat stress. There were 4565, 790 and 1117 specific genes altered their expression level in females, pseudomales and males, respectively. Of these, genes related to hsp gene family, cortisol synthesis and metabolism and epigenetic regulation were involved in early heat response. Furthermore, a total of 1048 DEGs were shared among females, pesudomales and males, which may represent the inherent difference between high and normal temperatures. Genes, such as eef1akmt3, eef1akmt4, pnmt and hsp family members, were found.

CONCLUSIONS

Our results depicted for the first time the gonadal gene expression under acute high temperature treatment in Chinese tongue sole. The findings may provide a clue for understanding the responses of genes and networks to environmental temperature.

Transcriptome Profiling and Analysis of Genes Associated with High Temperature-Induced Masculinization in Sex-Undifferentiated Nile Tilapia Gonad

Artificially high temperatures during critical thermosensitive periods (TSPs) can induce the sex reversal of Nile tilapia (Oreochromis niloticus) females into pseudomales; Nile tilapia is a GSD + TE (genotypic plus temperature effects) fish species. Previous studies have shown that water temperature affects the expression levels of many genes in the gonad or brain in various teleost species. However, few studies on the effect of temperature at the whole-gonad transcriptomic level in the early stage of sex differentiation have been reported in fish species exhibiting GSD + TE. In this study, RNA-Seq was performed to characterize the transcriptomic profile and identify genes exhibiting temperature- and sex-biased expressions in the Nile tilapia gonad at 21 dpf. A total of 42 genes were found to be associated with both high-temperature treatment and sex development, as the expression levels of these genes differed in both FC (female control) vs MC (male control) and FC vs FT (high temperature-treated females in the TSP). Among these genes, the transcriptional alterations of many male sex determination and differentiation genes, such as Dmrt1, Gsdf, and the DNA damage-inducible protein GADD45 alpha, suggested that the male pathway is initiated after high-temperature treatment and that its initiation may play a role in high temperature-induced masculinization in Nile tilapia. The qRT-PCR validation results for thirteen differentially expressed genes showed that the Pearson's correlation of the log₁₀ fold change values between the qPCR and RNA-Seq results was 0.70 (p < 0.001), indicating the accuracy and reliability of the RNA-Seq results. Our study provides insights into how high-temperature treatment induces the sex reversal of Nile tilapia females.

Transcriptome of Gonads From High Temperature Induced Sex Reversal During Sex Determination and Differentiation in Chinese Tongue Sole, Cynoglossus semilaevis

The sex of Chinese tongue sole (Cynoglossus semilaevis) is determined by both genetic sex determination (GSD) and environmental sex determination (ESD), making it an ideal model to study the relationship between sex-determination and temperature. In the present study, transcriptomes of undifferentiated gonads from genetic females and males, as well as differentiated gonads from males, females, and pseudomales under high and normal temperature treatments were generated for comparative transcriptomic analysis. A mean of 68.24 M high-quality clean reads was obtained for each library. Differentially expressed genes (DEGs) between different sexes and environmental treatments were identified, revealing that the heat shock protein gene family was involved in the high temperature induced sex reversal. The Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were enriched in pseudomale and genetic female comparison included neuroactive ligand-receptor interaction, cortisol synthesis and secretion, and steroid hormone biosynthesis. Furthermore, weighted gene co-expression network analyses were conducted on all samples, and two modules were positive correlated with pseudomale under high temperature. An illustrated protein-protein interaction map of the module identified a hub gene, hsc70. These findings provide insights into the genetic network that is involved in sex determination and sexual differentiation, and improve our understanding of genes involved in sex reversal under high temperature.

Profiling expression changes of genes associated with temperature and sex during high temperature-induced masculinization in the Nile tilapia brain

Fish sex-determining mechanisms can be classified as genotypic (GSD), temperature (TSD), or genotypic plus temperature effects (GSD+TE). Previous studies have shown that culturing water temperature during thermosensitive periods (TSP) could affect the expression of many genes in the gonad in some fish. However, few studies have focused on gene expression changes in the brain after temperature treatment during TSP in fish species. In this study, three families were developed by crossing XX neomales with XX females and one of them was used for transcriptome analysis. The results showed that a total of 105, 3164 and 4666 DEGs were respectively obtained in FC (female control) vs. FT (high temperature-treated females at TSP), FC vs. MC (male control), and MC vs. FT comparison groups. By profiling analysis, we show that the mRNA expression levels of 16 differentially expressed genes (DEGs) exhibited significant downregulation or upregulation after high temperature treatment and reached a similar level as that in MC. Among the 16 DEGs, LOC100699848 (lysine specific demethylase 6A) and Jarid2 contained JmjC domain, showing the possible important role of JmjC domain in response to temperature treatment in Nile tilapia. Kdm6b (lysine demethylase 6B) and Jarid2 have been shown to play important roles in reptile TSD, showing the relative conservation of underlying regulation mechanisms between TSD in reptile and TSD or GSD+TE in fish species. Finally, the transcriptome profiling was validated by quantitative real-time PCR in nine selected genes. These results provide a direction for investigating the GSD+TE molecular mechanism in fish species.

Transcriptome profiling reveals male- and female-specific gene expression pattern and novel gene candidates for the control of sex determination and gonad development in Xenopus laevis

Xenopus laevis is an amphibian (frog) species widely used in developmental biology and genetics. To unravel the molecular machinery regulating sex differentiation of Xenopus gonads, we analyzed for the first time the transcriptome of developing amphibian gonads covering sex determination period. We applied microarray at four developmental stages: (i) NF50 (undifferentiated gonad during sex determination), (ii) NF53 (the onset of sexual differentiation of the gonads), (iii) NF56 (sexual differentiation of the gonads), and (iv) NF62 (developmental progression of differentiated gonads). Our analysis showed that during the NF50, the genetic female (ZW) gonads expressed more sex-specific genes than genetic male (ZZ) gonads, which suggests that a robust genetic program is realized during female sex determination in Xenopus. However, a contrasting expression pattern was observed at later stages (NF56 and NF62), when the ZW gonads expressed less sex-specific genes than ZZ gonads, i.e., more genes may be involved in further development of the male gonads (ZZ). We identified sexual dimorphism in the expression of several functional groups of genes, including signaling factors, proteases, protease inhibitors, transcription factors, extracellular matrix components, extracellular matrix enzymes, cell adhesion molecules, and epithelium-specific intermediate filaments. In addition, our analysis detected a sexually dimorphic expression of many uncharacterized genes of unknown function, which should be studied further to reveal their identity and if/how they regulate gonad development, sex determination, and sexual differentiation. Comparison between genes sex-specifically expressed in developing gonads of Xenopus and available transcriptome data from zebrafish, two reptile species, chicken, and mouse revealed significant differences in the genetic control of sex determination and gonad development. This shows that the genetic control of gonad development is evolutionarily malleable.

Genomic analysis of Sparus aurata reveals the evolutionary dynamics of sex-biased genes in a sequential hermaphrodite fish

Sexual dimorphism is a fascinating subject in evolutionary biology and mostly results from sex-biased expression of genes, which have been shown to evolve faster in gonochoristic species. We report here genome and sex-specific transcriptome sequencing of Sparus aurata, a sequential hermaphrodite fish. Evolutionary comparative analysis reveals that sex-biased genes in S. aurata are similar in number and function, but evolved following strikingly divergent patterns compared with gonochoristic species, showing overall slower rates because of stronger functional constraints. Fast evolution is observed only for highly ovary-biased genes due to female-specific patterns of selection that are related to the peculiar reproduction mode of S. aurata, first maturing as male, then as female. To our knowledge, these findings represent the first genome-wide analysis on sex-biased loci in a hermaphrodite vertebrate species, demonstrating how having two sexes in the same individual profoundly affects the fate of a large set of evolutionarily relevant genes.