Genome sequence and genetic diversity of the common carp, Cyprinus carpio

The hybrid history of zebrafish

Since the description of zebrafish (Danio rerio) in 1822, the identity of its closest living relative has been unclear. To address this problem, we sequenced the exomes of 10 species in genus Danio, using the closely related Devario aequipinnatus as outgroup, to infer relationships across the 25 chromosomes of the zebrafish genome. The majority of relationships within Danio were remarkably consistent across all chromosomes. Relationships of chromosome segments, however, depended systematically upon their genomic location within zebrafish chromosomes. Regions near chromosome centers identified Danio kyathit and/or Danio aesculapii as the closest relative of zebrafish, while segments near chromosome ends supported only D. aesculapii as the zebrafish sister species. Genome-wide comparisons of derived character states revealed that danio relationships are inconsistent with a simple bifurcating species history but support an ancient hybrid origin of the D. rerio lineage by homoploid hybrid speciation. We also found evidence of more recent gene flow limited to the high recombination ends of chromosomes and several megabases of chromosome 20 with a history distinct from the rest of the genome. Additional insights gained from incorporating genome structure into a phylogenomic study demonstrate the utility of such an approach for future studies in other taxa. The multiple genomic histories of species in the genus Danio have important implications for comparative studies in these morphologically varied and beautiful species and for our understanding of the hybrid evolutionary history of zebrafish.

Genome-wide identification, evolution, and expression analysis of the bone morphogenetic protein gene family in Myxocyprinus asiaticus

Bone morphogenetic proteins (BMPs) are important growth factors belonging to the TGF-β superfamily. These factors not only play a vital role in skeleton formation in young fish but also regulate the morphological development of M. asiaticus, with Group II genes regulating morphology mainly during the juvenile stage. This study investigated how BMP genes regulate Myxocyprinus asiaticus development and function and explored the role of the BMP family in fish morphological development. In this study, 43 BMPs were identified and classified into five groups: BMP1/3/11/15 (Group I), BMP12/13/14 (Group II), BMP2/4/16 (Group III), BMP9/10 (Group IV), and BMP5/6/7/8 (Group V). Analyses of the gene structures and conserved motifs revealed the conservation of the BMP gene family in M. asiaticus. In M. asiaticus, gene fragmentation, duplication, and 4R whole-genome duplication events contributed to BMP gene family expansion. Furthermore, expression pattern analysis and qRT-PCR revealed that changes in M. asiaticus BMP gene expression during different developmental stages were due to body size alterations, highlighting the major impact of the BMP gene on body size variation in this species. Our study provides fundamental data for investigating the morphological development of M. asiaticus and lays the framework for understanding the genetic mechanisms of body size variation in scleractinian fishes, with potential applications in the artificial breeding and conservation of M. asiaticus.

Polyploidization-driven transcriptomic dynamics in Medicago sativa neotetraploids: mRNA, smRNA and allele-specific gene expression

Whole genome duplication (WGD) is a powerful evolutionary mechanism in plants. Autopolyploids have been comparatively less studied than allopolyploids, with sexual autopolyploidization receiving even less attention. In this work, we studied the transcriptomes of neotetraploids (2n = 4x = 32) obtained by crossing two diploid (2n = 2x = 16) plants of Medicago sativa that produce a significant percentage of either 2n eggs or pollen. Diploid progeny from the same cross allowed us to separate the transcriptional outcomes of hybridization from those of WGD. This material can help to elucidate events at the base of the domestication of cultivated 4x alfalfa, the world's most important leguminous forage. Three 2x and three 4x progeny plants and 2x parental plants were used for this study. The RNA-seq data revealed that WGD did not dramatically affect the transcription of leaf protein-coding genes. The two parental genotypes did not contribute equally to the progeny transcriptomes, and genome-wide expression level dominance of the male parent was observed. A large majority of the genes whose expression level changed due to WGD presented increased expression, indicating that the 4x state requires the upregulation of approximately 2.66% of the protein-coding genes. Overall, we estimated that 3.63% of the protein-coding genes were transcriptionally affected by WGD and may contribute to the phenotypic novelty of the neotetraploid plants. Pathway analysis suggested that WGD could affect secondary metabolite biosynthesis, which in turn may influence forage quality. We found four times as many transcription factor genes among the polyploidization-affected genes than among those affected only by hybridization. Several of these belong to classes involved in stress response. Small RNA-seq revealed that very few miRNAs were significantly associated with WGD, but they target several hundred genes, and their role in the WGD response may be relevant. Integrated network analysis led to the identification of putative miRNA: mRNA interactions potentially involved in transcriptome reprogramming. Allele-specific expression analysis indicated that parent-of-origin bias was not a significant outcome of WGD, but we found that parentally biased RNA editing may be a significant source of variation in neopolyploids.

Variation and Interaction of Distinct Subgenomes Contribute to Growth Diversity in Intergeneric Hybrid Fish

Intergeneric hybridization greatly reshapes regulatory interactions among allelic and non-allelic genes. However, their effects on growth diversity remain poorly understood in animals. In this study, we conducted whole-genome sequencing and RNA sequencing analyses in diverse hybrid varieties resulting from the intergeneric hybridization of goldfish (Carassius auratus red var.) and common carp (Cyprinus carpio). These hybrid individuals were characterized by distinct mitochondrial genomes and copy number variations. Through a weighted gene correlation network analysis, we identified 3693 genes as candidate growth-regulating genes. Among them, the expression of 3672 genes in subgenome R (originating from goldfish) displayed negative correlations with body weight, whereas 20 genes in subgenome C (originating from common carp) exhibited positive correlations. Notably, we observed intriguing expression patterns of solute carrier family 2 member 12 (slc2a12) in subgenome C, showing opposite correlations with body weight that changed with water temperatures, suggesting differential interactions between feeding activity and weight gain in response to seasonal changes for hybrid animals. In 40.30% of alleles, we observed dominant trans-regulatory effects in the regulatory interactions between distinct alleles from subgenomes R and C. Integrating analyses of allele-specific expression and DNA methylation data revealed that DNA methylation on both subgenomes shaped the relative contribution of allelic expression to the growth rate. These findings provide novel insights into the interactions of distinct subgenomes that underlie heterosis in growth traits and contribute to a better understanding of multiple allelic traits in animals.

A chromosome-level, haplotype-resolved genome assembly and annotation for the Eurasian minnow (Leuciscidae: Phoxinus phoxinus) provide evidence of haplotype diversity

BACKGROUND

In this study, we present an in-depth analysis of the Eurasian minnow (Phoxinus phoxinus) genome, highlighting its genetic diversity, structural variations, and evolutionary adaptations. We generated an annotated haplotype-phased, chromosome-level genome assembly (2n = 50) by integrating high-fidelity (HiFi) long reads and chromosome conformation capture data (Hi-C).

