Genomic incompatibilities in the diploid and tetraploid offspring of the goldfish × common carp cross

The key role of myostatin b in somatic growth in fishes derived from distant hybridization

The basic mechanism of heterosis has not been systematically and completely characterized. In previous studies, we obtained three economically important fishes that exhibit rapid growth, WR (WCC ♀ × RCC ♂), WR-II (WR ♀ × WCC ♂), and WR-III (WR-II ♀ × 4nAU ♂), through distant hybridization. However, the mechanism underlying this rapid growth remains unclear. In this study, we found that WR, WR-II, and WR-III showed muscle hypertrophy and higher muscle protein and fat contents compared with their parent species (RCC and WCC). Candidate genes responsible for this rapid growth were then obtained through an analysis of 12 muscle transcriptomes. Notably, the mRNA level of mstnb (myostatin b), which is a negative regulator of myogenesis, was significantly reduced in WR, WR-II, and WR-III compared with the parent species. To verify the function of mstnb, a mstnb-deficient mutant RCC line was generated using the CRISPR-Cas9 technique. The average body weight of mstnb-deficient RCC at 12 months of age was significantly increased by 29.57% compared with that in wild-type siblings. Moreover, the area and number of muscle fibers were significantly increased in mstnb-deficient RCC, indicating hypertrophy and hyperplasia. Furthermore, the muscle protein and fat contents were significantly increased in mstnb-deficient RCC. The molecular regulatory mechanism of mstnb was then revealed by transcription profiling, which showed that genes related to myogenesis (myod, myog, and myf5), protein synthesis (PI3K-AKT-mTOR), and lipogenesis (pparγ and fabp3) were highly activated in hybrid fishes and mstnb-deficient RCC. This study revealed that low expression or deficiency of mstnb regulates somatic growth by promoting myogenesis, protein synthesis, and lipogenesis in hybrid fishes and mstnb-deficient RCC, which provides evidence for the molecular mechanism of heterosis via distant hybridization.

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.

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.

Host Hybridization Dominates over Cohabitation in Affecting Gut Microbiota of Intrageneric Hybrid Takifugu Pufferfish

Microbial symbionts are of great importance for macroscopic life, including fish, and both collectively comprise an integrated biological entity known as the holobiont. Yet little is known as to how the normal balance within the fish holobiont is maintained and how it responds to biotic and/or abiotic influences. Here, through amplicon profiling, the genealogical relationship between artificial F1 hybrid pufferfish with growth heterosis, produced from crossing female Takifugu obscurus with male Takifugu rubripes and its maternal halfsibling purebred, was well recapitulated by their gut microbial community similarities, indicating an evident parallelism between host phylogeny (hybridity) and microbiota relationships therein. Interestingly, modest yet significant fish growth promotion and gut microbiota alteration mediated by hybrid-purebred cohabitation were observed, in comparison with their respective monoculture cohorts that share common genetic makeups, implying a certain degree of environmental influences. Moreover, the underlying assemblage patterns of gut microbial communities were found associated with a trade-off between variable selection and dispersal limitation, which are plausibly driven by the augmented social interactions between hybrid and purebred cohabitants differing in behaviors. Results from this study not only can enrich, from a microbial perspective, the sophisticated understanding of complex and dynamic assemblage of the fish holobiont, but will also provide deeper insights into the ecophysiological factors imposed on the diversity-function relationships thereof. Our findings emphasize the intimate associations of gut microbiota in host genetics-environmental interactions and would have deeper practical implications for microbial contributions to optimize performance prediction and to improve the production of farmed fishes. IMPORTANCE Microbial symbionts are of great importance for macroscopic life, including fish, and yet little is known as to how the normal balance within the fish holobiont is maintained and how it responds to the biotic and/or abiotic influences. Through gut microbiota profiling, we show that host intrageneric hybridization and cohabitation can impose a strong disturbance upon pufferfish gut microbiota. Moreover, marked alterations in the composition and function of gut microbiota in both hybrid and purebred pufferfish cohabitants were observed, which are potentially correlated with different metabolic priorities and behaviors between host genealogy. These results can enrich, from a microbial perspective, the sophisticated understanding of the complex and dynamic assemblage of the fish holobiont and would have deeper practical implications for microbial contributions to optimize performance prediction and to improve farmed fish production.

Effects of Sheep Sires on Muscle Fiber Characteristics, Fatty Acid Composition and Volatile Flavor Compounds in F1 Crossbred Lambs

