ksrates: positioning whole-genome duplications relative to speciation events in KS distributions

Telomere-to-telomere genome assembly and 3D chromatin architecture of Centella asiatica insight into evolution and genetic basis of triterpenoid saponin biosynthesis

Centella asiatica is renowned for its medicinal properties, particularly due to its triterpenoid saponins, such as asiaticoside and madecassoside, which are in excess demand for the cosmetic industry. However, comprehensive genomic resources for this species are lacking, which impedes the understanding of its biosynthetic pathways. Here, we report a telomere-to-telomere (T2T) C. asiatica genome. The genome size is 438.12 Mb with a contig N50 length of 54.12 Mb. The genome comprises 258.87 Mb of repetitive sequences and 25 200 protein-coding genes. Comparative genomic analyses revealed C. asiatica as an early-diverging genus within the Apiaceae family with a single whole-genome duplication (WGD, Apiaceae-ω) event following the ancient γ-triplication, contrasting with Apiaceae species that exhibit two WGD events (Apiaceae-α and Apiaceae-ω). We further constructed 3D chromatin structures, A/B compartments, and topologically associated domains (TADs) in C. asiatica leaves, elucidating the influence of chromatin organization on expression WGD-derived genes. Additionally, gene family and functional characterization analysis highlight the key role of CasiOSC03 in α-amyrin production while also revealing significant expansion and high expression of CYP716, CYP714, and UGT73 families involved in asiaticoside biosynthesis compared to other Apiaceae species. Notably, a unique and large UGT73 gene cluster, located within the same TAD, is potentially pivotal for enhancing triterpenoid saponin. Weighted gene coexpression network analysis (WGCNA) further highlighted the pathways modulated in response to methyl jasmonate (MeJA), offering insights into the regulatory networks governing saponin biosynthesis. This work not only provides a valuable genomic resource for C. asiatica but also sheds light on the molecular mechanisms driving the biosynthesis of pharmacologically important metabolites.

Comparative transcriptomics in ferns reveals key innovations and divergent evolution of the secondary cell walls

Ferns are essential for understanding plant evolution; however, their large and intricate genomes have kept their genetic landscape largely unexplored, with only a few genomes sequenced and limited transcriptomic data available. To bridge this gap, we generated extensive RNA-sequencing data across various organs from 22 representative fern species, resulting in high-quality transcriptome assemblies. These data enabled us to construct a time-calibrated phylogeny for ferns, encompassing all major clades, which revealed numerous instances of whole-genome duplication. We highlighted the distinctiveness of fern genetics, discovering that half of the identified gene families are unique to ferns. Our exploration of fern cell walls through biochemical and immunological analyses uncovered the presence of the lignin syringyl unit, along with evidence of its independent evolution in ferns. Additionally, the identification of an unusual sugar in fern cell walls suggests a divergent evolutionary trajectory in cell wall biochemistry, probably influenced by gene duplication and sub-functionalization. To facilitate further research, we have developed an online database that includes preloaded genomic and transcriptomic data for ferns and other land plants. We used this database to demonstrate the independent evolution of lignocellulosic gene modules in ferns. Our findings provide a comprehensive framework illustrating the unique evolutionary journey ferns have undertaken since diverging from the last common ancestor of euphyllophytes more than 360 million years ago.

Resolving floral development dynamics using genome and single-cell temporal transcriptome of Dendrobium devonianum

Dendrobium devonianum, a species of the Orchidaceae family, is notable for its unique floral characteristics, which include two yellow spots and purple tips on its labellum, as well as fringed edges. However, the molecular mechanisms underlying flower pattern formation in D. devonianum remain poorly understood, hindering advancements in its breeding process. Here, a chromosome-scale genome of D. devonianum was presented for the first time, revealing two significant polyploidization events. Additionally, a high-resolution single-cell transcriptomic atlas was constructed, capturing 11 distinct cell clusters. Expression patterns of MADS-box genes were identified through temporal and spatial bulk RNA-Seq, revealing alignment with the ABCDE model of flower formation. Meanwhile, mass spectrometry imaging and scRNA analyses showed that the yellow spots were primarily associated with carotenoid biosynthesis gene expression, while the purple colour is predominantly linked to anthocyanin biosynthesis gene expression. These genes were mainly expressed in the epidermis and vascular cells. Developmental trajectory analyses of epidermal cells further uncovered a gene regulatory network and several transcription factors likely responsible for fringes formation along the labellum margin. This study provides valuable insights into the molecular mechanisms driving floral colour differentiation and structural traits in D. devonianum, contributing to a deeper understanding of orchid evolution, diversification and breeding.

