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.

A Comprehensive Phylogenomic Platform for Exploring the Angiosperm Tree of Life.

The tree of life is the fundamental biological roadmap for navigating the evolution and properties of life on Earth, and yet remains largely unknown. Even angiosperms (flowering plants) are fraught with data gaps, despite their critical role in sustaining terrestrial life. Today, high-throughput sequencing promises to significantly deepen our understanding of evolutionary relationships. Here, we describe a comprehensive phylogenomic platform for exploring the angiosperm tree of life, comprising a set of open tools and data based on the 353 nuclear genes targeted by the universal Angiosperms353 sequence capture probes. The primary goals of this article are to (i) document our methods, (ii) describe our first data release, and (iii) present a novel open data portal, the Kew Tree of Life Explorer (https://treeoflife.kew.org). We aim to generate novel target sequence capture data for all genera of flowering plants, exploiting natural history collections such as herbarium specimens, and augment it with mined public data. Our first data release, described here, is the most extensive nuclear phylogenomic data set for angiosperms to date, comprising 3099 samples validated by DNA barcode and phylogenetic tests, representing all 64 orders, 404 families (96$\%$) and 2333 genera (17$\%$). A "first pass" angiosperm tree of life was inferred from the data, which totaled 824,878 sequences, 489,086,049 base pairs, and 532,260 alignment columns, for interactive presentation in the Kew Tree of Life Explorer. This species tree was generated using methods that were rigorous, yet tractable at our scale of operation. Despite limitations pertaining to taxon and gene sampling, gene recovery, models of sequence evolution and paralogy, the tree strongly supports existing taxonomy, while challenging numerous hypothesized relationships among orders and placing many genera for the first time. The validated data set, species tree and all intermediates are openly accessible via the Kew Tree of Life Explorer and will be updated as further data become available. This major milestone toward a complete tree of life for all flowering plant species opens doors to a highly integrated future for angiosperm phylogenomics through the systematic sequencing of standardized nuclear markers. Our approach has the potential to serve as a much-needed bridge between the growing movement to sequence the genomes of all life on Earth and the vast phylogenomic potential of the world's natural history collections. [Angiosperms; Angiosperms353; genomics; herbariomics; museomics; nuclear phylogenomics; open access; target sequence capture; tree of life.].

Microsatellites and petal morphology reveal new patterns of admixture in Orchis hybrid zones.

PREMISE: The genetic structure of hybrid zones provides insight into the potential for gene flow to occur between plant taxa. Four closely related European orchid species (Orchis anthropophora, O. militaris, O. purpurea, and O. simia) hybridize when they co-occur. We aimed to characterize patterns of hybridization in O. militaris-O. purpurea, O. purpurea-O. simia, and O. anthropophora-O. simia hybrid zones using molecular and morphological data. METHODS: We used 11 newly isolated nuclear microsatellites to genotype 695 individuals collected from seven hybrid zones and six allopatric parental populations in France. Geometric morphometric analysis was conducted using 15 labellum landmarks to capture the main aspects of petal shape. RESULTS: Backcrossing was asymmetric toward O. militaris in multiple O. militaris-O. purpurea hybrid zones. Hybrids in O. purpurea-O. simia and O. anthropophora-O. simia hybrid zones were largely limited to F1 and F2 generations, but further admixture had occurred. These patterns were reflected in labellum geometric morphometric data, which correlated strongly with nuclear microsatellite data in all three species combinations. CONCLUSIONS: The coexistence of parental and admixed individuals in these Orchis hybrid zones implies they are likely to be tension zones being maintained by a balance between gene flow into the hybrid zone and selection acting against admixed individuals. The pattern of admixture in the three species combinations suggests intrinsic selection acting on the hybrids is weaker in more closely related taxa.

A nuclear phylogenomic study of the angiosperm order Myrtales, exploring the potential and limitations of the universal Angiosperms353 probe set.

PREMISE: To further advance the understanding of the species-rich, economically and ecologically important angiosperm order Myrtales in the rosid clade, comprising nine families, approximately 400 genera and almost 14,000 species occurring on all continents (except Antarctica), we tested the Angiosperms353 probe kit. METHODS: We combined high-throughput sequencing and target enrichment with the Angiosperms353 probe kit to evaluate a sample of 485 species across 305 genera (76% of all genera in the order). RESULTS: Results provide the most comprehensive phylogenetic hypothesis for the order to date. Relationships at all ranks, such as the relationship of the early-diverging families, often reflect previous studies, but gene conflict is evident, and relationships previously found to be uncertain often remain so. Technical considerations for processing HTS data are also discussed. CONCLUSIONS: High-throughput sequencing and the Angiosperms353 probe kit are powerful tools for phylogenomic analysis, but better understanding of the genetic data available is required to identify genes and gene trees that account for likely incomplete lineage sorting and/or hybridization events.

Whole plastomes are not enough: phylogenomic and morphometric exploration at multiple demographic levels of the bee orchid clade Ophrys sect. Sphegodes.

Plastid sequences have long dominated phylogeny reconstruction at all time depths, predicated on a usually untested assumption that they accurately represent the evolutionary histories of phenotypically circumscribed species. We combined detailed in situ morphometrics (124 plants) and whole-plastome sequencing through genome skimming (71 plants) in order to better understand species-level diversity and speciation in arguably the most challenging monophyletic group within the taxonomically controversial, pseudo-copulatory bee orchid genus Ophrys. Using trees and ordinations, we interpreted the data at four nested demographic levels-macrospecies, mesospecies, microspecies, and local population-seeking the optimal level for bona fide species. Neither morphological nor molecular discontinuities are evident at any level below macrospecies, the observed overlap among taxa suggesting that both mesospecies and microspecies reflect arbitrary division of a continuum of variation. Plastomes represent geographic location more strongly than taxonomic assignment and correlate poorly with morphology, suggesting widespread plastid capture and possibly post-glacial expansion from multiple southern refugia. As they are rarely directly involved in the speciation process, plastomes depend on extinction of intermediate lineages to provide phylogenetic signal and so cannot adequately document evolutionary radiations. The popular 'ethological' evolutionary model recognizes as numerous 'ecological species' (microspecies) lineages perceived as actively diverging as a result of density-dependent selection on very few features that immediately dictate extreme pollinator specificity. However, it is assumed rather than demonstrated that the many microspecies are genuinely diverging. We conversely envisage a complex four-dimensional reticulate network of lineages, generated locally and transiently through a wide spectrum of mechanisms, but each unlikely to maintain an independent evolutionary trajectory long enough to genuinely speciate by escaping ongoing gene flow. The frequent but localized microevolution that characterizes the Ophrys sphegodes complex is often convergent and rarely leads to macroevolution. Choosing between the contrasting 'discontinuity' and 'ethology' models will require next-generation sequencing of nuclear genomes plus ordination of corresponding morphometric matrices, seeking the crucial distinction between retained ancestral polymorphism-consistent with lineage divergence-and polymorphisms reflecting gene flow through 'hybridization'-more consistent with lineage convergence.

Phylogeny and biogeography of Myrcia sect. Aguava (Myrtaceae, Myrteae) based on phylogenomic and Sanger data provide evidence for a Cerrado origin and geographically structured clades.

Myrcia is one of the largest exclusively Neotropical angiosperm genera, including ca. 800 species divided into nine sections. Myrcia sect. Aguava is one of most complex sections of Myrcia due to high morphological variation and wide distribution range of some species, including M. guianensis, with distribution throughout South America and a complex taxonomic history. We used complete plastid DNA sequences data generated using next-generation sequencing of 45 terminals, mostly from Myrcia sect. Aguava. These data were combined with five target DNA regions (ITS, psbA-trnH, trnL-trnF, trnQ-rps16, ndhF) of additional terminals to increase taxonomic coverage. Phylogenetic analyses were conducted using a maximum likelihood approach, and divergence times and ancestral range distributions were estimated. Myrcia sect. Aguava is monophyletic and exclusively comprises species with trilocular ovaries but has no relationship with other groups within Myrcia that possess trilocular ovaries. Three main lineages that correspond to geographical distribution are recognized within Myrcia sect. Aguava. Multiple accessions reveal a non-monophyletic Myrcia guianensis and stress the biogeographical structure inside the group. Myrcia sect. Aguava had a probable mid-Miocene origin in the Cerrado, but lineages that persisted there diversified only more recently, when the present-day vegetation started to stabilize. Posterior migrations to Atlantic Forest, Amazon and Caribbean occurred at the end of Miocene, evidencing transitions from open and dry to forested and more humid areas that are less frequent in the Neotropics. Overall, it is observed that related lineages remained in ecologically similar environments. Future perspectives on Myrcia and Myrteae in the phylogenomic era are also discussed.

Factors Affecting Targeted Sequencing of 353 Nuclear Genes From Herbarium Specimens Spanning the Diversity of Angiosperms.

The world's herbaria collectively house millions of diverse plant specimens, including endangered or extinct species and type specimens. Unlocking genetic data from the typically highly degraded DNA obtained from herbarium specimens was difficult until the arrival of high-throughput sequencing approaches, which can be applied to low quantities of severely fragmented DNA. Target enrichment involves using short molecular probes that hybridise and capture genomic regions of interest for high-throughput sequencing. In this study on herbariomics, we used this targeted sequencing approach and the Angiosperms353 universal probe set to recover up to 351 nuclear genes from 435 herbarium specimens that are up to 204 years old and span the breadth of angiosperm diversity. We show that on average 207 genes were successfully retrieved from herbarium specimens, although the mean number of genes retrieved and target enrichment efficiency is significantly higher for silica gel-dried specimens. Forty-seven target nuclear genes were recovered from a herbarium specimen of the critically endangered St Helena boxwood, Mellissia begoniifolia, collected in 1815. Herbarium specimens yield significantly less high-molecular-weight DNA than silica gel-dried specimens, and genomic DNA quality declines with sample age, which is negatively correlated with target enrichment efficiency. Climate, taxon-specific traits, and collection strategies additionally impact target sequence recovery. We also detected taxonomic bias in targeted sequencing outcomes for the 10 most numerous angiosperm families that were investigated in depth. We recommend that (1) for species distributed in wet tropical climates, silica gel-dried specimens should be used preferentially; (2) for species distributed in seasonally dry tropical climates, herbarium and silica gel-dried specimens yield similar results, and either collection can be used; (3) taxon-specific traits should be explored and established for effective optimisation of taxon-specific studies using herbarium specimens; (4) all herbarium sheets should, in future, be annotated with details of the preservation method used; (5) long-term storage of herbarium specimens should be in stable, low-humidity, and low-temperature environments; and (6) targeted sequencing with universal probes, such as Angiosperms353, should be investigated closely as a new approach for DNA barcoding that will ensure better exploitation of herbarium specimens than traditional Sanger sequencing approaches.

Conservation of the Threatened Species, Pulsatilla vulgaris Mill. (Pasqueflower), is Aided by Reproductive System and Polyploidy.

Population loss due to habitat disturbance is a major concern in biodiversity conservation. Here we investigate the genetic causes of the demographic decline observed in English populations of Pulsatilla vulgaris and the consequences for conservation. Using 10 nuclear microsatellite markers, we compare genetic variation in wild populations with restored and seed-regenerated populations (674 samples). Emergence of genetic structure and loss of allelic variation in natural populations are not as evident as expected from demographic trends. Restored populations show genetic variation comparable to their source populations and, in general, to the wild ones. Genetic homogeneity is observed in regeneration trials, although some alleles not captured in source populations are detected. We infer that polyploidy, longevity, and clonal reproduction have provided P. vulgaris with the standing genetic variation necessary to make the species resilient to the effects of demographic decline, suggesting that the use of multiple sources for reintroduction may be beneficial to mimic natural gene flow and the availability of multiple allele copies typical of polyploid species.

A Universal Probe Set for Targeted Sequencing of 353 Nuclear Genes from Any Flowering Plant Designed Using k-Medoids Clustering.

Sequencing of target-enriched libraries is an efficient and cost-effective method for obtaining DNA sequence data from hundreds of nuclear loci for phylogeny reconstruction. Much of the cost of developing targeted sequencing approaches is associated with the generation of preliminary data needed for the identification of orthologous loci for probe design. In plants, identifying orthologous loci has proven difficult due to a large number of whole-genome duplication events, especially in the angiosperms (flowering plants). We used multiple sequence alignments from over 600 angiosperms for 353 putatively single-copy protein-coding genes identified by the One Thousand Plant Transcriptomes Initiative to design a set of targeted sequencing probes for phylogenetic studies of any angiosperm group. To maximize the phylogenetic potential of the probes, while minimizing the cost of production, we introduce a k-medoids clustering approach to identify the minimum number of sequences necessary to represent each coding sequence in the final probe set. Using this method, 5-15 representative sequences were selected per orthologous locus, representing the sequence diversity of angiosperms more efficiently than if probes were designed using available sequenced genomes alone. To test our approximately 80,000 probes, we hybridized libraries from 42 species spanning all higher-order groups of angiosperms, with a focus on taxa not present in the sequence alignments used to design the probes. Out of a possible 353 coding sequences, we recovered an average of 283 per species and at least 100 in all species. Differences among taxa in sequence recovery could not be explained by relatedness to the representative taxa selected for probe design, suggesting that there is no phylogenetic bias in the probe set. Our probe set, which targeted 260 kbp of coding sequence, achieved a median recovery of 137 kbp per taxon in coding regions, a maximum recovery of 250 kbp, and an additional median of 212 kbp per taxon in flanking non-coding regions across all species. These results suggest that the Angiosperms353 probe set described here is effective for any group of flowering plants and would be useful for phylogenetic studies from the species level to higher-order groups, including the entire angiosperm clade itself.

Flower preferences and pollen transport networks for cavity-nesting solitary bees: Implications for the design of agri-environment schemes.

Floral foraging resources are valuable for pollinator conservation on farmland, and their provision is encouraged by agri-environment schemes in many countries. Across Europe, wildflower seed mixtures are widely sown on farmland to encourage pollinators, but the extent to which key pollinator groups such as solitary bees exploit and benefit from these resources is unclear. We used high-throughput sequencing of 164 pollen samples extracted from the brood cells of six common cavity-nesting solitary bee species (Osmia bicornis, Osmia caerulescens, Megachile versicolor, Megachile ligniseca, Megachile centuncularis and Hylaeus confusus) which are widely distributed across the UK and Europe. We documented their pollen use across 19 farms in southern England, UK, revealing their forage plants and examining the structure of their pollen transport networks. Of the 32 plant species included currently in sown wildflower mixes, 15 were recorded as present within close foraging range of the bees on the study farms, but only Ranunculus acris L. was identified within the pollen samples. Rosa canina L. was the most commonly found of the 23 plant species identified in the pollen samples, suggesting that, in addition to providing a nesting resource for Megachile leafcutter bees, it may be an important forage plant for these species. Higher levels of connectance and nestedness were characteristic of pollen transport networks on farms with abundant floral resources, which may increase resilience to species loss. Our data suggest that plant species promoted currently by agri-environment schemes are not optimal for solitary bee foraging. If a diverse community of pollinators is to be supported on UK and European farmland, additional species such as R. canina should be encouraged to meet the foraging requirements of solitary bees.

Isolation and characterization of microsatellite loci in Sorbus porrigentiformis and cross-amplification in S. aria and S. rupicola (Rosaceae).

PREMISE OF THE STUDY: Southwestern Britain is an emblematic hotspot of polyploid diversity of whitebeams (Sorbus aria agg.; Rosaceae) with ca. 30 polyploid endemic species. The tetraploid S. porrigentiformis is postulated as one of the parents of most of these endemics, along with the sexual diploid S. aria s. str. and the tetraploid S. rupicola. METHODS AND RESULTS: We isolated 16 nuclear microsatellite loci from S. porrigentiformis and characterized them on 45 trees representing the three putative parental species. Eleven loci were polymorphic, and eight of them exhibited species-specific alleles. Allele numbers ranged from one to 11, and observed heterozygosity ranged from 0.40 to 1.00. The intraspecific levels of variation were very low, in agreement with the facultative apomictic reproduction hypothesized for this species. CONCLUSIONS: The species-specific alleles will be useful for tracing the origin of the narrowly distributed Sorbus taxa. In addition, the assessment of diversity levels will help design a conservation strategy for the polyploid complex.

Key processes for Cheirolophus (Asteraceae) diversification on oceanic islands inferred from AFLP data.

The radiation of the genus Cheirolophus (Asteraceae) in Macaronesia constitutes a spectacular case of rapid diversification on oceanic islands. Twenty species - nine of them included in the IUCN Red List of Threatened Species - have been described to date inhabiting the Madeiran and Canarian archipelagos. A previous phylogenetic study revealed that the diversification of Cheirolophus in Macaronesia started less than 2 Ma. As a result of such an explosive speciation process, limited phylogenetic resolution was reported, mainly due to the low variability of the employed molecular markers. In the present study, we used highly polymorphic AFLP markers to i) evaluate species' boundaries, ii) infer their evolutionary relationships and iii) investigate the patterns of genetic diversity in relation to the potential processes likely involved in the radiation of Cheirolophus. One hundred and seventy-two individuals representing all Macaronesian Cheirolophus species were analysed using 249 AFLP loci. Our results suggest that geographic isolation played an important role in this radiation process. This was likely driven by the combination of poor gene flow capacity and a good ability for sporadic long-distance colonisations. In addition, we also found some traces of introgression and incipient ecological adaptation, which could have further enhanced the extraordinary diversification of Cheirolophus in Macaronesia. Last, we hypothesize that current threat categories assigned to Macaronesian Cheirolophus species do not reflect their respective evolutionary relevance, so future evaluations of their conservation status should take into account the results presented here.

