High-quality Chromosome-level Genome Research Articles

An improved chromosome-level genome assembly and annotation of Echeneis naucrates

Echeneis naucrates, as known as live sharksucker, is famous for the behavior of attaching to hosts using a highly modified dorsal fin with oval-shaped sucking disc. Here, we generated an improved high-quality chromosome-level genome assembly of E. naucrates using Illumina short reads, PacBio long reads and Hi-C data. Our assembled genome spans 572.85 Mb with a contig N50 of 23.19 Mb and is positioned to 24 pseudo-chromosomes. Additionally, at least one telomere was identified for 23 out of 24 chromosomes. Furthermore, we identified a total of 22,161 protein-coding genes, of which 21,402 genes (96.9%) were annotated successfully with functions. The combination of ab initio predictions and Repbase-based searches revealed that 15.57% of the assembled E. naucrates genome was identified as repetitive sequences. The completeness of the genome assembly and the gene annotation were estimated to be 97.5% and 95.4% with BUSCO analyses. This work enhances the utility of the live sharksucker genome and provides a valuable groundwork for the future study of genomics, biology and adaptive evolution in this species.

Chromosome-level genome assembly of the bethylid ectoparasitoid wasp Sclerodermus sp. ‘alternatusi’

The Bethylidae are the most diverse of Hymenoptera chrysidoid families. As external parasitoids, the bethylids have been widely adopted as biocontrol agents to control insect pests worldwide. Thus far, the genomic information of the family Bethylidae has not been reported yet. In this study, we crystallized into a high-quality chromosome-level genome of ant-like bethylid wasps Sclerodermus sp. ‘alternatusi’ (Hymenoptera: Bethylidae) using PacBio sequencing as well as Hi-C technology. The assembled S. alternatusi genome was 162.30 Mb in size with a contig N50 size of 3.83 Mb and scaffold N50 size of 11.10 Mb. Totally, 92.85% assembled sequences anchored to 15 pseudo-chromosomes. A total of 10,204 protein-coding genes were annotated, and 23.01 Mb repetitive sequences occupying 14.17% of genome were pinpointed. The BUSCO results showed that 97.9% of the complete core Insecta genes were identified in the genome, while 97.1% in the gene sets. The high-quality genome of S. alternatusi will not only provide valuable genomic information, but also show insights into parasitoid wasp evolution and bio-control application in future studies.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Whole-Genome Sequencing Analyses Reveal the Evolution Mechanisms of Typical Biological Features of Decapterus maruadsi.

Decapterus maruadsi is a typical representative of small pelagic fish characterized by fast growth rate, small body size, and high fecundity. It is a high-quality marine commercial fish with high nutritional value. However, the underlying genetics and genomics research focused on D. maruadsi is not comprehensive. Herein, a high-quality chromosome-level genome of a male D. maruadsi was assembled. The assembled genome length was 716.13 Mb with contig N50 of 19.70 Mb. Notably, we successfully anchored 95.73% contig sequences into 23 chromosomes with a total length of 685.54 Mb and a scaffold N50 of 30.77 Mb. A total of 22,716 protein-coding genes, 274.90 Mb repeat sequences, and 10,060 ncRNAs were predicted, among which 22,037 (97%) genes were successfully functionally annotated. The comparative genome analysis identified 459 unique, 73 expanded, and 52 contracted gene families. Moreover, 2804 genes were identified as candidates for positive selection, of which some that were related to the growth and development of bone, muscle, cardioid, and ovaries, such as some members of the TGF-β superfamily, were likely involved in the evolution of typical biological features in D. maruadsi. The study provides an accurate and complete chromosome-level reference genome for further genetic conservation, genomic-assisted breeding, and adaptive evolution research for D. maruadsi.

Chromosome-level genome provides insights into environmental adaptability and innate immunity in the common dolphin (Delphinus delphis)

