The genome sequence of lesser trefoil or Irish shamrock, Trifolium dubium Sibth. (Fabaceae).

We present a genome assembly from an individual Trifolium dubium (lesser trefoil; Tracheophyta; Magnoliopsida; Fabales; Fabaceae) as part of a collaboration between the Darwin Tree of Life and the European Reference Genome Atlas. The genome sequence is 679.1 megabases in span. Most of the assembly is scaffolded into 15 chromosomal pseudomolecules. The two mitochondrial genomes have lengths of 133.86 kb and 182.32 kb, and the plastid genome assembly has a length of 126.22 kilobases.

The genome sequence of the Lesser Skullcap, Scutellaria minor Huds., 1762 (Lamiaceae).

We present a genome assembly from an individual Scutellaria minor (Tracheophyta; Magnoliopsida; Lamiales; Lamiaceae). The genome sequence is 341.8 megabases in span. Most of the assembly is scaffolded into 14 chromosomal pseudomolecules. The mitochondrial and plastid genome assemblies have lengths of 376.64 kilobases and 152.59 kilobases in length, respectively.

A 160 Gbp fork fern genome shatters size record for eukaryotes.

Vascular plants are exceptional among eukaryotes due to their outstanding genome size diversity which ranges ∼2,400-fold, including the largest genome so far recorded in the angiosperm Paris japonica (148.89 Gbp/1C). Despite available data showing that giant genomes are restricted across the Tree of Life, the biological limits to genome size expansion remain to be established. Here, we report the discovery of an even larger eukaryotic genome in Tmesipteris oblanceolata, a New Caledonian fork fern. At 160.45 Gbp/1C, this record-breaking genome challenges current understanding and opens new avenues to explore the evolutionary dynamics of genomic gigantism.

Biased Retention of Environment-Responsive Genes Following Genome Fractionation.

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

Genome size is positively correlated with extinction risk in herbaceous angiosperms.

Angiosperms with large genomes experience nuclear-, cellular-, and organism-level constraints that may limit their phenotypic plasticity and ecological niche, which could increase their risk of extinction. Therefore, we test the hypotheses that large-genomed species are more likely to be threatened with extinction than those with small genomes, and that the effect of genome size varies across three selected covariates: life form, endemism, and climatic zone. We collated genome size and extinction risk information for a representative sample of angiosperms comprising 3250 species, which we analyzed alongside life form, endemism, and climatic zone variables using a phylogenetic framework. Genome size is positively correlated with extinction risk, a pattern driven by a signal in herbaceous but not woody species, regardless of climate and endemism. The influence of genome size is stronger in endemic herbaceous species, but is relatively homogenous across different climates. Beyond its indirect link via endemism and climate, genome size is associated with extinction risk directly and significantly. Genome size may serve as a proxy for difficult-to-measure parameters associated with resilience and vulnerability in herbaceous angiosperms. Therefore, it merits further exploration as a useful biological attribute for understanding intrinsic extinction risk and augmenting plant conservation efforts.

Genomic incongruence accompanies the evolution of flower symmetry in Eudicots: a case study in the poppy family (Papaveraceae, Ranunculales).

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events.

The genome sequence of rosebay willowherb Chamaenerion angustifolium (L.) Scop., 1771 (syn. Epilobium angustifolium L., 1753) (Onagraceae).

We present a genome assembly from an individual Chamaenerion angustifolium (fireweed; Tracheophyta; Magnoliopsida; Myrtales; Onagraceae). The genome sequence is 655.9 megabases in span. Most of the assembly is scaffolded into 18 chromosomal pseudomolecules. The mitochondrial and plastid genome assemblies have lengths of 495.18 kilobases and 160.41 kilobases in length, respectively.

The rise of baobab trees in Madagascar.

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Phylogeny, biogeography and ecological diversification of New Caledonian palms (Arecaceae).

