Fish Gonadogenesis. Part I: Genetic and Environmental Mechanisms of Sex Determination

Application of machine learning in the study of development, behavior, nerve, and genotoxicity of zebrafish

Machine learning (ML) as a novel model-based approach has been used in studying aquatic toxicology in the environmental field. Zebrafish, as an ideal model organism in aquatic toxicology research, has been widely used to study the toxic effects of various pollutants. However, toxicity testing on organisms may cause significant harm, consume considerable time and resources, and raise ethical concerns. Therefore, ML is used in related research to reduce animal experiments and assist researchers in conducting toxicological research. Although ML techniques have matured in various fields, research on ML-based aquatic toxicology is still in its infancy due to the lack of comprehensive large-scale toxicity databases for environmental pollutants and model organisms. Therefore, to better understand the recent research progress of ML in studying the development, behavior, nerve, and genotoxicity of zebrafish, this review mainly focuses on using ML modeling to assess and predict the toxic effects of zebrafish exposure to different toxic chemicals. Meanwhile, the opportunities and challenges faced by ML in the field of toxicology were analyzed. Finally, suggestions and perspectives were proposed for the toxicity studies of ML on zebrafish in future applications.

Proteome analysis provides insights into sex differences in Holothuria Scabra

Sex-determining mechanism is still ambiguous for sea cucumber Holothuria scabra which only manifests gonochorism in gonad. In this study, proteomic analysis was employed to delineate sex-related proteins and genes in gonads of H. scabra, subsequently validated through Quantitative real-time polymerase chain reaction (qRT-PCR). A total of 5,313 proteins were identified via proteome sequencing. Among these, 817 proteins exhibited expression in both the ovary and testis, with 445 proteins displaying up-regulation and 372 proteins showing down-regulation (ovary vs testis). Furthermore, 136 and 69 proteins were identified as ovary-specific and testis-specific Differentially Abundant Proteins (DAPs), respectively. And 9 DAP coding genes which play crucial role in ovary and testis were verified by qRT-PCR. Notably, 24 ovary-bias proteins enriched in ribosome pathway strongly indicated the crucial role of ribosome in ovary. This study serves to furnish novel evidence pertaining to sex differences in H. scabra.

Temperature regulates sex determination and growth in the paralichthid flatfish California halibut

California halibut (Paralichthys californicus) is a candidate species for aquaculture and stock enhancement. These applications rely on sex control, either to maximize the production of faster growing females or to match sex ratios in the wild. Other paralichthids exhibit temperature-dependent sex determination (TSD), but the presence and pattern of TSD is not well defined in California halibut. Juvenile California halibut were cultured at three distinct temperatures (15°C, 19°C, and 23°C) through the developmental period presumed to be thermosensitive based on findings from congeners. Sex ratios were quantified in each treatment using phenotypic sex identification techniques applied early (molecular biomarkers; 51-100 mm total length [TL]) and late (visual examination of the gonads; ≥100 mm TL) in the juvenile phase. Both techniques indicated similar sex determination trends at each temperature, with overall sex ratios assessed as 49.9% male at 15°C, 74.5% male at 19°C, and 98.2% male at 23°C. Growth rates were highest at 23°C and lowest at 15°C, with intrinsically fast- and slow-growing individuals at all temperatures. At 15°C and 19°C, females comprised a higher proportion among the fast growers than they did among the slow growers. These data show that California halibut exhibit TSD, with temperatures of 19°C and 23°C masculinizing fish while 15°C appears to produce a 1:1 sex ratio. This study will help optimize sex ratios and growth in hatcheries through thermal manipulation. Furthermore, the developed biomolecular tools and identified temperature thresholds will be important in future work to understand the influence of global warming on wild population demographics.

The Direct Effects of Climate Change on Tench (Tinca tinca) Sperm Quality under a Real Heatwave Event Scenario

Global aquaculture growth will most probably face specific conditions derived from climate change. In fact, the most severe impacts of these changes will be suffered by aquatic populations in restrictive circumstances, such as current aquaculture locations, which represent a perfect model to study global warming effects. Although the impact of temperature on fish reproduction has been characterized in many aspects, this study was focused on recreating more realistic models of global warming, particularly considering heatwave phenomena, in order to decipher its effects on male gametes (spermatozoa). For this purpose, thermal stress via a heatwave simulation (mimicking a natural occurring heatwave, from 24 to 30 °C) was induced in adult tench (Tinca tinca) males and compared with a control group (55.02 ± 16.44 g of average body wet weight). The impact of the thermal stress induced by this climate change event was assessed using cellular and molecular approaches. After the heatwave recreation, a multiparametric analysis of sperm quality, including some traditional parameters (such as sperm motility) and new ones (focus on redox balance and sperm quality biomarkers), was performed. Although sperm concentration and the volume produced were not affected, the results showed a significant deleterious effect on motility parameters (e.g., reduced progressive motility and total motility during the first minute post-activation). Furthermore, the sperm produced under the thermal stress induced by this heatwave simulation exhibited an increased ROS content in spermatic cells, confirming the negative effect that this thermal stress model (heatwave recreation) might have had on sperm quality. More importantly, the expression of some known sperm quality and fertilization markers was decreased in males exposed to thermal stress. This present study not only unveils the potential effects of climate change in contemporary and future fish farming populations (and their underlying mechanisms) but also provides insights on how to mitigate and/or avoid thermal stress due to heatwave events.

Effect of high temperatures on sex ratio and differential expression analysis (RNA-seq) of sex-determining genes in Oreochromis niloticus from different river basins in Benin

