Evolutionary history of two divergent Dmrt1 genes reveals two rounds of polyploidy origins in gibel carp

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

Efficiently Editing Multiple Duplicated Homeologs and Alleles for Recurrent Polyploids

Research on the evolutionary fate of duplicated genes in recurrent polyploids is scarce due to the difficulties in disentangling the different homeologs and alleles of duplicated genes. This chapter describes the detailed procedures to identify different homeologs and alleles of duplicated genes, to analyze their molecular characteristics, and to reveal their functional divergence by gene editing with CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9). Using the gene editing approach, we efficiently constructed multiple knockout mutant lines with single or simultaneously disrupted different homeologs or alleles in a recurrent polyploid fish, demonstrating its usability for targeting and mutating multiple divergent homeologs and alleles in recurrent duplicated genomes.

Equilibrated evolution of the mixed auto-/allopolyploid haplotype-resolved genome of the invasive hexaploid Prussian carp

Understanding genome evolution of polyploids requires dissection of their often highly similar subgenomes and haplotypes. Polyploid animal genome assemblies so far restricted homologous chromosomes to a 'collapsed' representation. Here, we sequenced the genome of the asexual Prussian carp, which is a close relative of the goldfish, and present a haplotype-resolved chromosome-scale assembly of a hexaploid animal. Genome-wide comparisons of the 150 chromosomes with those of two ancestral diploid cyprinids and the allotetraploid goldfish and common carp revealed the genomic structure, phylogeny and genome duplication history of its genome. It consists of 25 syntenic, homeologous chromosome groups and evolved by a recent autoploid addition to an allotetraploid ancestor. We show that de-polyploidization of the alloploid subgenomes on the individual gene level occurred in an equilibrated fashion. Analysis of the highly conserved actinopterygian gene set uncovered a subgenome dominance in duplicate gene loss of one ancestral chromosome set.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

Genome sequencing and haplotype assembly of two cyprinid teleosts, a sexual tetraploid and an unisexual hexaploid, reveal insights into the evolutionary mechanisms underpinning the reproductive success of unisexual polyploid vertebrates.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Phylogeographic relationships and the evolutionary history of the Carassius auratus complex with a newly born homodiploid raw fish (2nNCRC)

Background

An important aspect of studying evolution is to understand how new species are formed and their uniqueness is maintained. Hybridization can lead to the formation of new species through reorganization of the adaptive system and significant changes in phenotype. Interestingly, eight stable strains of 2nNCRC derived from interspecies hybridization have been established in our laboratory. To examine the phylogeographical pattern of the widely distributed genus Carassius across Eurasia and investigate the possible homoploid hybrid origin of the Carassius auratus complex lineage in light of past climatic events, the mitochondrial genome (mtDNA) and one nuclear DNA were used to reconstruct the phylogenetic relationship between the C. auratus complex and 2nNCRC and to assess how demographic history, dispersal and barriers to gene flow have led to the current distribution of the C. auratus complex.

Results

As expected, 2nNCRC had a very close relationship with the C. auratus complex and similar morphological characteristics to those of the C. auratus complex, which is genetically distinct from the other three species of Carassius. The estimation of divergence time and ancestral state demonstrated that the C. auratus complex possibly originated from the Yangtze River basin in China. There were seven sublineages of the C. auratus complex across Eurasia and at least four mtDNA lineages endemic to particular geographical regions in China. The primary colonization route from China to Mongolia and the Far East (Russia) occurred during the Late Pliocene, and the diversification of other sublineages of the C. auratus complex specifically coincided with the interglacial stage during the Early and Mid-Pleistocene in China.

Conclusion

Our results support the origin of the C. auratus complex in China, and its wide distribution across Eurasia was mainly due to natural Pleistocene dispersal and recent anthropogenic translocation. The sympatric distribution of the ancestral area for both parents of 2nNCRC and the C. auratus complex, as well as the significant changes in the structure of pharyngeal teeth and morphological characteristics between 2nNCRC and its parents, imply that homoploid hybrid speciation (HHS) for C. auratus could likely have occurred in nature. The diversification pattern indicated an independent evolutionary history of the C. auratus complex, which was not separated from the most recent common ancestor of C. carassius or C. cuvieri. Considering that the paleoclimate oscillation and the development of an eastward-flowing drainage system during the Pliocene and Pleistocene in China provided an opportunity for hybridization between divergent lineages, the formation of 2nNCRC in our laboratory could be a good candidate for explaining the HHS of C. auratus in nature.

Investigation of putative antimicrobial peptides in Carassius gibel, revealing a practical approach to screening antimicrobials

Antimicrobial peptides (AMPs) and their mimics are rapidly gaining attention as a new class of antimicrobials due to their clinical potential. AMPs are widely distributed throughout nature and participate in the innate host defense. In this study, 18 AMPs, including 3 β-defensins, 3 hepcidins, 4 liver-expressed antimicrobial peptide 2 (LEAP-2) compounds, 4 g-type lysozymes, 2 c-type lysozymes, and 2 NK-lysins, were identified from the genome of Carassius auratus by a homologous search and were further classified based on their fundamental structural features and molecular phylogeny. C. auratus AMPs were found to be ubiquitously distributed in all tested tissues and showed similar expression profiles, with the exception of β-defensins, when RT-qPCR was used to investigate the tissue distribution of AMPs in healthy Carassius gibel. In addition, the expression levels of NK-lysin genes in the tested tissues tended to be upregulated upon bacterial and viral infection when representative NK-lysins were chosen to examine their relative expression levels in various tissues. Importantly, the synthetic peptide caNKL2_102-119, which targets the functional domain of saposin B in caNK-lysins, could effectively counter Aeromonas hydrophila, Staphylococcus aureus, and Escherichia coli with minimum inhibitory concentration (MIC) values of 3-6 μg/mL, as well as inhibit the proliferation of spring viraemia of carp virus (SVCV). These results provide potential targets for antibiotic-free breeding in the aquaculture industry.

Two Duplicated Ptpn6 Homeologs Cooperatively and Negatively Regulate RLR-Mediated IFN Response in Hexaploid Gibel Carp

Src homology region 2 domain-containing phosphatase 1 (SHP1), encoded by the protein tyrosine phosphatase nonreceptor type 6 (ptpn6) gene, belongs to the family of protein tyrosine phosphatases (PTPs) and participates in multiple signaling pathways of immune cells. However, the mechanism of SHP1 in regulating fish immunity is largely unknown. In this study, we first identified two gibel carp (Carassius gibelio) ptpn6 homeologs (Cgptpn6-A and Cgptpn6-B), each of which had three alleles with high identities. Then, relative to Cgptpn6-B, dominant expression in adult tissues and higher upregulated expression of Cgptpn6-A induced by polyinosinic-polycytidylic acid (poly I:C), poly deoxyadenylic-deoxythymidylic (dA:dT) acid and spring viremia of carp virus (SVCV) were uncovered. Finally, we demonstrated that CgSHP1-A (encoded by the Cgptpn6-A gene) and CgSHP1-B (encoded by the Cgptpn6-B gene) act as negative regulators of the RIG-I-like receptor (RLR)-mediated interferon (IFN) response via two mechanisms: the inhibition of CaTBK1-induced phosphorylation of CaMITA shared by CgSHP1-A and CgSHP1-B, and the autophagic degradation of CaMITA exclusively by CgSHP1-A. Meanwhile, the data support that CgSHP1-A and CgSHP1-B have sub-functionalized and that CgSHP1-A overwhelmingly dominates CgSHP1-B in the process of RLR-mediated IFN response. The current study not only sheds light on the regulative mechanism of SHP1 in fish immunity, but also provides a typical case of duplicated gene evolutionary fates.

