Histological and transcriptomic effects of 17α-methyltestosterone on zebrafish gonad development

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.

Sub-acute exposure of male guppies (Poecilia reticulata) to environmentally relevant concentrations of PFOA and GenX induces significant changes in the testis transcriptome and reproductive traits

Poly- and perfluoroalkyl substances (PFAS) are frequently detected in the environment and are linked to adverse reproductive health outcomes in humans. Although legacy PFAS have been phased out due to their toxicity, alternative PFAS are increasingly used despite the fact that information on their toxic effects on reproductive traits is particularly scarce. Here, we exposed male guppies (Poecilia reticulata) for a short period (21 days) to an environmentally realistic concentration (1 ppb) of PFOA, a legacy PFAS, and its replacement compound, GenX, to assess their impact on reproductive traits and gene expression. Exposure to PFAS did not impair survival but instead caused sublethal effects. Overall, PFAS exposure caused changes in male sexual behaviour and had detrimental effects on sperm motility. Sublethal variations were also seen at the transcriptional level, with the modulation of genes involved in immune regulation, spermatogenesis, and oxidative stress. We also observed bioaccumulation of PFAS, which was higher for PFOA than for GenX. Our results offer a comprehensive comparison of these two PFAS and shed light on the toxicity of a newly emerging alternative to legacy PFAS. It is therefore evident that even at low concentrations and with short exposure, PFAS can have subtle yet significant effects on behaviour, fertility, and immunity. These findings underscore the potential ramifications of pollution under natural conditions and their impact on fish populations.

Effects of exogenous steroid hormones on growth, body color, and gonadal development in the Opsariichthys bidens

To investigate the effects of exogenous steroid hormones on growth, body color, and gonadal development in the Opsariichthys bidens (O. bidens), synthetic methyltestosterone (MT) and 17β-estradiol (E2) were used for 28 days' treatment of 4-month-old O. bidens before the breeding season. Our results suggested that MT had a significant growth-promoting effect (P < 0.05), whereas E2 played an inhibitory role. On the body surface, the females in the MT group showed gray stripes, and the fish in other groups showed no obvious stripes. The males with MT treatment displayed brighter blue-green stripes compared to the CK and E2 groups. The histological analysis showed that the MT significantly promoted testes development in males, blocked oocyte development, and caused massive apoptosis in females, whereas the E2 group promoted ovarian development and inhibited testes development. Based on qRT-PCR analysis, in females, the expression of igf-1, dmrt1, and cyp19a1a genes revealed that E2 treatment resulted in down-regulation of igf-1 expression and up-regulation of cyp19a1a expression. In males, igf-1 and dmrt1 were significantly up-regulated after MT treatment, and E2 treatment led to down-regulation of igf-1. Therefore, this study demonstrates that MT and E2 play an important role in reversing the morphological sex characteristics of females and males.

Interaction of 17α-ethinylestradiol and methyltestosterone in western mosquitofish (Gambusia affinis) across two generations

The interactions between estrogen and androgen in aquatic animals remain largely unknown. In this study, two generations (F0 and F1) of western mosquitofish (Gambusia affinis) were continuously exposed to 17α-ethinylestradiol (EE2, 10 ng/L), methyltestosterone (MT, 10 ng/L (MTL); 50 ng/L (MTH)), and mixtures (EE2+MTL and EE2+MTH). Various endpoints, including sex ratio (phenotypic and genetic), secondary sex characteristics, gonadal histology, and transcriptional profile of genes, were examined. The results showed that G. affinis exposed to MTH and EE2+MTH had a > 89.7 % of phenotypic males in F1 generation, with 34.5 and 50.0 % of these males originated from genetic females, respectively. Moreover, females from F0 and F1 generations exposed to MTH and EE2+MTH exhibited masculinized anal fins and skeletons. The combined effect of MT and EE2 on most endpoints was dependent on MT. Furthermore, significant transcriptional alterations in certain target genes were observed in both the F0 and F1 generations by EE2 and MT alone and by mixtures, showing some degree of interactions. These findings that the effects of EE2+MTH were primarily on the phenotypic sex of G. affinis in offspring generation suggest that G. affinis under chronic exposure to the binary mixture contaminated water could have sex-biased populations.

Identification of Insulin-like Growth Factor (IGF) Family Genes in the Golden Pompano, Trachinotus ovatus: Molecular Cloning, Characterization and Gene Expression

Insulin-like growth factors (IGFs) are hormones that primarily stimulate and regulate animal physiological processes. In this study, we cloned and identified the open reading frame (ORF) cDNA sequences of IGF family genes: the insulin-like growth factor 1 (IGF1), insulin-like growth factor 2 (IGF2), and insulin-like growth factor 3 (IGF3). We found that IGF1, IGF2, and IGF3 have a total length of 558, 648, and 585 base pairs (bp), which encoded a predicted protein with 185, 215, and 194 amino acids (aa), respectively. Multiple sequences and phylogenetic tree analysis showed that the mature golden pompano IGFs had been conserved and showed high similarities with other teleosts. The tissue distribution experiment showed that IGF1 and IGF2 mRNA levels were highly expressed in the liver of female and male fish. In contrast, IGF3 was highly expressed in the gonads and livers of male and female fish, suggesting a high influence on fish reproduction. The effect of fasting showed that IGF1 and mRNA expression had no significant difference in the liver but significantly decreased after long-term (7 days) fasting in the muscles and started to recover after refeeding. IGF2 mRNA expression showed no significant difference in the liver but had a significant difference in muscles for short-term (2 days) and long-term fasting, which started to recover after refeeding, suggesting muscles are more susceptible to both short-term and long-term fasting. In vitro incubation of 17β-estradiol (E2) was observed to decrease the IGF1 and IGF3 mRNA expression level in a dose- (0.1, 1, and 10 μM) and time- (3, 6, and 12 h) dependent manner. In addition, E2 had no effect on IGF2 mRNA expression levels in a time- and dose-dependent manner. The effect of 17α-methyltestosterone (MT) in vitro incubation was observed to significantly increase the IGF3 mRNA expression level in a time- and dose-dependent manner. MT had no effect on IGF2 mRNA but was observed to decrease the IGF1 mRNA expression in the liver. Taken together, these data indicate that E2 and MT may either increase or decrease IGF expression in fish; this study provides basic knowledge and understanding of the expression and regulation of IGF family genes in relation to the nutritional status, somatic growth, and reproductive endocrinology of golden pompano for aquaculture development.

Exploring the Effects of Rearing Densities on Epigenetic Modifications in the Zebrafish Gonads