BACKGROUND AND AIMS: The geographical origin and evolutionary mechanisms underpinning the rich and distinctive New Caledonian flora remain poorly understood. This is attributable to the complex geological past of the island and to the scarcity of well-resolved species-level phylogenies. Here, we infer phylogenetic relationships and divergence times of New Caledonian palms, which comprise 40 species. We use this framework to elucidate the biogeography of New Caledonian palm lineages and to explore how extant species might have formed. METHODS: A phylogenetic tree including 37 New Caledonian palm species and 77 relatives from tribe Areceae was inferred from 151 nuclear genes obtained by targeted sequencing. Fossil-calibrated divergence times were estimated and ancestral ranges inferred. Ancestral and extant ecological preferences in terms of elevation, precipitation and substrate were compared between New Caledonian sister species to explore their possible roles as drivers of speciation. KEY RESULTS: New Caledonian palms form four well-supported clades, inside which relationships are well resolved. Our results support the current classification but suggest that Veillonia and Campecarpus should be resurrected and fail to clarify whether Rhopalostylidinae is sister to or nested in Basseliniinae. New Caledonian palm lineages are derived from New Guinean and Australian ancestors, which reached the island through at least three independent dispersal events between the Eocene and Miocene. Palms then dispersed out of New Caledonia at least five times, mainly towards Pacific islands. Geographical and ecological transitions associated with speciation events differed across time and genera. Substrate transitions were more frequently associated with older events than with younger ones. CONCLUSIONS: Neighbouring areas and a mosaic of local habitats shaped the palm flora of New Caledonia, and the island played a significant role in generating palm diversity across the Pacific region. This new spatio-temporal framework will enable population-level ecological and genetic studies to unpick the mechanisms underpinning New Caledonian palm endemism.

Plant invasion and naturalization are influenced by genome size, ecology and economic use globally.

