From asymmetrical to balanced genomic diversification during rediploidization: Subgenomic evolution in allotetraploid fish

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

Evolutionary mechanisms and practical significance of reproductive success and clonal diversity in unisexual vertebrate polyploids

Unisexual reproduction is generally relevant to polyploidy, and unisexual vertebrates are often considered an evolutionary "dead end" due to the accumulation of deleterious mutations and absence of genetic diversity. However, some unisexual polyploids have developed strategies to avoid genomic decay, and thus provide ideal models to unveil unexplored evolutionary mechanisms, from the reproductive success to clonal diversity creation. This article reviews the evolutionary mechanisms for overcoming meiotic barrier and generating genetic diversity in unisexual vertebrates, and summarizes recent research advancements in the polyploid Carassius complex. Gynogenetic gibel carp (Carassius gibelio) is a unique amphitriploid that has undergone a recurrent autotriploidy and has overcome the bottleneck of triploid sterility via gynogenesis. Recently, an efficient strategy in which ploidy changes, including from amphitriploid to amphitetraploid, then from amphitetraploid to novel amphitriploid, drive unisexual-sexual-unisexual reproduction transition and clonal diversity has been revealed. Based on this new discovery, multigenomic reconstruction biotechnology has been used to breed a novel strain with superior growth and stronger disease resistance. Moreover, a unique reproduction mode that combines both abilities of ameiotic oogenesis and sperm-egg fusion, termed as ameio-fusiongensis, has been discovered, and it provides an efficient approach to synthesize sterile allopolyploids. In order to avoid ecological risks upon escape and protect the sustainable property rights of the aquaculture seed industry, a controllable fertility biotechnology approach for precise breeding is being developed by integrating sterile allopolyploid synthesis and gene-editing techniques. This review provides novel insights into the origin and evolution of unisexual vertebrates and into the attempts being made to exploit new breeding biotechnologies in aquaculture.

Asymmetric and parallel subgenome selection co-shape common carp domestication

BACKGROUND

The common carp (Cyprinus carpio) might best represent the domesticated allopolyploid animals. Although subgenome divergence which is well-known to be a key to allopolyploid domestication has been comprehensively characterized in common carps, the link between genetic architecture underlying agronomic traits and subgenome divergence is unknown in the selective breeding of common carps globally.

RESULTS

We utilized a comprehensive SNP dataset in 13 representative common carp strains worldwide to detect genome-wide genetic variations associated with scale reduction, vibrant skin color, and high growth rate in common carp domestication. We identified numerous novel candidate genes underlie the three agronomically most desirable traits in domesticated common carps, providing potential molecular targets for future genetic improvement in the selective breeding of common carps. We found that independently selective breeding of the same agronomic trait (e.g., fast growing) in common carp domestication could result from completely different genetic variations, indicating the potential advantage of allopolyploid in domestication. We observed that candidate genes associated with scale reduction, vibrant skin color, and/or high growth rate are repeatedly enriched in the immune system, suggesting that domestication of common carps was often accompanied by the disease resistance improvement.

CONCLUSIONS

In common carp domestication, asymmetric subgenome selection is prevalent, while parallel subgenome selection occurs in selective breeding of common carps. This observation is not due to asymmetric gene retention/loss between subgenomes but might be better explained by reduced pleiotropy through transposable element-mediated expression divergence between ohnologs. Our results demonstrate that domestication benefits from polyploidy not only in plants but also in animals.

The Goldfish Genome and Its Utility for Understanding Gene Regulation and Vertebrate Body Morphology

Goldfish, widely viewed as an ornamental fish, is a member of Cyprinidae family and has a very long history in research for both genetics and physiology studies. Among Cyprinidae, the chromosomal locations of orthologs and the amino acid sequences are usually highly conserved. Adult goldfish are 1000 times larger than adult zebrafish (who are in the same family of fishes), which can make it easier to perform several types of experiments compared to their zebrafish cousins. Comparing mutant phenotypes in orthologous genes between goldfish and zebrafish can often be very informative and provide a deeper insight into the gene function than studying the gene in either species alone. Comparative genomics and phenotypic comparisons between goldfish and zebrafish will provide new opportunities for understanding the development and evolution of body forms in the vertebrate lineage.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

PacBio Full-Length Transcriptome of a Tetraploid Sinocyclocheilus multipunctatus Provides Insights into the Evolution of Cavefish