RESULTS

We achieved a haploid size of 940 megabase pairs (Mbp) for haplome 1 and 929 Mbp for haplome 2 with high scaffold N50 values of 36.4 Mb and 36.6 Mb and BUSCO scores of 96.9% and 97.2%, respectively, indicating a highly complete genome assembly. We detected notable heterozygosity (1.43%) and a high repeat content (approximately 54%), primarily consisting of DNA transposons, which contribute to genome rearrangements and variations. We found substantial structural variations within the genome, including insertions, deletions, inversions, and translocations. These variations affect genes enriched in functions such as dephosphorylation, developmental pigmentation, phagocytosis, immunity, and stress response. In the annotation of protein-coding genes, 30,980 messenger RNAs and 23,497 protein-coding genes were identified with a high completeness score, which further underpins the high contiguity of our genome assemblies. We performed a gene family evolution analysis by comparing our proteome to 10 other teleost species, which identified immune system gene families that prioritize histone-based disease prevention over NB-LRR-related-based immune responses. Additionally, demographic analysis indicates historical fluctuations in the effective population size of P. phoxinus, likely correlating with past climatic changes.

CONCLUSIONS

This annotated, phased reference genome provides a crucial resource for resolving the taxonomic complexity within the genus Phoxinus and highlights the importance of haplotype-phased assemblies in understanding haplotype diversity in species characterized by high heterozygosity.

Viperin and Its Effect on SVCV Replication in Common Carp, Cyprinus carpio

Viperin is an interferon-stimulated gene (ISG) that plays an important role in the congenital antiviral immunity of vertebrates. In this study, the common carp viperin (cc-viperin) gene is characterized, and we determine whether it has the ability to inhibit spring viremia of carp virus (SVCV) replication in EPC cells. The results showed that the full-length cDNA of the cc-viperin gene was 1044 bp and it encoded 348 amino acids. The cc-viperin sequence contained a leucine zipper in the N-terminal, a CxxxCxxC motif in the SAM domain, and a conservative C-terminus. The cc-viperin gene's nucleotide and amino acid sequence alignment revealed that cc-viperin displayed relatively high sequence identity compared with other species. Phylogenetic analysis displayed the close relation of cc-viperin with Carassius auratus and Mylopharyngodon piceus. Subcellular localization analysis indicated that the cc-viperin protein was located in the cytoplasm. The gene expression results showed that cc-viperin was expressed in all of the tissues tested. Its expression level significantly increased in EPC cells after 24 h to 72 h compared to the control during SVCV infection. Moreover, cc-viperin significantly inhibited SVCV replication when it was overexpressed, whereas it increased SVCV replication when it had reduced expression in EPC cells, respectively. To summarize, the results obtained in this work show that cc-viperin shares similar sequence characteristics with other vertebrates, and it could inhibit SVCV replication in EPC cells, displaying an antiviral effect in common carp.

Molecular responses of paddy field carp (Cyprinus carpio) in the agricultural heritage to major environmental factors in paddy fields

As a core element of the Globally Important Agricultural Heritage System (GIAHS), the Qingtian paddy field carp (Cyprinus carpio, PF-carp) has been domesticated for over 1200 years in paddy field environments. This species has successfully adapted to shallow-water conditions in paddy fields. To reveal the adaptation mechanism, we conducted transcriptome sequencing on the hepatopancreas of PF-carp under two temperature conditions (28 °C and 38 °C) and concurrently analysed RNA-seq data from hypoxic conditions in the same tissue. By analysing high-temperature transcriptome data, 3154 differentially expressed genes (DEGs) were identified. KEGG analysis indicated that DEGs involved various pathways, including protein processing in endoplasmic reticulum, circadian rhythm, and HIF-1 signaling pathway. Notably, protein processing in endoplasmic reticulum was significantly enriched with key genes such as HSP70, HSP90, HSP40, CNX, CRT, and Bip. Through concurrent analysis of RNA-seq data from hypoxic conditions, we found that PF-carp regulate their metabolism through multiple pathways and produce almost opposite metabolic regulation to adapt to high temperature and hypoxic environments. The opposite activation state observed in the HIF-1 signaling pathway is particularly intriguing. In conclusion, PF-carp appear to rely on protein processing in endoplasmic reticulum to maintain cell homeostasis at high temperatures. The HIF-1 signaling pathway may emerged as a key player in adapting PF-carps to paddy fields. This study provides valuable insights into the adaptive mechanisms of domesticated fish in paddy fields.

Gene expression signatures between Limia perugiae (Poeciliidae) populations from freshwater and hypersaline habitats, with comparisons to other teleosts

Salinity gradients act as strong environmental barriers that limit the distribution of aquatic organisms. Changes in gene expression associated with transitions between freshwater and saltwater environments can provide insights into organismal responses to variation in salinity. We used RNA-sequencing (RNA-seq) to investigate genome-wide variation in gene expression between a hypersaline population and a freshwater population of the livebearing fish species Limia perugiae (Poeciliidae). Our analyses of gill gene expression revealed potential molecular mechanisms underlying salinity tolerance in this species, including the enrichment of genes involved in ion transport, maintenance of chemical homeostasis, and cell signaling in the hypersaline population. We also found differences in gene expression patterns associated with cell-cycle and protein-folding processes between the hypersaline and freshwater L. perugiae. Bidirectional freshwater-saltwater transitions have occurred repeatedly during the diversification of fishes, allowing for broad-scale examination of repeatable patterns in evolution. Therefore, we compared transcriptomic variation in L. perugiae with other teleosts that have made freshwater-saltwater transitions to test for convergence in gene expression. Among the four distantly related population pairs from high- and low-salinity environments that we included in our analysis, we found only ten shared differentially expressed genes, indicating little evidence for convergence. However, we found that differentially expressed genes shared among three or more lineages were functionally enriched for ion transport and immune functioning. Overall, our results-in conjunction with other recent studies-suggest that different genes are involved in salinity transitions across disparate lineages of teleost fishes.

The whole chromosome-level genome provides resources and insights into the endangered fish Percocypris pingi evolution and conservation

BACKGROUND

Percocypris pingi (Tchang) was classified as Endangered on the Red List of China's Vertebrates in 2015 and is widely distributed in the Upper Yangtze River. Although breeding and release into wild habitats have been performed for this commercially important fish in recent years, low genetic diversity has been found in wild populations. Genomic resources are strongly recommended before formulating and carrying out conservation strategies for P. pingi. Thus, there is an urgent need to conserve germplasm resources and improve the population diversity of P. pingi. To date, the whole genome of P. pingi has not been reported.

RESULTS

