Chromosome‐level genome assembly of a cyprinid fish Onychostoma macrolepis by integration of nanopore sequencing, Bionano and Hi‐C technology

Molecular characterization of an antimicrobial peptide LEAP-2 in Onychostoma macrolepis: Expression pattern, antimicrobial ability and immunomodulation function

Antimicrobial peptides (AMPs) as the host's antimicrobial source are promising candidates for the development of novel antibacterial agents. Our study identified and analyzed the LEAP-2 gene from the cavefish - Onychostoma macrolepis (OmLEAP-2), and antimicrobial activity of OmLEAP-2 was assessed through both in vivo and in vitro experiments. The OmLEAP-2 encodes 92 amino acids with a hydrophobic surface (LAMMPWY) and two disulfide bonds. Following A. hydrophila infection, OmLEAP-2 expression was up-regulated in the liver, spleen and intestine in Onychostoma macrolepis. Recombinant OmLEAP-2 protein exhibited dose-dependent antimicrobial activity against both Gram-negative and Gram-positive bacteria, particularly A. hydrophila and A. veronii. OmLEAP-2 overexpression not only reduced bacterial load and increased survival rate, but also suppressed expression of TNF-α and IL-1β in A. hydrophila-infected Onychostoma macrolepis. Furthermore, OmLEAP-2 regulated iron metabolism by increasing liver iron levels, altering serum iron content, and modulating iron-related genes (fpn1 and hepcidin) expression. Collectively, these findings demonstrate that OmLEAP-2 exhibits broad-spectrum antibacterial activity against aquatic pathogens and shows potential as a therapeutic peptide for treating bacterial disease.

Chromosome-scale assembly of the Xenocypris davidi using PacBio HiFi reads and Hi-C technologies

Xenocypris davidi is a benthic fish species widely distributed in the water systems south of the Yellow River in China, playing a significant role in aquatic ecosystems. Despite its ecological and economic importance, genomic resources for X. davidi are limited, hindering a comprehensive understanding of its evolutionary adaptations and genetic improvements. This study presents the first chromosome-level genome assembly of X. davidi, utilizing PacBio long-reads, Illumina short reads, and Hi-C sequencing data. The genome assembly spans 1.05 Gb with a scaffold N50 length of 33.99 Mb, and 95.12% of the genome sequence was successfully anchored onto 24 pseudochromosomes. We identified 27,360 protein-coding genes, of which 26,672 were functionally annotated. This genome sequence provides a valuable resource for exploring the molecular basis of agronomic traits in X. davidi and will facilitate its genetic enhancement.

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.

Chromosome-level genome assembly of the ratmouth barbel, Ptychidio jordani

The ratmouth barbel (Ptychidio jordani) is a critically endangered freshwater fish from the Cyprinidae family, primarily due to overfishing and habitat disruption. To address the challenges of its shrinking wild populations and the difficulties in artificial reproduction, we sequenced, assembled, and annotated a high-quality chromosome-level genome of P. jordani using next-generation short-read sequencing, third-generation long-read sequencing, and Hi-C sequencing. The final genome assembly was 1.14 Gb, consisting of 25 chromosomes with a contig N50 of 25.14 Mb and a scaffold N50 of 42.91 Mb. We identified 25,183 protein-coding genes, 751.75 Mb of repeats, and 19,373 ncRNAs. Methylation loci on most chromosomes ranged from 1,000 to 3,000 per 100 kb window. Gene expression levels across various tissues were analyzed, revealing 12,135 (caudal fin), 11,465 (liver), 14,438 (gill), 12,413 (heart), 8,301 (spleen), and 3,578 (kidney) differentially expressed genes compared to muscle. The comprehensive genomic and transcriptomic resources generated here will aid in understanding the ecology, adaptation, and environmental responses of P. jordani, supporting future research and conservation efforts.

Galectin-9 activates host immune response and improve immunoprotection of Onychostoma macrolepis against Aeromonas hydrophila infection