Deciphering reticulate evolution of the largest group of polyploid vertebrates, the subfamily cyprininae (Teleostei: Cypriniformes)

Despite the rarity of polyploidy in animals, some groups with polyploid species exhibit complicated and interesting patterns of reticulate evolution. Here we focus on fishes in the subfamily Cyprininae, the largest polyploid group of vertebrates. The large number of polyploid taxa poses significant challenges for phylogenetic and evolutionary studies on this subfamily. In this study, we cloned and sequenced three single-copy nuclear loci to investigate the evolution of polyploidy in the Cyprininae. Topologies of nuclear gene trees were compared with a newly reconstructed mitochondrial tree. The data provided herein corroborate the hybrid origins of the tribes Torini, Cyprinini, Spinibarbini, Barbini, and also Probarbini. Based on results from this study and previous studies, we hypothesize that at least 13 independent polyploidization events have occurred during the evolution of the Cyprininae. We offer hypotheses on the origin of each polyploid group and show that a diploid group or the diploid ancestor of a polyploid group might have served as progenitor of one or two other polyploid groups. The evolutionary history of Cyprinine (since its first divergence) can be divided into three stages: the "Diploid stage" (69.2-43.4 Ma or million years ago), the "Tetraploidization stage" (43.4-18.9 Ma), and the "Hexaploidization stage" (18.9 Ma to present). The second stage is when all tetraploidization events happened, while the last stage is when all hexaploidization events and most genus- or species-specific polyploidization events occurred. The post-polyploidization dynamics of polyploid groups are complicated and warrant more genomic level studies. We showed that the subfamily Cyprininae may represent a more complicated polyploid system than most, if not all, other vertebrates and some plants, if one or more of the following factors are considered: number of polyploid species, number of different ploidy levels, and number and type of independent polyploidization events.

Divergent Antiviral Mechanisms of Two Viperin Homeologs in a Recurrent Polyploid Fish

Polyploidy and subsequent diploidization provide genomic opportunities for evolutionary innovations and adaptation. The researches on duplicated gene evolutionary fates in recurrent polyploids have seriously lagged behind that in paleopolyploids with diploidized genomes. Moreover, the antiviral mechanisms of Viperin remain largely unclear in fish. Here, we elaborate the distinct antiviral mechanisms of two viperin homeologs (Cgviperin-A and Cgviperin-B) in auto-allo-hexaploid gibel carp (Carassius gibelio). First, Cgviperin-A and Cgviperin-B showed differential and biased expression patterns in gibel carp adult tissues. Subsequently, using co-immunoprecipitation (Co-IP) screening analysis, both CgViperin-A and CgViperin-B were found to interact with crucian carp (C. auratus) herpesvirus (CaHV) open reading frame 46 right (ORF46R) protein, a negative herpesvirus regulator of host interferon (IFN) production, and to promote the proteasomal degradation of ORF46R via decreasing K63-linked ubiquitination. Additionally, CgViperin-B also mediated ORF46R degradation through autophagosome pathway, which was absent in CgViperin-A. Moreover, we found that the N-terminal α-helix domain was necessary for the localization of CgViperin-A and CgViperin-B at the endoplasmic reticulum (ER), and the C-terminal domain of CgViperin-A and CgViperin-B was indispensable for the interaction with degradation of ORF46R. Therefore, the current findings clarify the divergent antiviral mechanisms of the duplicated viperin homeologs in a recurrent polyploid fish, which will shed light on the evolution of teleost duplicated genes.

Genomic anatomy of male-specific microchromosomes in a gynogenetic fish

Unisexual taxa are commonly considered short-lived as the absence of meiotic recombination is supposed to accumulate deleterious mutations and hinder the creation of genetic diversity. However, the gynogenetic gibel carp (Carassius gibelio) with high genetic diversity and wide ecological distribution has outlived its predicted extinction time of a strict unisexual reproduction population. Unlike other unisexual vertebrates, males associated with supernumerary microchromosomes have been observed in gibel carp, which provides a unique system to explore the rationales underlying male occurrence in unisexual lineage and evolution of unisexual reproduction. Here, we identified a massively expanded satellite DNA cluster on microchromosomes of hexaploid gibel carp via comparing with the ancestral tetraploid crucian carp (Carassius auratus). Based on the satellite cluster, we developed a method for single chromosomal fluorescence microdissection and isolated three male-specific microchromosomes in a male metaphase cell. Genomic anatomy revealed that these male-specific microchromosomes contained homologous sequences of autosomes and abundant repetitive elements. Significantly, several potential male-specific genes with transcriptional activity were identified, among which four and five genes displayed male-specific and male-biased expression in gonads, respectively, during the developmental period of sex determination. Therefore, the male-specific microchromosomes resembling common features of sex chromosomes may be the main driving force for male occurrence in gynogenetic gibel carp, which sheds new light on the evolution of unisexual reproduction.

Genotypic Males Play an Important Role in the Creation of Genetic Diversity in Gynogenetic Gibel Carp

Unisexual lineages are commonly considered to be short-lived in the evolutionary process as accumulation of deleterious mutations stated by Muller’s ratchet. However, the gynogenetic hexaploid gibel carp (Carassius gibelio) with existence over 0.5 million years has wider ecological distribution and higher genetic diversity than its sexual progenitors, which provides an ideal model to investigate the underlying mechanisms on countering Muller’s ratchet in unisexual taxa. Unlike other unisexual lineages, the wild populations of gibel carp contain rare and variable proportions of males (1–26%), which are determined via two strategies including genotypic sex determination and temperature-dependent sex determination. Here, we used a maternal gibel carp from strain F to be mated with a genotypic male from strain A⁺, a temperature-dependent male from strain A⁺, and a male from another species common carp (Cyprinus carpio), respectively. When the maternal individual was mated with the genotypic male, a variant of gynogenesis was initiated, along with male occurrence, accumulation of microchromosomes, and creation of genetic diversity in the offspring. When the maternal individual was mated with the temperature-dependent male and common carp, typical gynogenesis was initiated that all the offspring showed the same genetic information as the maternal individual. Subsequently, we found out that the genotypic male nucleus swelled and contacted with the female nucleus after fertilization although it was extruded from the female nucleus eventually, which might be associated with the genetic variation in the offspring. These results reveal that genotypic males play an important role in the creation of genetic diversity in gynogenetic gibel carp, which provides insights into the evolution of unisexual reproduction.

Functional Divergence of Multiple Duplicated Foxl2 Homeologs and Alleles in a Recurrent Polyploid Fish