The high temperature sex reversal process leading to functional phenotypic masculinization during development has been widely described in Nile tilapia (Oreochromis n iloticus) under laboratory or aquaculture conditions and in the wild. In this study, we selected five wild populations of O. niloticus from different river basins in Benin and produced twenty full-sib families of mixed-sex (XY and XX) by natural reproduction. Progenies were exposed to room temperature or high (36.5°C) temperatures between 10 and 30 days post-fertilization (dpf). In control groups, we observed sex ratios from 40% to 60% males as expected, except for 3 families from the Gobé region which showed a bias towards males. High temperature treatment significantly increased male rates in each family up to 88%. Transcriptome analysis was performed by RNA-sequencing (RNA-seq) on brains and gonads from control and treated batches of six families at 15 dpf and 40 dpf. Analysis of differentially expressed genes, differentially spliced genes, and correlations with sex reversal was performed. In 40 dpf gonads, genes involved in sex determination such as dmrt1, cyp11c1, amh, cyp19a1b, ara, and dax1 were upregulated. In 15 dpf brains, a negative correlation was found between the expression of cyp19a1b and the reversal rate, while at 40 dpf a negative correlation was found between the expression of foxl2, cyp11c1, and sf1 and positive correlation was found between dmrt1 expression and reversal rate. Ontology analysis of the genes affected by high temperatures revealed that male sex differentiation processes, primary male sexual characteristics, autophagy, and cilium organization were affected. Based on these results, we conclude that sex reversal by high temperature treatment leads to similar modifications of the transcriptomes in the gonads and brains in offspring of different natural populations of Nile tilapia, which thus may activate a common cascade of reactions inducing sex reversal in progenies.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Genome-wide SNP analysis provides insights into the XX/XY sex-determination system in silver barb (Barbonymus gonionotus)

Silver barb (Barbonymus gonionotus) is among the most economically important freshwater fish species in Thailand. It ranks fourth in economic value and third in production weight for fisheries and culture in Thailand. An XX/XY sex-determination system based on gynogenesis was previously reported for this fish. In this study, the molecular basis underlying the sex-determination system was further investigated. Genome-wide single-nucleotide polymorphism data were generated for 32 captive-bred silver barb individuals, previously scored by phenotypic sex, to identify sex-linked regions associated with sex determination. Sixty-three male-linked loci, indicating putative XY chromosomes, were identified. Male-specific loci were not observed, which indicates that the putative Y chromosome is young and the sex determination region is cryptic. A homology search revealed that most male-linked loci were homologous to the Mariner/Tc1 and Gypsy transposable elements and are probably the remnants of an initial accumulation of repeats on the Y chromosome from the early stages of sex chromosome differentiation. This research provides convincing insights into the mechanism of sex determination and reveals the potential sex determination regions in silver barb. The study provides the basic data necessary for increasing the commercial value of silver barbs through genetic improvements.

Validation of a male-specific DNA marker confirms XX/XY-type sex determination in several Hungarian strains of African catfish (Clarias gariepinus)

African catfish (Clarias gariepinus) is a promising food fish species with significant potential and growing mass of production in freshwater aquaculture. Male African catfish possess improved production characteristics over females, therefore the use of monosex populations could be advantageous for aquaculture production. However, our knowledge about the sex determination mechanism of this species is still limited and controversial. A previously isolated male-specific DNA marker (CgaY1) was validated using offspring groups from targeted crosses (n = 630) and it was found to predict the sex of 608 individuals correctly (96.43% accuracy). Using the proportion of recombinants, we estimated the average genetic distance between the potential sex determination locus and the sex-specific marker to be 3.57 cM. As an earlier study suggested that both XX/XY and ZZ/ZW systems coexist in this species, we tested the applicability of their putative 'moderately sex-linked loci' and found that no sex-specific amplification could be detected for any of them. In addition, temperature-induced masculinization suggested by others was also tested, but no such effect was detected in our stocks when the published parameters were used for heat treatment. Altogether, our results support an exclusive XX/XY sex determination system in our African catfish stock and indicate a good potential for the future use of this male-specific DNA marker in research and commercial production.

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

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

Use of male-to-female sex reversal as a welfare scoring system in the protandrous farmed gilthead sea bream (Sparus aurata)

Gilthead sea bream is a highly cultured marine fish throughout the Mediterranean area, but new and strict criteria of welfare are needed to assure that the intensification of production has no negative effects on animal farming. Most welfare indicators are specific to a given phase of the production cycle, but others such as the timing of puberty and/or sex reversal are of retrospective value. This is of particular relevance in the protandrous gilthead sea bream, in which the sex ratio is highly regulated at the nutritional level. Social and environmental factors (e.g., contaminant loads) also alter the sex ratio, but the contribution of the genetic component remains unclear. To assess this complex issue, five gilthead sea bream families representative of slow/intermediate/fast growth were grown out with control or a plant-based diet in a common garden system from early life to the completion of their sexual maturity in 3-year-old fish. The plant-based diet highly enhanced the male-to-female sex reversal. This occurred in parallel with the progressive impairment of growth performance, which was indicative of changes in nutrient requirements as the result of the different energy demands for growth and reproduction through development. The effect of a different nutritional and genetic background on the reproductive performance was also assessed by measurements of circulating levels of sex steroids during the two consecutive spawning seasons, varying plasma levels of 17β-estradiol (E₂) and 11-ketotestosterone (11-KT) with age, gender, diet, and genetic background. Principal component analysis (PCA) of 3-year-old fish displayed a gradual increase of the E₂/11-KT ratio from males to females with the improvement of nutritional/genetic background. Altogether, these results support the use of a reproductive tract scoring system for leading farmed fish toward their optimum welfare condition, contributing to improving the productivity of the current gilthead sea bream livestock.

Epigenetic Regulation of Phenotypic Sexual Plasticity Inducing Skewed Sex Ratio in Zebrafish