In our study, we constructed the first chromosome-level genome of P. pingi by high-throughput chromosome conformation capture (Hi-C) technology and PacBio long-read sequencing. The assembled genome was 1.7 Gb in size, with an N50 of 17,692 bp and a GC content from circular consensus sequencing of 37.67%. The Hi-C results again demonstrated that P. pingi was tetraploid (n = 98), with the genome consisting of 24-type and 25-type chromosomes. Chr.19 of the 24-type chromosomes in P. pingi resulted from the fusion of chr.19 and chr.22 in zebrafish. The divergence times between 24-type and 25-type chromosomes was around 6.1 million years ago. A total of 25,198 and 25,291 protein-coding genes were obtained from the 24-type and 25-type chromosomes, respectively. The ploidy of P. pingi is an allotetraploid. A total of 8,741 genes of P. pingi were clustered into 4,378 gene families that were shared with 14 other species, and the P. pingi genome had 68 unique gene families. Phylogenetic analyses indicated that P. pingi was most closely related to Schizothorax oconnori, and the genes were clustered on one branch. We identified 166 significantly expanded gene families and 173 significantly contracted gene families in P. pingi. The most enriched positive protein-coding genes, such as Bmp-4, Etfdh, homeobox protein HB9, and ATG3, were screened.

CONCLUSION

Our study provides a high-quality chromosome-anchored reference genome for P. pingi and provides sufficient information on the chromosomes, which will lead to valuable resources for genetic, genomic, and biological studies of P. pingi and for improving the genetic diversity, population size, and scientific conservation of endangered fish and other key cyprinid species in aquaculture.

Gap-free telomere-to-telomere haplotype assembly of the tomato hind (Cephalopholis sonnerati)

The tomato hind (Cephalopholis sonnerati) is an emerging economically important grouper in recent years. With the increasing maturity of sequencing technologies and assembly methodologies, a higher quality reference genome has become both accessible and necessary. In this study, we present two telomere-to-telomere (T2T) gap-free haplotype assemblies of the tomato hind with lengths of 1039.53 Mb (YSFRI_Csonn_HA_1.0, N50 43.83 Mb) and 1039.91 Mb (YSFRI_Csonn_HB_1.0, N50 44.09 Mb). Reads from next-generation sequencing, ONT ultra-long sequencing, and PacBio HiFi sequencing exhibited mapping rates exceeding 99.8% when aligned to these two assemblies. Evaluation using Merqury indicated high accuracy for both assemblies, with average quality values of 51.80 and 51.83, respectively. Percentages of 97.9% and 97.8% of complete BUSCOs were achieved, and a total of 23,270 and 23,184 protein-code genes were inferred in each assembly. Moreover, telomere identification, centromere prediction, and repetitive sequence annotation were also successfully performed. These two assemblies provide robust foundation for the genetic analysis and development of molecular genetic breeding technologies in C. sonnerati.

Characterization of Feeding Behaviors, Appetite Regulation and Growth Performance of All-Female (cyp17a1+/-;XX Genotype) Common Carp (Cyprinus carpio)

Genome editing has the potential to improve growth and traits of aquatic animals. Assessment of the feeding habits of the genetically modified farmed fish is necessary, as this is closely related to the assessment of their growth performance, which is one of the most important economic traits. Previously, we developed a novel strategy to produce all-female (AF) common carp (cyp17a1+/-;XX genotype) with genome editing, which exhibited a growth advantage compared to the control carp (including control male and female carp). However, the feeding behavior related to the growth performance of wild-type control and AF common carp remains elusive. The results of feeding and swimming behaviors showed that AF common carp exhibited a faster feeding activities and more active swimming activities, which probably enhanced its growth performance. Brain gene expression analysis revealed AF common carp had a significant upregulation of the orexigenic factors gene expression levels in the fed state, which would further promote the growth of AF carp. Here, AF carp exhibited higher growth performance with higher growth hormone (gh) gene expression. This study provided insight into the growth performance, feeding behaviors and appetite regulation of the genetically modified AF carp and the assessment of feeding behaviors in other genetically modified farmed fish.

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.

Comparative Analysis of Growth Performance, Morphological Development, and Physiological Condition in Three Romanian Cyprinus carpio Varieties and Koi: Implications for Aquaculture

This study investigates the influence of internal factors on growth dynamics in four Cyprinus carpio varieties, three Romanian strains (Frăsinet, Ineu, and Podu Iloaiei) and the Koi variety. Fish were measured for total length, maximum height, and weight at four ontogenic stages, namely 7 days post-hatch, 3 months (0+), 18 months old (1+), and 36 months (2+). Weight Gain (WG), Specific Growth Rate (SGR), Relative Growth Rate (RGR), Fulton's condition factor, and the profile index were calculated and analyzed. Results revealed significant intervariety differences in growth performance and physiological condition across life stages. At the 2+ stage, Podu Iloaiei exhibited the highest WG (849.73 ± 4.09 g), while Koi showed the lowest (403.99 ± 14.21 g). Koi demonstrated a unique growth pattern, with the highest SGR (0.18 ± 0.00% day-1) and RGR (0.98 ± 0.05 g day-1) at the 2+ stage. Fulton's condition factor varied markedly, with Frăsinet showing the highest value at 7 days post-hatch (149.57 ± 17.485) and Koi the lowest at the 1+ stage (0.63 ± 0.011). The profile index decreased with age in all varieties, with Podu Iloaiei showing the most dramatic change from 4.22 ± 0.149 at 7 days to 2.18 ± 0.004 at 2+. These findings highlight the complex interplay between genetic and developmental factors in carp varieties, offering new insights for tailored breeding programs and aquaculture practices.

Recent and Recurrent Autopolyploidization Fueled Diversification of Snow Carp on the Tibetan Plateau

Whole-genome duplication (WGD), or polyploidization, is a major contributor to biodiversity. However, the establishment and survival of WGDs are often considered to be stochastic, since elucidating the processes of WGD establishment remains challenging. In the current study, we explored the processes leading to polyploidy establishment in snow carp (Cyprinidae: Schizothoracinae), a predominant component of the ichthyofauna of the Tibetan Plateau and its surrounding areas. Using large-scale genomic data from isoform sequencing, we analyzed ohnolog genealogies and divergence in hundreds to thousands of gene families across major snow carp lineages. Our findings demonstrated that independent autopolyploidization subsequent to speciation was prevalent, while autopolyploidization followed by speciation also occurred in the diversification of snow carp. This was further supported by matrilineal divergence and drainage evolution evidence. Contrary to the long-standing hypothesis that ancient polyploidization preceded the diversification of snow carp, we determined that polyploidy in extant snow carp was established by recurrent autopolyploidization events during the Pleistocene. These findings indicate that the diversification of extant snow carp resembles a coordinated duet: first, the uplift of the Tibetan Plateau orchestrated the biogeography and diversification of their diploid progenitors; then, the extensive Pliocene-Pleistocene climate changes acted as relay runners, further fueling diversification through recurrent autopolyploidization. Overall, this study not only reveals a hitherto unrecognized recent WGD lineage in vertebrates but also advances current understanding of WGD processes, emphasizing that WGD establishment is a nonstochastic event, emerging from numerous adaptations to environmental challenges and recurring throughout evolutionary history rather than merely in plants.

