High male incidence and evolutionary implications of triploid form in northeast Asia Carassius auratus complex

A Fish of Multiple Faces, Which Show Us Enigmatic and Incredible Phenomena in Nature: Biology and Cytogenetics of the Genus Carassius

Sexual vs. asexual reproduction-unisexual vs. bisexual populations-diploid vs. polyploid biotypes-genetic vs. environmental sex determination: all these natural phenomena are associated with the genus of teleost fish, Carassius. This review places emphasis on two Carassius entities with completely different biological characteristics: one globally widespread and invasive Carassius gibelio, and the other C. carassius with a decreasing trend of natural occurrence. Comprehensive biological and cytogenetic knowledge of both entities, including the physical interactions between them, can help to balance the advantages of highly invasive and disadvantages of threatened species. For example, the benefits of a wide-ranged colonization can lead to the extinction of native species or be compensated by parasitic enemies and lead to equilibrium. This review emphasizes the comprehensive biology and cytogenetic knowledge and the importance of the Carassius genus as one of the most useful experimental vertebrate models for evolutionary biology and genetics. Secondly, the review points out that effective molecular cytogenetics should be used for the identification of various species, ploidy levels, and hybrids. The proposed investigation of these hallmark characteristics in Carassius may be applied in conservation efforts to sustain threatened populations in their native ranges. Furthermore, the review focuses on the consequences of the co-occurrence of native and non-native species and outlines future perspectives of Carassius research.

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.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

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.

Measurement of Chromosomal Arms and FISH Reveal Complex Genome Architecture and Standardized Karyotype of Model Fish, Genus Carassius

The widely distributed ray-finned fish genus Carassius is very well known due to its unique biological characteristics such as polyploidy, clonality, and/or interspecies hybridization. These biological characteristics have enabled Carassius species to be successfully widespread over relatively short period of evolutionary time. Therefore, this fish model deserves to be the center of attention in the research field. Some studies have already described the Carassius karyotype, but results are inconsistent in the number of morphological categories for individual chromosomes. We investigated three focal species: Carassius auratus, C. carassius and C. gibelio with the aim to describe their standardized diploid karyotypes, and to study their evolutionary relationships using cytogenetic tools. We measured length (q+plength) of each chromosome and calculated centromeric index (i value). We found: (i) The relationship between q+plength and i value showed higher similarity of C. auratus and C. carassius. (ii) The variability of i value within each chromosome expressed by means of the first quartile (Q1) up to the third quartile (Q3) showed higher similarity of C. carassius and C. gibelio. (iii) The fluorescent in situ hybridization (FISH) analysis revealed higher similarity of C. auratus and C. gibelio. (iv) Standardized karyotype formula described using median value (Q2) showed differentiation among all investigated species: C. auratus had 24 metacentric (m), 40 submetacentric (sm), 2 subtelocentric (st), 2 acrocentric (a) and 32 telocentric (T) chromosomes (24m+40sm+2st+2a+32T); C. carassius: 16m+34sm+8st+42T; and C. gibelio: 16m+22sm+10st+2a+50T. (v) We developed R scripts applicable for the description of standardized karyotype for any other species. The diverse results indicated unprecedented complex genomic and chromosomal architecture in the genus Carassius probably influenced by its unique biological characteristics which make the study of evolutionary relationships more difficult than it has been originally postulated.

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.

Coexistence of Three Divergent mtDNA Lineages in Northeast Asia Provides New Insights into Phylogeography of Goldfish (Carssius auratus)

Simple Summary

Goldfish (Carassius auratus) is a well-known fish as food and as a pet, which is also frequently used as experimental animal. A unique mtDNA sequence was detected in a sample from our experimental station, which motivated us to study genetic constitution of goldfish in Northeast Asia. Three divergent mtDNA lineages were confirmed to coexist in this region. Two of which corresponded to the known lineages (C2 and C6), which was consistent with the zoogeographical records that there were two sympatric subspecies in Amur river basin. However, the third one (lineage C7) was largely neglected in the previous studies. Our results suggested lineage C7 had a wide distribution from Central Asia to Northeast Asia.

Abstract

