Skip to main content

Comparative analyses of the Sox9a-Amh-Cyp19a1a regulatory Cascade in Autotetraploid fish and its diploid parent



Autotetraploid Carassius auratus (4nRCC, 4n = 200, RRRR) was derived from the whole genome duplication of diploid red crucian carp (Carassius auratus red var.) (2nRCC, 2n = 100, RR). To investigate the genetic effects of tetraploidization, we analyzed DNA variation, epigenetic modification and gene expression changes in the Sox9a-Amh-Cyp19a1a regulatory cascade between 4nRCC and 2nRCC.


We found that the Sox9a gene contained two variants in 2nRCC and four variants in 4nRCC. Compared with that in 2nRCC, DNA methylation in the promoter regions of the Amh and Cyp19a1a genes in 4nRCC was altered by single nucleotide polymorphism (SNP) mutations, which resulted in the insertions and deletions of CpG sites, and the methylation levels of the Sox9a, Amh and Cyp19a1a genes increased after tetraploidization. The gene expression level of the Sox9a-Amh-Cyp19a1a regulatory cascade was downregulated in 4nRCC compared with that in 2nRCC.


The above results demonstrate that tetraploidization leads to significant changes in the genome, epigenetic modification and gene expression in the Sox9a-Amh-Cyp19a1a regulatory cascade; these findings increase the extant knowledge regarding the effects of polyploidization.


Fertile polyploids play important roles in promoting the exchange of genetic material among species, enriching species diversity, and laying foundations for polyploid breeding [1,2,3,4]. Polyploidy can be classified into autopolyploidy and allopolyploidy [5, 6]. The former can be produced by the whole genome duplication (WGD) of a species (autopolyploidization), whereas the latter refers to the merging of genomes from different species after hybridization (allopolyploidization) [7]. Autopolyploidization is of great significance in the origin and evolution of vertebrates; however, it is rare in teleost fishes [8,9,10,11,12]. In previous studies, fertile allotetraploids (4nRB, 4n = 148) were successfully obtained in the first generation derived from the distant hybridization of Carassius auratus red var. (2nRCC, 2n = 100, female) × Megalobrama amblycephala (BSB, 2n = 48, male) [13]. Fertile autotetraploids (4nRCC, 4n = 200) were obtained in the second generation through crossing diploid sperm and diploid eggs produced by special chromosome behavior in germ cells of 4nRB [14]. The 4nRCC possessed four sets of chromosomes derived from 2nRCC, which produced diploid gametes [14, 15]. We successfully established the autotetraploid lineage (F2-F13) through continuous self-crossing. Phenotypic changes occurred both in 4nRB and 4nRCC, including the presence of the gray body color and the barbel, which were absent in 2nRCC [14]. Furthermore, the morphological traits of 4nRCC obviously differed from those of 4nRB (e.g. lateral scales, upper lateral scales, dorsal fins, and anal fins) [14]. Hence, this autotetraploid fish lineage provides not only abundant diploid gamete resources for polyploid genetic breeding but also an excellent model to study the consequences of inheritance and evolution in polyploid vertebrates.

The Sox9-Amh-Cyp19a1 regulatory cascade was identified in 2013, which played a role in sex differentiation and maintenance in Atlantic cod [16]. This regulatory cascade had a similar role in the gonad development of zebrafish [17, 18]. In this regulatory cascade, Sox9 (SRY-box containing gene 9) activated the expression of Amh (anti-Müllerian hormone), which induced the degeneration of the Müllerian ducts and inhibited the expression of Cyp19a1 (cytochrome P450 family 19 subfamily A member 1) [19,20,21]. Testicular expression of Sox9a in adult Atlantic cod and zebrafish was consistent with male-specific expression of Sox9 in mammal [16,17,18]. Amh had functions in both sexes in fish, especially in males, in which it was the basic factor underlying sexual differentiation [22]. In addition, P450 aromatase (Cyp19a1a), a key steroidogenic enzyme that converted androgen into estrogen, was widely accepted to be actively involved in ovarian formation, development, and maintenance in fish [23, 24]. Although the Sox9a-Amh-Cyp19ala regulatory cascade played a crucial role in the complex process of sex differentiation in some vertebrates, studies of conserved regulatory elements in 2nRCC and 4nRCC could be served as the starting point for a broader understanding of the impact of tetraploidization processes.