Alternative splicing dynamically regulates common carp embryogenesis under thermal stress

BACKGROUND

Thermal stress is a major environmental factor affecting fish development and survival. Common carp (Cyprinus carpio) are susceptible to heat stress in their embryonic and larval phases, but the thermal stress response of alternative splicing during common carp embryogenesis remains poorly understood.

RESULTS

Using RNA-seq data from eight developmental stages and four temperatures, we constructed a comprehensive profile of alternative splicing (AS) during the embryogenesis of common carp, and found that AS genes and events are widely distributed among all stages. A total of 5,835 developmental stage-specific AS (SAS) genes, 21,368 temperature-specific differentially expressed genes (TDEGs), and 2,652 temperature-specific differentially AS (TDAS) genes were identified. Hub TDAS genes in each developmental stage, such as taf2, hnrnpa1, and drg2, were identified through protein-protein interaction (PPI) network analysis. The early developmental stages may be more sensitive to temperature, with thermal stress leading to a massive increase in the number of expressed transcripts, TDEGs, and TDAS genes in the morula stage, followed by the gastrula stage. GO and KEGG analyses showed that from the morula stage to the neurula stage, TDAS genes were more involved in intracellular transport, protein modification, and localization processes, while from the optic vesicle stage to one day post-hatching, they participated more in biosynthetic processes. Further subgenomic analysis revealed that the number of AS genes and events in subgenome B was generally higher than that in subgenome A, and the homologous AS genes were significantly enriched in basic life activity pathways, such as mTOR signaling pathway, p53 signaling pathway, and MAPK signaling pathway. Additionally, lncRNAs can play a regulatory role in the response to thermal stress by targeting AS genes such as lmnl3, affecting biological processes such as apoptosis and axon guidance.

CONCLUSIONS

In short, thermal stress can affect alternative splicing regulation during common carp embryogenesis at multiple levels. Our work complemented some gaps in the study of alternative splicing at both levels of embryogenesis and thermal stress in C. carpio and contributed to the comprehension of environmental adaptation formation in polyploid fishes during embryogenesis.

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.

Deletion of Asb15b gene can lead to a significant decrease in zebrafish intermuscular bone

Intermuscular bones, which are present in numerous economically significant fish species, have a negative impact on the development of aquaculture. The Asb15b gene, primarily expressed in skeletal muscle, plays a crucial role in regulating protein turnover and the development of muscle fibers. It stimulates protein synthesis and controls the differentiation of muscle fibers. In this study, we employed CRISPR/Cas9 technology to generate homozygous zebrafish strains with 7 bp and 49 bp deletions in the Asb15b gene. Subsequent analyses using skeleton staining demonstrated a substantial reduction in the number of intermuscular bones in adult Asb15b-/- -7 bp and Asb15b-/- -49 bp mutants compared to the wild-type zebrafish, with decreases of 30 % (P < 0.001) and 40 % (P < 0.0001), respectively. Histological experiments further revealed that the diameter and number of muscle fibers in adult Asb15b-/- mutants did not exhibit significant changes when compared to wild-type zebrafish. Moreover, qRT-PCR experiments demonstrated significant differences in the expression of bmp6 and runx2b genes, which are key regulators of intermuscular bone development, during different stages of intermuscular bone development in Asb15b-/- mutants. This study strongly suggests that the Asb15b gene plays a crucial role in regulating intermuscular bone development in fish and lays the groundwork for further exploration of the role of the Asb15b gene in zebrafish intermuscular bone development.

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.

A high-quality chromosome-level genome assembly of the topmouth culter (Culter alburnus Basilewsky, 1855)

Culter alburnus is extensively distributed in various rivers and lakes across China. As a widely adaptive fish species, it has significant economic values and special ecological roles. To meet research demands and provide better genomic resources, in this research, a chromosome-level genome assembly was constructed using HiFi long-reads and Hi-C sequencing data. Compared with the published versions, our genome assembly is of higher quality with only 31 gaps and closer to its true structure and sequence. The genome size was 1.052 Gb, with a contig N50 of 32.92 Mb and a scaffold N50 of 43.09 Mb. 55 contigs were anchored to 24 chromosomes on the basis of Hi-C data. A total of 598.23 Mb of repetitive sequences were annotated and 28,228 protein-coding genes were predicted. Additionally, BUSCO assessment indicated assembly and annotation scores of 98.3% and 99.2%, respectively. This high-quality genome will provide scientific support for excavating the species characteristics of C. alburnus and exploring its molecular mechanisms in response to environmental changes and stress.

Evolution of two-pore domain potassium channels and their gene expression in zebrafish embryos

BACKGROUND

The two-pore domain potassium (K2P) channels are a major type of potassium channels that maintain the cell membrane potential by conducting passive potassium leak currents independent of voltage change. They play prominent roles in multiple physiological processes, including neuromodulation, perception of pain, breathing and mood control, and response to volatile anesthetics. Mutations in K2P channels have been linked to many human diseases, such as neuronal and cardiovascular disorders and cancers. Significant progress has been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and function during embryonic development remain largely unknown.

RESULTS

We employed the zebrafish model and identified 23 k2p genes using BLAST search and gene cloning. We first analyzed vertebrate K2P channel evolution by phylogenetic and syntenic analyses. Our data revealed that the six subtypes of the K2P genes have already evolved in invertebrates long before the emergence of vertebrates. Moreover, the vertebrate K2P gene number increased, most likely due to two whole-genome duplications. Furthermore, we examined zebrafish k2p gene expression during early embryogenesis by in situ hybridization. Each subgroup's genes showed similar but distinct gene expression domains with some exceptions. Most of them were expressed in neural tissues consistent with their known function of neural excitability regulation. However, a few k2p genes were expressed temporarily in specific tissues or organs, suggesting that these K2P channels may be needed for embryonic development.

CONCLUSIONS

Our phylogenetic and developmental analyses of K2P channels shed light on their evolutionary history and potential roles during embryogenesis related to their physiological functions and human channelopathies.

Genetic Diversity and Population Structure of Rhodeus uyekii in the Republic of Korea Revealed by Microsatellite Markers from Whole Genome Assembly