Crossbreeding significantly improves meat production performance in sheep; however, whether hybridization changes the meat quality characteristics of lambs is uncertain. We analyzed the effects of three different hybrid sires on muscle fiber characteristics (MFCs), fatty acid composition (FAC), and volatile flavor compounds (VFCs) in lambs under identical feeding conditions. Compared with those of purebred lambs, the muscle fiber diameter and cross-sectional areas of the crossbred lambs were significantly decreased (p < 0.05), and the collagen fiber content was significantly increased (p < 0.05). The numbers and area ratios of the fast and slow muscle fibers did not significantly differ between the purebred and crossbred lambs, but the expressions of four MyHC gene types differed significantly (p < 0.05). Twenty-three fatty acids were identified in both the purebred and crossbred lambs, of which thirteen were differentially expressed (p < 0.05). Saturated fatty acid (SFA) contents in the crossbred lambs were significantly increased (p < 0.05), whereas the monounsaturated fatty acid content was significantly decreased (p < 0.05). Polyunsaturated fatty acid/SFA and n-6/n-3 ratios were significantly lower in the crossbred lambs than in the purebred lambs (p < 0.05). Twenty-five VFCs were identified among the three hybrids, and aldehydes were the main VFCs. Eleven VFCs were differentially expressed in the crossbred lambs (p < 0.05). Hybrid sires affected the MFCs, FAC, and VFCs of the F1 lambs, thus providing a reference for high-quality mutton production.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10^-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.

Breaking the reproductive barrier of divergent species to explore the genomic landscape

Background: Climate change will have significant consequences for species. Species range shifts induce the emergence of new hybrid zones or the spatial displacement of pre-existing ones. These hybrid zones may become more porous as alleles are passed from one species to another. Currently, hybridization between highly divergent species living in sympatry seems extremely limited. Indeed, this phenomenon involves breaking two barriers. The first is the pre-mating barrier, related to the reproductive phenology of the two species. The second is the post-zygotic barrier, related to the genetic divergence between these species. Here, we were interested in identifying new hybridization patterns and potential implications, especially in the context of environmental modifications.

Methods: We sampled Telestes souffia and Parachondrostoma toxostoma wild specimens from different locations across France and genotyped them for SNP markers. We identified discriminant loci using F1-hybrid specimens and parental species and performed principal component analysis and Bayesian model-based clustering to analyze phylogenetic information. Furthermore, we assessed deviation in allele frequency from F1 to F2 and for Hardy–Weinberg equilibrium for F2 and assessed gene function associated with two F2 cohorts.

Results: We demonstrate that by breaking the ecological barrier, massive introgressive hybridization is possible between two endemic lineages of Cyprinidae belonging to two distinct genera. For both cohorts studied (=2 cm and >2 cm), a large majority of loci (>88%) presented no deviation in allele frequency and no departure from the Hardy–Weinberg equilibrium. For individuals beyond the 2 cm stage, two phenomena were observed. The first was an allelic imbalance in favor of P. toxostoma, for some genomic regions, with genes involved in developmental regulatory processes, cytoskeletal organization, and chromosome organization. The second was an excess of heterozygous loci coupled with an equilibrium of allelic frequencies for genes involved in immune response and kidney/liver development. Moreover, the 2 cm-sized specimens with high mortality yielded a particular genomic signature.

Conclusion: Our study displayed important results for understanding the early stages of hybridization between divergent lineages and predicting the emergence of future hybrid zones in the wild. Moreover, this hybridization generates a wide spectrum of hybrids that are a potential source of important evolutionary novelties.

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 analysis on the immunoregulatory roles of ferritin M in hybrid fish (Carassius cuvieri ♀ × Carassius auratus red var ♂) and its parental species after bacterial infection

Ferritin M is involved in the regulation of fish immunity. In this study, open reading frame (ORF) sequences of ferritin M from hybrid fish and its parental species were 534 bp. Tissue-specific analysis indicated that the highest level of ferritin M from red crucian carp was observed in kidney, while peaked expressions of ferritin M from white crucian carp and hybrid carp were observed in gill. Elevated levels of ferritin M from hybrid carp and its parental species were detected in immune-related tissues following Aeromonas hydrophila infection or in cultured fish cell lines after lipopolysaccharide (LPS) challenge. Ferritin M overexpression could attenuate NF-κB and TNFα promoter activity in their respective fish cells. Purified ferritin M fusion proteins elicited in vitro binding activity to A. hydrophila and Edwardsiella tarda, lowered bacterial dissemination to tissues and alleviated inflammatory response. Furthermore, treatment with ferritin M fusion proteins could mitigate bacteria-induced liver damage and rescue antioxidant activity. These results suggested that ferritin M in hybrid fish showed a similar immune defense against bacteria infection in comparison with those of its parental species.

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.

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.

Research Progress of the Gut Microbiome in Hybrid Fish

Fish, including hybrid species, are essential components of aquaculture, and the gut microbiome plays a vital role in fish growth, behavior, digestion, and immune health. The gut microbiome can be affected by various internal and/or external factors, such as host development, diet, and environment. We reviewed the effects of diet and dietary supplements on intestinal microorganisms in hybrid fish and the difference in the gut microbiome between the hybrid and their hybrids that originate. Then, we summarized the role of the gut microbiome in the speciation and ecological invasion of hybrid fish. Finally, we discussed possible future studies on the gut microbiome in hybrid fish, including the potential interaction with environmental microbiomes, the effects of the gut microbiome on population expansion, and fish conservation and management.

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.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

Blood cell characterization and transcriptome analysis reveal distinct immune response and host resistance of different ploidy cyprinid fish following Aeromonas hydrophila infection

