Whole genome sequencing of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) provide novel insights into their evolution and speciation

Comparative genomic analysis of immune-related genes and chemosensory receptors provides insights into the evolution and adaptation of four major domesticated Asian carps

BACKGROUND

Ctenopharyngodon idella (grass carp), Mylopharyngodon piceus (black carp), Hypophthalmichthys nobilis (bighead carp), and Hypophthalmichthys molitrix (silver carp), collectively known as the four major domesticated Asian carp, are freshwater fish species from the family Cyprinidae and are widely consumed in China. Current studies on these species primarily focus on immune system regulation and the growth and development of individual species. However, in-depth genomic investigations and comprehensive comparative analysis remained limited.

METHODS

The complete genomes of Ctenopharyngodon idella, Mylopharyngodon piceus and Hypophthalmichthys nobilis were assembled using a hybrid approach that integrated both next- and third-generation sequencing reads, followed by annotation using the MAKER2 pipeline. Based on the high-quality genomes of Ctenopharyngodon idella, Mylopharyngodon piceus Hypophthalmichthys nobilis, and Hypophthalmichthys molitrix, a comparative genomic analysis was conducted using bioinformatic tools to investigate gene family evolution in these four domesticated Asian carp species.

RESULTS

High-quality genomes of Ctenopharyngodon idella, Mylopharyngodon piceus, and Hypophthalmichthys nobilis were assembled, achieving over 90% completeness. Immune-related gene families, including MHC class I and NLRC3-like genes, have undergone rapid evolution, with Ctenopharyngodon idella exhibiting significant expansion of NLRC3-like genes. Massive tandem duplication events were identified in trace amine-associated receptors (TAARs), and rapid expansion was observed in TAAR16 and TAAR29. Additionally, a novel TAAR gene cluster was identified in all four Asian carp species. Comparative genomic analysis revealed the expansion of type 1 taste receptor genes, particularly in Ctenopharyngodon idella and Mylopharyngodon piceus.

CONCLUSION

This study has successfully constructed the high-quality genomes of Ctenopharyngodon idella, Mylopharyngodon piceus, and Hypophthalmichthys nobilis. The comparative genomic analysis revealed the evolution of immune-related genes and chemosensory receptors in the four major domesticated Asian carp species. These findings suggested the enhanced immunity and sensory perception in these species, providing valuable insights into their adaptation, survival and reproduction.

Multi-omics association pattern between gut microbiota and host metabolism of a filter-feeding fish in situ exposed to microplastics

Microplastics (MPs) are widespread in water environments and can affect gut microbiota and host metabolism of fish, but whether changes in host metabolism under MPs are mediated by gut microbiota remains unclear. Here, silver carp, a filter-feeding fish with important ecological functions, was in-situ exposure to environmentally relevant MPs. Multi-omics analysis and fecal microbiota transplantation were used to reveal the metabolic responses of carp along gut-liver-muscle axis. After three months of in situ exposure to MPs, community structure of gut microbiota of carp was reshaped, and five dominate phyla were significantly changed, including increased Cyanobacteria, Chloroflexi and Planctomycetota but decreased Firmicutes and Fusobacteriota. Weighted gene co-expression network analysis was further performed between these phyla and liver transcription spectrum, showing that the hub gene module contained up-regulated hppD, maiA and plg and activated ubiquinone and other terpenoid-quinone biosynthesis and phenylalanine metabolism. By fecal microbiota transplantation, the key gene module associated with core microbiota phyla of carp was verified in germ-free zebrafish. Interestingly, up-regulated hppD, maiA and plg and enriched phenylalanine metabolism were also observed in this module. Subsequently, metabolome performed in carp liver also shared activated phenylalanine metabolism, including increased trans-cinnamic acid and L-tyrosine. Furthermore, high-associated mapping showed that the differentially expressed metabolites (gamma-aminobutyric acid, ornithine and L-serine) related to amino acid metabolism in carp muscle were significantly accompanied with increased L-tyrosine in its liver. Overall, MPs exposure could change gut microbiome of silver carp and alter host metabolism especially amino acid metabolism along the gut-liver-muscle axis.

A chromosome-level genome assembly of the mud carp (Cirrhinus molitorella)

