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

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.

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.

Measurement(s)	whole genome sequencing
Technology Type(s)	Oxford Nanopore Sequencing
Sample Characteristic - Organism	Chanodichthys 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.