Skip to main content

Formation of autotriploid Carassius auratus and its fertility-related genes analysis



Formation of triploid organism is useful in genetics and breeding. In this study, autotriploid Carassius auratus (3nRR, 3n = 150) was generated from Carassius auratus red var. (RCC, 2n = 100) (♀) and autotetraploid Carassius auratus (4nRR, 4n = 200) (♂). The female 3nRR produced haploid, diploid and triploid eggs, whereas the male 3nRR was infertile. The aim of the present study was to explore fertility of potential candidate genes of 3nRR.


Gonadal transcriptome profiling of four groups (3 females RCC (FRCC), 3 males 4nRR (M4nRR), 3 males 3nRR (M3nRR) and 3 females 3nRR (F3nRR)) was performed using RNA-SEq. A total of 78.90 Gb of clean short reads and 24,262 differentially expressed transcripts (DETs), including 20,155 in F3nRR vs. FRCC and 4,107 in M3nRR vs. M4nRR were identified. A total of 106 enriched pathways were identified through KEGG enrichment analysis. Out of the enriched pathways, 44 and 62 signalling pathways were identified in F3nRR vs. FRCC and M3nRR vs. M4nRR, respectively. A total of 80 and 25 potential candidate genes for fertility-related in F3nRR and M3nRR were identified, respectively, through GO, KEGG analyses and the published literature. Moreover, protein-protein interaction (PPI) network construction of these fertility-associated genes were performed. Analysis of the PPI networks showed that 6 hub genes (MYC, SOX2, BMP4, GATA4, PTEN and BMP2) were involved in female fertility of F3nRR, and 2 hub genes (TP53 and FGF2) were involved in male sterility of M3nRR.


Establishment of autotriploid fish offers an ideal model to study reproductive traits of triploid fish. RNA-Seq data revealed 6 genes, namely, MYC, SOX2, BMP4, GATA4, PTEN and BMP2, involved in the female fertility of the F3nRR. Moreover, 2 genes, namely, TP53 and FGF2, were related to the male sterility of the M3nRR. These findings provide information on reproduction and breeding in triploid fish.

Peer Review reports


Polyploid organisms have three or more chromosome sets. Triploidy, an example of polyploids, plays a vital role in the process of biological evolution and can be divided into autotriploidy and allotriploidy [1]. Allotriploids have three chromosome sets from two or more different species, whereas autotriploids have three chromosome sets derived from a single taxon.

Development of gonads is critical to fertility in sexually reproducing organisms especially in triploids and is tightly regulated by complex processes [2]. Sex determination, sexual differentiation and gametogenesis are important processes during gonadal development. Any abnormality in these events can result in infertility. Several genes implicated in sexual determination and differentiation have been reported [3,4,5]. Gametogenesis, including oogenesis and spermatogenesis, are also regulated by complex mechanisms and several regulatory genes [6, 7]. A previous study explored regulation of early stages of oogenesis [8]. In addition, studies explored biological mechanisms that occur mid-oogenesis [9], and regulation of late oogenesis [10]. Spermatogenesis is divided into three steps: spermatogonial mitotic proliferation, two times of meiosis, and post-meiotic differentiation [11]. Previous studies explored the functional mechanisms of spermatocytogenesis [12], meiosis during spermatogenesis [13], and spermiogenesis [14].

Fertility of polyploids has important implications in fisheries and sustainable aquaculture. Artificial triploids of species such as Atlantic salmon, Oncorhynchus mykiss, Salmo trutta and Salvelinus fontinalis have been widely used in fish farming industry [15]. There has been a believe that triploidization causes infertility in fish. In addition, a previous study reports that the triploid can channel the energy required for gonad maturation to somatic growth, causing rapid growth rates compared with their diploid counterparts [16]. However, a different study reported that triploid fish can produce normal gametes [17]. In our previous study, Carassius auratus red var. (RCC) (female) and autotetraploid Carassius auratus (4nRR) (male) were artificially hybridized to produce hybrid autotriploid Carassius auratus (3nRR). After hybridization, the male 3nRR did not produce normal sperm, whereas the females generated dynamic eggs [18]. Analysis of meiosis-related gene expression showed that Dmc1 and Ph1 had higher expression level in female 3nRR compared with levels in the males, indicating that these genes are involved in regulating fertility of 3nRR [19]. Molecular mechanisms involving 3nRR fish in controlling fertility have not been explored fully.

RNA-Seq technology is utilised for analysis of the structure and function of genes at the organismal level, and for exploring a series of biological pathways [20]. RNA-Seq technique has been successfully used in studies various fishes in the past decade. In spotted scat species (Scatophagus argus), several candidate genes involved in reproduction and gonadal development were obtained by RNA-Seq [21]. Studies on Takifugu rubripes reported that sex-related genes play an important role at early sex differentiation stage [22]. Gonadal transcriptome profiling of triploid hybrid loaches (Misgurnus anguillicaudatus) and their diploid and tetraploid parents showed key genes implicated in low hybrid triploid fertility [17]. A study on Thunnus maccoyii reported sex and gonad-development-related genes in the gonads of Southern bluefin tuna through RNA-Seq [23]. In addition, RNA-Seq has been successfully used to analyze sex determination and differentiation related genes in tilapia [24]. In the present study, we successfully obtained triploid fish (3nRR) by crossing female RCC and male 4nRR. The diploid (2nF1), triploid (3nF1) and tetraploid (4nF1) hybrids were then generated by hybridization of female 3nRR and male RCC. In this study, we explored important biological traits and systematically compared gonadal transcriptome of the triploid fish (3nRR) with their parents. Further, the molecular mechanism of the low fertility of the autotriploid fish was explored. The findings of this study provide information on the biological characteristics of 3nRR and mechanisms associated with fertility regulation in triploid fish.


Fertility of autotriploid Carassius auratus