Galectin-9 (Gal-9) belongs to a family of the glycan-binding proteins (GBPs) and is known to restrict bacterial activity via interacting with pathogen associated molecular pattern (PAMPs). However, the underlying immune mechanism of endogenous Gal-9 on fish against bacterial infection is still unclear. In this study, effects of Gal-9 from Onychostoma macrolepis (OmGal-9) on expression of immune-related genes were measured by HEK293T. The immune response of O. macrolepis with OmGal-9 overexpression to Aeromonas hydrophila (A. hydrophila) infection (1.65 × 108 CFU/mL) was evaluated by tissue bacterial load, fish survival rate and transcriptome analysis. The results showed that OmGal-9 displayed a punctate distribution in the nucleus and cytoplasm of HEK293T cells. Compared to cells transfected with the empty vector (EV group), recombinant plasmid pEGFP-Gal9 treatment (Gal9 group) significantly down-regulated the expression of immune-related genes TNFα, STAT3, MyD88, LCK, and p52 of HEK293T cells stimulated with LPS at 24 h, while up-regulated IκBα and caspase-1 (P < 0.05). The activities of catalase (CAT), superoxide dismutase (SOD), the total antioxidant capacity (T-AOC), alkaline phosphatase (AKP), acid phosphatase (ACP), and lysozyme (LZM) of O. macrolepis were significantly increased on 7 days in Gal9 group compared to EV group (P < 0.05). The bacterial load of liver, spleen, and kidney of O. macrolepis infected with A. hydrophila in Gal9 group at 24 h was significantly lower than that in EV group (P < 0.05), and the survival rate had increased from 15 % to 35 %. A comparative transcriptome analysis between the Gal9 and EV group identified 305 differentially expressed genes (DEGs). The analysis showed that OmGal-9 might play an important regulatory role in glycolysis/gluconeogenesis, fatty acid degradation, and ascorbate and aldarate metabolism. Moreover, the immune-related DEGs were predominantly enriched in eleven pathways, with the most important three of them being linked to innate immunity: NOD-like, C-type lectin and Toll-like receptor signaling pathway. Taking together, OmGal-9 can enhance the resistance of fish to bacterial diseases by improving immune system function and activating immune-related pathways.

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.

Seasonal modulation of the testis transcriptome reveals insights into hibernation and reproductive adaptation in Onychostoma macrolepis

The Onychostoma macrolepis have a unique survival strategy, overwintering in caves and returning to the river for reproduction in summer. The current knowledge on the developmental status of its testes during winter and summer is still undiscovered. We performed RNA-seq analysis on O. macrolepis testes between January and June, using the published genome (NCBI, ASM1243209v1). Through KEGG and GO enrichment analysis, we were able to identify 2111 differentially expressed genes (DEGs) and demonstrate their functions in signaling networks associated with the development of organism. At the genomic level, we found that during the overwintering phase, genes associated with cell proliferation (ccnb1, spag5, hdac7) were downregulated while genes linked to testicular fat metabolism (slc27a2, scd, pltp) were upregulated. This indicates suppression of both mitosis and meiosis, thereby inhibiting energy expenditure through genetic regulation of testicular degeneration. Furthermore, in January, we observed the regulation of autophagy and apoptosis (becn1, casp13), which may have the function of protecting reproductive organs and ensuring their maturity for the breeding season. The results provide a basis for the development of specialized feed formulations to regulate the expression of specific genes, or editing of genes during the fish egg stage, to ensure that the testes of O. macrolepis can mature more efficiently after overwintering, thereby enhancing reproductive performance.

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.

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.

A screening for optimal selenium enrichment additives for selenium-enriched fish production: Application of a HPLC-ICP-MS method