Algal blooms, which have become increasingly prevalent worldwide over the past decade, significantly impact on water quality and aquatic organisms. Filter-feeding fish are used to control phytoplankton and improve the ecological quality of water bodies. Mud carp (Cirrhinus molitorella) is a freshwater cyprinid species that predominantly consumes algae. Here, we generated a high-quality chromosome-level assembly of C. molitorella by integrating PacBio and Hi-C sequencing strategies. The genome assembly is 1.05 Gb, with a contig N50 of 24.13 Mb and a scaffold N50 of 39.38 Mb. The Benchmarking Universal Single-Copy Orthologs (BUSCO) (v4.0.5) benchmark of genome assembly reached 97.4% (95.8% single-copy). The consensus quality value (QV) and k-mer completeness of the C. molitorella assembly evaluated by Merqury software were 30.35 and 92.16%, respectively. The construction of the C. molitorella genome provides a valuable genetic resource that will facilitate the investigation of the digestion mechanism of filter-feeding fish.

Chromosome-Level Genome Assembly of Discogobio brachyphysallidos (Teleostei, Cyprinidae) and Population Genomics of the D. brachyphysallidos Complex: Impacts of Geological and Climate Changes on Species Evolution in Southwest China

The genus Discogobio is distributed in the eastern three rivers on the Yunnan-Guizhou Plateau and its adjacent regions, located to the southeast of the Qinghai-Tibet Plateau. Its origin and evolution are likely influenced by the uplift of the Qinghai-Tibet Plateau. However, the historical impact of geological events on the divergence and distribution of this fish group has not been fully elucidated. In this study, we successfully assembled a chromosome-level genome for Discogobio brachyphysallidos, which is approximately 1.21 Gb in length with a contig N50 of 8.63 Mb. The completeness of the genome assembly was assessed with a BUSCO score of 94.78%. A total of 30,597 protein-coding genes were predicted, with 93.92% functionally annotated. Phylogenetic analysis indicated that D. brachyphysallidos was closely related to Labeo rohita, and the divergence of the subfamily Labeoninae coincided with the significant uplift events of the Qinghai-Tibet Plateau. Additionally, we analyzed 75 samples of D. brachyphysallidos and D. yunnanensis from five populations, yielding 1.82 Tb of clean data and identifying 891,303,336 high-quality SNP sites. Population structure analyses indicated that the populations were clustered into five distinct groups, demonstrating significant genetic differentiation among them and the presence of cryptic species within this genus. Analyses of linkage disequilibrium decay and selective sweep indicated that the Pearl River population exhibited relatively higher genetic diversity compared with the populations from other drainages, and none of the populations showed evidence of expansion. Notably, the two population declines coincided with the early Pleistocene and Quaternary glaciation. It can be assumed that the geological movements of the Qinghai-Tibet Plateau and the Quaternary glaciation contributed to the decline in Discogobio populations and shaped their current size. The population genomics results showed that the present distribution pattern of Discogobio was the outcome of a series of geological events following the uplift of the Qinghai-Tibet Plateau. This study reconstructed the geological evolutionary history of the region from the perspective of species evolution. Furthermore, our study presents the first genome-wide analysis of the genetic divergence of Discogobio.

Chromosome-level genome assembly of the mud carp (Cirrhinus molitorella) using PacBio HiFi and Hi-C sequencing

The mud carp (Cirrhinus molitorella) is an important economic farmed fish, mainly distributed in South China and Southeast Asia due to its strong adaptability and high yield. Despite its economic importance, the paucity of genomic information has constrained detailed genetic research and breeding efforts. In this study, we utilized PacBio HiFi long-read sequencing and Hi-C technologies to generate a meticulously assembled chromosome-level genome of the mud carp. This assembly spans 1,033.41 Mb, with an impressive 99.82% distributed across 25 chromosomes. The contig N50 and scaffold N50 are 33.29 Mb and 39.86 Mb, respectively. The completeness of the mud carp genome assembly is highlighted by a BUSCO score of 98.05%. We predict 25,865 protein-coding genes, with a BUSCO score of 96.54%, and functional annotations for 91.83% of these genes. Approximately 52.21% of the genome consists of repeat elements. This high-fidelity genome assembly is a vital resource for advancing molecular breeding, comparative genomics, and evolutionary studies of the mud carp and related species.

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Genetic improvement and genomic resources of important cyprinid species: status and future perspectives for sustainable production