Aeromonas hydrophila can pose a great threat to survival of freshwater fish. In this study, A. hydrophila infection could decrease blood cell numbers, promote blood cell damage as well as alter the levels of alkaline phosphatase (ALP), lysozyme (LZM), aspartate aminotransferase (AST), total antioxidant capacity (T-AOC), total superoxide dismutase (SOD), catalase (CAT) and malondialdehyde (MDA) in immune-related tissues of red crucian carp (RCC, 2 N = 100) and triploid cyprinid fish (3 N fish, 3 N = 150). In addition, the significant alternation of antioxidant status was observed in PBMCs isolated from RCC and 3 N following LPS stimulation. The core differential expression genes (DEGs) involved in apoptosis, immunity, inflammation and cellular signals were co-expressed differentially in RCC and 3 N following A. hydrophila challenge. NOD-like receptor (NLR) signals appeared to play a critical role in A. hydrophila-infected fish. DEGs of NLR signals in RCCah vs RCCctl were enriched in caspase-1-dependent Interleukin-1β (IL-1β) secretion, interferon (IFN) signals as well as cytokine activation, while DEGs of NLR signals in 3Nah vs 3Nctl were enriched in caspase-1-dependent IL-1β secretion and antibacterial autophagy. These results highlighted the differential signal regulation of different ploidy cyprinid fish to cope with bacterial infection.

Characteristics of SP600125 Induced Tetraploid Cells in Comparison With Diploid and Tetraploid Cells of Fish

In our previous research, SP600125 (Anthrapyrazolone) was used to induce autotetraploid of crucian carp cells (SP4N cells), and tetraploid fry was generated from the SP4N cells by somatic cell nuclear transfer technique. However, it is still unclear about biological characteristics of the SP4N cells. In this article, the cytological characteristic and gene expression profiles of the SP4N cells are investigated in comparison with the crucian carp cells (2N cells) and the tetraploid crucian carp cells (CC4N cells). The SP4N cells have tetraploid characteristics in terms of morphology and DNA ploidy levels, and their chromosome behavior is stable during the cell proliferation. The migration ability and the mtDNA copy number of SP4N cells are both lower than those in the CC4N cells and the 2N cells, but there exist giant mitochondria in the SP4N cells. The similar expression trends in the cell cycle regulation genes of the SP4N cells and 2N cells, while the corresponding expression profiles are clearly different between the SP4N cells and the CC4N cells. Moreover, the significant difference genes are associated with energy metabolism pathways among the SP4N cells, 2N cells and CC4N cells. These results can provide deeper understanding of SP600125 induction, as well as finding applications in polyploidization breeding of fish species.

Ferritin H can counteract inflammatory response in hybrid fish and its parental species after Aeromonas hydrophila infection

Ferritin H can participate in the regulation of fish immunity. Tissue-specific analysis revealed that the highest expressions of Ferritin H in parental species were observed in spleen, while peaked level of Ferritin H mRNA in hybrid fish was observed in liver. In addition, A. hydrophila challenge could sharply enhance their Ferritin H mRNA expression in liver, kidney and spleen. To further investigate their roles in immune regulation, their Ferritin H fusion proteins were produced in vitro. Ferritin H fusion proteins could exhibit a direct binding activity to A. hydrophila and endotoxin in a dose-dependent manner, restrict dissemination of A. hydrophila to tissues and abrogate inflammatory cascades. Moreover, treatment with Ferritin H fusion proteins could reduce A. hydrophila-induced lipid peroxidation. These results indicated that Ferritin H in hybrid fish elicited a similar immune regulation of A. hydrophila-induced inflammatory signals in comparison with those of its parents.

Interploidy gene flow involving the sexual-asexual cycle facilitates the diversification of gynogenetic triploid Carassius fish

Asexual vertebrates are rare and at risk of extinction due to their restricted adaptability through the loss of genetic recombination. We explore the mechanisms behind the generation and maintenance of genetic diversity in triploid asexual (gynogenetic) Carassius auratus fish, which is widespread in East Asian fresh waters and exhibits one of the most extensive distribution among asexual vertebrates despite its dependence on host sperm. Our analyses of genetic composition using dozens of genetic markers and genome-wide transcriptome sequencing uncover admixed genetic composition of Japanese asexual triploid Carassius consisting of both the diverged Japanese and Eurasian alleles, suggesting the involvement of Eurasian lineages in its origin. However, coexisting sexual diploid relatives and asexual triploids in Japan show regional genetic similarity in both mitochondrial and nuclear markers. These results are attributed to a unique unidirectional gene flow from diploids to sympatric triploids, with the involvement of occasional sexual reproduction. Additionally, the asexual triploid shows a weaker population structure than the sexual diploid, and multiple triploid lineages coexist in most Japanese rivers. The generated diversity via repeated interploidy gene flow as well as an increased establishment of immigrants is assumed to offset the cost of asexual reproduction and might contribute to the successful broad distribution of this asexual vertebrate.

Comparative analysis of erythrocyte hemolysis, plasma parameters and metabolic features in red crucian carp (Carassius auratus red var) and triploid hybrid fish following Aeromonas hydrophila challenge

