Chromosome-level Genome Assembly Research Articles

Chromosome genome assembly and annotation of Adzuki Bean (Vigna angularis)

Adzuki bean (Vigna angularis) is a significant dietary legume crop that is prevalent in East Asia. It also holds traditional medicinal importance in China. In this study, we report a high-quality, chromosome-level genome assembly of adzuki bean obtained by employing Illumina short-read sequencing, PacBio long-read sequencing, and Hi-C technology. The assembly spans 447.8 Mb, encompassing 96.32% of the estimated genome, with contig and scaffold N50 values of 16.5 and 41.0 Mb, respectively. More than 98.2% of the 1,614 BUSCO genes were fully identified, and 25,939 genes were annotated, with 98.23% of them being functionally identifiable. Vigna angularis was estimated to diverge successively from Vigna unguiculata and Vigna radiata about 15.3 and 8.7 million years ago (Ma), respectively. This chromosome-level reference genome of Vigna angularis provides a robust foundation for exploring the functional genomics and genome evolution of adzuki bean, thereby facilitating advancements in molecular breeding of adzuki bean.

A chromosome level reference genome of Diviner’s sage (Salvia divinorum) provides insight into salvinorin A biosynthesis

BackgroundDiviner’s sage (Salvia divinorum; Lamiaceae) is the source of the powerful hallucinogen salvinorin A (SalA). This neoclerodane diterpenoid is an agonist of the human Κ-opioid receptor with potential medical applications in the treatment of chronic pain, addiction, and post-traumatic stress disorder. Only two steps of the approximately twelve step biosynthetic sequence leading to SalA have been resolved to date.ResultsTo facilitate pathway elucidation in this ethnomedicinal plant species, here we report a chromosome level genome assembly. A high-quality genome sequence was assembled with an N50 value of 41.4 Mb and a BUSCO completeness score of 98.4%. The diploid (2n = 22) genome of ~ 541 Mb is comparable in size and ploidy to most other members of this genus. Two diterpene biosynthetic gene clusters were identified and are highly enriched in previously unidentified cytochrome P450s as well as crotonolide G synthase, which forms the dihydrofuran ring early in the SalA pathway. Coding sequences for other enzyme classes with likely involvement in downstream steps of the SalA pathway (BAHD acyl transferases, alcohol dehydrogenases, and O-methyl transferases) were scattered throughout the genome with no clear indication of clustering. Differential gene expression analysis suggests that most of these genes are not inducible by methyl jasmonate treatment.ConclusionsThis genome sequence and associated gene annotation are among the highest resolution in Salvia, a genus well known for the medicinal properties of its members. Here we have identified the cohort of genes responsible for the remaining steps in the SalA pathway. This genome sequence and associated candidate genes will facilitate the elucidation of SalA biosynthesis and enable an exploration of its full clinical potential.

Chromosome-level genome assembly and population genomic analysis provide insights into the genetic diversity and adaption of Schizopygopsis younghusbandi on the Tibetan Plateau.

The Yarlung Tsangpo River on the Tibetan Plateau provides a unique natural environment for studying fish evolution and ecology. However, the genomes and genetic diversity of plateau fish species have been rarely reported. Schizopygopsis younghusbandi, a highly specialized Schizothoracine species and economically important fish inhabiting the Yarlung Tsangpo River, is threatened by overfishing and biological invasion. Herein, we generated a chromosome-level genome of S. younghusbandi and whole-genome resequencing data for 59 individuals from six locations of the river. The results showed that the divergence time between S. younghusbandi and other primitive Schizothoracine species was ∼4.2 Mya, coinciding with the major phase of the Neogene Tibetan uplift. The expanded gene families enriched in DNA integration and replication, ion binding and transport, energy storage, and metabolism likely contribute to the adaption of this species. The S. younghusbandi may have diverged from other highly specialized Schizothoracine species in the Zanda basin during the Pliocene epoch, which underwent major population reduction possibly due to the drastic climate change during the last glacial period. Population analysis indicated that the ancient population might have originated upstream before gradually adapting to evolve into the populations inhabiting the mid-stream and downstream regions of the Yarlung Tsangpo River. In conclusion, the chromosome-level genome and population diversity of S. younghusbandi provide valuable genetic resources for the evolution, ecology, and conservation studies of endemic fishes on the Tibetan Plateau.