The common dolphin (Delphinus delphis) is widely distributed worldwide and well adapted to various habitats. Animal genomes store clues about their pasts, and can reveal the genes underlying their evolutionary success. Here, we report the first high-quality chromosome-level genome of D. delphis. The assembled genome size was 2.56 Gb with a contig N50 of 63.85 Mb. Phylogenetically, D. delphis was close to Tursiops truncatus and T. aduncus. The genome of D. delphis exhibited 428 expanded and 1,885 contracted gene families, and 120 genes were identified as positively selected. The expansion of the HSP70 gene family suggested that D. delphis has a powerful system for buffering stress, which might be associated with its broad adaptability, longevity, and detoxification capacity. The expanded IFN-α and IFN-ω gene families, as well as the positively selected genes encoding tripartite motif-containing protein 25, peptidyl-prolyl cis-trans isomerase NIMA-interacting 1, and p38 MAP kinase, were all involved in pathways for antiviral, anti-inflammatory, and antineoplastic mechanisms. The genome data also revealed dramatic fluctuations in the effective population size during the Pleistocene. Overall, the high-quality genome assembly and annotation represent significant molecular resources for ecological and evolutionary studies of Delphinus and help support their sustainable treatment and conservation.

Chromosome-level genome assembly of Hippophae tibetana provides insights into high-altitude adaptation and flavonoid biosynthesis

BackgroundAs an endemic shrub of the Qinghai-Tibetan Plateau (QTP), the distribution of Hippophae tibetana Schlecht. ranges between 2800 and 5200 m above sea level. As the most basal branch of the Hippophae genus, H. tibetana has an extensive evolutionary history. The H. tibetana is a valuable tree for studying the ecological evolution of species under extreme conditions.ResultsHere, we generated a high-quality chromosome-level genome of H. tibetana. The total size of the assembly genome is 917 Mb. The phylogenomic analysis of 1064 single-copy genes showed a divergence between 3.4 and 12.8 Mya for H. tibetana. Multiple gene families associated with DNA repair and disease resistance were significantly expanded in H. tibetana. We also identified many genes related to DNA repair with signs of positive selection. These results showed expansion and positive selection likely play important roles in H. tibetana’s adaptation to comprehensive extreme environments in the QTP. A comprehensive genomic and transcriptomic analysis identified 49 genes involved in the flavonoid biosynthesis pathway in H. tibetana. We generated transgenic sea buckthorn hairy root producing high levels of flavonoid.ConclusionsTaken together, this H. tibetana high-quality genome provides insights into the plant adaptation mechanisms of plant under extreme environments and lay foundation for the functional genomic research and molecular breeding of H. tibetana.

Genomic, transcriptomic, and metabolomic analyses provide insights into the evolution and development of a medicinal plant Saposhnikovia divaricata (Apiaceae)

Abstract Saposhnikovia divaricata, 2n = 2x = 16, as a perennial species, is widely distributed in China, Mongolia, Russia, etc. It is a traditional Chinese herb used to treat tetanus, rubella pruritus, rheumatic arthralgia and other diseases. Here, we assembled a 2.07 Gb and N50 scaffold length of 227.67 Mb high-quality chromosome-level genome of S. divaricata based on the PacBio Sequel II sequencing platform. The total number of genes identified was 42948, and 42456 of them were functionally annotated. A total of 85.07% of the genome was composed of repeat sequences, comprised mainly of long terminal repeats (LTRs) which represented 73.7% of the genome sequence. The genome size may have been affected by a recent whole-genome duplication event. Transcriptional and metabolic analyses revealed bolting and non-bolting S. divaricata differed in flavonoids, plant hormones, and some pharmacologically active components. The analysis of its genome, transcriptome, and metabolome helped to provide insights into the evolution of bolting and non-bolting phenotypes in wild and cultivated S. divaricata and lays the basis for genetic improvement of the species.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

A high-quality chromosome-level genome assembly of the Chinese medaka Oryzias sinensis

Oryzias sinensis, also known as Chinese medaka or Chinese ricefish, is a commonly used animal model for aquatic environmental assessment in the wild as well as gene function validation or toxicology research in the lab. Here, a high-quality chromosome-level genome assembly of O. sinensis was generated using single-tube long fragment read (stLFR) reads, Nanopore long-reads, and Hi-C sequencing data. The genome is 796.58 Mb, and a total of 712.17 Mb of the assembled sequences were anchored to 23 pseudo-chromosomes. A final set of 22,461 genes were annotated, with 98.67% being functionally annotated. The Benchmarking Universal Single-Copy Orthologs (BUSCO) benchmark of genome assembly and gene annotation reached 95.1% (93.3% single-copy) and 94.6% (91.7% single-copy), respectively. Furthermore, we also use ATAC-seq to uncover chromosome transposase-accessibility as well as related genome area function enrichment for Oryzias sinensis. This study offers a new improved foundation for future genomics research in Chinese medaka.

