Transcriptome Dynamics During Turbot Spermatogenesis Predicting the Potential Key Genes Regulating Male Germ Cell Proliferation and Maturation

Comparative transcriptome analysis identified genes involved in testicular development in Takifugu rubripes

To identify candidate genes and pathways involved in testicular development in Takifugu rubripes, a comparative transcription analysis was conducted across the various developmental stages of the testis (stages II to V). A total of 9520 differentially expressed genes (DEGs) were identified among the different stages, and they were significantly clustered into six clusters (P < 0.05). One thousand four hundred eleven DEGs such as gndf, wnt1, and cyp17b1 were found to be decreased from stage II to V. In contrast, 994 DEGs such as fn1, ift81, and cdc25a were found to be increased from stage II to V. Six thousand three hundred eighteen DEGs (e.g., dmrt1, sdk2, and chrna1) were identified as being expressed at similar levels at stages II and III. However, they were subsequently found to be decreased from stage III to IV. Four hundred one DEGs exhibited a significant upregulation trend from stage II to III. These genes were expressed at similar levels in stages III, IV, and V, including chrnd, wnt4a, and cyp7a1. The highest expression levels of 200 DEGs (e.g., ccnb2, cdk1, and sycp2) were observed in stage IV, while 196 DEGs (e.g., chmp1b, hsd17b3, and zp3) exhibited the highest expression level in stage III. Those DEGs were mainly enriched in the pathways (e.g., neuroactive ligand-receptor interaction, cell adhesion molecules, and calcium signaling pathways) associated with testicular development. Quantitative polymerase chain reaction of eight randomly selected genes validated the RNA sequencing results. This study may provide new insights into the molecular regulatory mechanisms governing testicular development and spermatogenesis in T. rubripes.

Dissecting the dynamic cellular transcriptional atlas of adult teleost testis development throughout the annual reproductive cycle

Teleost testis development during the annual cycle involves dramatic changes in cellular compositions and molecular events. In this study, the testicular cells derived from adult black rockfish at distinct stages - regressed, regenerating and differentiating - were meticulously dissected via single-cell transcriptome sequencing. A continuous developmental trajectory of spermatogenic cells, from spermatogonia to spermatids, was delineated, elucidating the molecular events involved in spermatogenesis. Subsequently, the dynamic regulation of gene expression associated with spermatogonia proliferation and differentiation was observed across spermatogonia subgroups and developmental stages. A bioenergetic transition from glycolysis to mitochondrial respiration of spermatogonia during the annual developmental cycle was demonstrated, and a deeper level of heterogeneity and molecular characteristics was revealed by re-clustering analysis. Additionally, the developmental trajectory of Sertoli cells was delineated, alongside the divergence of Leydig cells and macrophages. Moreover, the interaction network between testicular micro-environment somatic cells and spermatogenic cells was established. Overall, our study provides detailed information on both germ and somatic cells within teleost testes during the annual reproductive cycle, which lays the foundation for spermatogenesis regulation and germplasm preservation of endangered species.

Transcriptome Analysis Reveals Spermatogenesis-Related CircRNAs and LncRNAs in Goat Spermatozoa

Mammalian spermatozoa comprises both coding and non-coding RNAs, which are traditionally believed to be a residual of spermatogenesis. The differential expression level of spermatozoal RNAs is also observed between fertile and infertile human, thereby anticipated as potential molecular marker of male fertility. This study investigated the transcriptome profile of goat (Capra hircus) spermatozoa. The sperm transcriptome was analyzed by three different methods viz. RLM-RACE, long-read RNA sequencing (RNAseq) in Nanopore™ platform, and short-read RNAseq in Illumina™ platform. The Illumina™ sequencing discovered 16,604 transcripts with 357 mRNAs having FPKM (fragments per kilobase per million mapped reads) of more than five. The spermatozoal RNA suite included mRNA (94%), rRNA (3%), miscRNA (1%), circRNA (1%), miRNA (1%), etc. This study also predicted circRNAs (127), lncRNAs (655), and imprinted genes (160) that have potential role in male reproduction. The gene ontology analysis revealed the involvement of spermatozoal RNA in regulating male meiosis (TET3, STAT5B), capacitation (ACRBP, CATSPER4), sperm motility (GAS8, TEKT2), spermatogenesis (ADAMTS2, CREB3L4), etc. The spermatozoal RNA were also associated with different biological pathways viz. Wnt signaling pathway, cAMP signaling pathway, AMPK signaling pathway, and MAPK signaling pathways having potential role in spermatogenesis. Overall, this study enlightened the suite of spRNA transcripts in goat and their relevance in male fertility for diagnostic approach.

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Genome-Wide Identification of Aqp Family Related to Spermatogenesis in Turbot (Scophthalmus maximus)