An integrative strategy used by the aphid Uroleucon formosanum to counter host sesquiterpene lactone defense: Insights from combined genomic and transcriptomic analysis.

Insect herbivores adapt and develop strategies to counteract plant chemical defenses. The aphid Uroleucon formosanum is a serious sap-sucking pest that infests lettuces containing toxic sesquiterpene lactones (STLs). Herein, we employed a combination of genome sequencing and RNA-seq transcriptome profiling to understand the mechanisms underlying phytotoxin tolerance in U.formosanum. We generated the first chromosome-level genome assembly for U.formosanum, with a total size of 453.26Mb and a scaffold N50 of 33.22Mb. Comparative genomic analyses revealed an enrichment of signals for positive selection and gene family expansion in immune-related pathways. Specifically, the expanded set of heat shock protein70 (HSP70) genes showed upregulation after treatment with lactucin, suggesting that they may play a role in the immune response against STLs. The expression of takeout-like genes and cuticle-associated genes was also significantly increased in the lactucin-treated samples. Additionally, 53 cytochrome P450 monooxygenase, 30 carboxylesterase, 19 glutathione S-transferase, 32 uridine diphosphate glycosyltransferase and 63 ATP-binding cassette (ABC) transporter genes were identified in the U.formosanum genome. CYP4C1, CYP6A13 and 7 ABC genes were strongly upregulated in response to lactucin treatment, indicating the involvement of detoxifying enzymes in the tolerance of U.formosanum to STLs. Our findings suggest that the cuticle barrier, immune response and enzyme-mediated metabolic detoxification jointly enhance the tolerance of U.formosanum to phytotoxins and promote its adaptation to host plants. This study presents a valuable genomic resource and provides insights into insect adaptation to plant chemical challenges and future technological developments for pest management.

High-quality chromosome-level genome assembly of female Artemia franciscana reveals sex chromosome and Hox gene organization

Artemia is a crustacean genus belonging to the order Anostraca in the class Branchiopoda and lives in inland hypersaline lakes. Among the genus, A. franciscana is a valuable species as a fish food in the aquaculture industry or as an aquatic model organism for toxicity tests. However, genomic data for A. franciscana remains incomplete. In this study, high-quality genome assembly at the chromosome level of female A. franciscana was conducted by combining various sequencing and assembly technologies. The final A. franciscana assembled genome was 1.27 Gb in length, containing 21 chromosomal scaffolds (>10 Mb). The scaffold N50 was 45.3 Mb, with a complete BUSCO value of 91.0 %, thereby confirming that a high-quality genome was assembled. Gene annotation shows that the A. franciscana genome contained 67.26 % of repetitive sequences, and a total of 26,923 protein-coding genes were predicted. Among the 21 chromosome-scale scaffolds, chromosome 1 was identified as a sex chromosome Z. Additionally, five contigs of putative W chromosome fragments and the candidate sex-determining genes were suggested. Ten homeobox (Hox) genes were identified in A. franciscana on the chromosome 14, which were in two subclusters with a large gap. Hox gene organizations within 13 arthropods showed that four anostracans had conserved synteny. This study provides a new female Artemia genome with sex chromosome and the first complete genomic arrangement of the Hox cluster in Anostraca. This study will be a useful genomic and genetic reference for understanding the evolution and development of A. franciscana.

A new high-quality genome assembly and annotation for the threatened Florida Scrub-Jay (Aphelocoma coerulescens).