Evolutionary fates of duplicated genes have been widely investigated in many polyploid plants and animals, but research is scarce in recurrent polyploids. In this study, we focused on foxl2, a central player in ovary, and elaborated the functional divergence in gibel carp (Carassius gibelio), a recurrent auto-allo-hexaploid fish. First, we identified three divergent foxl2 homeologs (Cgfoxl2a-B, Cgfoxl2b-A, and Cgfoxl2b-B), each of them possessing three highly conserved alleles and revealed their biased retention/loss. Then, their abundant sexual dimorphism and biased expression were uncovered in hypothalamic-pituitary-gonadal axis. Significantly, granulosa cells and three subpopulations of thecal cells were distinguished by cellular localization of CgFoxl2a and CgFoxl2b, and the functional roles and the involved process were traced in folliculogenesis. Finally, we successfully edited multiple foxl2 homeologs and/or alleles by using CRISPR/Cas9. Cgfoxl2a-B deficiency led to ovary development arrest or complete sex reversal, whereas complete disruption of Cgfoxl2b-A and Cgfoxl2b-B resulted in the depletion of germ cells. Taken together, the detailed cellular localization and functional differences indicate that Cgfoxl2a and Cgfoxl2b have subfunctionalized and cooperated to regulate folliculogenesis and gonad differentiation, and Cgfoxl2b has evolved a new function in oogenesis. Therefore, the current study provides a typical case of homeolog/allele diversification, retention/loss, biased expression, and sub-/neofunctionalization in the evolution of duplicated genes driven by polyploidy and subsequent diploidization from the recurrent polyploid fish.

Abnormal meiosis in fertile and sterile triploid cyprinid fish

Meiosis is the key process for producing mature gametes. A natural fertile triploid Carassius auratus population (3nDTCC) and an artificially derived sterile triploid crucian carp (3nCC) have been previously observed, providing suitable model organisms for investigating meiosis characteristics in triploid fish. In the present study, the microstructures and ultrastructures of spermatogenesis were studied in these fishes. TdT-mediated dUTP nick end labeling detection was performed to investigate the apoptosis of spermatocytes. Fluorescence in situ hybridization was employed to trace chromatin pairing. In addition, the mRNA expressions of cell cycle-related genes (i.e., cell division control 2 and cell cycle protein B) were determined by quantitative realtime polymerase chain reaction to illustrate the molecular mechanism of abnormal meiosis in the 3nCC. The results showed that the 3nCC undergoes an irregular prophase I, with the chromosomes distributed in a unipolar radial manner and exhibiting partial pairing, hindered metaphase I, and degenerated cells in the subsequent stages. Meanwhile, the 3nDTCC presented a relatively regular meiotic prophase I with complete conjugate chromosome pairs and chromosomes distributed along the karyotheca, which were presented as a ring structure by slicing. Only the spreads with 130-150 irregular chromosomes can be easily detected in the 3nDTCC, suggesting that it may undergo an abnormal metaphase I. This study provides new insights into the meiosis of fertile and sterile triploid cyprinid fish.

The genomic structure of the highly-conserved dmrt1 gene in Solea senegalensis (Kaup, 1868) shows an unexpected intragenic duplication

Knowing the factors responsible for sex determination in a species has significant theoretical and practical implications; the dmrt1 gene (Doublesex and Mab-3 (DM)-related Transcription factor 1) plays this role in diverse animal species. Solea senegalensis is a commercially important flat fish in which females grow 30% faster than males. It has 2n = 42 chromosomes and an XX / XY chromosome system for sex determination, without heteromorph chromosomes but with sex proto-chromosome. In the present study, we are providing the genomic structure and nucleotide sequence of dmrt1 gene obtained from cDNA from male and female adult gonads. A cDNA of 2027 containing an open-reading frame (ORF) of 1206 bp and encoding a 402 aa protein it is described for dmrt1 gene of S. senegalensis. Multiple mRNA isoforms indicating a high variable system of alternative splicing in the expression of dmrt1 of the sole in gonads were studied. None isoforms could be related to sex of individuals. The genomic structure of the dmrt1 of S. senegalensis showed a gene of 31400 bp composed of 7 exons and 6 introns. It contains an unexpected duplication of more than 10399 bp, involving part of the exon I, exons II and III and a SINE element found in the sequence that it is proposed as responsible for the duplication. A mature miRNA of 21 bp in length was localized at 336 bp from exon V. Protein-protein interacting networks of the dmrt1 gene showed matches with dmrt1 protein from Cynoglossus semilaevis and a protein interaction network with 11 nodes (dmrt1 plus 10 other proteins). The phylogenetic relationship of the dmrt1 gene in S. senegalensis is consistent with the evolutionary position of its species. The molecular characterization of this gene will enhance its functional analysis and the understanding of sex differentiation in Solea senegalensis and other flatfish.

Dynamic and Differential Expression of Duplicated Cxcr4/Cxcl12 Genes Facilitates Antiviral Response in Hexaploid Gibel Carp

Chemokine receptor cxcr4 and its ligand cxcl12 have evolved two paralogs in the teleost lineage. In this study, we have identified four duplicated cxcr4 and cxcl12 genes from hexaploid gibel carp, Carassius gibelio, respectively. Cgcxcr4bs and Cgcxcl12as were dynamically and differentially expressed in immune-related tissues, and significantly up-regulated in head kidney and spleen after crucian carp herpesvirus (CaHV) infection. Blocking Cxcr4/Cxcl12 axis by injecting AMD3100 brought more severe bleeding symptom and lower survival rate in CaHV-infected fish. AMD3100 treatment also suppressed the up-regulation of key antiviral genes in head kidney and spleen, and resulted in more acute replication of CaHV in vivo. Consistently, the similar suppression of up-regulated expression of key antiviral genes were also observed in CAB cells treated by AMD3100 after poly(I:C) stimulation. Finally, MAPK3 and JAK/STAT were identified as the possible pathways that CgCxcr4s and CgCxcl12s participate in to promote the antiviral response in vitro.

Sox Gene Family Revealed Genetic Variations in Autotetraploid Carassius auratus

The Sox gene family encoded transcription factors that played key roles in developmental processes in vertebrates. To further understand the evolutionary fate of the Sox gene family in teleosts, the Sox genes were comprehensively characterized in fish of different ploidy levels, including blunt snout bream (2n = 48, Megalobrama amblycephala, BSB), goldfish (2n = 100, Carassius auratus red var., 2nRCC), and autotetraploid C. auratus (4n = 200, 4nRCC). The 4nRCC, which derived from the whole genome duplication (WGD) of 2nRCC, were obtained through the distant hybridization of 2nRCC (♀) × BSB (♂). Compared with the 26 Sox genes in zebrafish (2n = 50, Danio rerio), 26, 47, and 92 putative Sox genes were identified in the BSB, 2nRCC, and 4nRCC genomes, respectively, and classified into seven subfamilies (B1, B2, C, D, E, F, and K). Comparative analyses showed that 89.36% (42/47) of Sox genes were duplicated in 2nRCC compared with those in BSB, while 97.83% (90/92) of Sox genes were duplicated in 4nRCC compared with those in 2nRCC, meaning the Sox gene family had undergone an expansion in BSB, 2nRCC, and 4nRCC, respectively, following polyploidization events. In addition, potential gene loss, genetic variations, and paternal parent SNP locus insertion occurred during the polyploidization events. Our data provided new insights into the evolution of the Sox gene family in polyploid vertebrates after several rounds of WGD events.

Regain of sex determination system and sexual reproduction ability in a synthetic octoploid male fish