3nRR were generated by crossing female RCC and male 4nRR during the breeding season (Figs. 1a, b and c, 2 and 3a, b and c; Table 1). Testes of RCC and 4nRR (Fig. 4a, b) contained spermatogonia (SG), spermatocytes (SC) and a large number of mature spermatid (ST), whereas the mature sperm was not observed in 3nRR (Fig. 4c). Ovaries of RCC, 4nRR and 3nRR contained second, third and fourth phase oocytes (Fig. 4d, e, f). These results indicated that all ovaries, and RCC and 4nRR testes were fertile whereas 3nRR testes were sterile.

Fig. 1

Formation of polyploid fish

Fig. 2

DNA-content flow-cytometry histograms of RCC (a), 4nRR (b) and 3nRR (c)

Fig. 3

Chromosome spreads at metaphase in RCC, 4nRR, 3nRR, 2nF1, 3nF1 and  4nF1. a: The 100 chromosomes of RCC; b: The 200 chromosomes of 4nRR; c: The 150 chromosomes of 3nRR; d: The 100 chromosomes of 2nF1; e: The 150 chromosomes of 3nF1; f: The 200 chromosomes of 4nF1; bar = 5 μm

Table 1 Examination of chromosome number of RCC, 4nRR, 3nRR, 2nF1, 3nF1 and 4nF1
Fig. 4

Micrographs of the testes and ovaries of RCC, 3nRR and 4nRR. a: Micrographs of testis from RCC; b: Micrographs of testis from 4nRR; c: Micrographs of testis from 3nRR; d: Micrographs of ovary from RCC; e: Micrographs of ovary from 4nRR; f: Micrographs of ovary from 3nRR; SG: spermatogonia; SC: spermatocyte; ST: spermatid; II: stage II oocyte; III: stage III oocyte; IV: stage IV oocyte; Bars = 50 μm

Eggs and water-like semen were collected during the reproductive season from two years old males and females of 3nRR, respectively (Fig. 5). Ploidy levels of the offspring resulting from a cross of female 3nRR and male RCC (Fig. 1c, d, e, f, g) were determined by measuring the chromosome number (Fig. 3d, e, f; Table 1). These analyses showed that female 3nRR produced different sizes of eggs.

Fig. 5

Spermatozoa and eggs of 3nRR. a: Abnormal spermatozoa produced by 3nRR males; b: Different sizes of eggs collected from 2-year-old 3nRR females

Transcriptome sequencing and sequence alignment

Optical density (OD) ratio A260/A280 and RNA integrity numbers (RINs) of the RNA in 12 samples (Table 2) were 2.1 and 8.0-8.8, respectively (Additional file 1). These results indicate that all samples were free from contamination and their quality met the requirements for transcriptome sequencing.

Table 2 Sample information

RNA-seq from gonadal tissue samples of autotriploid fish and their parents was performed by Illumina. RNA-Seq results are presented in Tables 3 and 4. Number of clean reads from the 12 RNA-seq libraries ranged from 39,624,312 to 50,588,484. All clean reads were then aligned to the RCC genome sequences using HISAT2 software. Mapped genome reads ranged from 24,237,536 to 42,474,296, genome map rates ranged from 59.79 to 91.55 %, and unique match rates ranged from 57.84 to 85.93 %.

Table 3 Summary of the RNA-Seq data collected from FRCC, M4nRR, F3nRR and M3nRR
Table 4 Summary of clean reads mapped from FRCC, M4nRR, F3nRR and M3nRR to the reference genome

Identification of Differentially Expressed Transcripts (DETs)

Analysis of F3nRR and FRCC showed that a total of 13,467 DETs were downregulated whereas 6,688 DETs were up-regulated (Fig. 6a). DETs between F3nRR and FRCC included forkhead box L2 (FOXL2), LIM homeobox 8 (LHX8), lysine acetyltransferase 8 (KAT8), BCL2 apoptosis regulator (BCL2), doublesex and mab-3 related transcription factor 1 (DMRT1), ovarian serine protease (OSP) and CCM2 scaffold protein (CCM2). Analysis of M3nRR and M4nRR showed that a total of 1,886 DETs were downregulated and 2,221 DETs were up-regulated (Fig. 6b). DETs between M3nRR and M4nRR included septin 12 (SEPT12), ATPase copper transporting beta (ATP7B), CF transmembrane conductance regulator (CFTR), cAMP responsive element modulator (CREM), cytochrome P450 family 26 subfamily B member 1 (CYP26B1), EF-hand calcium binding domain 2 (EFCAB2) and inhibitor of kappa light polypeptide gene enhancer in B-cells and kinase complex-associated protein (IKBKAP).

Fig. 6

Volcano plot for transcript differential expression. a: F3nRR vs. FRCC; b: M3nRR vs. M4nRR. Transcripts with FDR < 0.05 and ratio of FPKMs of the two samples > 2 were considered to be differentially expressed transcripts. The red region shows significantly up-regulated transcripts, whereas the green region shows down-regulated transcripts

GO and KEGG enrichment analysis of DETs

GO enrichment analysis of the biological process, cellular component and molecular function categories yielded 242, 38 and 51 terms, respectively, for F3nRR vs. FRCC, and 223, 28 and 29 for M3nRR vs. M4nRR group. (Additional files 2 and 3). The most-enriched GO-terms for F3nRR vs. FRCC group were “induction of programmed cell death” in the biological process category, “neuron projection” in the cellular component category, and “channel activity” and “passive transmembrane transporter activity” in the molecular function category. The most-enriched GO-terms for M3nRR vs. M4nRR group were “extracellular region part” in the cellular component category; “kinase activity” and “transferase activity, transferring phosphorus-containing groups” in the molecular function category; and “response to osmotic stress” in the biological process category (Fig. 7).

Fig. 7

Gene Ontology (GO) functional classification of differentially expressed transcripts (DETs). a: F3nRR vs. FRCC; b: M3nRR vs. M4nRR. The three terms are presented on the x-axis shows and the proportion of DETs corresponding to each subcategory is presented on the y-axis