How to open the treasure chest? Optimising DNA extraction from herbarium specimens.

Herbarium collections are potentially an enormous resource for DNA studies, but the use of herbarium specimens in molecular studies has thus far been slowed down by difficulty in obtaining amplifiable DNA. Here we compare a set of commercially available DNA extraction protocols and their performance in terms of DNA purity and yield, and PCR amplification success as measured by using three differentially sized markers, the rbcL barcoding marker (cpDNA), the LEAFY exon 3 (nrDNA), and the trnL((UAA)) P6 loop (cpDNA). Results reveal large differences between extraction methods, where DNA purity rather than yield is shown to be strongly correlated with PCR success. Amplicon size shows similarly strong correlation with PCR success, with the shortest fragment showing the highest success rate (78%, P6 loop, 10-143 base pairs (bp)) and the largest fragment the lowest success (10%, rbcL, 670 bp). The effect of specimen preparation method on PCR success was also tested. Results show that drying method strongly affects PCR success, especially the availability of fragments longer than 250 bp, where longer fragments are more available for PCR amplification in air dried material compared to alcohol dried specimens. Results from our study indicate that projects relying on poor-quality starting material such as herbarium or scat samples should focus on extracting pure DNA and aim to amplify short target regions (<200-300 bp) in order to maximise outcomes. Development of shorter barcoding regions, or mini-barcodes within existing ones should be of high importance as only a few options are currently available; this is particularly important if we hope to incorporate the millions of herbarium samples available into barcoding initiatives and other molecular studies.