SOI: robust identification of orthologous synteny with the Orthology Index and broad applications in evolutionary genomics

With the explosive growth of whole-genome datasets, accurate detection of orthologous synteny has become crucial for reconstructing evolutionary history. However, current methods for identifying orthologous synteny face great limitations, particularly in scaling with varied polyploidy histories and accurately removing out-paralogous synteny. In this study, we developed a scalable and robust approach, based on the Orthology Index (OI), to effectively identify orthologous synteny. Our evaluation across a large-scale empirical dataset with diverse polyploidization events demonstrated the high reliability and robustness of the OI method. Simulation-based benchmarks further validated the accuracy of our method, showing its superior performance against existing methods across a wide range of scenarios. Additionally, we explored its broad applications in reconstructing the evolutionary histories of plant genomes, including the inference of polyploidy, identification of reticulation, and phylogenomics. In conclusion, OI offers a robust, interpretable, and scalable approach for identifying orthologous synteny, facilitating more accurate and efficient analyses in plant evolutionary genomics.

Chromosome-level de novo genome unveils the evolution of Gleditsia sinensis and thorns development

Gleditsia sinensis Lam. (G. sinensis) as an important species within the Leguminosae family, has been utilized in Chinese medicine for centuries, and its thorns serve as a chief medicinal ingredient. The absence of a comprehensive genome database has hindered its in-depth research. In this investigation, a chromosome-level de novo genome assembly of G. sinensis 'Yulin No.1' was achieved, which harbors a 786.13 Mb sized genome with 36,408 protein-coding genes and experiences two WGD events. The comparative and evolutionary analysis unveiled the close phylogenetic relationship between G. sinensis and eight other Leguminosae species. The WGCNA and gene family analysis further indicated that GsinMYB was involved in the development of thorns. This investigation offered a high-level genome of G. sinensis, facilitating comparisons in Leguminosae species evolution and functional elucidation. It also provided key insights for further research on the molecular regulation mechanisms of thorn development in plants and the molecular breeding of G. sinensis.

doubletrouble: an R/Bioconductor package for the identification, classification, and analysis of gene and genome duplications

SUMMARY

Gene and genome duplications are major evolutionary forces that shape the diversity and complexity of life. However, different duplication modes have distinct impacts on gene function, expression, and regulation. Existing tools for identifying and classifying duplicated genes are either outdated or not user-friendly. Here, we present doubletrouble, an R/Bioconductor package that provides a comprehensive and robust framework for analyzing duplicated genes from genomic data. doubletrouble can detect and classify gene pairs as derived from six duplication modes (segmental, tandem, proximal, retrotransposon-derived, DNA transposon-derived, and dispersed duplications), calculate substitution rates, detect signatures of putative whole-genome duplication events, and visualize results as publication-ready figures. We applied doubletrouble to classify the duplicated gene repertoire in 822 eukaryotic genomes, and results were made available through a user-friendly web interface.

AVAILABILITY AND IMPLEMENTATION

