Molecular players involved in temperature-dependent sex determination and sex differentiation in Teleost fish

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.

Identification and characterization of the Dmrt1B gene in the oriental river prawn, Macrobrachium nipponense

Dmrt (doublesex and mab-3 related transcription factor) is a protein family of transcription factors implicated in sexual regulation. Dmrt proteins are widely conserved and known for their involvement in sex determination and differentiation across species, from invertebrates to humans. In this study, we identified a novel gene with a DM (doublesex/Mab-3)-domain gene in the river prawn, Macrobrachium nipponense, which we named MniDmrt1B due to its similarities and close phylogenetic relationship with Dmrt1B in Macrobrachium rosenbergii. Through amino acid alignments and structural predictions, we observed conservation and identified putative active sites within the DM domain. qRT-PCR analysis revealed that MniDmrt1B exhibited high expression levels in the testis, with consistently higher expression in males compared to females during development. Additionally, similar to other sex-regulated genes, the MniDmrt1B gene exhibited high expression levels during the sex differentiation-sensitive periods in M. nipponense. These results strongly indicated that MniDmrt1B probably plays an important role in testis development and sex differentiation in M. nipponense.

Molecular and Physiological Effects of 17α-methyltestosterone on Sex Differentiation of Black Rockfish, Sebastes schlegelii

It is widely known that all-female fish production holds economic value for aquaculture. Sebastes schlegelii, a preeminent economic species, exhibits a sex dimorphism, with females surpassing males in growth. In this regard, achieving all-female black rockfish production could significantly enhance breeding profitability. In this study, we utilized the widely used male sex-regulating hormone, 17α-methyltestosterone (MT) at three different concentrations (20, 40, and 60 ppm), to produce pseudomales of S. schlegelii for subsequent all-female offspring breeding. Long-term MT administration severely inhibits the growth of S. schlegelii, while short term had no significant impact. Histological analysis confirmed sex reversal at all MT concentrations; however, both medium and higher MT concentrations impaired testis development. MT also influenced sex steroid hormone levels in pseudomales, suppressing E2 while increasing T and 11-KT levels. In addition, a transcriptome analysis revealed that MT down-regulated ovarian-related genes (cyp19a1a and foxl2) while up-regulating male-related genes (amh) in pseudomales. Furthermore, MT modulated the TGF-β signaling and steroid hormone biosynthesis pathways, indicating its crucial role in S. schlegelii sex differentiation. Therefore, the current study provides a method for achieving sexual reversal using MT in S. schlegelii and offers an initial insight into the underlying mechanism of sexual reversal in this species.

Innovative approach for high-throughput exploiting sex-specific markers in Japanese parrotfish Oplegnathus fasciatus

BACKGROUND

The use of sex-specific molecular markers has become a prominent method in enhancing fish production and economic value, as well as providing a foundation for understanding the complex molecular mechanisms involved in fish sex determination. Over the past decades, research on male and female sex identification has predominantly employed molecular biology methodologies such as restriction fragment length polymorphism, random amplification of polymorphic DNA, simple sequence repeat, and amplified fragment length polymorphism. The emergence of high-throughput sequencing technologies, particularly Illumina, has led to the utilization of single nucleotide polymorphism and insertion/deletion variants as significant molecular markers for investigating sex identification in fish. The advancement of sex-controlled breeding encounters numerous challenges, including the inefficiency of current methods, intricate experimental protocols, high costs of development, elevated rates of false positives, marker instability, and cumbersome field-testing procedures. Nevertheless, the emergence and swift progress of PacBio high-throughput sequencing technology, characterized by its long-read output capabilities, offers novel opportunities to overcome these obstacles.

FINDINGS

Utilizing male/female assembled genome information in conjunction with short-read sequencing data survey and long-read PacBio sequencing data, a catalog of large-segment (>100 bp) insertion/deletion genetic variants was generated through a genome-wide variant site-scanning approach with bidirectional comparisons. The sequence tagging sites were ranked based on the long-read depth of the insertion/deletion site, with markers exhibiting lower long-read depth being considered more effective for large-segment deletion variants. Subsequently, a catalog of bulk primers and simulated PCR for the male/female variant loci was developed, incorporating primer design for the target region and electronic PCR (e-PCR) technology. The Japanese parrotfish (Oplegnathus fasciatus), belonging to the Oplegnathidae family within the Centrarchiformes order, holds significant economic value as a rocky reef fish indigenous to East Asia. The criteria for rapid identification of male and female differences in Japanese parrotfish were established through agarose gel electrophoresis, which revealed 2 amplified bands for males and 1 amplified band for females. A high-throughput identification catalog of sex-specific markers was then constructed using this method, resulting in the identification of 3,639 (2,786 INS/853 DEL, ♀ as reference) and 3,672 (2,876 INS/833 DEL, ♂ as reference) markers in conjunction with 1,021 and 894 high-quality genetic sex identification markers, respectively. Sixteen differential loci were randomly chosen from the catalog for validation, with 11 of them meeting the criteria for male/female distinctions. The implementation of cost-effective and efficient technological processes would facilitate the rapid advancement of genetic breeding through expediting the high-throughput development of sex genetic markers for various species.