Challenges in the DNA barcoding of plant material.

DNA barcoding, using a short gene sequence from a standardized region of the genome, is a species identification tool which would not only aid species discovery but would also have applications ranging from large-scale biodiversity surveys through to identification of a single fragment of material in forensic contexts. To fulfill this vision a universal, relatively cheap, scalable system needs to be in place. The mitochondrial locus being used for many animal groups and algae is not suitable for use in land plants, and an appropriate alternative is needed.Progress has been made in the selection of two alternative regions for plant DNA barcoding. There are however many challenges in finding a solution that fulfills all the requirements of a successful, universally applicable barcode, and in the short term a pragmatic solution that achieves as much as possible and has payoffs in most areas has been chosen. Research continues in areas ranging from the technicalities of sequencing the regions to data analysis and the potential improvements that may result from the developing technology and data analysis systems.The ultimate success of DNA barcoding as a plant identification tool for all occasions depends on the building of a reference database and it fulfilling the requirements of potential users such that they are able to achieve valid results through its use, that would be more time consuming and costly, and less reliable using other techniques.

Towards writing the encyclopedia of life: an introduction to DNA barcoding.

An international consortium of major natural history museums, herbaria and other organizations has launched an ambitious project, the 'Barcode of Life Initiative', to promote a process enabling the rapid and inexpensive identification of the estimated 10 million species on Earth. DNA barcoding is a diagnostic technique in which short DNA sequence(s) can be used for species identification. The first international scientific conference on Barcoding of Life was held at the Natural History Museum in London in February 2005, and here we review the scientific challenges discussed during this conference and in previous publications. Although still controversial, the scientific benefits of DNA barcoding include: (i) enabling species identification, including any life stage or fragment, (ii) facilitating species discoveries based on cluster analyses of gene sequences (e.g. cox1 = CO1, in animals), (iii) promoting development of handheld DNA sequencing technology that can be applied in the field for biodiversity inventories and (iv) providing insight into the diversity of life.

