- Open Access
Formation of autotriploid Carassius auratus and its fertility-related genes analysis
BMC Genomics volume 22, Article number: 435 (2021)
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.
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 . 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 . 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 . In addition, studies explored biological mechanisms that occur mid-oogenesis , and regulation of late oogenesis . Spermatogenesis is divided into three steps: spermatogonial mitotic proliferation, two times of meiosis, and post-meiotic differentiation . Previous studies explored the functional mechanisms of spermatocytogenesis , meiosis during spermatogenesis , and spermiogenesis .
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 . 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 . However, a different study reported that triploid fish can produce normal gametes . 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 . 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 . 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 . 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 . Studies on Takifugu rubripes reported that sex-related genes play an important role at early sex differentiation stage . Gonadal transcriptome profiling of triploid hybrid loaches (Misgurnus anguillicaudatus) and their diploid and tetraploid parents showed key genes implicated in low hybrid triploid fertility . A study on Thunnus maccoyii reported sex and gonad-development-related genes in the gonads of Southern bluefin tuna through RNA-Seq . In addition, RNA-Seq has been successfully used to analyze sex determination and differentiation related genes in tilapia . 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.
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.
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.
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 %.
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).
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).
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.
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).
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.
Triploid animals are usually sterile and cannot form triploid populations. However, previous studies have been reporting contradicting results. Xiao et al.  reported that triploid Carassius auratus in Dongting water system produces normal gametes. Hu et al.  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 . 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 . In addition, TP53 codon 72 polymorphism in mice is involved in meiosis, implying that it plays a critical role in spermatogenesis . 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 . 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 . 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 . In mouse, knock out of FGF2 induces impaired sperm production and is associated with alterations in sperm morphology and function . 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 . In Xenopus, C-MYC was detected in oocytes, indicating that it plays a role in oogenesis . 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 . 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 . 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 . 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 . Sox genes have been explored in many organisms and can be classified into ten subgroups (A-J) . 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 . In Misgurnus anguillicaudatus, MaSOX3 is abundant in primary oocytes and previtellogenic oocyte cells, indicating that MaSox3 gene is involved in ovarian development . In Paramisgurnus dabryanus, SOX4 was detected in the ovary, showing that it plays an important role during ovarian development . In mouse, expression of SOX2 is required for establishment and maintenance of the oocyte cell . 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 . A previous study reports that BMPs play a role in regulation of ovarian follicular development . 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 . In mouse, BMP4 regulates the number of oocytes, suggesting its role in the process of oogenesis . Roles of BMP2 and BMP4 as important factors in survival and development of bovine secondary follicles were recently reported . 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 . GATA4 and GATA6 knockout female mice exhibited infertility due to disrupted formation of ovaries . GATA4 deletion resulted in a sterile female mice phenotype attributed to drastic reduction in number of developing follicles . The mRNA for GATA4 has been reported in human ovary implying that GATA4 plays a role in ovarian folliculogenesis . 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 . PTEN is a negative regulator of PI3K-Akt signaling pathway which is involved in growth of eggs . Deletion of PTEN from oocytes affects mouse fertility by interrupting oocyte growth . PTEN signaling pathway associated with ovarian follicle development has been reported in human . In Crassostrea gigas, PTEN is involved in insulin pathway in gonads and plays a critical role in reproduction . In Drosophila, loss of PTEN is related to IIS/mTORC1 signalling, which is important for oogenesis . 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 . 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 (https://bigd.big.ac.cn/search?dbId=gwh&q=GWHAAIA00000000) 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 (https://string-db.org/) 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 .
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) (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA694292) and the expression profiles of RNA-Seq data are included in Figshare (https://doi.org/10.6084/m9.figshare.14561670.v1).
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
- 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
Kyoto Encyclopedia of Genes and Genomes
- TP53 :
Tumor protein p53
- FGF2 :
Fibroblast growth factor 2
- MYC :
MYC proto-oncogene, bHLH transcription factor
- SOX2 :
SRY-box transcription factor 2
- BMP4 :
Bone morphogenetic protein 4
- GATA4 :
GATA binding protein 4
- PTEN :
Phosphatase and tensin homolog
- BMP2 :
Bone morphogenetic protein 2
Otto SP, Whitton J. Polyploid incidence and evolution. Annu Rev Genet. 2000;34(1):401–37.
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.
Nagahama Y. Molecular mechanisms of sex determination and gonadal sex differentiation in fish. Fish Physiol Biochem. 2005;31(2–3):105–9.
Quinn A, Koopman P. The molecular genetics of sex determination and sex reversal in mammals. Semin Reprod Med. 2012;30(5):351–63.
Windley SP, Wilhelm D. Signaling pathways involved in mammalian sex determination and gonad development. Sex Dev. 2015;9(6):297–315.
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.
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.
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.
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.
Kawaguchi S, Ueki M, Kai T. Drosophila MARF1 ensures proper oocyte maturation by regulating nanos expression. PLoS One. 2020;15(4):e0231114.
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.
Phillips BT, Gassei K, Orwig KE. Spermatogonial stem cell regulation and spermatogenesis. Phil Trans R Soc B. 2010;365:1663–78.
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.
Bettegowda A, Wilkinson MF. Transcription and post-transcriptional regulation of spermatogenesis. Phil Trans R Soc B. 2010;365:1637–51.
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.
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.
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.
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.
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.
Qian X, Ba Y, Zhuang QF, Zhong GF. RNA-Seq technology and its application in fish transcriptomics. OMICS. 2014;18(2):98–110.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Aughey GN, Grice SJ, Liu JL. The Interplay between Myc and CTP Synthase in Drosophila. PLoS Genet. 2016;12(2):e1005867.
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.
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.
Phochanukul N, Russell S. No backbone but lots of Sox: Invertebrate Sox genes. Int J Biochem Cell Biol. 2010;42(3):453–64.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Otsuka F. Interaction of melatonin and BMP-6 in ovarian steroidogenesis. Vitam Horm. 2018;107:137–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.
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.
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.
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.
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.
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.
Bennett-Toomey J, Stocco C. GATA regulation and function during the ovarian life cycle. Vitam Horm. 2018;107:193–225.
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.
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.
Adhikari D, Liu K. Molecular mechanisms underlying the activation of mammalian primordial follicles. Endocr Rev. 2009;30(5):438–64.
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.
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.
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.
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.
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.
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.
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).
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
The authors declare that there are no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Quality test results of RNA.
Transcript classification based on gene ontology (GO) for DETs in F3nRR vs FRCC.
Transcript classification based on gene ontology (GO) for DETs in M3nRR vs M4nRR.
Transcript classification based on Kyoto Encyclopedia of Genes and Genomes (KEGG) for DETs in F3nRR vs FRCC.
Transcript classification based on Kyoto Encyclopedia of Genes and Genomes (KEGG) for DETs in M3nRR vs M4nRR.
Female fertility-related genes.
Male fertility-related genes.
Table of F3nRR vs FRCC of fertility-related gene in protein-protein interaction network.
Table of M3nRR vs M4nRR of fertility-related gene in protein-protein interaction network.
Sequences of primers used in this study.
About this article
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). https://doi.org/10.1186/s12864-021-07753-5
- Autotriploid fish