The development and maturation of sperm entails intricate metabolic processes involving water molecules, amino acids, hormones, and various substances. Among these processes, the role of aquaporins (aqps) in the testis is crucial. Turbot (Scophthalmus maximus) is a significant marine flatfish species in China; however, natural egg laying in females is not feasible under cultured conditions. Consequently, artificial insemination becomes necessary, requiring the retrieval of sperm and eggs through artificial methods. In this study, we combined genomic, transcriptomics, RT-qPCR, computer-assisted sperm analysis (CASA), and immunohistochemistry to investigate the involvement of the aqp family in spermatogenesis in turbot. Through genomic data analysis, we identified 16 aqps genes dispersed across 13 chromosomes, each exhibiting the characteristic major intrinsic protein (MIP) domain associated with AQPs. The results from RNA-seq and RT-qPCR analysis revealed prominent expression of aqp4, 10, and 12 during the proliferative stage, whereas aqp1 showed primary expression during the mature stage. aqp11 displayed high expression levels during both MSII and MSV stages, potentially contributing significantly to the proliferation and maturation of male germ cells. Conversely, aqp8 showed elevated expression levels during the MSIII, MSIII-IV, and MSIV stages, suggesting its direct involvement in spermiogenesis. Immunohistochemical analysis unveiled the predominant localization of AQP1 protein in male germ cells rather than Sertoli cells, specifically concentrated in the head of sperm within cysts. Furthermore, a noteworthy decline in sperm motility was observed when sperm were subjected to treatment with either the AQP1-specific inhibitor (HgCl₂) or the AQP1 antibody. However, no direct correlation was found between the expression of Smaqp1 and sperm quality. Overall, these findings provide new insights into the involvement of aqps in teleost spermatogenesis. Moreover, they hold potential for improving techniques related to sperm activation and cryopreservation, offering valuable knowledge for future advancements in this field.

Dynamic cellular and molecular characteristics of spermatogenesis in the viviparous marine teleost Sebastes schlegelii†

Spermatogenesis is a dynamic cell developmental process that is essential for reproductive success. Vertebrates utilize a variety of reproductive strategies, including sperm diversity, and internal and external fertilization. Research on the cellular and molecular dynamic changes involved in viviparous teleost spermatogenesis, however, is currently lacking. Here, we combined cytohistology, 10 × genomic single-cell RNA-seq, and transcriptome technology to determine the dynamic development characteristics of the spermatogenesis of Sebastes schlegelii. The expressions of lhcgr (Luteinizing hormone/Choriogonadotropin receptor), fshr (follicle-stimulating hormone receptor), ar (androgen receptor), pgr (progesterone receptor), and cox (cyclo-oxygen-ase), as well as the prostaglandin E and F levels peaked during the maturation period, indicating that they were important for sperm maturation and mating. Fifteen clusters were identified based on the 10 × genomic single-cell results. The cell markers of the sub-cluster were identified by their upregulation; piwil, dazl, and dmrt1 were upregulated and identified as spermatogonium markers, and sycp1/3 and spo11 were identified as spermatocyte markers. For S. schlegelii, the sperm head nucleus was elongated (spherical to streamlined in shape), which is a typical characteristic for sperm involved in internal fertilization. We also identified a series of crucial genes associated with spermiogenesis, such as spata6, spag16, kif20a, trip10, and klf10, while kif2c, kifap3, fez2, and spaca6 were found to be involved in nucleus elongation. The results of this study will enrich our cellular and molecular knowledge of spermatogenesis and spermiogenesis in fish that undergo internal fertilization.

Whole Exome Sequencing and In Silico Analysis of Human Sertoli in Patients with Non-Obstructive Azoospermia

