Chromosome-scale Genome Assembly Research Articles

Chromosome-Scale Genome of the Fern Cibotium barometz Unveils a Genetic Resource of Medicinal Value

Ferns represent the second-largest group of vascular plants, yet their genomic resources lag far behind. Here, we present a chromosome-scale genome assembly of Cibotium barometz (L.) J. Sm., a medicinally important fern species. The 3.49 Gb genome, assembled into 66 chromosomes with 99.41% sequence anchorage, revealed an exceptionally high proportion (83.93%) of repetitive elements, dominated by recently expanded LTR retrotransposons. We identified 30,616 protein-coding genes, providing insights into fern-specific gene families. Genomic analyses uncover the evolutionary dynamics of 513 key biosynthetic genes, particularly those involved in terpenoid and flavonoid production. Expression profiling across tissues revealed tissue-specific regulation of these pathways, with notable upregulation of chalcone synthase genes in roots. Our structural analysis of 1-deoxy-d-xylulose-5-phosphate synthase, a key enzyme in terpenoid biosynthesis, demonstrated high conservation across land plants while highlighting fern-specific adaptations. The identification of multiple isoforms for key enzymes points to potential gene-duplication events or the evolution of fern-specific variants. This genome provides a foundation for understanding fern biology, evolution, and the molecular basis of their medicinal properties. It also offers valuable resources for conservation efforts and pharmacological research, paving the way for sustainable utilization of this valuable medicinal plant and advancing our understanding of plant diversity and natural product biosynthesis.

Repeat-based holocentromeres of the woodrush Luzula sylvatica reveal insights into the evolutionary transition to holocentricity

In most studied eukaryotes, chromosomes are monocentric, with centromere activity confined to a single region. However, the rush family (Juncaceae) includes species with both monocentric (Juncus) and holocentric (Luzula) chromosomes, where centromere activity is distributed along the entire chromosome length. Here, we combine chromosome-scale genome assembly, epigenetic analysis, immuno-FISH and super-resolution microscopy to study the transition to holocentricity in Luzula sylvatica. We report repeat-based holocentromeres with an irregular distribution of features along the chromosomes. Luzula sylvatica holocentromeres are predominantly associated with two satellite DNA repeats (Lusy1 and Lusy2), while CENH3 also binds satellite-free gene-poor regions. Comparative repeat analysis suggests that Lusy1 plays a crucial role in centromere function across most Luzula species. Furthermore, synteny analysis between L. sylvatica (n = 6) and Juncus effusus (n = 21) suggests that holocentric chromosomes in Luzula could have arisen from chromosome fusions of ancestral monocentric chromosomes, accompanied by the expansion of CENH3-associated satellite repeats.

Chromosome-scale genome assembly of Codonopsis pilosula and comparative genomic analyses shed light on its genome evolution

IntroductionCodonopsis pilosula is a significant plant in traditional Chinese medicine, valued for its edible and medicinal properties. However, the lack of available genomic resources has hindered further research.MethodsThis study presents the first chromosome-scale genome assembly of C. pilosula using PacBio CLR reads and Hi-C scaffolding technology. Additionally, Ks analysis and syntenic depth analysis were performed to elucidate its evolutionary history.ResultsThe final assembly yielded a high-quality genome of 679.20 Mb, which was anchored to 8 pseudo-chromosomes with an anchoring rate of 96.5% and a scaffold N50 of 80.50 Mb. The genome assembly showed a high completeness of 97.6% based on Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis. Repetitive elements constituted approximately 76.8% of the genome, with long terminal repeat retrotransposons (LTRs) accounting for about 39.17%. Ks and syntenic depth analyses revealed that the polyploidization history of three platycodonoid clade species involved only the γ-WGT event. Karyotype evolutionary analysis identified an ancestral karyotype with 9 protochromosomes for the three platycodonoid clade species. Moreover, non-WGD genes, particularly those arising from tandem duplications, were found to contribute significantly to gene family expansion.DiscussionThese findings provide essential insights into the genetic diversity and evolutionary biology of C. pilosula, aiding its conservation and sustainable use.

Whole-genome assembly and annotation of the firecracker penstemon (Penstemon eatonii).

The penstemons are ornamental annual flowering plants native to the Intermountain West and Rocky Mountains and commonly used for urban landscaping. Elite commercial penstemons are generally susceptible to abiotic stresses, including drought, root rot, cold, and high salinity. Firecracker penstemon (Penstemon eatonii), however, is much more tolerant to these stresses than most elite cultivars. Importantly, firecracker penstemon has been reported to hybridize with many other penstemons and therefore provides the opportunity to develop more tolerant elite cultivars through strategic crossing. To facilitate the study and utilization of firecracker penstemon, we sequenced and annotated the genome of a P. eatonii accession collected from Utah, USA. We also performed low-coverage, whole-genome sequencing of 26 additional accessions from three different varieties of P. eatonii. This chromosome-scale genome assembly is the most contiguous and complete Penstemon genome sequenced to date.