The production of selenium-enriched fish contributes to alleviating selenium deficiency for humans. In this study, selenium nanoparticles (SeNPs) comparable in bioavailability to selenomethionine (SeMet), increased SeMet content in O. macrolepis (Onychostoma macrolepis) muscle. Additionally, dietary SeNPs significantly enhanced selenocysteine (SeCys2) and methylselenocysteine (MeSeCys) levels in O. macrolepis muscle. The effect of SeNPs on selenium speciation in grass carp muscle was consistent with O. macrolepis results. SeCys2 and MeSeCys showed antioxidant capacity in HEK293T cells, indicating enhanced health benefits of Se-enriched fish produced using SeNPs. Furthermore, the addition of 0.3 mg/kg SeNPs significantly improved the flesh quality of O. macrolepis by reducing the content of crude fat and heavy metals, as well as increasing the levels of crude protein, eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA), and the ratio of n-3/n-6 polyunsaturated fatty acids (PUFAs). Therefore, selenium-enriched fish produced from SeNPs is a good source for improving human dietary selenium intake.

Evaluating effects of climate change on the spatial distribution of an atypical cavefish Onychostoma macrolepis

Comprehending endangered species' spatial distribution in response to global climate change (GCC) is of great importance for formulating adaptive management, conservation, and restoration plans. However, it is regrettable that previous studies mainly focused on geoclimatic species, while neglected climate-sensitive subterranean taxa to a large extent, which clearly hampered the discovery of universal principles. In view of this, taking the endemic troglophile riverine fish Onychostoma macrolepis (Bleeker, 1871) as an example, we constructed a MaxEnt (maximum-entropy) model to predict how the spatial distribution of this endangered fish would respond to future climate changes (three Global Climate Models × two Shared Socio-economic Pathways × three future time nodes) based on painstakingly collected species occurrence data and a set of bioclimatic variables, including WorldClim and ENVIREM. Model results showed that variables related to temperature rather than precipitation were more important in determining the geographic distribution of this rare and endemic fish. In addition, the suitable areas and their distribution centroids of O. macrolepis would shrink (average: 20,901.75 km2) and move toward the northeast or northwest within the study area (i.e. China). Linking our results with this species' limited dispersion potential and unique habitat requirements (i.e. karst landform is essential), we thus recommended in situ conservation to protect this relict.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Selenium nanoparticles improve fish sperm quality by enhancing glucose uptake capacity via AMPK activation

Appropriate additives can provide a suitable physiological environment for storage of fish sperm and facilitate the large-scale breeding of endangered species and commercial fish. Suitable additives for fish sperm storage in vitro are required for artificial insemination. This study evaluate the effects of 0.1, 0.5, 1.5, and 4.5 mg/L selenium nanoparticles (SeNPs) on the quality of Schizothorax prenanti and Onychostoma macrolepis sperm storage in vitro at 4 °C for 72 h. We found that 0.5 mg/L SeNPs was a suitable concentration for maintaining the normal physiological state of O. macrolepis sperm during storage at 4 °C (p < 0.05). Higher adenosine triphosphate (ATP) content of O. macrolepis sperm before and after activation was present at that concentration. To further explore the potential mechanism of action of SeNPs on O. macrolepis sperm, western blotting and glucose uptake analyses were performed. The results implied that after 24 h of in vitro preservation, 0.5 mg/L SeNPs significantly improved p-AMPK levels and glucose uptake capacity of O. macrolepis sperm, while compound C (CC), the inhibitor of activated AMP-activated protein kinase (p-AMPK), significantly restricted the function of SeNPs on stored sperm. Similar effects of 0.5 mg/L SeNPs were found on Schizothorax prenanti sperm. Our study demonstrates that SeNPs maintained ATP content and O. macrolepis and Schizothorax prenanti sperm function during storage in vitro for 72 h, possibly because SeNPs enhanced the glucose uptake capacity of sperm by maintaining the level of p-AMPK.

A high-quality chromosome-level genome assembly of rohu carp, Labeo rohita, and its utilization in SNP-based exploration of gene flow and sex determination

