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

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.

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

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

Genome size estimation and its associations with body length, chromosome number and evolution in teleost fishes

Precise estimation of genome size (GS) is vital for various genomic studies, such as deciding genome sequencing depth, genome assembly, biodiversity documentation, evolution, genetic disorders studies, duplication events etc. Animal Genome Size Database provides GS of over 2050 fish species, which ranges from 0.35 pg in pufferfish (Tetraodon nigroviridis) to 132.83 pg in marbled lungfish (Protopterus aethiopicus). The GS of majority of the fishes inhabiting waters of Indian subcontinent are still missing. In present study, we estimated GS of 51 freshwater teleost (31 commercially important, 7 vulnerable and 13 ornamental species) that ranged from 0.58 pg in banded gourami (Trichogaster fasciata) to 1.92 pg in scribbled goby (Awaous grammepomus). Substantial variation in GS was observed within the same fish orders (0.64-1.45 pg in cypriniformes, 0.70-1.41 pg in siluriformes and 0.58-1.92 pg in perciformes). We examined the relationship between the GS, chromosome number and body length across all the fishes. Body length was found to be associated with GS, whereas no relationship was noticed between the GS and the chromosome number. The analysis using ancestral information revealed haploid chromosome number 25, 27 and 24 for the most recent common ancestor of cypriniformes, siluriformes and perciformes, respectively. The study led to generation of new records on GS of 43 fish species and revalidated records for 8 species. The finding is valuable resource for further research in the areas of fish genomics, molecular ecology and evolutionary conservation genetics.

High production of transfer RNAs identifies the presence of developing oocytes in ovaries and intersex testes of teleost fish

5S rRNA is highly transcribed in fish oocytes and this transcription levels can be used to identify the presence of oocytes in the intersex testes of fish exposed to xenoestrogens. Similar to 5S rRNA, tRNAs are transcribed by RNA polymerase III (Pol-III) in eukaryotes, so this study focuses in the analysis of the levels of expression of tRNAs in the gonads (ovaries and testes) of eight teleost species as a possible new oocyte molecular marker. Total RNA extracted from gonads of six commercial teleost species in the Biscay Bay, from the pollution sentinel species thicklip grey mullet (Chelon labrosus) known present intersex testes in response to xenoestrogens in Gernika estuary and from the laboratory model species Danio rerio were analysed through capillary electrophoresis. Bioanalyzer electropherograms were used to quantify the concentrations of tRNAs, 5S and 5.8S rRNA. All studied ovaries expressed significantly higher levels of tRNAs and 5S rRNA than testes. A tRNA to 5.8S rRNA index was calculated which differentiates ovaries from testes, and identifies some intersex testes in between testes and ovaries in mullets. The tRNA/5.8S ratio was highest in ovaries in previtellogenic stage, decreasing towards maturity. Thus, strong oocyte expression of tRNAs is an additional proof of high activity levels of Pol-III during early stages of oocyte development in teleost ovaries. Incidentally, we observed that miRNA concentrations were always higher in testes than ovaries. The indexing approach developed in the present study could have multiple applications in teleost reproduction research and in the development of early molecular markers of intersex condition.

Comparative genomics studies on the stk gene family in vertebrates: From the bighead carp (Hypophthalmichthys nobilis) genome

The mammalian sterile 20-like (MST) family belongs to the serine/threonine protein kinase (STK) superfamily and participates in a variety of biological processes, such as cell apoptosis, polarity, migration, immune regulation, inflammatory responses, and cancer. In the economically important bighead carp (Hypophthalmichthys nobilis), the STK gene family and immune-related biological functions may be helpful in increasing its economic yield. However, the comprehensive role of STKs in the bighead carp remains unclear. In this study, the five stk sequences from the bighead carp were divided into two classes: stk3/4 and stk24/25/26. Gene structure and motif prediction analyses confirmed that stk is conserved in the bighead carp. Compared to 26 other vertebrate species, teleosts (including bighead carp) possess more stk members because of teleost-specific whole-genome duplication. Synteny analysis revealed that stk3, stk24, stk25, and stk26 have been relatively conserved in bighead carp during evolution. Meanwhile, stk4 was lost in most Cyprinid species, including bighead carp, during evolution. RNA-seq data revealed that STK expression was associated with various pathogens, and the expression of these STKs (Hnstk3, Hnstk24a, Hnstk24b, Hnstk25, and Hnstk26) was different in seven tissues of bighead carp. In addition, we showed that STK expression levels were dramatically altered in the head kidney and that stk24 was involved in defense against Aeromonas hydrophila. This study provides a molecular basis for the analysis of stk function in bighead carp, and can be used as a reference for further phylogenomics.