Cyprinid species are the most cultured aquatic species around the world in terms of quantity and total value. They account for 25% of global aquaculture production and significantly contribute to fulfilling the demand for fish food. The aquaculture of these species is facing severe concerns in terms of seed quality, rising feed costs, disease outbreaks, introgression of exotic species, environmental impacts, and anthropogenic activities. Numerous researchers have explored biological issues and potential methods to enhance cyprinid aquaculture. Selective breeding is extensively employed in cyprinid species to enhance specific traits like growth and disease resistance. In this context, we have discussed the efforts made to improve important cyprinid aquaculture practices through genetic and genomic approaches. The recent advances in DNA sequencing technologies and genomic tools have revolutionized the understanding of biological research. The generation of a complete genome and other genomic resources in cyprinid species has significantly strengthened molecular-level investigations into disease resistance, growth, reproduction, and adaptation to changing environments. We conducted a comprehensive review of genomic research in important cyprinid species, encompassing genome, transcriptome, proteome, metagenome, epigenome, etc. This review reveals that considerable data has been generated for cyprinid species. However, the seamless integration of this valuable data into genetic selection programs has yet to be achieved. In the upcoming years, genomic techniques, gene transfer, genome editing tools are expected to bring a paradigm shift in sustainable cyprinid aquaculture production. The comprehensive information presented here will offer insights for the cyprinid aquaculture research community.

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae.

RESULTS

Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity.

CONCLUSIONS

Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

Transcriptomic atlas for hypoxia and following re-oxygenation in Ancherythroculter nigrocauda heart and brain tissues: insights into gene expression, alternative splicing, and signaling pathways

Hypoxia is a mounting problem that affects the world's freshwaters, with severe consequence for many species, including death and large economical loss. The hypoxia problem has increased recently due to the combined effects of water eutrophication and global warming. In this study, we investigated the transcriptome atlas for the bony fish Ancherythroculter nigrocauda under hypoxia for 1.5, 3, and 4.5 h and its recovery to normal oxygen levels in heart and brain tissues. We sequenced 21 samples for brain and heart tissues (a total of 42 samples) plus three control samples and obtained an average of 32.40 million raw reads per sample, and 95.24% mapping rate of the filtered clean reads. This robust transcriptome dataset facilitated the discovery of 52,428 new transcripts and 6,609 novel genes. In the heart tissue, the KEGG enrichment analysis showed that genes linked to the Vascular smooth muscle contraction and MAPK and VEGF signaling pathways were notably altered under hypoxia. Re-oxygenation introduced changes in genes associated with abiotic stimulus response and stress regulation. In the heart tissue, weighted gene co-expression network analysis pinpointed a module enriched in insulin receptor pathways that was correlated with hypoxia. Conversely, in the brain tissue, the response to hypoxia was characterized by alterations in the PPAR signaling pathway, and re-oxygenation influenced the mTOR and FoxO signaling pathways. Alternative splicing analysis identified an average of 27,226 and 28,290 events in the heart and brain tissues, respectively, with differential events between control and hypoxia-stressed groups. This study offers a holistic view of transcriptomic adaptations in A. nigrocauda heart and brain tissues under oxygen stress and emphasizes the role of gene expression and alternative splicing in the response mechanisms.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Identification and Characterization of Immune-Associated MicroRNAs in Silver Carp (Hypophthalmichthys molitrix) Responding to Aeromonas veronii and LPS Stimulation

