Function of Foxl2 and Dmrt1 proteins during gonadal differentiation in the olive flounder Paralichthys olivaceus

Transcriptome analysis of four types of gonadal tissues in largemouth bass (Micropterus salmoides) to reveal its sex-related genes

The sex determination system of largemouth bass (Micropterus salmoides, LMB) is XX/XY; however, the underlying molecular mechanisms involved in early sex differentiation, gonadal development, and exogenous hormone-induced sex reversal remain unknown. In this study, LMB at 15 days post-hatching (dph) were fed diets containing 20 mg/kg of 17α-methyltestosterone (17α-MT) or 30 mg/kg of 17β-estradiol (17β-E2) for 60 days, respectively. Serum steroid levels, histological observations of the gonads, and identification of sex-specific markers were employed to screen the gonads of 60-day-old normal female fish (XX-F), normal male fish (XY-M), 17β-E2 induced pseudo-female fish (XY-F), and 17α-MT-induced pseudo-male fish (XX-M) for transcriptome sequencing in order to uncover genes and pathway involved in the process of sexual reversal. The results from histology and serum sex steroid hormone analysis showed that both 17α-MT and 17β-E2 were capable of inducing sex reversal of LMB at 15 dph. Transcriptome results revealed a total of 2,753 genes exhibiting differential expression, and the expression pattern of these genes in the gonads of XX-M or XY-F resembled that of normal females or males. The male sex-biased genes that are upregulated in XX-M and downregulated in XY-F are referred to as key genes for male reversal, while the female sex-biased genes that are upregulated in XY-F and downregulated in XX-M are referred to as key genes for female reversal. Finally, 12 differentially expressed genes (DEGs) related to male sex reversal were screened, including star2, cyp17a, cyp11b1, dmrt1, amh, sox9a, katnal1, spata4, spata6l, spata7, spata18 and foxl3. 2 DEGs (foxl2a and cyp19a1b) were found to be associated with female sex reversal. The changes in these genes collectively influence the direction of sex differentiation of LMB. Among them, star2, dmrt1 and cyp19a1b with significantly altered expression levels may play potentially crucial role in the process of gender reversal. The expression patterns of 21 randomly selected genes were verified using qRT-PCR which confirmed the reliability and accuracy of the RNA-seq results. These findings not only enhance our understanding of the molecular basis underlying sex reversal but also provide crucial data support for future breeding research on unisexual LMB.

A time-course transcriptome analysis revealing the potential molecular mechanism of early gonadal differentiation in the Chinese giant salamander

The Chinese giant salamander (CGS) Andrias davidianus is the largest extant amphibian and has recently become an important species for aquaculture with high economic value. Meanwhile, its wild populations and diversity are in urgent need of protection. Exploring the mechanism of its early gonadal differentiation will contribute to the development of CGS aquaculture and the recovery of its wild population. In this study, transcriptomic and phenotypic research was conducted on the critical time points of early gonadal differentiation of CGS. The results indicate that around 210 days post-hatching (dph) is the critical window for female CGS's gonadal differentiation, while 270 dph is that of male CGS. Besides, the TRPM1 gene may be the crucial gene among many candidates determining the sex of CGS. More importantly, in our study, key genes involved in CGS's gonadal differentiation and development are identified and their potential pathways and regulatory models at early stage are outlined. This is an initial exploration of the molecular mechanisms of CGS's early gonadal differentiation at multiple time points, providing essential theoretical foundations for its captive breeding and offering unique insights into the conservation of genetic diversity in wild populations from the perspective of sex development.

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Genome-wide investigation of the DMRT gene family sheds new insight into the regulation of sex differentiation in spotted knifejaw (Oplegnathus punctatus) with fusion chromosomes (Y)

The role of the DMRT family in male sex determination and differentiation is significant, but its regulatory role in spotted knifejaw with Y fusion chromosomes remains unclear. Through genome-wide scanning, transcriptome analysis, qPCR, FISH, and RNA interference (RNAi), we investigated the DMRT family and the dmrt1-based sex regulation network. Seven DMRTs were identified (DMRT1/2 (2a,2b)/6, DMRT4/5, DMRT3), and dmrt gene dispersion among chromosomes is possibly driven by three whole-genome duplications. Transcriptome analysis enriched genes were associated with sex regulation and constructed a network associated with dmrt1. qPCR and FISH results showed the expression dimorphism of sex-related genes in dmrt-related regulatory networks. RNAi experiments indicated a distinct sex regulation mode in spotted knifejaw. Dmrt1 knockdown upregulated male-related genes (sox9a, sox9b, dmrt1, amh, amhr2) and hsd11b2 expression, which is critical for androgen synthesis. Amhr2 is located on the heterozygous chromosome (Y) and is specifically localized in primary spermatocytes, and is extremely upregulated after dmrt1 knockdown which suggested besides the important role of dmrt1 in male differentiation, the amhr2 along with amhr2/amh system, also play important regulatory roles in maintaining high expression of the hsd11b2 and male differentiation. This study aims to further investigate sex regulatory mechanisms in species with fusion chromosomes.

