SubPhaser: a robust allopolyploid subgenome phasing method based on subgenome-specific k-mers

The complete genome assembly of Nicotiana benthamiana reveals the genetic and epigenetic landscape of centromeres

Nicotiana benthamiana is a model organism widely adopted in plant biology. Its complete assembly remains unavailable despite several recent improvements. To further improve its usefulness, we generate and phase the complete 2.85 Gb genome assembly of allotetraploid N. benthamiana. We find that although Solanaceae centromeres are widely dominated by Ty3/Gypsy retrotransposons, satellite-based centromeres are surprisingly common in N. benthamiana, with 11 of 19 centromeres featured by megabase-scale satellite arrays. Interestingly, the satellite-enriched and satellite-free centromeres are extensively invaded by distinct Gypsy retrotransposons which CENH3 protein more preferentially occupies, suggestive of their crucial roles in centromere function. We demonstrate that ribosomal DNA is a major origin of centromeric satellites, and mitochondrial DNA could be employed as a core component of the centromere. Subgenome analysis indicates that the emergence of satellite arrays probably drives new centromere formation. Altogether, we propose that N. benthamiana centromeres evolved via neocentromere formation, satellite expansion, retrotransposon enrichment and mtDNA integration.

The near-complete genome assembly of hexaploid wild oat reveals its genome evolution and divergence with cultivated oats

Avena sterilis, the ancestral species of cultivated oats, is a valuable genetic resource for oat improvement. Here we generated a near-complete 10.99 Gb A. sterilis genome and a high-quality 10.89 Gb cultivated oat genome. Genome evolution analysis revealed the centromeres dynamic and structural variations landscape associated with domestication between wild and cultivated oats. Population genetic analysis of 117 wild and cultivated oat accessions worldwide detected many candidate genes associated with important agronomic traits for oat domestication and improvement. Remarkably, a large fragment duplication from chromosomes 4A to 4D harbouring many agronomically important genes was detected during oat domestication and was fixed in almost all cultivated oats from around the world. The genes in the duplication region from 4A showed significantly higher expression levels and lower methylation levels than the orthologous genes located on 4D in A. sterilis. This study provides valuable resources for evolutionary and functional genomics and genetic improvement of oat.

The genome of Eleocharis vivipara elucidates the genetics of C3-C4 photosynthetic plasticity and karyotype evolution in the Cyperaceae

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C3 photosynthesis underwater and C4 photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi-C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C4 Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C4 pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C4 photosynthesis in E. vivipara, and hypothesize that the evolution of C4 photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C3 and C4 photosynthesis, offering promising avenues for crop improvement and breeding.

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.

Alternative splicing dynamically regulates common carp embryogenesis under thermal stress

BACKGROUND

Thermal stress is a major environmental factor affecting fish development and survival. Common carp (Cyprinus carpio) are susceptible to heat stress in their embryonic and larval phases, but the thermal stress response of alternative splicing during common carp embryogenesis remains poorly understood.

RESULTS

Using RNA-seq data from eight developmental stages and four temperatures, we constructed a comprehensive profile of alternative splicing (AS) during the embryogenesis of common carp, and found that AS genes and events are widely distributed among all stages. A total of 5,835 developmental stage-specific AS (SAS) genes, 21,368 temperature-specific differentially expressed genes (TDEGs), and 2,652 temperature-specific differentially AS (TDAS) genes were identified. Hub TDAS genes in each developmental stage, such as taf2, hnrnpa1, and drg2, were identified through protein-protein interaction (PPI) network analysis. The early developmental stages may be more sensitive to temperature, with thermal stress leading to a massive increase in the number of expressed transcripts, TDEGs, and TDAS genes in the morula stage, followed by the gastrula stage. GO and KEGG analyses showed that from the morula stage to the neurula stage, TDAS genes were more involved in intracellular transport, protein modification, and localization processes, while from the optic vesicle stage to one day post-hatching, they participated more in biosynthetic processes. Further subgenomic analysis revealed that the number of AS genes and events in subgenome B was generally higher than that in subgenome A, and the homologous AS genes were significantly enriched in basic life activity pathways, such as mTOR signaling pathway, p53 signaling pathway, and MAPK signaling pathway. Additionally, lncRNAs can play a regulatory role in the response to thermal stress by targeting AS genes such as lmnl3, affecting biological processes such as apoptosis and axon guidance.