The ubiquitous Gram-negative bacterial pathogen Aeromonas veronii (A. veronii) can easily cause inflammatory reactions in aquatic organisms, resulting in high mortality and huge economic losses. MicroRNAs (miRNAs) participate in immune regulation and have certain conserved properties. MiRNAs are involved in the immune responses of a variety of teleost fish infected with bacteria, whereas there is no related report in silver carp (Hypophthalmichthys molitrix). Therefore, we identified the expression profiles of miRNA in silver carp stimulated by A. veronii and LPS. Among them, the quantity of differentially expressed miRNAs (DEmiRNAs) obtained in the silver carp challenge group was 73 (A. veronii) and 90 (LPS). The GO enrichment and analysis of KEGG pathways have shown that the predicted target genes are mainly associated with lipid metabolism and the immune response in silver carp. This indicates the possibility that miRNAs play a role in regulating immune-related pathways. In addition, a total of eight DEmiRNAs validated the accuracy of the sequencing result via quantitative real-time PCR (qRT-PCR). Finally, we selected the silver carp head kidney macrophage cells (HKCs) as model cells and proved that miR-30b-5p can regulate the inflammatory response in silver carp HKCs. This study lays the foundation for exploring miRNA regulation in silver carp during pathogenic bacterial infection. In addition, it provides a reference for the future development of non-coding RNA antibacterial drugs.

Characteristics of conserved microRNAome and their evolutionary adaptation to regulation of immune defense functions in the spleen of silver carp and bighead carp

Immune defense functions of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) have shown obvious evolutionary divergence. MiRNAs participate in the fine regulation of immune function. However, the evolutionary adaptation of miRNAs in the regulation of immune defense function is still poorly understood in silver carp and bighead carp. Here, small RNA libraries were constructed from the spleen tissue of one-year-old and three-year-old healthy silver carp and bighead carp, 424 and 422 known conserved miRNAs were respectively identified from the spleen of silver carp and bighead carp by bioinformatic analysis, which 398 were shared between the two species. These conserved miRNAs showed highly similar expression patterns between silver carp and bighead carp, but the abundance in spleen varied greatly in different species. Family analysis showed that miRNA families including mir-8, mir-7, mir-23, mir-338, mir-30, mir-27, mir-221, mir-19, mir-181, mir-17, mir-15, mir-148, mir-130, mir-10 and let-7 were the main miRNAs in the spleen of silver carp and bighead carp. 27 and 51 significant differentially expressed (SDE) miRNAs were identified from silver carp and bighead carp, respectively. Evolution analysis for the predicted target genes of SDE-miRNAs showed that ten biological processes such as blood coagulation, cell adhesion mediated by integrin and adaptive immune response were positively selected. In addition, immune genes including TLR3, NFATC3, MALT1, B2M, GILT and MHCII were positively selected only in silver carp, and they were specifically targeted by the SDE-miRNAs including miR-9-5p, miR-196a-5p, miR-375, miR-122, miR-722, miR-132-3p, miR-727-5p, miR-724, miR-19d-5p and miR-138-5p, respectively. PLA2G4 in Fc epsilon RI signaling pathway was positively selected only in bighead carp and was specifically targeted by the SDE-miRNAs including miR-222b, miR-22b-5p, miR-15c, miR-146a, miR-125c-3p, miR-221-5p, miR-2188-5p, miR-142a-3p, miR-212, miR-138-5p and miR-15b-5p. In particular, SDE-miRNAs such as miR-144-3p, miR-2188-3p, miR-731, miR-363-3p and miR-218b could simultaneously target multiple evolutionarily differentiated immune-related genes. These results indicated that in the spleen of silver carp and bighead carp, conserved miRNAs have obvious evolutionary adaptations in the regulation of immune defense function. The results of this study can provide valuable resources for further revealing themechanism of miRNA in the formation of resistance traits evolution between silver carp and bighead carp.

Genetic Diversity and Signatures of Selection in the Roughskin Sculpin (Trachidermus fasciatus) Revealed by Whole Genome Sequencing

The roughskin sculpin (Trachidermus fasciatus) is an endangered fish species in China. In recent years, artificial breeding technology has made significant progress, and the population of roughskin sculpin has recovered in the natural environment through enhancement programs and the release of juveniles. However, the effects of released roughskin sculpin on the genetic structure and diversity of wild populations remain unclear. Studies on genetic diversity analysis based on different types and numbers of molecular markers have yielded inconsistent results. In this study, we obtained 2,610,157 high-quality SNPs and 494,698 InDels through whole-genome resequencing of two farmed populations and one wild population. Both farmed populations showed consistent levels of genomic polymorphism and a slight increase in linkage compared with wild populations. The population structure of the two farmed populations was distinct from that of the wild population, but the degree of genetic differentiation was low (overall average Fst = 0.015). Selective sweep analysis showed that 523,529 genes were selected in the two farmed populations, and KEGG enrichment analysis showed that the selected genes were related to amino acid metabolism, which might be caused by artificial feeding. The findings of this study provide valuable additions to the existing genomic resources to help conserve roughskin sculpin populations.

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.