Aeromonas hydrophila can pose a great threat to survival of freshwater fish. In this study, A. hydrophila challenge could promote the erythrocyte hemolysis, increase free hemoglobin (FHB) level and generate malondialdehyde (MDA) production in plasma but decrease the levels of total antioxidant capacity (T-AOC), total superoxide dismutase (SOD), catalase (CAT), alkaline phosphatase (ALP) and lysozyme (LZM) of red crucian carp (RCC, 2 N = 100) and triploid hybrid fish (3 N fish, 3 N = 150) following A. hydrophila challenge. Elevated expression levels of heat shock protein 90 alpha (HSP90α), matrix metalloproteinase 9 (MMP-9), free fatty acid receptor 3 (FFAR3), paraoxonase 2 (PON2) and cytosolic phospholipase A2 (cPLA2) were observed in A. hydrophila-infected fish. In addition, A. hydrophila challenge could significantly increase expressions of cortisol, leucine, isoleucine, glutamate and polyunsaturated fatty acids (PUFAs) in RCC and 3 N, while glycolysis and tricarboxylic acid cycle appeared to be inactive. We identified differential fatty acid derivatives and their metabolic networks as crucial biomarkers from metabolic profiles of different ploidy cyprinid fish subjected to A. hydrophila infection. These results highlighted the comparative metabolic strategy of different ploidy cyprinid fish against bacterial infection.

Effect of Lipopolysaccharide (LPS) stimulation on apoptotic process and oxidative stress in fibroblast cell of hybrid crucian carp compared with those of Carassius cuvieri and Carassius auratus red var

Bacterial LPS is a heat-stable endotoxin and wall components of gram negative bacteria, which can exhibit a toxicological effect on physiology and biochemical activities of fish. In this study, we investigated the effect of LPS exposure on cell viability, oxidative stress, caspase activity and immune-related gene expressions in cultured fin cell lines of red crucian carp, white crucian carp and their hybrid offspring. LPS stimulation could reduce fish cell viability, whereas gene expression levels and promoter activities in inflammatory signals increased dramatically. Moreover, enhanced levels of intracellular oxidative stress and decreased levels of mitochondrial membrane potential (MMP) were observed in LPS-induced fish cells. N-Acetyl-L-cysteine (NAC) could alleviate LPS-stimulated reactive oxygen species (ROS) generation and caspase-3 activity in fish cells. These results suggested that ROS-mediated cytotoxic stress was involved in LPS-induced inflammation and mitochondrial damage in cultured fish cells.

Localization of RNA Pol II CTD (S5) and Transcriptome Analysis of Testis in Diploid and Tetraploid Hybrids of Red Crucian Carp (♀) × Common Carp (♂)

Polyploidy occurs naturally in fish; however, the appearance of these species is an occasional and gradual process, which makes it difficult to trace the changes in phenotypes, genotypes, and regulation of gene expression. The allotetraploid hybrids (4nAT) of red crucian carp (RCC; ♀) × common carp (CC; ♂) generated from interspecies crossing are a good model to investigate the initial changes after allopolyploidization. In the present study, we focused on the changes in the active sites of the testicular transcriptome of the allotetraploid by localization of RNA Pol II CTD YSPTSPS (phospho S5) using immunofluorescence and RNA-seq data via bioinformatic analysis. The results showed that there was no significant difference in signal counts of the RNA Pol II CTD (S5) between the different types of fish at the same stages, including RCC, CC, 2nF₁, and 4nAT, which means that the number of transcriptionally active sites on germ cell chromosomes was not affected by the increase in chromosome number. Similarly, RNA-seq analysis indicated that in the levels of chromosomes and 10-kb regions in the genome, there were no significant changes in the highly active sites in RCC, 2nF₁, and 4nAT. These findings suggest that at the beginning of tetraploid origin, the active transcriptome site of 4nAT in the testis was conserved in the regions of the genome compared to that in RCC and 2nF₁. In conclusion, 4nAT shared a similar gene expression model in the regions of the genome with RCC and 2nF₁ with significantly different expression levels.

Diverse transcriptional patterns of homoeologous recombinant transcripts in triploid fish (Cyprinidae)

Homoeologous recombination (HR), the exchange of homoeologous chromosomes, contributes to subgenome adaptation to diverse environments by producing various phenotypes. However, the potential relevance of HR and innate immunity is rarely described in triploid cyprinid fish species. In our study, two allotriploid genotypes (R₂C and RC₂), whose innate immunity was stronger than their inbred parents (Carassius auratus red var. and Cyprinus carpio L.), were obtained from backcrossing between male allotetraploids of C. auratus red var.×C. carpio L. and females of their two inbred parents, respectively. The work detected 140 HRs shared between the two triploids at the genomic level. Further, transcriptions of 54 homoeologous recombinant genes (HRGs) in R₂C and 65 HRGs in RC₂ were detected using both Illumina and PacBio data. Finally, by comparing expressed recombinant reads to total expressed reads in each of the genes, a range of 0.1%-10% was observed in most of the 99-193 HRGs, of which six recombinant genes were classified as "response to stimulus". These results not only provide a novel way to predict HRs in allopolyploids based on cross prediction at both genomic and transcriptional levels, but also insight into the potential relationship between HRs related to innate immunity and adaptation of the triploids and allotetraploids.