Non-obstructive azoospermia (NOA) is a serious cause of male infertility. The Sertoli cell responds to androgens and takes on roles supporting spermatogenesis, which may cause infertility. This work aims to enhance the genetic diagnosis of NOA via the discovery of new and hub genes implicated in human NOA and to better assess the odds of successful sperm extraction according to the individual’s genotype. Whole exome sequencing (WES) was done on three NOA patients to find key genes involved in NOA. We evaluated genome-wide transcripts (about 50,000 transcripts) by microarray between the Sertoli of non-obstructive azoospermia and normal cells. The microarray analysis of three human cases with different non-obstructive azoospermia revealed that 32 genes were upregulated, and the expressions of 113 genes were downregulated versus the normal case. For this purpose, Enrich Shiny GO, STRING, and Cytoscape online evaluations were applied to predict the functional and molecular interactions of proteins and then recognize the master pathways. The functional enrichment analysis demonstrated that the biological process (BP) terms “inositol lipid-mediated signaling”, “positive regulation of transcription by RNA polymerase II”, and “positive regulation of DNA-templated transcription” significantly changed in upregulated differentially expressed genes (DEGs). The BP investigation of downregulated DEGs highlighted “mitotic cytokinesis”, “regulation of protein-containing complex assembly”, “cytoskeleton-dependent cytokinesis”, and the “peptide metabolic process”. Overrepresented molecular function (MF) terms in upregulated DEGs included “ubiquitin-specific protease binding”, “protease binding”, “phosphatidylinositol trisphosphate phosphatase activity”, and “clathrin light chain binding”. Interestingly, the MF analysis of the downregulated DEGs revealed overexpression in “ATPase inhibitor activity”, “glutathione transferase activity”, and “ATPase regulator activity”. Our findings suggest that these genes and their interacting hub proteins could help determine the pathophysiologies of germ cell abnormalities and infertility.

Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.

Comparative transcriptomic analysis of gonadal development and renewal in the ovoviviparous black rockfish (Sebastes schlegelii)

BACKGROUND

The black rockfish (Sebastes schlegelii) has an ovoviviparous reproductive pattern and long-term sperm storage, resulting in asynchronous gonadal development between the sexes. However, the comprehensive understanding of gonadal development in black rockfish has not yet been achieved. Here, we studied gonadal development and germ cell renewal using histology and RNA-seq.

RESULTS

In this study, RNA-seq was performed on testes and ovaries to characterize key pathways and genes that are active during development and gamete maturation in black rockfish. Differentially expressed genes (DEGs) were identified and annotated in 4 comparisons (F_III vs. F_IV, F_IV vs. F_V, M_III vs. M_IV and M_IV vs. M_V). Based on analysis of DEGs enriched in the testis, 11 and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mapped to the M_III vs. M_IV group and the M_IV vs. M_V group, respectively. DEGs in ovarian development were also classified into 10 groups according to their biological functions. The expression patterns of the selected genes determined by qPCR were significantly correlated with the RNA-Seq results, supporting the reliability and accuracy of the RNA-Seq analysis. E2 levels showed down regulation from previtellogenesis to mature stage in female and T level showed down regulation from spermatogenesis to regressed stage in the male.

CONCLUSIONS

The categories "intercellular interaction and cytoskeleton", "molecule amplification" and "repair in the cell cycle" were revealed to be crucial in testis development and spermatogenesis, as was the biosynthesis of a series of metabolites. Our results provide comprehensive insight into black rockfish gonadal development and provide a basis for further study of reproductive physiology and molecular biology in ovoviviparity teleosts.

A single intronic single nucleotide polymorphism in splicing site of steroidogenic enzyme hsd17b1 is associated with phenotypic sex in oyster pompano, Trachinotus anak

Teleosts show varied master sex determining (MSD) genes and sex determination (SD) mechanisms, with frequent turnovers of sex chromosomes. Tracing the origins of MSD genes and turnovers of sex chromosomes in a taxonomic group is of particular interest in evolutionary biology. Oyster pompano (Trachinotus anak), a marine fish, belongs to the family Carangidae, in which 17b-hydroxysteroid dehydrogenase 1 (hsd17b1) has repeatedly evolved to an MSD gene. Whole-genome resequencing identified a single nucleotide polymorphism (SNP) at chromosome 24 to be strictly associated with phenotypic sex, with females being the heterozygous sex. This SNP is located in a splicing site at the first exon/intron boundary of hsd17b1. The Z-linked SNP results in malfunction of all spliced isoforms, whereas the W-linked isoforms were predicted to have open reading frames that are conserved among vertebrates, suggesting that hsd17b1 is a female-determining gene. The differential alternative splicing patterns of ZZ and ZW genotypes were consistently observed both in undifferentiated stages and differentiated gonads. We observed elevated recombination around the SD locus and no differentiation between Z and W chromosomes. The extreme diversity of mutational mechanisms that hsd17b1 evolves to an MSD gene highlights frequent in situ turnovers between sex chromosomes in the Carangidae.