KEGG analysis of all DETs showed that 44 and 62 signaling pathways were enriched in the F3nRR vs. FRCC group and M3nRR vs. M4nRR group, respectively (Additional files 4 and 5). The top 20 most enriched KEGG pathways are shown in Fig. 8. The five most-enriched pathways in the F3nRR vs. FRCC group were “ion channels” (ko04040), “cAMP signaling pathway” (ko04024), “focal adhesion” (ko04510), “glycosaminoglycan binding proteins” (ko00536) and “glycosyltransferases” (ko01003). Moreover, several pathways implicated in female fertility of F3nRR were identified, including “MAPK signaling pathway - plant” (ko04016), and “p53 signaling pathway” (ko04115). The five most enriched pathways for the M3nRR vs. M4nRR group were “Ion channels” (ko04040), “rap1 signaling pathway” (ko04015), “ras signaling pathway” (ko04014), “alcoholism” (ko05034) and “axon guidance” (ko04360). Notably, four of the top 20 most-enriched pathways, “regulation of actin cytoskeleton” (ko04810), “calcium signaling pathway” (ko04020), “tight junction” (ko04530) and “cytokines and growth factors” (ko04052), play important roles in cellular processes such as differentiation, proliferation, migration and apoptosis, implying that they are potentially involved in male sterility of M3nRR.

Fig. 8

Statistics of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the functional significance of DETs. a: F3nRR vs. FRCC; b: M3nRR vs. M4nRR. The abscissa represents P value, which decreases with increase in significance in enrichment level of differentially expressed transcripts in the pathway. The ordinate represents log10 (Q value), which increases with increase in significance of differentially expressed transcripts in the pathway

Hub genes related to the fertility in 3nRR were identified

Eighty genes out of the DETs identified in the F3nRR vs. FRCC group related to female fertility were identified by literature supported searching (Additional file 6). On the other hand, 25 genes out of the DETs in the M3nRR vs. M4nRR group are implicated in male sterility (Additional file 7). To further identify hub genes associated with 3nRR fertility, PPI of the fertility-related genes was constructed using STRING tool and analysis was carried out using Cytoscape software. After analysis of PPI network of female fertility-related genes, 6 genes with the interaction degrees more than 15 were screened as hub genes (Fig. 9a, Additional file 8). Furthermore, PPI of male sterility-related genes showed that 2 hub genes, with degrees more than 5 showed strong interaction with other node proteins (Fig. 9b, Additional file 9).

Fig. 9

Protein-protein interaction network of fertility-related genes. a: F3nRR vs. FRCC; b: M3nRR vs. M4nRR. Genes with higher degree value are presented in deep red color

RT-qPCR verification

To verify RNA-Seq results, twenty-eight DETs were chosen for validation by RT-qPCR. Among the 28 DETs, 6 DETs and 7 DETs were up-regulated in the F3nRR vs. FRCC and M3nRR vs. M4nRR groups, respectively; whereas 10 DETs and 5 DETs were down-regulated in the F3nRR vs. FRCC and M3nRR vs. M4nRR groups, respectively (Fig. 10). Expression profiles of the twenty DEGs obtained by RT-qPCR and RNA-Seq were similar, implying that RNA-Seq results were reliable.

Fig. 10

qPCR analysis of selected DETs. a: F3nRR vs. FRCC; b: M3nRR vs. M4nRR. Data represent the means ± SD, n = 3 independent experiments. *p < 0.01 versus control


Triploid animals are usually sterile and cannot form triploid populations. However, previous studies have been reporting contradicting results. Xiao et al. [25] reported that triploid Carassius auratus in Dongting water system produces normal gametes. Hu et al. [26] reported that female autotriploid hybrids (3nAUT) generated by crossing females of Carassius auratus red var. with males of autotetraploid fish produced mature eggs. However, male 3nAUT showed abnormal gonadal development and could not produce mature sperm. In the present study, 3nRR was generated by crossing female RCC and male 4nRR. Three different ploidy offspring were then obtained by hybridization of female 3nRR and male RCC. The female 3nRR offspring were fertile, whereas male 3nRR were sterile. However, only a few studies have explored the molecular mechanisms modulating fertility of the autotriploid of Carassius auratus. In this study transcriptome analysis was used to explore the molecular mechanisms associated with poor fertility in 3nRR. Eight fertility-related hub genes of 3nRR were identified through GO and KEGG enrichment analyses, and previous published literature.

Candidate hub genes related to male sterility of 3nRR were identified

Hub genes identified in the M3nRR vs. M4nRR group included several genes involved in male sterility, such as the tumor protein p53 (TP53) and fibroblast growth factor 2 (FGF2).

TP53, also known as P53 is a transcriptional regulator and tumor suppressor implicated in spermatogenesis [27]. In vertebrates, partial or complete impairment of P53 expression causes disordered meiotic divisions, which in turn causes spermatogenesis defects [28, 29]. P53 mRNA and protein levels are downregulated in the testis of P53 promoter-chloramphenicol acetyltransferase (CAT)-harboring mice, indicating its important role in development of spermatocytes [30]. In addition, TP53 codon 72 polymorphism in mice is involved in meiosis, implying that it plays a critical role in spermatogenesis [31]. In human, P53 gene polymorphism is higher in infertile men compared with fertile men, implying that it may affect germ cell apoptosis and increase risk of male infertility [29, 32]. In the current study, analysis of expression levels of the TP53 gene showed significantly different expression levels between M3nRR and M4nRR. This finding implies that TP53 may disrupt meiosis during spermatogenesis in the male 3nRR fish causing sterility.

Fibroblast growth factor 2 (FGF2) plays essential functions in regulation of spermatogenesis and sperm physiology [33]. A study using a human model reported presence of FGF2 and FGFRs in testis and sperm, which are related with human spermatogenesis and sperm motility [34]. Furthermore, incubation of human sperm with recombinant FGF2 (rFGF2) causes an increase in number of motile cells, implying that the gene is involved in sperm motility [35]. In mouse, knock out of FGF2 induces impaired sperm production and is associated with alterations in sperm morphology and function [36]. In this study, FGF2 was significantly upregulated in M3nRR vs. 4nRR. High expression levels of the gene can cause abnormal shaping of the normal sperm, which resulted in male 3nRR sterility.

Candidate hub genes related to female fertility of 3nRR were identified

