Huan Li
Qihui Zhu
Ruiyi Chen2,3
Mingtao Liu
Dongdong Xu- 1School of Fisheries, Zhejiang Ocean University, Zhoushan, China
- 2Key Lab of Mariculture and Enhancement of Zhejiang Province, Zhejiang Marine Fisheries Research Institute, Zhoushan, China
- 3Ocean and Fisheries Research Institute, Zhejiang Ocean University, Zhoushan, China
The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.
Introduction
Teleost fish display an amazing diversity of sex-determination systems, such as XX/XY, ZZ/ZW, XX/XO, and ZZ/ZO. Systems that are more complicated can involve multiple sex chromosomes and multiple gene loci (influence of autosomal loci on sex determination and polyfactorial sex determination). The multiple ♀X1X1X2X2/♂X1X2Y sex chromosome system is characterized by a neo-Y sex chromosome, with males containing an exceptionally large heterotypic chromosome and possessing one less chromosome than females. To date, over 50 cases with multiple sex chromosomes have been reported in teleost phylogeny (Sember et al., 2015; Bitencourt et al., 2017; Zhang et al., 2018a; Krysanov and Demidova, 2018; Cai et al., 2019; Xu et al., 2019). However, little is known about the molecular sex determination mechanism of multiple sex chromosome system in fish. The rock bream (Oplegnathus fasciatus), a member of the Oplegnathidae family of the Centrarchiformes, is a commercially important rocky reef fish native to East Asia and has been reported as a typical fish with multiple X1X1X2X2/X1X2Y sex chromosome system (Xu et al., 2012; Xue et al., 2016; Xu et al., 2019). Owing to their high aquaculture value, unique sex-determining system, and susceptibility to widespread biotic diseases, rock bream fish have drawn increasing research interest. For example, studies have been conducted on the characterization of the differences between X and Y chromosomes by molecular markers and cytological genetic methods and preliminary clarification of the origin and differentiation of a multiple X1X2Y sex chromosome system (Xu et al., 2013a; Xu et al., 2013b; Xu et al., 2019). Moreover, the whole genomes of both male and female rock breams were sequenced, assembled, and analyzed, and their results suggested that the neo-Y chromosome probably originated from the centric fusion of acrocentric chromosomes (X1 and X2 chromosomes) (Shin et al., 2018; Xiao et al., 2019; Xiao et al., 2020). These studies have contributed to the understanding of sex determination in species with a unique multiple sex chromosome system. Nevertheless, to date, little is known about the molecular pathways and regulatory mechanisms of sex determination in fish with multiple sex chromosome system.
Genetic sexual determination has always been a hot topic in the field of aquatic research and involves a series of complex biochemical interactions that lead to gonadal sexual determination and differentiation. Genetic systems in vertebrates may operate by polygenic control mechanisms and dominant genes, along with autosomic controls and sex chromosomes. The genes involved in sexual determination may be distributed in the genome, localized in individual chromosomes, or restricted to a single genetic locus (Robert and Nagahama, 2002). In monogenic systems, sex is determined by a specific gene on the sex chromosome, which has strong specificity and functionality. Only a handful of sex-determining genes have been identified in fish, including the Y chromosome-linked anti-Müllerian hormone (Amhy) gene of Odontesthes hatcheri (Ricardo et al., 2012) and Nile tilapia (Oreochromis niloticus) (Sissao et al., 2019), Sdy of rainbow trout (Oncorhynchus mykiss) (Yano et al., 2012), and anti-Müllerian hormone type II receptor (AmhrII) of Fugu rubripes (Takifugu rubripes) (Kamiya et al., 2012). A number of genes were identified to be closely related to sexual determination and differentiation in the complex sexual determination system in fish. For example, the anti-Müllerian hormone (Amh), which is expressed in Sertoli cells, has been reported to inhibit the expression of ovary specific genes in Nile tilapia (O. niloticus) (Ijiri et al., 2008). Although studies on Nile tilapia (O. niloticus) have shown that Amh plays a key role in male sex determination, the regulatory mechanism of this gene remains unclear. A number of genes related to sexual determination and differentiation have been discovered in many fishes, such as Nanos2 (Tsuda et al., 2003), Sox9 (Chen et al., 2015), and Gsdf (Jiang et al., 2016). These studies have provided a comprehensive understanding of sex determination in fish, whereas information on the mechanism of sexual determination in the multiple sex chromosome system is still scarce.
RNA-Seq has been proven to be a powerful strategy for quantifying gene expression levels, even when the expression levels are low. Transcriptomic studies on gonads have been conducted in a relatively large number of fish species, including Nile tilapia (O. niloticus) (Tao et al., 2013; Tao et al., 2018), spotted knifejaw (Oplegnathus punctatus) (Li et al., 2020), fugu rubripes (T. rubripes) (Yan et al., 2018), and spotted scat (Scatophagus argus) (He et al., 2019). For example, in the channel catfish (Ictalurus punctatus), male-preferential genes, such as Gsdf and cxcl12, were identified by transcriptomic comparison of the testes and ovaries (Zeng et al., 2016). In the spotted knifejaw (O. punctatus), a set of unigenes related to gonadal development and gametogenesis was filtered by transcriptomic analysis and genes, such as bmp15, Nanos3, Sox9, and Amh, were found to have roles in the regulation of gonadal physiology and germ line cell maintenance in the ovaries and testes (Du et al., 2017).
In this study, transcriptional analyses were carried out on rock bream to investigate the male and female gonads at different developmental stages to filter potential genes and pathways associated with sex determination. Next, several genes were screened and further analyzed by quantitative PCR. The results of this study will increase our understanding of the molecular mechanisms involved in sex determination and differentiation in fish with multiple sex chromosomes.
Materials and Methods
Ethics Statement
The collection and handling of all animals used in this study were approved by the Animal Care and Use Committee of Zhejiang Ocean University and the Zhejiang Marine Fisheries Research Institute. All experimental procedures were performed in accordance with the prescribed guidelines.
Samples and Tissue Collection
The experimental fish were reared in the hatchery of the Zhejiang Marine Fisheries Research Institute, Xixuan Island, Zhoushan, China. The rock breams were randomly sampled at the age of 1 (245.50 ± 16.82 g), 1.5 (315.45 ± 18.5 g) and 2Y (year old) (397.68 ± 20.05 g). At each sampling time point, 10 fish were randomly chosen and dissected immediately. The brain samples were immediately stored in liquid nitrogen for RNA isolation and each gonad sample was divided into two aliquots: one was stored in liquid nitrogen and the other was fixed with Bouin’s solution for histological analysis. The brain and gonad tissues from the 1.5Y rock bream were also prepared for transcriptional analysis via RNA-seq. For the 1.5Y rock breams, in addition to the brain, gonad, heart, liver, kidney, gill, intestine, and muscle were collected and stored in liquid nitrogen. After sex detection by histological observation of the gonad tissue, three samples from each male and female from each sampling time were randomly chosen for RNA isolation and subsequent analyses.
Total RNA Extraction and Sequencing
Samples of the gonads, brain, and other tissues were stored in liquid nitrogen for gene expression analysis. Total RNA was extracted using a Takara RNA kit (RR047A, Dalian), and the quality was determined using 1.2% denatured agarose gel. The concentration and quality of RNA were detected by AB-Sorbance and agarose gel electrophoresis at 260 nm, respectively. With samples for transcriptional analysis (three samples of each gonad and brain tissues from both male and female fish), RNA integrity was assessed using an Agilent 2100 biological analyzer (Agilent technologies), and samples with RNA integrity numbers ≥7 were subsequently analyzed. Libraries were constructed and sequenced using the BGISEQ-500 platform.
Transcriptional Analysis and Functional Annotation
Raw data were processed using the SOAPnuke method. After removing the adaptor and low-quality sequences, the clean reads were subsequently aligned with Bowtie2 and mapped to the rock bream genome (NCBI accession number SRP220007) using HISAT. The FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) (Trapnell et al., 2010) and read count values of each gene were calculated using Bowtie2 (Langmead and Salzberg, 2012) and RSEM (Li and Dewey, 2011). Pearson correlation coefficients between two samples were calculated using the core function in R software. The DEGseq method was based on the Poisson distribution, and the present project performed differentially expressed gene (DEG) detection according to a previously described method (Anders and Huber, 2012). To improve the accuracy of DEGs, we defined the genes with more than two-fold difference and Q-value ≤ 0.001 (fold change ≥ 2 and adjusted p-value ≤ 0.001), which were screened as significant DEGs. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs were performed in R based on hypergeometric distribution.
Quantitative Polymerase Chain Reaction (qPCR) Validation
Quantitative PCR was performed on a StepOnePlus Real-Time PCR system (Thermo Fisher Scientific, United States) using SYBR Green Real-time PCR Mix (RR420A; Takara, Dalian), following the manufacturer’s protocol. An isolation step was conducted after each PCR run to verify amplification specificity. Three fish samples were randomly chosen as biological replicates in qPCR. For each sample, qPCR was performed in triplicate. Primer pairs for qPCR were designed using the Primer Premier V6.0. The details of these genes are listed in Supplementary Table S1. Target gene expression was analyzed using the 2−ΔΔCT method (Kenneth and Thomas, 2001). The copy number of a reference gene, β-actin, was used to normalize the expression values.
Histology Observation
The gonads were dissected and fixed in the Bouin’s solution, and the samples were embedded in paraffin using conventional techniques. Sections (5 μm) were cut, mounted on slides, and stained with hematoxylin and eosin. Images of the sections were obtained using a microscope and digital camera (Axio Imager A2; digital camera (Axiocam 506; Zeiss). The phenotypic sex of each fish was determined, and the development of testes and ovaries was observed.
Results
Transcriptome Data Analysis
A total of 12 samples were analyzed via RNA-Seq in this study, with three samples in each group, including ovary (F_sex group), testis (M_sex group), female brain (F_brain group), and male brain (M_brain group). As a result, an average output of 6.50 GB bases per sample was generated, with an average clean reads ratio of up to 89.93%. Detailed statistics of the transcriptome data are summarized in Supplementary Table S2. As the average genome mapping rate of samples was 90.96%, a total of 25,678 genes were obtained, of which 22,723 were known genes and 2,955 were predicted new genes. A total of 41,127 new transcripts were obtained, of which 29,536 were new alternative splicing subtypes of known protein-coding genes, 2,974 were new protein-coding genes, and the remaining 8,617 were long non-coding RNAs. The correlation coefficients between replicated samples were higher than 0.88 (Supplementary Figure S1), indicating the conserved similarity between the replicates, and the results also showed good correlation coefficients between all brain samples. The results of principal component analysis showed that the replicates of ovaries, testes, and brains were clustered closely in a region, further validating the results of the correlation coefficient analysis (Supplementary Figure S2). DEGs between males and females were identified (Supplementary Figure S3). The F_sex-vs-M_sex group contained 17,033 upregulated and 2024 downregulated DEGs, whereas the F_brain-vs-M_brain group contained only 1,325 upregulated and 523 downregulated DEGs.
In terms of the DEGs of F_brain-vs-M_brain group, GO classification results (Figure 1A) highlighted that GO terms of “biological regulation” (181 genes) and “cellular process” (203 genes) in the biological process category, “membrane” (269 genes) and “membrane part” (254 genes) in the cellular component category, and “catalytic activity” (204 genes) and “binding” (322 genes) in the molecular function category contained the highest number of DEGs. Further GO enrichment analysis showed that a total of 50 GO terms were significantly enriched (Q value < 0.05), including several GO terms relating to development or differentiation, such as “central nervous system development” and “spinal cord motor neuron differentiation” (Supplementary Figure S4). KEGG annotation results showed that for the DEGs between male and female brains, the terms of “signal transduction” (340 genes), and “signaling molecules and interaction” (240 genes), were classified with the highest number of genes (Figure 2A). KEGG enrichment analysis showed that 29 pathways were significantly enriched (Q value < 0.05), of which 10 belonged to the category of “organismal systems,” and the pathway of “neuroactive ligand-receptor interaction” (105 genes) contained the highest number of genes (Supplementary Figure S5).
FIGURE 1. GO classification based on the DEGs of the F_brain-vs-M_brain group (A) and the F_sex-vs-M_sex group (B). The horizontal axis shows the number of DEGs enriched in each GO term. GO, gene ontology; DEG, differentially expressed gene.
FIGURE 2. KEGG classification based on the DEGs of the F_brain-vs-M_brain group (A) and the F_sex-vs-M_sex group (B). The horizontal axis shows the number of DEGs enriched in each KEGG term. KEGG, Kyoto Encyclopedia of Genes and Genomes; DEG, differentially expressed gene.
For the DEGs of the F_sex-vs-M_sex group, which represented the molecular differences between testis and ovary in 1.5Y rock breams, the GO terms with the top number of DEGs were “cellular process” (2,250 genes), “biological regulation” (1,570 genes), “membrane” (2,632 genes), “membrane part” (2,484 genes), “binding” (3,619 genes), and “catalytic activity” (2,314 genes) (Figure 1B). Although these were the same terms as the result of the F_brain-vs-M_brain group, the GO enrichment results were quite different. A total of 30 GO terms were significantly enriched, including nine development-related terms, such as “nervous system development” and “multicellular organism development” (Supplementary Figure S6). Further KEGG annotation results showed that the terms of “signal transduction” (3,077 genes) and “signaling molecules and interaction” (240 genes) were classified with the highest number of genes (Figure 2B). In addition, 40 pathways were significantly enriched, and the pathways of “human papillomavirus infection” (652 genes) and “endocytosis” (529 genes) contained the highest number of genes (Supplementary Figure S7).
A Venn diagram representing the DEGs in the F_brain-vs-M_brain and the F_sex-vs-M_sex group (Supplementary Figure S8) illustrated 335 DEGs that were specific to the F_brain-vs-M_brain group, and 17,544 DEGs specific to the F_sex-vs-M_sex group, with 1,513 overlapping. Then, we set |log2Ratio| ≥ 2 (fold change ≥ 4) as a filtering threshold to identify significant DEGs. For the DEGs specific to the F_brain-vs-M_brain group, 103 significant DEGs were identified, and the GO terms “extracellular region” and “hormone activity” were significantly enriched, which included 11 and five genes, respectively. In terms of the DEGs specific to the F_sex-vs-M_sex group, 11,904 significant DEGs were identified, and 51 GO terms (Supplementary Figure S9) and 84 KEGG (Table S3) pathways were further significantly enriched, with the pathways of “PI3K-Akt signaling pathway” (388 genes), “MAPK signaling pathway” (351 genes), “Ras signaling pathway” (289 genes), and “Rap1 signaling pathway” (288 genes) in the top ranking of highest number of DEGs. For the overlapping DEGs, 284 significant DEGs were obtained, and 13 GO terms were further significantly enriched, which mostly were development-related, including “ventral spinal cord development,” “nervous system development,” and “forebrain development” (Supplementary Figure S10).
Validation by qPCR
To validate the reliability and accuracy of the RNA-Seq results, 18 DEGs were selected for qPCR analysis, including Amh, Amhr, Cyp19a, Dmrt1, Dmrtb1, Foxl2, Gsdf, Nanos2, Sox9, vasa, ZP4, Bouncer, CA4, FoxH1, H2A, RNF208, TOB, and somatostatin (Figure 3). The coefficient of determination (R2) between the Log2 (fold change) values of RNA-Seq and the qPCR results was 0.83, thus indicating good consistency.
FIGURE 3. Validation of RNA-Seq data by qPCR. Comparison of the relative log2 (Fold Change) values of 18 genes between RNA-Seq and qPCR results. qPCR, quantitative PCR.
Tissue Distribution Analysis by qPCR
According to the RNA-Seq results, we selected male-or female-biased expressed genes and analyzed their expression in different tissues by qPCR. The tissue distribution results (Figure 4) showed that for most of the male-biased genes, such as Dmrt1, Dmrtb1, Amh, Nanos2, and Sox9, their expression levels in the testis were significantly higher than in other tissues. However, we also detected high expression of Dmrtb1 in the liver, Sox9 in the liver and brain, and Amhr and Gsdf in the heart and liver. As for the female-biased genes, the expression levels of ZP4, CA4, RNF208, H2A, somatostatin, and TOB were significantly higher in the ovary than those in other tissues. In addition, Bouncer showed relatively high expression in the gills, and FoxH1 was highly expressed in the liver. Interestingly, although the expression levels of vasa peaked in the testis, it was also relatively highly expressed in the ovary, heart, liver, and both male and female brains (Figure 4A).
FIGURE 4. The expression patterns of selected genes in different tissues of the rock bream. (A) Heat map of the expression patterns for 18 selected genes; Red indicates relative high expression, blue indicates relative low expression, and color depth indicates relative expression levels; (B) Detailed expression patterns of selected genes (the histograms for the rest 12 genes are shown in Supplementary Figure S11).
Sex Dimorphic Gene Expression at Different Gonadal Development Stages
The expression patterns of 18 selected genes in the gonads of rock breams at different developmental stages were analyzed by qPCR (Figure 5). We first analyzed gonadal development in rock bream using histological staining. The ovary was in stage II at 1Y, with the main cell types in the ovary being oogonias and phase I and phase II oocytes (Figures 5Ca). The testes of the rock bream at this age had developed into stage III, a large number of sperm cells were observed, and the internal structure of the testis had developed and matured (Figures 5Cd). For the 1.5Y fish, the ovary structure had developed, and mainly contained phase II–IV oocytes and emptied follicles (Figures 5Cb). In addition, a small amount of sperm was observed in the testis, along with some spermatogonia and spermatocytes (Figures 5Ce). For the 2Y fish, the ovary was in stage VI–V ovary and mainly contained phase IV and V oocytes (Figures 5Cc). At this age, the seminiferous lobules were observed to be partially melted, and the whole testis was filled with sperm cells and mature sperm, and there were still some spermatocytes (Figures 5Cf).
FIGURE 5. The expression patterns of selected genes in the gonads and the histological observation of the gonads for the rock breams at different development stages (1, 1.5, and 2Y). (A) Heat map of the expression patterns for 18 selected genes; Red indicates relative high expression, blue indicates relative low expression, and color depth indicates relative expression levels; (B) Detailed expression patterns of selected genes (the histograms for the rest nine genes were shown in Supplementary Figure S12). (C) Histological observation of the gonads at different development stages. Scale bars = 50 μm. Y, year old; Oc, oocyte; Sp, sperm; Sc, spermatocyte; Sg, spermatogonia.
As shown in the heat map (Figure 5A), the expression patterns of the genes were mainly clustered into three branches: the expression levels of Amh, Cyp19a, Dmrt1, Dmrtb1, Nanos2, vasa, and Sox9 in the testis were significantly higher than those in the ovary, and ZP4, FoxH1, Bouncer, CA4, RNF208, H2A, somatostatin, and TOB were found to be highly expressed in the ovaries, whereas Amhr and Gsdf had relatively high expression in the gonads of both 2Y male and female fish. In particular, the expression of Dmrt1, Dmrtb1, Nanos2, and Sox9 increased with the age of the fish, but Amh showed the opposite pattern (Figure 5B). Interestingly, for most of the female biased expressed genes, such as ZP4, CA4, RNF208, and TOB, their expression in the ovaries reached peaks at 1.5Y fish instead of in 1Y or 2Y (Figure 5B).
Discussion
In this study, for a better understanding of sex determination and differentiation in the rock bream, we investigated the gene expression profiles in the gonads and brains of males and females by conducting a comparative transcriptome analysis and screened potential sex-related genes and pathways. We detected a low number of DEGs in the brains of males and females, which indicated fewer differences in the brain between sexes at the transcriptomic level. However, many DEGs related to nervous system development were significantly enriched. In particular, 103 significant DEGs were identified to be specific in the F_brain-vs-M_brain group, and these could be potential marker genes for exploring the differences between male and female brains. The number of DEGs between the ovary and testis was quite large (19,057 DEGs), and the identified significant DEGs specific to the F_sex-vs-M_sex group was also relatively large (11,904 DEGs). In addition, many signaling pathways were further enriched, including “PI3K-Akt signaling pathway,” “MAPK signaling pathway,” “Wnt signaling pathway,” and “Janus kinase-signal transducer and activator of transcription (Jak-STAT) signaling pathway.” Many of these factors have been reported to be related to gonadal development. For example, Jak-STAT signaling pathway is an evolutionarily conserved pathway that has been reported to play a pivotal role during development in both vertebrates and invertebrates, including the development of somatic and germ line cells (Luo and Dearolf, 2001). In the present study, the expression levels of Amh, Dmrt1, Sox9, and Nanos2 in testes were significantly higher than those in female ovaries and other tissues, and the expression levels of these genes in testes also varied significantly at different developmental stages, indicating that these genes might play important roles during testicular development in O. fasciatus.
Previous studies have reported that Amh, Dmrt1, vasa, Sox9, and Nanos2 genes are related to sex differentiation (Schmidt et al., 2004; Daniel et al., 2005; Rodríguez-Marí et al., 2005; Alam et al., 2008; Sun et al., 2017; Zhang et al., 2018b), reproductive system development (Tsuda et al., 2003; Kuiri-Hänninen et al., 2011), and growth of aquatic species (Li et al., 2012) and are closely related to the differentiation and development of testis (Adolfi et al., 2015; Cui et al., 2017), which is consistent with our results. In addition, the transcription factor Dmrtb1 was found to have a specific expression in testis. This gene is seldom reported in fish and currently, there are only a few reports of its occurrence in mice and humans, which showed that it might play a key role in the transition between mitosis and meiosis of mouse germ cells as well as in the process of male spermatogonial cells during meiosis (Hilbold et al., 2019). Amh has been studied in many vertebrates, including teleosts and mammals. Reports have shown that Amh is mainly expressed in Sertoli cells of the testis and granulosa cells of the ovary. It plays an important role in the process of sex reversal in medaka (Oryzias latpes) (Wang et al., 2018), maintaining the normal proliferation and differentiation of female and male germ cells in zebrafish (Danio rerio) (Zhang et al., 2020) and ensuring a balance between proliferation and differentiation. In this study, Amh was deteced to be significantly highly expressed in the testes of rock bream than other tissues including the ovaries, which was consistent with the results of previous studies. Although previous studies have shown that vasa is specifically expressed in the gonads in most teleost, where it plays an essential role in the formation of germ cells, several studies have also reported that vasa was expressed in tissues other than the gonads, such as the heart and brain of O. mykiss (Yoshizaki et al., 2000), the heart of tongue sole Cynoglossus semilaevis (Zhang et al., 2012). And such extragonadal expression was also detected in Mammalian (Zamboni and Upadhyay, 1983) and Amphibian (Ikenishi and Tanaka, 2000; Jia et al., 2009). Nevertheless, the functional significance of extragonadal expression of vasa remains to be determined. Dmrt1 is reported to be one of the most conserved genes involved in sex determination and differentiation. Studies on Nile tilapia showed that Dmrt1 was only expressed in mature testes, and it was also expressed in the gonads of XX tilapia after sex reversal by androgens (Wang et al., 2010). It has been reported that Dmrt1 maintains the normal development of testes in zebrafish but does not affect the development of ovaries (Bradley et al., 2011). All these studies suggest that Dmrt1 plays a role in testicular formation or function maintenance and can also be used as an effective molecular marker for sex identification in O. fasciatus. Sox9 also appears to be involved in sex determination in many species. Its high expression was identified in the testis of Nile tilapia (Melo et al., 2019), and the testis became abnormal when Sox9 was knocked out in male individuals. Our study showed that Sox9 expression was significantly high in the testis and its expression increased with the maturation of the testis, suggesting that Sox9 is closely associated with gonadal development. Although Nanos2 has been reported as a marker gene of germ stem cells, its tissue distribution varies among species. For example, Nanos2 was found to be specifically expressed in the male germ cells of Nile tilapia (O. niloticus) (Li et al., 2014; Jin et al., 2019). However, in orange-spotted grouper (Epinephelus coioides) and tongue sole (C. semilaevis), in addition to its high expression in the testis and ovaries, it was also expressed in small amounts in tissues, such as the brain and liver (Huang et al., 2017; Sun et al., 2017). In the present study, Nanos2 was detected to be significantly highly expressed in the testis and relatively low in the heart, gill, and brain, suggesting testis-specific expression.
Some female-biased expressed genes were also investigated in this study, including ZP4, Bouncer, RNF208, FoxH1, and TOB. The expression levels of these genes in ovaries were significantly higher than those in testes or other tissues. Moreover, the expression levels of these genes in the ovary were also significantly different at different developmental stages, indicating that they might play important roles in female gonad development in O. fasciatus. Based on the information available from a limited number of previous studies, ZP4 is mainly expressed in females in mammals, such as humans, pigs, and cows, and plays a regulatory role in sperm-egg binding, probably as a receptor or an auxiliary receptor (Meczekalski et al., 2015; Zou et al., 2019). Our results also showed a significant female-biased expression of ZP4. It has been reported that Bouncer is a key fertilization factor required for sperm-egg interaction in zebrafish (Herberg et al., 2018) and is a key determinant of species-specific fertilization in fish, which is mainly expressed in the ovary. In this study, the expression patterns of Bouncer were similar to those reported in previous studies, with a significantly lower expression level in the testis than in the ovary. There are few studies on this gene in fish, and no clear conclusion has been reached yet. FoxH1 was reported as a maternally expressed gene in zebrafish (Pogoda et al., 2000; Sirotkin et al., 2000; Christopher et al., 2010), and studies in tilapia (O. niloticus) showed that FoxH1 was expressed in phase I and II oocytes (Yuan et al., 2014). Our results also showed that the expression of FoxH1 was significantly higher in the ovary than in the testis and other tissues, suggesting it might play important roles in the ovary.
It has been reported that the Amh signals through Amh receptors (Amhr) to regulate differentiation and the growth of target cells. In medaka (Oryzias latipes), Amhr2 was detected to be expressed in somatic cells of the gonads and showed no sexually dimorphic expression during gonad development (Klüver et al., 2007). Similarly, in the spotted scat (S. argus), Amhr2 was also reported to be expressed in the both male and female gonads, suggesting that it could plays an important role in the development of gonads in both male and female (Wang et al., 2018). The gene Gsdf is the master sex determination gene in medaka (O. luzonensis) (Myosho et al., 2012) and Nile tilapia (Jiang et al., 2016), where it showed male-specific high expression. However, in the olive flounder (Paralichthys olivaceus) (Liu et al., 2017) and the Nile tilapia (Kaneko et al., 2015), Gsdf was detected to be predominantly expressed in the Sertoli cells and neighboring spermatogonia in testis, but also detectable in the somatic cells surrounding the oogonia cytoplasm of oocytes. In this study, Amhr and Gsdf were also detected to have relatively high expression in the gonads of both 2Y male and female fish, which suggested that they might be involved in gonad formation and maintenance in both male and female. Overall, the qPCR results of selected genes suggested that male-or female-biased genes might be suitable molecular makers during gonadal development of the rock bream.
In conclusion, the present study characterized and identified a number of potential genes and pathways associated with sex determination and gonad development in rock bream (O. fasciatus) using transcriptional analysis. Several sex-biased genes were screened, including male-biased genes, such as Dmrt1, Dmrtb1, Amh, Nanos2, and Sox9, and female-biased genes, such as ZP4, CA4, RNF208, H2A, somatostatin, and TOB. Most of these genes showed gonad-specific (ovary or testis) expression in tissue distribution analysis, and their expression patterns during different gonadal development stages confirmed the sex dimorphic expression, indicating that they are potential molecular markers for the investigation of sex determination and differentiation in rock bream. Further studies on the function of these genes and pathways and their roles during sex determination and differentiation are needed to reveal the complex sex determination and differentiation processes of the fish with the multiple sex chromosome system.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The animal study was reviewed and approved by the Animal Care and Use Committee of Zhejiang Ocean University and the Zhejiang Marine Fisheries Research Institute.
Author Contributions
QZ and DX designed the experiments. HL, QZ, ML, and RC performed the research, while HL, QZ, and DX wrote the manuscript. All the authors have read and approved the final manuscript.
Funding
This work was supported by grants from the Zhejiang Provincial Natural Science Foundation of China (No. LQ20C190008), National Key R&D Program of China (No. 2019YFD0901204), and Key R&D Program of Zhejiang Province (No. 2019C02056).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.791179/full#supplementary-material
References
Adolfi, M. C., Carreira, A. C., Jesus, L. W., Bogerd, J., Funes, R. M., Schartl, M., et al. (2015). Molecular Cloning and Expression Analysis of Dmrt1 and Sox9 during Gonad Development and Male Reproductive Cycle in the Lambari Fish, Astyanax Altiparanae. Reprod. Biol. Endocrinol. 13, 2. doi:10.1186/1477-7827-13-2
Alam, M. A., Kobayashi, Y., Horiguchi, R., Hirai, T., and Nakamura, M. (2008). Molecular Cloning and Quantitative Expression of Sexually Dimorphic Markers Dmrt1 and Foxl2 during Female-To-Male Sex Change in Epinephelus M. Gen. Comp. Endocrinol. 157, 75–85. doi:10.1016/j.ygcen.2008.03.018
Anders, S., and Huber, W. (2012). Differential Expression of RNA-Seq Data at the Gene Level–The DESeq Package. Heidelberg: European Molecular Biology Laboratory.
Bitencourt, J. A., Sampaio, I., Ramos, R. T., Vicari, M. R., Affonso, P. R., and Affonso, P. R. A. M. (2017). First Report of Sex Chromosomes in Achiridae (Teleostei: Pleuronectiformes) with Inferences About the Origin of the Multiple X1X1X2X2/X1X2Y System and Dispersal of Ribosomal Genes in Achirus Achirus. Zebrafish 14, 90–95. doi:10.1089/zeb.2016.1333
Bradley, K. M., Breyer, J. P., Melville, D. B., Broman, K. W., Knapik, E. W., and Smith, J. R. (2011). An SNP-Based Linkage Map for Zebrafish Reveals Sex Determination Loci. G3 (Bethesda) 1, 3–9. doi:10.1534/g3.111.000190
Cai, M., Zou, Y., Xiao, S., Li, W., Han, Z., Han, F., et al. (2019). Chromosome Assembly of Collichthys Lucidus, a Fish of Sciaenidae with a Multiple Sex Chromosome System. Sci. Data 6, 132. doi:10.1038/s41597-019-0139-x
Chen, J. H., He, M. X., Mu, X. J., Yan, B. L., Zhang, J. B., and Jin, S. C. (2015). cDNA Cloning and mRNA Expression Analysis of Sox9 in Scatophagus Argus. Chin. J. Zoolog. 50, 93–102. doi:10.13859/j.cjz.201501012
Christopher, E. S., Tsutomu, A., and Rebecca, D. B. (2010). FoxH1 and Nodal Signaling during Zebrafish Development. Develop. Biol. 344, 418. doi:10.1016/j.ydbio.2009.05.047
Cui, Z., Liu, Y., Wang, W., Wang, Q., Zhang, N., Lin, F., et al. (2017). Genome Editing Reveals Dmrt1 as an Essential Male Sex-Determining Gene in Chinese Tongue Sole (Cynoglossus Semilaevis). Sci. Rep. 7, 42213. doi:10.1038/srep42213
Daniel, B., Remi, H., Alexis, F., and Yann, G. (2005). Large-Scale Temporal Gene Expression Profiling during Gonadal Differentiation and Early Gametogenesis in Rainbow Trout. Biol. Reprod. 73, 959–966. doi:10.1095/biolreprod.105.041830
Du, X., Wang, B., Liu, X., Liu, X., He, Y., Zhang, Q., et al. (2017). Comparative Transcriptome Analysis of Ovary and Testis Reveals Potential Sex-Related Genes and Pathways in Spotted Knifejaw Oplegnathus Punctatus. Gene 637, 203–210. doi:10.1016/j.gene.2017.09.055
He, F.-X., Jiang, D.-N., Huang, Y.-Q., Mustapha, U. F., Yang, W., Cui, X.-F., et al. (2019). Comparative Transcriptome Analysis of Male and Female Gonads Reveals Sex-Biased Genes in Spotted Scat (Scatophagus Argus). Fish. Physiol. Biochem. 45, 1963–1980. doi:10.1007/s10695-019-00693-8
Herberg, S., Gert, K. R., Schleiffer, A., and Pauli, A. (2018). The Ly6/uPAR Protein Bouncer Is Necessary and Sufficient for Species-Specific Fertilization. Science 361, 1029–1033. doi:10.1126/science.aat7113
Hilbold, E., Bergmann, M., Fietz, D., Kliesch, S., Weidner, W., Langeheine, M., et al. (2019). Immunolocalization of Dmrtb1 in Human Testis with Normal and Impaired Spermatogenesis. Andrology 7, 428–440. doi:10.1111/andr.12617
Huang, J., Li, Y., Shao, C., Wang, N., and Chen, S. (2017). Identification, Characterization and Functional Analysis of Regulatory Region of Nanos Gene from Half-Smooth Tongue Sole ( Cynoglossus Semilaevis ). Gene 617, 8–16. doi:10.1016/j.gene.2017.03.033
Ijiri, S., Kaneko, H., Kobayashi, T., Wang, D.-S., Sakai, F., Paul-Prasanth, B., et al. (2008). Sexual Dimorphic Expression of Genes in Gonads during Early Differentiation of a Teleost Fish, the Nile Tilapia Oreochromis N. Biol. Reprod. 78, 333–341. doi:10.1095/biolreprod.107.064246
Ikenishi, K., and Tanaka, T. S. (2000). Spatio-Temporal Expression of Xenopus Vasa Homolog, XVLG1, in Oocytes and Embryos: The Presence of XVLG1 RNA in Somatic Cells as Well as Germline Cells. Dev. Growth Differ. 42, 95–103. doi:10.1046/j.1440-169x.2000.00493.x
Jia, R., Nie, L.-W., Wang, N., and Wang, J. (2009). Molecular Cloning and Expression Patterns of the Vasa Gene from Rana Nigromaculata (Amphibia: Anura). Zoologia (Curitiba, Impr.) 26, 316–322. doi:10.1590/s1984-46702009000200014
Jiang, D.-N., Yang, H.-H., Li, M.-H., Shi, H.-J., Zhang, X.-B., and Wang, D.-S. (2016). Gsdf Is a Downstream Gene of Dmrt1 that Functions in the Male Sex Determination Pathway of the Nile tilapia. Mol. Reprod. Dev. 83, 497–508. doi:10.1002/mrd.22642
Jin, Y. H., Davie, A., and Migaud, H. (2019). Expression Pattern of Nanos, Piwil, Dnd, Vasa and Pum Genes during Ontogenic Development in Nile tilapia Oreochromis Niloticus. Gene 688, 62–70. doi:10.1016/j.gene.2018.11.078
Kamiya, T., Kai, W., Tasumi, S., Oka, A., Matsunaga, T., Mizuno, N., et al. (2012). A Trans-species Missense SNP in Amhr2 Is Associated with Sex Determination in the Tiger Pufferfish, Takifugu Rubripes (Fugu). Plos Genet. 8, e1002798. doi:10.1371/journal.pgen.1002798
Kaneko, H., Ijiri, S., Kobayashi, T., Izumi, H., Kuramochi, Y., Wang, D.-S., et al. (2015). Gonadal Soma-Derived Factor (Gsdf), a TGF-Beta Superfamily Gene, Induces Testis Differentiation in the Teleost Fish Oreochromis Niloticus. Mol. Cell Endocrinol. 415, 87–99. doi:10.1016/j.mce.2015.08.008
Kenneth, J. L., and Thomas, D. S. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2 −ΔΔCT Method. Methods 25, 402–408. doi:10.1006/meth.2001.1262
Klüver, N., Pfennig, F., Pala, I., Storch, K., Schlieder, M., Froschauer, A., et al. (2007). Differential Expression of Anti-Müllerian Hormone (Amh) and Anti-Müllerian Hormone Receptor Type II (amhrII) in the Teleost Medaka. Dev. Dyn. 236, 271–281. doi:10.1002/dvdy.20997
Krysanov, E., and Demidova, T. (2018). Extensive Karyotype Variability of African Fish Genus Nothobranchius (Cyprinodontiformes). Comp. Cytogenet. 12, 387–402. doi:10.3897/compcytogen.v12i3.25092
Kuiri-Hänninen, T., Kallio, S., Seuri, R., Tyrväinen, E., Liakka, A., Tapanainen, J., et al. (2011). Postnatal Developmental Changes in the Pituitary-Ovarian axis in Preterm and Term Infant Girls. J. Clin. Endocrinol. Metab. 96, 3432–3439. doi:10.1210/jc.2011-1502
Langmead, B., and Salzberg, S. L. (2012). Fast Gapped-Read Alignment with Bowtie 2. Nat. Methods 9, 357–359. doi:10.1038/nmeth.1923
Li, B., and Dewey, C. N. (2011). RSEM: Accurate Transcript Quantification from RNA-Seq Data with or without a Reference Genome. BMC Bioinformatics 12, 323. doi:10.1186/1471-2105-12-323
Li, L., Liang, H. W., Li, Z., Luo, X. Z., Zhang, Z. W., Zhu, Y., et al. (2012). Full Length Cloning and Expression Analysis of Dmrt1 Gene cDNA of Pelteobagrus Fulvidraco. J. Huazhong Agric. Univ. 31, 220–226. doi:10.13300/j.cnki.hnlkxb.2012.02.019
Li, M., Xu, H., Xu, W., Zhou, Q., Xu, X., Zhu, Y., et al. (2020). Isolation of a Male-Specific Molecular Marker and Development of a Genetic Sex Identification Technique in Spotted Knifejaw (Oplegnathus Punctatus). Mar. Biotechnol. 22, 467–474. doi:10.1007/s10126-020-09966-3
Li, M., Yang, H., Zhao, J., Fang, L., Shi, H., Li, M., et al. (2014). Efficient and Heritable Gene Targeting in tilapia by CRISPR/Cas9. Genetics 197, 591–599. doi:10.1534/genetics.114.163667
Liu, Y., Zhang, W., Du, X., Zhao, J., Liu, X., Li, X., et al. (2017). Sexually Dimorphic Expression in Developing and Adult Gonads Shows an Important Role of Gonadal Soma-Derived Factor during Sex Differentiation in Olive Flounder ( Paralichthys O ). Comp. Biochem. Physiol. B: Biochem. Mol. Biol. 210, 1–8. doi:10.1016/j.cbpb.2017.05.003
Luo, H., and Dearolf, C. R. (2001). The JAK/STAT Pathway andDrosophila Development. BioEssays 23, 1138–1147. doi:10.1002/bies.10016
Meczekalski, B., Nawrot, R., Nowak, W., Czyzyk, A., Kedzia, H., and Gozdzicka-Jozefiak, A. (2015). Study on the Zona Pellucida 4 (ZP4) Gene Sequence and its Expression in the Ovaries of Patients with Polycystic Ovary Syndrome. J. Endocrinol. Invest. 38, 791–797. doi:10.1007/s40618-015-0260-4
Melo, L. H., Melo, R. M. C., Luz, R. K., Bazzoli, N., and Rizzo, E. (2019). Expression of Vasa, Nanos2 and Sox9 during Initial Testicular Development in Nile tilapia (Oreochromis Niloticus) Submitted to Sex Reversal. Reprod. Fertil. Dev. 31, 1637–1646. doi:10.1071/rd18488
Myosho, T., Otake, H., Masuyama, H., Matsuda, M., Kuroki, Y., Fujiyama, A., et al. (2012). Tracing the Emergence of a Novel Sex-Determining Gene in Medaka, Oryzias L. Genetics 191, 163–170. doi:10.1534/genetics.111.137497
Pogoda, H.-M., Solnica-Krezel, L., Driever, W., and Meyer, D. (2000). The Zebrafish Forkhead Transcription Factor FoxH1/Fast1 Is a Modulator of Nodal Signaling Required for Organizer Formation. Curr. Biol. 10, 1041–1049. doi:10.1016/s0960-9822(00)00669-2
Ricardo, S. H., Yu, M., Miho, O., Shuji, M., Sullip, K. M., Takashi, S., et al. (2012). A Y-Linked Anti-Mullerian Hormone Duplication Takes over a Critical Role in Sex Determination. Proc. Natl. Acad. Sci. 109, 2955–2959. doi:10.1073/pnas.1018392109
Robert, H. D., and Nagahama, Y. (2002). Sex Determination and Sex Differentiation in Fish: An Overview of Genetic, Physiological, and Environmental Influences. Aquaculture 208, 191–364. doi:10.1016/S0044-8486(02)00057-1
Rodríguez-Marí, A., Yan, Y. L., Ruth, A. B., Catherine, W., Cañestro, C., and John, H. P. (2005). Characterization and Expression Pattern of Zebrafish Anti-Mullerian Hormone (Amh) Relative to Sox9a, Sox9b, and Cyp19a1a during Gonad Development. Gene Expr. Patterns 5, 655–667. doi:10.1016/j.modgep.2005.02.008
Schmidt, K. L., Kryger-Baggesen, N., Byskov, A. G., and Andersen, C. Y. (2004). Anti-Müllerian Hormone Initiates Growth of Human Primordial Follicles In Vitro. Mol. Cel Endocrinol 234, 87–93. doi:10.1016/j.mce.2004.12.010
Sember, A., Bohlen, J., Šlechtová, V., Altmanová, M., Symonová, R., and Ráb, P. (2015). Karyotype Differentiation in 19 Species of River Loach Fishes (Nemacheilidae, Teleostei): Extensive Variability Associated with rDNA and Heterochromatin Distribution and its Phylogenetic and Ecological Interpretation. BMC Evol. Biol. 15, 251–272. doi:10.1186/s12862-015-0532-9
Shin, Y., Jung, M., Shin, G.-H., Jung, H.-J., Baek, S.-J., Lee, G.-Y., et al. (2018). First Draft Genome Sequence of the Rock Bream in the Family Oplegnathidae. Sci. Data. 5, 180234. doi:10.1038/sdata.2018.234
Sirotkin, H. I., Gates, M. A., Kelly, P. D., Schier, A. F., and Talbot, W. S. (2000). Fast1 Is Required for the Development of Dorsal Axial Structures in Zebrafish. Curr. Biol. 10, 1051–1054. doi:10.1016/s0960-9822(00)00679-5
Sissao, R., D’Cotta, H., Baroiller, J.-F., and Toguyeni, A. (2019). Mismatches between the Genetic and Phenotypic Sex in the Wild Kou Population of Nile tilapia Oreochromis N. PeerJ 7, e7709. doi:10.7717/peerj.7709
Sun, Z.-H., Wang, Y., Lu, W.-J., Li, Z., Liu, X.-C., Li, S.-S., et al. (2017). Divergent Expression Patterns and Function Implications of Four Nanos Genes in a Hermaphroditic Fish, Epinephelus Coioides. Int. J. Mol. Sci. 18, 685. doi:10.3390/ijms18040685
Tao, W., Chen, J., Tan, D., Yang, J., Sun, L., Wei, J., et al. (2018). Transcriptome Display during tilapia Sex Determination and Differentiation as Revealed by RNA-Seq Analysis. BMC genomics 19, 363. doi:10.1186/s12864-018-4756-0
Tao, W., Yuan, J., Zhou, L., Sun, L., Sun, Y., Yang, S., et al. (2013). Characterization of Gonadal Transcriptomes from Nile tilapia (Oreochromis Niloticus) Reveals Differentially Expressed Genes. PLoS One 8, e63604. doi:10.1371/journal.pone.0063604
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript Assembly and Quantification by RNA-Seq Reveals Unannotated Transcripts and Isoform Switching during Cell Differentiation. Nat. Biotechnol. 28, 511–515. doi:10.1038/nbt.1621
Tsuda, M., Sasaoka, Y., Kiso, M., Abe, K., Haraguchi, S., Kobayashi, S., et al. (2003). Conserved Role of Nanos Proteins in Germ Cell Development. Science 301, 1239–1241. doi:10.1126/science.1085222
Wang, D.-S., Zhou, L.-Y., Kobayashi, T., Matsuda, M., Shibata, Y., Sakai, F., et al. (2010). Doublesex- and Mab-3-Related Transcription Factor-1 Repression of Aromatase Transcription, a Possible Mechanism Favoring the Male Pathway in Tilapia. Endocrinology 151, 1331–1340. doi:10.1210/en.2009-0999
Wang, M., Deng, S. P., Chen, H. P., Zhu, C. H., Wu, T. Li., Li, G. L., et al. (2018). Molecular Cloning and Expression Analysis of Anti-Müllerian Hormone Receptor Type II (Amhr2) in Spotted Scat (Scatophagus argus). J. Guangdong Ocean Univ. 38, 17–24 [In Chinese]. doi:10.3969/j.issn.1673-9159.2018.03.003
Xiao, Y., Xiao, Z., Ma, D., Liu, J., and Li, J. (2019). Genome Sequence of the Barred Knifejaw Oplegnathus F (Temminck & Schlegel, 1844): The First Chromosome-Level Draft Genome in the Family Oplegnathidae. Gigascience 8, giz013. doi:10.1093/gigascience/giz013
Xiao, Y., Xiao, Z., Ma, D., Zhao, C., Liu, L., Wu, H., et al. (2020). Chromosome-level Genome Reveals the Origin of Neo-Y Chromosome in the Male Barred Knifejaw Oplegnathus Fasciatus. iScience 23, 101039. doi:10.1016/j.isci.2020.101039
Xu, D., Lou, B., Bertollo, L. A. C., and Cioffi, M. D. B. (2013). Chromosomal Mapping of Microsatellite Repeats in the Rock Bream Fish Oplegnathus Fasciatus, with Emphasis of Their Distribution in the Neo-Y Chromosome. Mol. Cytogenet. 6, 12. doi:10.1186/1755-8166-6-12
Xu, D., Lou, B., Xu, H., Li, S., and Geng, Z. (2013). Isolation and Characterization of Male-Specific DNA Markers in the Rock Bream Oplegnathus Fasciatus. Mar. Biotechnol. 15, 221–229. doi:10.1007/s10126-012-9480-1
Xu, D., Sember, A., Zhu, Q., Oliveira, E. A. d., Liehr, T., Al-Rikabi, A. B. H., et al. (2019). Deciphering the Origin and Evolution of the X1X2Y System in Two Closely-Related Oplegnathus Species (Oplegnathidae and Centrarchiformes). Int. J. Mol. Sci. 20, 3571. doi:10.3390/ijms20143571
Xu, D., You, F., Lou, B., Geng, Z., Li, J., and Xiao, Z. Z. (2012). Comparative Analysis of Karyotype and C-Banding in Male and Female Oplegnathus Fasciatus. Acta Hydrobiol. Sin. 36, 552–557. doi:10.3724/SP.J.1035.2012.00552
Xue, R., An, H., Liu, Q. H., Xiao, Z. Z., Wang, Y. F., and Li, J. (2016). Karyotype and Ag-NORs in Male and Female of Oplegnathus Punctatus. Ocean. Limnol. Sin. 47, 626–632. doi:10.11693/hyhz20160300042
Yan, H., Shen, X., Cui, X., Wu, Y., Wang, L., Zhang, L., et al. (2018). Identification of Genes Involved in Gonadal Sex Differentiation and the Dimorphic Expression Pattern in Takifugu Rubripes Gonad at the Early Stage of Sex Differentiation. Fish. Physiol. Biochem. 44, 1275–1290. doi:10.1007/s10695-018-0519-8
Yano, A., Guyomard, R., Nicol, B., Jouanno, E., Quillet, E., Klopp, C., et al. (2012). An Immune-Related Gene Evolved into the Master Sex-Determining Gene in Rainbow trout, Oncorhynchus Mykiss. Curr. Biol. 22, 1423–1428. doi:10.1016/j.cub.2012.05.045
Yoshizaki, G., Sakatani, S., Tominaga, H., and Takeuchi, T. (2000). Cloning and Characterization of Avasa-Like Gene in Rainbow trout and its Expression in the Germ Cell Lineage. Mol. Reprod. Dev. 55, 364–371. doi:10.1002/(sici)1098-2795(200004)55:4<364::aid-mrd2>3.0.co;2-8
Yuan, J., Tao, W., Cheng, Y., Huang, B., and Wang, D. (2014). Genome-Wide Identification, Phylogeny, and Gonadal Expression of Fox Genes in Nile tilapia, Oreochromis Niloticus. Fish. Physiol. Biochem. 40, 1239–1252. doi:10.1007/s10695-014-9919-6
Zamboni, L., and Upadhyay, S. (1983). Germ Cell Differentiation in Mouse Adrenal Glands. J. Exp. Zool. 228, 173–193. doi:10.1002/jez.1402280204
Zeng, Q., Liu, S., Yao, J., Zhang, Y., Yuan, Z., Jiang, C., et al. (2016). Transcriptome Display during Testicular Differentiation of Channel Catfish (Ictalurus Punctatus) as Revealed by RNA-Seq Analysis. Biol. Reprod. 95, 19. doi:10.1095/biolreprod.116.138818
Zhang, L., Bian, L., Chang, Q., Hou, J. L., Chen, S. Q., Zhao, Q., et al. (2018). Cloning and Expression of the Sox9 Gene in Spotted Halibut (Verasper Variegatus). Prog. Fish. Sci. 39, 72–80. doi:10.19663/j.issn2095-9869.20171010001
Zhang, S., Zheng, J., Zhang, J., Wang, Z., Wang, Y., and Cai, M. (2018). Cytogenetic Characterization and Description of an X1X1X2X2/X1X2Y Sex Chromosome System in Collichthys Lucidus (Richardson, 1844). Acta Oceanol. Sin. 37, 34–39. doi:10.1007/s13131-018-1152-1
Zhang, Y.-Q., Wen, H.-S., He, F., Li, J.-F., Shi, D., Hu, J., et al. (2012). Cloning of Vasa and Expression Analysis during the Reproductive Cycle inCynoglossus semilaevisGünther. J. Fish. China 36, 1–8. (In Chinese). doi:10.3724/sp.j.1231.2012.27584
Zhang, Z., Zhu, B., Chen, W., and Ge, W. (2020). Anti-Müllerian Hormone (Amh/amh) Plays Dual Roles in Maintaining Gonadal Homeostasis and Gametogenesis in Zebrafish. Mol. Cel Endocrinol 517, 110963. doi:10.1016/j.mce.2020.110963
Keywords: gonadal development, RNA-seq, histological observation, rock bream, sex-related genes
Citation: Li H, Zhu Q, Chen R, Liu M and Xu D (2021) Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes. Front. Genet. 12:791179. doi: 10.3389/fgene.2021.791179
Received: 08 October 2021; Accepted: 09 November 2021;
Published: 29 November 2021.
Edited by:
Chutian Ge, Zhejiang Wanli University, ChinaReviewed by:
Xiao Ling, Zhejiang Wanli University, ChinaJie Mei, Huazhong Agricultural University, China
Li Shuisheng, Sun Yat-sen University, China
Copyright © 2021 Li, Zhu, Chen, Liu and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dongdong Xu, xudong0580@163.com
†These authors have contributed equally to this work