CONCLUSIONS

Our study utilized assembled genome information from male and female individuals obtained from PacBio, in addition to data from short-read sequencing data survey and long-read PacBio sequencing data. We extensively employed genome-wide variant site scanning and identification, high-throughput primer design of target regions, and e-PCR batch amplification, along with statistical analysis and ranking of the long-read depth of the variant sites. Through this integrated approach, we successfully compiled a catalog of large insertion/deletion sites (>100 bp) in both male and female Japanese parrotfish.

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.

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.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

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.

High-temperature stress will put the thermo-sensitive teleost yellow catfish (Tachysurus fulvidraco) in danger through reducing reproductivity

Recently, concerns for species that sex differentiation is influenced by temperature in the context of global warming have increased because disrupted operational sex ratios could threaten population maintenance. In contrast, little attention has been given to the reproductive ability of populations that experienced elevated temperatures. In this study, we demonstrated that high temperature (HT) would decrease population size via three different aspects of reproductive ability for the first time. We show that, in a thermo-sensitive teleost yellow catfish, a short period of HT (+3 °C) exposure during the critical period of sex differentiation leads to a different percentage of masculinization of XX genotypic females (1-23%) in wet-lab and natural water bodies. Combining the results of gonadal appearance, histology, sperm parameters, and fertilization rate, we found that XX pseudo-males induced by HT display significantly discounted fertility and reproductive performance compared to XY normal males. We demonstrate that the survival of the XY genotype is lower than XX genotype under environmental stress, including HT, hypoxia, and parasite infection, and the differential survival seems unrelated to male-biased sexual size dimorphism. The mathematical model predicts that the phenotypic female percent will be stabilized at 50% and the population will be sustainably maintained when masculinizing force is less than 0.5, while HT will put the population in danger when the masculinizing force exceeds 0.5. However, when we combine the real-world data of reproductive ability and mathematic model, our results suggest the population size decreases and the long-term survival of the studied species are threatened under the projected pace of increasing temperature. These findings will be useful for understanding the long-term effects of increasing temperature on sex ratio, reproduction and population maintenance in teleost.

Timing of early gonadal differentiation and effects of estradiol-17β treatments on the sex differentiation in Largemouth bass (Micropterus salmoides)

In this study, an efficient estradiol-17β (E₂)-induced feminization method was established based on the timing of early gonadal differentiation in Largemouth bass (Micropterus salmoides). Histological section results showed that from 20 days post-hatch (dph) to 30 dph, the germ cells gradually differentiated into oogonium and spermatic deferent, respectively. Moreover, female-biased genes Foxl2 and Cyp19a1a were up-regulated to the first peak at 20 dph, while the male-biased genes Dmrt1 were up-regulated to the first peak at 30 dph. These results indicated that the timing of early gonadal differentiation in Largemouth bass was between 20 and 30 dph. Therefore, 15 dph Largemouth bass with a body length of 15.10 ± 0.09 mm were chosen, and four E₂-treated diets were set as 0 (E0, control), 50 mg/kg E₂ (E50), 100 mg/kg E₂ (E100), and 200 mg/kg E₂ (E200). After feeding with E₂-treated diets for 60 days, female ratios were 55%, 100%, 100%, and 100% in E0, E50, E100, and E200 groups, respectively. No intersex fish were observed in all the groups. However, 30% of females in the E200 group possessed thinner ovaries, with smaller ovary cavity structures and a decreased number of primary oocyte cells than those in other groups. Besides, the Largemouth bass in the E0 group grew more than those in E50, E100, and E200 groups during the E₂ treatments period (P < 0.05). In conclusion, our study suggested that 50-100 mg/kg E₂-treated diets could effectively induce the feminization of 15 dph Largemouth bass within 60 days duration time, which provided valuable information for the breeding of the all-male Largemouth bass population.

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.

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

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

A Multidimensional Matrix Model for Predicting the Effects of Male-Biased Sex Ratios on Fish Populations

Laboratory experiments have established that exposure to certain endocrine-active substances prior to and/or during the period of sexual differentiation can lead to skewed sex ratios in fish. However, the potential long-term population impact of biased sex ratio depends on multiple factors including the life history of the species and whether the ratio is male or female-biased. In the present study, we describe a novel multidimensional, density-dependent matrix model that analyzes age class-structure of both males and females over time, allowing for the quantitative evaluation of the effects of biased sex ratio on population status. This approach can be used in conjunction with field monitoring efforts and/or laboratory testing to link effects on sex ratio due to chemical and/or nonchemical stressors to adverse outcomes in whole organisms and populations. For demonstration purposes, we applied the model to evaluate population trajectories for fathead minnow (Pimephales promelas) exposed to prochloraz, an aromatase inhibitor, during sexual differentiation. The model also was used to explore the population impact in a more realistic exposure scenario in which both adult and early life stages of fish are exposed concurrently to prochloraz, which, in addition to altering sex ratio during development, can decrease vitellogenin and egg production in adult females. For each exposure scenario, the model was used to analyze total population size, numbers of females and of males, and sex specific recruitment of the F1 generation. The present study illustrates the utility of multidimensional matrix population models for ecological risk assessment in terms of integrating effects across a population of an organism even when chemical effects on individuals are manifested via different pathways depending on life stage. Environ Toxicol Chem 2022;41:1066-1077. Published 2022. This article is a U.S. Government work and is in the public domain in the USA.