A novel NK-lysin in hybrid crucian carp can exhibit cytotoxic activity in fish cells and confer protection against Aeromonas hydrophila infection in comparison with Carassius cuvieri and Carassius auratus red var

NK-lysin, an effector of natural killer (NK) cells and cytotoxic T lymphocytes (CTLs), not only exhibits cytotoxic effect in fish cells, but also participates in the immune defense against pathogenic infection. In this study, ORF sequences of RCC-NK-lysin, WCC-NK-lysin and WR-NK-lysin were 369 bp. Tissue-specific analysis revealed that the highest expressions of RCC-NK-lysin and WCC-NK-lysin were observed in gill, while the peaked level of WR-NK-lysin mRNA was observed in spleen. A. hydrophila infection sharply increased RCC-NK-lysin, WCC-NK-lysin and WR-NK-lysin mRNA expression in liver, trunk kidney and spleen. In addition, elevated levels of NK-lysin mRNA were observed in cultured fin cell lines of red crucian carp (RCC), white crucian carp (WCC) and their hybrid offspring (WR) after Lipopolysaccharide (LPS) challenge. RCC-NK-lysin, WCC-NK-lysin and WR-NK-lysin exerted regulatory roles in inducing ROS generation, modulating mitochondrial membrane potential, decreasing fish cell viability and antagonizing survival signalings, respectively. RCC/WCC/WR-NK-lysin-overexpressing fish could up-regulate expressions of inflammatory cytokines and decrease bacterial loads in spleen. These results indicated that NK-lysin in hybrid fish contained close sequence similarity to those of its parents, possessing the capacities of cytotoxicity and immune defense against bacterial infection.

Unusual chromosome numbers and polyploidy in invasive fire ant populations

Fire ants (Solenopsis invicta Buren in J Ga Entomol Soc 7:1-26, 1972), an invasive alien ant species, first spread from South America to the United States in the 1930s, the southern part of the United States by the end of the twentieth century, Oceania, Taiwan, and China in the twenty-first century, and finally to Japan and South Korea in 2017. As these ants have significant negative economic, human health, and environmental impacts, the purpose of this research was to accumulate cytogenetic information regarding fire ants and provide basic data for developing management strategies for their control. Fire ants were collected from invasive populations from Taiwan, Florida (USA), and Buenos Aires (Argentina), and a native population from Puerto Iguazu (Argentina), their point of origination, and analyzed with regard to chromosome number, morphology, and polyploidy, silver-stained nucleolar organizer regions (Ag-NORs), and 18S rDNA and telomere fluorescence in situ hybridization (FISH). The results showed that (1) fire ants from invaded populations differed in chromosome morphology compared to those from native populations; (2) the Florida and Taiwanese fire ant populations evinced greater variability in chromosome numbers and polyploidy variations; (3) the Taiwanese population exhibited significantly increased Ag-NOR signals in interphase cells, with signal number significantly positively correlating with distance from native populations; and (4) substantial diversity of signals was also apparent following 18S rDNA and telomere FISH analyses. Variation in these characteristics were hypothesized to be due to (1) the effect of hybridizations and interbreeding between closely related species or genetically distant populations, and (2) the potential effect of large amounts of insecticides sprayed for pest control.

A New Type of Allodiploid Hybrids Derived From Female Megalobrama amblycephala × Male Gobiocypris rarus

Distant hybridization can combine whole genomes from parent species and result in changes in the phenotypes and genotypes in hybrids. The characteristics of many hybrid fishes with even number of chromosomes have been reported, but the hybrids with odd number chromosomes are rarely reported. Blunt snout bream (Megalobrama amblycephala, BSB, 2n = 48) and rare gudgeon (Gobiocypris rarus, RG, 2n = 50) belong to two different subfamilies and have quite different biological characteristics. In this study, we obtain the hybrids (BR) derived from the inter-subfamily hybridization of female BSB and male RG. We investigate the fertilization rate, hatching rate, morphological traits, chromosomal numbers, DNA content, growth rates, and 5S rDNA in the BR. The results show that the BR is an allodiploid fish with 49 chromosomes, and all the measurable traits are significantly different (p < 0.05) among BR, BSB, and BR. Interestingly, the upper part of the BR body color is similar to BSB (gray), the lower part of the BR body color is similar to RG (light yellow), and the BR inherits a unique light yellow wide longitudinal band from the RG. Furthermore, the BR has a fast growth rate compared with RG. The 5S rDNA of the BR inherits the specific bands of its parental 5S rDNA respectively and has some mutations, which show obvious recombination, heredity, and variability in BR. This study will be of great significance in fish genetic breeding.

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.

Enhanced Immune Response Improves Resistance to Cadmium Stress in Triploid Crucian Carp