Goldfish (Carassius aurautus), which is a middle size cyprinid, widely distribute throughout Eurasia. Phylogeographic studies using mtDNA markers have revealed several divergent lineages within goldfish. In this study, mtDNA variations were determined to elucidate the phylogeographical pattern and genetic structure of goldfish in Northeast Asia. A total of 1054 individuals from Amur river basin were analyzed, which including five newly collected populations and four previously reported populations. Three distinct mtDNA lineages were identified in those samples, two of which corresponded to two known lineages C2 and C6, respectively. The third lineage referred to as C7, following six known lineages of goldfish in mainland Eurasia. AMOVA results suggested that most of the genetic variations were among lineages, rather than among populations or twice samplings. We noted that the control region (CR) and cytochrome b (cytb) sequences of lineage C7 have been reported in previous studies, respectively. However, the evolutionary position and distribution pattern of this lineage was not discussed in the context of the species. Our results showed that “odd” CR and “hidden” cytb sequences from Central Asia represent the same mtDNA lineage of goldfish. The known samples of C7 lineage were collected from Central Asia (Eastern Kazakhstan and Western Mongolia) to East Asia (Northeast China and Far East Russia), which suggested that it had a wider distribution, rather than limit in Central Asia.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Bisexual Fertile Triploid Zebrafish (Danio rerio): a Rare Case

Previous studies have suggested that artificially induced triploid zebrafish are exclusively male-biased. Owing to greatly inhibited gonadal development for the artificially induced triploid fish, they are regarded to be sterile in general. In this article, partially fertile bisexual triploid zebrafish are produced by suppressing extrusion of the second polar body by heat shock. Histological observation confirms that the early gonadal development of these triploid zebrafish is normal. Backcrossing and self-crossing are used to demonstrate that both the female and male triploid zebrafish have partial reproductive ability. Their dynamic of chromosomes during meiosis is revealed from the chromosome preparations of gonads. Examination of the expressed gonadal development-related genes shows some molecular evidence of the normal gonadal development in the triploid zebrafish. Clearly, these fertile bisexual triploid zebrafish can provide a unique system to study sex determination, as well as aneuploidy associated human diseases such as infertility and pregnancy loss.

Evolutionary dynamics of 18S and 5S rDNA in autotriploid Carassius auratus

The Carassius auratus (crucian carp) complex of the Dongting water system exhibits coexistence of diploid and triploid forms. As reported, triploid C. auratus is autotriploid origin. Ribosomal DNA (rDNA) with evolutionary conservation is widely used to study polyploidization. Here, we investigated genomic and transcribed rDNA sequences (18S and 5S) in diploid (2nCC, 2n = 100) and triploid (3nCC, 3n = 150) C. auratus. The results showed that the genetic traits and expression of 18S and 5S rDNA from 2nCC individuals were identified in 3nCC individuals. Moreover, pseudogenization of rDNA (18S and 5S) sequences were also observed in both 2nCC and 3nCC individuals, but expression of these variants was not detected. Based on the transcribed rDNA consensus sequence between 2nCC and 3nCC individuals, the functional secondary structures of 18S rRNA (expansion segments, ES6S) and 5S rRNA were predicted. These data demonstrated that complex evolutionary dynamics existed in the rDNA family of C. auratus. The evolutionary conservation of rDNA revealed that autotriploidization could not induce the divergence in Carassius taxa of the Dongting water system. These observations will expand our knowledge of the evolutionary dynamics of the rDNA family 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.

Differential expression of innate and adaptive immune genes in the survivors of three gibel carp gynogenetic clones after herpesvirus challenge

BACKGROUND

Accompanied with rapid growth and high density aquaculture, gibel carp has been seriously threatened by Carassius auratus herpesvirus (CaHV) since 2012. In previous study, distinct CaHV resistances and immune responses were revealed in the diseased individuals of three gibel carp gynogenetic clones (A+, F and H). However, little is known about the gene expression changes in the survivors after CaHV challenge, particularly their differences of innate and adaptive immune system between susceptible clone and resistant clone.

RESULTS