Polyploids in vertebrates are generally associated with unisexual reproduction, but the direct consequences of polyploidy on sex determination system and reproduction mode remain unknown. Here, we synthesized a group of artificial octoploids between unisexual gynogenetic hexaploid Carassius gibelio and sexual tetraploid Carassius auratus. The synthetic octoploids were revealed to have more than 200 chromosomes, in which 50 chromosomes including the X/Y sex determination system were identified to transfer from sexual tetraploid C. auratus into the unisexual gynogenetic hexaploid C. gibelio. Significantly, a few synthetic octoploid males were found to be fertile, and one octoploid male was confirmed to regain sexual reproduction ability, which exhibits characteristics that are the same to sexual reproduction tetraploid males, such as 1:1 sex ratio occurrence, meiosis completion and euploid sperm formation in spermatogenesis, as well as normal embryo development and gene expression pattern during embryogenesis. Therefore, the current finding provides a unique case to explore the effect of sex determination system incorporation on reproduction mode transition from unisexual gynogenesis to sexual reproduction along with genome synthesis of recurrent polyploidy in vertebrates.

Stable Genome Incorporation of Sperm-derived DNA Fragments in Gynogenetic Clone of Gibel Carp

How unisexual animals eliminate deleterious mutations to avoid dead ends is one of the most interesting puzzles in evolutionary genetics. Incorporation of microchromosomes derived from exogenous sperm had been observed in gynogenetic animals, but little is known about their detailed process and hereditary fate. Here, we show a stable genome incorporation case in an artificial clone F of gynogenetic gibel carp (Carassius gibelio). A total of 12 exogenous DNA fragments were screened through a read depth-dependent comparison strategy and confirmed to be specific to the clone F and the paternal blunt snout bream (Megalobrama amblycephala Yin) by SCAR (sequence characterized amplified regions) marker detection. Moreover, these sperm-derived DNA fragments were not detected in some samples in early gynogenetic generations, but they were found to exist in all examined individuals through artificial gynogenetic selections of 13 generations, implying that they might have stably incorporated into the genome of clone F. Furthermore, chromosome localization and sequence characterization indicate that the largest fragment CgA22_34 is derived from blunt snout bream non-LTR retrotransposon and durably incorporated into only one of three homologous chromosomes of gibel carp clone F. Our results suggest that the incorporated sperm-derived DNA fragments by allogynogenesis should increase genetic diversity and introduce new traits into unisexual animals which will benefit genetic breeding of gibel carp. During the process, transposable elements (TEs) may play significant roles in shaping the genome structures. Simultaneously, the incorporated DNA fragments are able to be used as genetic markers to perform selective breeding programs in aquaculture practices of gibel carp.

Molecular characterization, computational analysis and expression profiling of Dmrt1 gene in Indian major carp, Labeo rohita (Hamilton 1822)

Sexual dimorphism of fish morphology, physiology and behavior is diverse and complex in nature. Doublesex and mab-3 related transcription factor (Dmrt) is a large protein family whose function is sexual development and differentiation in vertebrates. Here, we report a full-length cDNA sequence of Labeo rohita (rohu) Dmrt1 of 907 bp length having 798 bp of open reading frame encoding 265 amino acids. The molecular weight of rohu DMRT1 protein was found to be 28.74 KDa and isoelectric point was 7.53. DMRT1 protein contains 23 positively and 24 negatively charged amino acids with a GRAVY score of -0.618. A characteristic DM domain was found in DMRT1 protein, which is a novel DNA-binding domain. Phylogenetic analysis showed maximum similarity with Cyprinus carpio when compared with DMRT1 of other vertebrates. Molecular docking study identified active sites to be targeted for drug designing. Rohu DMRT1 was observed to interact with other proteins such as FOXL2, CYP19a1a, AMH and SOX9a. Differential expression study revealed higher expression in testis tissue implying its role in male sex differentiation and testicular development. The information generated in the present work could facilitate further research to resolve the issues related to gonadal maturation and reproduction of commercially important aquaculture species.

Analysis on the Meiosis-Related Gene (Dmc1, Ph1) Expression in Autotriploid Carassius auratus

Triploid is usually considered to be unable to perform normal meiosis due to the abnormal behavior of the three sets of chromosomes. But autotriploid Carassius auratus in the Dongting water system (3n = 150, abbreviated as 3nCC) can perform normal meiosis. In artificial autotriploid Carassius auratus (3n = 150, abbreviated as 3nRR), female individuals undergo normal meiosis and produce mature gametes, while male individuals cannot. To better understand the effects of triploidization on meiosis in fish, we study the structure, methylation level, and expression level of meiosis-related genes (Dmc1, Ph1) in diploid Carassius auratus (2n = 100, abbreviated as 2nCC), Carassius auratus red var.(2n = 100, abbreviated as RCC), 3nCC and 3nRR. The results show that, compared with their diploid ancestors (2nCC and RCC), Dmc1 and Ph1 genes are hypomethylated in all 3nCC and female 3nRR, while are hypermethylated in male 3nRR. Correspondingly, Dmc1 and Ph1 genes are highly expressed in all 3nCC and female 3nRR, while are lowly expressed in male 3nRR. These results indicate that high expression of meiosis-related genes can contribute to restoration of bivalent pairing during meiosis in autotriploid Carassius auratus. This study provides new insights into the effect of DNA methylation on the fertility in triploid fish.

Female-to-male sex reversal in orange-spotted grouper (Epinephelus coioides) caused by overexpressing of Amh in vivo

A variety of mechanisms are involved in sex determination in vertebrates. The orange-spotted grouper (Epinephelus coioides), a teleost fish, functions first as females and later as a male and is an ideal model to investigate the regulation of sexual fate. Here, we report female-to-male sex reversal in juvenile orange-spotted groupers caused by overexpressing anti-Müllerian hormone (Amh). Tissue distribution analyses showed that amh and amhrII primarily expressed in the gonad, and expression level in the testis was much higher than that in the ovary. In gonads, the expression of amh was located in the Sertoli cells around spermatogonia of the testis and in the zona pellucida of the mature ovary, and the expression of amhrII was located in the Sertoli cells of the testis and in the oocytes of the ovary. Decrease in female-related genes and serum 17β-estradiol level, increase in male-related genes and serum 11-ketotestosterone, ovarian regression, and spermatogonia proliferation were observed during plasmid feeding experiment. These results illustrate that amh overexpression plasmid feeding can induce a female-to-male transition in grouper.

Diverse and variable sex determination mechanisms in vertebrates

Sex is prevalent in nature and sex determination is one of the most fundamental biological processes, while the way of initiating female and male development exhibits remarkable diversity and variability across vertebrates. The knowledge on why and how sex determination mechanisms evolve unusual plasticity remains limited. Here, we summarize sex determination systems, master sex-determining genes and gene-regulatory networks among vertebrates. Recent research advancements on sex determination system transition are also introduced and discussed in some non-model animals with multiple sex determination mechanisms. This review will provide insights into the origin, transition and evolutionary adaption of different sex determination strategies in vertebrates, as well as clues for future perspectives in this field.

Whole Genome Incorporation and Epigenetic Stability in a Newly Synthetic Allopolyploid of Gynogenetic Gibel Carp