Six hub genes associated with female fertility were identified in the F3nRR vs. F4nRR group including MYC proto-oncogene, bHLH transcription factor (MYC), SRY-box transcription factor 2 (SOX2), bone morphogenetic protein 4 (BMP4), GATA binding protein 4 (GATA4), phosphatase and tensin homolog (PTEN) and bone morphogenetic protein 2 (BMP2).

MYC gene encodes the MYC transcription factor which is involved in cell proliferation and gametogenesis [37]. In Xenopus, C-MYC was detected in oocytes, indicating that it plays a role in oogenesis [38]. In Drosophila, MYC was involved in distribution of CTPsyn in follicle cells, implying that it plays a role in synthesizing nutrients for the developing oocytes [39]. In addition, a previous study reported that MYC plays an important role in regulation of mitochondrial biogenesis in Drosophila ovary, indicating that the gene is involved in oocyte development [40]. In Larimichthys crocea, Lc-cMYC had different expression patterns in oocytes at various stages of development, implying that it plays an essential role in oogenesis [41]. In the current study, 3nRR ovary showed low expression level of MYC which may inhibited formation of oogenesis defects during triploidization, resulting in production mature eggs.

SOX is an ancient gene family involved in oogenesis [42]. Sox genes have been explored in many organisms and can be classified into ten subgroups (A-J) [43]. In Agasicles hygrophila, AhDichaete and AhSox3 expression levels are significantly high in ovary, indicating that it plays a vital regulatory role in during ovarian development and oogenesis [44]. In Misgurnus anguillicaudatus, MaSOX3 is abundant in primary oocytes and previtellogenic oocyte cells, indicating that MaSox3 gene is involved in ovarian development [45]. In Paramisgurnus dabryanus, SOX4 was detected in the ovary, showing that it plays an important role during ovarian development [46]. In mouse, expression of SOX2 is required for establishment and maintenance of the oocyte cell [47]. In our study, SOX2 was identified in F3nRR vs. FRCC, with lower expression level in F3nRR gonad compared with the expression level in FRCC gonad. This finding shows that SOX2 may be an important factor in normal ovarian development of F3nRR.

Bone morphogenetic proteins (BMPs) are belonging to the transforming growth factor-β superfamily of proteins, and they appear to be highly conserved [48]. A previous study reports that BMPs play a role in regulation of ovarian follicular development [49]. BMP1, BMP6 and BMP15 are implicated in ovarian development [50,51,52]. In Xenopus laevis, BMP2 gene is highly expressed during oogenesis, implying that it is an important factor in ovarian development [53]. In mouse, BMP4 regulates the number of oocytes, suggesting its role in the process of oogenesis [54]. Roles of BMP2 and BMP4 as important factors in survival and development of bovine secondary follicles were recently reported [55]. In this study, BMP2 and BMP4 genes were differentially expressed between F3nRR and FRCC, indicating that they may be involved in ovarian development.

GATA4, a member of the GATA-binding family, is highly expressed in ovarian granulosa cells [56, 57]. A previous study reported that the gene is involved in regulation of ovarian development [58]. GATA4 and GATA6 knockout female mice exhibited infertility due to disrupted formation of ovaries [59]. GATA4 deletion resulted in a sterile female mice phenotype attributed to drastic reduction in number of developing follicles [60]. The mRNA for GATA4 has been reported in human ovary implying that GATA4 plays a role in ovarian folliculogenesis [57]. GATA4 was identified through transcriptome analysis in this study. GATA4 was differentially expressed in F3nRR vs. FRCC, indicating that plays important roles in ovarian development of 3nRR.

Phosphatase and tensin homolog (PTEN) protein has phosphatase activity and belongs to protein-tyrosine phosphatase superfamily [61]. PTEN is a negative regulator of PI3K-Akt signaling pathway which is involved in growth of eggs [62]. Deletion of PTEN from oocytes affects mouse fertility by interrupting oocyte growth [63]. PTEN signaling pathway associated with ovarian follicle development has been reported in human [64]. In Crassostrea gigas, PTEN is involved in insulin pathway in gonads and plays a critical role in reproduction [65]. In Drosophila, loss of PTEN is related to IIS/mTORC1 signalling, which is important for oogenesis [66]. In the present study, F3nRR showed significantly low expression levels of PTEN, which may have caused F3nRR normal fertility through regulation of oogenesis.


The autotriploid Carassius auratus (3nRR, 3n = 150) is generated from Carassius auratus red var. (RCC, 2n = 100) (♀) and autotetraploid Carassius auratus (4nRR, 4n = 200) (♂), of which the female 3nRR can produce mature gametes, whereas the male 3nRR cannot. In addition, we produced diploid (2nF1, 2n = 100), triploid (3nF1, 3n = 150) and tetraploid (4nF1, 4n = 200) hybrids in the F1 generation by crossing females of 3nRR with males of RCC, which further indicated that female 3nRR were fertile. Gonadal transcriptome reveals 6 hub genes (MYC, SOX2, BMP4, GATA4, PTEN and BMP2) were involved in the female fertility of the female 3nRR, and 2 hub genes (TP53 and FGF2) were involved in the male sterility of the male 3nRR. The obtained data reveals novel candidate genes for the fertility in the autotriploid fish and also extends an understanding of the molecular aspects of fertility in triploid fish.


Animals and crosses

One year old RCC and one year old 4nRR (F11) were fed in the State Key Laboratory of Developmental Biology of Freshwater Fish, Hunan Normal University, China. Hybrids (3nRR) of RCC (♀) × 4nRR (♂) were generated in May 2018. All fish were maintained in open pools (0.067 ha) with suitable pH (7.0-8.5), water temperature (22–24℃), dissolved oxygen content (5.0–8.0 mg/L) and adequate forage. All dissections were performed under MS-222 anaesthesia (100 mg/L; Sigma-Aldrich).

Gonadal histologic analysis