Intensive masculinization caused by chronic heat stress in juvenile Cynoglossus semilaevis: Growth performance, gonadal histology and gene responses

The sea temperature has been observed to chronically increase during the past decades, leaving unpredictable influences to the marine biological resources. Thus, it is of vital significance to study the biological responses of ocean inhabited organisms with the artificially stimulated heat stress environment. Cynoglossus semilaevis provides us with an ideal model to study the influence of chronic heat stress on the sexual differentiation in marine teleosts for its genetic sex determination (GSD) + environmental effected (EE) sex determination system. In this study, the comparative experiment was conducted employing heated seawater (HT group) and ambient seawater (CT group) to cultivate juvenile C. semilaevis respectively. Significant differences were exhibited in growth performance and a delayed germ cell development effect was found in pseudomales formed under chronic heat stress. Using transcriptome analysis, the transcription profile of 55 days post fertilization (dpf) and 100 dpf juveniles' gonads were studied. A total of 47 libraries were constructed with an average mapping rate of 94.63% after assembling. GO and KEGG enrichment were proceeded using DEGs screened out between (1) pseudomale gonads at 55 dpf and 100 dpf in HT and CT group (2) pseudomale and female gonads at 55 dpf and 100 dpf in HT and CT group. Terms and pathways involved in steroid stimulation, reproduction ability, germ cell proliferation et al. were shed light on. The expression pattern of 29 DEGs including amh, hsp90b1, pgr et al. were also provided to supplement the results of functional enrichment. Weighted gene co-expression networks analysis (WGCNA) was constructed and hspb8-like, histone H2A.V were exhibited to play vital roles in the heat-induced masculinization. Our findings facilitate the understanding for transcriptional variations in intensive masculinization cause by chronic heat stress of C. semilaevis and provide referable study of the influences on the teleosts in elevated sea temperature.

Effects of global warming on species with temperature-dependent sex determination: Bridging the gap between empirical research and management

Global warming could threaten over 400 species with temperature-dependent sex determination (TSD) worldwide, including all species of sea turtle. During embryonic development, rising temperatures might lead to the overproduction of one sex and, in turn, could bias populations' sex ratios to an extent that threatens their persistence. If climate change predictions are correct, and biased sex ratios reduce population viability, species with TSD may go rapidly extinct unless adaptive mechanisms, whether behavioural, physiological or molecular, exist to buffer these temperature-driven effects. Here, we summarize the discovery of the TSD phenomenon and its still elusive evolutionary significance. We then review the molecular pathways underpinning TSD in model species, along with the hormonal mechanisms that interact with temperatures to determine an individual's sex. To illustrate evolutionary mechanisms that can affect sex determination, we focus on sea turtle biology, discussing both the adaptive potential of this threatened TSD taxon, and the risks associated with conservation mismanagement.

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

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

Epigenetic Regulation and Environmental Sex Determination in Cichlid Fishes

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

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

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

R-spondin1 in loach (Misgurnus anguillicaudatus): Identification, characterization, and analysis of its expression patterns and DNA methylation in response to high-temperature stress

With a well-understood function in mammals, R-spondin1 (Rspo1) is an important regulator of ovarian development via the Wnt/β-catenin pathway. Rspo1 deficiency causes retardation of ovarian development in XX fish, and increases Rspo1 function induces femininity and sex reversal in XY fish. In this study, Rspo1 was successfully cloned from loach (Misgurnus anguillicaudatus), and its expression profile was analyzed. The full-length cDNA of Misgurnus anguillicaudatus Rspo1 (MaRspo1) comprised 1322 bp and included an open reading frame (ORF) of 795 bp, which encoded a predicted polypeptide measuring 264 amino acids in length. Phylogenetic and gene structure analyses showed a highly conserved sequence of MaRspo1 (identical to the Rspo1 genes of other species), consisting of an N-terminal signal peptide (SP), two furin-like cysteine-rich domains (FU1 and FU2), a thrombospondin type 1 repeat (TSP1) and a C-terminal region. Real-time PCR revealed the female-biased expression profile of MaRspo1, with the highest expression level among tested tissues detected in ovary. Investigation of MaRspo1 expression levels throughout the early development stage (10-60 days post hatching) under three temperature treatments (25 °C, 28 °C, and 31 °C) revealed significantly differential expression of MaRspo1 among the three temperature groups, with decreased MaRspo1 expression in the high-temperature (31 °C) group. The results of DNA methylation analysis indicated that exposure to high temperature during early development can increase the average promoter methylation level of MaRspo1 in both females and males. Taken together, the results of this study provide the basis for the further investigation of the molecular mechanism of Rspo1 in response to temperature.

Epigenetic regulation of gonadal and brain aromatase expression in a cichlid fish with environmental sex determination