The plasticity of sexual phenotype in response to environmental conditions results in biased sex ratios, and their variation has an effect on population dynamics. Epigenetic modifications can modulate sex ratio variation in species, where sex is determined by genetic and environmental factors. However, the role of epigenetic mechanisms underlying skewed sex ratios is far from being clear and is still an object of debate in evolutionary developmental biology. In this study, we used zebrafish as a model animal to investigate the effect of DNA methylation on sex ratio variation in sex-biased families in response to environmental temperature. Two sex-biased families with a significant difference in sex ratio were selected for genome-wide DNA methylation analysis using reduced representation bisulfite sequencing (RRBS). The results showed significant genome-wide methylation differences between male-biased and female-biased families, with a greater number of methylated CpG sites in testes than ovaries. Likewise, pronounced differences between testes and ovaries were identified within both families, where the male-biased family exhibited a higher number of methylated sites than the female-biased family. The effect of temperature showed more methylated positions in the high incubation temperature than the control temperature. We found differential methylation of many reproduction-related genes (e.g., sox9a, nr5a2, lhx8a, gata4) and genes involved in epigenetic mechanisms (e.g., dnmt3bb.1, dimt1l, hdac11, h1m) in both families. We conclude that epigenetic modifications can influence the sex ratio variation in zebrafish families and may generate skewed sex ratios, which could have a negative consequence for population fitness in species with genotype-environment interaction sex-determining system under rapid environmental changes.

Sex-specific differences in zebrafish brains

In this systematic review, we highlight the differences between the male and female zebrafish brains to understand their differentiation and their use in studying sex-specific neurological diseases. Male and female brains display subtle differences at the cellular level which may be important in driving sex-specific signaling. Sex differences in the brain have been observed in humans as well as in non-human species. However, the molecular mechanisms of brain sex differentiation remain unclear. The classical model of brain sex differentiation suggests that the steroid hormones derived from the gonads are the primary determinants in establishing male and female neural networks. Recent studies indicate that the developing brain shows sex-specific differences in gene expression prior to gonadal hormone action. Hence, genetic differences may also be responsible for differentiating the brain into male and female types. Understanding the signaling mechanisms involved in brain sex differentiation could help further elucidate the sex-specific incidences of certain neurological diseases. The zebrafish model could be appropriate for enhancing our understanding of brain sex differentiation and the signaling involved in neurological diseases. Zebrafish brains show sex-specific differences at the hormonal level, and recent advances in RNA sequencing have highlighted critical sex-specific differences at the transcript level. The differences are also evident at the cellular and metabolite levels, which could be important in organizing sex-specific neuronal signaling. Furthermore, in addition to having one ortholog for 70% of the human gene, zebrafish also shares brain structural similarities with other higher eukaryotes, including mammals. Hence, deciphering brain sex differentiation in zebrafish will help further enhance the diagnostic and pharmacological intervention of neurological diseases.

Switches, stability and reversals in the evolutionary history of sexual systems in fish

Sexual systems are highly diverse and have profound consequences for population dynamics and resilience. Yet, little is known about how they evolved. Using phylogenetic Bayesian modelling and a sample of 4614 species, we show that gonochorism is the likely ancestral condition in teleost fish. While all hermaphroditic forms revert quickly to gonochorism, protogyny and simultaneous hermaphroditism are evolutionarily more stable than protandry. In line with theoretical expectations, simultaneous hermaphroditism does not evolve directly from gonochorism but can evolve slowly from sequential hermaphroditism, particularly protandry. We find support for the predictions from life history theory that protogynous, but not protandrous, species live longer than gonochoristic species and invest the least in male gonad mass. The distribution of teleosts' sexual systems on the tree of life does not seem to reflect just adaptive predictions, suggesting that adaptations alone may not fully explain why some sexual forms evolve in some taxa but not others (Williams' paradox). We propose that future studies should incorporate mating systems, spawning behaviours, and the diversity of sex determining mechanisms. Some of the latter might constrain the evolution of hermaphroditism, while the non-duality of the embryological origin of teleost gonads might explain why protogyny predominates over protandry in teleosts.

Proteome and Transcriptome Analysis of Gonads Reveals Intersex in Gigantidas haimaensis

Sex has proven to be one of the most intriguing areas of research across evolution, development, and ecology. Intersex or sex change occurs frequently in molluscs. The deep-sea mussel Gigantidas haimaensis often dominates within Haima cold seep ecosystems, but details of their reproduction remain unknown. Herein, we conducted a combined proteomic and transcriptomic analysis of G. haimaensis gonads to provide a systematic understanding of sexual development in deep-sea bivalves. A total of 2,452 out of 42,238 genes (5.81%) and 288 out of 7,089 proteins (4.06%) were significantly differentially expressed between ovaries and testes with a false discovery rate (FDR) <0.05. Candidate genes involved in sexual development were identified; among 12 differentially expressed genes between sexes, four ovary-biased genes (β-catenin, fem-1, forkhead box L2 and membrane progestin receptor α) were expressed significantly higher in males than females. Combining histological characteristics, we speculate that the males maybe intersex undergoing sex change, and implied that these genes may be involved in the process of male testis converting into female gonads in G. haimaensis. The results suggest that this adaptation may be based on local environmental factors, sedentary lifestyles, and patchy distribution, and sex change may facilitate adaptation to a changing environment and expansion of the population. The findings provide a valuable genetic resource to better understand the mechanisms of sex change and survival strategies in deep-sea bivalves.

Loss of Nobox prevents ovarian differentiation from juvenile ovaries in zebrafish

As a species without master sex-determining genes, zebrafish displays high plasticity in sex differentiation, making it an excellent model for studying the regulatory mechanisms underlying gonadal differentiation and gametogenesis. Despite being a gonochorist, zebrafish is a juvenile hermaphrodite that undergoes a special phase of juvenile ovary before further differentiation into functional testis and ovary. The mechanisms underlying juvenile ovary formation and subsequent gonadal differentiation remain largely unknown. In a recent study, we demonstrated an important role for Figla (factor in the germline alpha) in zebrafish oogenesis. In this study, we explored the role of Nobox/nobox (new born ovary homeobox protein), another oocyte-specific transcription factor in females, in early zebrafish gonadogenesis using CRISPR/Cas9 technology. As in mammals, nobox is specifically expressed in zebrafish gonads with a dimorphic pattern at juvenile stage. In contrast to the mutant of figla (another oocyte-specific transcription factor), the nobox mutants showed formation of typical perinucleolar (PN) follicles at primary growth (PG) stage in juvenile gonads, suggesting occurrence of follicle assembly from cystic oocytes (chromatin nucleolar stage, CN). These follicles, however, failed to develop further to form functional ovaries, resulting in all-male phenotype. Despite its expression in adult testis, the loss of nobox did not seem to affect testis development, spermatogenesis and male spawning. In summary, our results indicate an important role for Nobox in zebrafish ovarian differentiation and early folliculogenesis.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

