Molecular signaling pathways elicited by 17α-ethinylestradiol in Japanese medaka male larvae undergoing gonadal differentiation

Theoretical Analysis and Expression Profiling of 17β-Hydroxysteroid Dehydrogenase Genes in Gonadal Development and Steroidogenesis of Leopard Coral Grouper (Plectropomus leopardus)

The differentiation and developmental trajectory of fish gonads, significantly important for fish breeding, culture, and production, has long been a focal point in the fields of fish genetics and developmental biology. However, the mechanism of gonadal differentiation in leopard coral grouper (Plectropomus leopardus) remains unclear. This study investigates the 17β-Hydroxysteroid Dehydrogenase (Hsd17b) gene family in P. leopardus, with a focus on gene characterization, expression profiling, and functional analysis. The results reveal that the P. leopardus's Hsd17b gene family comprises 11 members, all belonging to the SDR superfamily. The amino acid similarity is only 12.96%, but conserved motifs, such as TGxxxGxG and S-Y-K, are present in these genes. Hsd17b12a and Hsd17b12b are unique homologs in fish, and chromosomal localization has confirmed that they are not derived from different transcripts of the same gene, but rather are two independent genes. The Hsd17b family genes, predominantly expressed in the liver, heart, gills, kidneys, and gonads, are involved in synthesizing or metabolizing sex steroid hormones and neurotransmitters, with their expression patterns during gonadal development categorized into three distinct categories. Notably, Hsd17b4 and Hsd17b12a were highly expressed in the testis and ovary, respectively, suggesting their involvement in the development of reproductive cells in these organs. Fluorescence in situ hybridization (FISH) further indicated specific expression sites for these genes, with Hsd17b4 primarily expressed in germ stem cells and Hsd17b12a in oocytes. This comprehensive study provides foundational insights into the role of the Hsd17b gene family in gonadal development and steroidogenesis in P. leopardus, contributing to the broader understanding of fish reproductive biology and aquaculture breeding.

Norgestrel causes digestive gland injury in the clam Mactra veneriformis: An integrated histological, transcriptomics, and metabolomics study

The potential adverse effects of progestins on aquatic organisms, especially non-target species, are of increasing concern worldwide. However, the effect and mechanism of progestin toxicity on aquatic invertebrates remain largely unexplored. In the present study, clams Mactra veneriformis were exposed to norgestrel (NGT, 0, 10, and 1000 ng/L), the dominant progestin detected in the aquatic environment, for 21 days. NGT accumulation, histology, transcriptome, and metabolome were assessed in the digestive gland. The bioconcentration factor (BCF) was 386 and 268 in the 10 ng/L NGT group and 1000 ng/L NGT group, respectively, indicating efficient accumulation of NGT in the clams. Histological analysis showed that NGT led to the swelling of epithelial cells and blurring of the basement membrane in the digestive gland. Differentially-expressed genes and KEGG pathway enrichment analysis using a transcriptomic approach suggested that NGT primarily disturbed the detoxification system, antioxidant defense, carbohydrate and amino acid metabolism, and steroid hormone metabolism, which was consistent with the metabolites analyzed using a metabolomic approach. Furthermore, we speculated that the oxidative stress caused by NGT resulted in histological damage to the digestive gland. This study showed that NGT caused adverse effects in the clams and sheds light on the mechanisms of progestin interference in aquatic invertebrates.

Tandem Mass Tag-Based Quantitative Proteomics Analysis of Gonads Reveals New Insight into Sexual Reversal Mechanism in Chinese Soft-Shelled Turtles

Chinese soft-shelled turtles display obvious sex dimorphism. The exogenous application of hormones (estradiol and methyltestosterone) can change the direction of gonadal differentiation of P. sinensis to produce sex reversed individuals. However, the molecular mechanism remains unclear. In this study, TMT-based quantitative proteomics analysis of four types of P. sinensis (female, male, pseudo-female, and pseudo-male) gonads were compared. Quantitative analysis of 6107 labeled proteins in the four types of P. sinensis gonads was performed. We identified 440 downregulated and 423 upregulated proteins between pseudo-females and males, as well as 394 downregulated and 959 upregulated proteins between pseudo-males and females. In the two comparisons, the differentially expressed proteins, including K7FKG1, K7GIQ2, COL4A6, K7F2U2, and K7FF80, were enriched in some important pathways, such as focal adhesion, endocytosis, apoptosis, extracellular matrix-receptor interaction, and the regulation of actin cytoskeleton, which were upregulated in pseudo-female vs. male and downregulated in pseudo-male vs. female. In pathways such as ribosome and spliceosome, the levels of RPL28, SRSF3, SNRNP40, and HNRNPK were increased from male to pseudo-female, while they decreased from female to pseudo-male. All differentially expressed proteins after sexual reversal were divided into six clusters, according to their altered levels in the four types of P. sinensis, and associated with cellular processes, such as embryonic development and catabolic process, that were closely related to sexual reversal. These data will provide clues for the sexual reversal mechanism in P. sinensis.

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