The Florida Scrub-Jay (Aphelocoma coerulescens), a Federally Threatened, cooperatively-breeding bird, is an emerging model system in evolutionary biology and ecology. Extensive individual-based monitoring and genetic sampling for decades has yielded a wealth of data, allowing for the detailed study of social behavior, demography, and population genetics of this natural population. Here, we report a linkage map and a chromosome-level genome assembly and annotation for a female Florida Scrub-Jay made with long-read sequencing technology, chromatin conformation data, and the linkage map. We constructed a linkage map comprising 4,468 SNPs that had 34 linkage groups and a total sex-averaged autosomal genetic map length of 2446.78 cM. The new genome assembly is 1.33 Gb in length, consisting of 33 complete or near-complete autosomes and the sex chromosomes (ZW). This highly contiguous assembly has an NG50 of 68 Mb and a Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score of 97.1% with respect to the Aves database. The annotated gene set has a BUSCO transcriptome completeness score of 95.5% and 17,964 identified protein-coding genes, 92.5% of which have associated functional annotations. This new, high-quality genome assembly and linkage map of the Florida Scrub-Jay provides valuable tools for future research into the evolutionary dynamics of small, natural populations of conservation concern.

High-quality Chromosomal-Level Genome Assembly of the Wasabi (Eutrema japonicum) ‘Magic’

Wasabi (Eutrema japonicum) is a plant belonging to the Brassicaceae family that produces its distinctive pungent taste through allyl isothiocyanate. This study achieved a high-quality chromosome-level genome assembly of the E. japonicum ‘Magic’ bred in Korea for its rapid growth cycle. The assembly was accomplished using a combination of Illumina, PacBio HIFI, Nanopore MinION, and Pore-C scaffolding technologies. The final assembled genome size is 794.6 Mb, anchored to 14 chromosomes. The genome comprises 67.56% repetitive elements and has a BUSCO score of 99.3%, indicating a high level of completeness. Compared to previously published assemblies with a different cultivar, the total length increased by approximately 48.08 Mb, while the number of Ns decreased from 89,000 to 49,000, and the assembly gaps (500 N padding) reduced from 178 to 98, resulting in a higher quality assembly. This genome will be a valuable resource for genetic and biological research on E. japonicum, aiding in its breeding and genetic improvement.

Limited population structure but signals of recent selection in introduced African Fig Fly (Zaprionus indianus) in North America.

Invasive species have devastating consequences for human health, food security, and the environment. Many invasive species adapt to new ecological niches following invasion, but little is known about the early steps of adaptation. Here we examine population genomics of a recently introduced drosophilid in North America, the African Fig Fly, Zaprionus indianus. This species is likely intolerant of subfreezing temperatures and recolonizes temperate environments yearly. We generated a new chromosome-level genome assembly for Z. indianus. Using resequencing of over 200 North American individuals collected over four years in temperate Virginia, plus a single collection from subtropical Florida, we tested for signatures of recolonization, population structure, and adaptation within invasive populations. We show founding populations are sometimes small and contain close genetic relatives, yet temporal population structure and differentiation of populations is mostly absent across recurrent recolonization events. Although we find limited signals of genome-wide spatial or temporal population structure, we identify haplotypes on the X chromosome that are repeatedly differentiated between Virginia and Florida populations. These haplotypes show signatures of natural selection and are not found in African populations. We also find evidence for several large structural polymorphisms segregating within North America populations and show X chromosome evolution in invasive populations is strikingly different from the autosomes. These results show that despite limited population structure, populations may rapidly evolve genetic differences early in an invasion. Further uncovering how these genomic regions influence invasive potential and success in new environments will advance our understanding of how organisms evolve in changing environments.

A deep learning-based method enables the automatic and accurate assembly of chromosome-level genomes.

The application of high-throughput chromosome conformation capture (Hi-C) technology enables the construction of chromosome-level assemblies. However, the correction of errors and the anchoring of sequences to chromosomes in the assembly remain significant challenges. In this study, we developed a deep learning-based method, AutoHiC, to address the challenges in chromosome-level genome assembly by enhancing contiguity and accuracy. Conventional Hi-C-aided scaffolding often requires manual refinement, but AutoHiC instead utilizes Hi-C data for automated workflows and iterative error correction. When trained on data from 300+species, AutoHiC demonstrated a robust average error detection accuracy exceeding 90%. The benchmarking results confirmed its significant impact on genome contiguity and error correction. The innovative approach and comprehensive results of AutoHiC constitute a breakthrough in automated error detection, promising more accurate genome assemblies for advancing genomics research.

High-quality haplotype-resolved chromosome assembly provides evolutionary insights and targeted steviol glycosides (SGs) biosynthesis in Stevia rebaudiana Bertoni.