Rearing density directly impacts fish welfare, which, in turn, affects productivity in aquaculture. Previous studies have indicated that high-density rearing during sexual development in fish can induce stress, resulting in a tendency towards male-biased sex ratios in the populations. In recent years, research has defined the relevance of the interactions between the environment and epigenetics playing a key role in the final phenotype. However, the underlying epigenetic mechanisms of individuals exposed to confinement remain elucidated. By using zebrafish (Danio rerio), the DNA methylation promotor region and the gene expression patterns of six genes, namely dnmt1, cyp19a1a, dmrt1, cyp11c1, hsd17b1, and hsd11b2, involved in the DNA maintenance methylation, reproduction, and stress were assessed. Zebrafish larvae were subjected to two high-density conditions (9 and 66 fish/L) during two periods of overlapping sex differentiation of this species (7 to 18 and 18 to 45 days post-fertilization, dpf). Results showed a significant masculinization in the populations of fish subjected to high densities from 18 to 45 dpf. In adulthood, the dnmt1 gene was differentially hypomethylated in ovaries and its expression was significantly downregulated in the testes of fish exposed to high-density. Further, the cyp19a1a gene showed downregulation of gene expression in the ovaries of fish subjected to elevated density, as previously observed in other studies. We proposed dnmt1 as a potential testicular epimarker and the expression of ovarian cyp19a1a as a potential biomarker for predicting stress originated from high densities during the early stages of development. These findings highlight the importance of rearing densities by long-lasting effects in adulthood conveying cautions for stocking protocols in fish hatcheries.

Germ cells do not progress through spermatogenesis in the infertile zebrafish testis

Vertebrate spermatogonial stem cells maintain sperm production over the lifetime of an animal but fertility declines with age. While morphological studies have greatly informed our understanding of typical spermatogenesis, the molecular and cellular mechanisms underlying spermatogenesis are not yet understood, particularly with respect to the onset of fertility. We used single-cell RNA sequencing to generate a developmental atlas of the zebrafish testis. Using 5 timepoints across the adult life of a zebrafish, we described cellular profiles in the testis during and after fertility. While all germ cell stages of spermatogenesis are detected in testes from fertile adult zebrafish, testes from older infertile males only contained spermatogonia and a reduced population of spermatocytes. These remaining germ cells are transcriptionally distinct from fertile spermatogonia. Immune cells including macrophages and lymphocytes drastically increase in abundance in infertile testes. Our developmental atlas reveals the cellular changes as the testis ages and defines a molecular roadmap for the regulation of male fertility.

Comparative Analysis of miRNA-mRNA Regulation in the Testes of Gobiocypris rarus following 17α-Methyltestosterone Exposure

17α-Methyltestosterone (17MT), a synthetic organic compound commonly found in sewage waters, can affect reproduction in aquatic animals, such as tilapia and yellow catfish. In the present study, male Gobiocypris rarus were exposed to 25, 50, and 100 ng/L of 17α-methyltestosterone (17MT) for 7 days. We first analyzed miRNA- and RNA-seq results to determine miRNA-target gene pairs and then developed miRNA-mRNA interactive networks after 17MT administration. Total weights, total lengths, and body lengths were not significantly different between the test groups and control groups. The paraffin slice method was applied to testes of G. rarus in the MT exposure and control groups. We found that there were more mature sperm (S) and fewer secondary spermatocytes (SSs) and spermatogonia (SGs) in the testes of control groups. As 17MT concentration increased, fewer and fewer mature sperm (S) were observed in the testes of male G. rarus. The results showed that FSH, 11-KT, and E2 were significantly higher in individuals exposed to 25 ng/L 17MT compared with the control groups. VTG, FSH, LH, 11-KT, and E2 were significantly lower in the 50 ng/L 17MT exposure groups compared to the control groups. VTG, FSH, LH, 11-KT, E2, and T were significantly lower in the groups exposed to 100 ng/L 17MT. High-throughput sequencing revealed 73,449 unigenes, 1205 known mature miRNAs, and 939 novel miRNAs in the gonads of G. rarus. With miRNA-seq, 49 (MT25-M vs. Con-M), 66 (MT50-M vs. Con-M), and 49 (MT100-M vs. Con-M) DEMs were identified in the treatment groups. Five mature miRNAs (miR-122-x, miR-574-x, miR-430-y, lin-4-x, and miR-7-y), as well as seven differentially expressed genes (soat2, inhbb, ihhb, gatm, faxdc2, ebp, and cyp1a1), which may be associated with testicular development, metabolism, apoptosis, and disease response, were assayed using qRT-PCR. Furthermore, miR-122-x (related to lipid metabolism), miR-430-y (embryonic development), lin-4-x (apoptosis), and miR-7-y (disease) were differentially expressed in the testes of 17MT-exposed G. rarus. This study highlights the role of miRNA-mRNA pairs in the regulation of testicular development and immune response to disease and will facilitate future studies on the miRNA-RNA-associated regulation of teleost reproduction.

Bisphenol A alters sexual dimorphism and gene expression in marine medaka Oryzias melastigma

Bisphenol A (BPA) is an endocrine disruptor that is present in freshwater and marine environments. However, conclusive evidence for the toxicity of chronic BPA exposure to marine fishes remains lacking. Therefore, we investigated the influence of BPA on male marine medaka (Oryzias melastigma). BPA exposure induced formation of testis-ova at 2610 µg/L, and male-type anal fins became more female type in a concentration-dependent manner. Some males with female-type anal fins had normal testes, indicating that anal fin shape is more sensitive to BPA. Gonadal soma-derived factor (gsdf) expression decreased after BPA exposure in the 746 and 2610 µg/L exposure groups, although the changes were not statistically significant. Additionally, liver vitellogenin (vtg) expression increased in a dose-dependent manner and was significantly higher in all exposure groups. vtg and gsdf are likely to be useful biomarkers for the impact of estrogenic endocrine disrupters in O. melastigma.

A Simple Method for Inducing Masculinization of Zebrafish Stocks Using 17α-Methyltestosterone

Severely skewed sex ratios in zebrafish stocks can pose significant hurdles for line propagation and sperm cryopreservation. To overcome female-biased sex ratios in stocks derived from imported sperm samples, the Zebrafish International Resource Center has implemented routine supplementation of larval food with 17α-methyltestosterone to skew gonadal sex differentiation toward masculinization. Resulting stocks averaged 80% males.

Environmentally-induced sex reversal in fish with chromosomal vs. polygenic sex determination

