Exploring Angiosperms353: Developing and applying a universal toolkit for flowering plant phylogenomics

Chloroplast genomic insights into adaptive evolution and rapid radiation in the genus Passiflora (Passifloraceae)

Chloroplasts are essential organelles in plants and eukaryotic algae, responsible for photosynthesis, fatty acid synthesis, amino acid production, and stress responses. The genus Passiflora, known for its species diversity and dynamic chloroplast (cp) genome evolution, serves as an excellent model for studying structural variations. This study investigates evolutionary relationships within Passiflora by sequencing 11 new chloroplast genomes, assessing selective pressures on cp genes, and comparing plastid and nuclear phylogenies. Passiflora cp genomes showed significant variations in size, gene content, and structure, ranging from 132,736 to 163,292 base pairs, especially in Decaloba. Structural rearrangements and species-specific repeat patterns were identified. Selective pressure tests revealed significant adaptive evolution in certain lineages, with several genes, including clpP and petL, under positive selection. Phylogenetic analyses confirmed the monophyly of subgenera Astrophea, Passiflora, and Decaloba, while Deidamioides appeared polyphyletic. Nuclear phylogenetic analysis based on 35S rDNA sequences supported the monophyly of Astrophea but showed inconsistencies within subgenus Passiflora compared to cp genome data. This study highlights the evolutionary complexity of Passiflora cp genomes, demonstrating significant structural variations and adaptive evolution. The findings underscore the effectiveness of plastid phylogenomics in resolving phylogenetic relationships and provide insights into adaptive mechanisms shaping cp genome diversity in angiosperms.

Toward a phylogenomic classification of magnoliids

PREMISE

Magnoliids are a strongly supported clade of angiosperms. Previous phylogenetic studies based primarily on analyses of a limited number of mostly plastid markers have led to the current classification of magnoliids into four orders and 18 families. However, uncertainty remains regarding the placement of several families.

METHODS

For the first comprehensive phylogenomic analysis of magnoliids as a whole, we sampled 235 species from 199 (74%) genera and representing all families and most subfamilies and tribes. We analyzed newly generated data from the Angiosperms353 probe set using both coalescent and concatenation analyses and testing the impact of multiple filtering and alignment strategies.

RESULTS

While our results generally provide further support for previously established phylogenetic relationships in both magnoliids as a whole and large families including Annonaceae and Lauraceae, they also provide new evidence for previously ambiguous relationships. In particular, we found support for the position of Hydnoraceae as sister to the remainder of Piperales, though evidence was conflicting, and resolved the backbone of relationships among most genera of Myristicaceae. Different analytical strategies tended to have rather small effects on branch support and topology.

CONCLUSIONS

Although some of our results are limited by low gene recovery for a number of taxa and significant gene tree conflict for some relationships, this study represents a significant step toward reconstructing the evolutionary history of a major lineage of angiosperms. Based on these results, we present an updated phylogenetic classification for magnoliids, recognizing 21 families, summarizing previously established subfamilies and tribes, and describing new tribes for Myristicaceae.

Phylogenomics of the pantropical Connaraceae: revised infrafamilial classification and the evolution of heterostyly

Connaraceae is a pantropical family of about 200 species containing lianas and small trees with remarkably diverse floral polymorphisms, including distyly, tristyly, homostyly, and dioecy. To date, relationships within the family have not been investigated using a targeted molecular phylogenetic treatment, severely limiting systematic understanding and reconstruction of trait evolution. Accordingly, their last infrafamilial classification was based only on morphological data. Here, we used phylogenomic data obtained using the Angiosperms353 nuclear target sequence capture probes, sampling all tribes and almost all genera, entirely from herbarium specimens, to revise infrafamilial classification and investigate the evolution of heterostyly. The backbone of the resulting molecular phylogenetic tree is almost entirely resolved. Connaraceae consists of two clades, one containing only the African genus Manotes (4 or 5 species), which we newly recognize at the subfamily level. Vegetative and reproductive synapomorphies are proposed for Manotoideae. Within Connaroideae, Connareae is expanded to include the former Jollydoreae. The backbone of Cnestideae, which contains more than half of the Connaraceae species, remains incompletely resolved. Reconstructions of reproductive system evolution are presented that tentatively support tristyly as the ancestral state for the family, with multiple parallel losses, in agreement with previous hypotheses, plus possible re-gains. However, the great diversity of stylar polymorphisms and their phylogenetic lability preclude a definitive answer. Overall, this study reinforces the usefulness of herbarium phylogenomics, and unlocks the reproductive diversity of Connaraceae as a model system for the evolution of complex biological phenomena.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00606-024-01909-y.

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.