Ploidy levels of the fish (RCC, 4nRR and 3nRR) were estimated using a flow cytometer (Gallios Flow Cytometer, Beckman Coulter). Blood was collected from the caudal vein using heparinized syringes. Samples were then resuspended in 4,6-Diamidino-2-Phenylindole solution (Sigma-Aldrich) for 10 min. DNA content was compared with that of RCC per sample. Gonadal tissues of two years old female RCC, male 4nRR and 3nRR were fixed in Bonn’s liquid and then dehydrated using graded series of alcohol, cleared with xylene, embedded in paraffin wax and cut into 5–8 μm sections. The sections were placed on slides, stained with hematoxylin and eosin, and viewed under a light microscope.

Gamete phenotypes and egg ploidy detection

The water-like semen and mature eggs of two years old 3nRR were sampled for morphological examination. The female 3nRR produced different sized eggs. To determine the egg ploidy, mature eggs were used in vitro fertilization of the RCC haploid sperm and then viable offspring (2nF1, 3nF1 and 4nF1) were generated. Ploidy of these offspring was detected by chromosome counts.

Preparation of chromosome spreads

For ploidy level analysis, chromosome counts were carried out using kidney tissues from 10 individuals each of RCC, 4nRR, 3nRR, 2nF1, 3nF1 and 4nF1 at eight months of age following a previously described method [67]. 200 metaphase chromosome spreads (20 spreads per sample) were analyzed for each type of fish. Each preparation was examined under 3330× magnification with an oil immersion lens.

Sample collection and preparation for transcriptomic sequencing

A total of 3 females RCC (FRCC), 3 males 4nRR (M4nRR), 3 males 3nRR (M3nRR) and 3 females 3nRR (F3nRR) were acquired at 24 months. Fish were anesthetized before surgical removal of tissues. Gonadal tissues were harvested from FRCC, M4nRR, F3nRR and M3nRR after euthanasia.

RNA extraction and sequencing

Total RNA was extracted from gonads of RCC, 4nRR and 3nRR using TRIzol reagent (Takara, Beijing, China) according to the instructions. RNA integrity (RNA integrity score ≥ 7.0) was checked on the bioanalyzer 2100 system (Agilent, Palo Alto, CA) and RNA quantity was measured using NanoDrop 2000 (Thermo, Waltham, MA, USA). A total of 12 libraries from FRCC, M4nRR, F3nRR and M3nRR groups were sequenced. In summary, mRNA was purified and broken into short fragments. Then, reverse transcription, cDNA synthesis and cluster generation were performed. RNA-Seq libraries were then sequenced on Illumina Hiseq2500 platform. The sequenced data are publicly available at the NCBI (PRJNA694292).

Differential expressed transcripts (DETs) and profiling of potential fertility-related genes