Stevia rebaudiana Bertoni is popular source of plant-derived low/no-calorie natural sweeteners (LNCSs), collectively known as steviol glycosides (SGs). Nevertheless, genetic predisposition for targeted biosynthesis of SGs is complex due to multi-substrate functionality of key uridine diphosphate glycosyltransferases (UGTs). Here, we created a high-quality monoploid assembly of 1.34 Gb with N50 value of 110 Mb, 55 551 predicted protein-coding genes, and ~80% repetitive regions in Rebaudioside-A (Reb-A) enriched cultivar of S. rebaudiana. Additionally, a haplotype-based chromosome assembly consisting of haplotype A and haplotype B with an overall genome size of 2.33Gb was resolved, harbouring 639 634 variants including single nucleotide polymorphisms (SNPs), indels and structural variants (SVs). Furthermore, a lineage-specific whole genome duplication analysis revealed that gene families encoding UGTs and Cytochrome-P450 (CYPs) were tandemly duplicated. Additionally, expression analysis revealed five tandemly duplicated gene copies of UGT76G1 having significant correlations with Reb-A content, and identified key residue (leu200val) in the glycosylation of Reb-A. Furthermore, missense variations identified in the acceptor region of UGT76G1 in haplotype resolve genome, transcriptional and molecular docking analysis were confirmed with resequencing of 10 diverse stevia genotypes (~25X). Gene regulatory network analysis identified key transcription factors (MYB, bHLH, bZIP and AP2-ERF) as potential regulators of SG biosynthesis. Overall, this study provides haplotype-resolved chromosome-level genome assembly for genome editing and enhancing breeding efforts for targeted biosynthesis of SGs in S. rebaudiana.

168 American mink Genome Sequencing Project

Abstract The long-term goal of American Mink (Neogale vison) Genome Sequencing Project is to develop novel genomics tools to study the genome biology and evolution of mink. This research agenda is organized around three major goals: 1) development of first chromosome-level genome assembly with high quality for American mink to help understanding the biology and evolution of the order Carnivora, 2) designing the first mink single nucleotide polymorphism (SNP) genotyping array for genomic discovery in mink, and 3) using genomics to understand the genetic architecture of economically important traits in mink. An initial draft of the genome assembly is generated using 2,884,047 PacBio long reads. Integration of Hi-C data into the initial draft led to an assembly with 183 scaffolds and scaffold N50 of 220 Mb. This gap-free genome assembly of American mink (ASM_NN_V1) has a length of 2.68 Gb in which about 98.6% of the whole genome is covered by 15 chromosomes. Significant improvements are observed in genome contiguity, the number of scaffolds, and annotation compared with the first draft of mink genome assembly. A custom genotyping array is developed using the developed mink reference genome and DNA panel composed of 100 whole-genome sequences data of American mink from five color-types. This Affymetrix Mink 70K SNP array created 70,000 SNPs, which is sufficient for conducting genomic studies in the mink research program and beyond. This project would facilitate the discovery of genomic regions, markers and genetic networks underlying economically important traits in American mink, which will contribute novel insights into the genetic architecture of key traits and is a critical step to implement genomic selection in mink breeding programs. The results of this research have the potential to impact the mink industry, especially in North America, and improve the overall performance of the mink industry.

An expanded odorant-binding protein mediates host cue detection in the parasitic wasp Baryscapus dioryctriae basis of the chromosome-level genome assembly analysis