Chromosome-Level Assembly and Comparative Genomic Analysis of Suillus bovinus Provides Insights into the Mechanism of Mycorrhizal Symbiosis.

Suillus bovinus is a wild edible ectomycorrhizal fungus with important economic and ecological value, which often forms an ectomycorrhiza with pine trees. We know little about the mechanisms associated with the metabolism and symbiosis of S. bovinus and its effects on the nutritional value. In this study, the whole-genome sequencing of S. bovinus was performed using Illumina, HiFi, and Hi-C technologies, and the sequencing data were subjected to genome assembly, gene prediction, and functional annotation to obtain a high-quality chromosome-level genome of S. bovinus. The final assembly of the S. bovinus genome includes 12 chromosomes, with a total length of 43.03 Mb, a GC content of 46.58%, and a contig N50 size of 3.78 Mb. A total of 11,199 coding protein sequences were predicted from genome annotation. The S. bovinus genome contains a large number of small secreted proteins (SSPs) and genes that encode enzymes related to carbohydrates, as well as genes related to terpenoids, auxin, and lipochitooligosaccharides. These genes may contribute to symbiotic processes. The whole-genome sequencing and genetic information provide a theoretical basis for a deeper understanding of the mechanism of the mycorrhizal symbiosis of S. bovinus and can serve as a reference for comparative genomics of ectomycorrhizal fungi.

The First High-Quality Genome Assembly of Freshwater Pearl Mussel Sinohyriopsis cumingii: New Insights into Pearl Biomineralization.

China leads the world in freshwater pearl production, an industry in which the triangle sail mussel (Sinohyriopsis cumingii) plays a pivotal role. In this paper, we report a high-quality chromosome-level genome assembly of S. cumingii with a size of 2.90 Gb-the largest yet reported among bivalves-and 89.92% anchorage onto 19 linkage groups. The assembled genome has 37,696 protein-coding genes and 50.86% repeat elements. A comparative genomic analysis revealed expansions of 752 gene families, mostly associated with biomineralization, and 237 genes under strong positive selection. Notably, the fibrillin gene family exhibited gene family expansion and positive selection simultaneously, and it also exhibited multiple high expressions after mantle implantation by transcriptome analysis. Furthermore, RNA silencing and an in vitro calcium carbonate crystallization assay highlighted the pivotal role played by one fibrillin gene in calcium carbonate deposition and aragonite transformation. This study provides a valuable genomic resource and offers new insights into the mechanism of pearl biomineralization.

A chromosome-level genome assembly of yellow stem borer (Scirpophaga incertulas)

The yellow stem borer Scirpophaga incertulas is the dominant pest of rice in tropical Asia. However, the lack of genomic resources makes it difficult to understand their invasiveness and ecological adaptation. A high-quality chromosome-level genome of S. incertulas, a monophagous rice pest, was assembled by combining Illumina short reads, PacBio HiFi long sequencing, and Hi-C scaffolding technology. The final genome size was 695.65 Mb, with a scaffold N50 of 28.02 Mb, and 93.50% of the assembled sequences were anchored to 22 chromosomes. BUSCO analysis demonstrated that this genome assembly had a high level of completeness, with 97.65% gene coverage. A total of 14,850 protein-coding genes and 366.98 Mb of transposable elements were identified. In addition, comparative genomic analyses indicated that chemosensory processes and detoxification capacity may play critical roles in the specialized host preference of S. incertulas. In summary, the chromosome-level genome assembly of S. incertulas provides a valuable genetic resource for understanding the biological characteristics of its invasiveness and developing an efficient management strategy.

High-quality chromosome-level genome assembly of the Northern Pacific sea star Asterias amurensis.

Asterias amurensis, a starfish species that is native to countries such as China and Japan, as well as non-native regions like Australia, has raised serious concerns in terms of its impact on ecology and economy. To gain a better understanding of its population genomics and dynamics, we successfully assembled a high-quality chromosome-level genome of A. amurensis using PacBio and Hi-C sequencing technologies. A total of 87 scaffolds assembly with contig N50 length of 10.85Mb and scaffold N50 length of 23.34Mb were obtained, with over 98.80% (0.48 Gb) of them anchored to 22 pseudochromosomes. We predicted 16,673 protein-coding genes, 95.19% of which were functionally annotated. Our phylogenetic analysis revealed that A. amurensis and Asterias rubens formed a clade, and their divergence time was estimated ~ 28 million years ago (Mya). The significantly enriched pathways and Gene Ontology terms related to the amplified gene family were mainly associated with immune response and energy metabolism, suggesting that these factors might have contributed to the adaptability of A. amurensis to its environment. This study provides valuable genomic resources for comprehending the genetics, dynamics, and evolution of A. amurensis, especially when population outbreaks or invasions occur.