After sequencing, clean reads were acquired by removing adapters and low-quality reads using fastp software (Version 0.20.0). The clean reads quality was assessed with FastQC software (Version 0.11.9). The clean reads of the libraries were aligned to the published RCC reference genome ( using HISAT2 tool (Version 2.1.0). To calculate gene expression level, we used fragments per kilobase per million mapped fragments (FPKM) method. Differentially expressed transcripts (DETs) analysis of F3nRR vs. FRCC, and M3nRR vs. M4nRR was performed using DEGSeq2 R package (Version 1.28.1). Transcripts having a fold change (FC) > 2 and a false discovery rate (FDR) < 0.05 were considered as DETs. To further explore these DETs, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed using clusterProfiler (Version 3.6.0) with p < 0.05. Moreover, fertility-related DETs were screened following GO, KEGG enrichment analyses and published literature. STRING database ( to construct protein-protein interaction (PPI) networks to explore protein relationships among the fertility-related DETs. Hub genes were obtained based on the ranking order of connectivity degree by Cytoscape software [68].

Quantitative real-time PCR verification

Ten significantly DETs (five up-regulated DETs and five down-regulated DETs) and eight important genes in this study were chosen for quantitative real-time (qRT) PCR to test the reliability of the F3nRR vs. FRCC and M3nRR vs. M4nRR transcriptome sequencing results. PrimeScript™ RT reagent kit (Takara, Dalian, China) was used to perform cDNA synthesis following the manufacturer’s instructions. Primer sequences for β-actin (the internal control gene) and these DETs are listed in Additional file 10. The 10-µl-volume qRT-PCR reaction mixture consisted of 5 µl SYBR Green qPCR Master Mix, 0.5 µl of 20 µM of each primer, 1 µl of cDNA (1:10 dilution) and 3 µL of nuclease-free water. qRT-PCR thermal cycle used was as follows: 95 °C for 2 min, 40 cycles of 95 °C for 15 s and annealing at 60 °C for 30 s. Three technical replicates were used for each biological sample in the qRT-PCR. Relative mRNA expression level was calculated by using the 2−ΔΔCt method. Data were analyzed statistically using SPSS (v22.0) software (SPSS Inc., Chicago, IL, USA). Statistical significance was determined using Student’s t-test analysis.

Availability of data and materials

Raw sequence reads are available from the NCBI (PRJNA694292) ( and the expression profiles of RNA-Seq data are included in Figshare (



Carassius auratus red var


Autotriploid Carassius auratus


Autotetraploid Carassius auratus


Female Carassius auratus red var


Male autotetraploid Carassius auratus


Female autotriploid Carassius auratus


Male autotriploid Carassius auratus


Protein-protein interaction

2nF1 :

Diploid hybrids of female autotriploid Carassius auratus × male Carassius auratus red var.;

3nF1 :

Triploid hybrids of female autotriploid Carassius auratus × male Carassius auratus red var

4nF1 :

Tetraploid hybrids of female autotriploid Carassius auratus × male Carassius auratus red var








Differentially expressed transcripts


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes

TP53 :

Tumor protein p53

FGF2 :

Fibroblast growth factor 2


MYC proto-oncogene, bHLH transcription factor

SOX2 :

SRY-box transcription factor 2

BMP4 :

Bone morphogenetic protein 4


GATA binding protein 4


Phosphatase and tensin homolog

BMP2 :

Bone morphogenetic protein 2


  1. 1.

    Otto SP, Whitton J. Polyploid incidence and evolution. Annu Rev Genet. 2000;34(1):401–37.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Phizicky DV, Berchowitz LE, Bell SP. Multiple kinases inhibit origin licensing and helicase activation to ensure reductive cell division during meiosis. Elife. 2018;7:e33309.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Nagahama Y. Molecular mechanisms of sex determination and gonadal sex differentiation in fish. Fish Physiol Biochem. 2005;31(2–3):105–9.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Quinn A, Koopman P. The molecular genetics of sex determination and sex reversal in mammals. Semin Reprod Med. 2012;30(5):351–63.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Windley SP, Wilhelm D. Signaling pathways involved in mammalian sex determination and gonad development. Sex Dev. 2015;9(6):297–315.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Schultz N, Hamra FK, Garbers DL. A multitude of genes expressed solely in meiotic or postmeiotic spermatogenic cells offers a myriad of contraceptive targets. Proc Natl Acad Sci. 2003;100:12201–6.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Lee CS, Lu T, Seydoux G. Nanos promotes epigenetic reprograming of the germline by down-regulation of the thap transcription factor lin-15b. Elife. 2017;6:e30201.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Marco N, Assa Y, Eckmann CR. Stage-specific combinations of opposing poly(a) modifying enzymes guide gene expression during early oogenesis. Nucleic Acids Res. 2019;20:20.

    Google Scholar 

  9. 9.

    Kong J, Han H, Bergalet J, Bouvrette LPB, Hernández G, Moon NS, et al. A ribosomal protein S5 isoform is essential for oogenesis and interacts with distinct RNAs in Drosophila melanogaster. Sci Rep. 2019;9(1):1–11.

    Google Scholar 

  10. 10.

    Kawaguchi S, Ueki M, Kai T. Drosophila MARF1 ensures proper oocyte maturation by regulating nanos expression. PLoS One. 2020;15(4):e0231114.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Snyder E, Chukrallah L, Seltzer K, Goodwin L, Braun RE. ADAD1 and ADAD2, testis-specific adenosine deaminase domain-containing proteins, are required for male fertility. Sci Rep. 2020;10(1):11536.

  12. 12.

    Phillips BT, Gassei K, Orwig KE. Spermatogonial stem cell regulation and spermatogenesis. Phil Trans R Soc B. 2010;365:1663–78.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Wolgemuth DJ, Roberts SS. Regulating mitosis and meiosis in the male germ line: critical functions for cyclins. Phil Trans R Soc B. 2010;365:1653–62.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Bettegowda A, Wilkinson MF. Transcription and post-transcriptional regulation of spermatogenesis. Phil Trans R Soc B. 2010;365:1637–51.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Piferrer F, Beaumont A, Falguiere JC, Flajshans M, Haffray P, Colombo L. Polyploid fish and shellfish: production, biology and applications to aquaculture for performance improvement and genetic containment. Aquaculture. 2009;293:125–56.

    Article  Google Scholar 

  16. 16.

    Cal RM, Vidal S, Gómez C, Álvarez-Blázquez B, Martínez P, Piferrer F. Growth and gonadal development in diploid and triploid turbot (Scophthalmus maximus). Aquaculture. 2006;251(1):99–108.

    CAS  Article  Google Scholar 

  17. 17.

    He Z, Xu QZ, Rui Z, Zhuang ZX, Ma YQ, Wei W, et al. Gonadal transcriptome analysis of hybrid triploid loaches (Misgurnus anguillicaudatus) and their diploid and tetraploid parents. PLoS One. 2018;13(5):e0198179.

    Article  CAS  Google Scholar 

  18. 18.

    Qin QB, Wang YD, Wang J, Dai J, Xiao J, Hu FZ, et al. The autotetraploid fish derived from hybridization of Carassius auratus red var. (female) × Megalobrama amblycephala (male). Biol Reprod. 2014;91(4):93.

    PubMed  Article  CAS  Google Scholar 

  19. 19.

    Qin QB, Zhou YW, Wang CQ, Zhang MH, Qin H, Zhao C, et al. Analysis on the meiosis-related gene (Dmc1, Ph1) expression in autotriploid Carassius auratus. Mar Biotechnol. 2019;21:753–61.

    CAS  Article  Google Scholar 

  20. 20.

    Qian X, Ba Y, Zhuang QF, Zhong GF. RNA-Seq technology and its application in fish transcriptomics. OMICS. 2014;18(2):98–110.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Yao Q, Dong YZ, Chen J, Quan LF, Zhang WQ, Chen BX. Transcriptome analysis of female and male Conopomorpha sinensis (Lepidoptera: Gracilariidae) adults with a focus on hormone and reproduction. J Econ Entomol. 2019;112(6):2966–75.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Yan HW, Shen XF, Cui X, Wu YW, Wang LS, Zhang L, et al. 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. 2018;44(5):1275–90.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Bar I, Cummins S, Elizur A. Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the Southern bluefin tuna (Thunnus maccoyii). BMC Genomics. 2016;17:217.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  24. 24.

    Tao WJ, Chen JL, Tan DJ, Yang J, Sun L, Wei J, et al. Transcriptome display during tilapia sex determination and differentiation as revealed by RNA-Seq analysis. BMC Genomics. 2018;19(1):363.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  25. 25.

    Xiao J, Zou TM, Chen YB, Chen L, Liu SJ, Tao M, et al. Coexistence of diploid, triploid and tetraploid crucian carp (Carassius auratus) in natural waters. BMC Genet. 2011;12(1):20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Hu F, Fan JJ, Qin QB, Huo YY, Wang YD, Wu C, et al. The sterility of allotriploid fish and fertility of female autotriploid fish. Front Genet. 2019;10:377.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Almon E, Goldfinger N, Kapon A, Schwartz D, Levine AJ, Rotter V. Testicular tissue-specific expression of the p53 suppressor gene. Dev Biol. 1993;156(1):107–16.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Harvey M, McArthur MJ, Montgomery CA, Bradley A, Donehower LA. Genetic background alters the spectrum of tumors that develop in p53-deficient mice. FASEB J. 1993;7:938–43.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Huang C, Liu W, Ji GX, Gu AH, Qu JH, Song L, et al. Genetic variants in TP53 and MDM2 associated with male infertility in Chinese population. Asian J Androl. 2012;14(5):691–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Rotter V, Schwartz D, Almon E, Goldfinger N, Kapon A, Meshorer A, et al. Mice with reduced levels of p53 protein exhibit the testicular giant-cell degenerative syndrome. P Proc Natl Acad Sci U S A. 1993;90(19):9075–9.

    CAS  Article  Google Scholar 

  31. 31.

    de Morais MP, Curado RF, Silva KS., Moura K, Arruda JT. Male idiopathic infertility and the TP53 polymorphism in codon 72. Genet Mol Res. 2016;15(4):gmr15048882.

  32. 32.

    Ebrahim Abadi Z, Khademi Bami M, Golzadeh M, Kalantar SM, Sheikhha MH. The frequency of TP53 R72P and MDM2 309T > G polymorphisms in Iranian infertile men with spermatogenetic failure: a case-control study. Int J Reprod Biomed. 2018;16(8):491–6.

    PubMed  PubMed Central  Google Scholar 

  33. 33.

    Saucedo L, Rumpel R, Sobarzo C, Schreiner D, Brandes G, Lustig L, et al. Deficiency of fibroblast growth factor 2 (FGF-2) leads to abnormal spermatogenesis and altered sperm physiology. J Cell Physiol. 2018;233(12):9640–51.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Garbarino Azúa DJ, Saucedo L, Giordana S, Magri ML, Buffone MG, Neuspiller F, et al. Fibroblast growth factor 2 (FGF2) is present in human spermatozoa and is related with sperm motility. The use of recombinant FGF2 to improve motile sperm recovery. Andrology. 2017;5(5):990–8.

    PubMed  Article  CAS  Google Scholar 

  35. 35.

    Saucedo L, Buffa GN, Rosso M, Guillardoy T, Góngora A, Munuce MJ, et al. Fibroblast growth factor receptors (FGFRs) in human dperm: expression, functionality and involvement in motility regulation. PLoS One. 2015;10(5):e0127297.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Saucedo L, Sobarzo C, Brukman NG, Guidobaldi HA, Lustig L, Giojalas LC, et al. Involvement of fibroblast growth factor 2 (FGF2) and its receptors in the regulation of mouse sperm physiology. Reproduction. 2018;156(2):163–72.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Hartl M, Glasauer S, Valovka T, Breuker K, Hobmayer B, Bister K. Hydra myc2, a unique pre-bilaterian member of the myc gene family, is activated in cell proliferation and gametogenesis. Biol Open. 2014;3(5):397–407.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Godeau F, Persson H, Gray HE, Pardee AB. C-myc expression is dissociated from DNA synthesis and cell division in Xenopus oocyte and early embryonic development. EMBO J. 1986;5(13):3571–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Aughey GN, Grice SJ, Liu JL. The Interplay between Myc and CTP Synthase in Drosophila. PLoS Genet. 2016;12(2):e1005867.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Wang ZH, Liu Y, Chaitankar V, Pirooznia M, Xu H. Electron transport chain biogenesis activated by a JNK-insulin-Myc relay primes mitochondrial inheritance in Drosophila. Elife. 2019;8:e49309.

    PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Jiang Y, Han K, Chen S, Wang Y, Zhang Z. Isolation, characterization, and expression of proto-oncogene cMyc in large yellow croaker Larimichthys crocea. Fish Physiol Biochem. 2017;43(5):1443–61.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Phochanukul N, Russell S. No backbone but lots of Sox: Invertebrate Sox genes. Int J Biochem Cell Biol. 2010;42(3):453–64.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Bowles J, Schepers G, Koopman P. Phylogeny of the SOX family of developmental transcription factors based on sequence and structural indicators. Dev Biol. 2000;227(2):239–55.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Dong WY, Wang Y, Zhou ZS, Guo JY. Sox genes in Agasicles hygrophila (Coleoptera: Chrysomelidae) are involved in ovarian development and oogenesis. Arch Insect Biochem Physiol. 2020;105(1):e21721.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Xia X, Huo W, Wan R, Zhang L, Xia X, Chang Z. Molecular cloning and expression analysis of Sox3 during gonad and embryonic development in Misgurnus anguillicaudatus. Int J Dev Biol. 2017;61(8–9):565–70.

  46. 46.

    Xia X, Wan R, Huo W, Zhang L, Xia X, Chang Z. Molecular cloning and mRNA expression pattern of Sox4 in Paramisgurnus dabryanus. Gene Expr Patterns. 2017;25–26:109–17.

  47. 47.

    Campolo F, Gori M, Favaro R, Nicolis S, Pellegrini M, Botti F, et al. Essential role of Sox2 for the establishment and maintenance of the germ cell line. Stem Cells. 2013;31(7):1408–21.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Dumic-Cule I, Peric M, Kucko L, Grgurevic L, Pecina M, Vukicevic S. Bone morphogenetic proteins in fracture repair. Int Orthop. 2018;42(11):2619–26.

    PubMed  Article  Google Scholar 

  49. 49.

    Juengel JL, McNatty KP. The role of proteins of the transforming growth factor-beta superfamily in the intraovarian regulation of follicular development. Hum Reprod Update. 2005;11(2):143–60.

    CAS  PubMed  Article  Google Scholar 

  50. 50.

    Peng C, Clelland E, Tan Q. Potential role of bone morphogenetic protein-15 in zebrafish follicle development and oocyte maturation. Comparative biochemistry and physiology. Comp Biochem Physiol A Mol Integr Physiol. 2009;153(1):83–7.

    PubMed  Article  CAS  Google Scholar 

  51. 51.

    Lei X, Cui K, Li Z, Su J, Jiang J, Zhang H, et al. BMP-1 participates in the selection and dominance of buffalo follicles by regulating the proliferation and apoptosis of granulosa cells. Theriogenology. 2016;85(5):999–1012.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Otsuka F. Interaction of melatonin and BMP-6 in ovarian steroidogenesis. Vitam Horm. 2018;107:137–53.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Clement JH, Fettes P, Knöchel S, Lef J, Knöchel W. Bone morphogenetic protein 2 in the early development of Xenopus laevis. Mech Dev. 1995;52(2–3):357–70.

  54. 54.

    Park ES, Woods DC, Tilly JL. Bone morphogenetic protein 4 promotes mammalian oogonial stem cell differentiation via Smad1/5/8 signaling. Fertil Steril. 2013;100(5):1468–75.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  55. 55.

    da Cunha EV, Melo L, Sousa GB, Araújo VR, Vasconcelos GL, Silva A, et al. Effect of bone morphogenetic proteins 2 and 4 on survival and development of bovine secondary follicles cultured in vitro. Theriogenology. 2018;110:44–51.

    PubMed  Article  CAS  Google Scholar 

  56. 56.

    Heikinheimo M, Ermolaeva M, Bielinska M, Rahman NA, Narita N, Huhtaniemi IT, et al. Expression and hormonal regulation of transcription factors GATA-4 and GATA-6 in the mouse ovary. Endocrinology. 1997;138(8):3505–14.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Laitinen MP, Anttonen M, Ketola I, Wilson DB, Ritvos O, Butzow R, et al. Transcription factors GATA-4 and GATA-6 and a GATA family cofactor, FOG-2, are expressed in human ovary and sex cord-derived ovarian tumors. J Clin Endocrinol Metab. 2000;85(9):3476–83.

    CAS  PubMed  Google Scholar 

  58. 58.

    Lowry JA, Atchley WR. Molecular evolution of the GATA family of transcription factors: conservation within the DNA-binding domain. J Mol Evol. 2000;50(2):103–15.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    Bennett-Toomey J, Stocco C. GATA regulation and function during the ovarian life cycle. Vitam Horm. 2018;107:193–225.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Efimenko E, Padua MB, Manuylov NL, Fox SC, Morse DA, Tevosian SG. The transcription factor GATA4 is required for follicular development and normal ovarian function. Dev Biol. 2013;381(1):144–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Maehama T, Dixon JE. The tumor suppressor, PTEN/MMAC1, dephosphorylates the lipid second messenger, phosphatidylinositol 3,4,5-trisphosphate. J Biol Chem. 1998;273(22):13375–8.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Adhikari D, Liu K. Molecular mechanisms underlying the activation of mammalian primordial follicles. Endocr Rev. 2009;30(5):438–64.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Reddy P, Liu L, Adhikari D, Jagarlamudi K, Rajareddy S, Shen Y, et al. Oocyte-specific deletion of Pten causes premature activation of the primordial follicle pool. Science. 2008;319(5863):611–3.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Ernst EH, Grøndahl M, Grund S, Hardy K, Heuck A, Sunde L, et al. Dormancy and activation of human oocytes from primordial and primary follicles: molecular clues to oocyte regulation. Hum Reprod. 2017;32(8):1684–700.

    CAS  PubMed  Article  Google Scholar 

  65. 65.

    Jouaux A, Franco A, Heude-Berthelin C, Sourdaine P, Blin JL, Mathieu M, et al. Identification of Ras, Pten and p70S6K homologs in the Pacific oyster Crassostrea gigas and diet control of insulin pathway. Gen Comp Endocrinol. 2012;176(1):28–38.

    CAS  PubMed  Article  Google Scholar 

  66. 66.

    Mensah LB, Goberdhan D, Wilson C. mTORC1 signalling mediates PI3K-dependent large lipid droplet accumulation in Drosophila ovarian nurse cells. Biol Open. 2017;6(5):563–70.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Qin QB, Wang J, Wang YD, Liu Y, Liu SJ. Organization and variation analysis of 5S rDNA in gynogenetic offspring of Carassius auratus red var. (♀) × Megalobrama amblycephala (♂). BMC Genetics. 2015;16:26.

    PubMed  PubMed Central  Article  Google Scholar 

  68. 68.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references


Not applicable.


This work was supported by the Natural Science Foundation for Distinguished Young Scholars of Hunan Province (Grant No. 2017JJ1022), the National Natural Science Foundation of China (Grant No. 31730098, U19A2040), the Earmarked Fund for China Agriculture Research System (Grant No. CARS-45), Hunan Provincial Science and Technology Department (2019RS5001) and the National Key Research and Development Program of China (2020YFD0900104).

Author information




SL and QQ have designed of the work. CW has contributed to this study for the design, in executing experiments and in writing manuscript. XL, HQ and CZ have made substantial contributions to the acquisition and analysis of data. LY, TY, YZ, XH and XX have substantively revised the work. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Qinbo Qin or Shaojun Liu.

Ethics declarations

Ethics approval and consent to participate

The study was approved by Ethics Committee of Hunan Normal University, all methods were carried out in accordance with relevant guidelines and regulations. This study was carried out in compliance with the ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that there are no competing financial interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Quality test results of RNA.

Additional file 2: Table S2.

Transcript classification based on gene ontology (GO) for DETs in F3nRR vs FRCC.

Additional file 3: Table S3.

Transcript classification based on gene ontology (GO) for DETs in M3nRR vs M4nRR.

Additional file 4: Table S4.

Transcript classification based on Kyoto Encyclopedia of Genes and Genomes (KEGG) for DETs in F3nRR vs FRCC.

Additional file 5: Table S5.

Transcript classification based on Kyoto Encyclopedia of Genes and Genomes (KEGG) for DETs in M3nRR vs M4nRR.

Additional file 6: Table S6.

Female fertility-related genes.

Additional file 7: Table S7.

Male fertility-related genes.

Additional file 8: Table S8.

Table of F3nRR vs FRCC of fertility-related gene in protein-protein interaction network.

Additional file 9: Table S9.

Table of M3nRR vs M4nRR of fertility-related gene in protein-protein interaction network.

Additional file 10: Table S10.

Sequences of primers used in this study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, C., Luo, X., Qin, H. et al. Formation of autotriploid Carassius auratus and its fertility-related genes analysis. BMC Genomics 22, 435 (2021).

Download citation


  • Fertility
  • Autotriploid fish
  • Transcriptome
  • Gonad