Phylogenomics and the rise of the angiosperms

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Melica as an emerging model system for comparative studies in temperate Pooideae grasses

BACKGROUND AND AIMS

Pooideae grasses contain some of the world's most important crop and forage species. Although much work has been conducted on understanding the genetic basis of trait diversification within a few annual Pooideae, comparative studies at the subfamily level are limited by a lack of perennial models outside 'core' Pooideae. We argue for development of the perennial non-core genus Melica as an additional model for Pooideae, and provide foundational data regarding the group's biogeography and history of character evolution.

METHODS

Supplementing available ITS and ndhF sequence data, we built a preliminary Bayesian-based Melica phylogeny, and used it to understand how the genus has diversified in relation to geography, climate and trait variation surveyed from various floras. We also determine biomass accumulation under controlled conditions for Melica species collected across different latitudes and compare inflorescence development across two taxa for which whole genome data are forthcoming.

KEY RESULTS

Our phylogenetic analyses reveal three strongly supported geographically structured Melica clades that are distinct from previously hypothesized subtribes. Despite less geographical affinity between clades, the two sister 'Ciliata' and 'Imperfecta' clades segregate from the more phylogenetically distant 'Nutans' clade in thermal climate variables and precipitation seasonality, with the 'Imperfecta' clade showing the highest levels of trait variation. Growth rates across Melica are positively correlated with latitude of origin. Variation in inflorescence morphology appears to be explained largely through differences in secondary branch distance, phyllotaxy and number of spikelets per secondary branch.

CONCLUSIONS

The data presented here and in previous studies suggest that Melica possesses many of the necessary features to be developed as an additional model for Pooideae grasses, including a relatively fast generation time, perenniality, and interesting variation in physiology and morphology. The next step will be to generate a genome-based phylogeny and transformation tools for functional analyses.

Horizon scan of DNA-based methods for quality control and monitoring of herbal preparations

Herbal medicines and preparations are widely used in healthcare systems globally, but concerns remain about their quality and safety. New herbal products are constantly being introduced to the market under varying regulatory frameworks, with no global consensus on their definition or characterization. These biologically active mixtures are sold through complex globalized value chains, which create concerns around contamination and profit-driven adulteration. Industry, academia, and regulatory bodies must collaborate to develop innovative strategies for the identification and authentication of botanicals and their preparations to ensure quality control. High-throughput sequencing (HTS) has significantly improved our understanding of the total species diversity within DNA mixtures. The standard concept of DNA barcoding has evolved over the last two decades to encompass genomic data more broadly. Recent research in DNA metabarcoding has focused on developing methods for quantifying herbal product ingredients, yielding meaningful results in a regulatory framework. Techniques, such as loop-mediated isothermal amplification (LAMP), DNA barcode-based Recombinase Polymerase Amplification (BAR-RPA), DNA barcoding coupled with High-Resolution Melting (Bar-HRM), and microfluidics-based methods, offer more affordable tests for the detection of target species. While target capture sequencing and genome skimming are considerably increasing the species identification resolution in challenging plant clades, ddPCR enables the quantification of DNA in samples and could be used to detect intended and unwanted ingredients in herbal medicines. Here, we explore the latest advances in emerging DNA-based technologies and the opportunities they provide as taxa detection tools for evaluating the safety and quality of dietary supplements and herbal medicines.

Phylogenomic analyses of Sapindales support new family relationships, rapid Mid-Cretaceous Hothouse diversification, and heterogeneous histories of gene duplication

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order's spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research.

Target sequence data shed new light on the infrafamilial classification of Araceae

PREMISE

Recent phylogenetic studies of the Araceae have confirmed the position of the duckweeds nested within the aroids, and the monophyly of a clade containing all the unisexual flowered aroids plus the bisexual-flowered Calla palustris. The main objective of the present study was to better resolve the deep phylogenetic relationships among the main lineages within the family, particularly the relationships between the eight currently recognized subfamilies. We also aimed to confirm the phylogenetic position of the enigmatic genus Calla in relation to the long-debated evolutionary transition between bisexual and unisexual flowers in the family.

METHODS

Nuclear DNA sequence data were generated for 128 species across 111 genera (78%) of Araceae using target sequence capture and the Angiosperms 353 universal probe set.

RESULTS