Chromosome-level genome assembly of the silver pomfret Pampus argenteus

Pampus argenteus (Euphrasen, 1788) is one of the major fishery species in coastal China. Pampus argenteus has a highly specialized morphology, and its declining fishery resources have encouraged massive research efforts on its aquacultural biology. In this study, we reported the first high-quality chromosome-level genome of P. argenteus obtained by integrating Illumina, PacBio HiFi, and Hi-C sequencing techniques. The final size of the genome was 518.06 Mb, with contig and scaffold N50 values of 20.47 and 22.86 Mb, respectively. The sequences were anchored and oriented onto 24 pseudochromosomes based on Hi-C data corresponding to the 24-chromatid karyotype of P. argenteus. A colinear relationship was observed between the P. argenteus genome and that of a closely related species (Scomber japonicus). A total of 24,696 protein-coding genes were identified from the genome, 98.9% of which were complete BUSCOs. This report represents the first case of high-quality chromosome-level genome assembly for P. argenteus and can provide valuable information for future evolutionary, conservation, and aquacultural research.

High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation.

High-quality genome of rosemary (Salvia rosmarinus) represents a valuable resource and tool for understanding genome evolution and environmental adaptation as well as its genetic improvement. However, the existing rosemary genome did not provide insights into the relationship between antioxidant components and environmental adaptability. In this study, by employing Nanopore sequencing and Hi-C technologies, a total of 1.17 Gb (97.96%) genome sequences were mapped to 12 chromosomes with 46 121 protein-coding genes and 1265 non-coding RNA genes. Comparative genome analysis reveals that rosemary had a closely genetic relationship with Salvia splendens and Salvia miltiorrhiza, and it diverged from them approximately 33.7 million years ago (MYA), and one whole-genome duplication occurred around 28.3 MYA in rosemary genome. Among all identified rosemary genes, 1918 gene families were expanded, 35 of which are involved in the biosynthesis of antioxidant components. These expanded gene families enhance the ability of rosemary adaptation to adverse environments. Multi-omics (integrated transcriptome and metabolome) analysis showed the tissue-specific distribution of antioxidant components related to environmental adaptation. During the drought, heat and salt stress treatments, 36 genes in the biosynthesis pathways of carnosic acid, rosmarinic acid and flavonoids were up-regulated, illustrating the important role of these antioxidant components in responding to abiotic stresses by adjusting ROS homeostasis. Moreover, cooperating with the photosynthesis, substance and energy metabolism, protein and ion balance, the collaborative system maintained cell stability and improved the ability of rosemary against harsh environment. This study provides a genomic data platform for gene discovery and precision breeding in rosemary. Our results also provide new insights into the adaptive evolution of rosemary and the contribution of antioxidant components in resistance to harsh environments.

Chromosome-level genome assembly and annotation of the yellow grouper, Epinephelus awoara

Epinephelus awoara, as known as yellow grouper, is a significant economic marine fish that has been bred artificially in China. However, the genetic structure and evolutionary history of yellow grouper remains largely unknown. Here, this work presents the high-quality chromosome-level genome assembly of yellow grouper using PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The 984.48 Mb chromosome-level genome of yellow grouper was assembled, with a contig N50 length of 39.77 Mb and scaffold N50 length of 41.39 Mb. Approximately 99.76% of assembled sequences were anchored into 24 pseudo-chromosomes with the assistance of Hi-C reads. Furthermore, approximately 41.17% of the genome was composed of repetitive elements. In total, 24,541 protein-coding genes were predicted, of which 22,509 (91.72%) genes were functionally annotated. The highly accurate, chromosome-level reference genome assembly and annotation are crucial to the understanding of population genetic structure, adaptive evolution and speciation of the yellow grouper.

A chromosome-level genome assembly of the Echiura Urechis unicinctus