Human factors and plant characteristics are important drivers of plant invasions, which threaten ecosystem integrity, biodiversity and human well-being. However, while previous studies often examined a limited number of factors or focused on a specific invasion stage (e.g., naturalization) for specific regions, a multi-factor and multi-stage analysis at the global scale is lacking. Here, we employ a multi-level framework to investigate the interplay between plant characteristics (genome size, Grime's adaptive CSR-strategies and native range size) and economic use and how these factors collectively affect plant naturalization and invasion success worldwide. While our findings derived from structural equation models highlight the substantial contribution of human assistance in both the naturalization and spread of invasive plants, we also uncovered the pivotal role of species' adaptive strategies among the factors studied, and the significantly varying influence of these factors across invasion stages. We further revealed that the effects of genome size on plant invasions were partially mediated by species adaptive strategies and native range size. Our study provides insights into the complex and dynamic process of plant invasions and identifies its key drivers worldwide.

The origin and speciation of orchids.

Orchids constitute one of the most spectacular radiations of flowering plants. However, their origin, spread across the globe, and hotspots of speciation remain uncertain due to the lack of an up-to-date phylogeographic analysis. We present a new Orchidaceae phylogeny based on combined high-throughput and Sanger sequencing data, covering all five subfamilies, 17/22 tribes, 40/49 subtribes, 285/736 genera, and c. 7% (1921) of the 29 524 accepted species, and use it to infer geographic range evolution, diversity, and speciation patterns by adding curated geographical distributions from the World Checklist of Vascular Plants. The orchids' most recent common ancestor is inferred to have lived in Late Cretaceous Laurasia. The modern range of Apostasioideae, which comprises two genera with 16 species from India to northern Australia, is interpreted as relictual, similar to that of numerous other groups that went extinct at higher latitudes following the global climate cooling during the Oligocene. Despite their ancient origin, modern orchid species diversity mainly originated over the last 5 Ma, with the highest speciation rates in Panama and Costa Rica. These results alter our understanding of the geographic origin of orchids, previously proposed as Australian, and pinpoint Central America as a region of recent, explosive speciation.

The global distribution of angiosperm genome size is shaped by climate.

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

Complementing model species with model clades.

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant Tree of Life continues to improve. The intersection of these two research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade". These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis thaliana and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

The genome sequence of the tree of heaven, Ailanthus altissima (Mill.) Swingle, 1916.

We present a genome assembly from an individual Ailanthus altissima (tree of heaven; Streptophyta; Magnoliopsida; Sapindales; Simaroubaceae). The genome sequence is 939 megabases in span. Most of the assembly is scaffolded into 31 chromosomal pseudomolecules. The mitochondrial and plastid genome assemblies are 661.1 kilobases and 161.1 kilobases long, respectively.

From genome size to trait evolution during angiosperm radiation.

Angiosperm diversity arises from trait flexibility and repeated evolutionary radiations, but the role of genomic characters in these radiations remains unclear. In this opinion article, we discuss how genome size can influence angiosperm diversification via its intricate link with cell size, tissue packing, and physiological processes which, in turn, influence the macroevolution of functional traits. We propose that integrating genome size, functional traits, and phylogenetic data across a wide range of lineages allows us to test whether genome size decrease consistently leads to increased trait flexibility, while genome size increase constrains trait evolution. Combining theories from molecular biology, functional ecology and macroevolution, we provide a framework to better understand the role of genome size in trait evolution, evolutionary radiations, and the global distribution of angiosperms.

"BIFloraExplorer": A Taxonomic, Genetic, and Ecological Data Resource for the Vascular Plants of Britain and Ireland.

The vascular flora of Britain and Ireland is a historically well-documented and clearly delimited study system that offers itself to large-scale analyses of ecology and species assemblages. However, such analyses require clean, curated, and taxonomically resolved data, which are often unavailable. In this chapter, we describe how to access and use a key data resource that combines a taxonomically stable species list with genetic data (genome size, chromosome counts, and DNA barcode information), ecological information (such as life-form, realized niche description, and geographic origin) and distribution records. The data resource enables and encourages the study of natural ecological and evolutionary patterns and processes within the vascular flora of Britain and Ireland.

The Plant DNA C-Values Database: A One-Stop Shop for Plant Genome Size Data.

Genome size is a plant character with far-reaching implications, ranging from impacts on the financial and computing feasibility of sequencing and assembling genomes all the way to influencing the very ecology and evolution of species. The increasing recognition of the role of genome size in plant science has led to a rising demand for comprehensive and easily accessible sources of genome size data. The Plant DNA C-values database has established itself as a trusted and widely used central hub for users needing to access available plant genome size data, complemented with related cytogenetic (ploidy level) and karyological (chromosome number) information where available. Since its inception in 2001, the database has undergone six major updates to incorporate newly available genome size information, leading to the most recent release (Release 7.1), which comprises data for 12,273 species across all the major land plant and some algal lineages. Here we describe how to use the database efficiently, making use of its different query and filtering settings.

Small genome size and variation in ploidy levels support the naturalization of vascular plants but constrain their invasive spread.

Karyological characteristics are among the traits underpinning the invasion success of vascular plants. Using 11 049 species, we tested the effects of genome size and ploidy levels on plant naturalization (species forming self-sustaining populations where they are not native) and invasion (naturalized species spreading rapidly and having environmental impact). The probability that a species naturalized anywhere in the world decreased with increasing monoploid genome size (DNA content of a single chromosome set). Naturalized or invasive species with intermediate monoploid genomes were reported from many regions, but those with either small or large genomes occurred in fewer regions. By contrast, large holoploid genome sizes (DNA content of the unreplicated gametic nucleus) constrained naturalization but favoured invasion. We suggest that a small genome is an advantage during naturalization, being linked to traits favouring adaptation to local conditions, but for invasive spread, traits associated with a large holoploid genome, where the impact of polyploidy may act, facilitate long-distance dispersal and competition with other species.

The Role of Chromatin Modifications in the Evolution of Giant Plant Genomes.

Angiosperm genome sizes (GS) range ~2400-fold and comprise genes and their regulatory regions, repeats, semi-degraded repeats, and 'dark matter'. The latter represents repeats so degraded that they can no longer be recognised as repetitive. In exploring whether the histone modifications associated with chromatin packaging of these contrasting genomic components are conserved across the diversity of GS in angiosperms, we compared immunocytochemistry data for two species whose GS differ ~286-fold. We compared published data for Arabidopsis thaliana with a small genome (GS = 157 Mbp/1C) with newly generated data from Fritillaria imperialis, which has a giant genome (GS = 45,000 Mbp/1C). We compared the distributions of the following histone marks: H3K4me1, H3K4me2, H3K9me1, H3K9me2, H3K9me3, H3K27me1, H3K27me2, and H3K27me3. Assuming these histone marks are associated with the same genomic features across all species, irrespective of GS, our comparative analysis enables us to suggest that while H3K4me1 and H3K4me2 methylation identifies genic DNA, H3K9me3 and H3K27me3 marks are associated with 'dark matter', H3K9me1 and H3K27me1 mark highly homogeneous repeats, and H3K9me2 and H3K27me2 mark semi-degraded repeats. The results have implications for our understanding of epigenetic profiles, chromatin packaging and the divergence of genomes, and highlight contrasting organizations of the chromatin within the nucleus depending on GS itself.

A Bioinformatic Pipeline to Estimate Ploidy Level from Target Capture Sequence Data Obtained from Herbarium Specimens.

Whole genome duplications (WGD) are frequent in many plant lineages; however, ploidy level variation is unknown in most species. The most widely used methods to estimate ploidy levels in plants are chromosome counts, which require living specimens, and flow cytometry estimates, which necessitate living or relatively recently collected samples. Newly described bioinformatic methods have been developed to estimate ploidy levels using high-throughput sequencing data, and these have been optimized in plants by calculating allelic ratio values from target capture data. This method relies on the maintenance of allelic ratios from the genome to the sequence data. For example, diploid organisms will generate allelic data in a 1:1 proportion, with an increasing number of possible allelic ratio combinations occurring in individuals with higher ploidy levels. In this chapter, we explain step-by-step this bioinformatic approach for the estimation of ploidy level.