The autotetraploid fish lineage provides good experimental material for the study of genomic DNA variation in polyploids, which often exhibited changes in gene structure, gene expression and DNA methylation [25,26,27,28,29,30]. To investigate the genetic effects of tetraploidization, we analyzed DNA variation, epigenetic modification and gene expression changes in the Sox9a-Amh-Cyp19a1a regulatory cascade between 4nRCC and 2nRCC. This is the first report concerning the Sox9a-Amh-Cyp19a1a regulatory cascade in autotetraploid fish, and its results provide a broad survey of the effects of polyploidization.


Cloning and homology analysis of CDS regions

By splicing the intermediate fragments and 3′ RACE sequences of Sox9a, Amh, and Cyp19a1a, we obtained complete CDS regions. For the Sox9a gene, two different forms in 2nRCC (GenBank Accession Nos. MK307791 and MK307773) and four different forms in 4nRCC (GenBank Accession Nos. MK307774-MK307777) with high homology were observed (namely, variants 1 and 2 in 2nRCC, variants 1 to 4 in 4nRCC); the sizes of the CDS region of these genes were 1398, 1329, 1398, 1344, 1383 and 1329 bp, respectively. For the Amh gene, the complete CDS lengths in 2nRCC and 4nRCC were 1707 bp and 1695 bp (GenBank Accession Nos. MK224515 and MK224514), respectively. The complete CDS length of Cyp19a1a in both 2nRCC and 4nRCC was 1557 bp (GenBank Accession Nos. MK224513 and MK224512). Homology analysis showed that 2nRCC and 4nRCC had high homology at the amino acid sequence level, with amino acid identities of 90% for Amh and 98% for Cyp19a1a. As shown in Fig. 1, CDS alterations of Sox9a involving SNP insertions, deletions and mutations were observed. Phylogenetic analysis was used to determine the relationships among the Sox9a variant sequences in 2nRCC and 4nRCC using the neighbor-joining method in the MEGA7 software (Fig. 2).

Fig. 1

Alignment of the Sox9a variant sequences of 2nRCC and 4nRCC. Deletions and insertions in all sequences are framed by black boxes. Dots represent sequence identity

Fig. 2

Phylogenetic analysis demonstrating the relationships among the Sox9a variant sequences of 2nRCC and 4nRCC using the neighbor-joining method in MEGA7

Expression of Sox9a, Amh, and Cyp19a1a in mature gonads

To investigate the expression of the target genes in mature gonads, qPCR was used to compare expression patterns between 2nRCC and 4nRCC. The results showed that the expression of Sox9a and Amh was predominant in the mature testes (Fig. 3a and b; P < 0.01), whereas the expression level of Cyp19a1a was higher in the mature ovaries (Fig. 3c; P < 0.01). Compared with 2nRCC, 4nRCC showed downregulated expression levels of Sox9a in the mature testes and ovaries (Fig. 3a; P < 0.01), and the expression pattern of Amh was observed to be same as that of Sox9a (Fig. 3b; P < 0.01). As shown in Fig. 3c, in the mature ovaries, 2nRCC showed high expression of Cyp19a1a compared with that of 4nRCC (P < 0.01), but there was no significant difference in the mature testes (P > 0.05). These results showed that the expression of the target genes in the Sox9a-Amh-Cyp19a1a regulatory cascade was down-regulation following polyploidization.

Fig. 3

Relative quantification of Sox9a (a), Amh (b), and Cyp19a1a (c) in 2nRCC and 4nRCC. Both *** and **** indicate significant differences between 2nRCC and 4nRCC (P < 0.01); ns indicates no significant difference between 2nRCC and 4nRCC

DNA methylation status of the Sox9a, Amh, and Cyp19a1a genes