This study is the first report to characterize the Rhodus uyekii genome and study the development of microsatellite markers and their markers applied to the genetic structure of the wild population. Genome assembly was based on PacBio HiFi and Illumina HiSeq paired-end sequencing, resulting in a draft genome assembly of R. uyekii. The draft genome was assembled into 2652 contigs. The integrity assessment of the assemblies indicates that the quality of the draft assemblies is high, with 3259 complete BUSCOs (97.2%) in the database of Verbrata. A total of 31,166 predicted protein-coding genes were annotated in the protein database. The phylogenetic tree showed that R. uyekii is a close but distinct relative of Onychostoma macrolepis. Among the 10 fish genomes, there were significant gene family expansions (8-2387) and contractions (16-2886). The average number of alleles amplified by the 21 polymorphic markers ranged from 6 to 23, and the average PIC value was 0.753, which will be useful for evolutionary and genetic analysis. Using population genetic analysis, we analyzed genetic diversity and the genetic structures of 120 individuals from 6 populations. The average number of alleles per population ranged from 7.6 to 9.9, observed heterozygosity ranged from 0.496 to 0.642, and expected heterozygosity ranged from 0.587 to 0.783. Discriminant analysis of principal components According to the analysis method, the population was divided into three populations (BS vs. DC vs. GG, GC, MS, DC). In conclusion, our study provides a useful resource for comparative genomics, phylogeny, and future population studies of R. uyekii.

The Expansion of Sirtuin Gene Family in Gilthead Sea Bream (Sparus aurata)-Phylogenetic, Syntenic, and Functional Insights across the Vertebrate/Fish Lineage

The Sirtuin (SIRT1-7) family comprises seven evolutionary-conserved enzymes that couple cellular NAD availability with health, nutrition and welfare status in vertebrates. This study re-annotated the sirt3/5 branch in the gilthead sea bream, revealing three paralogues of sirt3 (sirt3.1a/sirt3.1b/sirt3.2) and two of sirt5 (sirt5a/sirt5b) in this Perciform fish. The phylogeny and synteny analyses unveiled that the Sirt3.1/Sirt3.2 dichotomy was retained in teleosts and aquatic-living Sarcopterygian after early vertebrate 2R whole genome duplication (WGD). Additionally, only certain percomorphaceae and gilthead sea bream showed a conserved tandem-duplicated synteny block involving the mammalian-clustered sirt3.1 gene (psmd13-sirt3.1a/b-drd4-cdhr5-ctsd). Conversely, the expansion of the Sirt5 branch was shaped by the teleost-specific 3R WGD. As extensively reviewed in the literature, human-orthologues (sirt3.1/sirt5a) showed a high, conserved expression in skeletal muscle that increased as development advanced. However, recent sirt3.2 and sirt5b suffered an overall muscle transcriptional silencing across life, as well as an enhanced expression on immune-relevant tissues and gills. These findings fill gaps in the ontogeny and differentiation of Sirt genes in the environmentally adaptable gilthead sea bream, becoming a good starting point to advance towards a full understanding of its neo-functionalization. The mechanisms originating from these new paralogs also open new perspectives in the study of cellular energy sensing processes in vertebrates.

Database construction and comparative genomics analysis of genes involved in nutritional metabolic diseases in fish

Nutritional metabolic diseases in fish frequently arise in the setting of intensive aquaculture. The etiology and pathogenesis of these conditions involve energy metabolic disorders influenced by both internal genetic factors and external environmental conditions. The exploration of genes associated with nutritional and metabolic disorder has sparked considerable interest within both the aquaculture scientific community and the industry. High-throughput sequencing technology offers researchers extensive genetic information. Effectively mining, analyzing, and securely storing this data is crucial, especially for advancing disease prevention and treatment strategies. Presently, the exploration and application of gene databases concerning nutritional and metabolic disorders in fish are at a nascent stag. Therefore, this study focused on the model organism zebrafish and five primary economic fish species as the subjects of investigation. Using information from KEGG, OMIM, and existing literature, a novel gene database associated with nutritional metabolic diseases in fish was meticulously constructed. This database encompassed 4583 genes for Danio rerio, 6287 for Cyprinus carpio, 3289 for Takifugu rubripes, 3548 for Larimichthys crocea, 3816 for Oreochromis niloticus, and 5708 for Oncorhynchus mykiss. Through a comparative systems biology approach, we discerned a relatively high conservation of genes linked to nutritional metabolic diseases across these fish species, with over 54.9 % of genes being conserved throughout all six species. Additionally, the analysis pinpointed the existence of 13 species-specific genes within the genomes of large yellow croaker, tilapia, and rainbow trout. These genes exhibit the potential to serve as novel candidate targets for addressing nutritional metabolic diseases.

Genome-wide identification and expression analysis of Sox gene family in the Manila clam (Ruditapes philippinarum)

Sox transcription factors are vital in numerous fundamental biological processes. In this study, nine Sox gene family members were discovered in the Ruditapes philippinarum genome, classified into the SoxB1, SoxB2, SoxC, SoxD, SoxE, and SoxF groups, marking the first genome-wide identification of this gene family in R. philippinarum. Analyses of phylogeny, exon-intron structures, and domains bolster the support for their categorization and annotation. Furthermore, transcriptomic analyses across various developmental stages revealed that RpSox4, RpSox5, RpSox9, and RpSox11 were significantly expressed in the D-larval stage. Additionally, investigations into transcriptomes of clams with different shell colors indicated that most sox genes exhibited their highest expression levels in orange clams, followed by zebra, white zebra, and white clams, and the results of transcriptomes analysis in different tissues indicated that 8 Sox genes (except RpSox17) were highly expressed in the mantle tissue. Moreover, qPCR was used to detect the expression of Sox gene in R. philippinarum at different developmental periods, different shell colors and different tissues, and the results showed consistency with those of the transcriptomes. This study's findings lay the groundwork for additional exploration into the role of the Sox gene in melanin production in R. philippinarum shells.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

The reference genome of Macropodus opercularis (the paradise fish)

Amongst fishes, zebrafish (Danio rerio) has gained popularity as a model system over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species has a highly complex behavioral repertoire and has been the subject of many ethological investigations but lacks genomic resources. Here we report the reference genome assembly of M. opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of 483,077,705 base pairs (~483 Mb) on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ~90% of them to orthogroups.

Whole-genome Sequencing Reveals Autooctoploidy in Chinese Sturgeon and Its Evolutionary Trajectories

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Secreted novel AID/APOBEC-like deaminase 1 (SNAD1) - a new important player in fish immunology