Advances in Reproductive Endocrinology and Neuroendocrine Research Using Catfish Models

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

The Snakeskin Gourami (Trichopodus pectoralis) Tends to Exhibit XX/XY Sex Determination

The snakeskin gourami (Trichopodus pectoralis) has a high meat yield and is one of the top five aquaculture freshwater fishes in Thailand. The species is not externally sexually dimorphic, and its sex determination system is unknown. Understanding the sex determination system of this species will contribute to its full-scale commercialization. In this study, a cytogenetic analysis did not reveal any between-sex differences in chromosomal patterns. However, we used genotyping-by-sequencing to identify 4 male-linked loci and 1 female-linked locus, indicating that the snakeskin gourami tends to exhibit an XX/XY sex determination system. However, we did not find any male-specific loci after filtering the loci for a ratio of 100:0 ratio of males:females. This suggests that the putative Y chromosome is young and that the sex determination region is cryptic. This approach provides solid information that can help identify the sex determination mechanism and potential sex determination regions in the snakeskin gourami, allowing further investigation of genetic improvements in the species.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. 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)'.

A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus)

BACKGROUND

Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches.

RESULTS

The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD.

CONCLUSIONS

The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

Epigenetic differences in the innate response after immune stimulation during zebrafish sex differentiation

Infections are able to trigger epigenetic modifications; however, epigenetic-mediating infections in the immune system in fish is currently unavailable. Within this purpose, zebrafish were immune-stimulated with three lipopolysaccharides (LPS) during sex differentiation. Methylation patterns of three immune genes were studied by a candidate gene approach together with gene expression analysis, and in adulthood, sex ratios were determined. It was shown that the entrance of LPS was through the gills and accumulated in the pronephros. Significant hypomethylation levels of CASP9 and a significant CpG site for IL1β after Pseudomonas aeruginosa LPS exposure were found. No methylation difference was observed for TNFα. Gene expression and correlation data differed among studied genes. Sex ratios showed a feminization in dose and LPS strain-dependent manner. Here, it is provided epigenetic regulatory mechanisms derived by innate response and the first evidence of possible epigenetic interactions between the immune and reproductive systems.

Activation of stress response axis as a key process in environment-induced sex plasticity in fish

The determination of sex is an important hallmark in the life cycle of organisms, in which the fate of gonads and then the individual sex are defined. In gonochoristic teleost fish, this process is characterized by a high plasticity, considering that in spite of genotypic sex many environmental factors can cause shifts from one to another molecular pathway, resulting in organisms with mismatching genotypic and phenotypic sexes. Interestingly, in most instances, both female-to-male or male-to-female sex-reversed individuals develop functional gonads with normal gametogenesis and respective progenies with full viability. The study of these mechanisms is being spread to other non-model species or to those inhabiting more extreme environmental conditions. Although water temperature is an important mechanism involved in sex determination, there are other environmental stressors affected by the climate change which are also implicated in stress response-induced masculinization in fish. In this regard, the brain has emerged as the transducer of the environment input that can influence the gonadal fate. Furthermore, the evaluation of other environmental stressors or their synergic effect on sex determination at conditions that simulate the natural environments is growing gradually. Within such scope, the concerns related to climate change impacts rely on the fact that many of biotic and abiotic parameters reported to affect sex ratios are expected to increase concomitantly as a result of increased greenhouse gas emissions and, particularly worrying, many of them are related to male bias in the populations, such as high temperature, hypoxia, and acidity. These environmental changes can also generate epigenetic changes in sex-related genes affecting their expression, with implications on sex differentiation not only of exposed individuals but also in following generations. The co-analysis of multi-stressors with potential inter- and transgenerational effects is essential to allow researchers to perform long-term predictions on climate change impacts in wild populations and for establishing highly accurate monitoring tools and suitable mitigation strategies.

Are long-term growth responses to elevated pCO2 sex-specific in fish?

Whether marine fish will grow differently in future high pCO₂ environments remains surprisingly uncertain. Long-term and whole-life cycle effects are particularly unknown, because such experiments are logistically challenging, space demanding, exclude long-lived species, and require controlled, restricted feeding regimes—otherwise increased consumption could mask potential growth effects. Here, we report on repeated, long-term, food-controlled experiments to rear large populations (>4,000 individuals total) of the experimental model and ecologically important forage fish Menidia menidia (Atlantic silverside) under contrasting temperature (17°, 24°, and 28°C) and pCO₂ conditions (450 vs. ~2,200 μatm) from fertilization to ~ a third of this annual species’ life span. Quantile analyses of trait distributions showed mostly negative effects of high pCO₂ on long-term growth. At 17°C and 28°C, but not at 24°C, high pCO₂ fish were significantly shorter [17°C: -5 to -9%; 28°C: -3%] and weighed less [17°C: -6 to -18%; 28°C: -8%] compared to ambient pCO₂ fish. Reductions in fish weight were smaller than in length, which is why high pCO₂ fish at 17°C consistently exhibited a higher Fulton’s k (weight/length ratio). Notably, it took more than 100 days of rearing for statistically significant length differences to emerge between treatment populations, showing that cumulative, long-term CO₂ effects could exist elsewhere but are easily missed by short experiments. Long-term rearing had another benefit: it allowed sexing the surviving fish, thereby enabling rare sex-specific analyses of trait distributions under contrasting CO₂ environments. We found that female silversides grew faster than males, but there was no interaction between CO₂ and sex, indicating that males and females were similarly affected by high pCO₂. Because Atlantic silversides are known to exhibit temperature-dependent sex determination, we also analyzed sex ratios, revealing no evidence for CO₂-dependent sex determination in this species.