miR-301a-5p Regulates TGFB2 during Chicken Spermatogenesis

The process of spermatogenesis is complex and systemic, requiring the cooperation of many regulators. However, little is known about how micro RNAs (miRNAs) regulate spermatogenesis in poultry. In this study, we investigated key miRNAs and their target genes that are involved in spermatogenesis in chickens. Next-generation sequencing was conducted to determine miRNA expression profiles in five cell types: primordial germ cells (PGCs), spermatogonial stem cells (SSCs), spermatogonia (Spa), and chicken sperm. Next, we analyzed and identified several key miRNAs that regulate spermatogenesis in the four germline cell miRNA profiles. Among the enriched miRNAs, miRNA-301a-5p was the key miRNA in PGCs, SSCs, and Spa. Through reverse transcription quantitative PCR (RT-qPCR), dual-luciferase, and miRNA salience, we confirmed that miR-301a-5p binds to transforming growth factor-beta 2 (TGFβ2) and is involved in the transforming growth factor-beta (TGF-β) signaling pathway and germ cell development. To the best of our knowledge, this is the first demonstration of miR-301a-5p involvement in spermatogenesis by direct binding to TGFβ2, a key gene in the TGF-β signaling pathway. This finding contributes to the insights into the molecular mechanism through which miRNAs regulate germline cell differentiation and spermatogenesis in chickens.

Chronic Exposure to Cadmium Induces Differential Methylation in Mice Spermatozoa

Cadmium exposure is ubiquitous and has been linked to diseases including cancers and reproductive defects. Since cadmium is nonmutagenic, it is thought to exert its gene dysregulatory effects through epigenetic reprogramming. Several studies have implicated germline exposure to cadmium in developmental reprogramming. However, most of these studies have focused on maternal exposure, while the impact on sperm fertility and disease susceptibility has received less attention. In this study, we used reduced representation bisulfite sequencing to comprehensively investigate the impact of chronic cadmium exposure on mouse spermatozoa DNA methylation. Adult male C57BL/J6 mice were provided water with or without cadmium chloride for 9 weeks. Sperm, testes, liver, and kidney tissues were collected at the end of the treatment period. Cadmium exposure was confirmed through gene expression analysis of metallothionein-1 and 2, 2 well-known cadmium-induced genes. Analysis of sperm DNA methylation changes revealed 1788 differentially methylated sites present at regulatory regions in sperm of mice exposed to cadmium compared with vehicle (control) mice. Furthermore, most of these differential methylation changes positively correlated with changes in gene expression at both the transcription initiation stage as well as the splicing levels. Interestingly, the genes targeted by cadmium exposure are involved in several critical developmental processes. Our results present a comprehensive analysis of the sperm methylome in response to chronic cadmium exposure. These data, therefore, highlight a foundational framework to study gene expression patterns that may affect fertility in the exposed individual as well as their offspring, through paternal inheritance.

Plasma proteome profiling of freshwater and seawater life stages of rainbow trout (Oncorhynchus mykiss)

The sea-run phenotype of rainbow trout (Oncorhynchus mykiss), like other anadromous salmonids, present a juvenile stage fully adapted to life in freshwater known as parr. Development in freshwater is followed by the smolt stage, where preadaptations needed for seawater life are developed making fish ready to migrate to the ocean, after which event they become post-smolts. While these three life stages have been studied using a variety of approaches, proteomics has never been used for such purpose. The present study characterised the blood plasma proteome of parr, smolt and post-smolt rainbow trout using a gel electrophoresis liquid chromatography tandem mass spectrometry approach alone or in combination with low-abundant protein enrichment technology (combinatorial peptide ligand library). In total, 1,822 proteins were quantified, 17.95% of them being detected only in plasma post enrichment. Across all life stages, the most abundant proteins were ankyrin-2, DNA primase large subunit, actin, serum albumin, apolipoproteins, hemoglobin subunits, hemopexin-like proteins and complement C3. When comparing the different life stages, 17 proteins involved in mechanisms to cope with hyperosmotic stress and retinal changes, as well as the downregulation of nonessential processes in smolts, were significantly different between parr and smolt samples. On the other hand, 11 proteins related to increased growth in post-smolts, and also related to coping with hyperosmotic stress and to retinal changes, were significantly different between smolt and post-smolt samples. Overall, this study presents a series of proteins with the potential to complement current seawater-readiness assessment tests in rainbow trout, which can be measured non-lethally in an easily accessible biofluid. Furthermore, this study represents a first in-depth characterisation of the rainbow trout blood plasma proteome, having considered three life stages of the fish and used both fractionation alone or in combination with enrichment methods to increase protein detection.