The AID/APOBECs are a group of zinc-dependent cytidine deaminases that catalyse the deamination of bases in nucleic acids, resulting in a cytidine to uridine transition. Secreted novel AID/APOBEC-like deaminases (SNADs), characterized by the presence of a signal peptide are unique among all of intracellular classical AID/APOBECs, which are the central part of antibody diversity and antiviral defense. To date, there is no available knowledge on SNADs including protein characterization, biochemical characteristics and catalytic activity. We used various in silico approaches to define the phylogeny of SNADs, their common structural features, and their potential structural variations in fish species. Our analysis provides strong evidence of the universal presence of SNAD1 proteins/transcripts in fish, in which expression commences after hatching and is highest in anatomical organs linked to the immune system. Moreover, we searched published fish data and identified previously, "uncharacterized proteins" and transcripts as SNAD1 sequences. Our review into immunological research suggests SNAD1 role in immune response to infection or immunization, and interactions with the intestinal microbiota. We also noted SNAD1 association with temperature acclimation, environmental pollution and sex-based expression differences, with females showing higher level. To validate in silico predictions we performed expression studies of several SNAD1 gene variants in carp, which revealed distinct patterns of responses under different conditions. Dual sensitivity to environmental and pathogenic stress highlights its importance in the fish and potentially enhancing thermotolerance and immune defense. Revealing the biological roles of SNADs represents an exciting new area of research related to the role of DNA and/or RNA editing in fish biology.

Subgenomic T cell receptor alpha and delta (TRA/TRD) loci in common carp

In jawed vertebrates, the T cell receptor alpha (TRA) and delta (TRD) genes, which encode the TRα and TRδ chains, respectively, are located as a nested structure on a single chromosome. To date, no animal has been reported to harbor multiple TRA/TRD loci on different chromosomes. Therefore, herein, we describe the first full annotation of the TRA/TRD genomic regions of common carp, an allo-tetraploid fish species that experiences cyprinid-specific whole-genome duplication (WGD) in evolution. Fine genomic maps of TRA/TRD genomic regions 1 and 2, on LG30 and LG22, respectively, were constructed using the annotations of complete sets of TRA and TRD genes, including TRA/TRD variable (V), TRA junction (J), and constant (C), TRD diversity (D), and the J and C genes. The structure and synteny of the TRA/TRD genomic regions were highly conserved in zebrafish, indicating that these regions are on individual chromosomes. Furthermore, analysis of the variable regions of the TRA and TRD genes in a monoclonal T cell line revealed that both subgenomic regions 1 and 2 were indeed rearranged. Although carp TRAV and TRDV genes were phylogenetically divided into different lineages, they were mixed and organized into the TRA/TRD V gene clusters on the genome, similar to that in other vertebrates. Notably, 285 potential TRA/TRD V genes were detected in the TRA/TRD genomic regions, which is the most abundant number of genes in vertebrates and approximately two-fold that in zebrafish. The recombination signal sequences (RSSs) at the end of each V gene differed between TRAV and TRDV, suggesting that RSS variations might separate each V gene into a TRα or TRδ chain. This study is the first to describe subgenomic TRA/TRD loci in animals. Our findings provide fundamental insights to elucidate the impact of WGD on the evolution of immune repertoire.

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.

Genome-wide identification and structural analysis of the BMP gene family in Triplophysa dalaica

BACKGROUND

Bone morphogenetic proteins (BMPs) are part of the transforming growth factor beta (TGF-β) superfamily and play crucial roles in bone development, as well as in the formation and maintenance of various organs. Triplophysa dalaica, a small loach fish that primarily inhabits relatively high elevations and cooler water bodies, was the focus of this study. Understanding the function of BMP genes during the morphogenesis of T. dalaica helps to clarify the mechanisms of its evolution and serves as a reference for the study of BMP genes in other bony fishes. The data for the T. dalaica transcriptome and genome used in this investigation were derived from the outcomes of our laboratory sequencing.

RESULTS

This study identified a total of 26 BMP genes, all of which, except for BMP1, possess similar TGF-β structural domains. We conducted an analysis of these 26 BMP genes, examining their physicochemical properties, subcellular localization, phylogenetic relationships, covariance within and among species, chromosomal localization, gene structure, conserved motifs, conserved structural domains, and expression patterns. Our findings indicated that three BMP genes were associated with unstable proteins, while 11 BMP genes were located within the extracellular matrix. Furthermore, some BMP genes were duplicated, with the majority being enriched in the GO:0008083 pathway, which is related to growth factor activity. It was hypothesized that genes within the BMP1/3/11/15 subgroup (Group I) play a significant role in the growth and development of T. dalaica. By analyzing the expression patterns of proteins in nine tissues (gonad, kidney, gill, spleen, brain, liver, fin, heart, and muscle), we found that BMP genes play diverse regulatory roles during different stages of growth and development and exhibit characteristics of division of labor.

CONCLUSIONS

This study contributes to a deeper understanding of BMP gene family member expression patterns in high-altitude, high-salinity environments and provides valuable insights for future research on the BMP gene family in bony fishes.

De Novo Genome Assembly of the Whitespot Parrotfish (Scarus forsteni): A Valuable Scaridae Genomic Resource

Scarus forsteni, a whitespot parrotfish from the Scaridae family, is a herbivorous fish inhabiting coral reef ecosystems. The deterioration of coral reefs has highly affected the habitats of the parrotfish. The decline in genetic diversity of parrotfish emphasizes the critical importance of conserving their genetic variability to ensure the resilience and sustainability of marine ecosystems for future generations. In this study, a genome of S. forsteni was assembled de novo through using Illumina and Nanopore sequencing. The 1.71-Gb genome of S. forsteni, was assembled into 544 contigs (assembly level: contig). It exhibited an N50 length of 17.97 Mb and a GC content percentage of 39.32%. Our BUSCO analysis revealed that the complete protein of the S. forsteni genome had 98.10% integrity. Combined with structure annotation data, 34,140 (74.81%) genes were functionally annotated out of 45,638 predicted protein-coding genes. Upon comparing the genome size and TE content of teleost fishes, a roughly linear relationship was observed between these two parameters. However, TE content is not a decisive factor in determining the genome size of S. forsteni. Population history analysis results indicate that S. forsteni experienced two major population expansions, both of which occurred before the last interglacial period. In addition, through a comparative genomic analysis of the evolutionary relationship of other species, it was found that S. forsteni had the closest relationship with Cheilinus undulatus, another member of the Labridae family. Our expansion and contraction analysis of the gene family showed that the expansion genes were mainly associated with immune diseases, organismal systems, and cellular processes. At the same time, cell transcription and translation, sex hormone regulation, and other related pathways were also more prominent in the positive selection genes. The genomic sequence of S. forsteni offers valuable resources for future investigations on the conservation, evolution, and behavior of fish species.

Contribution of elovl5a to Docosahexaenoic Acid (DHA) Synthesis at the Transcriptional Regulation Level in Common Carp, Cyprinus carpio