Labeo rohita (rohu) is a carp important to aquaculture in South Asia, with a production volume close to Atlantic salmon. While genetic improvements to rohu are ongoing, the genomic methods commonly used in other aquaculture improvement programs have historically been precluded in rohu, partially due to the lack of a high-quality reference genome. Here we present a high-quality de novo genome produced using a combination of next-generation sequencing technologies, resulting in a 946 Mb genome consisting of 25 chromosomes and 2,844 unplaced scaffolds. Notably, while approximately half the size of the existing genome sequence, our genome represents 97.9% of the genome size newly estimated here using flow cytometry. Sequencing from 120 individuals was used in conjunction with this genome to predict the population structure, diversity, and divergence in three major rivers (Jamuna, Padma, and Halda), in addition to infer a likely sex determination mechism in rohu. These results demonstrate the utility of the new rohu genome in modernizing some aspects of rohu genetic improvement programs.

Comparative transcriptome provides insights into differentially expressed genes between testis and ovary of Onychostoma macrolepis in reproduction period

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered fishery species inhabiting the river of Qinling Mountains and some flowing freshwaters in China. The declining population of O. macrolepis caused by asynchrony of male and female development prompted us to focus on genetic regulation of its reproduction. In this study, high-throughput RNA-sequencing technology was applied to assemble and annotate the transcriptome of O. macrolepis testis and ovary. The results showed that a number of 338089335 (ovary:163216500, testis:174872835) raw sequences were obtained. After non-redundant analysis, a number of 207826065 (ovary:102334008, testis:105492057) high quality reads were obtained and predicted as unigenes, in which 201,038,682 unigenes were annotated with multiple databases. Taking the ovarian transcriptome as a control, comparative transcriptome analysis showed that 9918 differentially expressed genes (DEGs) up-regulated in the testis and 13,095 DEGs down-regulated. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, gonadal development, steroid biosynthesis pathways, MAPK signaling pathway and Wnt signaling pathway. Finally, the expression patterns of 19 unigenes were validated by using quantitative real-time polymerase chain reaction (qRT-PCR). This study illustrates a potential molecular mechanism on the unsynchronized male and female development of the O. macrolepis during the reproduction period in June and provides a theoretical basis for future artificial reproduction.

Equilibrated evolution of the mixed auto-/allopolyploid haplotype-resolved genome of the invasive hexaploid Prussian carp

Understanding genome evolution of polyploids requires dissection of their often highly similar subgenomes and haplotypes. Polyploid animal genome assemblies so far restricted homologous chromosomes to a 'collapsed' representation. Here, we sequenced the genome of the asexual Prussian carp, which is a close relative of the goldfish, and present a haplotype-resolved chromosome-scale assembly of a hexaploid animal. Genome-wide comparisons of the 150 chromosomes with those of two ancestral diploid cyprinids and the allotetraploid goldfish and common carp revealed the genomic structure, phylogeny and genome duplication history of its genome. It consists of 25 syntenic, homeologous chromosome groups and evolved by a recent autoploid addition to an allotetraploid ancestor. We show that de-polyploidization of the alloploid subgenomes on the individual gene level occurred in an equilibrated fashion. Analysis of the highly conserved actinopterygian gene set uncovered a subgenome dominance in duplicate gene loss of one ancestral chromosome set.

Chromosome-level genome of Tibetan naked carp (Gymnocypris przewalskii) provides insights into Tibetan highland adaptation

Gymnocypris przewalskii, a cyprinid fish endemic to the Qinghai-Tibetan Plateau, has evolved unique morphological, physiological and genetic characteristics to adapt to the highland environment. Herein, we assembled a high-quality G. przewalskii tetraploid genome with a size of 2.03 Gb and scaffold N50 of 44.93 Mb, which was anchored onto 46 chromosomes. The comparative analysis suggested that gene families related to highland adaptation were significantly expanded in G. przewalskii. According to the G. przewalskii genome, we evaluated the phylogenetic relationship of 13 schizothoracine fishes, and inferred that the demographic history of G. przewalskii was strongly associated with geographic and eco-environmental alterations. We noticed that G. przewalskii experienced whole-genome duplication, and genes preserved post duplication were functionally associated with adaptation to high salinity and alkalinity. In conclusion, a chromosome-scale G. przewalskii genome provides an important genomic resource for teleost fish, and will particularly promote our understanding of the molecular evolution and speciation of fish in the highland environment.