Sex ratio depends on sex determination mechanisms and is a key demographic parameter determining population viability and resilience to natural and anthropogenic stressors. There is increasing evidence that the environment can alter sex ratio even in genetically sex-determined species (GSD), as elevated temperature can cause female-to-male sex reversal (neomales). Alarmingly, neomales are being discovered in natural populations of several fish, amphibian and reptile species worldwide. Understanding the basis of neomale development is important for conservation biology. Among GSD species, it is unknown whether those with chromosomal sex determination (CSD), the most common system, will better resist the influence of high temperature than those with polygenic sex determination (PSD). Here, we compared the effects of elevated temperature in two wild zebrafish strains, Nadia (NA) and Ekkwill (EKW), which have CSD with a ZZ/ZW system, against the AB laboratory strain, which has PSD. First, we uncovered novel sex genotypes and the results showed that, at control temperature, the masculinization rate roughly doubled with the addition of each Z chromosome, while some ZW and WW fish of the wild strains became neomales. Surprisingly, we found that at elevated temperatures WW fish were just as likely as ZW fish to become neomales and that all strains were equally susceptible to masculinization. These results demonstrate that the Z chromosome is not essential for male development and that the dose of W buffers masculinization at the control temperature but not at elevated temperature. Furthermore, at the elevated temperature the testes of neomales, but not of normal males, contained more spermatozoa than at the control temperature. Our results show in an unprecedented way that, in a global warming scenario, CSD species may not necessarily be better protected against the masculinizing effect of elevated temperature than PSD species, and reveal genotype-by-temperature interactions in male sex determination and spermatogenesis.

Gonadal Soma-Derived Factor Expression is a Potential Biomarker for Predicting the Effects of Endocrine-Disrupting Chemicals on Gonadal Differentiation in Japanese Medaka (Oryzias Latipes)

Chemicals with androgenic or estrogenic activity induce the sex reversal and/or intersex condition in various teleost fish species. Previously, we reported that exposure to 17α-methyltestosterone, bisphenol A, or 4-nonylphenol induces changes in expression of the gonadal soma-derived factor (gsdf) gene accompanied by disruption of gonadal differentiation in Japanese medaka (Oryzias latipes). These findings suggest that gsdf expression might be a useful biomarker for predicting the potential effect of chemicals on gonadal differentiation. We examined the gsdf expression in Japanese medaka exposed to chemicals with estrogenic or androgenic activity. Exposure to the androgenic steroid 17β-trenbolone at 0.5-22.1 μg/L induced the development of ovotestis (presence of ovarian tissue with testicular tissue) and female-to-male sex reversal in XX embryos, and exposure at 6.32 and 22.1 μg/L significantly increased gsdf expression in XX embryos compared with controls at developmental stage 38 (1 day before hatching). In the present study, no statistically significant difference in gsdf mRNA expression was observed after exposure to 17β-estradiol, 17α-ethinylestradiol, and 4-t-octylphenol, which have estrogenic activity. In addition, antiandrogenic chemicals or chemicals without endocrine-disrupting activity did not induce changes in gsdf expression in XX or XY embryos. Thus, an increase in gsdf expression after androgen exposure was observed in XX embryos. Together, these findings indicate that gsdf expression might be useful for predicting the adverse effect of chemicals on gonadal differentiation. Environ Toxicol Chem 2022;41:1875-1884. © 2022 SETAC.

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 Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Transgenerational effects of androstadienedione and androstenedione at environmentally relevant concentrations in zebrafish (Danio rerio)

Androgens androstadienedione (ADD) and androstenedione (AED) are predominant steroid hormones in surface water, and can disrupt the endocrine system in fish. However, little is known about the transgenerational effects of ADD and AED in fish. In the present study, F0 generation was exposed to ADD and AED from 21 to 144 days post-fertilization (dpf) at nominal concentrations of 5 (L), 50 (M) and 500 (H) ng L^-1, and F1 generation was domesticated in clear water for 144 dpf. The sex ratio, histology and transcription in F0 and F1 generations were examined. In the F0 generation, ADD and AED tended to be estrogenic in zebrafish, resulting in female biased zebrafish populations. In the F1 generation, ADD at the H level caused 63.5% females, while AED at the H level resulted in 78.7% males. In brain, ADD and AED had similar effects on circadian rhythm in the F0 and F1 generations. In the F1 eleutheroembryos, transcriptomic analysis indicated that neuromast hair cell related biological processes (BPs) were overlapped in the ADD and AED groups. Taken together, ADD and AED at environmentally relevant concentrations had transgenerational effects on sex differentiation and transcription in zebrafish.

17β-Trenbolone binds to androgen receptor, decreases number of primordial germ cells, modulates expression of genes related to sexual differentiation, and affects sexual differentiation in zebrafish (Danio rerio)

Exposure to 17β-trenbolone caused a skewed sex ratio in fish. However, the molecular initiating event and key molecular event(s) remain unknown. In this study, zebrafish were exposed to 17β-trenbolone at nominal concentrations of 2 ng/L, 20 ng/L, 200 ng/L, and 2000 ng/L from fertilization to 60 days post fertilization (dpf). First, the sex ratio at 60 dpf was calculated to evaluate adverse outcomes on sexual differentiation. 17β-Trenbolone caused a skewed sex ratio toward males, with intersex individuals observed in the 20 ng/L group and all-male populations found in the 200 ng/L and 2000 ng/L groups. Then, the distribution and number of primordial germ cells, the expression of sex differentiation-related genes, and plasma vitellogenin concentrations were detected in wild-type zebrafish and the EGFP-nanos-3'UTR transgenic line using whole-mount in situ hybridization, real-time PCR, EGFP fluorescence quantification, and enzyme-linked immunosorbent assay. The results indicated that 17β-trenbolone exposure decreased the number of primordial germ cells at 1 dpf and 3 dpf, decreased expression of ovarian differentiation-related genes foxl2 and cyp19a1a at 60 dpf, increased expression of testis differentiation-related genes dmrt1, sox9a, and amh at 60 dpf, and decreased plasma vitellogenin levels at 60 dpf, revealing the key molecular events at different time points involved in affected sexual differentiation by 17β-trenbolone. Finally, molecular docking showed that 17β-trenbolone docked into ligand-binding domain of zebrafish androgen receptor with high binding energy (-3.72 kcal/mol), suggesting that binding to androgen receptor is the molecular initiating event affecting sexual differentiation by 17β-trenbolone. We found that 17β-trenbolone can bind to the zebrafish androgen receptor, decrease the number of primordial germ cells during the early embryonic stage, modulate the expression of genes related to sexual differentiation during gonadal differentiation, and eventually cause a skewed sex ratio toward males in adults.

Gene Expression of Takifugu rubripes Gonads During AI- or MT-induced Masculinization and E2-induced Feminization