We firstly confirmed the CaHV carrier state in the survivors of three gibel carp clones after CaHV challenge by evaluating the abundances of five CaHV genes. The assay of viral loads indicated the resistant clone H possessed not only stronger resistance but also higher tolerance to CaHV. Then, 2818, 4047 and 3323 differentially expressed unigenes (DEUs) were screened from the head-kidney transcriptome profiles of survivors compared with controls from clone A+, F and H. GO and KEGG analysis suggested that a persistent immune response might sustain in resistant clone H and F, while susceptible clone A+ had a long-term impact on the circulatory system which was consistent with the major symptoms of bleeding caused by CaHV. Among the top 30 enriched pathways of specifically up-regulated DEUs in respective clones, 26, 7 and 15 pathways in clone H, F and A+ were associated with infections, diseases, or immune-related pathways respectively. In addition, 20 pathways in clone F belonged to "metabolism" or "biogenesis", and 7 pathways involved in "circulatory system" were enriched in clone A+. Significantly, we revealed the differential expression changes of IFN system genes and immunoglobulin (Ig) genes among the survivors of three clones. Finally, myosins and Igs were identified as co-expression modules which were positively or negatively correlated to CaHV viral loads respectively.

CONCLUSIONS

Our results revealed the common and distinct gene expression changes in immune and circulatory system in the survivors of three gibel carp gynogenetic clones with different CaHV resistances. The current study represents a paradigm of differential innate and adaptive immune reactions in teleost, and will be beneficial to the disease-resistance breeding of gibel carp.

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.

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 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.

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.

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.

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.

Type-IV antifreeze proteins are essential for epiboly and convergence in gastrulation of zebrafish embryos

Many organisms in extremely cold environments such as the Antarctic Pole have evolved antifreeze molecules to prevent ice formation. There are four types of antifreeze proteins (AFPs). Type-IV antifreeze proteins (AFP4s) are present also in certain temperate and even tropical fish, which has raised a question as to whether these AFP4s have important functions in addition to antifreeze activity. Here we report the identification and functional analyses of AFP4s in cyprinid fish. Two genes, namely afp4a and afp4b coding for AFP4s, were identified in gibel carp (Carassius auratus gibelio) and zebrafish (Danio rerio). In both species, afp4a and afp4b display a head-to-tail tandem arrangement and share a common 4-exonic gene structure. In zebrafish, both afp4a and afp4b were found to express specifically in the yolk syncytial layer (YSL). Interestingly, afp4a expression continues in YSL and digestive system from early embryos to adults, whereas afp4b expression is restricted to embryogenesis. Importantly, we have shown by using afp4a-specific and afp4b-specifc morpholino knockdown and cell lineage tracing approaches that AFP4a participates in epiboly progression by stabilizing yolk cytoplasmic layer microtubules, and AFP4b is primarily related to convergence movement. Therefore, both AFP4 proteins are essential for gastrulation of zebrafish embryos. Our current results provide first evidence that AFP such as AFP4 has important roles in regulating developmental processes besides its well-known function as antifreeze factors.

The transcriptomes of the crucian carp complex (Carassius auratus) provide insights into the distinction between unisexual triploids and sexual diploids

Both sexual reproduction and unisexual reproduction are adaptive strategies for species survival and evolution. Unisexual animals have originated largely by hybridization, which tends to elevate their heterozygosity. However, the extent of genetic diversity resulting from hybridization and the genomic differences that determine the type of reproduction are poorly understood. In Carassius auratus, sexual diploids and unisexual triploids coexist. These two forms are similar morphologically but differ markedly in their modes of reproduction. Investigation of their genomic differences will be useful to study genome diversity and the development of reproductive mode. We generated transcriptomes for the unisexual and sexual populations. Genes were identified using homology searches and an ab initio method. Estimation of the synonymous substitution rate in the orthologous pairs indicated that the hybridization of gibel carp occurred 2.2 million years ago. Microsatellite genotyping in each individual from the gibel carp population indicated that most gibel carp genes were not tri-allelic. Molecular function and pathway comparisons suggested few gene expansions between them, except for the progesterone-mediated oocyte maturation pathway, which is enriched in gibel carp. Differential expression analysis identified highly expressed genes in gibel carp. The transcriptomes provide information on genetic diversity and genomic differences, which should assist future studies in functional genomics.

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