Disruption of T-box transcription factor eomesa results in abnormal development of median fins in Oujiang color common carp Cyprinus carpio

Median fins are thought to be ancestors of paired fins which in turn give rise to limbs in tetrapods. However, the developmental mechanisms of median fins remain largely unknown. Nonsense mutation of the T-box transcription factor eomesa in zebrafish results in a phenotype without dorsal fin. Compared to zebrafish, the common carp undergo an additional round of whole genome duplication, acquiring an extra copy of protein-coding genes. To verify the function of eomesa genes in common carp, we established a biallelic gene editing technology in this tetraploidy fish through simultaneous disruption of two homologous genes, eomesa1 and eomesa2. We targeted four sites located upstream or within the sequences encoding the T-box domain. Sanger sequencing data indicated the average knockout efficiency was around 40% at T1-T3 sites and 10% at T4 site in embryos at 24 hours post fertilization. The individual editing efficiency was high to about 80% at T1-T3 sites and low to 13.3% at T4 site in larvae at 7 days post fertilization. Among 145 mosaic F0 examined at four months old, three individuals (Mutant 1-3) showed varying degrees of maldevelopment in the dorsal fin and loss of anal fin. Genotyping showed the genomes of all three mutants were disrupted at T3 sites. The null mutation rates on the eomesa1 and eomesa2 loci were 0% and 60% in Mutant 1, 66.7% and 100% in Mutant 2, and 90% and 77.8% in Mutant 3, respectively. In conclusion, we demonstrated a role of eomesa in the formation and development of median fins in Oujiang color common carp and established an method that simultaneously disrupt two homologous genes with one gRNA, which would be useful in genome editing in other polyploidy fishes.

Navigating the CoGe Online Software Suite for Polyploidy Research

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloaded further. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.

Comparative transcriptome analysis of skin color-associated genes in leopard coral grouper (Plectropomus leopardus)

BACKGROUND

The leopard coral grouper (Plectropomus leopardus) is an important economic species in East Asia-Pacific countries. To meet the market demand, leopard coral grouper is facing overfishing and their population is rapidly declining. With the improvement of the artificial propagation technique, the leopard coral grouper has been successfully cultured by Fisheries Research Institute in Taiwan. However, the skin color of farmed individuals is often lacking bright redness. As such, the market price of farmed individuals is lower than wild-type.

RESULTS

To understand the genetic mechanisms of skin coloration in leopard coral grouper, we compared leopard coral grouper with different skin colors through transcriptome analysis. Six cDNA libraries generated from wild-caught leopard coral grouper with different skin colors were characterized by using the Illumina platform. Reference-guided de novo transcriptome data of leopard coral grouper obtained 24,700 transcripts, and 1,089 differentially expressed genes (DEGs) were found between red and brown skin color individuals. The results showed that nine candidate DEGs (epha2, sema6d, acsl4, slc7a5, hipk1, nol6, timp2, slc25a42, and kdf1) significantly associated with skin color were detected by using comparative transcriptome analysis and quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The findings may provide genetic information for further skin color research, and to boost the market price of farmed leopard coral grouper by selective breeding.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

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

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

Negative regulation of interferon regulatory factor 6 (IRF6) in interferon and NF-κB signalling pathways of common carp (Cyprinus carpio L.)