Elucidating the global molecular changes that occur during aromatase inhibitor (AI)- or 17α-methyltestosterone (MT)-induced masculinization and estradiol-17β (E2)-induced feminization is critical to understanding the roles that endocrine and genetic factors play in regulating the process of sex differentiation in fish. Here, fugu larvae were treated with AI (letrozole), MT, or E2 from 25 to 80 days after hatching (dah), and gonadal transcriptomic analysis at 80 dah was performed. The expression of dmrt1, gsdf, foxl2, and other key genes (star, hsd3b1, cyp11c1, cyp19a1a, etc.) involved in the steroid hormone biosynthesis pathway were found be altered. The expression of dmrt1, gsdf, cyp19a1a, and foxl2 was further verified by quantitative polymerase chain reaction. In the control group, the expression of dmrt1 and gsdf was significantly higher in XY larvae than in XX larvae, while the expression of foxl2 and cyp19a1a was significantly higher in XX larvae than in XY larvae (P < .05). AI treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 and gsdf in XX larvae. MT treatment suppressed the expression of foxl2, cyp19a1a, dmrt1, and gsdf in XX larvae. E2 treatment suppressed the expression of dmrt1 and gsdf, but did not restore the expression of foxl2 and cyp19a1a in XY larvae. The shared response following AI, MT, and E2 treatment suggested that these genes are essential for sex differentiation. This finding offers some insight into AI or MT-induced masculinization, and E2-induced femininization in fugu.

Integrated Analysis of mRNA- and miRNA-Seq in the Ovary of Rare Minnow Gobiocypris rarus in Response to 17α-Methyltestosterone

17α-Methyltestosterone (MT) is a synthetic androgen. The objective of this study was to explore the effects of exogenous MT on the growth and gonadal development of female rare minnow Gobiocypris rarus. Female G. rarus groups were exposed to 25–100 ng/L of MT for 7 days. After exposure for 7 days, the total weight and body length were significantly decreased in the 50-ng/L MT groups. The major oocytes in the ovaries of the control group were vitellogenic oocytes (Voc) and cortical alveolus stage oocytes (Coc). In the MT exposure groups, some fish had mature ovaries with a relatively lower proportion of mature oocytes, and the diameter of the perinucleolar oocytes (Poc) was decreased compared with those of the control group. Ovarian VTG, FSH, LH, 11-KT, E2, and T were significantly increased after exposure to 50 ng/L of MT for 7 days. Unigenes (73,449), 24 known mature microRNAs (miRNAs), and 897 novel miRNAs in the gonads of G. rarus were found using high-throughput sequencing. Six mature miRNAs (miR-19, miR-183, miR-203, miR-204, miR-205, and miR-96) as well as six differentially expressed genes (fabp3, mfap4, abca1, foxo3, tgfb1, and zfp36l1) that may be associated with ovarian development and innate immune response were assayed using qPCR. Furthermore, the miR-183 cluster and miR-203 were differentially expressed in MT-exposed ovaries of the different G. rarus groups. This study provides some information about the role of miRNA–mRNA pairs in the regulation of ovarian development and innate immune system, which will facilitate future studies of the miRNA–RNA-associated regulation of teleost reproduction.

Transcriptome analysis of five ovarian stages reveals gonad maturation in female Macrobrachium nipponense

Background

Macrobrachium nipponense is an economically important species of freshwater shrimp in China. Unlike other marine shrimps, the ovaries in adult female M. nipponense can mature rapidly and periodically during the reproductive period, but the resulting high stocking densities and environmental deterioration can negatively impact the harvest yield and economic benefits. To better understand ovary development in female M. nipponense, we performed systematic transcriptome sequencing of five different stages of ovarian maturation.

Results

We obtained 255,966 Gb of high quality transcriptome data from 15 samples. Of the 105,082 unigenes that were selected, 30,878 were successfully annotated. From these unigenes, we identified 17 differentially expressed genes and identified three distinct gene expression patterns related to different biological processes. We found that cathepins, legumains, and cystatin were enriched in the lysosome pathway, and they are related to vitellogenin hydrolysis. Additionally, we found that myosin heavy chain 67 participated in oocyte excretion.

Conclusions

We provide the first detailed transcriptome data relating to the ovarian maturation cycle in M. nipponense. Our results provide important reference information about the genomics, molecular biology, physiology, and population genetics of M. nipponense and other crustaceans. It is conducive to further solve the problem of M. nipponense rapid ovarian maturation from the aspects of energy supply and cell division.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07737-5.

Differential expression analysis and identification of sex-related genes by gonad transcriptome sequencing in estradiol-treated and non-treated Ussuri catfish Pseudobagrus ussuriensis

The Ussuri catfish (Pseudobagrus ussuriensis) has an XX/XY sex determination system but its sex determination gene(s) remain unknown. To better understand the molecular sex determination mechanism, transcriptome analysis was conducted to obtain sex-related gene expression profiles. Transcriptome analyses were made of male and female developing/differentiating gonads by high-throughput RNA sequencing, including gonads from fish given an estradiol-induced sex reversal treatment. A total of 81,569 unigenes were assembled and 39,904 were significantly matched to known unique proteins by comparison with public databases. Twenty specifically expressed and 142 differentially expressed sex-related genes were extracted from annotated data by comparing the treatment groups. These genes are involved in spermatogenesis (e.g., Dnali1, nectin3, klhl10, mybl1, Katnal1, Eno4, Mns1, Spag6, Tsga10, Septin7), oogenesis (e.g., Lagr5, Fmn2, Npm2, zar1, Fbxo5, Fbxo43, Prdx4, Nrip1, Lfng, Atrip), gonadal development/differentiation (e.g., Cxcr4b, Hmgb2, Cftr, Ch25h, brip1, Prdm9, Tdrd1, Star, dmrt1, Tut4, Hsd17b12a, gdf9, dnd, arf1, Spata22), and estradiol response (e.g., Mmp14, Lhcgr, vtg1, vtg2, esr2b, Piwil1, Aifm1, Hsf1, gdf9). Dmrt1 and gdf9 may play an essential role in sex determination in P. ussuriensis. The expression patterns of six random genes were validated by quantitative real-time PCR, which confirmed the reliability and accuracy of the RNA-seq results. These data provide a valuable resource for future studies of gene expression and for understanding the molecular mechanism of sex determination/differentiation and gonadal development/differentiation (including hormone-induced sexual reversal) in Ussuri catfish. This has the potential to assist in producing monosex Ussuri catfish to increase aquacultural productivity.

Transcriptomic and metabolomic insights into the variety of sperm storage in oviduct of egg layers

In birds, the sperm storage tubules (SST) are dispersed in uterovaginal junction (UVJ) and highly correlated with differential capacity of sperm storage (SS) in and among species with unspecified mechanisms. Here, the SS duration of 252 egg layer breeders was evaluated in 5 rounds with 3 phenotypic traits to screen high- and low-SS individuals, respectively, followed with transcriptome of UVJ tissues and metabolome of serum (high-SS vs. low-SS) to decipher the candidate genes and biochemical markers correlated with differential SS capacity. Histological characterization suggested slightly higher density of SST in UVJ (high-SS vs. low-SS). Transcriptome analyses identified 596 differentially expressed genes (336 upregulated vs. 260 downregulated), which were mainly enriched in gene ontology terms of homeostasis, steroid and lipid metabolism and hormone activity, and 12 significant pathways (P < 0.05) represented by calcium, steroid, and lipid metabolism. Immunohistochemical staining of GNAQ, ST6GAL1, ADFP, and PCNA showed similar distribution in UVJ tissues between 2 groups. Several candidates (HSD11B2, DIO2, AQP3, GNAQ, NANS, ST6GAL1) combined with 4 (11β-prostaglandin F2α, prostaglandin B1, 7α-hydroxytestosterone, and N-acetylneuraminic acid) of 40 differential metabolites enriched in serum metabolome were considered as regulators and biomarkers of SS duration in egg layer breeders. The integrated transcriptome and metabolome analyses of chicken breeder hens will provide novel insights for exploration and improvement of differential SS capacity in birds.