New chromosome-scale genome assembly for the Dwarf Goby, a River Rhine invader

New chromosome-scale genome assembly for the Dwarf Goby, a River Rhine invader

A chromosome-scale genome assembly of the peanut beetle, Ulomoides dermestoides (Coleptera: Tenebrionidae)

Ulomoides dermestoides, a member of the Tenebrionidae family, is known by various common names such as Chinese beetle, Chinese weevil, peanut beetle, cancer beetle, or asthma beetle. Despite its origin in Asia, it has achieved global distribution owing to its perceived medicinal properties. In this study, we employed a combination of Illumina short-read, PacBio long-read, and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to produce a high-quality, chromosome-scale genome assembly of U. dermestoides. The assembled genome spans 253.38 Mb across 11 chromosomes, including the two sex chromosomes, X and Y. Key metrics include contig and scaffold N50 values of 9.83 Mb and 28.62 Mb, respectively, with an estimated genome completeness of 99.62%. Annotation efforts identified 15,553 protein-coding genes and revealed that 48.11% of the genome comprises repetitive sequences. This meticulously assembled genome serves as a pivotal resource for elucidating the genetic determinants underlying the species’ ecological habits and potential medicinal properties.

Haplotype-resolved genome assembly and resequencing provide insights into the origin and breeding of modern rose.

Modern rose (Rosa hybrida) is a recently formed interspecific hybrid and has become one of the most important and widely cultivated ornamentals. Here we report the haplotype-resolved chromosome-scale genome assembly of the tetraploid R. hybrida 'Samantha' ('JACmantha') and a genome variation map of 233 Rosa accessions involving various wild species, and old and modern cultivars. Homologous chromosomes of 'Samantha' exhibit frequent homoeologous exchanges. Population genomic and genomic composition analyses reveal the contributions of wild Rosa species to modern roses and highlight that R. odorata and its derived cultivars are important contributors to modern roses, much like R. chinensis 'Old Blush'. Furthermore, selective sweeps during modern rose breeding associated with major agronomic traits, including continuous and recurrent flowering, double flower, flower senescence and disease resistance, are identified. This study provides insights into the genetic basis of modern rose origin and breeding history, and offers unprecedented genomic resources for rose improvement.

Two high quality chromosome-scale genome assemblies of female and male silver pomfret (Pampus argenteus)

Pampus argenteus is a highly commercial marine fish whose population is declining sharply. Here, we generated a female P. argenteus genome, spanning 536.33 Mb with contig N50 of 1.79 Mb; 24070 genes (99.50% of 24,182) were functionally annotated. To improve quality of it, we assembled a 553.79 Mb genome of male fish with contig N50 of 24.75 Mb through HiFi and ultra-long ONT sequence technologies; 550.82 Mb were anchored onto 24 gap-free chromosomes; 22,892 genes (98.1% of 23,346) were functionally annotated; the QV value was 51.55 with 98.9% of BUSCO and 99.39% coverage of Illumina reads. Finally, we compared this genome with previous published one, revealing 37,301 SVs. 52.82 Mb and 18.05 Mb SDs were characterized in our and published assemblies, respectively, and 48.96 Mb PURs were constructed. Thus, this genome assembly exhibits excellent completeness, continuity and accuracy comparing to the published one, which can be current preferred reference genome. Overall, these works help aquaculture and wild resources recovery of P. argenteus and provide a valuable genetic resource for study.

Chromosome-scale genome assembly of the hunt bumble bee, Bombus huntii Greene, 1860, a species of agricultural interest.

The Hunt bumble bee, Bombus huntii, is a widely distributed pollinator in western North America. The species produces large colony sizes in captive rearing conditions, experiences low parasite and pathogen loads, and has been demonstrated to be an effective pollinator of tomatoes grown in controlled environment agriculture systems. These desirable traits have galvanized producer efforts to develop commercial Bombus huntii colonies for growers to deliver pollination services to crops. To better understand Bombus huntii biology and support population genetic studies and breeding decisions, we sequenced and assembled the Bombus huntii genome from a single haploid male. High-fidelity sequencing of the entire genome using PacBio, along with HiC sequencing, led to a comprehensive contig assembly of high continuity. This assembly was further organized into a chromosomal arrangement, successfully identifying 18 chromosomes spread across the 317.4 Mb assembly with a BUSCO score indicating 97.6% completeness. Synteny analysis demonstrates shared chromosome number (n = 18) with Bombus terrestris, a species belonging to a different subgenus, matching the expectation that presence of 18 haploid chromosomes is an ancestral trait at least between the subgenera Pyrobombus and Bombus sensu stricto. In conclusion, the assembly outcome, alongside the minimal tissue sampled destructively, showcases efficient techniques for producing a comprehensive, highly contiguous genome.

