LTR_retriever: A Highly Accurate and Sensitive Program for Identification of Long Terminal Repeat Retrotransposons

Chromosome-Level Assembly and Annotation of the Endangered Red-Wing Fish (Distoechodon macrophthalmus)

Background/Objectives: The red-wing fish (Distoechodon macrophthalmus), an endangered species native to Yunnan, is endemic to Chenghai Lake. The natural population of this species has suffered a sharp decline due to the invasion of alien fish species. Fortunately, the artificial domestication and reproduction of D. macrophthalmus have been successful and this species has become an economic species locally. However, there is still little research on D. macrophthalmus. Methods: In this study, a high-quality genome of D. macrophthalmus was assembled and annotated. The genome was sequenced and assembled using the PacBio platform and Hi-C method. Results: The genome size is 1.01 Gb and N50 is 37.99 Mb. The assembled contigs were anchored into 24 chromosomes. BUSCO analysis revealed that the genome assembly has 95.6% gene coverage completeness. A total of 455.62 Mb repeat sequences (48.50% of the assembled genome) and 30,424 protein-coding genes were identified in the genome. Conclusions: This study provides essential genomic data for further research on the evolution and conservation of D. macrophthalmus. Meanwhile, the high-quality genome assembly also provides insights into the genomic evolution of the genus Distoechodon.

The tetraploid Camellia oleifera genome provides insights into evolution, agronomic traits, and genetic architecture of oil Camellia plants

Camellia oleifera is an economically important woody oil plant. Complex ploidy and lack of genomic information have seriously hindered the molecular breeding of C. oleifera. Here, we present an 11.43-Gb haplotype-resolved, chromosome-level genome assembly of tetraploid C. oleifera (COL-tetra). Methods employed in this study support the conclusion that COL-tetra is an autotetraploid and probably originates from genome doubling of the diploid C. brevistyla. In addition, DNA methylation plays a significant role in imbalanced allelic expression and seed development. Genetic divergence analyses reveal significant differentiation signals for flowering time between spring-flowering and autumn-flowering oil Camellia species, which probably account for reproductive isolation between species with distinct flowering times. Strong introgression signals are detected between COL-tetra and C. sasanqua and between C. vietnamensis and COL-hexa, which might affect the development of agronomic traits and environmental adaptability. This study provides valuable insights into the evolution, agronomic trait development, and genetic architecture of oil Camellia plants.

Genome assembly of an endemic butterfly (Minois Aurata) shed light on the genetic mechanisms underlying ecological adaptation to arid valley habitat

BACKGROUND

The Hengduan Mountains, one of the global biodiversity hotspots with exceptional species richness and high endemism, contains numerous arid valleys that create a distinctive geographical and ecological landscape. However, the adaptive evolutionary mechanisms of organism in the arid valley remain poorly understood. Minois aurata, an endemic butterfly species found exclusively in the arid valley of the upper Minjiang River, represents an attractive model system for studying adaptive evolutionary mechanisms to arid valley environments.

RESULTS

Here, we present the first chromosome-level genome assembly for Minois aurata, with a total size of approximately 609.17 Mb, and a scaffold N50 size of 23.88 Mb. These scaffolds were further clustered and anchored onto 29 chromosomes based on Hi-C data. A total of 16,163 protein-coding genes were predicted, of which 91.83% were functionally annotated. The expansion of transposable elements (TEs) accounts for the relatively large genome size of M. aurata, potentially aiding its adaptation to environmental conditions. Phylogenomic analyses based on 3,785 single-copy genes revealed that M. aurata is most closely related to Hipparchia semele. Further mitochondrial genome analysis of four Minois species placed M. aurata in a basal position within the genus, supporting it as an independent species. A total of 185 rapidly evolving and 232 specific gene families were identified in M. aurata. Functional enrichment analysis indicated that these gene families were mainly associated with ultraviolet radiation, heat and hypoxia responses. We also identified 234 positive selected genes in M. aurata, some of which are related to compound eye photoreceptor development, osmotic stress, and light stimulus response. Demographic analysis indicated that the effective population size of M. aurata decreased around 0.4 and 0.04 million years ago, respectively, coinciding with the localized sub-glaciation.

CONCLUSION

The chromosome-level genome offers a comprehensive genomic basis for understanding the evolutionary and adaptive strategies of Minois aurata in the unique arid valley environment of the Hengduan Mountains, while also providing valuable insights into the broader mechanisms of organism adaptation to such habitats.

Gap-free telomere-to-telomere haplotype assembly of the tomato hind (Cephalopholis sonnerati)

The tomato hind (Cephalopholis sonnerati) is an emerging economically important grouper in recent years. With the increasing maturity of sequencing technologies and assembly methodologies, a higher quality reference genome has become both accessible and necessary. In this study, we present two telomere-to-telomere (T2T) gap-free haplotype assemblies of the tomato hind with lengths of 1039.53 Mb (YSFRI_Csonn_HA_1.0, N50 43.83 Mb) and 1039.91 Mb (YSFRI_Csonn_HB_1.0, N50 44.09 Mb). Reads from next-generation sequencing, ONT ultra-long sequencing, and PacBio HiFi sequencing exhibited mapping rates exceeding 99.8% when aligned to these two assemblies. Evaluation using Merqury indicated high accuracy for both assemblies, with average quality values of 51.80 and 51.83, respectively. Percentages of 97.9% and 97.8% of complete BUSCOs were achieved, and a total of 23,270 and 23,184 protein-code genes were inferred in each assembly. Moreover, telomere identification, centromere prediction, and repetitive sequence annotation were also successfully performed. These two assemblies provide robust foundation for the genetic analysis and development of molecular genetic breeding technologies in C. sonnerati.

Chromosome scale genome assembly and annotation of coconut cultivar Chowghat Green Dwarf