Comparative transcriptome analysis of three gonadal development stages reveals potential genes involved in gametogenesis of the fluted giant clam (Tridacna squamosa)

BACKGROUND

Gonad development and differentiation is an essential function for all sexually reproducing species, and many aspects of these developmental processes are highly conserved among the metazoa. However, the mechanisms underlying gonad development and gametogenesis remain unclear in Tridacna squamosa, a large-size bivalve of great ecological value. They are protandrous simultaneous hermaphrodites, with the male gonad maturing first, eventually followed by the female gonads. In this study, nine gonad libraries representing resting, male and hermaphrodite stages in T. squamosa were performed to identify the molecular mechanisms.

RESULTS

Sixteen thousand four hundred ninety-one unigenes were annotated in the NCBI non-redundant protein database. Among the annotated unigenes, 5091 and 7328 unigenes were assigned to Gene Ontology categories and the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway database, respectively. A total of 4763 differentially expressed genes (DEGs) were identified by comparing male to resting gonads, consisting of 3499 which were comparatively upregulated in males and 1264 which were downregulated in males. Six hundred-ninteen DEGs between male and hermaphroditic gonads were identified, with 518 DEGs more strongly expressed in hermaphrodites and 101 more strongly expressed in males. GO (Gene Ontology) and KEGG pathway analyses revealed that various biological functions and processes, including functions related to the endocrine system, oocyte meiosis, carbon metabolism, and the cell cycle, were involved in regulating gonadal development and gametogenesis in T. squamosa. Testis-specific serine/threonine kinases 1 (TSSK1), TSSK4, TSSK5, Doublesex- and mab-3-related transcription factor 1 (DMRT1), SOX, Sperm surface protein 17 (SP17) and other genes were involved in male gonadal development in Tridacna squamosal. Both spermatogenesis- (TSSK4, spermatogenesis-associated protein 17, spermatogenesis-associated protein 8, sperm motility kinase X, SP17) and oogenesis-related genes (zona pellucida protein, Forkhead Box L2, Vitellogenin, Vitellogenin receptor, 5-hydroxytryptamine, 5-hydroxytryptamine receptor) were simultaneously highly expressed in the hermaphroditic gonad to maintain the hermaphroditism of T. squamosa.

CONCLUSION

All these results from our study will facilitate better understanding of the molecular mechanisms underlying giant clam gonad development and gametogenesis, which can provided a base on obtaining excellent gametes during the seed production process for giant clams.

Transcriptomic effects of tributyltin (TBT) in zebrafish eleutheroembryos. A functional benchmark dose analysis

Exposure to the antifouling tributyltin (TBT) has been related to imposex in mollusks and to obesogenicity, adipogenesis and masculinization in fish. To understand the underlying molecular mechanisms, we evaluated dose-response effects of TBT (1.7-56 nM) in zebrafish eleutheroembryos transcriptome exposed from 2 to 5 days post-fertilization. RNA-sequencing analysis identified 3238 differentially expressed transcripts in eleutheroembryos exposed to TBT. Benchmark dose analyses (BMD) showed that the point of departure (PoD) for transcriptomic effects (9.28 nM) was similar to the metabolomic PoD (11.5 nM) and about one order of magnitude lower than the morphometric PoD (67.9 nM) or the median lethal concentration (LC₅₀: 93.6 nM). Functional analysis of BMD transcriptomic data identified steroid metabolism and cholesterol and vitamin D₃ biosynthesis as the most sensitive pathways to TBT (<50% PoD). Conversely, transcripts related to general stress and DNA damage became affected only at doses above the PoD. Therefore, our results indicate that transcriptomes can act as early molecular indicators of pollutant exposure, and illustrates their usefulness for the mechanistic identification of the initial toxic events. As the estimated molecular PoDs are close to environmental levels, we concluded that TBT may represent a substantial risk in some natural environments.

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.

Gonadal response of juvenile protogynous grouper (Epinephelus fuscoguttatus) to long-term recombinant follicle-stimulating hormone administration†

The role of follicle-stimulating hormone (FSH) in the gonadal development of protogynous hermaphroditic grouper (Epinephelus fuscoguttatus) was investigated. Recombinant giant grouper (E. lanceolatus) FSH (rggFSH) was produced in yeast. Its receptor-binding capacity and steroidogenic potency were confirmed in vitro. Weekly injections of rggFSH to juvenile tiger grouper for 8 weeks (100 μg/kg body weight, BW) resulted in significantly larger and more advanced oocytes (cortical alveolar stage vs primary growth stage in control). Sustained treatment with rggFSH (20 to 38 weeks at 200 μg/kg BW) resulted in significant reduction in gonad size, degeneration of oocytes, and proliferation of spermatogonial cells, indicative of female to male sex change. Gene expression analysis showed that, while initiating female to male sex change, the rggFSH significantly suppressed the steroidogenic genes cyp11b, cyp19a1a, and foxl2 which restrained the endogenous production of sex steroid hormones and thus prevented the differentiation of spermatogonial cells. Expression profile of sex markers dmrt1, amh, figla, and bmp15 suggests that the observed sex change was restricted at the initiation stage. Based on these results, we propose that the process of female to male sex change in the protogynous grouper is initiated by FSH, rather than sex steroids, and likely involves steroid-independent pathway. The cortical alveolar stage in oocyte development is the critical point after which FSH-induced sex change is possible in grouper.

Mummichog (Fundulus heteroclitus) are less sensitive to 17α-ethinylestradiol (EE2) than other common model teleosts: A comparative review of reproductive effects