Chromosome-scale genome assembly of acerola (Malpighia emarginata DC.).

Acerola (Malpighia emarginata DC.) is a tropical evergreen shrub that produces vitamin C-rich fruits. Increasing fruit nutrition is one of the main targets of acerola breeding programs. Genomic tools have been shown to accelerate plant breeding even in fruiting tree species, which generally have a long life cycle; however, the availability of genomic resources in acerola, so far, has been limited. In this study, as a first step toward developing an efficient breeding technology for acerola, we established a chromosome-scale genome assembly of acerola using high-fidelity long-read sequencing and genetic mapping. The resultant assembly comprises 10 chromosome-scale sequences that span a physical distance of 1,032.5Mb and contain 35,892 predicted genes. Phylogenetic analysis of genome-wide SNPs in 60 acerola breeding materials revealed three distinct genetic groups. Overall, the genomic resource of acerola developed in this study, including its genome and gene sequences, genetic map, and phylogenetic relationship among breeding materials, will not only be useful for acerola breeding but will also facilitate genomic and genetic studies on acerola and related species.

Chromosome-scale genome assembly of Astragalus membranaceus using PacBio and Hi-C technologies

Astragalus membranaceus (Fisch.) Bge (AM) is a medicinal herb plant belonging to the Leguminosae family. In this study, we present a chromosome-scale genome assembly of AM, aiming to enhance the molecular biology and functional studies of Astragali Radix. The genome size of AM is about 1.43 Gb, with a contig N50 value of 1.67 Mb. A total of 98.16% of the assembly anchored to 9 pseudochromosomes using Hi-C technology. The assembly completeness was estimated to be 97.27% using BUSCO with the long terminal repeat assembly index (LAI) of 16.22 and quality value (QV) of 48.58. Additionally, the genome contained 67.98% repetitive sequences. Genome annotation predicted 29,914 protein-coding genes, including 73 genes involved in the flavonoid biosynthetic pathway and 2,048 transcription factors. The high-quality genome assembly and gene annotation resources will greatly facilitate future functional genomic studies in Leguminosae species.

A comparative analysis of planarian genomes reveals regulatory conservation in the face of rapid structural divergence

The planarian Schmidtea mediterranea is being studied as a model species for regeneration, but the assembly of planarian genomes remains challenging. Here, we report a high-quality haplotype-phased, chromosome-scale genome assembly of the sexual S2 strain of S. mediterranea and high-quality chromosome-scale assemblies of its three close relatives, S. polychroa, S. nova, and S. lugubris. Using hybrid gene annotations and optimized ATAC-seq and ChIP-seq protocols for regulatory element annotation, we provide valuable genome resources for the planarian research community and a first comparative perspective on planarian genome evolution. Our analyses reveal substantial divergence in protein-coding sequences and regulatory regions but considerable conservation within promoter and enhancer annotations. We also find frequent retrotransposon-associated chromosomal inversions and interchromosomal translocations within the genus Schmidtea and, remarkably, independent and nearly complete losses of ancestral metazoan synteny in Schmidtea and two other flatworm groups. Overall, our results suggest that platyhelminth genomes can evolve without syntenic constraints.

A reference genome for the Harpy Eagle reveals steady demographic decline and chromosomal rearrangements in the origin of Accipitriformes

The Harpy Eagle (Harpia harpyja) is an iconic species that inhabits forested landscapes in Neotropical regions, with decreasing population trends mainly due to habitat loss, and currently classified as vulnerable. Here, we report on a chromosome-scale genome assembly for a female individual combining long reads, optical mapping, and chromatin conformation capture reads. The final assembly spans 1.35 Gb, with N50scaffold equal to 58.1 Mb and BUSCO completeness of 99.7%. We built the first extensive transposable element (TE) library for the Accipitridae to date and identified 7,228 intact TEs. We found a burst of an unknown TE ~ 13–22 million years ago (MYA), coincident with the split of the Harpy Eagle from other Harpiinae eagles. We also report a burst of solo-LTRs and CR1 retrotransposons ~ 31–33 MYA, overlapping with the split of the ancestor to all Harpiinae from other Accipitridae subfamilies. Comparative genomics with other Accipitridae, the closely related Cathartidae and Galloanserae revealed major chromosome-level rearrangements at the basal Accipitriformes genome, in contrast to a conserved ancient genome architecture for the latter two groups. A historical demography reconstruction showed a rapid decline in effective population size over the last 20,000 years. This reference genome serves as a crucial resource for future conservation efforts towards the Harpy Eagle.