Using bisulfite sequencing and BiQ Analyzer software, we evaluated the differential DNA methylation status of target gene promoters in the Sox9a-Amh-Cyp19a1a cascade in the mature gonads of 2nRCC and 4nRCC. Among the 3 CpGs found between − 1120 bp and − 764 bp of the Sox9a promoter sequence (Fig. 4a), the average methylation levels of 4nRCC were higher than those of 2nRCC in the mature testes and ovaries (Fig. 5). Numbers with a plus or minus sign show CpG positions relative to the transcription start site ( Compared with the 2nRCC sequence, the 4nRCC Amh gene had deletions of three CpG sites, and the Cyp19a1a gene had three deleted CpG sites and two newly inserted CpG sites (Fig. 5). In 2nRCC and 4nRCC, 7 CpGs at − 1363 bp and − 915 bp, 4 CpGs at − 1372 bp and − 915 bp were found in the Amh promoter sequence, respectively (Fig. 4b). In addition, 10 CpGs at − 359 bp and + 36 bp, 9 CpGs at − 365 bp and + 36 bp were found in the Cyp19a1a promoter sequences of 2nRCC and 4nRCC, respectively (Fig. 4c). These results revealed that SNP mutations resulted in the insertions and deletions of CpG sites. We demonstrated the average methylation levels (Table 1) and differential CpG methylation status (Fig. 5) of the Sox9a, Amh and Cyp19a1a promoters in the mature testes and ovaries of 2nRCC and 4nRCC. As shown in Fig. 6, the methylation levels of the target gene promoters showed an obvious negative correlation with gene expression levels.

Fig. 4

Amplified regions of the predicted promoter sequences relative to the transcription start site. Numbers with a plus or minus sign show CpG positions relative to the transcription start site. The yellow boxes represent CpG dinucleotide sites in the target gene promoter. a Three CpG sites in the Sox9a gene promoters of 2nRCC and 4nRCC are marked with yellow boxes. b Seven CpG sites in the Amh gene promoter of 2nRCC are marked with yellow boxes, whereas only four CpG sites are present in the 4nRCC promoter. c Ten CpG sites in the Cyp19a1a gene promoter of 2nRCC are marked with yellow boxes, and nine CpG sites in the promoter of 4nRCC are marked. The last CpG site in both 4nRCC and 2nRCC was not included due to the primer position

Fig. 5

Differential CpG methylation status of Sox9a (a), Amh (b) and Cyp19a1a (c) promoters in mature testes and ovaries of 2nRCC and 4nRCC. Each box corresponds to one CpG position in the genomic sequence. Ten clones per sample were analyzed. Yellow boxes indicate methylation status, blue boxes indicate unmethylated status, and gray boxes indicate the non-present CpG positions

Table 1 Average methylation levels of all CpG sites in 2nRCC and 4nRCC
Fig. 6

Correlation between the expression of the Sox9a (a), Amh (b) and Cyp19a1a (c) genes and CpG methylation in 2nRCC and 4nRCC


After genome duplication, polyploids were thought to display changes in their genome structure, gene expression, and epigenetic modification [1, 28, 31]. The structural changes to the genome in polyploids consisted of deletions, insertions, duplications, translocations and transpositions [32]. Some duplicated genes might undergo pseudogenization, neofunctionalization or subfunctionalization during polyploidization [16, 33]. In previous studies, obvious DNA variation, pseudogenization and sequence divergence were observed in allotetraploids and autotetraploids [33, 34]. Although 4nRCC was originated from 2nRCC, the CDS regions and amino acid sequences of Sox9a gene presented large differences, which provided direct evidence of genome duplication and alterations following tetraploidization. First, the newly synthetic polyploidy required to solve the problem of survival and reproductive fitness, meaning they needed to overcome the genomic incompatibility and transcriptome shock by genomic changes. Importantly, a variety of genomic alterations at structural level (e.g. extensive intra- and interchromosomal rearrangements ensuring normal chromosome pairing in meiosis) had been shown to re-establish the genome balance after WGDs [10]. Thus, we speculated that the genomic alteration was used to adapt the tetraploidization, which not only increased genomic flexibility and the rate of gene differentiation but also had great advantages for long-term adaptation and persistence of autopolyploids.

Genetic redundancy and intergenomic interactions were among the general features of polyploids that induced genetic and epigenetic changes (e.g. epigenetic remodeling to silence parts of the duplicated genome), leading to recombination of the genome and regulatory networks [10, 35]. Compared with that of 2nRCC, the methylation level of the Sox9a-Amh-Cyp19a1a regulatory cascade increased in 4nRCC after genome duplication, and a similar result was found previously in the autopolyploid Chrysanthemum lavandulifolium [6]. Meanwhile, the genomic alterations in the promoter region of target genes led to insertions and deletions of CpG sites after an autopolyploidization event, producing an important effect on the genomic methylation status. It is possible that genetic and epigenetic alterations may be related to each other after a polyploidization event. Further, epigenetic changes may regulate the level of gene expression, which had the potential to add the stability of the genome after polyploidization. Notably, the same sex-specific methylation patterns were retained in 4nRCC and 2nRCC. In the previous study, DNA methylation may be a mechanism used for natural sex maintenance, as in the Chinese sea perch Lateolabrax maculatus, and a potential regulatory mechanism during temperature-induced sex change, as in Nile tilapia and European sea bass (Dicentrarchus labrax) [36,37,38]. Although the Sox9a, Amh and Cyp19a1a genes presented significant sex correlations in fishes, further experiments are required to verify whether the Sox9a-Amh-Cyp19a1a regulatory cascade is related to sex determination or sex differentiation in 4nRCC and 2nRCC. The results so far suggest that 2nRCC and 4nRCC may have similar sex differentiation mechanisms.

Similarly, to re-establish the transcriptional balance, the changes of gene expression patterns (increase or decrease) were required across the genome [10]. Previous studies of gene expression changes in polyploids had highlighted multiple potential mechanisms, including homologous gene silencing, altered regulatory networks, DNA methylation, and chromatin remodeling [26, 27, 39]. Because allopolyploidy encompasses different genomes and increases genome complexity, autopolyploids are thought to be the best material to evaluate the possibility of ploidy-dependent regulatory changes [26]. Here, the expression levels of the Sox9a, Amh, and Cyp19a1a genes were significantly downregulated (P < 0.01) after tetraploidization, which suggested that tetraploidization has the potential to lead to changes in gene expression through the combined effects of genomic variation and DNA methylation. Autopolyploids were generally expected to maintain genomes highly similar to those of their diploid parents and not exhibited significant changes in gene expression [28, 40], whereas our results clearly exhibited such changes. These changes in gene expression may be very important to help polyploids adapt to the consequences of WGD and maintain genome balance.


Our data demonstrated that tetraploidization was associated with significant changes in genome sequences, epigenetic modifications and gene expression by surveying the Sox9a-Amh-Cyp19a1a regulatory cascade. This study provides important evidence on genomic DNA variation and new information about the effects of polyploidization.



Both 2nRCC and 4nRCC used in this study were cultivated in ponds with expanded feed at the Protection Station of Polyploid Fish of Hunan Normal University, Hunan, China. During the reproductive season (April–June), both 2nRCC and 4nRCC reached sexual maturity at 1 year-old and produced a large number of mature eggs and white sperm normally, respectively. All fish used as samples were anesthetized with 100 mg/L MS-222 (Sigma-Aldrich, St Louis, MO, USA) prior to dissection. Gonads were excised from three adult males and three adult females under sterile conditions, frozen quickly in liquid nitrogen, and stored at − 80 °C for further use.

RNA isolation and RT-PCR

Total RNA was isolated from the mature gonad tissues of 2nRCC and 4nRCC using RNAiso reagent (TaKaRa, Japan) following the manufacturer’s protocol. The RNA purity and concentration were measured using Synergy 2 Multi-Mode Microplate Reader (BioTek,, and their integrity were evaluated using 1% agarose gel electrophoresis [34]. The first-strand cDNA was synthesized using the Maxima H Minus First Strand cDNA Synthesis Kit with dsDNase (Thermo Scientific, USA) in a 20 μl reaction volume. The synthesized cDNA was stored at − 20 °C for further use.

Coding sequence (CDS) cloning of Sox9a, Amh and Cyp19a1a

Based on the cDNA sequences of the Carassius auratus × Cyprinus carpio × Carassius cuvieri Sox9a gene (DQ201318), the goldfish Amh gene (KF640083) and the Carassius auratus Cyp19a1a gene (KC147009) available at NCBI, special primers for each gene were designed using Primer Premier 5 software to amplify the intermediate fragments. Subsequently, 3′ RACE was performed with forward primers designed from each intermediate fragment and a reverse adaptor primer. The 3′ RACE products of the target genes from 2nRCC and 4nRCC were sequenced, and their complete CDS regions were obtained after splicing the partial CDS regions. All primers were presented in Table 2, and the PCR products were cloned as described in [41] and sequenced by Tsingke (Beijing, China).

Table 2 Nucleotide sequences of primers used for CDS cloning, 3′ RACE and RT-qPCR in this study

Quantitative real-time PCR analysis

The cDNAs prepared from the gonads of 2nRCC and 4nRCC were used to detect relative gene expression using quantitative real-time PCR (qPCR). Species-specific primers for a housekeeping gene (β-actin) and three target genes (Sox9a, Amh and Cyp19a1a) were designed with Primer Premier 5 software and presented in Table 2. The cDNAs from the adult testes and ovaries of 2nRCC and 4nRCC were used as templates. All of the experiments in this stage were performed three samples and three replicas to improve the accuracy of the results. The qPCR amplification, using ABI Prism 7500 Sequence Detection System (Applied Biosystems, USA), was performed in a 10 μl reaction mixture containing 5 μl of SYBR Green PCR Master Mix (ABI), 1 μl of cDNAs (in a dilution of 1:15), 0.4 μl of each primer (10 μM) and 3.2 μl of sterilized water. The amplification conditions were 50 °C for 5 min, 95 °C for 10 min, and 40 cycles at 95 °C for 15 s and 60 °C for 45 s. After that, the melting curve was performed at the end of the assay to verify the generation of a single product. The dissociation conditions were 95 °C for 1 min, 60 °C for 30 s and 95 °C for 30 s. The relative expression was determined with GraphPad Prism 6 software using the 2-ΔΔCT method [42].

Genomic DNA isolation and DNA bisulfite modification

Genomic DNA was extracted from all samples of adult testes and ovaries using the TaKaRa MiniBEST Universal Genomic DNA Extraction Kit Ver. 5.0 (TaKaRa, Japan). The methods for detecting DNA purity and concentration were similar to those used after RNA isolation. Then, 400 ng of each sample was sodium bisulfite-modified using the EZ DNA Methylation-Gold™ Kit (Zymo Research, USA) according to the manufacturer’s protocol. The bisulfite-treated DNA was stored at − 20 °C.

Amplifying promoter regions and predicting CpG-rich regions

The first exon of each of the three target genes was aligned to the genome of 2nRCC (NCBI project accession no. PRJNA289059) using Blast (V2.5.0) software, and the promoter positions were considered to be within 2 kb upstream of the first exon. The promoters of Sox9a, Amh and Cyp19a1a were amplified using untreated genomic DNA as the template, and the primers were listed in Table 3. The CpG-rich regions of the Sox9a, Amh and Cyp19a1a promoters and exons were predicted by the online MethPrimer design software (

Table 3 Nucleotide sequences of primers used for promoter amplification and BS-PCR in this study

Bisulfite PCR (BS-PCR) and analysis

BS-PCR primers were designed with Primer Premier 5 software based on the sense strand of the bisulfite-modified DNA (Table 3). PCR was performed in a final volume of 50 μl using LA Taq™ (TaKaRa). The amplification conditions were as follows: 4 min at 94 °C; 35 cycles of 30 s at 94 °C, 30 s at 56 °C for Sox9a/ Cyp19a1a or 54 °C for Amh, and 30 s at 72 °C; and 5 min at 72 °C for extension [43]. All PCR products were detected and separated on a 2% agarose gel, purified using the Gel Extraction Kit (Tsingke), cloned into the pMD18-T vector (TaKaRa) and transferred into E. coli DH5α (Sangon, China). To investigate the status of each methylation site, at least 10 positive clones obtained after screening by PCR amplification were sequenced by Tsingke (Beijing, China).

Statistical analysis

Analyses of variance and pairwise comparisons of the data were analyzed by SPSS 17.0 software.

Availability of data and materials

The datasets supporting the conclusions of this article were available in the GenBank repository with access No. MK307791, MK307773-MK307777, and MK224512-MK224515.



Autotetraploid Carassius auratus


Diploid red crucian carp (Carassius auratus red var.)


Single nucleotide polymorphism


Coding sequence


Rapid amplification of cDNA ends


Real-time quantitative polymerase chain reaction


Bisulfite polymerase chain reaction


Megalobrama amblycephala


Allotetraploid Carassius auratus


  1. 1.

    Song C, Liu S, Xiao J, He W, Zhou Y, Qin Q, Zhang C, Liu Y. Polyploid organisms. Sci China Life Sci. 2012;55(4):301–11.

    Article  Google Scholar 

  2. 2.

    Liu S. Fish distant hybridization: China social sciences press; 2014.

    Google Scholar 

  3. 3.

    Chen J, Luo M, Li S, Tao M, Ye X, Duan W, Zhang C, Qin Q, Xiao J, Liu S. A comparative study of distant hybridization in plants and animals[J]. Science China Life Sciences, 2018;61(3):285–309.

  4. 4.

    Wang S, Tang C, Tao M, Qin Q, Zhang C, Luo K, Zhao R, Wang J, Ren L, Xiao J, et al. Establishment and application of distant hybridization technology in fish. Sci China Life Sci. 2018;62(1):22–45.

    Article  Google Scholar 

  5. 5.

    Qin Q, He W, Liu S, Wang J, Xiao J, Liu Y. Analysis of 5S rDNA organization and variation in polyploid hybrids from crosses of different fish subfamilies. J Exp Zool B Mol Dev Evol. 2010;314(5):403–11.

    Article  Google Scholar 

  6. 6.

    Gao R, Wang H, Dong B, Yang X, Chen S, Jiang J, Zhang Z, Liu C, Zhao N, Chen F. Morphological, Genome and Gene Expression Changes in Newly Induced Autopolyploid Chrysanthemum lavandulifolium (Fisch. ex Trautv.) Makino. Int J Mol Sci. 2016;17(10).

  7. 7.

    Julca I, Marcethouben M, Vargas P, Gabaldón T. Phylogenomics of the olive tree ( Olea europaea ) reveals the relative contribution of ancient Allo- and autopolyploidization events. BMC Biol. 2018;16(1):15.

    Article  Google Scholar 

  8. 8.

    Jenczewski E, Alix K. From diploids to allopolyploids: the emergence of efficient pairing control genes in plants. Crit Rev Plant Sci. 2004;23(1):21–45.

  9. 9.

    Parisod C, Holderegger R, Brochmann C. Evolutionary consequences of autopolyploidy. New Phytol. 2010;186(1):5–17.

    CAS  Article  Google Scholar 

  10. 10.

    Wertheim B, Beukeboom LW, van de Zande L. Polyploidy in animals: effects of gene expression on sex determination, evolution and ecology. Cytogenet Genome Res. 2013;140(2–4):256–69.

    CAS  Article  Google Scholar 

  11. 11.

    Liu S, Luo J, Chai J, Ren L, Zhou Y, Huang F, Liu X, Chen Y, Zhang C, Tao M, et al. Genomic incompatibilities in the diploid and tetraploid offspring of the goldfish x common carp cross. Proc Natl Acad Sci U S A. 2016;113(5):1327–32.

    CAS  Article  Google Scholar 

  12. 12.

    Qin Q, Lai Z, Cao L, Xiao Q, Wang Y, Liu S. Rapid genomic changes in allopolyploids of Carassius auratus red var. (♀) × Megalobrama amblycephala (♂)[J]. Sci Rep. 2016;6:34417.

  13. 13.

    Liu S, Qin Q, Xiao J, Lu W, Shen J, Li W, Liu J, Duan W, Zhang C, Tao M, et al. The formation of the polyploid hybrids from different subfamily fish crossings and its evolutionary significance. Genetics. 2007;176(2):1023–34.

    CAS  Article  Google Scholar 

  14. 14.

    Qin Q, Wang Y, Wang J, Dai J, Xiao J, Hu F, Luo K, Tao M, Zhang C, Liu Y, et al. The autotetraploid fish derived from hybridization of Carassius auratus red var. (female) x Megalobrama amblycephala (male). Biol Reprod. 2014;91(4):93.

    Article  Google Scholar 

  15. 15.

    Qin Q, Huo Y, Liu Q, Wang C, Zhou Y, Liu S. Induced gynogenesis in autotetraploids derived from Carassius auratus red var. (♀) × Megalobrama amblycephala (♂). Aquaculture. 2018;495:710–4.

    Article  Google Scholar 

  16. 16.

    Johnsen H, Tveiten H, Torgersen JS, Andersen O. Divergent and sex-dimorphic expression of the paralogs of the Sox9-Amh-Cyp19a1 regulatory cascade in developing and adult Atlantic cod (Gadus morhua L.). Mol Reprod Dev. 2013;80(5):358–70.

    CAS  Article  Google Scholar 

  17. 17.

    Rodriguez-Mari A, Yan YL, Bremiller RA, Wilson C, Canestro C, Postlethwait JH. Characterization and expression pattern of zebrafish anti-Mullerian hormone (Amh) relative to sox9a, sox9b, and cyp19a1a, during gonad development. Gene Expr Patterns. 2005;5(5):655–67.

    CAS  Article  Google Scholar 

  18. 18.

    Sun D, Zhang Y, Wang C, Hua X, Zhang XA, Yan J. Sox9-related signaling controls zebrafish juvenile ovary-testis transformation. Cell Death Dis. 2013;4:e930.

    CAS  Article  Google Scholar 

  19. 19.

    Cutting A, Chue J, Smith CA. Just how conserved is vertebrate sex determination? Dev Dyn. 2013;242(4):380–7.

    CAS  Article  Google Scholar 

  20. 20.

    Nagabhushana A, Mishra RK. Finding clues to the riddle of sex determination in zebrafish. J Biosci. 2016;41(1):145–55.

    CAS  Article  Google Scholar 

  21. 21.

    Yang Y-J, Wang Y, Li Z, Zhou L, Gui J-F. Sequential, divergent, and cooperative requirements of Foxl2a and Foxl2b in ovary development and maintenance of Zebrafish. Genetics. 2017;205(4):1551–72.

    CAS  Article  Google Scholar 

  22. 22.

    Pfennig F, Standke A, Gutzeit HO. The role of Amh signaling in teleost fish--multiple functions not restricted to the gonads. Gen Comp Endocrinol. 2015;223:87–107.

    CAS  Article  Google Scholar 

  23. 23.

    Guiguen Y, Fostier A, Piferrer F, Chang CF. Ovarian aromatase and estrogens: a pivotal role for gonadal sex differentiation and sex change in fish. General Comp Endocrinol. 2010;165(3):352–66.

    CAS  Article  Google Scholar 

  24. 24.

    Chen W, Liu L, Ge W. Expression analysis of growth differentiation factor 9 (Gdf9/gdf9), anti-mullerian hormone (Amh/amh) and aromatase (Cyp19a1a/cyp19a1a) during gonadal differentiation of the zebrafish, Danio rerio. Biol Reprod. 2017;96(2):401–13.

    Article  Google Scholar 

  25. 25.

    Soltis DE, Soltis PS, Tate JA. Advances in the study of polyploidy since plant speciation. New Phytol. 2004;161(1):173–91.

    CAS  Article  Google Scholar 

  26. 26.

    Comai L. The advantages and disadvantages of being polyploid. Nat Rev Genet. 2005;6(11):836–46.

    CAS  Article  Google Scholar 

  27. 27.

    Hegarty MJ, Barker GL, Wilson ID, Abbott RJ, Edwards KJ, Hiscock SJ. Transcriptome shock after interspecific hybridization in senecio is ameliorated by genome duplication. Curr Biol. 2006;16(16):1652–9.

    CAS  Article  Google Scholar 

  28. 28.

    Otto SP. The evolutionary consequences of polyploidy. Cell. 2007;131(3):452–62.

    CAS  Article  Google Scholar 

  29. 29.

    Kassahn KS, Dang VT, Wilkins SJ, Perkins AC, Ragan MA. Evolution of gene function and regulatory control after whole-genome duplication: comparative analyses in vertebrates. Genome Res. 2009;19(8):1404–18.

    CAS  Article  Google Scholar 

  30. 30.

    Li Z, Lu X, Gao Y, Liu S, Tao M, Xiao H, Qiao Y, Zhang Y, Luo J. Polyploidization and epigenetics. Chin Sci Bull. 2011;56(3):245–52.

    CAS  Article  Google Scholar 

  31. 31.

    Salmon A, Ainouche ML, Wendel JF. Genetic and epigenetic consequences of recent hybridization and polyploidy in Spartina (Poaceae). Mol Ecol. 2005;14(4):1163–75.

    CAS  Article  Google Scholar 

  32. 32.

    Chen ZJ, Ni Z. Mechanisms of genomic rearrangements and gene expression changes in plant polyploids. BioEssays. 2006;28(3):240–52.

    Article  Google Scholar 

  33. 33.

    Wang YD, Qin QB, Yang R, Sun WZ, Liu QW, Huo YY, Huang X, Tao M, Zhang C, Li T, et al. Hox genes reveal genomic DNA variation in tetraploid hybrids derived from Carassius auratus red var. (female) x Megalobrama amblycephala (male). BMC genetics. 2017;18(1):86.

    CAS  Article  Google Scholar 

  34. 34.

    Qin QB, Liu QW, Zhou YW, Wang CQ, Qin H, Zhao C, Liu SJ. Differential expression of HPG-axis genes in autotetraploids derived from red crucian carp Carassius auratus red var., female symbol x blunt snout bream Megalobrama amblycephala, male symbol. J Fish Biol. 2018;93(6):1082–9.

    CAS  Article  Google Scholar 

  35. 35.

    Song Q, Chen ZJ. Epigenetic and developmental regulation in plant polyploids. Curr Opin Plant Biol. 2015;24:101–9.

    CAS  Article  Google Scholar 

  36. 36.

    Navarro-Martin L, Vinas J, Ribas L, Diaz N, Gutierrez A, Di Croce L, Piferrer F. DNA methylation of the gonadal aromatase (cyp19a) promoter is involved in temperature-dependent sex ratio shifts in the European sea bass. PLoS Genet. 2011;7(12):e1002447.

    CAS  Article  Google Scholar 

  37. 37.

    Sun LX, Wang YY, Zhao Y, Wang H, Li N, Ji XS. Global DNA methylation changes in Nile Tilapia gonads during high temperature-induced masculinization. PLoS One. 2016;11(8):e0158483.

    Article  Google Scholar 

  38. 38.

    Chen X, He Y, Wang Z, Li J. Expression and DNA methylation analysis of cyp19a1a in Chinese sea perch Lateolabrax maculatus. Comp Biochem Physiol B Biochem Mol Biol. 2018;226:85–90.

    CAS  Article  Google Scholar 

  39. 39.

    Adams KL, Wendel JF. Novel patterns of gene expression in polyploid plants. Trends Genet. 2005;21(10):539–43.

    CAS  Article  Google Scholar 

  40. 40.

    Yu Z, Haberer G, Matthes M, Rattei T, Mayer KF, Gierl A, Torres-Ruiz RA. Impact of natural genetic variation on the transcriptome of autotetraploid Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2010;107(41):17809–14.

    CAS  Article  Google Scholar 

  41. 41.

    Zhang C, Ye L, Chen Y, Xiao J, Wu Y, Tao M, Xiao Y, Liu S. The chromosomal constitution of fish hybrid lineage revealed by 5S rDNA FISH. BMC Genet. 2015;16:140.

    Article  Google Scholar 

  42. 42.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    CAS  Article  Google Scholar 

  43. 43.

    Ding Y, He F, Wen H, Li J, Ni M, Chi M, Qian K, Bu Y, Zhang D, Si Y, et al. DNA methylation status of cyp17-II gene correlated with its expression pattern and reproductive endocrinology during ovarian development stages of Japanese flounder (Paralichthys olivaceus). Gene. 2013;527(1):82–8.

    CAS  Article  Google Scholar 

Download references


The authors thank all the researchers who helped this study. They are Shi Wang, Minghe Zhang, Huan Qin, Chun Zhao.


This research was financially supported by grants from the National Natural Science Foundation of China (Grant Nos. 31430088 and 31730098), the earmarked fund for China Agriculture Research System (Grant No. CARS-45), the Key Research and Development Program of Hunan Province (Grant No. 2018NK2072), Hunan Provincial Natural Science and Technology Major Project (Grant No. 2017NK1031), and the Cooperative Innovation Center of Engineering and New Products for Developmental Biology of Hunan Province (Grant No. 20134486). The funding bodies only provided the financial means to allow the authors to carry out the study, and the funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




XH, QBQ, SJL and CW1 designed the experiments, performed the analyses, performed the technical discussions; XH and KJG performed the experiments and drafted the manuscript; QBQ, SJL and CW1 modified the manuscript; XH, CW1, YWZ, QC participated in the sequence alignment; WJF, YYX, CW2 participated in phylogenetic analysis and discussions; XH, YDW, LC, MT collected the experimental materials. All authors read and approved the final manuscript. The CW1 corresponded to Chang Wu, and CW2 corresponded to Chongqing Wang.

Corresponding author

Correspondence to Shaojun Liu.

Ethics declarations

Ethics approval and consent to participate

All experiments were approved by the Animal Care Committee of Hunan Normal University and followed the guidelines statement of the Administration of Affairs Concerning Animal Experimentation of China. All samples were raised in natural ponds, all dissections were performed under MS-222 anaesthesia, and all efforts were made to minimize suffering.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

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

Huang, X., Qin, Q., Gong, K. et al. Comparative analyses of the Sox9a-Amh-Cyp19a1a regulatory Cascade in Autotetraploid fish and its diploid parent. BMC Genet 21, 35 (2020).

Download citation


  • Autopolyploidization
  • Sox9a-Amh-Cyp19a1a
  • DNA methylation
  • Genetic/epigenetic diversity