CONCLUSIONS

In short, thermal stress can affect alternative splicing regulation during common carp embryogenesis at multiple levels. Our work complemented some gaps in the study of alternative splicing at both levels of embryogenesis and thermal stress in C. carpio and contributed to the comprehension of environmental adaptation formation in polyploid fishes during embryogenesis.

Novel centromeric repetitive DNA elements reveal karyotype dynamics in polyploid sainfoin (Onobrychis viciifolia)

Polyploidy is a common feature in eukaryotes with one of paramount consequences leading to better environmental adaptation. Heterochromatin is often located at telomeres and centromeres and contains repetitive DNA sequences. Sainfoin (Onobrychis viciifolia) is an important perennial forage legume for sustainable agriculture. However, there are only a few studies on the sainfoin genome and chromosomes. In this study, novel tandem repetitive DNA sequences of the sainfoin genome (OnVi180, OnVi169, OnVi176 and OnVidimer) were characterized using bioinformatics, molecular and cytogenetic approaches. The OnVi180 and OnVi169 elements colocalized within functional centromeres. The OnVi176 and OnVidimer elements were localized in centromeric, subtelomeric and interstitial regions. We constructed a sainfoin karyotype that distinguishes the seven basic chromosome groups. Our study provides the first detailed description of heterochromatin and chromosome structure of sainfoin and proposes an origin of heterozygous ancestral genomes, possibly from the same ancestral diploid species, not necessarily from different species, or for chromosome rearrangements after polyploidy. Overall, we discuss our novel and complementary findings in a polyploid crop with unique and complex chromosomal features.

Unravelling the novel sex determination genotype with 'ZY' and a distinctive 2.15-2.95 Mb inversion among poplar species through haplotype-resolved genome assembly and comparative genomics analysis

Populus tomentosa, an indigenous tree species, is widely distributed and cultivated over 1,000,000 km2 in China, contributing significantly to forest production, ecological conservation and urban-rural greening. Although a reference genome is available for P. tomentosa, the intricate interspecific hybrid origins, chromosome structural variations (SVs) and sex determination mechanisms remain confusion and unclear due to its broad and even overlapping geographical distribution, extensive morphological variations and cross infiltration among white poplar species. We conducted a haplotype-resolved de novo assembly of P. tomentosa elite individual GM107, which comprises subgenomes a and b with a total genome size of 714.9 Mb. We then analysed the formation of hybrid species and the phylogenetic evolution and sex differentiation across the entire genus. Phylogenomic analyses suggested that GM107 likely originated from a hybridisation event between P. alba (♀) and P. davidiana (♂) which diverged at approximately 3.8 Mya. A total of 1551 chromosome SVs were identified between the two subgenomes. More noteworthily, a distinctive inversion structure spanning 2.15-2.95 Mb was unveiled among Populus, Tacamahaca, Turaga, Aigeiros poplar species and Salix, highlighting a unique evolutionary feature. Intriguingly, a novel sex genotype of the ZY type, which represents a crossover between XY and ZW systems, was identified and confirmed through both natural and artificial hybrids populations. These novel insights offer significant theoretical value for the study of the species' evolutionary origins and serve as a valuable resource for ecological genetics and forest biotechnology.

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

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

Genome architecture of the allotetraploid wild grass Aegilops ventricosa reveals its evolutionary history and contributions to wheat improvement