Simple Summary

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

Abstract

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

Genome-Wide Association Study Identifies Genomic Loci of Sex Determination and Gonadosomatic Index Traits in Large Yellow Croaker (Larimichthys crocea)

Larimichthys crocea is one of the traditional marine culture fishes in China, widely distributed in South China Sea, East Sea, and southern Yellow Sea. Sex dimorphism is evident in this species that females present a substantial growth strength than males, suggesting breeding females could obtain more economic benefits in L. crocea aquaculture industry. With the continuous expansion of aquaculture industry, both identifying sex-associated genome region and understanding the genetic basis underlying gonad differentiation and development matter to not only sex control aquaculture but also breeding industry. Thus, genome-wide association analysis (GWAS) of sex determination was conducted with a random breeding population of 905 individuals (including 463 females and 442 males) by ddRAD sequencing. For sex determination, 21 significant single nucleotide polymorphisms (SNPs) in chromosome (Chr) 22 were identified. Surrounding these SNPs, we founded 14 candidate genes, including dmrt1, dmrt3, and piwil2, fam102a, and odf2. The sex-associated region was narrowed down further to 2.4 Mb on Chr22 through F_st scanning and insertion-deletion (InDel) analysis. Besides, 3 SNPs in the supposed sex-determining region on Chr22 were identified as highly associated with gonad differentiation through GWAS on gonadosomatic index (GSI) in 350 males and 231 females. Because of the significant difference of GSI between females and males of L. crocea, GWAS on GSI of different genders was also conducted independently. Finally, we identified a SNP in Chr18 showing genome-wide significant association with male GSI (MGSI) and three genes axl, cyp2a10, and cyp2g1 involved in the gonadal development regulation process of aromatase. Overall, this study explored the genetic basis of sex determination mechanism and provided novel insights into gonad differentiation and development, offering solid genetic support for sex control breeding, marker-assisted selection, and marine resources conservation.

Identification of estrogen receptor target genes involved in gonadal feminization caused by estrogen in Xenopus laevis

Estrogens and estrogenic endocrine disrupting chemicals can cause gonadal feminization in some vertebrates mainly through estrogen receptor (ER), but the underlying molecular mechanisms are unclear. The present study aimed to identify ER target genes involved in estrogen-caused gonadal feminization in Xenopus laevis. Based on our recent transcriptomic data that 10 nM 17β-estradiol (E2) altered gene transcription in feminizing gonads of male X. laevis at NF stages 48, 50, and 52, we searched estrogen response element (ERE) using the Dragon ERE Finder software in the promoter region of all the E2-regulated genes. As a result, 163 genes containing ERE sequence were identified as predicted ER target genes at NF stage 50 (on the 14th day postfertilization), a crucial stage for gonadal feminization. Then, some of these predicted ER target genes were further investigated, mainly including the genes that were suggested to be involved in E2-caused gonadal feminization and genes being dramatically up or down-regulated by E2. Fifteen genes were demonstrated to be responsive to E2, in turn ER antagonist blocked the E2-regulated transcription. Finally, we identified 10 genes that can bind to ERα by a chromatin immunoprecipitation-qPCR. Taken together, we identified the 10 genes that contain predicted ERE sequences, are responsive to estrogen and ER antagonist, and have ability to bind to ER as ER target genes, including pglyrp2, apoa1, fgb, tdo2, ca6, nags, cpb2, tmprss6, nudc, zwilch. Our results could help to improve the understanding of the molecular mechanisms for gonadal feminization caused by estrogenic endocrine disrupting chemicals in X. laevis, and even in other species.

Effects of early life stage exposure of largemouth bass to atrazine or a model estrogen (17α-ethinylestradiol)