Allopolyploidization plays an important role in speciation, and some natural or synthetic allopolyploid fishes have been extensively applied to aquaculture. Although genetic and epigenetic inheritance and variation associated with plant allopolyploids have been well documented, the relative research in allopolyploid animals is scarce. In this study, the genome constitution and DNA methylation inheritance in a newly synthetic allopolyploid of gynogenetic gibel carp were analyzed. The incorporation of a whole genome of paternal common carp sperm in the allopolyploid was confirmed by genomic in situ hybridization, chromosome localization of 45S rDNAs, and sequence comparison. Pooled sample-based methylation sensitive amplified polymorphism (MSAP) revealed that an overwhelming majority (98.82%) of cytosine methylation patterns in the allopolyploid were inherited from its parents of hexaploid gibel carp clone D and common carp. Compared to its parents, 11 DNA fragments in the allopolyploid were proved to be caused by interindividual variation, recombination, deletion, and mutation through individual sample-based MSAP and sequencing. Contrast to the rapid and remarkable epigenetic changes in most of analyzed neopolyploids, no cytosine methylation variation was detected in the gynogenetic allopolyploid. Therefore, the newly synthetic allopolyploid of gynogenetic gibel carp combined genomes from its parents and maintained genetic and epigenetic stability after its formation and subsequently seven successive gynogenetic generations. Our current results provide a paradigm for recurrent polyploidy consequences in the gynogenetic allopolyploid animals.

Differential interferon system gene expression profiles in susceptible and resistant gynogenetic clones of gibel carp challenged with herpesvirus CaHV

Interferon (IFN) system plays a vital role in the first line of defense against viruses. In this study, we first identified multiple transcripts of 15 IFN system genes, including PRRs (TLR2, TLR3, RIG-I, and LGP2), PRR-mediated IFN signal pathway (MyD88, MITA, and MAVS), IFN regulatory factors (IRF1, IRF3, IRF7, and IRF9), IFNs (IFNφ1 and IFNφ3), and ISGs (Mx and viperin), and one transcript of TLR9 in de novo transcriptome assembly data of gibel carp head-kidney. Multiple nucleotide alignments and phylogenetic analysis of common region showed that the transcripts of every of the 15 IFN system genes were classified into two homologs with distinctly divergent sequences, indicating that hexaploid gibel carp may be an allopolyploid. During Carassius auratus herpesvirus (CaHV) infection, gibel carp resistant clone H significantly suppressed CaHV replication with markedly less viral loads than those in highly susceptible clone A⁺ and moderately resistant clone F. Then, qPCR analyses were performed to reveal their differential and dynamic expression changes during CaHV infection in head kidney, spleen and liver among three gibel carp gynogenetic clones. Through qPCR and hierarchical clustering analysis, 8 genes, such as RIG-Is, LGP2s, IRF1-B, IRF3s, IRF7s, IRF9-B, Mxs, and viperins, were identified as candidate resistant-related genes. They remarkably increased their expression in immune tissues of three clones after CaHV infection. Significantly, the up-regulation folds of these genes in clone A⁺, F and H were related to their resistance ability to CaHV, progressively increasing from susceptible clone to resistant clone at 1 dpi. The positive correlation to the resistance ability suggested that resistant clone H immediately triggered stronger IFN response. IFNφ3 showed a different dynamic change and was sharply induced in moderately resistant clone F at 3 dpi. The other 5 IFN system genes (TLR2, TLR3, TLR9, MyD88, and MITA) maintained a low expression level after CaHV challenge. Interestingly, the A or B copies/homologs of almost these IFN system genes exhibited differential transcript abundance in immune tissue after CaHV challenge, suggesting A or B homologs might occur dominant or biased expression of homeologs during gibel carp evolution. These data provide candidate resistant-related genes for disease-resistance breeding of gibel carp.

Differential expression and functional diversification of diverse immunoglobulin domain-containing protein (DICP) family in three gynogenetic clones of gibel carp

Diverse immunoglobulin (Ig) domain-containing protein (DICP) family is a novel bony fish-specific multi-gene family encoding diversified immune receptors. However, their function and the implication of binding partners remain unknown. In this study, we first identified 28 DICPs from three gibel carp gynogenetic clones and revealed their high variability and clone-specific feature. After crucian carp herpesvirus (CaHV) infection, these DICPs were significantly upregulated in head kidney, kidney and spleen. The up-regulation folds in clone A⁺, F and H were related to the susceptibility to CaHV, progressively increasing from resistant clone to susceptible clone. Overexpression of gibel carp DICPs inhibited interferon (IFN) and viperin promoter-driven luciferase activity. The additions of E. coli extracts and lipid A significantly enhanced the inhibition effect. In addition, gibel carp DICPs can interact with SHP-1 and SHP-2. These findings suggest that gible carp DICPs, as inhibitory receptors, might specifically recognize lipid A, and then interact with SHP-1 and SHP-2 to inhibit the induction of IFN and ISGs.

Origin and transition of sex determination mechanisms in a gynogenetic hexaploid fish

Most vertebrates reproduce sexually, and plastic sex determination mechanisms including genotypic sex determination (GSD) and environmental sex determination (ESD) have been extensively revealed. However, why sex determination mechanisms evolve diversely and how they correlate with diverse reproduction strategies remain largely unclear. Here, we utilize the superiority of a hexaploid gibel carp (Carassius gibelio) that is able to reproduce by unisexual gynogenesis and contains a rare but diverse proportion of males to investigate these puzzles. A total of 2248 hexaploid specimens were collected from 34 geographic wild populations throughout mainland China, in which 24 populations were revealed to contain 186 males with various incidences ranging from 1.2 to 26.5%. Subsequently, the proportion of temperature-dependent sex determination (TSD) was revealed to be positively correlated to average annual temperature in wild populations, and male incidence in lab gynogenetic progenies was demonstrated to increase with the increasing of larval rearing temperature. Meanwhile, extra microchromosomes were confirmed to play genotypic male determination role as previously reported. Thereby, GSD and TSD were found to coexist in gibel carp, and the proportions of GSD were observed to be much higher than that of TSD in sympatric wild populations. Our findings uncover a potential new mechanism in the evolution of sex determination system in polyploid vertebrates with unisexual gynogenesis ability, and also reveal a possible association of sex determination mechanism transition between TSD and GSD and reproduction mode transition between unisexual gynogenesis and bisexual reproduction.

Distinct sperm nucleus behaviors between genotypic and temperature-dependent sex determination males are associated with replication and expression-related pathways in a gynogenetic fish

BACKGROUND

Coexistence and transition of diverse sex determination strategies have been revealed in some ectothermic species, but the variation between males caused by different sex determination strategies and the underlying mechanism remain unclear. Here, we used the gynogenetic gibel carp (Carassius gibelio) with both genotypic sex determination (GSD) and temperature-dependent sex determination (TSD) strategies to illustrate this issue.

RESULTS

We found out that males of GSD and TSD in gibel carp had similar morphology, testicular histology, sperm structure and sperm vitality. However, when maternal individuals were mated with males of GSD, sperm nucleus swelling and fusing with the female pronucleus were observed in the fertilized eggs. On the contrary, when maternal individuals were mated with males of TSD, sperm nucleus remained in the condensed status throughout the whole process. Subsequently, semen proteomics analysis unveiled that DNA replication and gene expression-related pathways were inhibited in the sperm from males of TSD compared to males of GSD, and most differentially expressed proteins associated with DNA replication, transcription and translation were down-regulated. Moreover, via BrdU incorporation and immunofluorescence detection, male nucleus replication was revealed to be present in the fertilized eggs by the sperm from males of GSD, but absent in the fertilized eggs by the sperm from males of TSD.