A fit animal must develop testes or ovaries, with brain and physiology to match. In species with alternative male morphs this coordination of development across tissues operates within sexes as well as between. For Pelvicachromis pulcher, an African cichlid in which early pH exposure influences both sex and alternative male morph, we sequence both copies of aromatase (cyp19a1), a key gene for sex determination. We analyze gene expression and epigenetic state, comparing gonad and brain tissue from females, alternative male morphs, and fry. Relative to brain, we find elevated expression of the A-copy in the ovaries but not testes. Methylation analysis suggests strong epigenetic regulation, with one region specifying sex and another specifying tissue. We find elevated brain expression of the B-copy with no sex or male morph differences. B-copy methylation follows that of the A-copy rather than corresponding to B-copy expression. In 30-day old fry, we see elevated B-copy expression in the head, but we do not see the expected elevated A-copy expression in the trunk that would reflect ovarian development. Interestingly, the A-copy epialleles that distinguish ovaries from testes are among the most explanatory patterns for variation among fry, suggesting epigenetic marking of sex prior to differentiation and thus laying the groundwork for mechanistic studies of epigenetic regulation of sex and morph differentiation.

Molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with XX/XY sex determination

Phenotypic sex of an organism is determined by molecular changes in the gonads, so-called molecular sex differentiation, which should precede the rise of cellular or anatomical sex-distinguishing features. This study characterized molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with established XX/XY genotypic sex determination. Next generation sequencing was conducted on sablefish ovarian and testicular mRNAs to obtain sequences for transcripts associated with vertebrate sex determination and differentiation and early reproductive development. Gene-specific PCRs were developed to determine the distribution and ontogenetic gonadal expression of transcription, growth, steroidogenic and germline factors, as well as gonadotropin and steroid receptors. Molecular changes associated with sex differentiation were first apparent in both XY- and XX-genotype sablefish at ~ 60 mm in body length and prior to histological signs of sex differentiation. The earliest and most robust markers of testicular differentiation were gsdf, amh, dmrt1, cyp11b, star, sox9a, and fshr. Markedly elevated mRNA levels of several steroidogenesis-related genes and ar2 in differentiating testes suggested that androgens play a role in sablefish testicular differentiation. The earliest markers of ovarian differentiation were cyp19a1a, lhcgr, foxl2, nr0b1, and igf3. Other transcripts such as figla, zp3, and pou5f3 were expressed predominantly in XX-genotype fish and significantly increased with the first appearance and subsequent development of primary oocytes. This study provides valuable insight to the developmental sequence of events associated with gonadal sex differentiation in marine teleosts with XX/XY sex determination. It also implicates particular genes in processes of male and female development and establishes robust molecular markers for phenotypic sex in sablefish, useful for ongoing work related to sex control and reproductive sterilization.

Gonad-Specific Transcriptomes Reveal Differential Expression of Gene and miRNA Between Male and Female of the Discus Fish (Symphysodon aequifasciatus)

The discus fish (Symphysodon aequifasciatus) is an ornamental fish that is well-known around the world. In artificial reproduction, they must be matched by one male and one female, whereas phenotype investigation indicated that there are no significant differences in appearance between males and females, which causes great difficulties in the mating during artificial reproduction. So, it is of great importance to establish artificial sex identification methods for the discus fish. The molecular mechanism of the sexual dimorphism of the discus fish was previously unknown. In this study, we constructed six cDNA libraries from three adult testes and three adult ovaries and performed RNA sequencing for identifying sex-biased candidate genes and microRNAs (miRNAs). A total of 50,082 non-redundant genes (unigenes) were identified, of which 18,570 unigenes were significantly overexpressed in testes, and 11,182 unigenes were significantly overexpressed in ovaries. A total of 551 miRNAs were identified, of which 47 miRNAs were differentially expressed between testes and ovaries. Eight differentially expressed unigenes, seven differentially expressed miRNAs and one non-differential miRNA were validated by quantitative real-time polymerase chain reaction. Twenty-four of these differentially expressed miRNAs and their 15 predicted target genes constituted 41 miRNA–mRNA interaction pairs, and some of vital sex-related metabolic pathways were also identified. These results revealed these differentially expressed genes and miRNAs between ovary and testis might be involved in regulating gonadal development, sex determination, gametogenesis, and physiological function maintenance, and there are complex regulatory networks between genes and miRNAs. It can help us understand the molecular mechanism of the sexual dimorphism and obtain a high-efficiency sex identification method in the artificial reproduction process of the discus fish.

Roles of forkhead box protein L2 (foxl2) during gonad differentiation and maintenance in a fish, the olive flounder (Paralichthys olivaceus)