Chamaeleo calyptratus (veiled chameleon) chromosome-scale genome assembly and annotation provides insights into the evolution of reptiles and developmental mechanisms.

The family Chamaeleonidae comprises 228 species, boasting an extensive geographic spread and an array of evolutionary novelties and adaptations, but a paucity of genetic and molecular analyses. Veiled chameleon ( Chamaeleo calyptratus ) has emerged as a tractable research organism for the study of squamate early development and evolution. Here we report a chromosomal-level assembly and annotation of the veiled chameleon genome. We note a remarkable chromosomal conservation across squamates, but comparisons to more distant genomes reveal GC peaks correlating with ancestral chromosome fusion events. We subsequently identified the XX/XY region on chromosome 5, confirming environmental-independent sex determination in veiled chameleons. Furthermore, our analysis of the Hox gene family indicates that veiled chameleons possess the most complete set of 41 Hox genes, retained from an amniote ancestor. Lastly, the veiled chameleon genome has retained both ancestral paralogs of the Nodal gene, but is missing Dand5 and several other genes, recently associated with the loss of motile cilia during the establishment of left-right patterning. Thus, a complete veiled chameleon genome provides opportunities for novel insights into the evolution of reptilian genomes and the molecular mechanisms driving phenotypic variation and ecological adaptation.

Chromosome-scale genome assembly and annotation of the tetraploid potato cultivar Diacol Capiro adapted to the Andean region.

Potato (Solanum tuberosum) is an essential crop for food security and is ranked as the third most important crop worldwide for human consumption. The Diacol Capiro cultivar holds the dominant position in Colombian cultivation, primarily catering to the food processing industry. This highly heterozygous, autotetraploid cultivar belongs to the Andigenum group and it stands out for its adaptation to a wide variety of environments spanning altitudes from 1,800 to 3,200 meters above sea level. Here, a chromosome-scale assembly, referred to as DC, is presented for this cultivar. The assembly was generated by combining circular consensus sequencing with proximity ligation Hi-C for the scaffolding and represents 2.369 Gb with 48 pseudochromosomes covering 2,091 Gb and an anchor rate of 88.26%. The reference genome metrics, including an N50 of 50.5 Mb, a BUSCO (Benchmarking Universal Single-Copy Orthologue) score of 99.38%, and an Long Terminal Repeat Assembly Index score of 13.53, collectively signal the achieved high assembly quality. A comprehensive annotation yielded a total of 154,114 genes, and the associated BUSCO score of 95.78% for the annotated sequences attests to their completeness. The number of predicted NLR (Nucleotide-Binding and Leucine-Rich-Repeat genes) was 2107 with a large representation of NBARC (for nucleotide binding domain shared by Apaf-1, certain R gene products, and CED-4) containing domains (99.85%). Further comparative analysis of the proposed annotation-based assembly with high-quality known potato genomes, showed a similar genome metrics with differences in total gene numbers related to the ploidy status. The genome assembly and annotation of DC presented in this study represent a valuable asset for comprehending potato genetics. This resource aids in targeted breeding initiatives and contributes to the creation of enhanced, resilient, and more productive potato varieties, particularly beneficial for countries in Latin America.

Chromosome-scale genome assembly and annotation of Paspalum notatum Flüggé var. saurae

Paspalum notatum Flüggé is an economically important subtropical fodder grass that is widely used in the Americas. Here, we report a new chromosome-scale genome assembly and annotation of a diploid biotype collected in the center of origin of the species. Using Oxford Nanopore long reads, we generated a 557.81 Mb genome assembly (N50 = 56.1 Mb) with high gene completeness (BUSCO = 98.73%). Genome annotation identified 320 Mb (57.86%) of repetitive elements and 45,074 gene models, of which 36,079 have a high level of confidence. Further characterisation included the identification of 59 miRNA precursors together with their putative targets. The present work provides a comprehensive genomic resource for P. notatum improvement and a reference frame for functional and evolutionary research within the genus.

Chromosome-scale genome assembly of Docynia delavayi provides new insights into the α-farnesene biosynthesis