CONCLUSIONS

These findings indicate that DNA replication and gene expression-related pathways are associated with the distinct sperm nucleus development behaviors in fertilized eggs in response to the sperm from males of GSD and TSD. And this study is the first attempt to screen the differences between males determined via GSD and TSD in gynogenetic species, which might give a hint for understanding evolutionary adaption of diverse sex determination mechanisms in unisexual vertebrates.

Numerous mitochondrial DNA haplotypes reveal multiple independent polyploidy origins of hexaploids in Carassius species complex

Evolutionary trajectory and occurrence history of polyploidy have been extensively studied in plants, but they remain quite elusive in vertebrates. Here, we sampled and gathered 4,159 specimens of polyploid Carassius species complex including 1,336 tetraploids and 2,823 hexaploids from a large geographic scale (49 localities) across East Asia, and identified a huge number of 427 diverse haplotypes of mitochondrial control region, in which 74 haplotypes with total occurrence frequency up to 75.498% were shared by hexaploids and tetraploids. Significantly, these diverse haplotypes were clustered into four major lineages, and many haplotypes of hexaploids and tetraploids were intermixed in every lineage. Moreover, the evolutionary trajectory and occurrence history of four different lineages were revealed by a simplified time-calibrated phylogenetic tree, and their geographic distribution frequencies and haplotype diversity were also analyzed. Furthermore, lineage C and D were revealed to undergo population expansion throughout mainland China. Therefore, our current data indicate that hexaploids should undergo multiple independent polyploidy origins from sympatric tetraploids in the polyploid Carassius species complex across East Asia.

Distinct and Cooperative Roles of amh and dmrt1 in Self-Renewal and Differentiation of Male Germ Cells in Zebrafish

Spermatogenesis is a fundamental process in male reproductive biology and depends on precise balance between self-renewal and differentiation of male germ cells. However, the regulative factors for controlling the balance are poorly understood. In this study, we examined the roles of amh and dmrt1 in male germ cell development by generating their mutants with Crispr/Cas9 technology in zebrafish. Amh mutant zebrafish displayed a female-biased sex ratio, and both male and female amh mutants developed hypertrophic gonads due to uncontrolled proliferation and impaired differentiation of germ cells. A large number of proliferating spermatogonium-like cells were observed within testicular lobules of the amh-mutated testes, and they were demonstrated to be both Vasa- and PH3-positive. Moreover, the average number of Sycp3- and Vasa-positive cells in the amh mutants was significantly lower than in wild-type testes, suggesting a severely impaired differentiation of male germ cells. Conversely, all the dmrt1-mutated testes displayed severe testicular developmental defects and gradual loss of all Vasa-positive germ cells by inhibiting their self-renewal and inducing apoptosis. In addition, several germ cell and Sertoli cell marker genes were significantly downregulated, whereas a prominent increase of Insl3-positive Leydig cells was revealed by immunohistochemical analysis in the disorganized dmrt1-mutated testes. Our data suggest that amh might act as a guardian to control the balance between proliferation and differentiation of male germ cells, whereas dmrt1 might be required for the maintenance, self-renewal, and differentiation of male germ cells. Significantly, this study unravels novel functions of amh gene in fish.

Analyzing homoeolog expression provides insights into the rediploidization event in gynogenetic hybrids of Carassius auratus red var. × Cyprinus carpio

Rediploidization is considered to be a part of the evolutionary history of allotetraploids, and resulted in the emergence of novel epigenetic regulatory activities. To study the changing patterns of gene expression following the reduction of a genome by 50%, we used RNA-seq and quantitative real-time PCR (qPCR) to investigate total gene expression and homoeolog expression in three hybrids of a C. auratus red var. (2n = 100, ♀) (R) and C. carpio (2n = 100, ♂) (C) (i.e., F1, F18, and G4) and their original parents. A comparison of homoeolog expression between G4 and F18 identified 7 genes (0.22%) that exhibited novel R/C homoeolog expression patterns in G4, while 4 genes (0.12%) were affected by R/C homoeolog silencing. We determined the direction and extent of the homoeolog expression bias (HEB). The C-HEB genes (i.e., nrp1a and igf1) and R-HEB genes (i.e., fgf23 and esm1) provided insights into the effects of the dominance of one parental homoeolog expression on growth regulation. This dominance may contribute to the rapid growth of G4 fish. Our findings may be relevant for clarifying the relationship between growth heterosis and differences in homoeolog expression patterns.

Genomic and transcriptomic approaches to study immunology in cyprinids: What is next?

Accelerated by the introduction of Next-Generation Sequencing (NGS), a number of genomes of cyprinid fish species have been drafted, leading to a highly valuable collective resource of comparative genome information on cyprinids (Cyprinidae). In addition, NGS-based transcriptome analyses of different developmental stages, organs, or cell types, increasingly contribute to the understanding of complex physiological processes, including immune responses. Cyprinids are a highly interesting family because they comprise one of the most-diversified families of teleosts and because of their variation in ploidy level, with diploid, triploid, tetraploid, hexaploid and sometimes even octoploid species. The wealth of data obtained from NGS technologies provides both challenges and opportunities for immunological research, which will be discussed here. Correct interpretation of ploidy effects on immune responses requires knowledge of the degree of functional divergence between duplicated genes, which can differ even between closely-related cyprinid fish species. We summarize NGS-based progress in analysing immune responses and discuss the importance of respecting the presence of (multiple) duplicated gene sequences when performing transcriptome analyses for detailed understanding of complex physiological processes. Progressively, advances in NGS technology are providing workable methods to further elucidate the implications of gene duplication events and functional divergence of duplicates genes and proteins involved in immune responses in cyprinids. We conclude with discussing how future applications of NGS technologies and analysis methods could enhance immunological research and understanding.

Distinct herpesvirus resistances and immune responses of three gynogenetic clones of gibel carp revealed by comprehensive transcriptomes

Background

Gibel carp is an important aquaculture species in China, and a herpesvirus, called as Carassius auratus herpesvirus (CaHV), has hampered the aquaculture development. Diverse gynogenetic clones of gibel carp have been identified or created, and some of them have been used as aquaculture varieties, but their resistances to herpesvirus and the underlying mechanism remain unknown.

Results

To reveal their susceptibility differences, we firstly performed herpesvirus challenge experiments in three gynogenetic clones of gibel carp, including the leading variety clone A⁺, candidate variety clone F and wild clone H. Three clones showed distinct resistances to CaHV. Moreover, 8772, 8679 and 10,982 differentially expressed unigenes (DEUs) were identified from comparative transcriptomes between diseased individuals and control individuals of clone A⁺, F and H, respectively. Comprehensive analysis of the shared DEUs in all three clones displayed common defense pathways to the herpesvirus infection, activating IFN system and suppressing complements. KEGG pathway analysis of specifically changed DEUs in respective clones revealed distinct immune responses to the herpesvirus infection. The DEU numbers identified from clone H in KEGG immune-related pathways, such as “chemokine signaling pathway”, “Toll-like receptor signaling pathway” and others, were remarkably much more than those from clone A⁺ and F. Several IFN-related genes, including Mx1, viperin, PKR and others, showed higher increases in the resistant clone H than that in the others. IFNphi3, IFI44-like and Gig2 displayed the highest expression in clone F and IRF1 uniquely increased in susceptible clone A⁺. In contrast to strong immune defense in resistant clone H, susceptible clone A⁺ showed remarkable up-regulation of genes related to apoptosis or death, indicating that clone A⁺ failed to resist virus offensive and evidently induced apoptosis or death.