As an important maricultured fish, the olive flounder Paralichthys olivaceus shows sex-dimorphic growth. Thus, the molecular mechanisms involved in sex control in P. olivaceus have attracted researchers' attention. Among the sex-related genes, forkhead box protein L2 (foxl2) exhibits significant sex-dimorphic expression patterns and plays an important role in fish gonad differentiation and development. The present study first investigated the expression levels and promoter methylation dynamics of foxl2 during flounder gonad differentiation under treatments of high temperature and exogenous 17β-oestradiol (E2). During high temperature treatment, the expression of flounder foxl2 may be repressed via maintenance of DNA methylation. Then, flounder with differentiated testis at Stages I-II were treated with exogenous 5ppm E2 or 5ppm E2+150ppm trilostane (TR) to investigate whether exogenous sex hormones could induce flounder sex reversal. The differentiated testis exhibited phenotypic variations of gonadal dysgenesis with upregulation of female-related genes (foxl2 and cytochrome P450 family 19 subfamily A (cyp19a)) and downregulation of male-related genes (cytochrome P450 family 11 subfamily B member 2 (cyp11b2), doublesex- and mab-3 related transcription factor 1 (dmrt1), anti-Mullerian hormone (amh) and SRY-box transcription factor 9 (sox9)). Furthermore, a cotransfection assay of the cells of the flounder Sertoli cell line indicated that Foxl2 was able alone or with nuclear receptor subfamily 5 group A member 2 (Nr5a2) jointly to upregulate expression of cyp19a. Moreover, Foxl2 and Nr5a2 repressed the expression of dmrt1. In summary, Foxl2 may play an important role in ovarian differentiation by maintaining cyp19a expression and antagonising the expression of dmrt1. However, upregulation of foxl2 is not sufficient to induce the sex reversal of differentiated testis.

Identification of male heterogametic sex-determining regions on the Atlantic herring Clupea harengus genome

The sex determination system of Atlantic herring Clupea harengus L., a commercially important fish, was investigated. Low coverage whole-genome sequencing of 48 females and 55 males and a genome-wide association study revealed two regions on chromosomes 8 and 21 associated with sex. The genotyping data of the single nucleotide polymorphisms associated with sex showed that 99.4% of the available female genotypes were homozygous, whereas 68.6% of the available male genotypes were heterozygous. This is close to the theoretical expectation of homo/heterozygous distribution at low sequencing coverage when the males are factually heterozygous. This suggested a male heterogametic sex determination system in C. harengus, consistent with other species within the Clupeiformes group. There were 76 protein coding genes on the sex regions but none of these genes were previously reported master sex regulation genes, or obviously related to sex determination. However, many of these genes are expressed in testis or ovary in other species, but the exact genes controlling sex determination in C. harengus could not be identified.

Amh dominant expression in Sertoli cells during the testicular differentiation and development stages in the olive flounder Paralichthys olivaceus

The olive flounder Paralichthys olivaceus, an important marine fish, shows gender differences in growth. The mechanism on its gonadal differentiation direction affected with exogenous factors still needs to be clarified. The anti-Müllerian hormone (amh) gene is involved in fish testicular differentiation and maintenance. The aim of this study was to investigate the expression of the flounder amh in tissues and the gonads. The quantitative expression analysis results showed that it was highly expressed in the testis, especially in the testis at stages I - IV (P < 0.05). Also, amh was detected in Sertoli cells surrounding spermatogonia and peripheral seminiferous lobule of the testis with in situ hybridization (ISH) and immunohistochemistry (IHC). During the differentiation period, the amh expression in the testis of the tamoxifen treatment group (100 ppm) was higher than that in the ovary of the 17β-estradiol (E2, 5 ppm) group, and the expression levels of amh during process of the male differentiation in the tamoxifen group were higher than those in the 17ɑ-methyltestosterone (MT, 5 ppm) group (P < 0.05). ISH results also exhibited that amh was expressed in the somatic cells that surrounded the germ cells of juvenile flounder similar to adult ones. Furthermore, the flounder gonads in the tamoxifen group maintained more germ cells and somatic cells than those in the MT group from 20 to 80 mm total length (TL). Especially, at 60 and 80 mm TL, the numbers of germ and somatic cells in the tamoxifen group were significantly higher than those in the MT group (P < 0.05). In summary, amh might initiate the process of testicular differentiation, and is involved in the early development and maintenance of testis.

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

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

High-Density Linkage Map and Mapping for Sex and Growth-Related Traits of Largemouth Bass (Micropterus salmoides)

The largemouth bass is an important species, and its culture has risen sharply with the surge in fish aquaculture in China. Due to the lack of selective breeding technology for the largemouth bass, the growth rate and disease resistance are low, its sexual maturation is slow, and other serious problems are contributing to a sharp decline in the safety and quality of largemouth bass products in recent decades. Therefore, comprehensive breeding programs to improve the economic performance and promote the modern industrial development of largemouth bass must be considered a priority. Here, a total of 152 adult largemouth bass, including two parents and 150 progenies, were selected to produce the genetic mapping family. Then, a high-density linkage map was constructed based on restriction site–associated DNA sequencing using 6,917 single-nucleotide polymorphisms (SNPs) located in 24 linkage groups (LGs). The total genetic length of the linkage map was 1,261.96 cM, and the length of each LG varied from 24.72 cM for LG02 to 117.53 cM for LG16, with an average length of 52.58 cM and an average SNP number of 286. Thirteen significant quantitative trait loci (QTLs) for sex determination were located on LG04, LG05, LG08, LG12, LG15, LG21, and LG23. An informative QTL cluster that included six QTLs was detected on LG12. However, one notable QTL, which accounted for 71.48% of the total phenotypic variation, was located in the region of 1.85 cM on LG05. In addition, 32 identified QTLs were related to growth, including body weight, body length, body height, and head length. The QTLs for these growth-related traits are located in 13 LG regions and have little effect on phenotypic variation. This high-density genetic linkage map will enable the fine-mapping of economic traits and support the future genome assembly of the largemouth bass. Additionally, our study will be useful for future selective culture of largemouth bass and could potentially be used in molecular-assisted breeding of largemouth bass for aquaculture.

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.

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

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