Sinocyclocheilus multipunctatus is a second-class nationally protected wild animal in China. As one of the cavefish, S. multipunctatus has strong adaptability to harsh subterranean environments. In this study, we used PacBio SMRT sequencing technology to generate a first representative full-length transcriptome for S. multipunctatus. Sequence clustering analysis obtained 232,126 full-length transcripts. Among all transcripts, 40,487 were annotated in public databases, while 70,300 microsatellites, 2384 transcription factors, and 16,321 long non-coding RNAs were identified. The phylogenetic tree showed that S. multipunctatus shows a closer relationship to Carassius auratus and Cyprinus carpio, phylogenetically diverging from the common ancestor ~14.74 million years ago (Mya). We also found that between 15.6 and 17.5 Mya, S. multipunctatus also experienced an additional whole-genome duplication (WGD) event, which may have promoted the species evolution of S. multipunctatus. Meanwhile, the overall rates of evolutionary of polyploid S. multipunctatus were significantly higher than those of the other cyprinids, and 220 positively selected genes (PSGs) were identified in two sub-genomes of S. multipunctatus. These PSGs are likely to fulfill critical roles in the process of adapting to diverse cave environments. This study has the potential to facilitate future investigations into the genomic characteristics of S. multipunctatus and provide valuable insights into revealing the evolutionary history of polyploid S. multipunctatus.

New reports on iron related proteins: Molecular characterization of two ferroportin genes in common carp (Cyprinus carpio L.) and its expression pattern

Iron uptake, transport, and storage require the involvement of several proteins, including ferroportin (fpn), the sole known iron efflux transporter. Due to its critical function fpn has been studied, particularly in humans. Here, we characterized the ferroportin gene in common carp (Cyprinus carpio L.) and performed RNA-seq analysis to evaluate its constitutive transcription levels across different tissues. Our results indicate that C. carpio possesses two functional fpns with distinct expression patterns, highlighting the potential for functional divergence and expression differentiation among fpns in this species.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Transposon signatures of allopolyploid genome evolution

Hybridization brings together chromosome sets from two or more distinct progenitor species. Genome duplication associated with hybridization, or allopolyploidy, allows these chromosome sets to persist as distinct subgenomes during subsequent meioses. Here, we present a general method for identifying the subgenomes of a polyploid based on shared ancestry as revealed by the genomic distribution of repetitive elements that were active in the progenitors. This subgenome-enriched transposable element signal is intrinsic to the polyploid, allowing broader applicability than other approaches that depend on the availability of sequenced diploid relatives. We develop the statistical basis of the method, demonstrate its applicability in the well-studied cases of tobacco, cotton, and Brassica napus, and apply it to several cases: allotetraploid cyprinids, allohexaploid false flax, and allooctoploid strawberry. These analyses provide insight into the origins of these polyploids, revise the subgenome identities of strawberry, and provide perspective on subgenome dominance in higher polyploids.

An efficient approach to synthesize sterile allopolyploids through the combined reproduction mode of ameiotic oogenesis and sperm-egg fusion in the polyploid Carassius complex

The association between polyploidy and reproduction transition, which is an intriguing issue in evolutionary genetics, can also be exploited as an approach for genetic improvement in agriculture. Recently, we generated novel amphitriploids (NA3n) by integrating the genomes of the gynogenetic Carassius gibelio and sexual C. auratus, and found gynogenesis was recovered in most NA3n females (NA3n♀I). Here, we discovered a unique reproduction mode, termed ameio-fusiongenesis, which combines the abilities of both ameiotic oogenesis and sperm-egg fusion, in a few NA3n females (NA3n♀II). These females inherited ameiotic oogenesis to produce unreduced eggs from gynogenetic C. gibelio and sperm-egg fusion from sexual C. auratus. Subsequently, we utilized this unique reproduction mode to generate a group of synthetic alloheptaploids by crossing NA3n♀II with Megalobrama amblycephala. They contained all chromosomes of maternal NA3n♀II and a chromosomal set of paternal M. amblycephala. Intergenomic chromosome translocations between NA3n♀II and M. amblycephala were also observed in a few somatic cells. Primary oocytes of the alloheptaploid underwent severe apoptosis owing to incomplete double-strand break repair at prophase I. Although spermatocytes displayed similar chromosome behavior at prophase I, they underwent apoptosis due to chromosome separation failure at metaphase I. Therefore, the alloheptaploid females and males were all sterile. Finally, we established a sustainable clone for the large-scale production of NA3n♀II and developed an efficient approach to synthesize diverse allopolyploids containing genomes of different cyprinid species. These findings not only broaden our understanding of reproduction transition but also offer a practical strategy for polyploidy breeding and heterosis fixing.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Monsoon boosted radiation of the endemic East Asian carps