Molecular and functional characterization of ribosome protein S24 in ovarian development of Macrobrachium nipponense

Ribosomal proteins (RPs) have mang extraribosomal functions including regulation of ovarian development in some organisms. In order to solve the problem of rapid ovarian maturation in Macrobrachium nipponense aquaculture, this study identified a RPS24 (MnRPS24) gene from M. nipponense, which encodes a protein of ββαβαααα folding structure type. MnRPS24 exhibited the greatest expressions in the female adult stage among the six growth stages, in the ovary among the nine tissues, and in the stage I ovary among the six ovarian development stages. The MnRPS24 protein located in the cytoplasm of oogonia, previtellogenic and early-vitellogenic oocytes, and the follicular cells surrounding the oocytes. The expression of the vitellogenin (MnVg), vitellogenin receptor (MnVgr), cell cycle protein B (MnCyclin B) and cell division cyclin 2 (MnCdc2) genes were increased by recombinant MnRPS24 protein incubation. Conversely, the expression of the Wee1 kinase (MnWee1) gene was decreased. MnRPS24 gene silencing downregulated the expression for MnVg, MnVgr, MnCyclin B and MnCdc2 and upregulated the expression for MnWee1. Furthermore, MnRPS24 gene silencing delayed the vitellogenesis of oocytes, halting the progression of ovarian development. The findings of this research demonstrate that MnRPS24 could potentially function as a stimulator in promoting the development of ovaries in M. nipponense.

Cloning, Identification, and Functional Analysis of the Foxl2 Gene in Procambarus clarkii

Procambarus clarkii is the most widely distributed freshwater shrimp in China, with important economic value and great potential for development. The forkheadboxL2 (Foxl2) gene has been found to be involved in the reproductive development of many crustaceans. To understand the role of the Foxl2 gene in the gonad development of P. clarkii, we designed CDS-specific primers for the P. clarkii Foxl2 (PcFoxl2) gene and cloned its CDS sequence using RT-PCR. The nucleotide and protein sequence information was then analyzed through bioinformatics analysis. The expression and subcellular localization of PcFoxl2 in various tissues were detected using qRT-PCR and in situ hybridization. The effects of PcFoxl2 knockdown on gonad development were investigated using RNA interference. The results showed that the CDS length of the PcFoxl2 gene was 1614 bp and encoded 537 amino acids. Protein sequence comparison and phylogenetic analysis showed that PcFoxl2 was the closest relative to Crayfish. qRT-PCR analysis indicated that the expression level of PcFoxl2 in the testis was significantly higher (>40 fold) than that in the ovary (p < 0.01). The in situ hybridization results showed that PcFoxl2 was expressed in both the cytoplasm and the nucleus of egg cells, and that the expression was strongest in egg cells at the early stage of yolk synthesis, while weak in the secondary oocytes. The positive signal was strongest in the spermatocyte nucleolus, while only a trace signal was observed in the cytoplasm. After interfering with the PcFoxl2 gene using dsRNA, the expression of PcFoxl2 in the RNA interference group was significantly lower than that in the control group, and this interference effect lasted for one week. Moreover, the gonad index of the experimental group was significantly lower than that of the control group (p < 0.05) after 10 days of P. clarkii cultivation following PcFoxl2 knockdown. The expression levels of the nanos and S3a genes, which are related to gonad development, decreased significantly after PcFoxl2 gene interference. The results suggest that the Foxl2 gene is involved in the growth and development of gonads, particularly in the development of testis, and is related to the early development of oocytes. This study provides a theoretical basis for the artificial breeding of P. clarkii.

Expression characteristics of Hsd3b7 in the gonads of Paralichthys olivaceus