BackgroundBaryscapus dioryctriae (Chalcidodea: Eulophidae) is a parasitic wasp that parasitizes the pupae of many Pyralidae members and has been used as a biological control agent against Dioryctria pests of pinecones.ResultsThis B. dioryctriae assembly has a genome size of 485.5 Mb with a contig N50 of 2.17 Mb, and scaffolds were assembled onto six chromosomes using Hi-C analysis, significantly increasing the scaffold N50 to 91.17 Mb, with more than 96.13% of the assembled bases located on chromosomes, and an analysis revealed that 94.73% of the BUSCO gene set. A total of 54.82% (279.27 Mb) of the assembly was composed of repetitive sequences and 24,778 protein-coding genes were identified. Comparative genomic analysis demonstrated that the chemosensory perception, genetic material synthesis, and immune response pathways were primarily enriched in the expanded genes. Moreover, the functional characteristics of an odorant-binding protein (BdioOBP45) with ovipositor-biased expression identified from the expanded olfactory gene families were investigated by the fluorescence competitive binding and RNAi assays, revealing that BdioOBP45 primarily binds to the D. abietella-induced volatile compounds, suggesting that this expanded OBP is likely involved in locating female wasp hosts and highlighting a direction for future research.ConclusionsTaken together, this work not only provides new genomic sequences for the Hymenoptera systematics, but also the high-quality chromosome-level genome of B. dioryctriae offers a valuable foundation for studying the molecular, evolutionary, and parasitic processes of parasitic wasps.

Genomic changes and stabilization following homoploid hybrid speciation of the Oxford ragwort Senecio squalidus

Oxford ragwort (Senecio squalidus) is one of only two homoploid hybrid species known to have originated very recently, so it is a unique model for determining genomic changes and stabilization following homoploid hybrid speciation. Here, we provide a chromosome-level genome assembly of S.squalidus with 95% of the assembly contained in the 10 longest scaffolds, corresponding to its haploid chromosome number. We annotated30,249 protein-coding genes and estimated that ∼62% of the genome consists of repetitive elements. We then characterized genome-wide patterns of linkage disequilibrium, polymorphism, and divergence in S.squalidus and its two parental species, finding that (1) linkage disequilibrium is highly heterogeneous, with a region on chromosome 4 showing increased values across all three species but especially in S.squalidus; (2) regions harboring genetic incompatibilities between the two parental species tend to be large, show reduced recombination, and have lower polymorphism in S.squalidus; (3) the two parental species have an unequal contribution (70:30) to the genome of S.squalidus, with long blocks of parent-specific ancestry supporting a very rapid stabilization of the hybrid lineage after hybrid formation; and (4) genomic regions with major parent ancestry exhibit an overrepresentation of loci with evidence for divergent selection occurring between the two parental species on Mount Etna. Our results show that both genetic incompatibilities and natural selection play a role in determining genome-wide reorganization following hybrid speciation and that patterns associated with homoploid hybrid speciation-typically seen in much older systems-can evolve very quickly following hybridization.

A chromosome-level genome assembly of Drosophila madeirensis, a fruit fly species endemic to the island of Madeira.

Drosophila subobscura is distributed across Europe, the Near East, and the Americas, while its sister species, Drosophila madeirensis, is endemic to the island of Madeira in the Atlantic Ocean. D. subobscura is known for its strict light-dependence in mating and its unique courtship displays, including nuptial gift-giving. D. subobscura has also attracted the interest of researchers because of its abundant variations in chromosomal polymorphisms correlated to the latitude and season, which have been used as a tool to track global climate warming. Although D. madeirensis can be an important resource for understanding the evolutionary underpinning of these genetic characteristics of D. subobscura, little work has been done on the biology of this species. Here, we used a HiFi long-read sequencing data set to produce a de novo genome assembly for D. madeirensis. This assembly comprises a total of 111 contigs spanning 135.5 Mb and has an N50 of 24.2 Mb and a BUSCO completeness score of 98.6%. Each of the 6 chromosomes of D. madeirensis consisted of a single contig except for some centromeric regions. Breakpoints of the chromosomal inversions between D. subobscura and D. madeirensis were characterized using this genome assembly, updating some of the previously identified locations.

The first Chromosomal-level genome assembly of Sageretia thea using Nanopore long reads and Pore-C technology