Genome-wide association analyses based on whole-genome sequencing of Protosalanx hyalocranius provide insights into sex determination of Salangid fishes

Identification of sex determination system and sex-determining genes have important implications in conservation, ecology and evolution. However, much remains to be discovered about the evolution of different sexual determination systems in teleost fishes, of which the mechanisms of sex determination are remarkably variable. In the present study, the whole genomes of 20 males and 20 females of a Salangid fish, Protosalanx hyalocranius, were sequenced and genome wide association analyses were conducted to uncover its sex determination system and putative sex-determining genes. A total of 150 SNPs were significantly associated with sex, which showed high differentiation between sexes (F_ST ranged from 0.245 to 0.556). Of the 150 sex-associated SNPs, 76 SNPs displayed sex specificity with even coverage of depth and were female heterogametic, which suggested a ZZ/ZW sex determination system. Interestingly, one scaffold containing sex-specific SNPs displayed synteny to the sex chromosome of medaka. Annotations of sex-associated loci suggested that both transcriptional regulators (e.g., FOX genes) and secreted hormones and their receptors might be involved in the sex determination/differentiation of P. hyalocranius. More strikingly, we found a nonsense mutation in one copy of GALNT homology gene of all females, which suggested that "Z dosage" effect might play a vital role in the processes of sex determination/differentiation. These sex-specific loci could be a valuable resource for further research on sex determination of Salangid fishes and the results could contribute to the understanding of sex determination mechanisms and the evolution of sex chromosome in teleost fishes.

Climate change impacts on fish reproduction are mediated at multiple levels of the brain-pituitary-gonad axis

Anthropogenic emissions of carbon dioxide in the atmosphere have generated rapid variations in atmospheric composition which drives major climate changes. Climate change related effects include changes in physico-chemical proprieties of sea and freshwater, such as variations in water temperature, salinity, pH/pCO₂ and oxygen content, which can impact fish critical physiological functions including reproduction. In this context, the main aim of the present review is to discuss how climate change related effects (variation in water temperature and salinity, increases in duration and frequency of hypoxia events, water acidification) would impact reproduction by affecting the neuroendocrine axis (brain-pituitary-gonad axis). Variations in temperature and photoperiod regimes are known to strongly affect sex differentiation and the timing and phenology of spawning period in several fish species. Temperature mainly acts at the level of gonad by interfering with steroidogenesis, (notably on gonadal aromatase activity) and gametogenesis. Temperature is also directly involved in the quality of released gametes and embryos development. Changes in salinity or water acidification are especially associated with reduction of sperm quality and reproductive output. Hypoxia events are able to interact with gonad steroidogenesis by acting on the steroids precursor cholesterol availability or directly on aromatase action, with an impact on the quality of gametes and reproductive success. Climate change related effects on water parameters likely influence also the reproductive behavior of fish. Although the precise mechanisms underlying the regulation of these effects are not always understood, in this review we discuss different hypothesis and propose future research perspectives.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23

Background

Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins.

Results

A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284).

Conclusions

Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

2,4-Dichlorophenol induces feminization of zebrafish (Danio rerio) via DNA methylation

2,4-Dichlorophenol (2,4-DCP) is a ubiquitous contaminant of aquatic environments with an estrogenic effect on fish. However, the molecular mechanism underlying this effect remains elusive. To this end, the present study aimed to explore the effect of 2,4-DCP on sex differentiation and its relevant mechanism in zebrafish (Danio rerio). The results showed that a female-biased sex ratio was induced after exposing larval zebrafish to 2,4-DCP (0-160 μg/L) from 20 to 50 days post fertilization (dpf). The feminization of zebrafish was accompanied by decreased expression of male-related genes (sox9a, amh and dmrt1) under 2,4-DCP from 20 to 50 dpf. However, the expression of female-related genes (cyp19a1a, foxl2 and esr1) was also suppressed. Nevertheless, it is noteworthy that the methylation level of sox9a promoter was significantly increased, which may result in the significantly decreased expression of sox9a and ultimately the feminization effect of 2,4-DCP on zebrafish. In addition, 5-aza-2'-deoxycytidine (5-AZA), a methyltransferase inhibitor, significantly reduced the methylation level, increased the expression of sox9a, and partly impaired the feminization effect caused by 2,4-DCP, which further confirmed the importance of DNA methylation of sox9a in 2,4-DCP-induced feminization. These findings provide novel insights into the epigenetic mechanisms of DCP-induced estrogenic effect in fish.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Analyses of mRNA-seq and miRNA-seq of the brain reveal the sex differences of gene expression and regulation before and during gonadal differentiation in 17β-estradiol or 17α-methyltestosterone-induced olive flounder (Paralichthys olivaceus)

Olive flounder (Paralichthys olivaceus) is a commercially important flatfish species cultured in East Asia. Female flounders generally grow more rapidly than males, therefore control of the sex ratio seems to be a proposed way to increase production. However, the sex determination gene and sex determination mechanism have yet been elucidated. The brain is an important organ that is involved in gonadal development. To explore the sex differences of gene expression in the brain before and during the flounder gonadal differentiation, we used messenger RNA (mRNA)-seq technology to investigate transcriptomes of male and female brains. Between female and male brains, 103 genes were differentially expressed before ovarian differentiation, 16 genes were differentially expressed before testicular differentiation, and 64 genes were differentially expressed during gonadal differentiation. According to annotation and Kyoto Encyclopedia of Genes and Genomes information, the differentially expressed genes (DEGs) were involved in circadian rhythm, circadian rhythm-fly, circadian entrainment, dopaminergic synapse, calcium signaling, glutamatergic synapse, taste transduction, herpes simplex infection, long-term depression, retrograde endocannabinoid signaling, and the synaptic vesicle cycle pathways. MicroRNA (miRNA)-seq was performed during the gonadal differentiation and the target genes of miRNAs were predicted. Integrated analysis of mRNA-seq and miRNA-seq showed that 29 of the 64 DEGs were regulated by the differentially expressed miRNAs during the gonadal differentiation. Our study provides a basis for further studies of brain sex differentiation and the molecular mechanism of sex determination in olive flounder.