Major historical events often trigger the rapid flourishing of a few lineages, which in turn shape established biodiversity patterns. How did this process occur and develop? This study provides a window into this issue. The endemic East Asian carps (EEAC) dominated the ichthyofauna of East Asia and exhibited a high degree of adaptation to monsoonal river-lake ecosystems. A series of evidence, including ecogeography, phylogenetics, and macroevolution, suggests that the EEAC is a lineage that arose with the East Asian monsoon and thrived intimately with subsequent monsoon activities. We further deduce the evolution of the EEAC and find that a range of historical events in the monsoon setting (e.g., marine transgression and regression and glacial-interglacial cycle) have further reshaped the distribution patterns of EEAC's members. Comparative genomics analyses reveal that introgressions during the initial period of EEAC radiation and innovations in the regulation of the brain and nervous system may have aided their adaptation to river-lake ecosystems in a monsoon setting, which boosted radiation. Overall, this study strengthens knowledge of the evolutionary patterns of freshwater fishes in East Asia and provides a model case for understanding the impact of major historical events on the evolution of biota.

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.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Genome-wide identification, evolutionary analysis, and antimicrobial activity prediction of CC chemokines in allotetraploid common carp, Cyprinus carpio

Chemokines are a group of secreted small molecules which are essential for cell migration in physiological and pathological conditions by binding to specific chemokine receptors. They are structurally classified into five groups, namely CXC, CC, CX3C, XC and CX. CC chemokine group is the largest one among them. In this study, we identified and characterized 61 CC chemokines from allotetraploid common carp (Cyprinus carpio). The sequence analyses showed that the majority of CC chemokines had an N-terminal signal peptide, and an SCY domain, and all CC chemokines were located in the extracellular region. Phylogenetic, evolutionary and syntenic analyses confirmed that CC chemokines were annotated as 11 different types (CCL19, CCL20, CCL25, CCL27, CCL32, CCL33, CCL34, CCL35, CCL36, CCL39, and CCL44), which exhibited unique gene arrangement pattern and chromosomal location respectively. Furthermore, genome synteny analyses between common carp and four representative teleost species indicated expansion of common carp CC chemokines resulted from the whole genome duplication (WGD) event. Additionally, the continuous evolution of gene CCL25s in teleost afforded a novel viewpoint to explain the WGD event in teleost. Then, we predicted the three-dimensional structures and probable function regions of common carp CC chemokines. All the CC chemokines core structures were constituted of an N-loop, a three-stranded β-sheet, and a C-terminal helix. Finally, 43 CC chemokines were predicted to have probable general antimicrobial activity. Their tertiary structures, cationic and amphiphilic physicochemical property supported the viewpoint. To verify the prediction, six recombinant CCL19s proteins were prepared and the antibacterial activity against Escherichia coli and Aeromonas hydrophila were verified. The results supported our prediction that rCCL19a.1s (rCCL19a.1_a, rCCL19a.1_b) and rCCL19bs (rCCL19b_a, rCCL19b_b), especially rCCL19bs, exhibited extremely significant inhibition to the growth of both E. coli and A. hydrophila. On the contrary, two rCCL19a.2s had no significant inhibitory effect. These studies suggested that CC chemokines were essential in immune system evolution and not monofunctional during pathogen infection.

Immunoglobulin heavy-chain loci in ancient allotetraploid goldfish

As an ancient allotetraploid species, goldfish (Carassius auratus) have two sets of subgenomes. In this study, immunoglobulin heavy-chain (IGH) genes were cloned from the red crucian carp (Carassius auratus red var.), and the corresponding loci were identified in the gynogenetic diploid red crucian carp (GRCC) genome as well as the genomes of three other goldfish strains (Wakin, G-12, and CaTCV-1). Examination showed that each goldfish strain possessed two sets of parallel IGH loci: a complete IGHA locus and a degenerated IGHB locus that was nearly 40 × smaller. In the IGHA locus, multiple τ chain loci were arranged in tandem between the μ&δ chain locus and the variable genes, but no τ-like genes were found in the IGHB locus.

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