Sageretia thea, a notable species within the mock buckthorn genus, is recognized for its intriguing biogeographical distribution and diverse medicinal properties. Despite this significance, genomic studies on S. thea are still in the nascent stages. We present the first chromosome-level genome assembly of S. thea that was generated using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies complemented by Pore-C chromatin conformation capture. The genome assembly had a size of 197.8 Mb with 12 chromosomal scaffolds and a scaffold N50 length of 15.9 Mb. A total of 25,434 protein-coding genes were identified and functionally annotated, and the gene model indicated 96.5% complete eukaryotic BUSCOs. Additionally, orthologous gene profiling and synteny analysis were performed to elucidate the evolutionary relationships within the Rhamnaceae family and Rosales. This high-quality chromosomal genome is the first genomic view of S. thea, which will serve as the basis for future studies on its biological and medicinal properties, and evolutionary history.

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila.

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Pan-genome analyses of 11 Fraxinus species provide insights into salt adaptation in ash trees

Ash trees (Fraxinus) exhibit rich genetic diversity and wide adaptation to various ecological environments, and several species are highly salt tolerant. Dissecting the genomic basis of salt adaptation in Fraxinus is vital for its resistance breeding. Here, we present 11 high-quality chromosome-level genome assemblies for Fraxinus species, which reveal two unequal subgenome compositions and two recent whole-genome triplication events in their evolutionary history. A Fraxinus pan-genome was constructed on the basis of structural variations and revealed that presence–absence variations (PAVs) of transmembrane transport genes have likely contributed to salt adaptation in Fraxinus. Through whole-genome resequencing of an F1 population from an interspecies cross of F. velutina ‘Lula 3’ (salt tolerant) with F. pennsylvanica ‘Lula 5’ (salt sensitive), we mapped salt-tolerance PAV-based quantitative trait loci (QTLs) and pinpointed two PAV-QTLs and candidate genes associated with Fraxinus salt tolerance. Mechanistically, FvbHLH85 enhances salt tolerance by mediating reactive oxygen species and Na+/K+ homeostasis, whereas FvSWEET5 enhances salt tolerance by mediating osmotic homeostasis. Collectively, these findings provide valuable genomic resources for Fraxinus salt-resistance breeding and the research community.

The chromosome-level genome assembly of Cananga odorata provides insights into its evolution and terpenoid biosynthesis.

Cananga odorata is known as a natural perfume tree of the Annonaceae family in Magnoliales. However, its phylogenetic position and the molecular mechanisms involved in the biosynthesis of the floral volatile organic compounds (VOCs) remain unclear. Here, by combining a variety of sequencing platforms, we present a telomere-to-telomere (T2T) genome of C. odorata with 735.83 Mb, which represents the highest integrity and assembly quality of genome in magnoliid plants reported to date. Phylogenetic analysis based on multiple datasets and approaches showed that C. odorata, as a member of magnoliids, is sister to eudicots, after their divergence from monocots. Metabolomic of VOCs in the essential oil and flowers scent showed that sesquiterpenes, especially β-caryophyllene, were the major compounds. Two CoTPS21 homologues derived from tandem duplication events were highly expressed during flower development and were identified as the key sesquiterpene synthases for the production of β-caryophyllene. In addition, CoSPL3 and CoSPL9 were considered as potential transcription factors for activating the expression of CoTPS21 homologues. Our results shed light on the molecular mechanisms underlying the biosynthesis of the unique floral fragrance in C. odorata and provide new insights into the phylogenetic position of magnoliids.

Chromosome-level genome assemblies of vulnerable male and female elongate loach (Leptobotia elongata)

Endemic to the upper and middle reaches of the Yangtze River in China, elongate loach (Leptobotia elongata) has become a vulnerable species mainly due to overfishing and habitat destruction. Thus far, no genome data of this species are reported. As a result, lacking of such genomic information has restricted practical conservation and utilization of this economic fish. Here, we constructed chromosome-level genome assemblies for both male and female elongate loach by integration of MGI, PacBio HiFi and Hi-C sequencing technologies. Two primary genome assemblies (586-Mb and 589-Mb) were obtained for female and male fishes, respectively. Indeed, 98.22% and 98.61% of the contig sequences were anchored onto 25 chromosomes, with identification of 26.22% and 25.92% repeat contents in both assembled genomes. Meanwhile, a total of 25,215 and 25,253 protein-coding genes were annotated, of which 97.41% and 98.8% could be predicted with functions. Taken together, our genome data presented here provide a valuable genomic resource for in-depth evolutionary and functional research, as well as molecular breeding and conservation of this economic fish species.

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.