Warmer waters masculinize wild populations of a fish with temperature-dependent sex determination

Southern flounder (Paralichthys lethostigma) exhibit environmental sex determination (ESD), where environmental factors can influence phenotypic sex during early juvenile development but only in the presumed XX female genotype. Warm and cold temperatures masculinize fish with mid-range conditions producing at most 50% females. Due to sexually dimorphic growth, southern flounder fisheries are dependent upon larger females. Wild populations could be at risk of masculinization from ESD due to globally increasing water temperatures. We evaluated the effects of habitat and temperature on wild populations of juvenile southern flounder in North Carolina, USA. While northern habitats averaged temperatures near 23 °C and produced the greatest proportion of females, more southerly habitats exhibited warmer temperatures (>27 °C) and consistently produced male-biased sex ratios (up to 94% male). Rearing flounder in the laboratory under temperature regimes mimicking those of natural habitats recapitulated sex ratio differences observed across the wild populations, providing strong evidence that temperature is a key factor influencing sex ratios in nursery habitats. These studies provide evidence of habitat conditions interacting with ESD to affect a key demographic parameter in an economically important fishery. The temperature ranges that yield male-biased sex ratios are within the scope of predicted increases in ocean temperature under climate change.

High temperature increases the gsdf expression in masculinization of genetically female Japanese flounder (Paralichthys olivaceus)

In teleosts, sex is plastic and is influenced by environmental factors. Elevated temperatures have masculinizing effects on the phenotypic sex of certain sensitive species. In this study, we reared genetic XX Japanese flounder at a high temperature (27.5 ± 0.5 °C) and obtained a population of sex-reversal XX males (male ratio, 95.24%). We comparatively analyzed the dynamic characteristics of germ cells and gsdf (gonadal soma-derived factor) expression during sexual differentiation for the experimental (27.5 ± 0.5 °C) and control (18 °C ± 0.5 °C) groups. The results revealed that the germ cell proliferation inhibited and gsdf expression up-regulated in the experimental group, and the gsdf mRNA and proteins expressed in somatic cells that had direct contact with germline stem cells (with Nanos 2 protein expression) including spermatogonia and oogonia by ISH (in situ hybridization) and IHC (immunohistochemistry). In addition, we also overexpressed the gsdf in XX flounders, and the germ cell number of XX flounders bearing gsdf gene significantly decreased and sometimes disappeared completely, which was consistent with the results from high-temperature induction. Therefore, based on all the results, we speculated that the high expression of gsdf might inhibit germ cell proliferation during sex differentiation, and eventually cause sex reversal in the high-temperature induced masculinization of XX Japanese flounder.

The effects of bifenthrin and temperature on the endocrinology of juvenile Chinook salmon

The San Francisco Bay delta (USA) is experiencing seasonally warmer waters attributable to climate change and receives rainstorm runoff containing pyrethroid pesticides. Chinook salmon (Oncorhynchus tshawytscha) inhabit the affected waterways from hatch through smoltification, and thus juvenile fish may experience both pyrethroid and warmer water exposures. The effects of higher temperatures and pesticide exposure on presmolt Chinook are unknown. To improve understanding of the potential interaction between temperature and pesticide exposure on salmonid development, juvenile alevin and fry were reared in 11, 16.4, and 19 °C freshwater for 11 d and 2 wk, respectively, and exposed to nominal concentrations of 0, 0.15, and 1.5 µg/L bifenthrin for the final 96 h of rearing. Estradiol-17β (E2), testosterone, triiodothyronine, and thyroxine levels were measured in whole-body homogenates using hormone-specific enzyme-linked immunosorbent assays. Brain gonadotropin-releasing hormone receptor (GnRH2), dopamine receptor 2A, and growth hormone 1 (GH1) mRNA levels were measured using quantitative PCR. Results showed significantly decreased survival and condition factors observed with increasing temperature in alevin. Alevin thyroid hormones increased significantly with temperature, but fry thyroid hormones trended toward a decrease at lower temperatures with increasing bifenthrin exposure. There were significant reductions in fry testosterone and E2 at 11 °C with increasing bifenthrin treatments and significant changes in GnRH2 and GH1 gene expression in both alevin and fry, indicating potential disruption of hormonal and signaling pathways. Environ Toxicol Chem 2019;38:852-861. © 2019 SETAC.

Caste-specific development of the dopaminergic system during metamorphosis in female honey bees