Polyploidy lineages, despite very rare in vertebrates, have been proposed to play significant role in speciation and evolutionary success, but the occurrence history and consequences are still largely unknown. In this study, we used the conserved Dmrt1 to analyze polyploidy occurrence and evolutionary process in polyploid gibel carp. We identified two divergent Dmrt1 genes and respectively localized the two genes on three homologous chromosomes. Subsequently, the corresponding full-length cDNAs and genomic sequences of Dmrt1 genes were also characterized from the closely related species including Carassius auratus auratus and Cyprinus carpio, and their two Dmrt1 genes were respectively localized on two homologous chromosomes. Significantly, the evolutionary relationship analyses among cDNA and genomic DNA sequences of these Dmrt1 genes revealed two rounds of polyploidy origins in the gibel carp: an early polyploidy might result in an common tetraploid ancestor of Carassius auratus gibelio, Carassius auratus auratus and Cyprinus carpio before 18.49 million years ago (Mya), and an late polyploidy might occur from evolutionary branch of Carassius auratus at around 0.51 Mya, which lead to the occurrence of the hexaploid gibel carp. Therefore, this study provides clear genetic evidence for understanding occurrence time and historical process of polyploidy in polyploid vertebrates.

C1q-like factor, a target of miR-430, regulates primordial germ cell development in early embryos of Carassius auratus

C1q-like is a significant maternal factor of TNF/C1q super-family, and the abundant protein has been observed in both mature eggs of Carassius auratus and Carassius auratus gibelio, but its biological function in early embryo development has remained unclear. In this study, we firstly revealed a high level of maternal C1q-like transcript existence only in mature eggs of Carassius auratus, whereas no any maternal C1q-like transcript was observed in that of Carassius auratus gibelio. During embryonic development, the C1q-like zygotic expression begins around cardiopalmus stage in embryos of both Carassius auratus and Carassius auratus gibelio. Then, we examined the biological role of C1q-like by morpholino-mediated knockdown in early embryo development. Knockdown of CaOC1q resulted in a significant reduction of primordial germ cells (PGCs) in Carassius auratus, as shown by whole mount in situ hybridization with vasa-specific RNA probe, fluorescence immunostaining of vasa protein, and GFP imaging of the GFP-nanos1-3'UTR mRNA reporter. In vitro and in vivo evidence indicated that a microRNA, miR-430 could repress the C1q-like expression and PGC development. These data suggest that C1q-like should be a direct target of miR-430 and play an essential role in PGC development of Carassius auratus.

Genetic differentiation and efficient sex-specific marker development of a pair of Y- and X-linked markers in yellow catfish

Pf62-Y and Pf62-X is a pair of allelic Y chromosome-linked and X chromosome-linked markers, and have been used to identify YY super-males, XY males and XX females for commercial production of all-male populations in yellow catfish (Pelteobagrus fulvidraco). However, the SCAR primers used previously have only two nucleotide difference, which restricts the wide utility because of nucleotide polymorphism. In this study, a continuous 8102 bp Pf62-Y sequence and a 5362 bp Pf62-X sequence have been cloned by genome walking, and significant genetic differentiation has been revealed between the corresponding X and Y chromosome allele sequences. Moreover, three pairs of primers were designed to efficiently identify YY super-males, XY males and XX females in an artificial breeding population, and to distinguish XY males and XX females in various wild populations. Together, the three new sex-specific genetic markers develop a highly stable and efficient method for genetic sex identification and sex control application in sustainable aquaculture of all-male yellow catfish.

Oocyte-specific H2A variant H2af1o is required for cell synchrony before midblastula transition in early zebrafish embryos

Oocyte-specific histone variants have been expected to play significant roles in early embryonic development, but the exact evidence and the biological function have remained unclear. Here, we present evidence that H2af1o, an oocyte-specific H2A variant, is required for cell synchrony before midblastula transition in early zebrafish embryos. The H2A variant is oocyte specific, peaks in mature eggs, and is supplied to early embryos. We constructed a series of deletion plasmids of the zebrafish h2af1o tagged with EGFP and determined the main key function regions including nuclear localization signal of N-terminal 25 amino acids and nucleosome binding region of 110-122 amino acid sequence in the C-terminus by microinjecting them into one-cell-stage zebrafish embryos. In comparison with ubiquitous H2A.X, the H2af1o was revealed to confer a more open structure than canonical H2A in the nucleosomes. Furthermore, we conducted the h2af1o-specific morpholino knockdown analysis in early embryos of zebrafish and revealed its biological function for maintaining cell synchrony division because the H2af1o deficiency disturbed cell synchrony in early cleavages before midblastula transition. Therefore, our current findings provided the first case to understand the biological function of maternal oocyte-specific histone variants in vertebrates.