Changes in DNA methylation during epigenetic-associated sex reversal under low temperature in Takifugu rubripes

DNA methylation has frequently been implicated in sex determination and differentiation in teleost species. In order to detect the DNA methylation patterns established during sexual differentiation in tiger pufferfish T. rubripes, we performed comprehensive whole genome methylation sequencing and analyses of the gonads of male, female, and pseudo male. We obtained a total of 33.12, 32.44, and 31.60 Gb clean data for male, female, and pseudo male, with a sequencing depth of 66.44×, 60.47× and 54.86×, respectively. The methylation level of cytosine (C) residues in the genomic DNA from gonads was 11.016%, 10.428%, and 11.083% in male, female, and pseudo male, respectively. More than 65% of C methylation was at CpG sites, and less than 1% was at CHG and CHH sites. In each regulatory element, there were low methylation levels on both sides of the transcription start site, and higher methylation levels in exons, introns, and downstream of genes. The highest mCpG was on chromosome 8 and the lowest mCpG was on chromosome 5. Comparisons of whole-genome DNA methylation between pairs of samples revealed that there were 3,173 differentially methylated regions (DMRs) between female and male, and 3,037 DMRs between male and pseudo male, corresponding to 0.232% and 0.223% of the length of the genome, respectively. There were only 1,635 DMRs between female and pseudo male, representing 0.127% of the length of the genome. A number of differentially methylated genes (DMGs) related to sex determination and differentiation were selected, such as amhr2 and pfcyp19a. After Bisulfite Sequencing PCR (BSP) verification, amhr2 was exhibited low methylation level in normal males and pseudo male, and high methylation level in normal females but pfcyp19a showed low methylation level in normal females and high methylation level in normal males and pseudo males. These results provide information about the molecular epigenetic mechanisms of DNA methylation during low-temperature induced masculinization of tiger pufferfish, and increase our understanding of the mechanisms of sex determination and differentiation in this important aquaculture fish species.

Inhibition of ovarian development and instances of sex reversal in genotypic female sablefish (Anoplopoma fimbria) exposed to elevated water temperature

This study determined high temperature effects on ovarian development in a marine groundfish species, sablefish (Anoplopoma fimbria), with potential application in sex reversal or sterilization for aquaculture. Monosex female (XX-genotype) sablefish larvae (∼30 mm) were randomly divided into three groups and exposed to control (15.6 °C ± 0.8 °C), moderate (20.4 °C ± 0.5 °C), or high (21.7 °C ± 0.5 °C) temperatures for 19 weeks. Treated fish were then tagged and transferred to ambient seawater (11.2 °C ± 2.3 °C) for one year to determine whether temperature effects on reproductive development were maintained post-treatment. Fish were periodically sampled for gonadal histology, gene expression and plasma 17β-estradiol (E2) analyses to assess gonadal development. Short-term (4-week) exposure to elevated temperatures had only minor effects, whereas longer exposure (12-19 weeks) markedly inhibited ovarian development. Fish from the moderate and high treatment groups had significantly less developed ovaries relative to controls, and mRNA levels for germ cell (vasa, zpc) and apoptosis-associated genes (p53, casp8) generally indicated gonadal degeneration. The high treatment group also had significantly reduced plasma E2 levels and elevated gonadal amh gene expression. After one year at ambient temperatures, however, ovaries of moderate and high treatment fish exhibited compensatory recovery and were indistinguishable from controls. Two genotypic females possessing immature testes (neomales) were observed in the high treatment group, indicating sex reversal had occurred (6% rate). These results demonstrate that extreme elevated temperatures may inhibit ovarian development or trigger sex reversal. High temperature treatment is likely not an effective sterilization method but may be preferable for sablefish neomale broodstock production.

Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome

The accumulation of repetitive DNA sequences on the sex-limited W or Y chromosomes is a well-known process that is likely triggered by the suppression of recombination between the sex chromosomes, which leads to major differences in their sizes and genetic content. Here, we report an analysis conducted on the satellitome of Megaleporinus macrocephalus that focuses specifically on the satDNAs that have been shown to have higher abundances in females and are putatively located on the W chromosome in this species. We characterized 164 satellite families in M. macrocephalus, which is, by far, the most satellite-rich species discovered to date. Subsequently, we mapped 30 satellites, 22 of which were located on the W chromosome, and 14 were shown to exist only on the W chromosome. Finally, we report two simple, quick and reliable methods that can be used for sex identification in M. macrocephalus individuals using fin clips or scales, which could be applicable to future studies conducted in the field of aquaculture.

Sex determination in Neotropical fish: Implications ranging from aquaculture technology to ecological assessment

The high biodiversity of fish in the Neotropical region contrasts with scarce or biased studies on the mechanisms involved in the sex determination in members of this fauna. In this review, we attempted to compile the information available on determination, differentiation, and manipulation of sex for Neotropical species, with special focus on silversides and other two speciose groups, known as characins (Characiformes) and catfishes (Siluriformes). Currently, there is plenty of information available on chromosomal sex determination systems, which includes both male and female heterogamety with many variations, and sex chromosomes evolution at the macro chromosomal level. However, there is hitherto a blank in information at micro, gene/molecule levels and in research related to the effects of environmental cues on sex determination; most of reported studies are limited to silversides and guppies. In view of such a high diversity, it is critically necessary to establish key model species for relevant Neotropical fish taxa and also multi-disciplinary research groups in order to uncover the main patterns and trends that dictate the mechanisms of sex determination and gonadal differentiation in this icthyofauna. By increasing our knowledge on sex determination/differentiation with the identification of sex chromosome-linked markers or sex-determining genes, characterization of the onset timing of morphological gonadal differentiation, and determination of the environmental-hormonal labile period of gonadal sex determination in reference species, it will be possible to use those information as guidelines for application in other related groups. Overall, the strategic advance in this research field will be crucial for the development of biotechnological tools for aquaculture industry and for conservation of fish fauna from the Neotropical Region.