The phylogenomic data confirmed the monophyly of the eight Araceae subfamilies, but the phylogenetic position of subfamily Lasioideae remains uncertain. The genus Calla is included in subfamily Aroideae, which has also been expanded to include Zamioculcadoideae. The tribe Aglaonemateae is newly defined to include the genera Aglaonema and Boycea.

CONCLUSIONS

Our results strongly suggest that new research on African genera (Callopsis, Nephthytis, and Anubias) and Calla will be important for understanding the early evolution of the Aroideae. Also of particular interest are the phylogenetic positions of the isolated genera Montrichardia, Zantedeschia, and Anchomanes, which remain only moderately supported here.

A target Capture Probe Set Useful for Deep- and Shallow-Level Phylogenetic Studies in Cactaceae

The molecular phylogenies of Cactaceae have enabled us to better understand their systematics, biogeography, and diversification ages. However, most of the phylogenetic relationships within Cactaceae major groups remain unclear, largely due to the lack of an appropriate set of molecular markers to resolve its contentious relationships. Here, we explored the genome and transcriptome assemblies available for Cactaceae and identified putative orthologous regions shared among lineages of the subfamily Cactoideae. Then we developed a probe set, named Cactaceae591, targeting both coding and noncoding nuclear regions for representatives from the subfamilies Pereskioideae, Opuntioideae, and Cactoideae. We also sampled inter- and intraspecific variation to evaluate the potential of this panel to be used in phylogeographic studies. We retrieved on average of 547 orthologous regions per sample. Targeting noncoding nuclear regions showed to be crucial to resolving inter- and intraspecific relationships. Cactaceae591 covers 13 orthologous genes shared with the Angiosperms353 kit and two plastid regions largely used in Cactaceae studies, enabling the phylogenies generated by our panel to be integrated with angiosperm and Cactaceae phylogenies, using these sequences. We highlighted the importance of using coalescent-based species tree approaches on the Cactaceae591 dataset to infer accurate phylogenetic trees in the presence of extensive incomplete lineage sorting in this family.

A New Approach Using Targeted Sequence Capture for Phylogenomic Studies across Cactaceae

Relationships within the major clades of Cactaceae are relatively well known based on DNA sequence data mostly from the chloroplast genome. Nevertheless, some nodes along the backbone of the phylogeny, and especially generic and species-level relationships, remain poorly resolved and are in need of more informative genetic markers. In this study, we propose a new approach to solve the relationships within Cactaceae, applying a targeted sequence capture pipeline. We designed a custom probe set for Cactaceae using MarkerMiner and complemented it with the Angiosperms353 probe set. We then tested both probe sets against 36 different transcriptomes using Hybpiper preferentially retaining phylogenetically informative loci and reconstructed the relationships using RAxML-NG and Astral. Finally, we tested each probe set through sequencing 96 accessions, representing 88 species across Cactaceae. Our preliminary analyses recovered a well-supported phylogeny across Cactaceae with a near identical topology among major clade relationships as that recovered with plastome data. As expected, however, we found incongruences in relationships when comparing our nuclear probe set results to plastome datasets, especially at the generic level. Our results reveal great potential for the combination of Cactaceae-specific and Angiosperm353 probe set application to improve phylogenetic resolution for Cactaceae and for other studies.

A customised target capture sequencing tool for molecular identification of Aloe vera and relatives

Plant molecular identification studies have, until recently, been limited to the use of highly conserved markers from plastid and other organellar genomes, compromising resolution in highly diverse plant clades. Due to their higher evolutionary rates and reduced paralogy, low-copy nuclear genes overcome this limitation but are difficult to sequence with conventional methods and require high-quality input DNA. Aloe vera and its relatives in the Alooideae clade (Asphodelaceae, subfamily Asphodeloideae) are of economic interest for food and health products and have horticultural value. However, pressing conservation issues are increasing the need for a molecular identification tool to regulate the trade. With > 600 species and an origin of ± 15 million years ago, this predominantly African succulent plant clade is a diverse and taxonomically complex group for which low-copy nuclear genes would be desirable for accurate species discrimination. Unfortunately, with an average genome size of 16.76 pg, obtaining high coverage sequencing data for these genes would be prohibitively costly and computationally demanding. We used newly generated transcriptome data to design a customised RNA-bait panel targeting 189 low-copy nuclear genes in Alooideae. We demonstrate its efficacy in obtaining high-coverage sequence data for the target loci on Illumina sequencing platforms, including degraded DNA samples from museum specimens, with considerably improved phylogenetic resolution. This customised target capture sequencing protocol has the potential to confidently indicate phylogenetic relationships of Aloe vera and related species, as well as aid molecular identification applications.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.