Background

Interspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.

Results

We firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F₁ and allotetraploid (F₂₂ and F₂₄) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.

Conclusions

Our results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01401-4.

Correlations of expression of nuclear and mitochondrial genes in triploid fish

The expression of nuclear and mitochondrial genes, as well as their coordinated control, regulates cell proliferation, individual development, and disease in animals. However, the potential coregulation between nuclear and mitochondrial genes is unclear in triploid fishes. The two triploids (R2C and RC2) with distinct mitochondrial genomes but similar nuclear genomes exhibit different embryonic development times and growth rates. They are an excellent model for studying how nuclear and mitochondrial genes coordinate. Here, we performed the mRNA-seq of four stages of embryonic development (blastula, gastrula, segmentation, and hatching periods) in the two triploids (R2C and RC2) and their diploid inbred parents (red crucian carp and common carp). After establishing the four patterns of mitochondrial and nuclear gene expression, 270 nuclear genes regulated by mitochondrial genes were predicted. The expression levels of APC16 and Trim33 were higher in RC2 than in R2C, suggesting their potential effects on regulating embryonic development time. In addition, 308 differentially expressed genes filtered from the list of nuclear-encoded mitochondrial genes described by Mercer et al. in 2011 were considered potential genes for which nuclear genes regulate mitochondrial function. The findings might aid in our understanding of the correlation between mitochondrial and nuclear genomes as well as their synergistic effects on embryonic development.

Comparative transcriptome and adaptive evolution analysis on the main liver and attaching liver of Pareuchiloglanis macrotrema

Pareuchiloglanis macrotrema is a glyptosternoid fish belonging to the Siluriform family and is endemic to the Qinghai-Tibet Plateau tributaries. P. macrotrema is an ideal model for studying the adaptive evolution of fish at high altitudes. P. macrotrema has two attaching livers connected to the main liver, a common feature in most Sisoridae fishes but is a special phenomenon relative to other vertebrates. Using RNA-Seq, 42 differentially expressed genes were found between the main liver and attaching liver, of which 31 were upregulated and 11 were downregulated in the main liver. The major differentially expressed genes between the main liver and attaching liver of P. macrotrema are related to metabolism, immunity, and digestive processes. Meanwhile, a comparative transcriptome analysis was carried out on P. macrotrema fish and six non-plateau Siluriformes fishes. We found 268 positively selected genes in P. macrotrema that are related to energy metabolism, immunity, and hypoxic responses. The findings of this study highlight the gene expression differences between the main liver and attaching livers of Sisoridae fishes and provide greater insight into the evolution of Tibetan fishes.

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.

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.

The neuropeptide Y receptor gene repository, phylogeny and comparative expression in allotetraploid common carp

NPY-family receptors belong to G protein-coupled receptors (GPCR), which lays a physiological foundation for the transmembrane transport of an endogenous appetite-stimulating factor neuropeptide Y and related peptides. In this study, we investigated the npyr genes in ten representative species, and twelve npyr genes were identified from allotetraploid C. carpio, the npyr gene number of C. carpio was twice the number of its subgenome B progenitor-like diploid Poropuntius huangchuchieni. Phylogenetic analysis showed that all npyr genes were divided into three subgroups, and they underwent strong purifying selection according to selection pressure analysis. Subsequently, synteny analysis showed that most npyr genes were evenly distributed on the homologous chromosomes of two subgenomes in allotetraploid C. carpio, in which npy1r and npy2r were tandem duplicated, respectively. In addition, the global expression of npyr genes during embryonic development in allotetraploid C. carpio suggested the potential function of npyr genes in immunity and reproduction. In adult tissues, npyr genes were mainly distributed in the brain, gonad, and skin, which displayed a similar expression pattern between the C. carpio B subgenome and P. huangchuchieni. In general, our research could provide reference information for future exploration of the NPY receptors and neuroendocrine system of allotetraploid C. carpio and vertebrates.

Heterosis of growth trait regulated by DNA methylation and miRNA in allotriploid fish

Background

Heterosis of growth traits in allotriploid fish has benefited the production of aquaculture for many years, yet its genetic and molecular basis has remained obscure. Now, an allotriploid complex, including two triploids and their diploid inbred parents, has provided an excellent model for investigating the potential regulatory mechanisms of heterosis.