Caste-specific differences in the dopaminergic systems of social insects assist in maintaining caste-specific behavior. To determine how caste differences in the honey bee occur during metamorphosis, a number of comparative analyses between castes were performed including comprehensive quantification of: levels of dopamine and its metabolite in the brain, the gene expression levels of enzymes involved in dopamine biosynthesis and conversion as well as expression levels of dopamine receptors and a dopamine transporter. Dopamine levels standardized to the protein contents of a whole brain at the day of eclosion in queens were 3.6-fold higher than those in workers. Dopamine levels increased until eclosion (7 days) in queens, whereas those in workers increased until 5–6 days before eclosion and then maintained until eclosion (10 days). These caste-specific dopamine dynamics in the brain were supported by the higher expression of genes (Amddc and Amth) encoding enzymes involved in dopamine synthesis in queens. The distribution of cells expressing Amddc in the brain revealed that soma clusters of dopaminergic cells were similar between castes at 7–8 days after pupation, suggesting the upregulation of Amddc expression in some cells in queens rather than addition of cell clusters. In contrast, genes encoding dopamine receptors were downregulated in queens or showed similar expression levels between castes. The expression of genes encoding dopamine transporters did not differ between castes. These results reveal the developmental process of caste-specific dopaminergic systems during metamorphosis in the honey bee, suggesting caste-specific behavior and division of reproduction in this highly eusocial species.

Effects of Temperature on the Expression of Two Ovarian Differentiation-Related Genes foxl2 and cyp19a1a

Exposure to stress induces a series of responses and influences a wide range of biological processes including sex differentiation in fish. The present work investigated the molecular and physiological response to thermal stress throughout the early development stage covering the whole period of sex differentiation of bluegill, Lepomis macrochirus. Larvae were treated using three temperatures, 17, 24, and 32°C from 6 to 90 days posthatching (dph) in 30-L round tanks. There is no significant difference of the sex ratio and survival among the three temperature groups in the geographic population used in this study. Two ovarian differentiation-related genes foxl2 and cyp19a1a were detected at 7 dph suggesting that these genes have already played a role prior to sex differentiation. The expression of foxl2 reached the peak and was thermosensitive just prior to the onset of ovarian differentiation at 27 dph. Histological examination displayed that the proliferation of germ cells and ovarian differentiation were delayed at the low-temperature treatment (17°C) at 97 dph compared with higher temperatures. In conclusion, the water temperature regulates the sex differentiation of bluegill through modulation of the expression of foxl2 and cyp19a1a. A comparative study of the expression profile of sex differentiation-related genes in species will shed light on the evolution of sex-determination mechanisms and the impact of stress on sex differentiation.

Distinct sperm nucleus behaviors between genotypic and temperature-dependent sex determination males are associated with replication and expression-related pathways in a gynogenetic fish

BACKGROUND

Coexistence and transition of diverse sex determination strategies have been revealed in some ectothermic species, but the variation between males caused by different sex determination strategies and the underlying mechanism remain unclear. Here, we used the gynogenetic gibel carp (Carassius gibelio) with both genotypic sex determination (GSD) and temperature-dependent sex determination (TSD) strategies to illustrate this issue.

RESULTS

We found out that males of GSD and TSD in gibel carp had similar morphology, testicular histology, sperm structure and sperm vitality. However, when maternal individuals were mated with males of GSD, sperm nucleus swelling and fusing with the female pronucleus were observed in the fertilized eggs. On the contrary, when maternal individuals were mated with males of TSD, sperm nucleus remained in the condensed status throughout the whole process. Subsequently, semen proteomics analysis unveiled that DNA replication and gene expression-related pathways were inhibited in the sperm from males of TSD compared to males of GSD, and most differentially expressed proteins associated with DNA replication, transcription and translation were down-regulated. Moreover, via BrdU incorporation and immunofluorescence detection, male nucleus replication was revealed to be present in the fertilized eggs by the sperm from males of GSD, but absent in the fertilized eggs by the sperm from males of TSD.

CONCLUSIONS

These findings indicate that DNA replication and gene expression-related pathways are associated with the distinct sperm nucleus development behaviors in fertilized eggs in response to the sperm from males of GSD and TSD. And this study is the first attempt to screen the differences between males determined via GSD and TSD in gynogenetic species, which might give a hint for understanding evolutionary adaption of diverse sex determination mechanisms in unisexual vertebrates.

Expression profiles of amh and foxl2 in Schizothorax kozlovi, and their response to temperature during the early developmental stage

To elucidate the role of amh and foxl2 in sex differentiation of the teleost fish Schizothorax kozlovi, the full-length cDNAs were cloned from the mature testis and ovary by rapid amplification of cDNA ends (RACE), and their relative mRNA expression levels were determined by quantitative real-time polymerase chain reaction among tissues and temperature groups. The complete amh and foxl2 cDNAs of S. kozlovi were 2060 bp and 1750 bp, which encoded 568 and 306 amino acids, respectively. The amh were expressed only in gonads, while foxl2 was expressed in the gills, brain and gonads, both exhibiting relatively high tissue specificity. The amh exhibited sex-specific expression pattern in the gonads. No sex differences in the foxl2 expression were observed in the brain and gonads, but significant sex differences were found in the gills. No significant differences were found in the foxl2 expression, from the larval to the juvenile stage, and also between different temperature groups. However, significant differences were found in the expression levels of amh from the larval (12-63 days posthatching (dph)) to the juvenile stage (190 dph), and also among the 18°C and 10°C groups at 31 dph. This result suggested that amh plays an important role in male sex differentiation of S. kozlovi during the early developmental stage, but no similar effect was observed in foxl2.