Docosahexaenoic acid (DHA) is an essential nutrient for humans and plays a critical role in human development and health. Freshwater fish, such as the common carp (Cyprinus carpio), have a certain degree of DHA biosynthesis ability and could be a supplemental source of human DHA needs. The elongase of very-long-chain fatty acid 5 (Elovl5) is an important enzyme affecting polyunsaturated fatty acid (PUFA) biosynthesis. However, the function and regulatory mechanism of the elovl5 gene related to DHA synthesis in freshwater fish is not clear yet. Previous studies have found that there are two copies of the elovl5 gene, elovl5a and elovl5b, which have different functions. Our research group found significant DHA content differences among individuals in Yellow River carp (Cyprinus carpio var.), and four candidate genes were found to be related to DHA synthesis through screening. In this study, the expression level of elovl5a is decreased in the high-DHA group compared to the low-DHA group, which indicated the down-regulation of elovl5a in the DHA synthesis pathways of Yellow River carp. In addition, using a dual-luciferase reporter gene assay, we found that by targeting the 3'UTR region of elovl5a, miR-26a-5p could regulate DHA synthesis in common carp. After CRISPR/Cas9 disruption of elovl5a, the DHA content in the disrupted group was significantly higher than in the wildtype group; meanwhile, the expression level of elovl5a in the disrupted group was significantly reduced compared with the wildtype group. These results suggest that elovl5a may be down-regulating DHA synthesis in Yellow River carp. This study could provide useful information for future research on the genes and pathways that affect DHA synthesis.

Asymmetric and parallel subgenome selection co-shape common carp domestication

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Effect of aqueous extract of Millettia speciosa Champ on intestinal health maintenance and immune enhancement of Cyprinus carpio

Millettia speciosa Champ (MSP) is a natural Chinese herb that improves gastrointestinal health and enhances animal immunity. An 8-week feeding trial with different MSP levels (0, 150, 300, and 600 mg/kg) was conducted to evaluate the promotive effects of MSP in Cyprinus carpio. Results indicate that MSP improved intestinal immunity to some extent evidenced by the immuno-antioxidant parameters and the 16S rRNA in the Illumina MiSeq platform. With the analysis of transcriptome sequencing, 4685 differentially expressed genes (DEGs) were identified, including 2149 up-regulated and 2536 down-regulated. According to the GO and KEGG enrichments, DEGs were mainly involved in the immune system. Transcriptional expression of the NOD-like signaling pathway and key genes retrieved from the transcriptome database confirmed that innate immunity was improved in response to dietary MSP administration. Therefore, MSP could be used as a feed supplement that enhances immunity. This may provide insight into Chinese herb additive application in aquaculture production.

Role of HIF in fish inflammation

The hypoxia-inducing factor (HIF) is a central transcription factor in cellular oxygen sensing and regulation. It is common that the inflammation always appears in many diseases, like infectious diseases in fishes, and the inflammation is often accompanied by hypoxia, as a hallmark of inflammation. Besides coordinating cellular responses to low oxygen, HIF-mediated hypoxia signaling pathway is also crucial for immune responses such as the regulations of innate immune cell phenotype and function, as well as metabolic reprogramming under the inflammation. However, the understanding of the molecular mechanisms by which HIFs regulate the inflammatory response in fish is still very limited. Here, we review the characteristics of HIF as well as its roles in innate immune cells and the infections caused by bacteria and viruses. The regulatory effects of HIF on the metabolic reprogramming of innate immune cells are also discussed and the future research directions are outlooked. This paper will serve as a reference for elucidating the molecular mechanism of HIF regulating inflammation and identifying treatment strategies to target HIF for fish disease.

Differential Expression of miRNAs, lncRNAs, and circRNAs between Ovaries and Testes in Common Carp (Cyprinus carpio)

Female common carp grow faster than male individuals, implying that rearing females could be more profitable in aquaculture. Non-coding RNAs (ncRNAs) serve as versatile regulators with multiple functions in diverse biological processes. However, the roles of ncRNAs in the sex differentiation of common carp are less studied. In this study, we investigated the expression profiles of ncRNAs, including miRNAs, lncRNAs, and circRNAs, in the gonads to comprehend the roles of ncRNAs in sex differentiation in common carp. A substantial number of differentially expressed (DE) ncRNAs in ovaries and testes were identified. Some miRNAs, notably miR-205, miR-214, and miR-460-5p, might modulate hormone synthesis and thus maintain sex. A novel miRNA, novel_158, was predicted to suppress the expression of foxl3. DE lncRNAs were associated with oocyte meiosis, GnRH signaling pathways, and steroid biosynthesis, while DE circRNA target genes were enriched in the WNT signaling pathway and MAPK signaling pathway. We also analyzed ncRNA-mRNA interactions to shed light on the crosstalk between competing endogenous RNAs (ceRNAs), which is the critical mechanism by which lncRNAs and circRNAs function. Some lncRNAs and circRNAs may be able to competitively bind novel_313, a new miRNA, and thus regulate hsd17β3. Our research will provide a valuable resource for understanding the genetic basis of gonadal differentiation and development in common carp.

Transcriptome analysis in hepatopancreases reveals the response of domesticated common carp to a high-temperature environment in the agricultural heritage rice-fish system

Qingtian paddy field carp (PF-carp) is a local carp cultivated in the paddy field of Qingtian, Zhejiang. This rice-fish co-culture system has been recognized as one of the Globally Important Agriculture Heritage Systems (GIAHS). PF-carp has been acclimatized to the high-temperature environment of shallow paddy fields after several centuries of domestication. To reveal the physiological and molecular regulatory mechanisms of PF-carp, we chose to use 28°C as the control group and 34°C as the treatment group. We measured biochemical parameters in their serum and hepatopancreases and also performed transcriptome sequencing analysis. Compared with the control group, the serum levels of malondialdehyde (MDA), glucose (GLU), glutathione peroxidase (GSH-Px), catalase (CAT), alanine aminotransferase (ALT), and aspartate aminotransferase (AST) show no significant change. In addition, superoxide dismutase (SOD), GSH-Px, and CAT also show no significant change in hepatopancreases. We identified 1,253 differentially expressed genes (DEGs), and their pathway analysis revealed that heat stress affected AMPK signaling pathway, protein export, and other biological processes. It is worth noting that protein processing in the endoplasmic reticulum (ER) was the most significantly enriched pathway identified by the Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis (GSEA). Significantly higher levels of HSP40, HSP70, HSP90, and other ubiquitin ligase-related genes were upregulated. In summary, heat stress did not lead to tissue damage, inflammation, oxidative stress, and ER stress in the hepatopancreases of PF-carp. This study provides valuable insights into the adaptation mechanism of this species to the high-temperature environment of paddy fields.

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

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

Chromosome-level genome assembly and whole-genome resequencing of topmouth culter (Culter alburnus) provide insights into the intraspecific variation of its semi-buoyant and adhesive eggs