The initiation of puberty in Atlantic salmon brings about large changes in testicular gene expression that are modulated by the energy status

Background

When puberty starts before males reach harvest size, animal welfare and sustainability issues occur in Atlantic salmon (Salmo salar) aquaculture. Hallmarks of male puberty are an increased proliferation activity in the testis and elevated androgen production. Examining transcriptional changes in salmon testis during the transition from immature to maturing testes may help understanding the regulation of puberty, potentially leading to procedures to modulate its start. Since differences in body weight influence, via unknown mechanisms, the chances for entering puberty, we used two feed rations to create body weight differences.

Results

Maturing testes were characterized by an elevated proliferation activity of Sertoli cells and of single undifferentiated spermatogonia. Pituitary gene expression data suggest increased Gnrh receptor and gonadotropin gene expression, potentially responsible for the elevated circulating androgen levels in maturing fish. Transcriptional changes in maturing testes included a broad variety of signaling systems (e.g. Tgfβ, Wnt, insulin/Igf, nuclear receptors), but also, activation of metabolic pathways such as anaerobic metabolism and protection against ROS. Feed restriction lowered the incidence of puberty. In males maturing despite feed restriction, plasma androgen levels were higher than in maturing fish receiving the full ration. A group of 449 genes that were up-regulated in maturing fully fed fish, was up-regulated more prominently in testis from fish maturing under caloric restriction. Moreover, 421 genes were specifically up-regulated in testes from fish maturing under caloric restriction, including carbon metabolism genes, a pathway relevant for nucleotide biosynthesis and for placing epigenetic marks.

Conclusions

Undifferentiated spermatogonia and Sertoli cell populations increased at the beginning of puberty, which was associated with the up-regulation of metabolic pathways (e.g. anaerobic and ROS pathways) known from other stem cell systems. The higher androgen levels in males maturing under caloric restriction may be responsible for the stronger up-regulation of a common set of (449) maturation-associated genes, and the specific up-regulation of another set of (421) genes. The latter opened regulatory and/or metabolic options for initiating puberty despite feed restriction. As a means to reduce the incidence of male puberty in salmon, however, caloric restriction seems unsuitable.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5869-9) contains supplementary material, which is available to authorized users.

A SNP in a Steroidogenic Enzyme Is Associated with Phenotypic Sex in Seriola Fishes

Vertebrate sex development consists largely of two processes: "sex determination," the initial bifurcation of sexual identity, and "sex differentiation," which subsequently facilitates maleness or femaleness according to the sex determination signal. Steroid hormones promote multiple types of sexual dimorphism in eutherian mammals and avians [1-3], in which they are indispensable for proper sex differentiation. By contrast, in many poikilothermic vertebrates, steroid hormones have been proposed to be key players in sex determination as well as sex differentiation [4-8]. This hypothesis was introduced more than 50 years ago but has never been rigorously tested due to difficulties in discriminating the roles of steroids in sex determination and differentiation. We found that a missense SNP in the gene encoding the steroidogenic enzyme 17β-hydroxysteroid dehydrogenase 1 (Hsd17b1) was perfectly associated with ZZ/ZW sex determination in Seriola fishes. Biochemical analyses revealed that a glutamate residue present specifically in Z-type HSD17B1 attenuated interconversion between 17-keto and 17β-hydroxy steroids relative to the allelic product from the W chromosome, which harbors glycine at that position, by disrupting the hydrogen bond network between the steroid and the enzyme's catalytic residues. Hsd17b1 mRNA is constitutively expressed in undifferentiated and differentiating gonads of both genotypic sexes, whereas W-type mRNA is expressed only in genotypic females. Meanwhile, Cyp19a1 is predominantly expressed in differentiating ovary. We conclude that the combination of Hsd17b1 alleles determines sex by modulating endogenous estrogen levels in Seriola species. These findings strongly support the long-standing hypothesis on steroids in sex determination.