Chromosome-Level Assembly of Male Opsariichthys bidens Genome Provides Insights into the Regulation of the GnRH Signaling Pathway and Genome Evolution

The hook snout carp Opsariichthys bidens is an important farmed fish in East Asia that shows sexual dimorphism in growth, with males growing faster and larger than females. To understand these complex traits and improve molecular breeding, chromosome-level genome assembly of male O. bidens was performed using Illumina, Nanopore, and Hi-C sequencing. The 992.9 Mb genome sequences with a contig N50 of 5.2 Mb were anchored to 38 chromosomes corresponding to male karyotypes. Of 30,922 functionally annotated genes, 97.5% of BUSCO genes were completely detected. Genome evolution analysis showed that the expanded and contracted gene families in the male O. bidens genome were enriched in 76 KEGG pathways, and 78 expanded genes were involved in the GnRH signaling pathway that regulates the synthesis and secretion of luteinizing hormone and glycoprotein hormones, further acting on male growth by inducing growth hormone. Compared to the released female O. bidens genome, the number of annotated genes in males was much higher (23,992). The male chromosome LG06 exhibited over 97% identity with the female GH14/GH38. Male-specific genes were identified for LG06, where structural variation, including deletions and insertions, occurred at a lower rate, suggesting a centric fusion of acrocentric chromosomes GH14 and GH38. The genome-synteny analysis uncovered significant inter-chromosome conservation between male O. bidens and grass carp, the former originating from ancestral chromosome breakage to increase the chromosome number. Our results provide a valuable genetic resource for studying the regulation of sexual dimorphism, sex-determining mechanisms, and molecular-guided breeding of O. bidens.

A chromosome-level genome assembly of the redfin culter (Chanodichthys erythropterus)

Chanodichthys erythropterus is a fierce carnivorous fish widely found in East Asian waters. It is not only a popular food fish in China, it is also a representative victim of overfishing. Genetic breeding programs launched to meet market demands urgently require high-quality genomes to facilitate genomic selection and genetic research. In this study, we constructed a chromosome-level reference genome of C. erythropterus by taking advantage of long-read single-molecule sequencing and de novo assembly by Oxford Nanopore Technology (ONT) and Hi-C. The 1.085 Gb C. erythropterus genome was assembled from 132 Gb of Nanopore sequence. The assembled genome represents 98.5% completeness (BUSCO) with a contig N50 length of 23.29 Mb. The contigs were clustered and ordered onto 24 chromosomes covering roughly 99.49% of the genome assembly with Hi-C data. Additionally, 33,041 (98.0%) genes were functionally annotated from a total of 33,706 predicted protein-coding sequences by combining transcriptome data from seven tissues. This high-quality assembled genome will be a precious resource for future molecular breeding and functional genomics research of C. erythropterus.

Chromosome-level genome assembly of grass carp (Ctenopharyngodon idella) provides insights into its genome evolution

Background

The grass carp has great economic value and occupies an important evolutionary position. Genomic information regarding this species could help better understand its rapid growth rate as well as its unique body plan and environmental adaptation.

Results

We assembled the chromosome-level grass carp genome using the PacBio sequencing and chromosome structure capture technique. The final genome assembly has a total length of 893.2 Mb with a contig N50 of 19.3 Mb and a scaffold N50 of 35.7 Mb. About 99.85% of the assembled contigs were anchored into 24 chromosomes. Based on the prediction, this genome contained 30,342 protein-coding genes and 43.26% repetitive sequences. Furthermore, we determined that the large genome size can be attributed to the DNA-mediated transposable elements which accounted for 58.9% of the repetitive sequences in grass carp. We identified that the grass carp has only 24 pairs of chromosomes due to the fusion of two ancestral chromosomes. Enrichment analyses of significantly expanded and positively selected genes reflected evolutionary adaptation of grass carp to the feeding habits. We also detected the loss of conserved non-coding regulatory elements associated with the development of the immune system, nervous system, and digestive system, which may be critical for grass carp herbivorous traits.