Testicular miRNAs and tsRNAs provide insight into gene regulation during overwintering and reproduction of Onychostoma macrolepis

The late overwintering period and breeding period are two important developmental stages of testis in Onychostoma macrolepis. Small non-coding RNAs (sncRNAs) are well-known regulators of biological processes associated with numerous biological processes. This study aimed to elucidate the roles of four sncRNA classes (microRNAs [miRNAs], Piwi-interacting RNAs [piRNAs], tRNA-derived small RNAs [tsRNAs], and rRNA-derived small RNAs [rsRNAs]) across testes in the late overwintering period (in March) and breeding period (in June) by high-throughput sequencing. The testis of O. macrolepis displayed the highest levels of piRNAs and lowest levels of rsRNAs. Compared with miRNAs and tsRNAs in June, tsRNAs in March had a higher abundance, while miRNAs in March had a much lower abundance. Bioinformatics analysis identified 1,362 and 1,340 differentially expressed miRNAs and tsRNAs, respectively. Further analysis showed that miR-200-1, miR-143-1, tRFi-Lys-CTT-1, and tRFi-Glu-CTC-1 could play critical roles during the overwintering and breeding periods. Our findings provided an unprecedented insight to reveal the epigenetic mechanism underlying the overwintering and reproduction process of male O. macrolepis.

A comprehensive overview of ovarian small non-coding RNAs in the late overwintering and breeding periods of Onychostoma macrolepis

The development of the ovary of Onychostoma macrolepis undergoes distinct annual cyclic changes in which small non-coding RNAs (sncRNAs) could play vital roles. In this study, four sncRNA classes in the ovary, including miRNA, piRNAs, tsRNA, and rsRNAs, were systematically profiled by high-throughput sequencing. In adult ovaries of O. macrolepis, 247 miRNAs and 235 tsRNAs were identified as differentially expressing in the late overwintering period (in March) and breeding period (in June). Some up-regulated sncRNAs in March, such as miR-125-1 and tRFi-Lys-CTT-1, could be involved in inhibiting biomolecule metabolism and enhancing stress tolerance during the overwintering period. Compared with the level expression of sncRNAs in March, some sncRNAs were up-regulated in June, such as miR-146-1 and tRFi-Gly-GCC-1, and could be involved in influencing molecular synthesis and metabolism, enhancing oocyte proliferation and maturation, accelerating ovarian development, and increasing fertilization of oocytes by regulating related target mRNAs. The results suggested that sncRNAs in the ovary of Onychostoma macrolepis not only reflect characteristics of the fish's physiology at different developmental periods, but also directly affect ovarian development and oocyte maturation during the breeding period. In conclusion, these results significantly advance our understanding of the roles of sncRNA during overwintering and reproduction periods, and provide a novel perspective for uncovering characteristics of the special overwintering ecology and reproductive physiology of an atypical cavefish.

MOSGA 2: Comparative genomics and validation tools

Due to the highly growing number of available genomic information, the need for accessible and easy-to-use analysis tools is increasing. To facilitate eukaryotic genome annotations, we created MOSGA. In this work, we show how MOSGA 2 is developed by including several advanced analyses for genomic data. Since the genomic data quality greatly impacts the annotation quality, we included multiple tools to validate and ensure high-quality user-submitted genome assemblies. Moreover, thanks to the integration of comparative genomics methods, users can benefit from a broader genomic view by analyzing multiple genomic data sets simultaneously. Further, we demonstrate the new functionalities of MOSGA 2 by different use-cases and practical examples. MOSGA 2 extends the already established application to the quality control of the genomic data and integrates and analyzes multiple genomes in a larger context, e.g., by phylogenetics.