Docynia delavayi is an economically significant fruit species with a high market potential due to the special aroma of its fruit. Here, a 653.34 Mb high-quality genome of D. delavayi was first reported, of which 93.8 % of the sequences (612.98 Mb) could be anchored to 17 chromosomes, containing 48,325 protein-coding genes. Ks analysis proved that two whole genome duplication (WGD) events occurred in D. delavayi, resulting in the expansion of genes associated with terpene biosynthesis, which promoted its fruit-specific aroma production. Combined multi-omics analysis, α-farnesene was detected as the most abundant aroma substance emitted by D. delavayi fruit during storage, meanwhile one α-farnesene synthase gene (AFS) and 15 transcription factors (TFs) were identified as the candidate genes potentially involved in α-farnesene biosynthesis. Further studies for the regulation network of α-farnesene biosynthesis revealed that DdebHLH, DdeERF1 and DdeMYB could activate the transcription of DdeAFS. To our knowledge, it is the first report that MYB TF plays a regulatory role in α-farnesene biosynthesis, which will greatly facilitate future breeding programs for D. delavayi.

The Chrysosplenium sinicum genome provides insights into adaptive evolution of shade plants

Chrysosplenium sinicum, a traditional Tibetan medicinal plant, can successfully thrive in low-light environments for long periods of time. To investigate the adaptive evolution of shade plants in low-light environments, we generated a chromosome-scale genome assembly (~320 Mb) for C. sinicum by combining PacBio sequencing and Hi-C technologies. Based on our results, gene families related to photosynthesis and cell respiration greatly expanded and evolved in C. sinicum genome due to intracellular DNA transfer from organelle genome to nuclear genome. Under positive selective pressure, adaptive evolution of light-harvesting complex II (LHCII) component protein CsLhcb1s resulted in the expansion of threonine residues at the phosphorylation site of STN7 kinase, potentially establishing a crucial genomic foundation for enhancing C. sinicum’s adaptability in low-light environments. Through transcriptome and metabolome analysis, we identified chrysosplenol and chrysosplenoside as predominant flavonoid metabolites of C. sinicum and predicted their synthesis pathways. In addition, analysis of alternative splicing (AS) revealed that AS events help regulate state transition and flavonoid biosynthesis. The present study provides new insights into the genomes of shade plants exposed to low-light conditions and adaptive evolution of these genomes; in addition, the results improve our current knowledge on the biosynthetic and regulatory processes of chrysosplenol and chrysosplenoside.

Chromosome-scale pearl millet genomes reveal CLAMT1b as key determinant of strigolactone pattern and Striga susceptibility

The yield of pearl millet, a resilient cereal crop crucial for African food security, is severely impacted by the root parasitic weed Striga hermonthica, which requires host-released hormones, called strigolactones (SLs), for seed germination. Herein, we identify four SLs present in the Striga-susceptible line SOSAT-C88-P10 (P10) but absent in the resistant 29Aw (Aw). We generate chromosome-scale genome assemblies, including four gapless chromosomes for each line. The Striga-resistant Aw lacks a 0.7 Mb genome segment containing two putative CARLACTONOIC ACID METHYLTRANSFERASE1 (CLAMT1) genes, which may contribute to SL biosynthesis. Functional assays show that P10CLAMT1b produces the SL-biosynthesis intermediate methyl carlactonoate (MeCLA) and that MeCLA is the precursor of P10-specific SLs. Screening a diverse pearl millet panel confirms the pivotal role of the CLAMT1 section for SL diversity and Striga susceptibility. Our results reveal a reason for Striga susceptibility in pearl millet and pave the way for generating resistant lines through marker-assisted breeding or direct genetic modification.

Genomic and transcriptomic analyses of the elite rice variety Huizhan provide insight into disease resistance and heat tolerance

The indica rice variety Huizhan shows elite traits of disease resistance and heat tolerance. However, the underlying genetic basis of these traits is not fully understood due to limited genomic resources. Here, we used Nanopore long-read and next-generation sequencing technologies to generate a chromosome-scale genome assembly of Huizhan. Comparative genomics analysis uncovered a large chromosomal inversion and expanded gene families that are associated with plant growth, development and stress responses. Functional rice blast resistance genes, including Pi2, Pib and Ptr, and bacterial blight resistance gene Xa27, contribute to disease resistance of Huizhan. Furthermore, integrated genomics and transcriptomics analyses showed that OsHIRP1, OsbZIP60, the SOD gene family, and various transcription factors are involved in heat tolerance of Huizhan. The high-quality genome assembly and comparative genomics results presented in this study facilitate the use of Huizhan as an elite parental line in developing rice varieties adapted to disease pressure and climate challenges.