Sex chromosome and sex locus characterization in goldfish, Carassius auratus (Linnaeus, 1758)

Cytogenetic Analysis of the Fish Genus Carassius Indicates Divergence, Fission, and Segmental Duplication as Drivers of Tandem Repeat and Microchromosome Evolution

Fishes of the genus Carassius are useful experimental vertebrate models for the study of evolutionary biology and cytogenetics. Carassius demonstrates diverse biological characteristics, such as variation in ploidy levels and chromosome numbers, and presence of microchromosomes. Those Carassius polyploids with ≥150 chromosomes have microchromosomes, but the origin of microchromosomes, especially in European populations, is unknown. We used cytogenetics to study evolution of tandem repeats (U1 and U2 small nuclear DNAs and H3 histone) and microchromosomes in Carassius from the Czech Republic. We tested the hypotheses whether the number of tandem repeats was affected by polyploidization or divergence between species and what mechanism drives evolution of microchromosomes. Tandem repeats were found in tetraploid and hexaploid Carassius gibelio, and tetraploid Carassius auratus and Carassius carassius in conserved numbers, with the exception of U1 small nuclear DNA in C. auratus. This conservation indicates reduction and/or loss in the number of copies per locus in hexaploids and may have occurred by divergence rather than polyploidization. To study the evolution of microchromosomes, we used the whole microchromosome painting probe from hexaploid C. gibelio and hybridized it to tetraploid and hexaploid C. gibelio, and tetraploid C. auratus and C. carassius. Our results revealed variation in the number of microchromosomes in hexaploids and indicated that the evolution of the Carassius karyotype is governed by macrochromosome fissions followed by segmental duplication in pericentromeric areas. These are potential mechanisms responsible for the presence of microchromosomes in Carassius hexaploids. Differential efficacy of one or both of these mechanisms in different tetraploids could ensure variability in chromosome number in polyploids in general.

Polygenic sex determination in vertebrates - is there any such thing?

Genetic sex determination (SD) in most vertebrates is controlled by a single master sex gene, which ensures a 1:1 sex ratio. However, more complex systems abound, and several have been ascribed to polygenic SD (PSD), in which many genes at different loci interact to produce the sexual phenotype. Here we examine claims for PSD in vertebrates, finding that most constitute transient states during sex chromosome turnover, or aberrant systems in species hybrids. To avoid confusion about terminology, we propose a consistent nomenclature for genetic SD systems.

Origin of scarlet gynogenetic triploid Carassius fish: Implications for conservation of the sexual-gynogenetic complex

Conservation of sperm-dependent asexual (gynogenetic) species is challenging due to their complicated ecological dynamics, which requires the stable coexistence with their sperm-providing sexual relatives, who often share similar niches. A symbolic but vulnerable gynogenetic animal is the scarlet Carassius fish, or Hibuna, which is mainly found in Lake Harutori on Hokkaido, Japan. Although Hibuna in Lake Harutori has been protected as a symbol of the Natural Monument of Japan, it has recently suffered population decline. To establish an effective conservation strategy for Hibuna, we investigated its origin, reproductive mode, and genetic diversity, with reference to the surrounding wild populations, using nuclear microsatellites and mitochondrial gene sequences. Our genetic analyses revealed that the main ploidy of Hibuna was triploid or tetraploid, and it reproduces gynogenetically. However, no co-existing sexual diploid Carassius was detected among our samples, suggesting that the sexual diploids and the gynogenetic population including Hibuna would be at risk of co-extirpation. In addition, Hibuna showed high genetic/clonal diversity and most Hibuna had nonindigenous mitochondrial haplotypes that are mostly identical to those reported from goldfish. These results indicate that Hibuna most probably originated from hybridization between indigenous gynogenetic triploids and goldfish introduced about 100 years ago, involving rare sexual reproduction. This spontaneous long-term field experiment exemplifies the recently documented diversification process of gynogenetic Carassius via complex interploidy gene flow. Although the priority to be placed on the conservation of Hibuna is controversial, the maintenance of gynogenetic Carassius, including Hibuna, requires strategic conservation of sexual populations.

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.

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.

The haplotype-resolved assembly of the asexual invasive Prussian carp shows six genome copies (AAABBB), evolved from two ancestral species by a recent self-addition (AB) to its hybrid-tetraploid (AABB) goldfish ancestor. Equilibrated gene loss led to subgenome dominance.

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.

Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome

Pseudobagrus ussuriensis is an aquaculture catfish with significant sexual dimorphism. In this study, a chromosome-level genome with a size of 741.97 Mb was assembled for female P. ussuriensis. A total of 26 chromosome-level contigs covering 97.34% of the whole-genome assembly were obtained with an N50 of 28.53 Mb and an L50 of 11. A total of 24,075 protein-coding genes were identified, with 91.54% (22,039) genes being functionally annotated. Based on the genome assembly, four chromosome evolution clusters of catfishes were identified and the formation process of P. ussuriensis chromosomes was predicted. A total of 55 sex-related quantitative trait loci (QTLs) with a phenotypic variance explained value of 100% were located on chromosome 8 (chr08). The QTLs and other previously identified sex-specific markers were located in a sex-determining region of 16.83 Mb (from 6.90 to 23.73 Mb) on chr08, which was predicted as the X chromosome. The sex-determining region comprised 554 genes, with 135 of which being differently expressed between males and females/pseudofemales, and 16 candidate sex-determining genes were screened out. The results of this study provided a useful chromosome-level genome for genetic, genomic and evolutionary studies of P. ussuriensis, and also be useful for further studies on sex-determination mechanism analysis and sex-control breeding of this fish.

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.

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

Identification of Sex-Specific Markers Through 2b-RAD Sequencing in the Sea Urchin (Mesocentrotus nudus)

Sex-specific markers play an important role in revealing sex-determination mechanism. Sea urchin (Mesocentrotus nudus) is an economically important mariculture species in several Asian countries and its gonads are the sole edible parts for people. However, growth rate and immunocompetence differ by sex in this species, sex-specific markers have not been identified, and the sex-determination mechanism of sea urchin remains undetermined. In this study, type IIB endonuclease restriction-site associated DNA sequencing (2b-RAD-seq) and a genome survey of M. nudus were performed, and three female-specific markers and three female heterogametic single nucleotide polymorphism (SNP) loci were identified. We validated these sex-specific markers via PCR amplification in a large number of individuals, including wild and artificially bred populations. Several open reading frames (ORFs) were predicted, although there are no potential genes known for sex determination and sex differentiation within the scaffold in which the sex-specific markers are located. Importantly, the female-specific sequences and female heterozygous SNP loci indicate that a female heterogametic and male homogametic ZW/ZZ sex-determination system should exist in M. nudus. The results provide a solid basis for revealing the sex-determination mechanism of this species, and open up new possibilities for developing sex-control breeding in sea urchin.

Is the use of recombinant cGnRH may be a future alternative to control the fish spawning? Let us go with the goldfish example

The use of recombinant gonadotropin-releasing hormone (rGnRH) has very rarely been tested in fish to promote spawning. This study evaluated the impact of recombinant chicken gonadotropin-releasing hormone (rcGnRH) with metoclopramide on the release of sex steroids and final maturation induction in goldfish (Carassius auratus) broodstock. For this purpose, goldfish broodstock was divided into four groups and treated with 0.9% NaCl with 20 mg/kg metoclopramide (Met) (C); 10 μg/kg body weight (BW) rcGnRH with 20 mg/kg metoclopramide (rcGn10); 15 μg/kg BW rcGnRH with 20 mg/kg metoclopramide (rcGn15); and 20 μg/kg BW rcGnRH with 20 mg/kg metoclopramide (rcGn20). The capability of the rcGnRH for eliciting biological response was tested in vivo by evaluating the changes of 17β estradiol (E2), testosterone (T), and 17α, 20β-dihydroxy-4-pregnen-3-one (DHP) and the induced spawning. Blood samples were obtained at 0 h, 12 h, and 24 h after injection. The rcGn10, rcGn15, and rcGn20 treatments induced lower E2 concentration, especially 24 h post-injection. T levels were significantly higher in rcGn10, rcGn15, and rcGn20 treatments 12 h post-injection than at 0 h and then decreased at 24 h post-injection. Furthermore, the rcGnRH tested significantly enhanced DHP secretion in rcGn10, rcGn15, and rcGn20 treatments 12 h post-injection before a decline at 24 h post-injection. No significant difference between the sampling times was found in the C treatment for the 3 sex steroids tested. The results also displayed that rcGnRH at 10-20 µg/kg of body weight can trigger spawning with the highest speed and efficiency of spawning at 20 µg/kg. The obtained results represent a possible strategy for enhancing the artificial reproduction and ovulation of broodstock fish by rGnRH and further support the use of recombinant hormones to promote reproduction in aquaculture.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.