Transcriptomic analysis of the differentiating ovary of the protogynous ricefield eel Monopterus albus

BACKGROUND

The ricefield eel is a protogynous hermaphroditic Synbranchiform species that changes sex naturally from female to male, which offers an interesting model for studying gonadal (particularly ovarian) differentiation in vertebrates. In the present study, transcriptome sequencing of the gonad of ricefield eel larvae was performed to explore the molecular mechanisms underlying the ovarian differentiation and development.

RESULTS

A total of 301,267,988 clean reads were generated from cDNA libraries of gonadal tissues of ricefield eel larvae at 6, 9, 12, and 20 days post hatching (dph), which contained undifferentiated gonads, differentiating ovaries, ovaries with oogonia, and ovaries with meiotic oocytes, respectively. De-novo assembly of all the clean reads generated a total of 265,896 unigenes with a mean size of 720 bp and a N50 of 1107 bp. RT-qPCR analysis of the developmental expression of 13 gonadal development-related functional genes indicated that RNA-seq data are reliable. Transcriptome data suggest that high expression of female development-related genes and low expression of male development-related genes in the early gonads of ricefield eel larvae participate in the cascade of sex differentiation leading to the final female phenotype. The contrasting expression patterns of genes involved in retinoid acid (RA) synthesis and degradation might result in peak production of RA at 12 dph in the gonad of ricefield eel larvae, and induce molecular events responsible for the initiation of meiosis before the meiotic signs could be observed at 20 dph. In addition, only stra6 but not stra8 could be identified in gonadal transcriptome data of ricefield eel larvae, and the expression pattern of stra6 paralleled those of genes involved in RA synthesis, suggesting that stra6 may be a downstream target of RA and play a role in RA metabolism and/or meiotic initiation in the gonad of ricefield eel larvae.

CONCLUSIONS

The present study depicted the first large-scale RNA sequencing of the gonad of ricefield eel larvae, and identified many important functional genes, GO terms and KEGG pathways involved in gonadal development and germ cell meiosis. Results of the present study will facilitate future study on the ovarian differentiation of ricefield eels and other teleosts as well.

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.

Human-mediated secondary contact of two tortoise lineages results in sex-biased introgression

Human-mediated secondary contact of recently diverged taxa offers valuable opportunities for studying the evolutionary mechanisms involved in the establishment and maintenance of genetic boundaries between taxa. We used mitochondrial and microsatellite markers to examine a recently introduced population of the spur-thighed tortoise (Testudo graeca) of mixed origin in the Doñana National Park (SW Spain). The earliest records of tortoises in Doñana trace back to the 18th century, but several population reinforcements in the 20th century with animals from Morocco are well-documented. Consequently, different genetic lineages, which represent distinct subspecies, are thought to co-exist there. Our results confirmed the presence of distinct lineages by revealing that tortoises of the subspecies T. g. marokkensis were introduced into a local allochthonous T. g. graeca population. Unexpectedly, T. g. marokkensis haplotypes exclusively appeared in males, and admixture levels were statistically sex-biased toward males. The sex ratio of the population deviated from parity, with males being 2.36-fold more abundant than females. Our results indicated that population reinforcements had a strong effect on the genetic composition of this population and aggravated its sex ratio deviation. We predict that this sex-biased pattern of introgression is ephemeral and advocated to the near loss of T. g. marokkensis haplotypes.

Evolution of the sex-determining gene in the teleostean genus Oryzias

In the genetic sex determination of vertebrates, the gonadal sex depends on the combination of sex chromosomes that a zygote possesses. Despite the discovery of the sex-determining gene (SRY/Sry) in mammals in 1990s, the sex-determining gene in non-mammalian vertebrates remained an enigma for over a decade. In most mammals, the male-inducing master sex-determining gene is located on the Y chromosome and is therefore absent from XX females. A second sex-determining gene, Dmy, was described in the Oryzias latipes in 2002 and has a DNA-binding motif that is different from the motif in the mammalian sex-determining gene SRY or Sry. Dmy is also located on the Y chromosome and is therefore absent in XX females. Seven other sex-determining genes, including candidate genes, are now known in birds, a frog species, and 5 fish species. These findings over the past twenty years have increased our knowledge of sex-determining genes and sex chromosomes among vertebrates. Here, we review recent advances in our understanding of sex-determining genes and genetic sex determination systems in fish, especially those of the Oryzias species, which are described in detail. The facts suggest some patterns of how new sex-determining genes emerged and evolved. We believe that these facts are common not only in Oryzias but also in other fish species. This knowledge will help to elucidate the conserved mechanisms from which various sex-determining mechanisms have evolved.