The allotetraploid wild grass Aegilops ventricosa (2n = 4x = 28, genome DvDvNvNv) has been recognized as an important germplasm resource for wheat improvement owing to its ability to tolerate biotic stresses. In particular, the 2NvS segment from Ae. ventricosa, as a stable and effective resistance source, has contributed greatly to wheat improvement. The 2NvS/2AS translocation is a prevalent chromosomal translocation between common wheat and wild relatives, ranking just behind the 1B/1R translocation in importance for modern wheat breeding. Here, we assembled a high-quality chromosome-level reference genome of Ae. ventricosa RM271 with a total length of 8.67 Gb. Phylogenomic analyses revealed that the progenitor of the Dv subgenome of Ae. ventricosa is Ae. tauschii ssp. tauschii (genome DD); by contrast, the progenitor of the D subgenome of bread wheat (Triticum aestivum L.) is Ae. tauschii ssp. strangulata (genome DD). The oldest polyploidization time of Ae. ventricosa occurred ∼0.7 mya. The Dv subgenome of Ae. ventricosa is less conserved than the D subgenome of bread wheat. Construction of a graph-based pangenome of 2AS/6NvL (originally known as 2NvS) segments from Ae. ventricosa and other genomes in the Triticeae enabled us to identify candidate resistance genes sourced from Ae. ventricosa. We identified 12 nonredundant introgressed segments from the Dv and Nv subgenomes using a large winter wheat collection representing the full diversity of the European wheat genetic pool, and 29.40% of European wheat varieties inherit at least one of these segments. The high-quality RM271 reference genome will provide a basis for cloning key genes, including the Yr17-Lr37-Sr38-Cre5 resistance gene cluster in Ae. ventricosa, and facilitate the full use of elite wild genetic resources to accelerate wheat improvement.

Chromosome-level genome assemblies reveal genome evolution of an invasive plant Phragmites australis

Biological invasions pose a significant threat to ecosystems, disrupting local biodiversity and ecosystem functions. The genomic underpinnings of invasiveness, however, are still largely unknown, making it difficult to predict and manage invasive species effectively. The common reed (Phragmites australis) is a dominant grass species in wetland ecosystems and has become particularly invasive when transferred from Europe to North America. Here, we present a high-quality gap-free, telomere-to-telomere genome assembly of Phragmites australis consisting of 24 pseudochromosomes and a B chromosome. Fully phased subgenomes demonstrated considerable subgenome dominance and revealed the divergence of diploid progenitors approximately 30.9 million years ago. Comparative genomics using chromosome-level scaffolds for three other lineages and a previously published draft genome assembly of an invasive lineage revealed that gene family expansions in the form of tandem duplications may have contributed to the invasiveness of the lineage. This study sheds light on the genome evolution of Arundinoideae grasses and suggests that genetic drivers, such as gene family expansions and tandem duplications, may underly the processes of biological invasion in plants. These findings provide a crucial step toward understanding and managing the genetic basis of invasiveness in plant species.