Endocrine disrupting contaminants are of continuing concern for potentially contributing to reproductive dysfunction in largemouth and smallmouth bass in the Chesapeake Bay watershed (CBW) and elsewhere. Exposures to atrazine (ATR) have been hypothesized to have estrogenic effects on vertebrate endocrine systems. The incidence of intersex in male smallmouth bass from some regions of CBW has been correlated with ATR concentrations in water. Fish early life stages may be particularly vulnerable to ATR exposure in agricultural areas, as a spring influx of pesticides coincides with spawning and early development. Our objectives were to investigate the effects of early life stage exposure to ATR or the model estrogen 17α-ethinylestradiol (EE2) on sexual differentiation and gene expression in gonad tissue. We exposed newly hatched largemouth bass (LMB, Micropterus salmoides) from 7 to 80 days post-spawn to nominal concentrations of 1, 10, or 100 µg ATR/L or 1 or 10 ng EE2/L and monitored histological development and transcriptomic changes in gonad tissue. We observed a nearly 100% female sex ratio in LMB exposed to EE2 at 10 ng/L, presumably due to sex reversal of males. Many gonad genes were differentially expressed between sexes. Multidimensional scaling revealed clustering by gene expression of the 1 ng EE2/L and 100 µg ATR/L-treated male fish. Some pathways responsive to EE2 exposure were not sex-specific. We observed differential expression in male gonad in LMB exposed to EE2 at 1 ng/L of several genes involved in reproductive development and function, including star, cyp11a2, ddx4 (previously vasa), wnt5b, cyp1a and samhd1. Expression of star, cyp11a2 and cyp1a in males was also responsive to ATR exposure. Overall, our results confirm that early development is a sensitive window for estrogenic endocrine disruption in LMB and are consistent with the hypothesis that ATR exposure induces some estrogenic responses in the developing gonad. However, ATR-specific and EE2-specific responses were also observed.

Transcriptomic analysis identifies early cellular and molecular events by which estrogen disrupts testis differentiation and causes feminization in Xenopus laevis

Extensive studies have shown that estrogenic endocrine-disrupting chemicals (EDCs) can disrupt testis differentiation and even cause feminization in vertebrates. However, little is known about the mechanisms by which estrogenic EDCs disrupt testis differentiation. Here, we employed Xenopus laevis, a model amphibian species sensitive to estrogenic EDCs, to explore the molecular and cellular events by which 17β-estradiol (E2) disrupts testis differentiation and causes feminization. Following waterborne exposure to E2 from stage 45/46, genetically male X. laevis were confirmed to undergo testis differentiation inhibition and ovary differentiation activation at stages 52 and 53, ultimately displaying gonadal feminization at stage 66. Using a time-course RNA sequencing approach, we then identified thousands of differentially expressed transcripts (DETs) in genetically male gonad-mesonephros complexes at stages 48, 50 and 52 (the window for testis differentiation) between E2 treatment and the control. Enrichment analysis suggests alterations in cell proliferation, extracellular matrix, and cell motility following E2 exposure. Further verification by multiple methods demonstrated that E2 inhibited cell proliferation, disrupted extracellular matrix, and altered cell motility in the genetically male gonads compared with controls, implying that these events together contributed to testis differentiation disruptions and feminization in X. laevis. This study for the first time uncovered some of the early molecular and cellular events by which estrogen disrupts testicular differentiation and causes feminization in X. laevis. These new findings improve our understanding of the mechanisms by which estrogenic EDCs disrupt testicular differentiation in vertebrates.

Comparative transcriptome analysis on four types of gonadal tissues of blotched snakehead (Channa maculata)

Blotched snakehead (Channa maculata) is an economically important freshwater fish in China, of which males grow much faster than females. To illuminate the molecular mechanism of sex differentiation and gonad development, RNA-Sequencing was performed to identify sex-related genes and pathway in gonads of 6-month-old normal XX females (XX-F), normal XY males (XY-M), XY sex reversal females (XY-F) and YY super-males (YY-M). The analysis showed that many differentially expressed genes (DEGs) had similar expression patterns in XY-F and XX-F, which were different from XY-M and YY-M. qRT-PCR indicated that Amh, Dmrt1, and Sox9 had relatively high expression in testes of XY-M and YY-M. Taking Amh as an example, there was a relative fold change of 1.0 in XX-F, 2.1 fold change in XY-F, 36.1 fold change in XY-M, and 26.0 fold change in YY-M. Cyp19a1a, Figla, and Foxl2 were highly expressive in ovaries of XX-F and XY-F. Taking Figla as an example, there was a relative fold change of 557 in XX-F, 304.5 fold change in XY-F, 5.6 fold change in XY-M, and 4.4 fold change in YY-M. KEGG analysis revealed many DEGs distributed in pathways related to sex differentiation, steroid hormone synthesis and growth, etc. Significant variation and trends in relative expression levels tested by qRT-PCR were consistent with those recorded by RNA-Sequencing. This is the first time that transcriptome of snakehead has been investigated systematically and in an integrated way. Large quantities of candidate genes involved in sex differentiation, gonad development and growth dimorphism were identified. The study provides useful resources for understanding sex differentiation and growth dimorphism, potentially assisting mono-sex production of snakehead in aquaculture.