The high-quality genome of coconut (Cocos nucifera L.) is a crucial resource for enhancing agronomic traits and studying genome evolution within the Arecaceae family. We sequenced the Chowghat Green Dwarf cultivar, which is resistant to the root (wilt) disease, utilizing Illumina, PacBio, ONT, and Hi-C technologies to produce a chromosome-level genome of ~ 2.68 Gb with a scaffold N50 of 174 Mb; approximately 97% of the genome could be anchored to 16 pseudo-molecules (2.62 Gb). In total, 34,483 protein-coding genes were annotated; the BUSCO completeness score was 96.80%, while the k-mer completeness was ~ 87%. The assembled genome includes 2.19 Gb (81.64%) of repetitive sequences, with long terminal repeats (LTRs) constituting the most abundant class at 53.76%. Additionally, our analysis confirms two whole-genome duplication (WGD) events in the C. nucifera lineage. A genome-wide analysis of LTR insertion time revealed ancient divergence and proliferation of copia and gypsy elements. In addition, 1368 RGAs were discovered in the CGD genome. We also developed a web server 'Kalpa Genome Resource' ( http://210.89.54.198:3000/ ), to manage and store a comprehensive array of genomic data, including genome sequences, genetic markers, structural and functional annotations like metabolic pathways, and transcriptomic profiles. The web server has an embedded genome browser to analyze and visualize the genome, its genomics elements, and transcriptome data. The in-built BLAST server allows sequence homology searches against genome, annotated transcriptome & proteome sequences. The genomic dataset and the database will support comparative genome analysis and can expedite genome-driven breeding and enhancement efforts for tapping genetic gains in coconut.

A first look at the genome structure of hexaploid "Mitcham" peppermint (Mentha × piperita L.)

Peppermint, Mentha × piperita L., is a hexaploid (2n = 6x = 72) and the predominant cultivar of commercial mint oil production in the US. This cultivar is threatened because of high susceptibility to the fungal disease verticillium wilt, caused by Verticillium dahliae. This report details the first draft polyploid chromosome-level genome assembly for this mint species. The "Mitcham" genome resource will broaden comparative studies of disease resistance, essential oil biosynthesis, and hybridization events within the genus Mentha. It will also be a valuable contribution to the body of phylogenetic studies involving Mentha and other genera that contain species with varying ploidy levels.

Genome evolution and diversity of wild and cultivated rice species

Wild species of crops serve as a valuable germplasm resource for breeding of modern cultivars. Rice (Oryza sativa L.) is a vital global staple food. However, research on genome evolution and diversity of wild rice species remains limited. Here, we present nearly complete genomes of 13 representative wild rice species. By integrating with four previously published genomes for pangenome analysis, a total of 101,723 gene families are identified across the genus, including 9834 (9.67%) core gene families. Additionally, 63,881 gene families absent in cultivated rice species but present in wild rice species are discovered. Extensive structural rearrangements, sub-genomes exchanges, widespread allelic variations, and regulatory sequence variations are observed in wild rice species. Interestingly, expanded but less diverse disease resistance genes in the genomes of cultivated rice, likely due to the loss of some resistance genes and the fixing and amplification of genes encoding resistance genes to specific diseases during domestication and artificial selection. This study not only reveals natural variations valuable for gene-level studies and breeding selection but also enhances our understanding on rice evolution and domestication.

Chromosome-level genome assembly and annotation of the Patagonian toothfish Dissostichus eleginoides

The Patagonian toothfish (Dissostichus eleginoides) belongs to the Actinopterygii class, and the suborder Notothenioidei, which lives in cold waters in the Southern Hemisphere. We performed assembly and annotation, and we integrated the Illumina short-read sequencing for polishinng, PacBio long-read sequencing for contig-level assembly, and Hi-C sequencing technology to obtain high-quality of chromosome-level genome assembly. The final assembly analysis resulted in a total of 495 scaffolds, a genome size of 844.7 Mbp and an N50 length of 36 Mbp. Among these data, we confirmed 24 scaffolds exceeded 10 Mbp and classified as chromosome-level. The completeness of BUSCO rate was over 97%. A total gene set of 32,224 was identified. Furthermore, we analyzed the presence of AFGP genes, classified into Antarctic and sub-Antarctic categories through phylogenetic analysis. This study provides a useful resource for the genomic analysis of Patagonian toothfish and genetic insights into the comparison with Antarctic fishes.

The chromosome-level genome assembly and annotation of an invasive forest pest Obolodiplosis robiniae

Biological invasion is a major global problem, leading to the loss of biodiversity and species extinction, and causing huge economic losses to countries. Obolodiplosis robiniae is a major invasive forest pest that has caused economic losses in Asia and Europe. Here, the chromosome- level genome of O. robiniae was assembled using the PacBio platform and Hi-C technology. A contig-level genome with a length of 199.49 Mb and a contig N50 of 4.66 Mb was assembled. Approximately 98.05% of contigs were successfully anchored to four chromosomes using Hi-C assisted genome assembly. The genome integrity was assessed to be 90.3% based on BUSCOs analysis. The high-quality genome provides valuable data for the study of invasive species, and a foundation for the understanding the biology and ecology of O. robiniae.

Comparative genome analysis and the genome-shaping role of long terminal repeat retrotransposons in the evolutionary divergence of fungal pathogens Blastomyces dermatitidis and Blastomyces gilchristii

Blastomyces dermatitidis and Blastomyces gilchristii are cryptic species of fungi that cause blastomycosis, an often severe disease involving pulmonary infection capable of systemic dissemination. While these species appear morphologically identical, differences exist in the genetic makeup, geographical range, and possibly the clinical presentation of infection. Here, we show genetic divergence between the cryptic species through both a Blastomyces species tree constructed from orthologous protein sequences and whole genome single-nucleotide variant phylogenomic analysis. Following linked-read sequencing and de novo genome assembly, we characterized and compared the genomes of 3 B. dermatitidis and 3 B. gilchristii isolates. The B. gilchristii genomes (73.25-75.4 Mb) were ∼8 Mb larger than the B. dermatitidis genomes (64.88-66.61 Mb). Average nucleotide identity was lower between genomes of different species than genomes of the same species, yet functional classification of genes suggested similar proteomes. The most striking difference involved long terminal repeat retrotransposons. Although the same retrotransposon elements were detected in the genomes, the quantity of elements differed between the 2 species. Gypsy retrotransposon content was significantly higher in B. gilchristii (38.04-39.26 Mb) than in B. dermatitidis (30.85-32.40 Mb), accounting for the majority of genome size difference between species. Age estimation and phylogenetic analysis of the reverse transcriptase domains suggested that these retrotransposons are relatively ancient, with genome insertion predating the speciation of B. dermatitidis and B. gilchristii. We postulate that different trajectories of genome contraction led to genetic incompatibility, reproductive isolation, and speciation, highlighting the role of transposable elements in fungal evolution.

Chromosome-level genome assembly of the ivory shell Babylonia areolata

The ivory shell Babylonia areolata is an economically important marine benthic gastropod known for its rapid growth and high nutritional value. B. areolata is distributed in Southeast Asia and the southeast coastal areas of China. In this study, we constructed a high-quality genome for B. areolata using PacBio, Illumina, and Hi-C sequencing technologies. The genome assembly comprised 35 chromosomal sequences with a total length of 1.65 Gb. The scaffold and contig N50 lengths were 53.17 Mb and 2.64 Mb, respectively, with repeat sequences constituting 64.46% of the genome. Furthermore, 26,130 protein-coding genes and 96.75% of the genome's BUSCOs were identified. This inaugural report of a B. areolata genome provides crucial foundational information for further investigations into the biology, genomics, and genetic improvement of economic traits of this species.

Chromosome-scale genome assembly of Apocynum pictum, a drought-tolerant medicinal plant from the Tarim Basin

Apocynum pictum Schrenk is a semishrub of the Apocynaceae family with a wide distribution throughout the Tarim Basin that holds significant ecological, medicinal, and economic values. Here, we report the assembly of its chromosome-level reference genome using Nanopore long-read, Illumina HiSeq paired-end, and high-throughput chromosome conformation capture sequencing. The final assembly is 225.32 Mb in length with a scaffold N50 of 19.64 Mb. It contains 23,147 protein-coding genes across 11 chromosomes, 21,148 of which (91.36%) have protein functional annotations. Comparative genomics analysis revealed that A. pictum diverged from the closely related species Apocynum venetum approximately 2.2 million years ago and has not undergone additional polyploidizations after the core eudicot WGT-γ event. Karyotype evolution analysis was used to characterize interchromosomal rearrangements in representative Apocynaceae species and revealed that several A. pictum chromosomes were derived entirely from single chromosomes of the ancestral eudicot karyotype. Finally, we identified 50 members of the well-known stress-responsive WRKY transcription factor family and used transcriptomic data to document changes in their expression at 2 stages of drought stress, identifying a number of promising candidate genes. Overall, this study provides high-quality genomic resources for evolutionary and comparative genomics of the Apocynaceae, as well as initial molecular insights into the drought adaptation of this valuable desert plant.

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

Introduction

Codonopsis 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.

Methods

This 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.

Results

The 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.

Discussion

These findings provide essential insights into the genetic diversity and evolutionary biology of C. pilosula, aiding its conservation and sustainable use.

Retrotransposon involves in photoperiodic spermatogenesis in Brandt's voles (Lasiopodomys brandtii) by co-transcription with flagellar genes

Photoperiod is a pivotal factor in affecting spermatogenesis in seasonal-breeding animals. Transposable elements have regulatory functions during spermatogenesis. However, whether it also functions in photoperiodic spermatogenesis in seasonal breeding animals is unknown. To explore this, we first annotated 5,501,822 transposons in the whole genome of Brandt's voles (Lasiopodomys brandtii), and revealed that LINEs were the most abundant, comprising 16.61 % of the genome. Following closely, SINEs accounted for 10.13 %, LTRs for 7.54 %, and DNA transposons for 0.70 %. Subsequently, we exposed male Brandt's voles to long-photoperiod (LP, 16 h/day) and short-photoperiod (SP, 8 h/day) from their embryonic stages, and obtained testes transcriptome at 4 and 10 weeks after birth. Differential expression and Pearson analysis indicated strongly positive correlations between the expression of differentially expressed retrotransposons and the adjacent genes. KO, KEGG and GSEA results showed that sperm flagellar genes were most enriched nearby the retrotransposons such as Dnah1, Dnah2, Dnah17, Dnali1. RT-PCR results showed that SINE/Alu_1213291 co-transcripted with Dnali1 gene. Our findings first reveal the regulatory function of transposons in photoperiodic spermatogenesis, providing insights into the role of photoperiod in seasonal reproduction in wild animals.

Integrated genomic, transcriptomic, and metabolomic analyses of Ilex hylonoma provide insights into the triterpenoid saponin biosynthesis

Ilex is known for its rich content of secondary metabolites, particularly triterpenoid saponins. These compounds hold significant value in natural remedies and herbal medicine. However, the molecular mechanisms responsible for triterpenoid biosynthesis in plants of this genus remain largely unexplored. In this study, we successfully generated the first chromosome-scale genome of Ilex hylonoma. The assembly, comprising 20 anchored chromosomes, has an N50 contig size of 2.13 Mb and a scaffold size of 33.68 Mb. Comparative genome analyses with two other congeners with available chromosome-level genomes suggested that an end-to-end chromosome fusion event likely contributed to the reduction in chromosome number from n = 20 to n = 19 within this genus. By integrating transcriptomic and metabolomic data, we identified the gene expression patterns and metabolite profiles of I. hylonoma across three commonly utilized medicinal tissues. We subsequently pinpointed candidate genes involved in the regulation of triterpenoid saponin biosynthesis, including CYP450 genes, UGT genes, and associated transcription factors. Furthermore, yeast heterologous expression analysis revealed that ihyl08363 catalyzed the conversion of β-amyrin into oleanolic acid, while ihyl04303 catalyzed the C-2α hydroxylation of oleanolic acid to produce maslinic acid. This integrated analysis provides valuable insights into the biosynthesis of important triterpenoid saponins in medicinal Ilex plants.

Convergence and divergence of diploid and tetraploid cotton genomes

Polyploidy is an important driving force in speciation and evolution; however, the genomic basis for parallel selection of a particular trait between polyploids and ancestral diploids remains unexplored. Here we construct graph-based pan-genomes for diploid (A2) and allotetraploid (AD1) cotton species, enabled by an assembly of 50 genomes of genetically diverse accessions. We delineate a mosaic genome map of tetraploid cultivars that illustrates genomic contributions from semi-wild forms into modern cultivars. Pan-genome comparisons identify syntenic and hyper-divergent regions of continued variation between diploid and tetraploid cottons, and suggest an ongoing process of sequence evolution potentially linked to the contrasting genome size change in two subgenomes. We highlight 43% of genetic regulatory relationships for gene expression in diploid encompassing sequence divergence after polyploidy, and specifically characterize six underexplored convergent genetic loci contributing to parallel selection of fiber quality. This study offers a framework for pan-genomic dissection of genetic regulatory components underlying parallel selection of desirable traits in organisms.

The chromosome-level Elaeagnus mollis genome and transcriptomes provide insights into genome evolution, glycerolipid and vitamin E biosynthesis in seeds

Elaeagnus mollis, which has seeds with high lipid and vitamin E contents, is a valuable woody oil plant with potential for utilization. Currently, the biosynthesis and regulation mechanism of glycerolipids and vitamin E are still unknown in E. mollis. Here, we present the chromosome-level reference genome of E. mollis (scaffold N50: ~40.66Mbp, genome size: ~591.48Mbp) by integrating short-read, long-read, and Hi-C sequencing platforms. A total of 36,796 protein-coding sequences, mainly located on 14 proto-chromosomes, were predicted. Additionally, two whole genome duplication (WGD) events were suggested to have occurred ~54.07 and ~35.06 million years ago (MYA), with Elaeagnaceae plants probably experiencing both WGD events. Furthermore, the long terminal retrotransposons in E. mollis were active ~0.23MYA, and one of them was inferred to insert into coding sequence of the negative regulatory lipid synthesis gene, EMF2. Through transcriptomic and metabonomic analysis, key genes contributing to the high lipid and vitamin E levels of E. mollis seeds were identified, while miRNA regulation was also considered. This comprehensive work on the E. mollis genome not only provides a solid theoretical foundation and experimental basis for the efficient utilization of seed lipids and vitamin E, but also contributes to the exploration of new genetic resources.

T2T genomes of carrot and Alternaria dauci and their utility for understanding host-pathogen interactions during carrot leaf blight disease

Carrot (Daucus carota) is one of the most popular and nutritious vegetable crops worldwide. However, significant yield losses occur every year due to leaf blight, a disease caused by a fungal pathogen (Alternaria dauci). Past research on resistance to leaf blight disease in carrots has been slow because of the low-quality genome assemblies of both carrot and the pathogen. Here, we report the greatly improved assemblies and annotations of telomere-to-telomere (T2T) reference genomes of carrot DH13M14 (451.04 Mb) and A. dauci A2016 (34.91 Mb). Compared with the previous carrot genome versions, our assembly featured notable improvements in genome size, continuity, and completeness of centromeres and telomeres. In addition, we generated a time course transcriptomic atlas during the infection of carrots by A. dauci and captured their dynamic gene expression reprogramming during the interaction process. During infection, A. dauci genes encoding effectors and enzymes responsible for the degradation of plant cell wall components, e.g., cellulose and pectin, were identified, which appeared to increase pathogenic ability through upregulation. In carrot, the coordinated gene expression of components of pattern- and effector-triggered immunity (PTI and ETI) in response to A. dauci attack was characterized. The biosynthesis or signal transduction of plant hormones, including JA, SA, and ethylene, was also involved in the carrot response to A. dauci. This work provides a foundation for understanding A. dauci pathogenic progression and carrot defense mechanisms to improve carrot resistance to leaf blight disease. The Carrot Database (CDB) developed also provides a useful resource for the carrot community.

Haplotype-resolved genome assembly and resequencing analysis provide insights into genome evolution and allelic imbalance in Pinus densiflora

Haplotype-level allelic characterization facilitates research on the functional, evolutionary and breeding-related features of extremely large and complex plant genomes. We report a 21.7-Gb chromosome-level haplotype-resolved assembly in Pinus densiflora. We found genome rearrangements involving translocations and inversions between chromosomes 1 and 3 of Pinus species and a proliferation of specific long terminal repeat (LTR) retrotransposons (LTR-RTs) in P. densiflora. Evolutionary analyses illustrated that tandem and LTR-RT-mediated duplications led to an increment of transcription factor (TF) genes in P. densiflora. The haplotype sequence comparison showed allelic imbalances, including presence-absence variations of genes (PAV genes) and their functional contributions to flowering and abiotic stress-related traits in P. densiflora. Allele-aware resequencing analysis revealed PAV gene diversity across P. densiflora accessions. Our study provides insights into key mechanisms underlying the evolution of genome structure, LTR-RTs and TFs within the Pinus lineage as well as allelic imbalances and diversity across P. densiflora.

Genome-wide characterization and evolution analysis of miniature inverted-repeat transposable elements in Barley (Hordeum vulgare)

Miniature inverted-repeat transposable elements (MITEs) constitute a class of class II transposable elements (TEs) that are abundant in plant genomes, playing a crucial role in their evolution and diversity. Barley (Hordeum vulgare), the fourth-most important cereal crop globally, is widely used for brewing, animal feed, and human consumption. However, despite their significance, the mechanisms underlying the insertion or amplification of MITEs and their contributions to barley genome evolution and diversity remain poorly understood. Through our comprehensive analysis, we identified 32,258 full-length MITEs belonging to 2,992 distinct families, accounting for approximately 0.17% of the barley genome. These MITE families can be grouped into four well-known superfamilies (Tc1/Mariner-like, PIF/Harbinger-like, hAT-like, and Mutator-like) and one unidentified superfamily. Notably, we observed two major expansion events in the barley MITE population, occurring approximately 12-13 million years ago (Mya) and 2-3 Mya. Our investigation revealed a strong preference of MITEs for gene-related regions, particularly in promoters, suggesting their potential involvement in regulating host gene expression. Additionally, we discovered that 7.73% miRNAs are derived from MITEs, thereby influencing the origin of certain miRNAs and potentially exerting a significant impact on post-transcriptional gene expression control. Evolutionary analysis demonstrated that MITEs exhibit lower conservation compared to genes, consistent with their dynamic mobility. We also identified a series of MITE insertions or deletions associated with domestication, highlighting these regions as promising targets for crop improvement strategies. These findings significantly advance our understanding of the fundamental characteristics and evolutionary patterns of MITEs in the barley genome. Moreover, they contribute to our knowledge of gene regulatory networks and provide valuable insights for crop improvement endeavors.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Chromosome-scale genome assembly of the mangrove climber species Dalbergia candenatensis

Consisting of trees, climbers and herbs exclusively in the intertidal environments, mangrove forest is one of the most extreme and vulnerable ecosystems of our planet and has long been of great interest for biologists and ecologists. Here, we first assembled the chromosome-scale genome of a climber mangrove plant, Dalbergia candenatensis. The assembled genome size is approximately 474.55 Mb, with a scaffold N50 of 48.1 Mb, a complete BUSCO score of 98.4%, and a high LTR Assembly Index value of 21. The genome contained 283.46 Mb (59.74%) repetitive sequences, and 29,554 protein-coding genes were predicted, of which 87.54% were functionally annotated in five databases. The high-quality genome assembly and annotation presented herein provide a valuable genomic resource that will expedite genomic and evolutionary studies of mangrove plants and facilitate the elucidation of molecular mechanisms underlying the salt- and water-logging-tolerance of mangrove plants.

Chromosome level genome assembly of giant freshwater prawn (Macrobrachium rosenbergii)

The giant freshwater prawn (Macrobrachium rosenbergii) has many advantages in aquaculture, such as fast growth rate, short breeding cycle and good nutrition, which makes it a freshwater shrimp with high economic value. Herein, high-quality chromosome-level genome of both female and male prawns were obtained by combining Illumina paired-end sequencing, PacBio single molecule sequencing technique and High-through chromosome conformation capture (Hi-C) technologies. In ZZ male prawn, a final contig assembly of 3118.58 Mb with a N50 length of 956,237 bp was obtained. In WW female prawn, a final contig assembly of 3333.31 Mb with a N50 length of 1,143,555 bp was obtained. The assembled genome sequences from prawns were anchored to 59 chromosomes. Moreover, the sex chromosomes including W chromosome and Z chromosome were generated in prawn with the length of 36.23 Mb and 27.33 Mb, respectively. The sequence similarity of Z chromosome and W chromosome reached to 74.90%. The high-quality genome resource will be useful for further molecular breeding and functional genomic research of giant freshwater prawns.

Inversions encounter relaxed genetic constraints and balance birth and death of TPS genes in Curcuma

Evolutionary dynamics of inversion and its impact on biochemical traits are a puzzling question. Here, we show abundance of inversions in three Curcuma species (turmeric, hidden ginger and Siam tulip). Genes within inversions display higher long terminal repeat content and lower expression level compared with genomic background, suggesting inversions in Curcuma experience relaxed genetic constraints. It is corroborated by depletion of selected SNPs and enrichment of deleterious mutations in inversions detected among 56 Siam tulip cultivars. Functional verification of tandem duplicated terpene synthase (TPS) genes reveals that genes within inversions become pseudogenes, while genes outside retain catalytic function. Our findings suggest that inversions act as a counteracting force against tandem duplication in balancing birth and death of TPS genes and modulating terpenoid contents in Curcuma. This study provides an empirical example that inversions are likely not adaptive but affect biochemical traits.

Global characterization of somatic mutations and DNA methylation changes during vegetative propagation in strawberries

Somatic mutations arise and accumulate during tissue culture and vegetative propagation, potentially affecting various traits in horticultural crops, but their characteristics are still unclear. Here, somatic mutations in regenerated woodland strawberry derived from tissue culture of shoot tips under different conditions and 12 cultivated strawberry individuals are analyzed by whole genome sequencing. The mutation frequency of single nucleotide variants is significantly increased with increased hormone levels or prolonged culture time in the range of 3.3 × 10-8-3.0 × 10-6 mutations per site. CG methylation shows a stable reduction (0.71%-8.03%) in regenerated plants, and hypoCG-DMRs are more heritable after sexual reproduction. A high-quality haplotype-resolved genome is assembled for the strawberry cultivar "Beni hoppe." The 12 "Beni hoppe" individuals randomly selected from different locations show 4731-6005 mutations relative to the reference genome, and the mutation frequency varies among the subgenomes. Our study has systematically characterized the genetic and epigenetic variants in regenerated woodland strawberry plants and different individuals of the same strawberry cultivar, providing an accurate assessment of somatic mutations at the genomic scale and nucleotide resolution in plants.

Comparative Genomics Reveals Species-Specific Genes and Symbiotic Adaptations in Tricholoma matsutake

Tricholoma matsutake, a highly valued ectomycorrhizal fungus, requires a symbiotic relationship with pine trees for growth, complicating its cultivation. This study presents a comprehensive comparative genomic analysis of Tricholoma species, with a focus on T. matsutake. Genomic data from 19 assemblies representing 13 species were analyzed to identify genus-, species-, and strain-specific genes, revealing significant evolutionary adaptations. Notably, T. matsutake exhibits a higher proportion of repetitive elements compared to other species, with retrotransposons like LTR Gypsy dominating its genome. Phylogenomic analyses showed that T. matsutake forms a monophyletic group closely related to T. bakamatsutake. Gene family expansion and contraction analyses highlighted the unique evolutionary pressures on T. matsutake, particularly the loss of tryptophan-related metabolic pathways and the gain of genes related to iron ion homeostasis, which may be crucial for its adaptation to nutrient-limited environments. Additionally, the reduction in secreted proteins and carbohydrate-active enzymes reflects the host-dependent lifestyle of T. matsutake and related species. These findings enhance our understanding of the genetic and evolutionary mechanisms underlying the complex symbiotic relationships of T. matsutake, offering potential avenues for optimizing its cultivation and commercial value.

Chromosome-level genome assembly of Cyamophila willieti (Hemiptera: Psyllidae)

Cyamophila willieti (Hemiptera: Psyllidae) is a significant pest that adversely affects the growth of Styphnolobium japonicum and its variant, Styphnolobium japonicum f. pendula. Despite its impact, research on this species remains limited. In this study, we successfully assembled a chromosome-level genome for Cyamophila willieti using a comprehensive approach that integrated Illumina sequencing, PacBio sequencing, and Hi-C technology. The genome size was determined to be 361.61 Mb, with a scaffold N50 length of 28.90 Mb. Additionally, the genome was mapped to 13 chromosomes (N = 12 A + X), and a total of 15,841 genes were predicted, with 90.38% of them functionally annotated. In summary, the high-quality genome of Cyamophila willieti provides valuable data to support further research, including pest management strategies.

Genomic insights into adaptive evolution of the species-rich cosmopolitan plant genus Rhododendron

The species-rich cosmopolitan genus Rhododendron offers a good system for exploring the genomic mechanisms underlying adaptation to diverse habitats. Here, we report high-quality chromosomal-level genome assemblies of nine species, representing all five subgenera, different altitudinal distributions, and all flower color types of this genus. Further comprehensive genomic analyses indicate diverse adaptive strategies employed by Rhododendron, particularly adaptation to alpine and subalpine habitats by expansion/contraction of gene families involved in pathogen defense and oxidative phosphorylation, genomic convergent evolution, and gene copy-number variation. The convergent adaptation to high altitudes is further shown by population genomic analysis of R. nivale from the Himalaya-Hengduan Mountains. Moreover, we identify the genes involved in the biosynthesis of anthocyanins and carotenoids, which play a crucial role in shaping flower color diversity and environmental adaptation. Our study is significant for comprehending plant adaptive evolution and the uneven distribution of species diversity across different geographical regions.

SesamumGDB: a comprehensive platform for Sesamum genetics and genomics analysis

Sesame (Sesamum indicum L., 2n = 26) is a crucial oilseed crop cultivated worldwide. The ancient evolutionary position of the Sesamum genus highlights its value for genomics and molecular genetics research among the angiosperms of other genera. However, Sesamum is considered a small orphan genus with only a few genomic databases for cultivated sesame to date. The urgent need to construct comprehensive, curated genome databases that include genus-specific gene resources for both cultivated and wild Sesamum species is being recognized. In response, we developed Sesamum Genomics Database (SesamumGDB), a user-friendly genomic database that integrates extensive genomic resources from two cultivated sesame varieties (S. indicum) and seven wild Sesamum species, covering all three chromosome groups (2n = 26, 32, and 64). This database showcases a total of 352 471 genes, including 6026 related to lipid metabolism and 17 625 transcription factors within Sesamum. Equipped with an array of bioinformatics tools such as BLAST (basic local alignment search tool) and JBrowse (the Javascript browser), SesamumGDB facilitates data downloading, screening, visualization, and analysis. As the first centralized Sesamum genome database, SesamumGDB offers extensive insights into the genomics and genetics of sesame, potentially enhancing the molecular breeding of sesame and other oilseed crops in the future. Database URL: http://www.sgbdb.com/sgdb/.

Genome sequence of a European Diplocarpon coronariae strain and in silico structure of the mating-type locus

Diplocarpon coronariae is a fungal pathogen that is prevalent in low-input apple production. Over the past 15 years, it has become increasingly distributed in Europe. However, comprehensive insights into its biology and pathogenicity remain limited. One particular aspect is the rarity of the sexual morph of this pathogen, a phenomenon hitherto unobserved in Europe. Diplocarpon coronariae reproduces through a heterothallic mating system requiring at least two different mating types for sexual reproduction. Genes determining the mating types are located on the mating-type locus. In this study, D. coronariae strain DC1_JKI from Dresden, Germany, was sequenced and used to unravel the structure of the mating type locus. Using short-read and long-read sequencing methods, the first gapless and near-complete telomere-to-telomere genome assembly of D. coronariae was achieved. The assembled genome spans 51.2 Mbp and comprises 21 chromosome-scale contigs of high completeness. The generated genome sequence was used to in silico elucidate the structure of the mating-type locus, identified as MAT1-2. Furthermore, an examination of MAT1-1 and MAT1-2 frequency across a diverse set of samples sourced from Europe and Asia revealed the exclusive presence of MAT1-2 in European samples, whereas both MAT loci were present in Asian counterparts. Our findings suggest an explanation for the absence of the sexual morph, potentially linked to the absence of the second mating idiomorph of D. coronariae in European apple orchards.

Chromosome-level genome assembly of two cultivated Jujubes

Jujube (Ziziphus jujuba Mill.) is a valuable tree species with economic, edible, medicinal, and ecological conservation benefits. In this study, we used PacBio HiFi sequencing and Hi-C technology to assemble chromosome-level genomes of two cultivated Jujubes, namely 'Lingwuchangzao' and 'Shiguang'. The genome sizes of 'Lingwuchangzao' and 'Shiguang' were 385.66 Mb and 394.12 Mb, respectively, with contig N50 sizes of 30.62 Mb and 4.30 Mb. These genomes contained 31,082 and 31,015 protein-coding genes, with repeat element contents of 42.11% and 42.33%, respectively. Phylogenetic analysis revealed that 'Lingwuchangzao' was sister to 'Shiguang' and followed by 'Dongzao'. Additionally, comparative analysis of gene families among 'Lingwuchangzao', 'Shiguang', 'Dongzao', 'Junzao', and 'Suanzao' identified 15,988 (57.98%) core gene families, 11,191 (40.59%) disposable gene families, and 394 (1.43%) private gene families. Overall, the assembly of the genomes of these two cultivated Jujube varieties provides valuable genetic information for Jujube evolution, functional genomics research, and molecular breeding.

T2T genome assemblies of Fallopia multiflora (Heshouwu) and F. multiflora var. angulata

The traditional Chinese medicinal plant Fallopia multiflora (hereafter AYY) is well known for its anti-hyperlipidaemia, immunomodulating, and hepatoprotective effects, attributed to its abundance of anthraquinones and stilbene glycosides, which are distinct to its variety F. multiflora var. angulata (hereafter CYY) in proportion and composition. In this study, telomere-to-telomere (T2T) genomes were assembled for AYY and CYY using PacBio HiFi reads and Hi-C data. The genome sizes, percentages of repetitive sequences, and numbers of protein-coding genes of AYY and CYY assemblies were 1,458.37 Mb/70.48%/84,768 and 1,174.38 Mb/67.36%/69,100, respectively. Comprehensive assessments confirmed high continuity (contig N50: 112.58 Mb and 94.83 Mb; number of gaps: 9 and 5), completeness (BUSCOs: 97.30% and 97.60%; LAI: 16.93 and 16.77), and correctness (QV: 51.42 and 52.60) of AYY and CYY assemblies. These T2T genomes of F. multiflora provide valuable resources for studying the biosynthesis of specialized metabolites and facilitating precise genetic improvement.

Genomic diversity and evolution of the Hawaiian Islands endemic Kokia (Malvaceae)

Island species are highly vulnerable due to habitat destruction and their often small population sizes with reduced genetic diversity. The Hawaiian Islands constitute the most isolated archipelago on the planet, harboring many endemic species. Kokia is an endangered flowering plant genus endemic to these islands, encompassing 3 extant and 1 extinct species. Recent studies provided evidence of unexpected genetic diversity within Kokia. Here, we provide high-quality genome assemblies for all 3 extant Kokia species, including an improved genome for Kokia drynarioides. All 3 Kokia genomes contain 12 chromosomes exhibiting high synteny within and between Kokia and the sister taxon Gossypioides kirkii. Gene content analysis revealed a net loss of genes in K. cookei compared to other species, whereas the gene complement in K. drynarioides remains stable and that of Kokia kauaiensis displays a net gain. A dated phylogeny estimates the divergence time from the last common ancestor for the 3 Kokia species at ∼1.2 million years ago (mya), with the sister taxa (K. cookei + K. drynarioides) diverging ∼0.8 mya. Kokia appears to have followed a stepping-stone pattern of colonization and diversification of the Hawaiian archipelago, likely starting on low or now submerged older islands. The genetic resources provided may benefit conservation efforts of this endangered endemic genus.

Haplotype-resolved and near-T2T genome assembly of the African catfish (Clarias gariepinus)

Airbreathing catfish are stenohaline freshwater fish capable of withstanding various environmental conditions and farming practices, including breathing atmospheric oxygen. This unique ability has enabled them to thrive in semi-terrestrial habitats. However, the genomic mechanisms underlying their adaptation to adverse ecological environments remain largely unexplored, primarily due to the limited availability of high-quality genomic resources. Here, we present a haplotype-resolved and near telomere-to-telomere (T2T) genome assembly of the African catfish (Clarias gariepinus), utilizing Oxford Nanopore, PacBio HiFi, Illumina and Hi-C sequencing technologies. The primary assembly spans 969.62 Mb with only 47 contigs, achieving a contig N50 of 33.71 Mb. Terminal telomeric signals were detected in 22 of 47 contigs, suggesting T2T assembled chromosomes. BUSCO analysis confirmed gene space completeness of 99% against the Actinopterygii dataset, highlighting the high quality of the assembly. Genome annotation identified 25,655 protein-coding genes and estimated 43.94% genome-wide repetitive elements. This data provides valuable genomic resources to advance aquaculture practices and to explore the genomic underpinnings of the ecological resilience of airbreathing catfish and related teleosts.

Chromosome-level genome assembly of the sacoglossan sea slug Elysia timida (Risso, 1818)

BACKGROUND

Sequencing and annotating genomes of non-model organisms helps to understand genome architecture, the genetic processes underlying species traits, and how these genes have evolved in closely-related taxa, among many other biological processes. However, many metazoan groups, such as the extremely diverse molluscs, are still underrepresented in the number of sequenced and annotated genomes. Although sequencing techniques have recently improved in quality and quantity, molluscs are still neglected due to difficulties in applying standardized protocols for obtaining genomic data.

RESULTS

In this study, we present the chromosome-level genome assembly and annotation of the sacoglossan sea slug species Elysia timida, known for its ability to store the chloroplasts of its food algae. In particular, by optimizing the long-read and chromosome conformation capture library preparations, the genome assembly was performed using PacBio HiFi and Arima HiC data. The scaffold and contig N50s, at 41.8 Mb and 1.92 Mb, respectively, are approximately 30-fold and fourfold higher compared to other published sacoglossan genome assemblies. Structural annotation resulted in 19,904 protein-coding genes, which are more contiguous and complete compared to publicly available annotations of Sacoglossa with respect to metazoan BUSCOs. We found no evidence for horizontal gene transfer (HGT), i.e. no photosynthetic genes encoded in the sacoglossan nucleus genome. However, we detected genes encoding polyketide synthases in E. timida, indicating that polypropionates are produced. HPLC-MS/MS analysis confirmed the presence of a large number of polypropionates, including known and yet uncharacterised compounds.

CONCLUSIONS

We can show that our methodological approach helps to obtain a high-quality genome assembly even for a "difficult-to-sequence" organism, which may facilitate genome sequencing in molluscs. This will enable a better understanding of complex biological processes in molluscs, such as functional kleptoplasty in Sacoglossa, by significantly improving the quality of genome assemblies and annotations.

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.

Biparental graph strategy to represent and analyze hybrid plant genomes

Hybrid plants are found extensively in the wild, and they often demonstrate superior performance of complex traits over their parents and other selfing plants. This phenomenon, known as heterosis, has been extensively applied in plant breeding for decades. However, the process of decoding hybrid plant genomes has seriously lagged due to the challenges associated with genome assembly and the lack of appropriate methodologies for their subsequent representation and analysis. Here, we present the assembly and analysis of 2 hybrids, an intraspecific hybrid between 2 maize (Zea mays ssp. mays) inbred lines and an interspecific hybrid between maize and its wild relative teosinte (Z. mays ssp. parviglumis), utilizing a combination of PacBio High Fidelity sequencing and chromatin conformation capture sequencing data. The haplotypic assemblies are well phased at chromosomal scale, successfully resolving the complex loci with extensive parental structural variations (SVs). By integrating into a biparental genome graph, the haplotypic assemblies can facilitate downstream short-read-based SV calling and allele-specific gene expression analysis, demonstrating outstanding advantages over a single linear genome. Our work offers a comprehensive workflow that aims to facilitate the decoding of numerous hybrid plant genomes, particularly those with unknown or inaccessible parentage, thereby enhancing our understanding of genome evolution and heterosis.

De Novo Assembly and Annotation of the Siganus fuscescens (Houttuyn, 1782) Genome: Marking a Pioneering Advance for the Siganidae Family

This study presents the first draft genome of Siganus fuscescens, and thereby establishes the first whole-genome sequence for a species in the Siganidae family. Leveraging both long and short read sequencing technologies, i.e., Oxford Nanopore and Illumina sequencing, we successfully assembled a mitogenome spanning 16.494 Kb and a first haploid genome encompassing 498 Mb. The assembled genome accounted for a 99.6% of the estimated genome size and was organized into 164 contigs with an N50 of 7.2 Mb. This genome assembly showed a GC content of 42.9% and a high Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness score of 99.5% using actinopterygii_odb10 lineage, thereby meeting stringent quality standards. In addition to its structural aspects, our study also examined the functional genomics of this species, including the intricate capacity to biosynthesize long-chain polyunsaturated fatty acids (LC-PUFAs) and secrete venom. Notably, our analyses revealed various repeats elements, which collectively constituted 17.43% of the genome. Moreover, annotation of 28,351 genes uncovered both shared genetic signatures and those that are unique to S. fuscescens. Our assembled genome also displayed a moderate prevalence of gene duplication compared to other fish species, which suggests that this species has a distinctive evolutionary trajectory and potentially unique functional constraints. Taken altogether, this genomic resource establishes a robust foundation for future research on the biology, evolution, and the aquaculture potential of S. fuscescens.

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological, and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb 'Qinjiao'. Herein, we present its first chromosome-level genome sequence assembly and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5 Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70 Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

Differences in activity and stability drive transposable element variation in tropical and temperate maize

Much of the profound interspecific variation in genome content has been attributed to transposable elements (TEs). To explore the extent of TE variation within species, we developed an optimized open-source algorithm, panEDTA, to de novo annotate TEs in a pangenome context. We then generated a unified TE annotation for a maize pangenome derived from 26 reference-quality genomes, which reveals an excess of 35.1 Mb of TE sequences per genome in tropical maize relative to temperate maize. A small number (n = 216) of TE families, mainly LTR retrotransposons, drive these differences. Evidence from the methylome, transcriptome, LTR age distribution, and LTR insertional polymorphisms reveals that 64.7% of the variability is contributed by LTR families that are young, less methylated, and more expressed in tropical maize, whereas 18.5% is driven by LTR families with removal or loss in temperate maize. Additionally, we find enrichment for Young LTR families adjacent to nucleotide-binding and leucine-rich repeat (NLR) clusters of varying copy number across lines, suggesting TE activity may be associated with disease resistance in maize.

Genome of the most noxious weed water hyacinth (Eichhornia crassipes) provides insights into plant invasiveness and its translational potential

The invasive character of Eichhornia crassipes (water hyacinth) is a major threat to global biodiversity and ecosystems. To investigate the genomic basis of invasiveness, we performed the genome and transcriptome sequencing of E. crassipes and reported the genome of 1.11 Gbp size with 63,299 coding genes and N50 of 1.98 Mb. We confirmed a recent whole genome duplication event in E. crassipes that resulted in high intraspecific collinearity and significant expansion in gene families. Further, the orthologs gene clustering analysis and comparative evolutionary analysis with 14 other aquatic invasive and non-invasive angiosperm species revealed adaptive evolution in genes associated with plant-pathogen interaction, hormone signaling, abiotic stress tolerance, heavy metals sequestration, photosynthesis, and cell wall biosynthesis with highly expanded gene families, which contributes toward invasive characteristics of the water hyacinth. However, these characteristics also make water hyacinth an excellent candidate for biofuel production, phytoremediation, and other translational applications.

Chromosome-level genome assembly of Chinese water Scorpion Ranatra chinensis (Heteroptera: Nepidae)

Heteroptera (the true bugs), one of the most diverse lineages of insects, diversified in feeding strategies and living habitats, and thus become an ideal lineage for studies on adaptive evolution. Chinese water scorpion Ranatra chinensis (Heteroptera: Nepidae) is a predaceous bug living in lentic water systems, representing an ideal model for studying habitat transition and adaptation to water environment. However, genetic studies on this water bug remain limited. Here, we obtained a chromosome-level genome of R. chinensis using PacBio HiFi long reads and Hi-C sequencing reads. The total assembly size of genome is 867.89 Mb, with a scaffold N50 length of 26.48 Mb and the GC content of 39.50%. All contigs were assembled into 23 pseudo-chromosomes (N = 19 A + X1X2X3X4), and we predicted 18,424 protein-coding genes in this genome. This study will provide valuable genomic resources for future studies on the biology, water adaptation, and genome evolution of water bugs.

Chromosome-level genome assembly of cotton thrips Thrips tabaci (Thysanoptera: Thripidae)

Cotton thrip, Thrips tabaci is a major polyphagous pest widely distributed on a variety of crops around the world, causing huge economic losses to agricultural production. Due to its biological and genomic characteristics, this pest can reproduce quickly and develop resistance to various pesticides in a very short time. However, the lack of high-quality reference genomes has hindered deeper gene function exploration and slows down the development of new management strategies. Here, we assembled a high-quality genome of T. tabaci at the chromosome level for the first time by using Illumina, PacBio long reads, and Hi-C technologies. The 329.59 Mb genome was obtained from 320 contigs, with a contig N50 of 1.53 Mb, and 94.21% of the assembly was anchored to 18 chromosomes. In total, 17,816 protein-coding genes were annotated, and 96.78% of BUSCO genes were fully represented. In conclusion, this high-quality genome provides a valuable genetic basis for our understanding of the biology of T. tabaci and contributes to the development of management strategies for cotton thrip.

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion

BACKGROUND

Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions.

RESULTS

Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP".

CONCLUSION

The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

The genomes of seven economic Caesalpinioideae trees provide insights into polyploidization history and secondary metabolite biosynthesis

The Caesalpinioideae subfamily contains many well-known trees that are important for economic sustainability and human health, but a lack of genomic resources has hindered their breeding and utilization. Here, we present chromosome-level reference genomes for the two food and industrial trees Gleditsia sinensis (921 Mb) and Biancaea sappan (872 Mb), the three shade and ornamental trees Albizia julibrissin (705 Mb), Delonix regia (580 Mb), and Acacia confusa (566 Mb), and the two pioneer and hedgerow trees Leucaena leucocephala (1338 Mb) and Mimosa bimucronata (641 Mb). Phylogenetic inference shows that the mimosoid clade has a much higher evolutionary rate than the other clades of Caesalpinioideae. Macrosynteny comparison suggests that the fusion and breakage of an unstable chromosome are responsible for the difference in basic chromosome number (13 or 14) for Caesalpinioideae. After an ancient whole-genome duplication (WGD) shared by all Caesalpinioideae species (CWGD, ∼72.0 million years ago [MYA]), there were two recent successive WGD events, LWGD-1 (16.2-19.5 MYA) and LWGD-2 (7.1-9.5 MYA), in L. leucocephala. Thereafter, ∼40% gene loss and genome-size contraction have occurred during the diploidization process in L. leucocephala. To investigate secondary metabolites, we identified all gene copies involved in mimosine metabolism in these species and found that the abundance of mimosine biosynthesis genes in L. leucocephala largely explains its high mimosine production. We also identified the set of all potential genes involved in triterpenoid saponin biosynthesis in G. sinensis, which is more complete than that based on previous transcriptome-derived unigenes. Our results and genomic resources will facilitate biological studies of Caesalpinioideae and promote the utilization of valuable secondary metabolites.

A chromosome-level genome assembly of the disco clam, Ctenoides ales

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales (C. ales), which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and reported a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein-coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, is a valuable resource for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

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.

Chromosome-scale Elaeis guineensis and E. oleifera assemblies: comparative genomics of oil palm and other Arecaceae

Elaeis guineensis and E. oleifera are the two species of oil palm. E. guineensis is the most widely cultivated commercial species, and introgression of desirable traits from E. oleifera is ongoing. We report an improved E. guineensis genome assembly with substantially increased continuity and completeness, as well as the first chromosome-scale E. oleifera genome assembly. Each assembly was obtained by integration of long-read sequencing, proximity ligation sequencing, optical mapping, and genetic mapping. High interspecific genome conservation is observed between the two species. The study provides the most extensive gene annotation to date, including 46,697 E. guineensis and 38,658 E. oleifera gene predictions. Analyses of repetitive element families further resolve the DNA repeat architecture of both genomes. Comparative genomic analyses identified experimentally validated small structural variants between the oil palm species and resolved the mechanism of chromosomal fusions responsible for the evolutionary descending dysploidy from 18 to 16 chromosomes.

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.

A near-complete chromosome-level genome assembly of looseleaf lettuce (Lactuca sativa var. crispa)

Lettuce (Lactuca sativa L., Asteraceae) is one of the most important vegetable crops, known for its various horticultural types and significant morphological variation. The first reference genome of lettuce, a crisphead type (L. sativa var. capitata cv. Salinas), was previously released. Here, we reported a near-complete chromosome-level reference genome for looseleaf lettuce (L. sativa var. crispa). PacBio high-fidelity sequencing, Oxford Nanopore, and Hi-C technologies were employed to produce genome assembly. The final assembly is 2.59 Gb in length with a contig N50 of 205.47 Mb, anchored onto nine chromosomes, containing 14 recognizable telomeres and only 11 gaps. Repetitive sequences account for 77.11% of the genome, and 41,375 protein-coding genes were predicted, with 99.10% of these assigned functional annotations. This chromosome-level genome enriched genomic resources for various horticultural types of lettuce and will facilitate the characterization of morphological variation and genetic improvement in lettuce.