Previous research has indicated that triploid crucian carp (3n fish) have preferential resistance to cadmium (Cd) compared to Carassius auratas red var. (2n fish). In this article, comparative research is further conducted between the 2n and 3n fish in terms of the immune response to Cd-induced stress. Exposure to 9 mg/L Cd for 96 h changed the hepatic function indexes remarkably in the 2n fish, but not in the 3n fish. In the serum of Cd-treated 2n fish, the levels of alanine amino transferase, aspartate aminotransferase, adenosine deaminase, and total bilirubin significantly increased, while the levels of total protein, albumin, lysozyme, and anti-superoxide anion radicals decreased demonstrating hepatotoxicity. By analysis of transcriptome profiles, many immune-related pathways were found to be involved in the response of 3n fish to the Cd-induced stress. Expression levels of the immune genes, including the interleukin genes, tumor necrosis factor super family member genes, chemokine gene, toll-like receptor gene, and inflammatory marker cyclooxygenase 2 gene were significantly enhanced in the hepatopancreas of the Cd-treated 3n fish. In contrast, the expression levels of these genes decreased in the 2n fish. This research provides a theoretical basis for polyploid fish breeding and is helpful for the ecological restoration of water due to pollution.

Asymmetric expression of homoeologous genes contributes to dietary adaption of an allodiploid hybrid fish derived from Megalobrama amblycephala (♀) × Culter alburnus (♂)

BACKGROUND

Hybridization, which can quickly merge two or more divergent genomes and form new allopolyploids, is an important technique in fish genetic breeding. However, the merged subgenomes must adjust and coexist with one another in a single nucleus, which may cause subgenome interaction and dominance at the gene expression level and has been observed in some allopolyploid plants. In our previous studies, newly formed allodiploid hybrid fish derived from herbivorous Megalobrama amblycephala (♀) × carnivorous Culter alburnus (♂) had herbivorous characteristic. It is thus interesting to further characterize whether the subgenome interaction and dominance derive dietary adaptation of this hybrid fish.

RESULTS

Differential expression, homoeolog expression silencing and bias were investigated in the hybrid fish after 70 days of adaptation to carnivorous and herbivorous diets. A total of 2.65 × 108 clean reads (74.06 Gb) from the liver and intestinal transcriptomes were mapped to the two parent genomes based on specific SNPs. A total of 2538 and 4385 differentially expressed homoeologous genes (DEHs) were identified in the liver and intestinal tissues between the two groups of fish, respectively, and these DEHs were highly enriched in fat digestion and carbon metabolism, amino acid metabolism and steroid biosynthesis. Furthermore, subgenome dominance were observed in tissues, with paternal subgenome was more dominant than maternal subgenome. Moreover, subgenome expression dominance controlled functional pathways in metabolism, disease, cellular processes, environment and genetic information processing during the two dietary adaptation processes. In addition, few but sturdy villi in the intestine, significant fat accumulation and a higher concentration of malondialdehyde in the liver were observed in fish fed carnivorous diet compared with fish fed herbivorous diet.

CONCLUSIONS

Our results indicated that diet drives phenotypic and genetic variation, and the asymmetric expression of homoeologous genes (including differential expression, expression silencing and bias) may play key roles in dietary adaptation of hybrid fish. Subgenome expression dominance may contribute to uncovering the mechanistic basis of heterosis and also provide perspectives for fish genetic breeding and application.

Chimeric ferritin H in hybrid crucian carp exhibits a similar down-regulation in lipopolysaccharide-induced NF-κB inflammatory signal in comparison with Carassius cuvieri and Carassius auratus red var

Ferritin H can participate in the regulation of teleostean immunity. ORF sequences of RCC/WCC/WR-ferritin H were 609 bp, while WR-ferritin H gene possessed chimeric fragments or offspring-specific mutations. In order to elucidate regulation of immune-related signal transduction, three fibroblast-like cell lines derived from caudal fin of red crucian carp (RCC), white crucian carp (WCC) and their hybrid offspring (WR) were characterized and designated as RCCFCs, WCCFCs and WRFCs. A sharp increase of ferritin H mRNA was observed in RCCFCs, WCCFCs and WRFCs following lipopolysaccharide (LPS) challenge. Overexpression of RCC/WCC/WR-ferritin H can decrease MyD88-IRAK4 signal and antagonize NF-κB, TNFα promoter activity in RCCFCs, WCCFCs and WRFCs, respectively. These results indicated that ferritin H in hybrid offspring harbors highly-conserved domains with a close sequence similarity to those of its parents, playing a regulatory role in inflammatory signals.

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.

Species and tissue specific analysis based on quantitative proteomics from allotetraploid and the parents