Conclusions

Our study is the first attempt to screen distinct resistances and immune responses of three gynogenetic gibel carp clones to herpesvirus infection by comprehensive transcriptomes. These differential DEUs, immune-related pathways and IFN system genes identified from susceptible and resistant clones will be beneficial to marker-assisted selection (MAS) breeding or molecular module-based resistance breeding in gibel carp.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3945-6) contains supplementary material, which is available to authorized users.

Wider geographic distribution and higher diversity of hexaploids than tetraploids in Carassius species complex reveal recurrent polyploidy effects on adaptive evolution

Polyploidy roles on adaptive evolution and ecological novelty have been extensively studied in plants but remained unclear in vertebrates owing to the rare polyploidy incidences. Here, a huge number of 3105 specimens in Carassius species complex including 2211 hexaploids and 894 tetraploids were sampled from 34 locations through mainland China. And hexaploids had wider geographic distribution than tetraploids especially in the areas with high altitude, high latitude and low annual precipitation. Then, an approximate 1050 bp transferrin (tf) fragments were amplified from all the samples, and 526 tf alleles were identified from a total of 37260 sequences at last. Intriguingly, higher nucleotide diversity of tf alleles in hexaploids than in tetraploids was revealed. Moreover, via phylogenetic analysis of tf alleles, potential origin center of Carassius species complex was deduced to be Yangtze River basin and hexaploids should undergo multiple independent polyploidy origins from sympatric tetraploids. These findings indicate that the hexaploids might possess stronger environmental adaptation and ecological novelty than the tetraploids, which provide an association paradigm of recurrent polyploidy and ecological context in polyploid vertebrates.

A novel male-specific SET domain-containing gene setdm identified from extra microchromosomes of gibel carp males

Various genes have been screened on extra chromosomes, but their molecular characterization, expression pattern and biological function are still unclear. Here, we utilized a male-specific sequence of polyploid gibel carp (Carassius gibelio) to identify a novel male-specific SET (Su(var)3-9, Enhancer-of-zeste, Trithorax) domain-containing gene setdm on extra microchromosomes of gibel carp males. And setdm was characterized in molecule and expression aspects, in which its expression was specific to testis, and had relative high transcription during middle/late stages of testis development. Moreover, prominent expression of Setdm in spermatogenic cells was observed in testis through immunofluorescence co-localization analysis. These results suggest that biological function of setdm might be related to testis development and spermatogenesis of gibel carp. Additionally, the homeologous gene setdmf of setdm, was also characterized, and its expression was gonad-specific, in which its expressed product was detected to mainly distribute in gametogenic cells of testis and ovary, and to have dynamic expression pattern similar to that of setdm. Based on the current results, we propose that the novel male-specific setdm on extra microchromosomes might be functional divergence gene of the gonad-specific setdmf. Therefore, these findings will help us to further understand evolutionary fate and functional role of genes on extra microchromosomes.

The genetic map of goldfish (Carassius auratus) provided insights to the divergent genome evolutions in the Cyprinidae family

A high-density linkage map of goldfish (Carassius auratus) was constructed using RNA-sequencing. This map consists of 50 linkage groups with 8,521 SNP markers and an average resolution of 0.62 cM. Approximately 84% of markers are in protein-coding genes orthologous to zebrafish proteins. We performed comparative genome analysis between zebrafish and medaka, common carp, grass carp, and goldfish to study the genome evolution events in the Cyprinidae family. The comparison revealed large synteny blocks among Cyprinidae fish and we hypothesized that the Cyprinidae ancestor undergone many inter-chromosome rearrangements after speciation from teleost ancestor. The study also showed that goldfish genome had one more round of whole genome duplication (WGD) than zebrafish. Our results illustrated that most goldfish markers were orthologous to genes in common carp, which had four rounds of WGD. Growth-related regions and genes were identified by QTL analysis and association study. Function annotations of the associated genes suggested that they might regulate development and growth in goldfish. This first genetic map enables us to study the goldfish genome evolution and provides an important resource for selective breeding of goldfish.

Molecular characterization and expression pattern of a germ cell marker gene dnd in gibel carp (Carassius gibelio)

As a germ cell marker gene, Dead end (dnd) has been identified and characterized in many vertebrates. Recently, we created a complete germ cell-depleted gonad model by the dnd-specific morpholino-mediated knockdown approach, and revealed sex-biased gene expression alteration through utilizing unisexual gynogenetic superiority in polyploid gibel carp. However, dnd and its expression pattern are still unclear in the gibel carp. In this study, we further analyzed molecular characterization of gibel carp dnd and its dynamic expression pattern during gametogenesis and embryogenesis. Similar to other homologs in vertebrates, gibel carp dnd contains a conserved RRM motif and five other motifs, and is highly evolutionary conserved in genomic organization and neighborhood gene synteny. RT-PCR and Western blot analyses showed its gonad-specific expression intensively in testis and ovary. Section in situ hybridization (SISH) and immunofluorescence localization revealed its dynamic expression pattern specific to oogenic cells and spermatogenetic cells during oogenesis and spermatogenesis. Moreover, its temporal and spatial distribution specific to PGCs were also demonstrated by RT-PCR and whole mount in situ hybridization (WISH) during embryogenesis. Therefore, gibel carp Dnd is a conserved germ cell marker during gametogenesis, and its maternal transcript is also a useful marker for tracing PGC specification and migration.

Extra Microchromosomes Play Male Determination Role in Polyploid Gibel Carp

Sex is generally determined by sex chromosomes in vertebrates, and sex chromosomes exhibit the most rapidly-evolving traits. Sex chromosome evolution has been revealed previously in numerous cases, but the association between sex chromosome origin and the reproduction mode transition from unisexual to sexual reproduction remains unclear. Here, we have isolated a male-specific sequence via analysis of amplified fragment length polymorphism from polyploid gibel carp (Carassius gibelio), a species that not only has the ability to reproduce unisexually but also contains males in wild populations. Subsequently, we have found through FISH analysis that males have several extra microchromosomes with repetitive sequences and transposable elements when compared to females. Moreover, we produced sex-reversed physiological females with a male-specific marker by using estradiol hormone treatment, and two gynogenetic families were established from them. In addition, the male incidence rates of two gynogenetic families were revealed to be closely associated with the extra microchromosome number of the sex-reversed physiological females. These results suggest that the extra microchromosomes in males might resemble a common feature of sex chromosomes and might play a significant role in male determination during the evolutionary trajectory of the reproduction mode transition from unisexual to sexual reproduction in the polyploid fish.

Complete depletion of primordial germ cells in an All-female fish leads to Sex-biased gene expression alteration and sterile All-male occurrence

Background