doubletrouble is available on Bioconductor (https://bioconductor.org/packages/doubletrouble), and the source code is available in a GitHub repository (https://github.com/almeidasilvaf/doubletrouble). doubletroubledb is available online at https://almeidasilvaf.github.io/doubletroubledb/.

Lonicera caerulea genome reveals molecular mechanisms of freezing tolerance and anthocyanin biosynthesis

INTRODUCTION

Lonicera caerulea L. (blue honeysuckle) is a noteworthy fleshy-fruited tree and a prominent medicinal plant, which possesses notable characteristics such as exceptional resilience to winter conditions and early maturation, and the richest source of functional anthocyanins, particularly cyanidin-3-glucoside. The molecular mechanisms responsible for its freezing tolerance and anthocyanin biosynthesis remain largely unknown.

OBJECTIVES

Here, a chromosome-scale genome of L. caerulea was presented, aiming to examine the genetic foundations that underlie these characteristics of blue honeysuckle.

METHODS

The PacBio HiFi reads and Hi-C data were used to construct high-quality genome of blue honeysuckle. Comparative genomic and transcriptomic analyses were conducted to elucidate the molecular mechanisms of freezing tolerance and anthocyanin biosynthesis.

RESULTS

Comparative genomics analysis between L. caerulea and L. japonica revealed that the dynamic changes of duplicated genes contributed to their phytochemical reconstruction and environmental adaptation. Moreover, the ABA and ICE-CBF-COR signaling pathways were closely correlated to the freezing tolerance of L. caerulea. Genome-wide identification and biochemical function indicated that three anthocyanin 3',5'-O-methyltransferases (LcOMT2, LcOMT14, and LcOMT20) and two 3'-O-glycosyltransferases (LcUGT78X1 and LcUGT95P1) were responsible for anthocyanin biosynthesis. In addition, LcUGT78X1 was regarded as the potent glycosyltransferase for the accumulation of cyanidin-3-glucoside in L. caerulea.

CONCLUSION

This research elucidates the crucial roles of the ABA and ICE-CBF-COR signaling pathways in enhancing freezing tolerance, while also identifying highly efficient anthocyanin biosynthetic enzymes in L. caerulea. These findings advance the understanding of environmental adaptation and phytochemical production in Lonicera species.

Accurate Inference of the Polyploid Continuum Using Forward-Time Simulations

Multiple rounds of whole-genome duplication (WGD) followed by diploidization have occurred throughout the evolutionary history of angiosperms. Much work has been done to model the genomic consequences and evolutionary significance of WGD. While researchers have historically modeled polyploids as either allopolyploids or autopolyploids, the variety of natural polyploids span a continuum of differentiation across multiple parameters, such as the extent of polysomic versus disomic inheritance, and the degree of genetic differentiation between the ancestral lineages. Here we present a forward-time polyploid genome evolution simulator called SpecKS. SpecKS models polyploid speciation as originating from a 2D continuum, whose dimensions account for both the level of genetic differentiation between the ancestral parental genomes, as well the time lag between ancestral speciation and their subsequent reunion in the derived polyploid. Using extensive simulations, we demonstrate that changes in initial conditions along either dimension of the 2D continuum deterministically affect the shape of the Ks histogram. Our findings indicate that the error in the common method of estimating WGD time from the Ks histogram peak scales with the degree of allopolyploidy, and we present an alternative, accurate estimation method that is independent of the degree of allopolyploidy. Lastly, we use SpecKS to derive tests that infer both the lag time between parental divergence and WGD time, and the diversity of the ancestral species, from an input Ks histogram. We apply the latter test to transcriptomic data from over 200 species across the plant kingdom, the results of which are concordant with the prevailing theory that the majority of angiosperm lineages are derived from diverse parental genomes and may be of allopolyploid origin.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

The link between ancient whole-genome duplications and cold adaptations in the Caryophyllaceae

PREMISE

The Caryophyllaceae (the carnation family) have undergone multiple transitions into colder climates and convergence on cushion plant adaptation, indicating that they may provide a natural system for cold adaptation research. Previous research has suggested that putative ancient whole-genome duplications (WGDs) are correlated with niche shifts into colder climates across the Caryophyllales. Here, we explored the genomic changes potentially involved in one of these discovered shifts in the Caryophyllaceae.

METHODS

We constructed a data set combining 26 newly generated transcriptomes with 45 published transcriptomes, including 11 cushion plant species across seven genera. With this data set, we inferred a dated phylogeny for the Caryophyllaceae and mapped ancient WGDs and gene duplications onto the phylogeny. We also examined functional groups enriched for gene duplications related to the climatic shift.

RESULTS

The ASTRAL topology was mostly congruent with the current consensus of relationships within the family. We inferred 15 putative ancient WGDs in the family, including eight that have not been previously published. The oldest ancient WGD (ca. 64.4-56.7 million years ago), WGD1, was found to be associated with a shift into colder climates by previous research. Gene regions associated with ubiquitination were overrepresented in gene duplications retained after WGD1 and those convergently retained by cushion plants in Colobanthus and Eremogone, along with other functional annotations.

CONCLUSIONS

Gene family expansions induced by ancient WGDs may have contributed to the shifts to cold climatic niches in the Caryophyllaceae. Transcriptomic data are crucial resources that help unravel heterogeneity in deep-time evolutionary patterns in plants.

Doubling down on polyploid discoveries: Global advances in genomics and ecological impacts of polyploidy

All flowering plants are now recognized as diploidized paleopolyploids (Jiao et al., 2011; One Thousand Plant Transcriptomes Initiative, 2019), and polyploid species comprise approximately 30% of contemporary plant species (Wood et al., 2009; Barker et al., 2016a). A major implication of these discoveries is that, to appreciate the evolution of plant diversity, we need to understand the fundamental biology of polyploids and diploidization. This need is broadly recognized by our community as there is a continued, growing interest in polyploidy as a research topic. Over the past 25 years, the sequencing and analysis of plant genomes has revolutionized our understanding of the importance of polyploid speciation to the evolution of land plants.

Species-tree topology impacts the inference of ancient whole-genome duplications across the angiosperm phylogeny

PREMISE

The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.

METHODS

We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.

RESULTS

The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.

CONCLUSIONS

Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

The chromosome-scale assembly of the Notopterygium incisum genome provides insight into the structural diversity of coumarins

Coumarins, derived from the phenylpropanoid pathway, represent one of the primary metabolites found in angiosperms. The alignment of the tetrahydropyran (THP) and tetrahydrofuran (THF) rings with the lactone structure results in the formation of at least four types of complex coumarins. However, the mechanisms underlying the structural diversity of coumarin remain poorly understood. Here, we report the chromosome-level genome assembly of Notopterygium incisum, spanning 1.64 Gb, with a contig N50 value of 22.7 Mb and 60,021 annotated protein-coding genes. Additionally, we identified the key enzymes responsible for shaping the structural diversity of coumarins, including two p-coumaroyl CoA 2'-hydroxylases crucial for simple coumarins basic skeleton architecture, two UbiA prenyltransferases responsible for angular or linear coumarins biosynthesis, and five CYP736 cyclases involved in THP and THF ring formation. Notably, two bifunctional enzymes capable of catalyzing both demethylsuberosin and osthenol were identified for the first time. Evolutionary analysis implies that tandem and ectopic duplications of the CYP736 subfamily, specifically arising in the Apiaceae, contributed to the structural diversity of coumarins in N. incisum. Conclusively, this study proposes a parallel evolution scenario for the complex coumarin biosynthetic pathway among different angiosperms and provides essential synthetic biology elements for the heterologous industrial production of coumarins.

wgd v2: a suite of tools to uncover and date ancient polyploidy and whole-genome duplication

MOTIVATION

Major improvements in sequencing technologies and genome sequence assembly have led to a huge increase in the number of available genome sequences. In turn, these genome sequences form an invaluable source for evolutionary, ecological, and comparative studies. One kind of analysis that has become routine is the search for traces of ancient polyploidy, particularly for plant genomes, where whole-genome duplication (WGD) is rampant.

RESULTS

Here, we present a major update of a previously developed tool wgd, namely wgd v2, to look for remnants of ancient polyploidy, or WGD. We implemented novel and improved previously developed tools to (a) construct KS age distributions for the whole-paranome (collection of all duplicated genes in a genome), (b) unravel intragenomic and intergenomic collinearity resulting from WGDs, (c) fit mixture models to age distributions of gene duplicates, (d) correct substitution rate variation for phylogenetic placement of WGDs, and (e) date ancient WGDs via phylogenetic dating of WGD-retained gene duplicates. The applicability and feasibility of wgd v2 for the identification and the relative and absolute dating of ancient WGDs is demonstrated using different plant genomes.

AVAILABILITY AND IMPLEMENTATION

wgd v2 is open source and available at https://github.com/heche-psb/wgd.

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

Crocus genome reveals the evolutionary origin of crocin biosynthesis

Crocus sativus (saffron) is a globally autumn-flowering plant, and its stigmas are the most expensive spice and valuable herb medicine. Crocus specialized metabolites, crocins, are biosynthesized in distant species, Gardenia (eudicot) and Crocus (monocot), and the evolution of crocin biosynthesis remains poorly understood. With the chromosome-level Crocus genome assembly, we revealed that two rounds of lineage-specific whole genome triplication occurred, contributing important roles in the production of carotenoids and apocarotenoids. According to the kingdom-wide identification, phylogenetic analysis, and functional assays of carotenoid cleavage dioxygenases (CCDs), we deduced that the duplication, site positive selection, and neofunctionalization of Crocus-specific CCD2 from CCD1 members are responsible for the crocin biosynthesis. In addition, site mutation of CsCCD2 revealed the key amino acids, including I143, L146, R161, E181, T259, and S292 related to the catalytic activity of zeaxanthin cleavage. Our study provides important insights into the origin and evolution of plant specialized metabolites, which are derived by duplication events of biosynthetic genes.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Chromosome-level genome assembly of Niphotrichum japonicum provides new insights into heat stress responses in mosses

With a diversity of approximately 22,000 species, bryophytes (hornworts, liverworts, and mosses) represent a major and diverse lineage of land plants. Bryophytes can thrive in many extreme environments as they can endure the stresses of drought, heat, and cold. The moss Niphotrichum japonicum (Grimmiaceae, Grimmiales) can subsist for extended periods under heat and drought conditions, providing a good candidate for studying the genetic basis underlying such high resilience. Here, we de novo assembled the genome of N. japonicum using Nanopore long reads combined with Hi-C scaffolding technology to anchor the 191.61 Mb assembly into 14 pseudochromosomes. The genome structure of N. japonicum's autosomes is mostly conserved and highly syntenic, in contrast to the sparse and disordered genes present in its sex chromosome. Comparative genomic analysis revealed the presence of 10,019 genes exclusively in N. japonicum. These genes may contribute to the species-specific resilience, as demonstrated by the gene ontology (GO) enrichment. Transcriptome analysis showed that 37.44% (including 3,107 unique genes) of the total annotated genes (26,898) exhibited differential expression as a result of heat-induced stress, and the mechanisms that respond to heat stress are generally conserved across plants. These include the upregulation of HSPs, LEAs, and reactive oxygen species (ROS) scavenging genes, and the downregulation of PPR genes. N. japonicum also appears to have distinctive thermal mechanisms, including species-specific expansion and upregulation of the Self-incomp_S1 gene family, functional divergence of duplicated genes, structural clusters of upregulated genes, and expression piggybacking of hub genes. Overall, our study highlights both shared and species-specific heat tolerance strategies in N. japonicum, providing valuable insights into the heat tolerance mechanism and the evolution of resilient plants.

A first draft genome of holm oak (Quercus ilex subsp. ballota), the most representative species of the Mediterranean forest and the Spanish agrosylvopastoral ecosystem "dehesa"

The holm oak (Quercus ilex subsp. ballota) is the most representative species of the Mediterranean Basin and the agrosylvopastoral Spanish "dehesa" ecosystem. Being part of our life, culture, and subsistence since ancient times, it has significant environmental and economic importance. More recently, there has been a renewed interest in using the Q. ilex acorn as a functional food due to its nutritional and nutraceutical properties. However, the holm oak and its related ecosystems are threatened by different factors, with oak decline syndrome and climate change being the most worrying in the short and medium term. Breeding programs informed by the selection of elite genotypes seem to be the most plausible biotechnological solution to rescue populations under threat. To achieve this and other downstream analyses, we need a high-quality and well-annotated Q. ilex reference genome. Here, we introduce the first draft genome assembly of Q. ilex using long-read sequencing (PacBio). The assembled nuclear haploid genome had 530 contigs totaling 842.2 Mbp (N50 = 3.3 Mbp), of which 448.7 Mb (53%) were repetitive sequences. We annotated 39,443 protein-coding genes of which 94.80% were complete and single-copy genes. Phylogenetic analyses showed no evidence of a recent whole-genome duplication, and high synteny of the 12 chromosomes between Q. ilex and Quercus lobata and between Q. ilex and Quercus robur. The chloroplast genome size was 142.3 Kbp with 149 protein-coding genes successfully annotated. This first draft should allow for the validation of omics data as well as the identification and functional annotation of genes related to phenotypes of interest such as those associated with resilience against oak decline syndrome and climate change and higher acorn productivity and nutraceutical value.

Inference of Ancient Polyploidy from Genomic Data

Whole-genome sequence data have revealed that numerous eukaryotic organisms derive from distant polyploid ancestors, even when these same organisms are genetically and karyotypically diploid. Such ancient whole-genome duplications (WGDs) have been important for long-term genome evolution and are often speculatively associated with important evolutionary events such as key innovations, adaptive radiations, or survival after mass extinctions. Clearly, reliable methods for unveiling ancient WGDs are key toward furthering understanding of the long-term evolutionary significance of polyploidy. In this chapter, we describe a set of basic established comparative genomics approaches for the inference of ancient WGDs from genomic data based on empirical age distributions and collinearity analyses, explain the principles on which they are based, and illustrate a basic workflow using the software "wgd," geared toward a typical exploratory analysis of a newly obtained genome sequence.

Revisiting ancient polyploidy in leptosporangiate ferns

Ferns, and particularly homosporous ferns, have long been assumed to have experienced recurrent whole-genome duplication (WGD) events because of their substantially large genome sizes, surprisingly high chromosome numbers, and high degrees of polyploidy among many extant members. As the number of sequenced fern genomes is limited, recent studies have employed transcriptome data to find evidence for WGDs in ferns. However, they have reached conflicting results concerning the occurrence of ancient polyploidy, for instance, in the lineage of leptosporangiate ferns. Because identifying WGDs in a phylogenetic context is the foremost step in studying the contribution of ancient polyploidy to evolution, we here revisited earlier identified WGDs in leptosporangiate ferns, mainly the core leptosporangiate ferns, by building K_S -age distributions and applying substitution rate corrections and by conducting statistical gene tree-species tree reconciliation analyses. Our integrative analyses not only identified four ancient WGDs in the sampled core leptosporangiate ferns but also identified false positives and false negatives for WGDs that recent studies have reported earlier. In conclusion, we underscore the significance of substitution rate corrections and uncertainties in gene tree-species tree reconciliations in calling WGD events and advance an exemplar workflow to overcome such often-overlooked issues.

A genome resource for Acacia, Australia's largest plant genus

Acacia (Leguminosae, Caesalpinioideae, mimosoid clade) is the largest and most widespread genus of plants in the Australian flora, occupying and dominating a diverse range of environments, with an equally diverse range of forms. For a genus of its size and importance, Acacia currently has surprisingly few genomic resources. Acacia pycnantha, the golden wattle, is a woody shrub or tree occurring in south-eastern Australia and is the country's floral emblem. To assemble a genome for A. pycnantha, we generated long-read sequences using Oxford Nanopore Technology, 10x Genomics Chromium linked reads, and short-read Illumina sequences, and produced an assembly spanning 814 Mb, with a scaffold N50 of 2.8 Mb, and 98.3% of complete Embryophyta BUSCOs. Genome annotation predicted 47,624 protein-coding genes, with 62.3% of the genome predicted to comprise transposable elements. Evolutionary analyses indicated a shared genome duplication event in the Caesalpinioideae, and conflict in the relationships between Cercis (subfamily Cercidoideae) and subfamilies Caesalpinioideae and Papilionoideae (pea-flowered legumes). Comparative genomics identified a suite of expanded and contracted gene families in A. pycnantha, and these were annotated with both GO terms and KEGG functional categories. One expanded gene family of particular interest is involved in flowering time and may be associated with the characteristic synchronous flowering of Acacia. This genome assembly and annotation will be a valuable resource for all studies involving Acacia, including the evolution, conservation, breeding, invasiveness, and physiology of the genus, and for comparative studies of legumes.

The genome of homosporous maidenhair fern sheds light on the euphyllophyte evolution and defences

Euphyllophytes encompass almost all extant plants, including two sister clades, ferns and seed plants. Decoding genomes of ferns is the key to deep insight into the origin of euphyllophytes and the evolution of seed plants. Here we report a chromosome-level genome assembly of Adiantum capillus-veneris L., a model homosporous fern. This fern genome comprises 30 pseudochromosomes with a size of 4.8-gigabase and a contig N50 length of 16.22 Mb. Gene co-expression network analysis uncovered that homospore development in ferns has relatively high genetic similarities with that of the pollen in seed plants. Analysing fern defence response expands understanding of evolution and diversity in endogenous bioactive jasmonates in plants. Moreover, comparing fern genomes with those of other land plants reveals changes in gene families important for the evolutionary novelties within the euphyllophyte clade. These results lay a foundation for studies on fern genome evolution and function, as well as the origin and evolution of euphyllophytes.

Phylotranscriptomics Illuminates the Placement of Whole Genome Duplications and Gene Retention in Ferns

Ferns are the second largest clade of vascular plants with over 10,000 species, yet the generation of genomic resources for the group has lagged behind other major clades of plants. Transcriptomic data have proven to be a powerful tool to assess phylogenetic relationships, using thousands of markers that are largely conserved across the genome, and without the need to sequence entire genomes. We assembled the largest nuclear phylogenetic dataset for ferns to date, including 2884 single-copy nuclear loci from 247 transcriptomes (242 ferns, five outgroups), and investigated phylogenetic relationships across the fern tree, the placement of whole genome duplications (WGDs), and gene retention patterns following WGDs. We generated a well-supported phylogeny of ferns and identified several regions of the fern phylogeny that demonstrate high levels of gene tree-species tree conflict, which largely correspond to areas of the phylogeny that have been difficult to resolve. Using a combination of approaches, we identified 27 WGDs across the phylogeny, including 18 large-scale events (involving more than one sampled taxon) and nine small-scale events (involving only one sampled taxon). Most inferred WGDs occur within single lineages (e.g., orders, families) rather than on the backbone of the phylogeny, although two inferred events are shared by leptosporangiate ferns (excluding Osmundales) and Polypodiales (excluding Lindsaeineae and Saccolomatineae), clades which correspond to the majority of fern diversity. We further examined how retained duplicates following WGDs compared across independent events and found that functions of retained genes were largely convergent, with processes involved in binding, responses to stimuli, and certain organelles over-represented in paralogs while processes involved in transport, organelles derived from endosymbiotic events, and signaling were under-represented. To date, our study is the most comprehensive investigation of the nuclear fern phylogeny, though several avenues for future research remain unexplored.

Chromosome Level Genome Assembly and Annotation of Highly Invasive Japanese Stiltgrass (Microstegium vimineum)

The invasive Japanese stiltgrass (Microstegium vimineum) affects a wide range of ecosystems and threatens biodiversity across the eastern USA. However, the mechanisms underlying rapid adaptation, plasticity, and epigenetics in the invasive range are largely unknown. We present a chromosome-level assembly for M. vimineum to investigate genome dynamics, evolution, adaptation, and the genomics of phenotypic plasticity. We generated a 1.12-Gb genome with scaffold N50 length of 53.44 Mb respectively, taking a de novo assembly approach that combined PacBio and Dovetail Genomics Omni-C sequencing. The assembly contains 23 pseudochromosomes, representing 99.96% of the genome. BUSCO assessment indicated that 80.3% of Poales gene groups are present in the assembly. The genome is predicted to contain 39,604 protein-coding genes, of which 26,288 are functionally annotated. Furthermore, 66.68% of the genome is repetitive, of which unclassified (35.63%) and long-terminal repeat (LTR) retrotransposons (26.90%) are predominant. Similar to other grasses, Gypsy (41.07%) and Copia (32%) are the most abundant LTR-retrotransposon families. The majority of LTR-retrotransposons are derived from a significant expansion in the past 1-2 Myr, suggesting the presence of relatively young LTR-retrotransposon lineages. We find corroborating evidence from Ks plots for a stiltgrass-specific duplication event, distinct from the more ancient grass-specific duplication event. The assembly and annotation of M. vimineum will serve as an essential genomic resource facilitating studies of the invasion process, the history and consequences of polyploidy in grasses, and provides a crucial tool for natural resource managers.