BACKGROUND

Interferon (IFN) regulatory factors (IRFs) is a kind of transcription factors, which play an important role in regulating the expression of type I IFN and related genes. In mammals, IRF6 is not relevant with IFN expression, while zebrafish IRF6 was reported to be a positive regulator of IFN expression and could be phosphorylated by both MyD88 and TBK1. However, the role of IRF6 in the immune response and IFN transcription of common carp is unknown.

RESULTS

In the present study, the cDNA of IRF6 gene (CcIRF6) was cloned from common carp using RACE technique, with a total length of 1905 bp, encoding 471 amino acid residues, which possesses two functional domains of DBD and IAD. Similarity analysis showed that CcIRF6 had more than 50% similarity with IRFs of other vertebrates, and had the highest similarity with grass carp and zebrafish, among which the DBD domain was much more conserved. The phylogenetic analysis showed that CcIRF6 is in the branch of Osteichthyes and has the closest relationship with grass carp. In healthy common carp, the CcIRF6 was expressed in all the examined tissues, with the highest level in the oral epithelium, and the lowest level in the head kidney. After intraperitoneal injection of poly(I:C) or Aeromonas hydrophila, the expression of CcIRF6 increased in spleen, head kidney, foregut and hindgut of common carp. Moreover, poly(I:C), LPS, PGN and flagellin induced the expression of CcIRF6 in peripheral leukocytes and head kidney leukocytes of common carp in vitro. In EPC cells, CcIRF6 inhibited the expression of some IFN-related genes and pro-inflammatory cytokines, and dual luciferase reporter assay showed that CcIRF6 reduced the activity of IFN and NF-κB reporter genes.

CONCLUSIONS

The present study suggests that CcIRF6 is involved in the antiviral and antibacterial immune response of common carp, and negatively regulate the expression of IFN and NF-κB signalling pathways, which provides a theoretical basis for the study and prevention of fish disease pathogenesis.

A high-quality chromosome-level genome assembly of Pelteobagrus vachelli provides insights into its environmental adaptation and population history

Pelteobagrus vachelli is a freshwater fish with high economic value, but the lack of genome resources has severely restricted its industrial development and population conservation. Here, we constructed the first chromosome-level genome assembly of P. vachelli with a total length of approximately 662.13 Mb and a contig N50 was 14.02 Mb, and scaffolds covering 99.79% of the assembly were anchored to 26 chromosomes. Combining the comparative genome results and transcriptome data under environmental stress (high temperature, hypoxia and Edwardsiella. ictaluri infection), the MAPK signaling pathway, PI3K-Akt signaling pathway and apelin signaling pathway play an important role in environmental adaptation of P. vachelli, and these pathways were interconnected by the ErbB family and involved in cell proliferation, differentiation and apoptosis. Population evolution analysis showed that artificial interventions have affected wild populations of P. vachelli. This study provides a useful genomic information for the genetic breeding of P. vachelli, as well as references for further studies on fish biology and evolution.

Genome-wide survey reveals the phylogenomic relationships of Chirolophisjaponicus Herzenstein, 1890 (Stichaeidae, Perciformes)

Fish are the largest vertebrate group, consisting of more than 30 000 species with important ecological and economical value, while less than 3% of fish genomes have been published. Herein, a fish, Chirolophisjaponicus, was sequenced using the next-generation sequencing. Approximately 595.7 megabase pair of the C.japonicus genome was assembled (49 901 contigs with 42.61% GC contents), leading to a prediction of 46 729 protein-coding gene models. A total of 554 136 simple sequence repeats was identified in the whole genome of C.japonicus, and dinucleotide microsatellite motifs were the most abundant, accounting for 59.49%. Phylogenomic analysis of 16 genomes based on the 694 single-copy genes suggests that C.japonicus is closely related with Anarrhichthysocellatus, Cebidichthysviolaceus, and Pholisgunnellus. The results provide more thorough genetic information of C.japonicus and a theoretical basis and reference for further genome-wide analysis.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.