Allotetraploids play central roles in the field of polyploid breeding of freshwater fishes. The molecular basis underlying distant hybridization and individual differences between allotetraploids and their parents is largely unknown. In this study, we performed quantitative proteomics profiling in gonad and liver tissues between allotetraploid and its parents Red Crucian Carp (♀) and Common Carp (♂) respectively. Thousands of proteins were identified and quantified. Species and tissue specific analysis revealed that a large number of causal proteins are specifically regulated in gonad or liver tissue in allotetraploid and its parents respectively. Subsequently, integrative bioinformatics analyses including functional enrichment, pathway and network analyses were conducted. The results suggested a series of gonad and liver specifically regulated proteins such as LSM3, LSM7, PABPC1B and ALDH3A1, EHHADHB, ACAT2 play crucial roles in reproduction-related and metabolism-related pathways including "DNA replication", "Spliceosome" and "Metabolic pathways", "Biosynthesis of antibiotics". Meanwhile, species specifically regulated proteins such as FMR1, MAO are involved in "RNA transport", "Glycine, serine and threonine metabolism". Herein, we established the first comprehensive proteomics knowledgebase for particular freshwater fishes. It may shed light on the molecular mechanisms underlying polyploidy breeding and individual differences and serve as an indispensable reference for further studies. SIGNIFICANCE: The molecular basis underlying distant hybridization and individual differences between allotetraploid and their parents Red Crucian Carp and Common Carp is largely unknown. Quantitative proteomics profiling integrated with multiple bioinformatics analysis revealed that a large number of causal proteins are specifically expressed in gonad or liver tissue in allotetraploid and its parents. Herein, we established the first comprehensive proteomics knowledgebase for particular freshwater fishes. It may shed light on the molecular mechanisms underlying polyploidy breeding and individual differences and serve as an indispensable reference for further associated studies.

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.

Optimal reference genes for gene expression analysis in polyploid of Cyprinus carpio and Carassius auratus

BACKGROUND

Reference genes are usually stably expressed in various cells and tissues. However, it was reported that the expression of some reference genes may be distinct in different species. In this study, we intend to answer whether the expression of reported traditional reference genes changes or not in the polyploid fish RESULTS: By retrieving the mRNA sequencing data of three different ploidy fish from the NCBI SRA database, we selected 12 candidate reference genes, and examined their expression levels in the 10 tissues and in the four cell lines of three different ploidy fish by real-time PCR. Then, the expression profiles of these 12 candidate reference genes were systematically evaluated by using the software platforms: BestKeeper, NormFinder and geNorm.

CONCLUSION

The 28S ribosomal protein S5 gene (RPS5) and the ribosomal protein S18 gene (RPS18) are the most suitable reference genes for the polyploid of Cyprinus carpio and Carassius auratus, demonstrated by both of the tissues and the cultured cells.

Sox Gene Family Revealed Genetic Variations in Autotetraploid Carassius auratus

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

3nLcn2, a teleost lipocalin 2 that possesses antimicrobial activity and inhibits bacterial infection in triploid crucian carp

Lipocalin 2 (Lcn2) has been identified in mammals, however, the in vivo function of fish Lcn2 is essentially unknown. Triploid crucian carp (3n = 150) of red crucian carp (female, 2n = 100) and allotetraploid (male, 4n = 200) shows better resistance to pathogenic infections. To elucidate the antimicrobial mechanism of triploid crucian carp, we examined the function of a novel Lcn2 from triploid crucian carp (3nLcn2). 3nLcn2 is 183 residues in length and contains a conserved lipocalin domain. Quantitative real time reverse transcription PCR (qRT-PCR) analysis showed that 3nLcn2 expression occurred in multiple tissues and was upregulated by bacterial infection in a time-dependent manner. We found that purified recombinant 3nLcn2 (r3nLcn2) exerted bactericidal activity to Aeromonas hydrophila and Escherichia coli. qRT-PCR detected increased expression of pro-inflammatory cytokines and tight junctions in fish with 3nLcn2 overexpression. Fish administered with 3nLcn2 exhibited enhanced intestinal barrier and resistance against bacterial infection. These results provide the first evidence that 3nLcn2 is a functional lipocalin with antimicrobial activity and plays a positive role in the immune defense during bacterial infection.

Further evidence for paternal DNA transmission in gynogenetic grass carp

Gynogenesis is an important breeding method in aquaculture and has been widely applied to many fish species. If gynogenetic progenies are to inherit paternal partial genomic DNA, this will increase genetic variation and will provide a useful outcome for breeding. In this study, we investigated the genetic variation in homeobox (Hox) gene clusters (HoxA4a, HoxA9a, HoxA11b, HoxB1b, HoxC4a, HoxC6b, and HoxD10a) among koi carp (Cyprinus carpio haematopterus, KOC; the stimulation sperm source), grass carp (Ctenopharyngodon idellus), and gynogenetic grass carp (GGC). We found paternal DNA (a special DNA fragment and HoxC6b) derived from KOC and a recombinant gene belonging to HoxC6b in GGC. We are the first to report the recombinant HoxC6b in GGC. Our study provides further evidence for paternal DNA transmission to gynogenetic progenies, which is a finding with great significance for the genetic breeding of fish.

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

A persistent enigma is the rarity of polyploidy in animals, compared to its prevalence in plants. Although animal polyploids are thought to experience deleterious genomic chaos during initial polyploidization and subsequent rediploidization processes, this hypothesis has not been tested. We provide an improved reference-quality de novo genome for allotetraploid goldfish whose origin dates to ~15 million years ago. Comprehensive analyses identify changes in subgenomic evolution from asymmetrical oscillation in goldfish and common carp to diverse stabilization and balanced gene expression during continuous rediploidization. The homoeologs are coexpressed in most pathways, and their expression dominance shifts temporally during embryogenesis. Homoeolog expression correlates negatively with alternation of DNA methylation. The results show that allotetraploid cyprinids have a unique strategy for balancing subgenomic stabilization and diversification. Rediploidization process in these fishes provides intriguing insights into genome evolution and function in allopolyploid vertebrates.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Expression Profile Analysis of the Cell Cycle in Diploid and Tetraploid Carassius auratus red var