Conclusions

The high-quality reference genome reported here provides a valuable resource for the genetic improvement and molecular-guided breeding of the grass carp.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08503-x.

Updated Genome Assembly of Bighead Carp (Hypophthalmichthys nobilis) and Its Differences Between Male and Female on Genomic, Transcriptomic, and Methylation Level

Cyprinidae is one of the largest family in freshwater fishes, and it is most intensively cultured fish taxon of the world. However, studies about sex determination in this large family is still rear, and one of the reasons is lack of high quality and complete genome. Here, we used nanopore to sequence the genome of a male bighead carp, obtaining contig N50 = 24.25 Mb, which is one of the best assemblies in Cyprinidae. Five males and five females were re-sequenced, and a male-specific region on LG19 was confirmed. We find this region holds many male-specific markers in other Cyprinidae fishes, such as grass carp and silver carp. Transcriptome analyses of hypothalamus and pituitary tissues showed that several sex-specific differentially expressed genes were associated with steroid biosynthesis. The UCH64E gene, located in the male-specific region on LG19, showed higher expression levels in male than female tissues of bighead carp. The methyl-RAD of hypothalamus tissues between males and females indicated that the sexual methylation differences are significant in bighead carp. We also compared the methylation sites recognized using methyl-RAD and nanopore raw reads and found that approximately 73% of the methylation sites identified using methyl-RAD were within nanopore CpG sites.

Chromosome-level assembly of the Hypophthalmichthys molitrix (Cypriniformes: Cyprinidae) genome provides insights into its ecological adaptation

Hypophthalmichthys molitrix (silver carp) is phytoplanktivorous and is an economically and ecologically important fish species. As a well-known invasive species, a number of factors associated with the ecological adaptations of this species are largely unknown. Here, we present a chromosomal-level assembly of the species based on the PacBio Sequel II platform and Hi-C scaffolding technology. Based on the high-quality genome sequences and previous genome sequencing projects, a number of genes that were probably subject to positive selection reside in the genome of H. molitrix, and the last common ancestors of H. molitrix and H. nobilis were identified. Some of these genes may partially explain the mechanisms of H. molitrix for surviving damage due to toxic algae. Demographic history estimation suggests that the effective population size (EPS) of the species may have constantly increased along with the uplift of the Qinghai-Tibet Plateau, started to decline when quaternary glaciation started, and further declined during the Younger Dryas Period. Moreover, the introgression from H. nobilis to H. molitrix in North America was corroborated based on the whole-genome sequencing data, and the proportion of introgressed regions was estimated to be approximately 5.8%. Based on the high-quality assembly, the possible mechanisms by which H. molitrix adapts to its endemic and invaded locations were profiled.

Genome of the butterfly hillstream loach provides insights into adaptations to torrential mountain stream life

Butterfly hillstream loach (Beaufortia kweichowensis), a benthic fish in the torrential mountain streams, possesses a totally flat ventrum, flattened craniofacial and body skeletons, and enlarged paired fins covered by substantially small keratinous structures. However, little is known about the genetic basis of these specialized morphological adaptations. Here we present a 448.52-Mb genome assembly with contig N50 length of 5.53 Mb by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin map. Demographic history reconstruction of the butterfly hillstream loach reveals that the population dynamics is correlated with the different stages of uplifting of the Tibetan Plateau. Comparative genomic analysis finds evidence of six keratin genes in butterfly hillstream loach evolving under positive selection. Within these genes, two keratin genes exhibit species-specific and divergent amino acid changes, suggesting a role in the formation of the unculi. Additionally, a series of positively selected genes, rapid evolving genes, specific variant genes and expanded gene families are found, including genes related to Hedgehog, Notch and BMP pathways, which may be involved in craniofacial development. These findings may have important implications for understanding the genetic basis of phenotypic adaptation to torrential mountain stream life.