The environmental estrogen 17α-ethinylestradiol (EE₂) will depress or completely inhibit egg production in many common model teleosts at low concentrations (≤0.5 ng/L; Runnalls et al., 2015). This inhibition is not seen in the estuarine killifish, or mummichog (Fundulus heteroclitus), even when exposed to 100 ng/L EE₂. This relative insensitivity to EE₂ exposure indicates species-specific mechanisms for compensating for exogenous estrogenic exposure. This review compares various reproductive responses elicited by EE₂ in mummichog to other common model teleosts, such as zebrafish (Danio rerio) and fathead minnow (Pimephales promelas), identifying key endpoints where mummichog differ from other studied fish. For example, EE₂ accumulates primarily in the liver/gall bladder of mummichog, which is different than zebrafish and fathead minnow in which accumulation is predominantly in the carcass. Despite causing species-specific differences in fecundity, EE₂ has been shown to consistently induce hepatic vitellogenin in males and cause feminization/sex reversal during gonadal differentiation in larval mummichog, similar to other species. In addition, while gonadal steroidogenesis and plasma steroid levels respond to exogenous EE_2, it is generally at higher concentrations than observed in other species. In mummichog, production of 17β-estradiol (E₂) by full grown ovarian follicles remains high; unlike other teleost models where E₂ synthesis decreases as 17α,20β-dihydroxy-4-prenen-3-on levels increase to induce oocyte maturation. New evidence in mummichog indicates some dissimilarity in gonadal steroidogenic gene expression responses compared to gene expression responses in zebrafish and fathead minnow exposed to EE₂. The role of ovarian physiology continues to warrant investigation regarding the tolerance of mummichog to exogenous EE₂ exposure. Here we present a comprehensive review, highlighting key biological differences in response to EE₂ exposure between mummichog and other commonly used model teleosts.

Immune genes, IL1β and Casp9, show sexual dimorphic methylation patterns in zebrafish gonads

There is crosstalk between the immune and reproductive systems in which sexual dimorphism is a common pattern in vertebrates. In recent years, epigenetics has emerged as a way to study the molecular mechanisms involved in gonadal development, those responsible for integrating environmental information that contribute to assigning a specific sexual phenotype (either an ovary or a testis). The knowledge of epigenetic mechanisms in certain molecular processes allows the development of epigenetic markers. In fish gonads, the existence of reproduction-immune system interactions is known, although the epigenetic mechanisms involved are far from clear. Here, we used the zebrafish (Danio rerio) as a model to study the DNA methylation patterns in gonads of two well-known innate immune genes: IL1β and Casp9. DNA methylation levels were studied by a candidate gene approach at single nucleotide resolution and gene expression analyses were also carried out. Results showed that there was clear sexual dimorphism in the DNA methylation levels of the two immune genes studied, being significantly higher in the testes when compared to the ovaries. In summary, and although further research is needed, this paper presents sexual dimorphic methylation patterns of two immune-related genes, thus sex-biased differences in methylation profiles should considered when analyzing immune responses in fish. Data showed here can help to develop epimarkers with forthcoming applications in livestock and fish farming production, for example, in immune fish diseases or sexual control programs as epigenetic molecular tools to predict environmental pressure in the gonads.

Genotoxic biomarkers and histological changes in marine medaka (Oryzias melastigma) exposed to 17α-ethynylestradiol and 17β-trenbolone

Endocrine-disrupting pollutants in marine environments have aroused great concern for their adverse effects on the reproduction of marine organisms. This study aimed to seek promising biomarkers for estrogenic/androgenic chemicals. First, two possible male-specific genes, SRY-box containing gene 9a2 (sox9a2) and gonadal soma-derived factor (gsdf), were cloned from marine medaka (Oryzias melastigma). Then the responses of sox9a2, gsdf, choriogenin (chgH and chgL), vitellogenin (vtg1 and vtg2), and cytochrome P450 aromatase (cyp19a and cyp19b) were investigated after exposure to 17α-ethynylestradiol (EE₂) and 17β-trenbolone (TB) at 2, 10, and 50 ng/L. The results showed that gsdf was specifically expressed in the testes and easily induced in the ovaries after TB exposure, indicating that gsdf was a potential biomarker of environmental androgens. ChgL was a useful biomarker of weak estrogen pollution for its high sensitivity to low levels of EE₂. In addition, both EE₂ and TB exposure damaged gonadal structures and inhibited gonadal development.

Identification of sex differences in zebrafish (Danio rerio) brains during early sexual differentiation and masculinization using 17α-methyltestoterone

Sexual behavior in teleost fish is highly plastic. It can be attributed to the relatively few sex differences found in adult brain transcriptomes. Environmental and hormonal factors can influence sex-specific behavior. Androgen treatment stimulates behavioral masculinization. Sex dimorphic gene expression in developing teleost brains and the molecular basis for androgen-induced behavioral masculinization are poorly understood. In this study, juvenile zebrafish (Danio rerio) were treated with 100 ng/L of 17 alpha-methyltestosterone (MT) during sexual development from 20 days post fertilization to 40 days and 60 days post fertilization. We compared brain gene expression patterns in MT-treated zebrafish with control males and females using RNA-Seq to shed light on the dynamic changes in brain gene expression during sexual development and how androgens affect brain gene expression leading to behavior masculinization. We found modest differences in gene expression between juvenile male and female zebrafish brains. Brain aromatase (cyp19a1b), prostaglandin 3a synthase (ptges3a), and prostaglandin reductase 1 (ptgr1) were among the genes with sexually dimorphic expression patterns. MT treatment significantly altered gene expression relative to both male and female brains. Fewer differences were found among MT-treated brains and male brains compared to female brains, particularly at 60 dpf. MT treatment upregulated the expression of hydroxysteroid 11-beta dehydrogenase 2 (hsd11b2), deiodinase, iodothyronine, type II (dio2), and gonadotrophin releasing hormones (GnRH) 2 and 3 (gnrh2 and gnrh3) suggesting local synthesis of 11-ketotestosterone, triiodothyronine, and GnRHs in zebrafish brains which are influenced by androgens. Androgen, estrogen, prostaglandin, thyroid hormone, and GnRH signaling pathways likely interact to modulate teleost sexual behavior.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

Morphometric signatures of exposure to endocrine disrupting chemicals in zebrafish eleutheroembryos

Understanding the mode of action of the different pollutants in human and wildlife health is a key step in environmental risk assessment. The aim of this study was to determine signatures that could link morphological phenotypes to the toxicity mechanisms of four Endocrine Disrupting Chemicals (EDCs): bisphenol A (BPA), perfluorooctanesulfonate potassium salt (PFOS), tributyltin chloride (TBT), and 17-ß-estradiol (E2). Zebrafish (Danio rerio) eleutheroembryos were exposed from 2 to 5 dpf to a wide range of BPA, PFOS, TBT and E2 concentrations. At the end of the exposures several morphometric features were assessed. Common and non-specific effects on larvae pigmentation or swim bladder area were observed after exposures to all compounds. BPA specifically induced yolk sac malabsorption syndrome and altered craniofacial parameters, whereas PFOS had specific effects on the notochord formation presenting higher rates of scoliosis and kyphosis. The main effect of E2 was an increase in the body length of the exposed eleutheroembryos. In the case of TBT, main alterations on the morphological traits were related to developmental delays. When integrating all morphometrical parameters, BPA showed the highest rates of malformations in terms of equilethality, followed by PFOS and, distantly, by TBT and E2. In the case of BPA and PFOS, we were able to relate our results with effects on the transcriptome and metabolome, previously reported. We propose that methodized morphometric analyses in zebrafish embryo model can be used as an inexpensive and easy screening tool to predict modes of action of a wide-range number of contaminants.