Polyploidization often leads to “transcriptome shock,” and is considered an important factor in evolution of species. Analysis of the cell cycle, which is associated with survival in polyploidy, has proved useful in investigating polyploidization. Here, we used mRNA sequencing to investigate global expression in vitro (in cultured cells) and in vivo (in fin and liver tissues) in both the diploid and tetraploid Carassius auratus red var.. Differential expression (DE) of genes in diploid (7482, 36.0%) and tetraploid (3787, 18.2%) states suggested that in vitro and in vivo conditions dramatically change mRNA expression levels. However, of the 20,771 total shared expressed genes, 18,050 (87.0%), including 17,905 (86.2%) non-differentially expressed genes (DEGs) and 145 (0.7%) DEGs between diploids and tetraploids, showed the same expression trends in both cultured cells and liver tissues. Of the DEGs, four of seven genes in the cell cycle pathway had the same expression trends (upregulated in diploids and tetraploids) in both cultured cells and liver tissues. Quantitative PCR analysis confirmed the same expression trends in the nine DEGs associated with regulation of the cell cycle. This research on common characteristics between diploids and tetraploids provides insights into the potential molecular regulatory mechanisms of polyploidization. The steady changes that occur between diploids and tetraploids in vitro and in vivo show the potential value of studying polyploidy processes using cultured cell lines, especially with respect to cell cycle regulation.

Hox genes reveal variations in the genomic DNA of allotetraploid hybrids derived from Carassius auratus red var. (female) × Cyprinus carpio L. (male)

Background

Hox transcription factors are master regulators of animal development. Although highly conserved, they can contribute to the formation of novel biological characteristics when modified, such as during the generation of hybrid species, thus potentially serving as species-specific molecular markers. Here, we systematically studied the evolution of genomic sequences of Hox loci in an artificial allotetraploid lineage (4nAT, 4n = 200) derived from a red crucian carp (♀, RCC, 2n = 100) × common carp (♂, CC, 2n = 100) cross and its parents (RCC and CC).

Results

PCR amplification yielded 23 distinct Hox gene fragments from 160 clones in 4nAT, 22 fragments from 90 clones in RCC, and 19 fragments from 90 clones in CC. Sequence alignment of the HoxA3a and HoxC10a genes indicated both the inheritance and loss of paternal genomic DNA in 4nAT. The HoxA5a gene from 4nAT consisted of two subtypes from RCC and two subtypes from CC, indicating that homologous recombination occurred in the 4nAT hybrid genome. Moreover, 4nAT carried genomic pseudogenization in the HoxA10b and HoxC13a loci. Interestingly, a new type of HoxC9a gene was found in 4nAT as a hybrid sequence of CC and RCC by recombination in the intronic region.

Conclusion

The results revealed the influence of Hox genes during polyploidization in hybrid fish. The data provided insight into the evolution of vertebrate genomes and might be benefit for artificial breeding programs.

The Hybrid Genome of a New Goldfish-Like Fish Lineage Provides Insights Into the Origin of the Goldfish

Distant hybridization leads to obvious changes in genotypes and phenotypes, giving rise to species with novel capabilities. However, the fusion of distinct genomes also polymerizes the DNA or gene variations that occur during the course of evolution. Knowledge of the early stages of post-hybridization evolution is particularly important. Here, we investigated the full-length (FL) transcriptomes and the sequences resulting from the genome resequencing of the red crucian carp-like homodiploid fish (RCC-L) and goldfish-like homodiploid fish (GF-L) derived from the interspecific hybridization of koi carp (KOC) and blunt snout bream (BSB) to provide molecular evidence for the hybrid origin of the goldfish (GF). We compared the orthologous genes in the transcriptomes of RCC-L and GF-L with those of KOC and BSB. We also mapped the orthologous genes to the common carp (CC) and BSB genomes and classified them into eight gene patterns in three categories (chimaera, mutant, and biparental origin genes). The results showed that 48.20% and 46.50% of the genes were chimaera and that 3.70% and 8.30% of the genes were mutations of orthologous genes in RCC-L and GF-L, respectively. In RCC-L and GF-L, 63.70% and 68.20% of the genetic materials were from KOC, and 12.30% and 11.90% of the genetic materials were from BSB. The sequences from the genome resequencing of RCC-L and GF-L were mapped to the genome sequences of CC and BSB, revealing that the similarities of both RCC-L and GF-L to the CC genome (92.52%, 90.18%) were obviously higher than to the BSB genome (50.33%, 49.18%), supporting the suggestion that the genomes of both RCC-L and GF-L were mainly inherited from KOC but had some DNA fragments from BSB. Overall, our results provide molecular biological evidence for the hybrid origin of red crucian carp (RCC) and GF.