Twenty years of transcriptomics, 17alpha-ethinylestradiol, and fish

In aquatic toxicology, perhaps no pharmaceutical has been investigated more intensely than 17alpha-ethinylestradiol (EE2), the active ingredient of the birth control pill. At the turn of the century, the fields of comparative endocrinology and endocrine disruption research witnessed the emergence of omics technologies, which were rapidly adapted to characterize potential hazards associated with exposures to environmental estrogens, such as EE2. Since then, significant advances have been made by the scientific community, and as a result, much has been learned about estrogen receptor signaling in fish from environmental xenoestrogens. Vitellogenin, the egg yolk precursor protein, was identified as a major estrogen-responsive gene, establishing itself as the premier biomarker for estrogenic exposures. Omics studies have identified a plethora of estrogen responsive genes, contributing to a wealth of knowledge on estrogen-mediated regulatory networks in teleosts. There have been ~40 studies that report on transcriptome responses to EE2 in a variety of fish species (e.g., zebrafish, fathead minnows, rainbow trout, pipefish, mummichog, stickleback, cod, and others). Data on the liver and testis transcriptomes dominate in the literature and have been the subject of many EE2 studies, yet there remain knowledge gaps for other tissues, such as the spleen, kidney, and pituitary. Inter-laboratory genomics studies have revealed transcriptional networks altered by EE2 treatment in the liver; networks related to amino acid activation and protein folding are increased by EE2 while those related to xenobiotic metabolism, immune system, circulation, and triglyceride storage are suppressed. EE2-responsive networks in other tissues are not as comprehensively defined which is a knowledge gap as regulated networks are expected to be tissue-specific. On the horizon, omics studies for estrogen-mediated effects in fish include: (1) Establishing conceptual frameworks for incorporating estrogen-responsive networks into environmental monitoring programs; (2) Leveraging in vitro and computational toxicology approaches to identify chemicals associated with estrogen receptor-mediated effects in fish (e.g., male vitellogenin production); (3) Discovering new tissue-specific estrogen receptor signaling pathways in fish; and (4) Developing quantitative adverse outcome pathway predictive models for estrogen signaling. As we look ahead, research into EE2 over the past several decades can serve as a template for the array of hormones and endocrine active substances yet to be fully characterized or discovered.

Sex-specific endocrine-disrupting effects of three halogenated chemicals in Japanese medaka

Several halogenated chemicals are found in an array of products that can cause endocrine disruption. Human studies have shown that endocrine responses are sex specific, with females more likely to develop hypothyroidism and males more likely to have reproductive impairment. The objective of this study was to assess sex differences on thyroid and estrogenic effects after exposure of Japanese medaka (Oryzias latipes, SK2MC) to halogenated compounds. This strain is an excellent model for these studies as sex can be determined non-destructively a few hours postfertilization. Medaka embryos were exposed to sublethal concentrations of Tris(1,3-dichloro-2-propyl) phosphate (TDCPP, 0.019 mg/L), perfluorooctanoic acid (PFOA, 4.7 mg/L) and its next generation alternative, perfluorobutyric acid (PFBA, 137 mg/L). Methimazole (inhibits thyroid hormone synthesis) and the thyroid hormone triiodothyronine served as reference controls. Fish were exposed throughout embryo development until 10 days postfertilization. Females displayed significantly larger swim bladders (which are under thyroid hormone control) after exposure to all chemicals with the exception of triiodothyronine, which caused the opposite effect. Females exposed to TDCPP and PFOA had increased expression of vitellogenin and exposure to PFOA upregulated expression of multiple thyroid-related genes. Upregulation of estrogenic-regulated genes after exposure to TDCPP, PFOA and methimazole was only observed in males. Overall, our results suggest that females and males show an estrogenic response when exposed to these halogenated chemicals and that females appear more susceptible to thyroid-induced swim bladder dysfunction compared with males. These results further confirm the importance of considering sex effects when assessing the toxicity of endocrine-disrupting compounds.