Steroidogenesis is an important biological process for gonadal differentiation and development. In mammals, 3β-hydroxysteroid dehydrogenase 7 (HSD3B7) could convert 3β-hydroxy of 7α-hydroxycholesterol into a ketone and form 7α-hydroxy-4-cholesten-3-one, which may affect steroidogenesis. However, in fish, the study of Hsd3b7 is still lacking. In this study, Hsd3b7 was identified in the olive flounder Paralichthys olivaceus, an important mariculture fish. According to bioinformatics analysis, Hsd3b7 belongs to a Rossmann-fold NAD(P)(+)-binding protein and can interact in a predictable manner with Hsd17b2, -3, and - 4, which play a role in steroidogenesis. In the adult flounder, Hsd3b7 was expressed in various tissues, at particularly high level in male muscle. The expression levels of Hsd3b7 at gonadal development stages I-V initially increased and then decreased, with an inflection point in the ovary at stage III and in the testis at stage IV. At stage III, the expression level of Hsd3b7 was significantly higher in the ovary than in the testis (P < 0.01). The results of in situ hybridization (ISH) revealed that it was mainly expressed in oocytes of phases I-IV or around oocytes of phases IV-V in the ovaries and around spermatid lobules at stages IV-V in the testes. Three regulatory sites of SRY-box transcription factor 9 (Sox9), a transcription factor involved in steroidogenesis and gonadal differentiation, were predicted in the promoter of Hsd3b7. After intraperitoneal injection with the recombination flounder Sox9a, the expression of Hsd3b7 was significantly up-regulated (P < 0.01). During the flounder gonadal differentiation, 17β-estradiol (E2, 5 μg/g feed) and 17α-methyltestosterone (T, 5 μg/g feed) were used to obtain the phenotypic female or male flounder, and the results showed that in the E2 group, Hsd3b7 expression was highest at 2 cm TL, the primordial gonad stage, which was significantly higher than that at 12 cm TL (P < 0.05). In the T group, Hsd3b7 expression level was also highest at 2 cm TL and significantly higher than at 10 and 12 cm TL (P < 0.05). Moreover, Hsd3b7 was detected to be localized mainly around oogonia and spermatogonia during the differentiated gonads with ISH. These findings first introduce the expression characteristics of Hsd3b7 and the effect of Sox9a on its expression, which contribute to our understanding of the function of Hsd3b7 in fish gonads.

Rxrs and their partner receptor genes inducing masculinization plausibly mediated by endocrine disruption in Paralichthys olivaceus

Retinoid X receptors (RXRs) can form homo- or heterodimers with orphan receptors involved in multiple intertwined signaling pathways. However, there is limited study on the formation of sex phenotypes and the regulation of steroidogenesis by RXRs in fish. Here, in Paralichthys olivaceus, we first indicated that PPARγ::RXRα was predictably a transcription factor for steroidogenesis genes, and Foxl2 and Dmrt1 were also transcription factors for rxrs and their partner receptor genes. When the flounder fry were exposed to LG100268 (LG, RXRs agonist, 50 mg/kg diet), the percentage of males increased from 50% to 71.4%. Before histological differentiation of the flounder ovary (3 cm TL) and testis (6 cm TL), the transcripts of rar β and rar γ (P < 0.05) were activated, and the steroidogenesis gene Hsd3b1 was down-regulated (P < 0.05). The ratios of testosterone (T)/17β-estradiol (E2) were all greatly increased (P < 0.05), and the ratio of 11-ketotestosterone (11-KT)/E2 was elevated at 3 cm TL. Moreover, LG was used to treat the cultured gonads in vitro (10 μM) and the fish with intraperitoneal injection in vivo (12 mg/kg body weight), respectively. LG was able to up-regulate rxr γ, rar γ, and ppar δ, and Hsd3b1 was significantly up-regulated (P < 0.05). The ratios of 11-KT/E2 in the culture medium and in the ovaries of the fish were decreased. Furthermore, the recombinant flounder Foxl2 protein was able to significantly down-regulate ppar γ (P < 0.05) and tr β (P < 0.01) in the ovaries in vitro, and the result of the Dmrt1 in the testes was opposite to that of the Foxl2, probably indicating a feedback loop between RXRs' partner receptors and Foxl2/Dmrt1. These findings introduce for the first time the mode of action of RXRs on the flounder steroidogenesis and provide important data to learn the potential function of RXRs in fish sex differentiation and the potential role of RXRs in aquatic animals in the presence of water pollutants.

Dynamic transcriptome analysis reveals the gene network of gonadal development from the early history life stages in dwarf surfclam Mulinia lateralis

BACKGROUND

Gonadal development is driven by a complex genetic cascade in vertebrates. However, related information remains limited in molluscs owing to the long generation time and the difficulty in maintaining whole life cycle in the lab. The dwarf surfclam Mulinia lateralis is considered an ideal bivalve model due to the short generation time and ease to breed in the lab.

RESULTS

To gain a comprehensive understanding of gonadal development in M. lateralis, we conducted a combined morphological and molecular analysis on the gonads of 30 to 60 dpf. Morphological analysis showed that gonad formation and sex differentiation occur at 35 and 40-45 dpf, respectively; then the gonads go through gametogenic cycle. Gene co-expression network analysis on 40 transcriptomes of 35-60 dpf gonads identifies seven gonadal development-related modules, including two gonad-forming modules (M6, M7), three sex-specific modules (M14, M12, M11), and two sexually shared modules (M15, M13). The modules participate in different biological processes, such as cell communication, glycan biosynthesis, cell cycle, and ribosome biogenesis. Several hub transcription factors including SOX2, FOXZ, HSFY, FOXL2 and HES1 are identified. The expression of top hub genes from sex-specific modules suggests molecular sex differentiation (35 dpf) occurs earlier than morphological sex differentiation (40-45 dpf).

CONCLUSION

This study provides a deep insight into the molecular basis of gonad formation, sex differentiation and gametogenesis in M. lateralis, which will contribute to a comprehensive understanding of the reproductive regulation network in molluscs.