The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P. setigerum (2n = 44), P. somniferum (2n = 22), P. rhoeas (2n = 14), and P. bracteatum (2n = 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P. somniferum and P. setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

Polyploids of Brassicaceae: Genomic Insights and Assembly Strategies

The Brassicaceae family is distinguished by its inclusion of high-value crops such as cabbage, broccoli, mustard, and wasabi, all noted for their glucosinolates. In this family, many polyploidy species are distributed and shaped by numerous whole-genome duplications, independent genome doublings, and hybridization events. The evolutionary trajectory of the family is marked by enhanced diversification and lineage splitting after paleo- and meso-polyploidization, with discernible remnants of whole-genome duplications within their genomes. The recent neopolyploidization events notably increased the proportion of polyploid species within the family. Although sequencing efforts for the Brassicaceae genome have been robust, accurately distinguishing sub-genomes remains a significant challenge, frequently complicating the assembly process. Assembly strategies include comparative analyses with ancestral species and examining k-mers, long terminal repeat retrotransposons, and pollen sequencing. This review comprehensively explores the unique genomic characteristics of the Brassicaceae family, with a particular emphasis on polyploidization events and the latest strategies for sequencing and assembly. This review will significantly improve our understanding of polyploidy in the Brassicaceae family and assist in future genome assembly methods.

Chromosome-scale assembly of the wild cereal relative Elymus sibiricus

Elymus species, belonging to Triticeae tribe, is a tertiary gene pool for improvement of major cereal crops. Elymus sibiricus, a tetraploid with StH genome, is a typical species in the genus Elymus, which is widely utilized as a high-quality perennial forage grass in template regions. In this study, we report the construction of a chromosome-scale reference assembly of E. sibiricus line Gaomu No. 1 based on PacBio HiFi reads and chromosome conformation capture. Subgenome St and H were well phased by assisting with kmer and subgenome-specific repetitive sequence. The total assembly size was 6.929 Gb with a contig N50 of 49.518 Mb. In total, 89,800 protein-coding genes were predicted. The repetitive sequences accounted for 82.49% of the genome in E. sibiricus. Comparative genome analysis confirmed a major species-specific 4H/6H reciprocal translocation in E. sibiricus. The E. sibiricus assembly will be much helpful to exploit genetic resource of StH species in genus Elymus, and provides an important tool for E. sibiricus domestication.

Homoeologs in Allopolyploids: Navigating Redundancy as Both an Evolutionary Opportunity and a Technical Challenge-A Transcriptomics Perspective

Allopolyploidy in plants involves the merging of two or more distinct parental genomes into a single nucleus, a significant evolutionary process in the plant kingdom. Transcriptomic analysis provides invaluable insights into allopolyploid plants by elucidating the fate of duplicated genes, revealing evolutionary novelties and uncovering their environmental adaptations. By examining gene expression profiles, scientists can discern how duplicated genes have evolved to acquire new functions or regulatory roles. This process often leads to the development of novel traits and adaptive strategies that allopolyploid plants leverage to thrive in diverse ecological niches. Understanding these molecular mechanisms not only enhances our appreciation of the genetic complexity underlying allopolyploidy but also underscores their importance in agriculture and ecosystem resilience. However, transcriptome profiling is challenging due to genomic redundancy, which is further complicated by the presence of multiple chromosomes sets and the variations among homoeologs and allelic genes. Prior to transcriptome analysis, sub-genome phasing and homoeology inference are essential for obtaining a comprehensive view of gene expression. This review aims to clarify the terminology in this field, identify the most challenging aspects of transcriptome analysis, explain their inherent difficulties, and suggest reliable analytic strategies. Furthermore, bulk RNA-seq is highlighted as a primary method for studying allopolyploid gene expression, focusing on critical steps like read mapping and normalization in differential gene expression analysis. This approach effectively captures gene expression from both parental genomes, facilitating a comprehensive analysis of their combined profiles. Its sensitivity in detecting low-abundance transcripts allows for subtle differences between parental genomes to be identified, crucial for understanding regulatory dynamics and gene expression balance in allopolyploids.

Hijacking a rapid and scalable metagenomic method reveals subgenome dynamics and evolution in polyploid plants

Premise

The genomes of polyploid plants archive the evolutionary events leading to their present forms. However, plant polyploid genomes present numerous hurdles to the genome comparison algorithms for classification of polyploid types and exploring genome dynamics.

Methods

Here, the problem of intra- and inter-genome comparison for examining polyploid genomes is reframed as a metagenomic problem, enabling the use of the rapid and scalable MinHashing approach. To determine how types of polyploidy are described by this metagenomic approach, plant genomes were examined from across the polyploid spectrum for both k-mer composition and frequency with a range of k-mer sizes. In this approach, no subgenome-specific k-mers are identified; rather, whole-chromosome k-mer subspaces were utilized.

Results

Given chromosome-scale genome assemblies with sufficient subgenome-specific repetitive element content, literature-verified subgenomic and genomic evolutionary relationships were revealed, including distinguishing auto- from allopolyploidy and putative progenitor genome assignment. The sequences responsible were the rapidly evolving landscape of transposable elements. An investigation into the MinHashing parameters revealed that the downsampled k-mer space (genomic signatures) produced excellent approximations of sequence similarity. Furthermore, the clustering approach used for comparison of the genomic signatures is scrutinized to ensure applicability of the metagenomics-based method.

Discussion

The easily implementable and highly computationally efficient MinHashing-based sequence comparison strategy enables comparative subgenomics and genomics for large and complex polyploid plant genomes. Such comparisons provide evidence for polyploidy-type subgenomic assignments. In cases where subgenome-specific repeat signal may not be adequate given a chromosomes' global k-mer profile, alternative methods that are more specific but more computationally complex outperform this approach.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

Chromosome-level assembly of the synthetic hexaploid wheat-derived cultivar Chuanmai 104

Synthetic hexaploid wheats (SHWs) are effective genetic resources for transferring agronomically important genes from wild relatives to common wheat (Triticum aestivum L.). Dozens of reference-quality pseudomolecule assemblies of hexaploid wheat have been generated, but none is reported for SHW-derived cultivars. Here, we generated a chromosome-scale assembly for the SHW-derived cultivar 'Chuanmai 104' based on PacBio HiFi reads and chromosome conformation capture sequencing. The total assembly size was 14.81 Gb with a contig N50 length of 58.25 Mb. A BUSCO analysis yielded a completeness score of 99.30%. In total, repetitive elements comprised 81.36% of the genome and 122,554 high-confidence protein-coding gene models were predicted. In summary, the first chromosome-level assembly for a SHW-derived cultivar presents a promising outlook for the study and utilization of SHWs in wheat improvement, which is essential to meet the global food demand.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B. racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77 266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

Pangenome analysis reveals transposon-driven genome evolution in cotton

BACKGROUND

Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution.

RESULTS

Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes.

CONCLUSIONS

Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

Chromosome-level genome of putative autohexaploid Actinidia deliciosa provides insights into polyploidisation and evolution

Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

The genomes of Dahlia pinnata, Cosmos bipinnatus, and Bidens alba in tribe Coreopsideae provide insights into polyploid evolution and inulin biosynthesis

BACKGROUND

The Coreopsideae tribe, a subset of the Asteraceae family, encompasses economically vital genera like Dahlia, Cosmos, and Bidens, which are widely employed in medicine, horticulture, ecology, and food applications. Nevertheless, the lack of reference genomes hinders evolutionary and biological investigations in this tribe.

RESULTS

Here, we present 3 haplotype-resolved chromosome-level reference genomes of the tribe Coreopsideae, including 2 popular flowering plants (Dahlia pinnata and Cosmos bipinnatus) and 1 invasive weed plant (Bidens alba), with assembled genome sizes 3.93 G, 1.02 G, and 1.87 G, respectively. We found that Gypsy transposable elements contribute mostly to the larger genome size of D. pinnata, and multiple chromosome rearrangements have occurred in tribe Coreopsideae. Besides the shared whole-genome duplication (WGD-2) in the Heliantheae alliance, our analyses showed that D. pinnata and B. alba each underwent an independent recent WGD-3 event: in D. pinnata, it is more likely to be a self-WGD, while in B. alba, it is from the hybridization of 2 ancestor species. Further, we identified key genes in the inulin metabolic pathway and found that the pseudogenization of 1-FEH1 and 1-FEH2 genes in D. pinnata and the deletion of 3 key residues of 1-FFT proteins in C. bipinnatus and B. alba may probably explain why D. pinnata produces much more inulin than the other 2 plants.

CONCLUSIONS

Collectively, the genomic resources for the Coreopsideae tribe will promote phylogenomics in Asteraceae plants, facilitate ornamental molecular breeding improvements and inulin production, and help prevent invasive weeds.

A reference genome of Commelinales provides insights into the commelinids evolution and global spread of water hyacinth (Pontederia crassipes)

Commelinales belongs to the commelinids clade, which also comprises Poales that includes the most important monocot species, such as rice, wheat, and maize. No reference genome of Commelinales is currently available. Water hyacinth (Pontederia crassipes or Eichhornia crassipes), a member of Commelinales, is one of the devastating aquatic weeds, although it is also grown as an ornamental and medical plant. Here, we present a chromosome-scale reference genome of the tetraploid water hyacinth with a total length of 1.22 Gb (over 95% of the estimated size) across 8 pseudochromosome pairs. With the representative genomes, we reconstructed a phylogeny of the commelinids, which supported Zingiberales and Commelinales being sister lineages of Arecales and shed lights on the controversial relationship of the orders. We also reconstructed ancestral karyotypes of the commelinids clade and confirmed the ancient commelinids genome having 8 chromosomes but not 5 as previously reported. Gene family analysis revealed contraction of disease-resistance genes during polyploidization of water hyacinth, likely a result of fitness requirement for its role as a weed. Genetic diversity analysis using 9 water hyacinth lines from 3 continents (South America, Asia, and Europe) revealed very closely related nuclear genomes and almost identical chloroplast genomes of the materials, as well as provided clues about the global dispersal of water hyacinth. The genomic resources of P. crassipes reported here contribute a crucial missing link of the commelinids species and offer novel insights into their phylogeny.

Telomere-to-telomere genome of the allotetraploid legume Sesbania cannabina reveals transposon-driven subgenome divergence and mechanisms of alkaline stress tolerance

Alkaline soils pose an increasing problem for agriculture worldwide, but using stress-tolerant plants as green manure can improve marginal land. Here, we show that the legume Sesbania cannabina is very tolerant to alkaline conditions and, when used as a green manure, substantially improves alkaline soil. To understand genome evolution and the mechanisms of stress tolerance in this allotetraploid legume, we generated the first telomere-to-telomere genome assembly of S. cannabina spanning ∼2,087 Mb. The assembly included all centromeric regions, which contain centromeric satellite repeats, and complete chromosome ends with telomeric characteristics. Further genome analysis distinguished A and B subgenomes, which diverged approximately 7.9 million years ago. Comparative genomic analysis revealed that the chromosome homoeologs underwent large-scale inversion events (>10 Mb) and a significant, transposon-driven size expansion of the chromosome 5A homoeolog. We further identified four specific alkali-induced phosphate transporter genes in S. cannabina; these may function in alkali tolerance by relieving the deficiency in available phosphorus in alkaline soil. Our work highlights the significance of S. cannabina as a green tool to improve marginal lands and sheds light on subgenome evolution and adaptation to alkaline soils.

Phased gap-free genome assembly of octoploid cultivated strawberry illustrates the genetic and epigenetic divergence among subgenomes

The genetic and epigenetic mechanisms underlying the coexistence and coordination of the four diverged subgenomes (ABCD) in octoploid strawberries (Fragaria × ananassa) remains poorly understood. In this study, we have assembled a haplotype-phased gap-free octoploid genome for the strawberry, which allowed us to uncover the sequence, structure, and epigenetic divergences among the subgenomes. The diploid progenitors of the octoploid strawberry, apart from subgenome A (Fragaria vesca), have been a subject of public controversy. Phylogenomic analyses revealed a close relationship between diploid species Fragaria iinumae and subgenomes B, C, and D. Subgenome A, closely related to F. vesca, retains the highest number of genes, exhibits the lowest content of transposable elements (TEs), experiences the strongest purifying selection, shows the lowest DNA methylation levels, and displays the highest expression level compared to the other three subgenomes. Transcriptome and DNA methylome analyses revealed that subgenome A-biased genes were enriched in fruit development biological processes. In contrast, although subgenomes B, C, and D contain equivalent amounts of repetitive sequences, they exhibit diverged methylation levels, particularly for TEs located near genes. Taken together, our findings provide valuable insights into the evolutionary patterns of subgenome structure, divergence and epigenetic dynamics in octoploid strawberries, which could be utilized in strawberry genetics and breeding research.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

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

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Telomere-to-telomere and haplotype-resolved genome assembly of the Chinese cork oak (Quercus variabilis)

The Quercus variabilis, a deciduous broadleaved tree species, holds significant ecological and economical value. While a chromosome-level genome for this species has been made available, it remains riddled with unanchored sequences and gaps. In this study, we present a nearly complete comprehensive telomere-to-telomere (T2T) and haplotype-resolved reference genome for Q. variabilis. This was achieved through the integration of ONT ultra-long reads, PacBio HiFi long reads, and Hi-C data. The resultant two haplotype genomes measure 789 Mb and 768 Mb in length, with a contig N50 of 65 Mb and 56 Mb, and were anchored to 12 allelic chromosomes. Within this T2T haplotype-resolved assembly, we predicted 36,830 and 36,370 protein-coding genes, with 95.9% and 96.0% functional annotation for each haplotype genome. The availability of the T2T and haplotype-resolved reference genome lays a solid foundation, not only for illustrating genome structure and functional genomics studies but also to inform and facilitate genetic breeding and improvement of cultivated Quercus species.

Genome evolution and initial breeding of the Triticeae grass Leymus chinensis dominating the Eurasian Steppe

Leymus chinensis, a dominant perennial grass in the Eurasian Steppe, is well known for its remarkable adaptability and forage quality. Hardly any breeding has been done on the grass, limiting its potential in ecological restoration and forage productivity. To enable genetic improvement of the untapped, important species, we obtained a 7.85-Gb high-quality genome of L. chinensis with a particularly long contig N50 (318.49 Mb). Its allotetraploid genome is estimated to originate 5.29 million years ago (MYA) from a cross between the Ns-subgenome relating to Psathyrostachys and the unknown Xm-subgenome. Multiple bursts of transposons during 0.433-1.842 MYA after genome allopolyploidization, which involved predominantly the Tekay and Angela of LTR retrotransposons, contributed to its genome expansion and complexity. With the genome resource available, we successfully developed a genetic transformation system as well as the gene-editing pipeline in L. chinensis. We knocked out the monocot-specific miR528 using CRISPR/Cas9, resulting in the improvement of yield-related traits with increases in the tiller number and growth rate. Our research provides valuable genomic resources for Triticeae evolutionary studies and presents a conceptual framework illustrating the utilization of genomic information and genome editing to accelerate the improvement of wild L. chinensis with features such as polyploidization and self-incompatibility.

A high-quality chromosome-level wild rice genome of Oryza coarctata

Oryza coarctata (2n = 4X = 48, KKLL) is an allotetraploid, undomesticated relative of rice and the only species in the genus Oryza with tolerance to high salinity and submergence. Therefore, it contains important stress and tolerance genes/factors for rice. The initial draft genome published was limited by data and technical restrictions, leading to an incomplete and highly fragmented assembly. This study reports a new, highly contiguous chromosome-level genome assembly and annotation of O. coarctata. PacBio high-quality HiFi reads generated 460 contigs with a total length of 573.4 Mb and an N50 of 23.1 Mb, which were assembled into scaffolds with Hi-C data, anchoring 96.99% of the assembly onto 24 chromosomes. The genome assembly comprises 45,571 genes, and repetitive content contributes 25.5% of the genome. This study provides the novel identification of the KK and LL genome types of the genus Oryza, leading to valuable insights into rice genome evolution. The chromosome-level genome assembly of O. coarctata is a valuable resource for rice research and molecular breeding.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

Haplotype-resolved genomes of wild octoploid progenitors illuminate genomic diversifications from wild relatives to cultivated strawberry

Strawberry is an emerging model for studying polyploid genome evolution and rapid domestication of fruit crops. Here we report haplotype-resolved genomes of two wild octoploids (Fragaria chiloensis and Fragaria virginiana), the progenitor species of cultivated strawberry. Substantial variation is identified between species and between haplotypes. We redefine the four subgenomes and track the genetic contributions of diploid species by additional sequencing of the diploid F. nipponica genome. We provide multiple lines of evidence that F. vesca and F. iinumae, rather than other described extant species, are the closest living relatives of these wild and cultivated octoploids. In response to coexistence with quadruplicate gene copies, the octoploid strawberries have experienced subgenome dominance, homoeologous exchanges and coordinated expression of homoeologous genes. However, some homoeologues have substantially altered expression bias after speciation and during domestication. These findings enhance our understanding of the origin, genome evolution and domestication of strawberries.

The genome of okra (Abelmoschus esculentus) provides insights into its genome evolution and high nutrient content

Okra (Abelmoschus esculentus) is an important vegetable crop with high nutritional value. However, the mechanism underlying its high nutrient content remains poorly understood. Here, we present a chromosome-scale genome of okra with a size of 1.19 Gb. Comparative genomics analysis revealed the phylogenetic status of A. esculentus, as well as whole-genome duplication (WGD) events that have occurred widely across the Malvaceae species. We found that okra has experienced three additional WGDs compared with the diploid cotton Gossypium raimondii, resulting in a large chromosome number (2n = 130). After three WGDs, okra has undergone extensive genomic deletions and retained substantial numbers of genes related to secondary metabolite biosynthesis and environmental adaptation, resulting in significant differences between okra and G. raimondii in the gene families related to cellulose synthesis. Combining transcriptomic and metabolomic analysis, we revealed the relationship between gene expression and metabolite content change across different okra developmental stages. Furthermore, the sinapic acid/S-lignin biosynthesis-related gene families have experienced remarkable expansion in okra, and the expression of key enzymes involved in the sinapic acid/S-lignin biosynthesis pathway vary greatly across developmental periods, which partially explains the differences in metabolite content across the different stages. Our study gains insights into the comprehensive evolutionary history of Malvaceae species and the genetic basis that underlies the nutrient content changes in okra, which will facilitate the functional study and genetic improvement of okra varieties.

The allotetraploid horseradish genome provides insights into subgenome diversification and formation of critical traits

Polyploidization can provide a wealth of genetic variation for adaptive evolution and speciation, but understanding the mechanisms of subgenome evolution as well as its dynamics and ultimate consequences remains elusive. Here, we report the telomere-to-telomere (T2T) gap-free reference genome of allotetraploid horseradish (Armoracia rusticana) sequenced using a comprehensive strategy. The (epi)genomic architecture and 3D chromatin structure of the A and B subgenomes differ significantly, suggesting that both the dynamics of the dominant long terminal repeat retrotransposons and DNA methylation have played critical roles in subgenome diversification. Investigation of the genetic basis of biosynthesis of glucosinolates (GSLs) and horseradish peroxidases reveals both the important role of polyploidization and subgenome differentiation in shaping the key traits. Continuous duplication and divergence of essential genes of GSL biosynthesis (e.g., FMO_GS-OX, IGMT, and GH1 gene family) contribute to the broad GSL profile in horseradish. Overall, the T2T assembly of the allotetraploid horseradish genome expands our understanding of polyploid genome evolution and provides a fundamental genetic resource for breeding and genetic improvement of horseradish.

Differential subgenome expression underlies biomass accumulation in allotetraploid Pennisetum giganteum

BACKGROUND

Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation.

RESULTS

Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum.

CONCLUSIONS

Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes

Sour cherry (Prunus cerasus L.) is a valuable fruit crop in the Rosaceae family and a hybrid between progenitors closely related to extant Prunus fruticosa (ground cherry) and Prunus avium (sweet cherry). Here we report a chromosome-scale genome assembly for sour cherry cultivar Montmorency, the predominant cultivar grown in the USA. We also generated a draft assembly of P. fruticosa to use alongside a published P. avium sequence for syntelog-based subgenome assignments for 'Montmorency' and provide compelling evidence P. fruticosa is also an allotetraploid. Using hierarchal k-mer clustering and phylogenomics, we show 'Montmorency' is trigenomic, containing two distinct subgenomes inherited from a P. fruticosa-like ancestor (A and A') and two copies of the same subgenome inherited from a P. avium-like ancestor (BB). The genome composition of 'Montmorency' is AA'BB and little-to-no recombination has occurred between progenitor subgenomes (A/A' and B). In Prunus, two known classes of genes are important to breeding strategies: the self-incompatibility loci (S-alleles), which determine compatible crosses, successful fertilization, and fruit set, and the Dormancy Associated MADS-box genes (DAMs), which strongly affect dormancy transitions and flowering time. The S-alleles and DAMs in 'Montmorency' and P. fruticosa were manually annotated and support subgenome assignments. Lastly, the hybridization event 'Montmorency' is descended from was estimated to have occurred less than 1.61 million years ago, making sour cherry a relatively recent allotetraploid. The 'Montmorency' genome highlights the evolutionary complexity of the genus Prunus and will inform future breeding strategies for sour cherry, comparative genomics in the Rosaceae, and questions regarding neopolyploidy.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.

High-quality haplotype-resolved genome assembly of cultivated octoploid strawberry

Cultivated strawberry (Fragaria × ananassa), a perennial herb belonging to the family Rosaceae, is a complex octoploid with high heterozygosity at most loci. However, there is no research on the haplotype of the octoploid strawberry genome. Here we aimed to obtain a high-quality genome of the cultivated strawberry cultivar, "Yanli", using single molecule real-time sequencing and high-throughput chromosome conformation capture technology. The "Yanli" genome was 823 Mb in size, with a long terminal repeat assembly index of 14.99. The genome was phased into two haplotypes, Hap1 (825 Mb with contig N50 of 26.70 Mb) and Hap2 (808 Mb with contig N50 of 27.51 Mb). Using the combination of Hap1 and Hap2, we obtained for the first time a haplotype-resolved genome with 56 chromosomes for the cultivated octoploid strawberry. We identified a ~ 10 Mb inversion and translocation on chromosome 2-1. 104 957 and 102 356 protein-coding genes were annotated in Hap1 and Hap2, respectively. Analysis of the genes related to the anthocyanin biosynthesis pathway revealed the structural diversity and complexity in the expression of the alleles in the octoploid F. × ananassa genome. In summary, we obtained a high-quality haplotype-resolved genome assembly of F. × ananassa, which will provide the foundation for investigating gene function and evolution of the genome of cultivated octoploid strawberry.