Gynogenesis is one of unisexual reproduction modes in vertebrates, and produces all-female individuals with identical genetic background. In sexual reproduction vertebrates, the roles of primordial germ cells on sexual dimorphism and gonadal differentiation have been largely studied, and two distinct functional models have been proposed. However, the role of primordial germ cells remains unknown in unisexual animals, and it is also unclear whether the functional models in sexual reproduction animals are common in unisexual animals.

Results

To solve these puzzles, we attempt to utilize the gynogenetic superiority of polyploid Carassius gibelio to create a complete germ cell-depleted gonad model by a similar morpholino-mediated knockdown approach used in other examined sexual reproduction fishes. Through the germ cell-depleted gonad model, we have performed comprehensive and comparative transcriptome analysis, and revealed a complete alteration of sex-biased gene expression. Moreover, the expression alteration leads to up-regulation of testis-biased genes and down-regulation of ovary-biased genes, and results in the occurrence of sterile all-males with testis-like gonads and secondary sex characteristics in the germ cell-depleted gynogenetic Carassius gibelio.

Conclusions

Our current results have demonstrated that unisexual gynogenetic embryos remain keeping male sex determination information in the genome, and the complete depletion of primordial germ cells in the all-female fish leads to sex-biased gene expression alteration and sterile all-male occurrence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2130-z) contains supplementary material, which is available to authorized users.

Identification of Sex-Specific Markers Reveals Male Heterogametic Sex Determination in Pseudobagrus ussuriensis

Comprehending sex determination mechanism is a first step for developing sex control breeding biotechnologies in fish. Pseudobagrus ussuriensis, one of bagrid catfishes in Bagridae, had been observed to have about threefold size dimorphism between males and females, but its sex determination mechanism had been unknown. In this study, we firstly used the amplified fragment length polymorphism (AFLP)-based screening approach to isolate a male-specific DNA fragment and thereby identified a 10,569 bp of male-specific sequence and a 10,365 bp of female-related sequence by genome walking in the bagrid catfish, in which a substantial genetic differentiation with 96.35 % nucleotide identity was revealed between them. Subsequently, a high differentiating region of 650 bp with only 70.26 % nucleotide identity was found from the corresponding two sequences, and three primer pairs of male-specific marker, male and female-shared marker with different length products in male and female genomes, and female-related marker were designed. Significantly, when these markers were used to identify genetic sex of the bagrid catfish, only male individuals was detected to amplify the male-specific marker fragment, and female-related marker was discovered to produce dosage association in females and in males. Our current data provide significant genetic evidence that P. ussuriensis has heterogametic XY sex chromosomes in males and homogametic XX sex chromosomes in females. Therefore, sex determination mechanism of P. ussuriensis is male heterogametic XX/XY system.

Comparative Evolution of Duplicated Ddx3 Genes in Teleosts: Insights from Japanese Flounder, Paralichthys olivaceus

Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication event occurred in the stem lineage of ray-finned fishes. This teleost-specific genome duplication is thought to be responsible for the biological diversification of ray-finned fishes. DEAD-box polypeptide 3 (DDX3) belongs to the DEAD-box RNA helicase family. Although their functions in humans have been well studied, limited information is available regarding their function in teleosts. In this study, two teleost Ddx3 genes were first identified in the transcriptome of Japanese flounder (Paralichthys olivaceus). We confirmed that the two genes originated from teleost-specific genome duplication through synteny and phylogenetic analysis. Additionally, comparative analysis of genome structure, molecular evolution rate, and expression pattern of the two genes in Japanese flounder revealed evidence of subfunctionalization of the duplicated Ddx3 genes in teleosts. Thus, the results of this study reveal novel insights into the evolution of the teleost Ddx3 genes and constitute important groundwork for further research on this gene family.

Meiosis completion and various sperm responses lead to unisexual and sexual reproduction modes in one clone of polyploid Carassius gibelio

Unisexual polyploid vertebrates are commonly known to reproduce by gynogenesis, parthenogenesis, or hybridogenesis. One clone of polyploid Carassius gibelio has been revealed to possess multiple modes of unisexual gynogenesis and sexual reproduction, but the cytological and developmental mechanisms have remained unknown. In this study, normal meiosis completion was firstly confirmed by spindle co-localization of β-tubulin and Spindlin. Moreover, three types of various nuclear events and development behaviors were revealed by DAPI staining and BrdU-incorporated immunofluorescence detection during the first mitosis in the fertilized eggs by three kinds of different sperms. They include normal sexual reproduction in response to sperm from the same clone male, typical unisexual gynogenesis in response to sperm from the male of another species Cyprinus carpio, and an unusual hybrid-similar development mode in response to sperm from another different clone male. Based on these findings, we have discussed cytological and developmental mechanisms on multiple reproduction modes in the polyploid fish, and highlighted evolutionary significance of meiosis completion and evolutionary consequences of reproduction mode diversity in polyploid vertebrates.

A comprehensive transcriptome provides candidate genes for sex determination/differentiation and SSR/SNP markers in yellow catfish

Sex dimorphic growth pattern has significant theory and application implications in fish. Recently, a Y- and X-specific allele marker-assisted sex control technique has been developed for mass production of all-male population in yellow catfish (Pelteobagrus fulvidraco), but the genetic information for sex determination and sex control breeding has remained unclear. Here, we attempted to provide the first insight into a comprehensive transcriptome covering multiple tissues from XX females, XY males, and YY super-males of yellow catfish by using 454 GS-FLX platform, for a better assembly and gene coverage. A total of 1,202,933 high quality reads (about 540 Mbp) were obtained and assembled into 28,297 contigs and 141,951 singletons. BLASTX searches against the NCBI non-redundant protein database (nr) led a total of 52,564 unique sequences including 18,748 contigs and 33,816 singletons to match 25,669 known or predicted unique proteins. All of them with annotated function were categorized by gene ontology (GO) analysis, and 712 were assigned to reproduction and reproductive process. Some potential genes relevant to reproductive system including steroid hormone biosynthesis and GnRH (gonadotropin-releasing hormone) signaling pathway were further identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis; and at least 21 sex determination and differentiation-related genes, such as Dmrt1, Sox9a/b, Cyp19b, WT1, and AMH were identified and characterized. Additionally, a total of 82,794 simple sequence repeats (SSRs), 26,450 single nucleotide polymorphisms (SNPs), and 4,145 insertions and deletions (INDELs) were revealed from the transcriptome data. Therefore, the current transcriptome resources highlight further studies on sex-control breeding in yellow catfish and will benefit future studies on reproduction and sex determination in teleost fish.

Genetic basis and biotechnological manipulation of sexual dimorphism and sex determination in fish

Aquaculture has made an enormous contribution to the world food production, especially to the sustainable supply of animal proteins. The utility of diverse reproduction strategies in fish, such as the exploiting use of unisexual gynogenesis, has created a typical case of fish genetic breeding. A number of fish species show substantial sexual dimorphism that is closely linked to multiple economic traits including growth rate and body size, and the efficient development of sex-linked genetic markers and sex control biotechnologies has provided significant approaches to increase the production and value for commercial purposes. Along with the rapid development of genomics and molecular genetic techniques, the genetic basis of sexual dimorphism has been gradually deciphered, and great progress has been made in the mechanisms of fish sex determination and identification of sex-determining genes. This review summarizes the progress to provide some directive and objective thinking for further research in this field.