High-Quality de novo Chromosome-Level Genome Assembly of a Single Bombyx mori With BmNPV Resistance by a Combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing

The reference genomes of Bombyx mori (B. mori), Silkworm Knowledge-based database (SilkDB) and SilkBase, have served as the gold standard for nearly two decades. Their use has fundamentally shaped model organisms and accelerated relevant studies on lepidoptera. However, the current reference genomes of B. mori do not accurately represent the full set of genes for any single strain. As new genome-wide sequencing technologies have emerged and the cost of high-throughput sequencing technology has fallen, it is now possible for standard laboratories to perform full-genome assembly for specific strains. Here we present a high-quality de novo chromosome-level genome assembly of a single B. mori with nuclear polyhedrosis virus (BmNPV) resistance through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. In addition, regular bioinformatics analyses, such as gene family, phylogenetic, and divergence analyses, were performed. The sample was from our unique B. mori species (NB), which has strong inborn resistance to BmNPV. Our genome assembly showed good collinearity with SilkDB and SilkBase and particular regions. To the best of our knowledge, this is the first genome assembly with BmNPV resistance, which should be a more accurate insect model for resistance studies.

Comparative genomic analysis of different sexes and diet-specific amino acid mutation identification in Ancherythroculter nigrocauda

Determining the sex and controlling the sex ratio are essential aspects of fish genetics that can assist in developing successful fish breeding programs. High quality genome assembly and annotations are prerequisites to determine sex-specific genes and their expression. In addition, analysis of resequencing data can identify genomic difference between male and female fishes. In this study, we performed chromosome-level genome assembly in female Ancherythroculter nigrocauda fish having low heterozygosity using PacBio reads. High-throughput chromatin conformation capture (HiC) yielded a genome of size 1054.05 Mb, with a contig N50 length of 3.40 Mb and a scaffold N50 length of 42.68 Mb. In addition, we sequenced 5 female and 5 male A. nigrocauda samples and identified sex-specific regions on LG20 Furthermore, diet-specific amino acid mutation were found on 582 genes between herbivorous and carnivorous fishes, with 26 of them exhibiting significantly different expression patterns in the liver tissue of these two types of fishes. The chromosome-level genome assembly of A. nigrocauda provides valuable resources for conducting in-depth comparative genomic studies with immense applications in fish genetic breeding and farming. Similarly, the diet-specific amino acid mutations are useful in the breeding of new strains of carnivorous fishes with an herbivorous diet.

Chromosome-level genome of Poropuntius huangchuchieni provides a diploid progenitor-like reference genome for the allotetraploid Cyprinus carpio

The diploid Poropuntius huangchuchieni in the cyprinid family, which is widely distributed in the Mekong and Red River basins, is one of the most closely related diploid progenitor-like species of allotetraploid common carp, which was generated by merging of two diploid genomes during evolution. Therefore, the P. huangchuchieni genome is essential for polyploid evolution studies in Cyprinidae. Here, we report a high-quality chromosome-level genome assembly of P. huangchuchieni by integrating Oxford Nanopore and Hi-C technologies. The assembled genome size was 1,021.38 Mb, 895.66 Mb of which was anchored onto 25 chromosomes with a N50 of 32.93 Mb. The genome contained 486.28 Mb repetitive elements and 24,099 protein-coding genes. Approximately 95.9% of the complete BUSCOs were detected, suggesting a high completeness of the genome. Evolutionary analysis revealed that P. huangchuchieni diverged from Cyprinus carpio at approximately 12 Mya. Genome comparison between P. huangchuchieni and the B subgenome of C. carpio provided insights into chromosomal rearrangements during the allotetraploid speciation. With the complete gene set, 17,474 orthologous genes were identified between P. huangchuchieni and C. carpio, providing a broad view of the gene component in the allotetraploid genome, which is critical for future genetic analyses. The high-quality genomic data set created for P. huangchuchieni provides a diploid progenitor-like reference for the evolution and adaptation of allotetraploid carps.