Zebrafish preserve global germline DNA methylation while sex-linked rDNA is amplified and demethylated during feminisation

The germline is the only cellular lineage capable of transferring genetic information from one generation to the next. Intergenerational transmission of epigenetic memory through the germline, in the form of DNA methylation, has been proposed; however, in mammals this is largely prevented by extensive epigenetic erasure during germline definition. Here we report that, unlike mammals, the continuously-defined ‘preformed’ germline of zebrafish does not undergo genome-wide erasure of DNA methylation during development. Our analysis also uncovers oocyte-specific germline amplification and demethylation of an 11.5-kb repeat region encoding 45S ribosomal RNA (fem-rDNA). The peak of fem-rDNA amplification coincides with the initial expansion of stage IB oocytes, the poly-nucleolar cell type responsible for zebrafish feminisation. Given that fem-rDNA overlaps with the only zebrafish locus identified thus far as sex-linked, we hypothesise fem-rDNA expansion could be intrinsic to sex determination in this species.

Germline cells transfer genetic information to offspring, and in zebrafish, drive sex determination. Here the authors report that, unlike mammals, the germline of zebrafish does not undergo genome-wide DNA methylation erasure, while amplifying and demethylating sex-linked rDNA during feminisation.

Endocrine disrupting effects of tebuconazole on different life stages of zebrafish (Danio rerio)

Tebuconazole is a widely used fungicide that has been detected in water ecosystems, of which the concentrations may affect the endocrine function of aquatic organisms. At present study, tissue-specific bioaccumulation of tebuconazole was found in ovary of adult zebrafish, indicating a potential risk of endocrine disruption. In order to evaluate the potential endocrine disrupting effects, three life stages (2 hpf (hours post-fertilization) -60 dpf (days post-fertilization), Stage I; 60-120 dpf, Stage II; 180-208 dpf, Stage III) of zebrafish (Danio rerio) were chronically exposed to tebuconazole at the concentrations ranging from 0.05 mg/L to 1.84 mg/L. Result showed that exposed to tebuconazole could lead to a male-biased sex differentiation in juvenile zebrafish and significant decrease of the percentage of germ cells in sexually-mature zebrafish. Egg production was significantly inhibited by 57.8% and 19.2% after Stage II- and Stage III-exposures, respectively. The contents of 17β-estradiol in gonad decreased by 63.5% when exposed to 0.20 mg/L tebuconazole at Stage II and by 49.5% after exposed to 0.18 mg/L tebuconazole at Stage III, respectively. For all stages exposure, reductions in 17β-estradiol/testosterone ratio were observed, indicating an imbalance in steroids synthesis. Additionally, tebuconazole reduced the expression of cyp19a, which was consistent with the decrease of E2 level. In overall, the present findings indicated that, playing as an anti-estrogen-like chemical, tebuconazole inhibited the expression of Cyp19, thereby impairing steroid hormones biosynthesis, leading to a diminished fecundity of zebrafish.

Molecular characterization and expression profiles provide new insights into GATA5 functions in tongue sole (Cynoglossus semilaevis)

GATA5 is a member of the GATA transcription factor family, which serves essential roles in varieties of cellular functions and biological processes. In this study, we have accomplished the molecular cloning, bioinformatic analysis and preliminary function study of C. semilaevis GATA5. The full-length cDNA nucleotide sequence is 1955 bp, with a coding sequence of 1167 bp, which encodes a polypeptide of 388 amino acids. Homology, phylogenetic, gene structure and synteny analysis showed that C. semilaevis GATA5 was highly conserved among vertebrates. Tissue distribution pattern exhibited that C. semilaevis GATA5 was significantly expressed in heart, intestine, liver, kidney and gonad, with a sexual dimorphic feature observed in testis and ovary. Embryonic development expression profiles showed that C. semilaevis GATA5 transcripts increased at the blastula stage, and peaked at the heat-beating period. Strong signals were detected at spermatids of male testis and stage III oocytes of female ovary by ISH. The expression of C. semilaevis GATA5 was regulated by 17α-MT and E2 after hormone stimulation to the ovary. Together, all the results pointed out that GATA5 might play a vital role during gonadal maturation and the reproductive cycle of C. semilaevis. This study lays the foundation for further researches on the sex control breeding in tongue sole.

Genetic mechanism underlying sexual plasticity and its association with colour patterning in zebrafish (Danio rerio)

BACKGROUND

Elevated water temperature, as is expected through climate change, leads to masculinization in fish species with sexual plasticity, resulting in changes in population dynamics. These changes are one important ecological consequence, contributing to the risk of extinction in small and inbred fish populations under natural conditions, due to male-biased sex ratio. Here we investigated the effect of elevated water temperature during embryogenesis on sex ratio and sex-biased gene expression profiles between two different tissues, namely gonad and caudal fin of adult zebrafish males and females, to gain new insights into the molecular mechanisms underlying sex determination (SD) and colour patterning related to sexual attractiveness.

RESULTS

Our study demonstrated sex ratio imbalances with 25.5% more males under high-temperature condition, resulting from gonadal masculinization. The result of transcriptome analysis showed a significantly upregulated expression of male SD genes (e.g. dmrt1, amh, cyp11c1 and sept8b) and downregulation of female SD genes (e.g. zp2.1, vtg1, cyp19a1a and bmp15) in male gonads compared to female gonads. Contrary to expectations, we found highly differential expression of colour pattern (CP) genes in the gonads, suggesting the 'neofunctionalisation' of those genes in the zebrafish reproduction system. However, in the caudal fin, no differential expression of CP genes was identified, suggesting the observed differences in colouration between males and females in adult fish may be due to post-transcriptional regulation of key enzymes involved in pigment synthesis and distribution.

CONCLUSIONS

Our study demonstrates male-biased sex ratio under high temperature condition and support a polygenic SD (PSD) system in laboratory zebrafish. We identify a subset of pathways (tight junction, gap junction and apoptosis), enriched for SD and CP genes, which appear to be co-regulated in the same pathway, providing evidence for involvement of those genes in the regulation of phenotypic sexual dimorphism in zebrafish.

Expression and transcriptional regulation of gsdf in spotted scat (Scatophagus argus)