The effects of nuclear DNA content (C-value) on the quality and utility of AFLP fingerprints.

BACKGROUND AND AIMS: Nuclear DNA content (C-value) varies approximately 1000-fold across the angiosperms, and this variation has been reported to have an effect on the quality of AFLP fingerprints. Various methods have been proposed for circumventing the problems associated with small and large genomes. Here we investigate the range of nuclear DNA contents across which the standard AFLP protocol can be used. METHODS: AFLP fingerprinting was conducted on an automated platform using the standard protocol (with 3 + 3 selective bases) in which DNA fragments are visualized as bands. Species with nuclear DNA contents ranging from 1C = 0.2 to 32.35 pg were included, and the total number of bands and the number of polymorphic bands were counted. For the species with the smallest C-value (Bixa orellana) and for one of the species with a large C-value (Damasonium alisma), alternative protocols using 2 + 3 and 3 + 4 selective bases, respectively, were also used. KEY RESULTS: Acceptable AFLP traces were obtained using the standard protocol with 1C-values of 0.30-8.43 pg. Below this range, the quality was improved by using 2 + 3 selective bases. Above this range, the traces were generally characterized by a few strongly amplifying bands and noisy baselines. Damasonium alisma, however, gave more even traces, probably due to it being a tetraploid. CONCLUSIONS: We propose that for known polyploids, genome size is a more useful indicator than the 1C-value in deciding which AFLP protocol to use. Thus, knowledge of ploidy (allowing estimation of genome size) and C-value are both important. For small genomes, the number of interpretable bands can be increased by decreasing the number of selective bases. For larger genomes, increasing the number of bases does not necessarily decrease the number of bands as predicted. The presence of a small number of strongly amplifying bands is likely to be linked to the presence of repetitive DNA sequences in high copy number in taxa with large genomes.