Echiura is a distinctive family of unsegmented sausage-shaped marine worms whose phylogenetic relationship still needs strong evidence from the phylogenomic analysis. In this family, Urechis unicinctus is known for its high nutritional and medicinal value and adaptation to harsh intertidal conditions. Herein, we combined PacBio long-read, short-read Illumina and Hi-C sequencing, generating a high-quality chromosome-level genome assembly of U. unicinctus. The assembled genome spans ~1,138.6 Mb with a scaffold N50 of 68.3 Mb, of which 1,113.8 Mb (97.82%) were anchored into 17 pseudo-chromosomes. The BUSCO analysis demonstrated the completeness of the genome assembly and gene model prediction are 93.5% and 91.5%, respectively. A total of 482.1 Mb repetitive sequences, 21,524 protein-coding genes, 1,535 miRNAs, 3,431 tRNAs, 124 rRNAs, and 348 snRNAs were annotated. This study significantly improves the quality of U. unicinctus genome assembly, sets the footsteps for molecular breeding and further study in genome evolution, genetic and molecular biology of U. unicinctus.

A chromosome-level genome assembly of a free-living white-crowned sparrow (Zonotrichia leucophrys gambelii)

The white-crowned sparrow, Zonotrichia leucophrys, is a passerine bird with a wide distribution and it is extensively adapted to environmental changes. It has historically acted as a model species in studies on avian ecology, physiology and behaviour. Here, we present a high-quality chromosome-level genome of Zonotrichia leucophrys using PacBio and OmniC sequencing data. Gene models were constructed by combining RNA-seq and Iso-seq data from liver, hypothalamus, and ovary. In total a 1,123,996,003 bp genome was generated, including 31 chromosomes assembled in complete scaffolds along with other, unplaced scaffolds. This high-quality genome assembly offers an important genomic resource for the research community using the white-crowned sparrow as a model for understanding avian genome biology and development, and provides a genomic basis for future studies, both fundamental and applied.

A chromosome-level genome assembly for Onobrychis viciifolia reveals gene copy number gain underlying enhanced proanthocyanidin biosynthesis

Sainfoin (Onobrychis viciifolia), which belongs to subfamily Papilionoideae of Leguminosae, is a vital perennial forage known as “holy hay” due to its high contents of crude proteins and proanthocyanidins (PAs, also called condensed tannins) that have various pharmacological properties in animal feed, such as alleviating rumen tympanic disease in ruminants. In this study, we select an autotetraploid common sainfoin (2n = 4x = 28) and report its high-quality chromosome-level genome assembly with 28 pseudochromosomes and four haplotypes (~1950.14 Mb, contig N50 = 10.91 Mb). The copy numbers of genes involved in PA biosynthesis in sainfoin are significantly greater than those in four selected Fabales species, namely, autotetraploid Medicago sativa and three other diploid species, Lotus japonicus, Medicago truncatula, and Glycine max. Furthermore, gene expansion is confirmed to be the key contributor to the increased expression of these genes and subsequent PA enhancement in sainfoin. Transcriptomic analyses reveal that the expression of genes involved in the PA biosynthesis pathway is significantly increased in the lines with high PA content compared to the lines with medium and low PA content. The sainfoin genome assembly will improve our understanding of leguminous genome evolution and biosynthesis of secondary metabolites in sainfoin.

Improved chromosome-level genome assembly of Indian sandalwood (Santalum album)

Santalum album is a well-known aromatic and medicinal plant that is highly valued for the essential oil (EO) extracted from its heartwood. In this study, we present a high-quality chromosome-level genome assembly of S. album after integrating PacBio Sequel, Illumina HiSeq paired-end and high-throughput chromosome conformation capture sequencing technologies. The assembled genome size is 207.39 M with a contig N50 of 7.33 M and scaffold N50 size of 18.31 M. Compared with three previously published sandalwood genomes, the N50 length of the genome assembly was longer. In total, 94.26% of the assembly was assigned to 10 pseudo-chromosomes, and the anchor rate far exceeded that of a recently released value. BUSCO analysis yielded a completeness score of 94.91%. In addition, we predicted 23,283 protein-coding genes, 89.68% of which were functionally annotated. This high-quality genome will provide a foundation for sandalwood functional genomics studies, and also for elucidating the genetic basis of EO biosynthesis in S. album.