Topmouth culter (Culter alburnus) is an ecologically and economically important species belonging to the subfamily Culterinae that is native to and widespread in East Asia. Intraspecific variation of semi-buoyant and adhesive eggs in topmouth culter provides an ideal opportunity to investigate the genetic mechanisms of spawning habits underlying the adaptive radiation of cyprinids in East Asia. In this study, we present a chromosome-level genome assembly of topmouth culter and re-sequenced 158 individuals from six locations in China covering three geographical groups and two egg type variations. The topmouth culter genome size was 1.05 Gb, with a contig N50 length of 17.8 Mb and anchored onto 24 chromosomes. Phylogenetic analysis showed that the divergence time of the Culterinae was coinciding with the time of initiation of the Asian monsoon intensification. Gene family evolutionary analysis indicated that the expanded gene families in topmouth culter were associated with dietary adaptation. Population-level genetic analysis indicated clear differentiation among the six populations, which were clustered into three distinct clusters, consistent with their geographical divergence. The historical effective population size of topmouth culter correlated with the Tibetan Plateau uplifting according to the demographic history reconstruction. A selective sweep analysis between adhesive and semi-buoyant egg populations revealed the genes associated with the hydration and adhesiveness of eggs, indicating divergent selection towards different hydrological environments. This study offers a high-resolution genetic resource for further studies on evolutionary adaptation, genetic breeding and conservation of topmouth culter, providing insights into the molecular mechanisms for egg type variation of East Asian cyprinids.

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

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

The food safety assessment of all-female common carp (Cyprinus carpio) (cyp17a1+/-;XX genotype) generated using genome editing technology

There are several technical challenges and public issues concerning genome editing applications before they become viable in commercial aquaculture. Recently, we developed a novel strategy to generate all-female (AF) common carp, which exhibited a growth advantage over the control carp, using genetic editing through single gene-targeting manipulation. Here, we found that the body weight of the AF common carp was higher by 22.58% than that of the control common carp. Because the genotype of the AF common carp was cyp17a1+/-;XX, the contents of sex steroids were normally synthesized, as they were comparable to that of the control female carp. To evaluate the food safety of the AF carp, Wistar rats were fed a diet containing control female carp (control, C) or all-female (AF) carp at an incorporation rate of 5, 10 and 20% (w/w) for 90 days. Compared with those fed control carp, the rats fed AF common carp exhibited no significant difference in body weight, food intake, feed conversion ratio, hematology, serum biochemistry, urine test, relative organ weight, gross necropsy, and histopathological examination. This is the first food safety assessment of the farmed fish strain cultured using CRISPR/Cas9, which will further advance the fishery development of genome-edited animals.

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

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

Morphological, histological and phylogenetic identification of three species of Myxobolus (Myxosporea: Myxobolidae) infecting different carp lineages in China

The global cultivation of common carp Cyprinus carpio is developed primarily for either food or recreational purposes and myxosporean infections causing significant economic losses are regularly troublesome for carp farmers. However, most myxosporean species are poorly characterized, making it difficult to correctly elaborate on their parasitism mode and determine pathogenicity. Here, based on an integrative approach, we supplemented fundamental information on three Myxobolus species infecting different carp lineages in China. Myxobolus elliptoides Wu et Chen, 1987; Nine new species of Myxosporida from freshwater fish of Wuhu lake, Hubei, China. Acta Hydrobiologica Sinica, 11, 161 formed yellowish pseudocysts in the anal fin of common carp and were observed with ellipsoidal myxospores and two polar capsules of distinct dimensions. The dermal tissue within the fin ray was regarded as the type site for growing plasmodia of M. elliptoides. Myxobolus basilamellaris Lom et Molnar, 1983; Myxobolus basilamellaris sp. n. (myxozoa, Myxosporea), a parasite of the gills of common carp (Cyprinus carpio L.). Folia Parasitologica, 30, 1 aggregated ellipsoidal myxospores at the base of the gill filaments of koi carp juvenile, causing significant tissue damage. Myxobolus artus Achmerov, 1960, enwrapped in diffuse pseudocysts, was observed throughout the trunk muscle of mirror carp. Its oblate myxospores were asymmetrical and contained two unequal-sized polar capsules. According to the localization of large plasmodia, inter- and intramuscular preferences were displayed by M. artus. For precise species identification, we provided the SSU rDNA sequences for each species of Myxobolus. Among them, M. elliptoides was molecularly characterized for the first time, showing the highest 94.21% identity to Thelohanellus sinensis (KY469292). Phylogenetically, the affinity of both M. artus and M. basilamellaris to their conspecific species derived from different carp lineages was highly supported. Moreover, the intermixed cluster of Myxobolus species, including M. elliptoides, with those of Thelohanellus species provided evidence for querying the monophyletic evolution history of these taxa. This work updates the appreciation of the concerned species and enhances our understanding of the parasite fauna of common carp.

The reference genome of the paradise fish (Macropodus opercularis)

Over the decades, a small number of model species, each representative of a larger taxa, have dominated the field of biological research. Amongst fishes, zebrafish (Danio rerio) has gained popularity over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species from Southeast Asia, has a highly complex behavioral repertoire and has been the subject of many ethological investigations, but lacks genomic resources. Here we report the reference genome assembly of Macropodus opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of ≈483 Mb on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ≈90% of them to orthogroups. Completeness analysis showed that 98.5% of the Actinopterygii core gene set (ODB10) was present as a complete ortholog in our reference genome with a further 1.2 % being present in a fragmented form. Additionally, we cloned multiple genes important during early development and using newly developed in situ hybridization protocols, we showed that they have conserved expression patterns.

The vertebrate sialylation machinery: structure-function and molecular evolution of GT-29 sialyltransferases

Every eukaryotic cell is covered with a thick layer of complex carbohydrates with essential roles in their social life. In Deuterostoma, sialic acids present at the outermost positions of glycans of glycoconjugates are known to be key players in cellular interactions including host-pathogen interactions. Their negative charge and hydrophilic properties enable their roles in various normal and pathological states and their expression is altered in many diseases including cancers. Sialylation of glycoproteins and glycolipids is orchestrated by the regulated expression of twenty sialyltransferases in human tissues with distinct enzymatic characteristics and preferences for substrates and linkages formed. However, still very little is known on the functional organization of sialyltransferases in the Golgi apparatus and how the sialylation machinery is finely regulated to provide the ad hoc sialome to the cell. This review summarizes current knowledge on sialyltransferases, their structure-function relationships, molecular evolution, and their implications in human biology.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information available, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprising 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system (immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes) exhibit expansion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum when compared to mammals. Our functional studies suggest that lack of the gc gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi, as well as the genome evolution in cartilaginous fishes.

Transposon signatures of allopolyploid genome evolution

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