Results

Here, we performed a series of analyses on DNA methylation modification and miRNA expression in combination with gene expression in the allotriploid complex. We first established a model of cis- and trans-regulation related to DNA methylation and miRNA in allotriploids. Then, comparative analyses showed that DNA methylation contributed to the emergence of a dosage compensation effect, which reduced gene expression levels in the triploid to the diploid state. We detected 31 genes regulated by DNA methylation in the subgenomes of the allotriploids. Finally, the patterns of coevolution between small RNAs and their homoeologous targets were classified and used to predict the regulation of miRNA expression in the allotriploids.

Conclusions

Our results uncovered the regulatory network between DNA methylation and miRNAs in allotriploids, which not only helps us understand the regulatory mechanisms of heterosis of growth traits but also benefits the study and application of epigenetics in aquaculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00455-6.

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.

Parallel subgenome structure and divergent expression evolution of allo-tetraploid common carp and goldfish

How two subgenomes in allo-tetraploids adapt to coexistence and coordinate through structure and expression evolution requires extensive studies. In the present study, we report an improved genome assembly of allo-tetraploid common carp, an updated genome annotation of allo-tetraploid goldfish and the chromosome-scale assemblies of a progenitor-like diploid Puntius tetrazona and an outgroup diploid Paracanthobrama guichenoti. Parallel subgenome structure evolution in the allo-tetraploids was featured with equivalent chromosome components, higher protein identities, similar transposon divergence and contents, homoeologous exchanges, better synteny level, strong sequence compensation and symmetric purifying selection. Furthermore, we observed subgenome expression divergence processes in the allo-tetraploids, including inter-/intrasubgenome trans-splicing events, expression dominance, decreased expression levels, dosage compensation, stronger expression correlation, dynamic functionalization and balancing of differential expression. The potential disorders introduced by different progenitors in the allo-tetraploids were hypothesized to be alleviated by increasing structural homogeneity and performing versatile expression processes. Resequencing three common carp strains revealed two major ecotypes and uncovered candidate genes relevant to growth and survival rate.

Sequencing an F1 hybrid of Silurus asotus and S. meridionalis enabled the assembly of high-quality parental genomes

Genome complexity such as heterozygosity may heavily influence its de novo assembly. Sequencing somatic cells of the F1 hybrids harboring two sets of genetic materials from both of the paternal and maternal species may avoid alleles discrimination during assembly. However, the feasibility of this strategy needs further assessments. We sequenced and assembled the genome of an F1 hybrid between Silurus asotus and S. meridionalis using the SequelII platform and Hi-C scaffolding technologies. More than 300 Gb raw data were generated, and the final assembly obtained 2344 scaffolds composed of 3017 contigs. The N50 length of scaffolds and contigs was 28.55 Mb and 7.49 Mb, respectively. Based on the mapping results of short reads generated for the paternal and maternal species, each of the 29 chromosomes originating from S. asotus and S. meridionalis was recognized. We recovered nearly 94% and 96% of the total length of S. asotus and S. meridionalis. BUSCO assessments and mapping analyses suggested that both genomes had high completeness and accuracy. Further analyses demonstrated the high collinearity between S. asotus, S. meridionalis, and the related Pelteobagrus fulvidraco. Comparison of the two genomes with that assembled only using the short reads from non-hybrid parental species detected a small portion of sequences that may be incorrectly assigned to the different species. We supposed that at least part of these situations may have resulted from mitotic recombination. The strategy of sequencing the F1 hybrid genome can recover the vast majority of the parental genomes and may improve the assembly of complex genomes.

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.

Chromosome-level genome of Poropuntius huangchuchieni provides a diploid progenitor-like reference genome for the allotetraploid Cyprinus carpio