Characterization of a male specific region containing a candidate sex determining gene in Atlantic cod

The genetic mechanisms determining sex in teleost fishes are highly variable and the master sex determining gene has only been identified in few species. Here we characterize a male-specific region of 9 kb on linkage group 11 in Atlantic cod (Gadus morhua) harboring a single gene named zkY for zinc knuckle on the Y chromosome. Diagnostic PCR test of phenotypically sexed males and females confirm the sex-specific nature of the Y-sequence. We identified twelve highly similar autosomal gene copies of zkY, of which eight code for proteins containing the zinc knuckle motif. 3D modeling suggests that the amino acid changes observed in six copies might influence the putative RNA-binding specificity. Cod zkY and the autosomal proteins zk1 and zk2 possess an identical zinc knuckle structure, but only the Y-specific gene zkY was expressed at high levels in the developing larvae before the onset of sex differentiation. Collectively these data suggest zkY as a candidate master masculinization gene in Atlantic cod. PCR amplification of Y-sequences in Arctic cod (Arctogadus glacialis) and Greenland cod (Gadus macrocephalus ogac) suggests that the male-specific region emerged in codfishes more than 7.5 million years ago.

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.

Roles of Figla/figla in Juvenile Ovary Development and Follicle Formation During Zebrafish Gonadogenesis

Sex determination and differentiation are complex processes. As a juvenile hermaphrodite or undifferentiated gonochorist, zebrafish undergo a special juvenile ovarian phase during sex differentiation, making it an excellent model for studying early oogenesis and folliculogenesis. We provide lines of evidence at morphological, molecular, and genetic levels for roles of factor in the germline α (Figla), an oocyte-specific transcription factor, in early zebrafish gonadogenesis. As in mammals, Figla/figla was also expressed in the gonads and its expression in the ovary was also restricted to early oocytes. Disruption of figla gene by CRISPR/Cas9 led to an all-male phenotype in the mutant. Detailed analysis of early gonadal development showed that the germ cells in the mutant were clustered in cysts and underwent meiosis, forming oocytes at prefollicular chromatin nucleolar (CN) stage (stage IA). However, the subsequent transition from cystic CN oocytes to individual follicular perinucleolar oocytes (stage IB) was blocked, resulting in an all-male phenotype in the mutant. The phenotype of figla mutant could not be rescued by estrogen treatment, in contrast to cyp19a1a mutant, and introduction of tp53 mutation also had no effect, unlike in fancd1 and fancl mutants. Transcriptome analysis revealed that many biological processes and pathways related to germ cell development, especially oogenesis, were upregulated in the presence of Figla and that the regulation of figla expression may involve heat shock proteins. Our results strongly suggest important roles for Figla in juvenile ovary development, especially the formation of individual follicles from cystic oocytes.

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

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

FOXL2 and DMRT1L Are Yin and Yang Genes for Determining Timing of Sex Differentiation in the Bivalve Mollusk Patinopecten yessoensis

Sex determination and differentiation have long been a research hotspot in metazoans. However, little is known about when and how sex differentiation occurs in most mollusks. In this study, we conducted a combined morphological and molecular study on sex differentiation in the Yesso scallop Patinopecten yessoensis. Histological examination on gonads from 5- to 13-month-old juveniles revealed that the morphological sex differentiation occurred at 10 months of age. To determine the onset of molecular sex differentiation, molecular markers were screened for early identification of sex. The gonadal expression profiles of eight candidate genes for sex determination or differentiation showed that only two genes displayed sexually dimorphic expression, with FOXL2 being abundant in ovaries and DMRT1L in testes. In situ hybridization revealed that both of them were detected in germ cells and follicle cells. We therefore developed LOG₁₀(DMRT1L/FOXL2) for scallop sex identification and confirmed its feasibility in differentiated individuals. By tracing its changes in 5- to 13-month-old juveniles, molecular sex differentiation time was determined: some scallops differentiate early in September when they are 7 months old, and some do late in December when they are 10 months old. Two kinds of coexpression patterns were found between FOXL2 and DMRT1L: expected antagonism after differentiation and unexpected coordination before differentiation. Our results revealed that scallop sex differentiation co-occurs with the formation of follicles, and molecular sex differentiation is established prior to morphological sex differentiation. Our study will assist in a better understanding of the molecular mechanism underlying bivalve sex differentiation.

Expression analysis of growth differentiation factor 9 (Gdf9/gdf9), anti-müllerian hormone (Amh/amh) and aromatase (Cyp19a1a/cyp19a1a) during gonadal differentiation of the zebrafish, Danio rerio

In the zebrafish, no sex-determining gene has been identified, while some sex-related genes, such as cyp19a1a and amh, show sexually dimorphic expression. Interestingly, most of these genes are expressed in the somatic cells. With increasing evidence suggesting roles of germ cells in gonadal differentiation, there is an increasing interest in the factors released by the germ cells for the bidirectional communication between the two compartments. We have reported that Gdf9/gdf9 is an oocyte-specific factor in the zebrafish, similar to that of mammals. Whether and how Gdf9 is involved in gonadal differentiation is unknown. In this study, we compared the expression levels of gdf9, cyp19a1a, and amh among several other sex-related genes in the gonads before, during, and after sex differentiation. The expression of gdf9 started in the gonads before sex differentiation, and its level surged in the differentiated ovary. Its expression pattern was similar to that of cyp19a1a, but reciprocal to amh expression. Using recombinant zebrafish Gdf9 (rzfGdf9), we further showed that Gdf9 significantly suppressed the expression of amh while increased that of activin beta subunits (inhbaa and inhbb) in vitro. Although gdf9 and cyp19a1a showed co-expression during gonadal differentiation, we only observed a slight but not significant response of cyp19a1a to rzfGdf9. Knocking down the expression of gdf9 and cyp19a1a with vivo-morpholinos caused a male-skewed sex ratio. Our data suggested that Gdf9 is likely involved in promoting oocyte/ovary differentiation in the zebrafish and it may act by suppressing amh expression, at least partly, in the somatic cells.

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.