Gonadal soma-derived factor (Gsdf) is critical for testicular differentiation and early germ cell development in teleosts. The spotted scat (Scatophagus argus), with a stable XX-XY sex-determination system and the candidate sex determination gene dmrt1, provides a good model for understanding the mechanism of sex determination and differentiation in teleosts. In this study, we analyzed spotted scat gsdf tissue distribution and gene expression patterns in gonads, as well as further analysis of transcriptional regulation. Tissue distribution analysis showed that gsdf was only expressed in testis and ovary. Real-time PCR showed that both gsdf and dmrt1 were expressed significantly higher in testes at different phases (phase III, IV and V) compared to ovaries at phase II, III and IV, while gsdf was expressed significantly higher in phase II ovaries than those of phase III and IV. Western blot analysis also showed that Gsdf was more highly expressed in the testis than ovary. Immunohistochemistry analysis showed that Gsdf was expressed in Sertoli cells surrounding spermatogonia in the testis, while it was expressed in the somatic cells surrounding the oogonia of the ovary. Approximately 2.7 kb of the 5' upstream region of gsdf was cloned from the spotted scat genomic DNA and in silico promoter analysis revealed the putative transcription factor binding sites of Dmrt1 and Sf1. The luciferase reporter assay, using the human embryonic kidney cells, demonstrated that Dmrt1 activated gsdf expression in a dose-dependent manner in the presence of Sf1 in spotted scat. These results suggest that Gsdf could play a role in regulating the development of spermatogonia and oogonia, and also participate in male sex differentiation by acting as a downstream gene of Dmrt1 in spotted scat.

Testicular transcriptome alterations in zebrafish (Danio rerio) exposure to 17β-estradiol

The hormone 17β-estradiol (E₂) can be found in rivers, effluents, and even drinking water. Researches have demonstrated that E₂ affects various metabolic pathways through gene activation and may cause reproductive toxicity in fish. Therefore, the aim of this study was to evaluate E₂-induced toxicity via testicular transcriptome of zebrafish (Danio rerio) exposed to different concentrations (10 ng L^-1, and 100 ng L^-1) of E₂. A total of >600 significant differentially expressed genes (DEGs) were enriched among the three treatments. Short time-series expression miner analysis revealed five KEGG pathways including drug metabolism, other enzymes, calcium signaling pathway, ECM-receptor interaction, gap junction, and cell adhesion molecules. Twenty genes were selected to verify the accuracy of RNA-Seq. Other reported genes related to sex differentiation, development, energy metabolism, and other processes were found. One set of genes significantly increased/decreased/fluctuated over time, especially 12 h after E₂ exposure. Genes associated with ovaries (zp3c), and development (bmp15, gdf9, and sycp2l) were significantly upregulated with increasing E₂ concentration. E₂ and testosterone was significantly decreased by 10 (except for T) and 100 ng L^-1 E₂ exposure at 12 h. The current study demonstrated that sex differentiation, development, energy metabolism, immunity, and ribosome biogenesis in male zebrafish were all significantly affected by 17β-estradiol exposure through transcriptional alterations.

Natural sex change in fish

Sexual fate can no longer be considered an irreversible deterministic process that once established during early embryonic development, plays out unchanged across an organism's life. Rather, it appears to be a dynamic process, with sexual phenotype determined through an ongoing battle for supremacy between antagonistic male and female developmental pathways. That sexual fate is not final and is actively regulated via the suppression or activation of opposing genetic networks creates the potential for flexibility in sexual phenotype in adulthood. Such flexibility is seen in many fish, where sex change is a usual and adaptive part of the life cycle. Many fish are sequential hermaphrodites, beginning life as one sex and changing sometime later to the other. Sequential hermaphrodites include species capable of female-to-male (protogynous), male-to-female (protandrous), or bidirectional (serial) sex change. These natural forms of sex change involve coordinated transformations across multiple biological systems, including behavioral, anatomical, neuroendocrine and molecular axes. Here we review the biological processes underlying this amazing transformation, focusing particularly on the molecular aspects, where new genomic technologies are beginning to help us understand how sex change is initiated and regulated at the molecular level.

Phenotypic plasticity induced using high ambient temperature during embryogenesis in domesticated zebrafish, Danio rerio

Ambient temperature during early stages of life has a substantial effect on physiological processes, eliciting phenotypic plasticity during zebrafish developmental stages. Zebrafish are known to possess a noteworthy ability to modify their phenotype in dependence of environmental factors. However, there is a poor understanding of the effects of temperature during embryogenesis, which influences the biological functions such as survival ability and masculinization in later developmental stages. Since the middle embryonic phase (pharyngula period) is genetically the most conserved stage in embryogenesis, it is very susceptible to embryonic lethality in developmental processes of vertebrates. Here, we tested the effect of transient perturbations (heat shock) during early development (5–24 hr post‐fertilization; hpf) at 35°C compared to control group at 28°C, on survival ability of zebrafish to study the embryonic and post‐embryonic mortality. We studied the variation of heat‐induced masculinization among and across the families in response to high temperature. Furthermore, morphometric traits of adult zebrafish at different developmental time points were measured in order to estimate the temperature × sex interaction effect. We found the highest embryonic mortality around the gastrula and segmentation periods in both experimental groups, with significantly lower survival ability in the temperature‐treated group (73.30% ± 0.58% vs. 70.19% ± 0.57%, respectively). A higher hatching success was observed in the control group (71.08% ± 0.61%) compared to the heat‐induced group (67.95% ± 0.60%). A distinct reduction in survival ability was also observed in both experimental groups during the first two weeks after hatching, followed by a reduced level of changes thereafter. We found sex ratio imbalances across all families, with 25.2% more males under temperature treatment. Our study on growth performance has shown a positive effect of increased temperature on growth plasticity, with a greater impact on female fish in response to high ambient temperature.

Roles of Two Sox9 Genes during Gonadal Development in Japanese Flounder: Sex Differentiation, Spermatogenesis and Gonadal Function Maintenance

The transcription factor sox9 has been implicated in cartilage formation and testis determination in mammals. Here, two duplicates of sox9 were found in Japanese flounder (Paralichthys olivaceus) named Posox9a and Posox9b, respectively. Phylogenetic and gene structure analyses revealed that Posox9a and Posox9b were homologous to that of teleosts and tetrapods. Quantitative real-time polymerase chain reaction (qRT-PCR) showed that both Posox9a and Posox9b expressed higher in testis than in ovary of adult tissues. The in situ hybridization (ISH) analysis of gonads showed that Posox9a and Posox9b mRNA were both detected in oocytes, Sertoli cells and spermatocytes. During sex differentiation, the expression of Posox9a exhibited obvious sexual dimorphic expression from 60 days after hatch (dah) with higher expression in male preferred individuals than female preferred individuals and increased gradually from 30 to 100 dah. A similar pattern was detected in Posox9b expression. After injection of androgen (17α-methyltestosterone) of different concentrations, the expression level of Posox9b increased significantly, whereas Posox9a did not change obviously. These results indicated that the two sox9 genes of Japanese flounder had converse functions in sex differentiation, whereas their differences in 17α-methyltestosterone administration were obvious and worthwhile for exploring evolutionary and adaptive significance. This study provided a foundation for further exploration of the roles of Posox9 genes during the sex determination and differentiation, spermatogenesis and gonadal function maintenance of Japanese flounder.