The diploid Poropuntius huangchuchieni in the cyprinid family, which is widely distributed in the Mekong and Red River basins, is one of the most closely related diploid progenitor-like species of allotetraploid common carp, which was generated by merging of two diploid genomes during evolution. Therefore, the P. huangchuchieni genome is essential for polyploid evolution studies in Cyprinidae. Here, we report a high-quality chromosome-level genome assembly of P. huangchuchieni by integrating Oxford Nanopore and Hi-C technologies. The assembled genome size was 1,021.38 Mb, 895.66 Mb of which was anchored onto 25 chromosomes with a N50 of 32.93 Mb. The genome contained 486.28 Mb repetitive elements and 24,099 protein-coding genes. Approximately 95.9% of the complete BUSCOs were detected, suggesting a high completeness of the genome. Evolutionary analysis revealed that P. huangchuchieni diverged from Cyprinus carpio at approximately 12 Mya. Genome comparison between P. huangchuchieni and the B subgenome of C. carpio provided insights into chromosomal rearrangements during the allotetraploid speciation. With the complete gene set, 17,474 orthologous genes were identified between P. huangchuchieni and C. carpio, providing a broad view of the gene component in the allotetraploid genome, which is critical for future genetic analyses. The high-quality genomic data set created for P. huangchuchieni provides a diploid progenitor-like reference for the evolution and adaptation of allotetraploid carps.

Comparative analyses of hypothalamus transcriptomes reveal fertility-, growth-, and immune-related genes and signal pathways in different ploidy cyprinid fish

Triploid crucian carp (TCC) is obtained by hybridization of female diploid red crucian carp (Carassius auratus red var., RCC) and male allotetraploid hybrids. In this study, high-throughput sequencing was used to conduct the transcriptome analysis of the female hypothalamus of diploid RCC, diploid common carp (Cyprinus carpio L., CC) and TCC. The key functional expression genes of the hypothalamus were obtained through functional gene annotation and differential gene expression screening. A total of 71.56 G data and 47,572 genes were obtained through sequencing and genome mapping, respectively. The Fuzzy Analysis Clustering assigned the differentially expressed genes (DEGs) into eight groups, two of which, overdominance expression (6005, 12.62%) and underdominance expression (3849, 8.09%) in TCC were further studied. KEGG enrichment analysis showed that the DEGs in overdominance were mainly enriched in four pathways. The expression of several fertility-related genes was lower levels in TCC, whereas the expression of several growth-related genes and immune-related genes was higher levels in TCC. Besides, 15 DEGs were verified by quantitative real-time PCR (qPCR). The present study can provide a reference for breeding sterility, fast-growth, and disease-resistant varieties by distant hybridization.

Genomic mosaicism due to homoeologous exchange generates extensive phenotypic diversity in nascent allopolyploids

Allopolyploidy is an important process in plant speciation, yet newly formed allopolyploid species typically suffer from extreme genetic bottlenecks. One escape from this impasse might be homoeologous meiotic pairing, during which homoeologous exchanges (HEs) generate phenotypically variable progeny. However, the immediate genome-wide patterns and resulting phenotypic diversity generated by HEs remain largely unknown. Here, we analyzed the genome composition of 202 phenotyped euploid segmental allopolyploid individuals from the fourth selfed generation following chromosomal doubling of reciprocal F1 hybrids of crosses between rice subspecies, using whole-genome sequencing. We describe rampant occurrence of HEs that, by overcoming incompatibility or conferring superiority of hetero-cytonuclear interactions, generate extensive and individualized genomic mosaicism across the analyzed tetraploids. We show that the resulting homoeolog copy number alteration in tetraploids affects known-function genes and their complex genetic interactions, in the process creating extraordinary phenotypic diversity at the population level following a single initial hybridization. Our results illuminate the immediate genomic landscapes possible in a tetraploid genomic environment, and underscore HE as an important mechanism that fuels rapid phenotypic diversification accompanying the initial stages of allopolyploid evolution.

The evolutionary origin and domestication history of goldfish (Carassius auratus)

Goldfish have been subjected to over 1,000 y of intensive domestication and selective breeding. In this report, we describe a high-quality goldfish genome (2n = 100), anchoring 95.75% of contigs into 50 pseudochromosomes. Comparative genomics enabled us to disentangle the two subgenomes that resulted from an ancient hybridization event. Resequencing 185 representative goldfish variants and 16 wild crucian carp revealed the origin of goldfish and identified genomic regions that have been shaped by selective sweeps linked to its domestication. Our comprehensive collection of goldfish varieties enabled us to associate genetic variations with a number of well-known anatomical features, including features that distinguish traditional goldfish clades. Additionally, we identified a tyrosine-protein kinase receptor as a candidate causal gene for the first well-known case of Mendelian inheritance in goldfish-the transparent mutant. The goldfish genome and diversity data offer unique resources to make goldfish a promising model for functional genomics, as well as domestication.