Proteome and Transcriptome Analysis of Ovary, Intersex Gonads, and Testis Reveals Potential Key Sex Reversal/Differentiation Genes and Mechanism in Scallop Chlamys nobilis

Bivalve mollusks exhibit hermaphroditism and sex reversal/differentiation. Studies generally focus on transcriptional profiling and specific genes related to sex determination and differentiation. Few studies on sex reversal/differentiation have been reported. A combination analysis of gonad proteomics and transcriptomics was conducted on Chlamys nobilis to provide a systematic understanding of sex reversal/differentiation in bivalves. We obtained 4258 unique peptides and 93,731 unigenes with good correlation between messenger RNA and protein levels. Candidate genes in sex reversal/differentiation were found: 15 genes differentially expressed between sexes were identified and 12 had obvious sexual functions. Three novel genes (foxl2, β-catenin, and sry) were expressed highly in intersex individuals and were likely involved in the control of gonadal sex in C. nobilis. High expression of foxl2 or β-catenin may inhibit sry and activate 5-HT receptor and vitellogenin to maintain female development. High expression of sry may inhibit foxl2 and β-catenin and activate dmrt2, fem-1, sfp2, sa6, Amy-1, APCP4, and PLK to maintain male function. High expression of sry, foxl2, and β-catenin in C. nobilis may be involved in promoting and maintaining sex reversal/differentiation. The downstream regulator may not be dimorphic expressed genes, but genes expressed in intersex individuals, males and females. Different expression patterns of sex-related genes and gonadal histological characteristics suggested that C. nobilis may change its sex from male to female. These findings suggest highly conserved sex reversal/differentiation with diverged regulatory pathways during C. nobilis evolution. This study provides valuable genetic resources for understanding sex reversal/differentiation (intersex) mechanisms and pathways underlying bivalve reproductive regulation.

Presence of gynogenetic males suggests a female heterogamety in sterlet Acipenser ruthenus L

Investigation of the heterogametic sex in sterlet Acipenser ruthenus L. was performed using meiotic gynogenesis and gonadal histology. Eggs from the albino females were activated by UV irradiated sperm of wild colored males and exposed to a heat shock. The resultant fish were all albino and exhibited exclusively maternal inheritance of the microsatellite DNA markers. Cytogenetic analysis indicated that gynogenetic progeny were diploids with 120 chromosomes. Based on the histological analysis, more than 86% of the gynogenetic individuals were found to be females. Moreover, some males (7%), sterile speciemens (3.5%) and fish with unidentified gonads (3.5%) were observed among the gynogenetic fish. Presence of both females and males in the gynogenetic offspring is indicative that the heterogametic sex in sterlets is female.

Characterization of a low-density lipoprotein receptor, Lrp13, in Chinese tongue sole (Cynoglossus semilaevis) and medaka (Oryzias latipes)

As an important economic marine species cultured in China, Chinese tongue sole (Cynoglossus semilaevis) has interested us due to its sexual dimorphism and ZW/ZZ sex determination system. In a previous study, dmrt1 was identified as a dosage-dependent male-determining gene. In the present study, a female-specific expressed gene, cse0440, initially annotated as lrp1b-like, was identified from chromosome W of C. semilaevis. In view of the differences between cse0440 and lrp1b in terms of expression pattern, a phylogenetic analysis containing 85 LRP proteins was constructed and provided an evidence to re-annotate cse0440 as cseLRP13. In addition, two orthologues of cseLRP13 were separately identified from W and Z chromosomes: cseLRP13-W and cseLRP13-Z. The subsequent multiple sequence alignment and syntenic arrangements of LRP13 in C. semilaevis, Japanese medaka (Oryzias latipes), large yellow croaker (Larimichthys crocea), striped bass (Morone saxatilis), white perch (Morone americana) and Fugu rubripes (Takifugu rubripes) further supported this re-annotation. RT-PCR and in situ hybridization revealed that cselrp13 was exclusively expressed in the oocytes and follicles of ovaries. These results suggested that lrp13 may play important roles in female reproduction. In future, with the advancement of micromanipulation in flatfish, the detailed function of two lrp13 orthologues in C. semilaevis will be elucidated.

Ribosome RNA Profiling to Quantify Ovarian Development and Identify Sex in Fish

Terminologies of ovary development, by somewhat subjective describing and naming main changes of oocytes, have been criticized for confusing and inconsistency of terms and classifications, and the incurred consequences impede communication among researchers. In the present work, we developed regression between ovary development and three ribosome RNA (rRNA) indexes, namely 5S rRNA percent, 18S rRNA percent, and 5S–18S rRNA ratio, using close relationship between volume percent of primary growth stage oocytes or gonadosomatic index and rRNA content, demonstrating species-specific quantification of ovary development can be established in species with either synchronous and asynchronous oogenesis. This approach may be extended to any species with primary growth oocytes, e.g. anurans and reptiles, to predict maturity stages in females. We further confirmed that 5S rRNA percent and 5S/18S rRNA ratio can serve as markers to distinguish sexes unambiguously. A micro-invasive sampling method may be invented for non-lethal prediction of ovary development and sex because only a small amount of ovary sample (<50 mg) is needed for the approach established in the current work. Researchers who work with ovary RNA-seq in these taxa should realize that insufficient depletion of rRNA will probably lead to incorrect quantification of gene expression and inaccurate conclusions.

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

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