Phylogenomics of the Major Tropical Plant Family Annonaceae Using Targeted Enrichment of Nuclear Genes

Dispersal from Africa to the Neotropics was followed by multiple transitions across Neotropical biomes facilitated by frugivores

BACKGROUND AND AIM

Plant disjunctions have fascinated biogeographers and ecologists for a long time. We use tribe Bocageeae (Annonaceae), a predominantly Neotropical plant group distributed across several present-day Neotropical biomes and with an African-American disjunction, to investigate long-distance dispersal mediated by frugivorous animals at both intercontinental and intracontinental scales.

METHODS

We reconstructed a species-level phylogeny of tribe Bocageeae with a dataset composed of 116 nuclear markers. We sampled 70 % of Bocageeae species, covering their geographical range and representing all eight genera. We estimated divergence times using BEAST, inferred ancestral range distributions and reconstructed ancestral states for fruit traits related to long-distance dispersal in a Bayesian framework.

KEY RESULTS

The ancestral Bocageeae date to the Early Eocene and were inferred to occur in Africa and proto-Amazonia. Their ancestral fruits were large and dehiscent. The first lineage split gave rise to an exclusively Neotropical clade during the Middle Eocene, in proto-Amazonia. Range exchange between the Amazon and the Atlantic Forest occurred at least once during the Miocene, and from Amazonia to Central America and Mexico during the Early Miocene. Transitions in different sets of fruit morphologies were inferred to be related to dispersal events across South American regions/biomes.

CONCLUSIONS

In Bocageeae, mammals might have been responsible for long-distance dispersal through the Boreotropics. In the Neotropics, proto-Amazonia is proposed to be the source for dispersal to other tropical American biomes. Long-distance dispersal might have happened via a wide range of dispersal guilds, depending on frugivore radiations, diversity and abundance in particular time periods and places. Hence, inter- and intracontinental dispersal might not rely on a single dispersal syndrome or guild, but more on the availability of frugivorous lineages for seed dispersal.

Sequential diversification with Miocene extinction and Pliocene speciation linked to mountain uplift explains the diversity of the African rain forest clade Monodoreae (Annonaceae)

BACKGROUND AND AIMS

Throughout the Cenozoic, Africa underwent several climatic and geological changes impacting the evolution of tropical rain forests (TRFs). African TRFs are thought to have extended from east to west in a 'pan-African' TRF, followed by several events of fragmentation during drier climate periods. During the Miocene, climate cooling and mountain uplift led to the aridification of tropical Africa and open habitats expanded at the expense of TRFs, which probably experienced local extinctions. However, in plants, these drivers were previously inferred using limited taxonomic and molecular data. Here, we tested the impact of climate and geological changes on diversification within the diverse clade Monodoreae (Annonaceae) composed of 90 tree species restricted to African TRFs.

METHODS

We reconstructed a near-complete phylogenetic tree, based on 32 nuclear genes, and dated using relaxed clocks and fossil calibrations in a Bayesian framework. We inferred the biogeographical history and the diversification dynamics of the clade using multiple birth-death models.

KEY RESULTS

Monodoreae originated in East African TRFs ~25 million years ago (Ma) and expanded toward Central Africa during the Miocene. We inferred range contractions during the middle Miocene and document important connections between East and West African TRFs after 15-13 Ma. Our results indicated a sudden extinction event during the late Miocene, followed by an increase in speciation rates. Birth-death models suggested that African elevation change (orogeny) is positively linked to speciation in this clade.

CONCLUSION

East Africa is inferred as an important source of Monodoreae species, and possibly for African plant diversity in general. Our results support a 'sequential scenario of diversification' in which increased aridification triggered extinction of TRF species in Monodoreae. This was quickly followed by fragmentation of rain forests, subsequently enhancing lagged speciation resulting from vicariance and improved climate conditions. In contrast to previous ideas, the uplift of East Africa is shown to have played a positive role in Monodoreae diversification.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

An updated and extended version of the Melastomataceae probe set for target capture

Premise

A probe set was previously designed to target 384 nuclear loci in the Melastomataceae family; however, when trying to use it, we encountered several practical and conceptual problems, such as the presence of sequences in reverse complement, intronic regions with stop codons, and other issues. This raised concerns regarding the use of this probe set for sequence recovery in Melastomataceae.

Methods

In order to correct these issues, we cleaned the Melastomataceae probe set, extended it with additional sequences, and compared its performance with the original version.

Results

The final probe set targets 396 putative nuclear loci represented by 6009 template sequences. The probe set has been made available, along with details on the cleaning process, for reproducibility. We show that the new probe set performs better than the original version in terms of sequence recovery.

Discussion

This updated, extended, and cleaned probe set will improve the availability of phylogenomic resources across the Melastomataceae family. It is fully compatible with sequence recovery and extraction pipelines. The cleaning process can also be applied to any plant-targeting probe set that would need to be cleaned or updated if new genomic resources for the targeted taxa become available.

Phylogenomic inference of the African tribe Monodoreae (Annonaceae) and taxonomic revision of Dennettia, Uvariodendron and Uvariopsis

Monodoreae (Annonaceae) is a tribe composed of 11 genera and 90 species restricted to the tropical African rain forests. All the genera are taxonomically well circumscribed except the species rich genera Uvariodendron and Uvariopsis which lack a recent taxonomic revision. Here, we used a robust phylogenomic approach, including all the 90 currently accepted species, with several specimens per species, and based on more than 300 Annonaceae-specific nuclear genes, to infer the phylogenetic tree of the Monodoreae and test the limits between the genera and species. We recover all the genera as monophyletic, except the genus Uvariopsis for which the species Uvariopsistripetala falls outside this clade. We thus reinstate the monotypic genus Dennettia for its single species Dennettiatripetala. We also erect a new tribe, Ophrypetaleae trib. nov., to accommodate the genera Ophrypetalum and Sanrafaelia, as we recover them excluded from the Monodoreae tribe with good support. Below the genus level, the genera Isolona, Monodora, Uvariastrum, Uvariodendron and Uvariopsis show weakly supported nodes and phylogenetic conflicts, suggesting that population level processes of evolution might occur in these clades. Our results also support, at the molecular level, the description of several new species of Uvariodendron and Uvariopsis, as well as several new synonymies. Finally, we present a taxonomic revision of the genera Dennettia, Uvariodendron and Uvariopsis, which contain one, 18 and 17 species respectively. We provide a key to the 11 genera of the Monodoraeae and describe four new species to science: Uvariodendronkimbozaense Dagallier & Couvreur, sp. nov., Uvariodendronmossambicense Robson ex Dagallier & Couvreur, sp. nov., Uvariodendronpilosicarpum Dagallier & Couvreur, sp. nov. and Uvariopsisoligocarpa Dagallier & Couvreur, sp. nov., and provide provisional descriptions of three putatively new species. We also present lectotypifications and nomenclatural changes implying synonymies and new combinations (Uvariodendroncitriodorum (Le Thomas) Dagallier & Couvreur, comb. et stat. nov., Uvariodendronfuscumvar.magnificum (Verdc.) Dagallier & Couvreur, comb. et stat. nov., Uvariopsiscongensisvar.angustifolia Dagallier & Couvreur, var. nov., Uvariopsisguineensisvar.globiflora (Keay) Dagallier & Couvreur, comb. et stat. nov., and Uvariopsissolheidiivar.letestui (Pellegr.) Dagallier & Couvreur, comb. et stat. nov.).

Desmopsisterriflora, an extraordinary new species of Annonaceae with flagelliflory

Flagelliflory refers to the production of inflorescences exclusively on long, whip-like branches which emerge from the main trunk and extend along the ground or below it. It is the rarest type of cauliflory and only a few cases have been reported in the world. Here, a new species of Annonaceae with flagelliflory is described and illustrated. The phylogenetic relationships of the new species were inferred using a hybrid-capture phylogenomic approach and we present some notes on its reproductive ecology and pollen characteristics. The new species, namely Desmopsisterriflorasp. nov., is part of a clade composed of Mexican species of Stenanona with long, awned petals. Desmopsisterriflora is distinguished by its flageliflorous inflorescences, basely fused sepals, thick red petals, reduced number of ovules per carpel, pollen grains with a weakly rugulate to fossulate exine ornamentation, and its globose, apiculate fruits with a woody testa. The morphological characteristics of the flagella suggest that these are specialized branches rather than inflorescences, and the absence of ramiflory implies an exclusively reproductive function. The flowers are infrequently visited by insects, their potential pollinators being flies and ants.

New insights into polyploid evolution and dynamic nature of Ludwigia section Isnardia (Onagraceae)

BACKGROUND

While polyploids are common in plants, the evolutionary history and natural dynamics of most polyploid groups are still unclear. Owing to plentiful earlier systematic studies, Ludwigia sect. Isnardia (comprising 22 wetland taxa) is an ideal allopolyploid complex to investigate polyploid evolution and natural dynamics within and among taxa. With a considerable sampling, we concentrated on revisiting earlier phylogenies of Isnardia, reevaluating the earlier estimated age of the most recent common ancestor (TMRCA), exploring the correlation between infraspecific genetic diversity and ploidy levels, and inspecting interspecific gene flows among taxa.

RESULTS

Phylogenetic trees and network concurred with earlier phylogenies and hypothesized genomes by incorporating 192 atpB-rbcL and ITS sequences representing 91% of Isnardia taxa. Moreover, we detected three multi-origin taxa. Our findings on L. repens and L. sphaerocarpa were consistent with earlier studies; L. arcuata was reported as a multi-origin taxon here, and an additional evolutionary scenario of L. sphaerocarpa was uncovered, both for the first time. Furthermore, estimated Isnardia TMRCA ages based on our data (5.9 or 8.9 million years ago) are in accordance with earlier estimates, although younger than fossil dates (Middle Miocene). Surprisingly, infraspecific genetic variations of Isnardia taxa did not increase with ploidy levels as anticipated from many other polyploid groups. In addition, the exuberant, low, and asymmetrical gene flows among Isnardia taxa indicated that the reproductive barriers may be weakened owing to allopolyploidization, which has rarely been reported.

CONCLUSIONS

The present research gives new perceptions of the reticulate evolution and dynamic nature of Isnardia and points to gaps in current knowledge about allopolyploid evolution.

Flora of Cameroon - Annonaceae Vol 45

Annonaceae is a major pantropical family with 113 genera and about 2550 species. Cameroon is one of the most biodiverse countries in Africa but its flora remains incompletely known. In this volume of the Flora of Cameroon, we describe 166 native taxa representing 163 species in 28 native genera within the family Annonaceae. A total of 22 species (about 13%) are endemic to the country. We provide keys to all native genera, species, and infraspecific taxa. For each species a detailed morphological description and a map of its distribution in Cameroon are provided. Distribution maps and diversity analyses are based on a taxonomically verified database of 2073 collections. Across Africa, Cameroon is a center of diversity for Annonaceae harboring one of the highest numbers of species and genera. For example, Cameroon harbors the highest number of African species for the only pantropical genus of Annonaceae, Xylopia. Annonaceae are found across all 10 administrative regions of Cameroon but diversity is concentrated within the tropical rain forest areas situated in the south and South-West. The areas around Bipindi and Mount Cameroon show the highest levels of diversity, but this is correlated with collection effort. Line drawings and/or photographs accompany most species. One species new to science Uvariopsisetugeana Dagallier & Couvreur sp. nov. is described. We also undertake a number of nomenclatural changes such as lectotypifications, six new synonymies and two new combinations (Uvariaanisotricha (Le Thomas) Couvreur, comb. nov.; Uvariodendronfuscumvar.giganteum (Engl.) Dagallier & Couvreur, comb. nov.).

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.

Oxygenated Cyclohexene Derivatives from the Stem and Root Barks of Uvaria pandensis

Five new cyclohexene derivatives, dipandensin A and B (1 and 2) and pandensenols A-C (3-5), and 16 known secondary metabolites (6-21) were isolated from the methanol-soluble extracts of the stem and root barks of Uvaria pandensis. The structures were characterized by NMR spectroscopic and mass spectrometric analyses, and that of 6-methoxyzeylenol (6) was further confirmed by single-crystal X-ray crystallography, which also established its absolute configuration. The isolated metabolites were evaluated for antibacterial activity against the Gram-positive bacteria Bacillus subtilis and Staphylococcus epidermidis and the Gram-negative bacteria Enterococcus raffinosus, Escherichia coli, Paraburkholderia caledonica, Pectobacterium carotovorum, and Pseudomonas putida, as well as for cytotoxicity against the MCF-7 human breast cancer cell line. A mixture of uvaretin (20) and isouvaretin (21) exhibited significant antibacterial activity against B. subtilis (EC50 8.7 μM) and S. epidermidis (IC50 7.9 μM). (8'α,9'β-Dihydroxy)-3-farnesylindole (12) showed strong inhibitory activity (EC50 9.8 μM) against B. subtilis, comparable to the clinical reference ampicillin (EC50 17.9 μM). None of the compounds showed relevant cytotoxicity against the MCF-7 human breast cancer cell line.

Diversity of Insect Flower Visitors of Xylopia aromatica (Magnoliales, Annonaceae) in a Brazilian Savanna

Small beetles are important pollinators of Annonaceae whose flower chambers are small and have diurnal and/or nocturnal anthesis. The pollinators of these flowers belong to the families Nitidulidae, Staphylinidae, Chrysomelidae, and Curculionidae. In this study, the first conducted in the Cerrado of Chapada dos Guimarães, Mato Grosso, Brazil, the behavior of the insect flower visitors of Xylopia aromatica was observed, in both the field and the laboratory. The chambers of 253 flowers were collected from 11 plants, and the biological aspects of their visitors were analyzed quantitatively and qualitatively. The most abundant visitors were thrips and beetles. Coleoptera was represented by four morphospecies occurring frequently in the floral chambers (>70% of individuals). Among beetles, one species belonged to Nitidulidae (Cillaeinae, Conotelus sp. 1) and two belonged to Staphylinidae (Aleocharinae sp. 1 and Aleocharinae sp. 2). These three morphospecies of small elongate beetles have setae where pollen may adhere. In addition, they were present on both male and female phases of the flowers, indicating potential cross-pollination. In the study area, X. aromatica possesses mixed pollination promoted by Thysanoptera and small Nitidulidae and Staphylinidae beetles. This study brings the first record of Lamprosomatinae (Chrysomelidae) and, especially, of Conotelus (Nitidulidae) in the flower chambers of X. aromatica, with new information on behavior of floral visitors coupled with their morphological traits that may promote cross-pollination in this plant species.

Admixture may be extensive among hyperdominant Amazon rainforest tree species

Admixture is a mechanism by which species of long-lived plants may acquire novel alleles. However, the potential role of admixture in the origin and maintenance of tropical plant diversity is unclear. We ask whether admixture occurs in an ecologically important clade of Eschweilera (Parvifolia clade, Lecythidaceae), which includes some of the most widespread and abundant tree species in Amazonian forests. Using target capture sequencing, we conducted a detailed phylogenomic investigation of 33 species in the Parvifolia clade and investigated specific hypotheses of admixture within a robust phylogenetic framework. We found strong evidence of admixture among three ecologically dominant species, E. coriacea, E. wachenheimii and E. parviflora, but a lack of evidence for admixture among other lineages. Accepted species were largely distinguishable from one another, as was geographic structure within species. We show that hybridization may play a role in the evolution of the most widespread and ecologically variable Amazonian tree species. While admixture occurs among some species of Eschweilera, it has not led to widespread erosion of most species' genetic or morphological identities. Therefore, current morphological based species circumscriptions appear to provide a useful characterization of the clade's lineage diversity.

Phylogenomics of Salvia L. subgenus Calosphace (Lamiaceae)

The evolutionary relationships of Salvia have been difficult to estimate. In this study, we used the Next Generation Sequencing method Hyb-Seq to evaluate relationships among 90 Lamiaceae samples, including representatives of Mentheae, Ocimeae, Salvia subgenera Audibertia, Leonia, Salvia, and 69 species of subgenus Calosphace, representing 32 of Epling's sections. A bait set was designed in MarkerMiner using available transcriptome data to enrich 119 variable nuclear loci. Nuclear and chloroplast loci were assembled with hybphylomaker (HPM), followed by coalescent approach analyses for nuclear data (ASTRAL, BEAST) and a concatenated Maximum Likelihood analysis of chloroplast loci. The HPM assembly had an average of 1,314,368 mapped reads for the sample and 527 putative exons. Phylogenetic inferences resolved strongly supported relationships for the deep-level nodes, agreeing with previous hypotheses which assumed that subgenus Audibertia is sister to subgenus Calosphace. Within subgenus Calosphace, we recovered eight monophyletic sections sensu Epling, Cardinalis, Hastatae, Incarnatae, and Uricae in all the analyses (nDNA and cpDNA), Biflorae, Lavanduloideae, and Sigmoideae in nuclear analyses (ASTRAL, BEAST) and Curtiflorae in ASTRAL trees. Network analysis supports deep node relationships, some of the main clades, and recovers reticulation within the core Calosphace. The chloroplast phylogeny resolved deep nodes and four monophyletic Calosphace sections. Placement of S. axillaris is distinct in nuclear evidence and chloroplast, as sister to the rest of the S. subg. Calosphace in chloroplast and a clade with "Hastatae clade" sister to the rest of the subgenus in nuclear evidence. We also tested the monophyly of S. hispanica, S. polystachia, S. purpurea, and S. tiliifolia, including two samples of each, and found that S. hispanica and S. purpurea are monophyletic. Our baits can be used in future studies of Lamiaceae phylogeny to estimate relationships between genera and among species. In this study, we presented a Hyb-Seq phylogeny for complex, recently diverged Salvia, which could be implemented in other Lamiaceae.

Phylogenomic relationships and historical biogeography in the South American vegetable ivory palms (Phytelepheae)

The palm tribe Phytelepheae form a clade of three genera and eight species whose phylogenetic relationships and historical biogeography are not fully understood. Based on morphological similarities and phylogenetic relatedness, it has been suggested that Phytelephas seemannii and Phytelephas schottii are synonyms of Phytelephas macrocarpa, implying the existence of only six species within the Phytelepheae. In addition, uncertainty in their phylogenetic relationships in turn results in blurred biogeographic history. We inferred the phylogenomic relationships in the Phytelepheae by target-capturing 176 nuclear genes and estimated divergence times by using four fossils for time calibration. We lastly explored the biogeographic history of the tribe by inferring its ancestral range evolution. Our phylogenomic trees showed that P. seemannii and P. schottii are not closely related with P. macrocarpa, and therefore, support the existence of eight species in the Phytelepheae. The ancestor of the tribe was widely-distributed in the Chocó, Magdalena, and Amazonia during the Miocene at 19.25 Ma. Early diversification in Phytelephas at 5.27 Ma could have occurred by trans-Andean vicariance after the western Andes uplifted rapidly at ∼ 10 Ma. Our results show the utility of phylogenomic approaches to shed light on species relationships and their biogeographic history.

Coleopteran Pollinators of Annonaceae in the Brazilian Cerrado—A Review

The evolutionary trajectories of insects and angiosperms appear to be intimately interconnected. Increases in the diversity of phytophagous beetles and angiosperms co-occur in the Mesozoic fossil record, and there is fossil evidence of pollinivory and pollination by insects, both in flowering plants and in gymnosperms. The oldest records of angiosperm pollination indicate flies as pollen vectors. A basal group of angiosperms, the order Magnoliales, has retained plesiomorphic characters such as dozens of pistils and stamens spiraling around the receptacle. In a family of this order, Annonaceae, over 90% of species are pollinated by beetles. In many Annonaceae species, flowers display wide spaces, referred to as floral chambers, where beetles can find shelter from weather conditions and predators, food in the form of pollen and tissues, and a mating site. Two basic types of floral chambers can be distinguished: small chambers visited by small beetles (Nitidulidae, Staphylinidae, Chrysomelidae, and Curculionidae) with diurnal and/or nocturnal activity and large and thermogenic floral chambers visited by beetles of the tribe Cyclocephalini (Scarabaeoidea, Melolonthidae). In the latter case, the heat that the flowers produce may serve as a resource for the beetles that visit them, resulting in smaller endothermy costs for the scarabs. This study reviewed the literature including PhD and MSc theses on cantharophilous Annonaceae in the Cerrado. In this biome, both types of associations are found, although cantharophilous Annonaceae represent a small portion of the plant species (<5%). Cantharophilous Annonaceae in the Cerrado share attributes according to the beetles that pollinate them: species pollinated by small beetles, for instance, may flower throughout the year, whereas Annonaceae pollinated by Cyclocephalini normally flower in the beginning of the rainy season (October/November), in synchrony with the phenological patterns of their pollinators. Cantharophilous Annonaceae flowers, regardless of their size, tend to have light colors and sweet and fruity odors. In addition to the lack of studies on the attraction of beetles by these floral characters, the taxonomic composition of the beetles that pollinate Annonaceae in the Cerrado is poorly known. This review attempts to discuss, in light of what has already been published, potential fields of investigation concerning pollinating beetles’ behavior and evolution.

The best of both worlds: Combining lineage-specific and universal bait sets in target-enrichment hybridization reactions

PREMISE

Researchers adopting target-enrichment approaches often struggle with the decision of whether to use universal or lineage-specific probe sets. To circumvent this quandary, we investigate the efficacy of a simultaneous enrichment by combining universal probes and lineage-specific probes in a single hybridization reaction, to benefit from the qualities of both probe sets with little added cost or effort.

METHODS AND RESULTS

Using 26 Brassicaceae libraries and standard enrichment protocols, we compare results from three independent data sets. A large average fraction of reads mapping to the Angiosperms353 (24-31%) and Brassicaceae (35-59%) targets resulted in a sizable reconstruction of loci for each target set (x̄ ≥ 70%).

CONCLUSIONS

High levels of enrichment and locus reconstruction for the two target sets demonstrate that the sampling of genomic regions can be easily extended through the combination of probe sets in single enrichment reactions. We hope that these findings will facilitate the production of expanded data sets that answer individual research questions and simultaneously allow wider applications by the research community as a whole.

New targets acquired: Improving locus recovery from the Angiosperms353 probe set

PREMISE

Universal target enrichment kits maximize utility across wide evolutionary breadth while minimizing the number of baits required to create a cost-efficient kit. The Angiosperms353 kit has been successfully used to capture loci throughout the angiosperms, but the default target reference file includes sequence information from only 6-18 taxa per locus. Consequently, reads sequenced from on-target DNA molecules may fail to map to references, resulting in fewer on-target reads for assembly, and reducing locus recovery.

METHODS

We expanded the Angiosperms353 target file, incorporating sequences from 566 transcriptomes to produce a 'mega353' target file, with each locus represented by 17-373 taxa. This mega353 file is a drop-in replacement for the original Angiosperms353 file in HybPiper analyses. We provide tools to subsample the file based on user-selected taxon groups, and to incorporate other transcriptome or protein-coding gene data sets.

RESULTS

Compared to the default Angiosperms353 file, the mega353 file increased the percentage of on-target reads by an average of 32%, increased locus recovery at 75% length by 49%, and increased the total length of the concatenated loci by 29%.

DISCUSSION

Increasing the phylogenetic density of the target reference file results in improved recovery of target capture loci. The mega353 file and associated scripts are available at: https://github.com/chrisjackson-pellicle/NewTargets.

Chromosome-level reference genome of the soursop (Annona muricata): A new resource for Magnoliid research and tropical pomology

The flowering plant family Annonaceae includes important commercially grown tropical crops, but development of promising species is hindered by a lack of genomic resources to build breeding programs. Annonaceae are part of the magnoliids, an ancient lineage of angiosperms for which evolutionary relationships with other major clades remain unclear. To provide resources to breeders and evolutionary researchers, we report a chromosome-level genome assembly of the soursop (Annona muricata). We assembled the genome using 444.32 Gb of DNA sequences (676× sequencing depth) from PacBio and Illumina short-reads, in combination with 10× Genomics and Bionano data (v1). A total of 949 scaffolds were assembled to a final size of 656.77 Mb, with a scaffold N50 of 3.43 Mb (v1), and then further improved to seven pseudo-chromosomes using Hi-C sequencing data (v2; scaffold N50: 93.2 Mb, total size in chromosomes: 639.6 Mb). Heterozygosity was very low (0.06%), while repeat sequences accounted for 54.87% of the genome, and 23,375 protein-coding genes with an average of 4.79 exons per gene were annotated using de novo, RNA-seq and homology-based approaches. Reconstruction of the historical population size showed a slow continuous contraction, probably related to Cenozoic climate changes. The soursop is the first genome assembled in Annonaceae, supporting further studies of floral evolution in magnoliids, providing an essential resource for delineating relationships of ancient angiosperm lineages. Both genome-assisted improvement and conservation efforts will be strengthened by the availability of the soursop genome. As a community resource, this assembly will further strengthen the role of Annonaceae as model species for research on the ecology, evolution and domestication potential of tropical species in pomology and agroforestry.

Joining forces in Ochnaceae phylogenomics: a tale of two targeted sequencing probe kits

PREMISE

Both universal and family-specific targeted sequencing probe kits are becoming widely used for reconstruction of phylogenetic relationships in angiosperms. Within the pantropical Ochnaceae, we show that with careful data filtering, universal kits are equally as capable in resolving intergeneric relationships as custom probe kits. Furthermore, we show the strength in combining data from both kits to mitigate bias and provide a more robust result to resolve evolutionary relationships.

METHODS

We sampled 23 Ochnaceae genera and used targeted sequencing with two probe kits, the universal Angiosperms353 kit and a family-specific kit. We used maximum likelihood inference with a concatenated matrix of loci and multispecies-coalescence approaches to infer relationships in the family. We explored phylogenetic informativeness and the impact of missing data on resolution and tree support.

RESULTS

For the Angiosperms353 data set, the concatenation approach provided results more congruent with those of the Ochnaceae-specific data set. Filtering missing data was most impactful on the Angiosperms353 data set, with a relaxed threshold being the optimum scenario. The Ochnaceae-specific data set resolved consistent topologies using both inference methods, and no major improvements were obtained after data filtering. Merging of data obtained with the two kits resulted in a well-supported phylogenetic tree.

CONCLUSIONS

The Angiosperms353 data set improved upon data filtering, and missing data played an important role in phylogenetic reconstruction. The Angiosperms353 data set resolved the phylogenetic backbone of Ochnaceae as equally well as the family specific data set. All analyses indicated that both Sauvagesia L. and Campylospermum Tiegh. as currently circumscribed are polyphyletic and require revised delimitation.

Development and phylogenetic utilities of a new set of single-/low-copy nuclear genes in Senecioneae (Asteraceae), with new insights into the tribal position and the relationships within subtribe Tussilagininae

The tribe Senecioneae is one of the largest tribes in Asteraceae, with a nearly cosmopolitan distribution. Despite great efforts devoted to elucidate the evolution of Senecioneae, many questions still remain concerning the systematics of this group, from the tribal circumscription and position to species relationships in many genera. The hybridization-based target enrichment method of next-generation sequencing has been accepted as a promising approach to resolve phylogenetic problems. We herein develop a set of single-/low-copy genes for Senecioneae, and test their phylogenetic utilities. Our results demonstrate that these genes work highly efficiently for Senecioneae, with a high average gene recovery of 98.8% across the tribe and recovering robust phylogenetic hypotheses at different levels. In particular, the delimitation of the Senecioneae has been confirmed to include Abrotanella and exclude Doronicum, with the former sister to core Senecioneae and the latter shown to be more closely related to Calenduleae. Moreover, Doronicum and Calenduleae are inferred to be the closest relatives of Senecioneae, which is a new hypothesis well supported by statistical topology tests, morphological evidence, and the profile of pyrrolizidine alkaloids, a special kind of chemical characters generally used to define Senecioneae. Furthermore, this study suggests a complex reticulation history in the diversification of Senecioneae, accounting for the prevalence of polyploid groups in the tribe. With subtribe Tussilagininae s.str. as a case study showing a more evident pattern of gene duplication, we further explored reconstructing the phylogeny in the groups with high ploidy levels. Our results also demonstrate that tree topologies based on sorted paralogous copies are stable across different methods of phylogenetic inference, and more congruent with the morphological evidence and the results of previous phylogenetic studies.

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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document}) and 2333 genera (17\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document}). 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.]

Phylogenomics of Gesneriaceae using targeted capture of nuclear genes

Gesneriaceae (ca. 3400 species) is a pantropical plant family with a wide range of growth form and floral morphology that are associated with repeated adaptations to different environments and pollinators. Although Gesneriaceae systematics has been largely improved by the use of Sanger sequencing data, our understanding of the evolutionary history of the group is still far from complete due to the limited number of informative characters provided by this type of data. To overcome this limitation, we developed here a Gesneriaceae-specific gene capture kit targeting 830 single-copy loci (776,754 bp in total), including 279 genes from the Universal Angiosperms-353 kit. With an average of 557,600 reads and 87.8% gene recovery, our target capture was successful across the family Gesneriaceae and also in other families of Lamiales. From our bait set, we selected the most informative 418 loci to resolve phylogenetic relationships across the entire Gesneriaceae family using maximum likelihood and coalescent-based methods. Upon testing the phylogenetic performance of our baits on 78 taxa representing 20 out of 24 subtribes within the family, we showed that our data provided high support for the phylogenetic relationships among the major lineages, and were able to provide high resolution within more recent radiations. Overall, the molecular resources we developed here open new perspectives for the study of Gesneriaceae phylogeny at different taxonomical levels and the identification of the factors underlying the diversification of this plant group.

Complete chloroplast sequence of Fenerivia ghesquiereana (Annonaceae): a rare and endemic tree from Madagascar

Fenerivia ghesquiereana (Cavaco & Keraudren) R.M.K. Saunders (Annaonaceae) is a rare and endemic tree restricted to the warm tropical forests of Madagascar, a major global hotspot for biodiversity. Species in the genus are mostly restricted to a thin belt along the eastern edge of the island and remain under intense pressure due to deforestation for logging, mining and slash-and-burn agriculture (‘tavy’), despite national pledges to curb biodiversity losses and increase park protection to combat illegal logging and poaching. Here we report the complete chloroplast genome sequence of this species to support ongoing efforts to complete the (sub-)tribal classification of the family. The chloroplast sequence of F. ghesquiereana was 160,194 bp in length, including two inverted repeat regions of 26,093 bp, a large single-copy region of 89,041 bp and a small single-copy region of 18,967 bp. A total of 160 genes were annotated, of which 115 are coding, 37 are tRNA genes, and eight are rRNA genes. The overall GC content was 39%; this was higher in the IRs (43.4%) when compared to the LSC (30.7%) and the SSC (33.9%) regions. A Maximum Likelihood phylogenetic analysis with a selection of other plastomes in Annonaceae placed F. ghesquiereana as sister to Meiogyne hainanensis (Merr.) in subfamily Malmeoideae.

The Effects of Herbarium Specimen Characteristics on Short-Read NGS Sequencing Success in Nearly 8000 Specimens: Old, Degraded Samples Have Lower DNA Yields but Consistent Sequencing Success

Phylogenetic datasets are now commonly generated using short-read sequencing technologies unhampered by degraded DNA, such as that often extracted from herbarium specimens. The compatibility of these methods with herbarium specimens has precipitated an increase in broad sampling of herbarium specimens for inclusion in phylogenetic studies. Understanding which sample characteristics are predictive of sequencing success can guide researchers in the selection of tissues and specimens most likely to yield good results. Multiple recent studies have considered the relationship between sample characteristics and DNA yield and sequence capture success. Here we report an analysis of the relationship between sample characteristics and sequencing success for nearly 8,000 herbarium specimens. This study, the largest of its kind, is also the first to include a measure of specimen quality ("greenness") as a predictor of DNA sequencing success. We found that taxonomic group and source herbarium are strong predictors of both DNA yield and sequencing success and that the most important specimen characteristics for predicting success differ for DNA yield and sequencing: greenness was the strongest predictor of DNA yield, and age was the strongest predictor of proportion-on-target reads recovered. Surprisingly, the relationship between age and proportion-on-target reads is the inverse of expectations; older specimens performed slightly better in our capture-based protocols. We also found that DNA yield itself is not a strong predictor of sequencing success. Most literature on DNA sequencing from herbarium specimens considers specimen selection for optimal DNA extraction success, which we find to be an inappropriate metric for predicting success using next-generation sequencing technologies.

Resolving Recalcitrant Clades in the Pantropical Ochnaceae: Insights From Comparative Phylogenomics of Plastome and Nuclear Genomic Data Derived From Targeted Sequencing

Plastid DNA sequence data have been traditionally widely used in plant phylogenetics because of the high copy number of plastids, their uniparental inheritance, and the blend of coding and non-coding regions with divergent substitution rates that allow the reconstruction of phylogenetic relationships at different taxonomic ranks. In the present study, we evaluate the utility of the plastome for the reconstruction of phylogenetic relationships in the pantropical plant family Ochnaceae (Malpighiales). We used the off-target sequence read fraction of a targeted sequencing study (targeting nuclear loci only) to recover more than 100 kb of the plastid genome from the majority of the more than 200 species of Ochnaceae and all but two genera using de novo and reference-based assembly strategies. Most of the recalcitrant nodes in the family's backbone were resolved by our plastome-based phylogenetic inference, corroborating the most recent classification system of Ochnaceae and findings from a phylogenomic study based on nuclear loci. Nonetheless, the phylogenetic relationships within the major clades of tribe Ochnineae, which comprise about two thirds of the family's species diversity, received mostly low support. Generally, the phylogenetic resolution was lowest at the infrageneric level. Overall there was little phylogenetic conflict compared to a recent analysis of nuclear loci. Effects of taxon sampling were invoked as the most likely reason for some of the few well-supported discords. Our study demonstrates the utility of the off-target fraction of a target enrichment study for assembling near-complete plastid genomes for a large proportion of samples.

Individualistic evolutionary responses of Central African rain forest plants to Pleistocene climatic fluctuations

Understanding the evolutionary dynamics of genetic diversity is fundamental for species conservation in the face of climate change, particularly in hyper-diverse biomes. Species in a region may respond similarly to climate change, leading to comparable evolutionary dynamics, or individualistically, resulting in dissimilar patterns. The second-largest expanse of continuous tropical rain forest (TRF) in the world is found in Central Africa. Here, present-day patterns of genetic structure are thought to be dictated by repeated expansion and contraction of TRFs into and out of refugia during Pleistocene climatic fluctuations. This refugia model implies a common response to past climate change. However, given the unrivalled diversity of TRFs, species could respond differently because of distinct environmental requirements or ecological characteristics. To test this, we generated genome-wide sequence data for >700 individuals of seven codistributed plants from Lower Guinea in Central Africa. We inferred species' evolutionary and demographic histories within a comparative phylogeographic framework. Levels of genetic structure varied among species and emerged primarily during the Pleistocene, but divergence events were rarely concordant. Demographic trends ranged from repeated contraction and expansion to continuous growth. Furthermore, patterns in genetic variation were linked to disparate environmental factors, including climate, soil, and habitat stability. Using a strict refugia model to explain past TRF dynamics is too simplistic. Instead, individualistic evolutionary responses to Pleistocene climatic fluctuations have shaped patterns in genetic diversity. Predicting the future dynamics of TRFs under climate change will be challenging, and more emphasis is needed on species ecology to better conserve TRFs worldwide.

Hybrid capture of 964 nuclear genes resolves evolutionary relationships in the mimosoid legumes and reveals the polytomous origins of a large pantropical radiation

PREMISE

Targeted enrichment methods facilitate sequencing of hundreds of nuclear loci to enhance phylogenetic resolution and elucidate why some parts of the "tree of life" are difficult (if not impossible) to resolve. The mimosoid legumes are a prominent pantropical clade of ~3300 species of woody angiosperms for which previous phylogenies have shown extensive lack of resolution, especially among the species-rich and taxonomically challenging ingoids.

METHODS

We generated transcriptomes to select low-copy nuclear genes, enrich these via hybrid capture for representative species of most mimosoid genera, and analyze the resulting data using de novo assembly and various phylogenomic tools for species tree inference. We also evaluate gene tree support and conflict for key internodes and use phylogenetic network analysis to investigate phylogenetic signal across the ingoids.

RESULTS

Our selection of 964 nuclear genes greatly improves phylogenetic resolution across the mimosoid phylogeny and shows that the ingoid clade can be resolved into several well-supported clades. However, nearly all loci show lack of phylogenetic signal for some of the deeper internodes within the ingoids.

CONCLUSIONS

Lack of resolution in the ingoid clade is most likely the result of hyperfast diversification, potentially causing a hard polytomy of six or seven lineages. The gene set for targeted sequencing presented here offers great potential to further enhance the phylogeny of mimosoids and the wider Caesalpinioideae with denser taxon sampling, to provide a framework for taxonomic reclassification, and to study the ingoid radiation.

Enabling evolutionary studies at multiple scales in Apocynaceae through Hyb-Seq

PREMISE

Apocynaceae is the 10th largest flowering plant family and a focus for study of plant-insect interactions, especially as mediated by secondary metabolites. However, it has few genomic resources relative to its size. Target capture sequencing is a powerful approach for genome reduction that facilitates studies requiring data from the nuclear genome in non-model taxa, such as Apocynaceae.

METHODS

Transcriptomes were used to design probes for targeted sequencing of putatively single-copy nuclear genes across Apocynaceae. The sequences obtained were used to assess the success of the probe design, the intrageneric and intraspecific variation in the targeted genes, and the utility of the genes for inferring phylogeny.

RESULTS

From 853 candidate nuclear genes, 835 were consistently recovered in single copy and were variable enough for phylogenomics. The inferred gene trees were useful for coalescent-based species tree analysis, which showed all subfamilies of Apocynaceae as monophyletic, while also resolving relationships among species within the genus Apocynum. Intraspecific comparison of Elytropus chilensis individuals revealed numerous single-nucleotide polymorphisms with potential for use in population-level studies.

DISCUSSION

Community use of this Hyb-Seq probe set will facilitate and promote progress in the study of Apocynaceae across scales from population genomics to phylogenomics.

Larvicidal effect from different Annonaceae species on Culex quinquefasciatus

The recent outbreaks of mosquito-borne diseases highlighted the pivotal importance of mosquito vector control in tropical areas worldwide. Several strategies have been developed to control vector populations and disease transmission in endemic areas. The steps to obtain natural active compounds involve the pre-selection in a biological model and subsequently evaluation on specific models. The present study reports the evaluation of 35 extracts, fractions, and essential oils obtained from five species from the Annonaceae family on Artemia salina and Culex quinquefasciatus. The A. salina results were used as a pre-screening for larvicidal test about mosquitoes. A correlation of biological activity in both bioassays was observed for the hydroethanolic extracts and their respective hexane and chloroform fractions of the leaves of Annona species, except A. nutans. The same correlation was also observed for all tested essential oils and petroleum ether extracts from Duguetia species. It was possible to limit an interval of lethality about A. salina, which has a corresponding range to the larvicidal test against the mosquito. The main components present in D. lanceolata essential oil or enriched fraction were α-selinene, aristolochene, (E)-caryophyllene, and (E)-calamenene. For D. furfuracea, the main components present of the underground parts were (E)-asarone, 2,4,5-trimethoxystyrene, spathulenol, and bicyclogermacrene for aerial parts. The A. salina test could be used as a model for the pre-screening of larvicidal activity.

Phylogenomic relationships and character evolution of the grape family (Vitaceae)

The grape family consists of 16 genera and ca. 950 species. It is best known for the economically important fruit crop - the grape Vitis vinifera. The deep phylogenetic relationships and character evolution of the grape family have attracted the attention of researchers in recent years. We herein reconstruct the phylogenomic relationships within Vitaceae using nuclear and plastid genes based on the Hyb-Seq approach and test the newly proposed classification system of the family. The five tribes of the grape family, including Ampelopsideae, Cayratieae, Cisseae, Parthenocisseae, and Viteae, are each robustly supported by both nuclear and chloroplast genomic data and the backbone relationships are congruent with previous reports. The cupular floral disc (raised above and free from ovary at the upper part) is an ancestral state of Vitaceae, with the inconspicuous floral disc as derived in the tribe Parthenocisseae, and the state of adnate to the ovary as derived in the tribe Viteae. The 5-merous floral pattern was inferred to be the ancestral in Vitaceae, with the 4-merous flowers evolved at least two times in the family. The compound dichasial cyme (cymose with two secondary axes) is ancestral in Vitaceae and the thyrse inflorescence (a combination of racemose and cymose branching) in tribe Viteae is derived. The ribbon-like trichome only evolved once in Vitaceae, as a synapomorphy for the tribe Viteae.

Phylogenomic approaches reveal how climate shapes patterns of genetic diversity in an African rain forest tree species

The world's second largest expanse of tropical rainforest is in Central Africa, and it harbours enormous species diversity. Population genetic studies have consistently revealed significant structure across Central African rainforest plants. In particular, previous studies have repeatedly demonstrated a north-south genetic discontinuity around the equatorial line, in a continuous expanse of rainforest where a climatic inversion is documented. Here, we took a phylogeographic approach by sequencing 351 nuclear markers in 112 individuals across the distribution of the African rainforest tree species Annickia affinis (Annonaceae). We showed for the first time that the north-south divide is the result of a single, major colonization event across the climatic inversion from an ancestral population located in Gabon. We suggested that differences in ecological niche of populations located on either side of this inversion may have contributed to this phylogenetic discontinuity. We found evidence for inland dispersal, predominantly in northern areas, and variable demographic histories among genetic clusters, indicating that populations responded differently to past climate change. We show how newly developed genomic tools can provide invaluable insights into our understanding of tropical rainforest evolutionary dynamics.

Museomics for reconstructing historical floristic exchanges: Divergence of stone oaks across Wallacea

Natural history collections and tropical tree diversity are both treasure troves of biological and evolutionary information, but their accessibility for scientific study is impeded by a number of properties. DNA in historical specimens is generally highly fragmented, complicating the recovery of high-grade genetic material. Furthermore, our understanding of hyperdiverse, wide-spread tree assemblages is obstructed by extensive species ranges, fragmented knowledge of tropical tree diversity and phenology, and a widespread lack of species-level diagnostic characters, prohibiting the collecting of readily identifiable specimens which can be used to build, revise or strengthen taxonomic frameworks. This, in turn, delays the application of downstream conservation action. A sizable component of botanical collections are sterile-thus eluding identification and are slowing down progress in systematic treatments of tropical biodiversity. With rapid advances in genomics and bioinformatic approaches to biodiversity research, museomics is emerging as a new field breathing life into natural collections that have been built up over centuries. Using MIGseq (multiplexed ISSR genotyping by sequencing), we generated 10,000s of short loci, for both freshly collected materials and museum specimens (aged >100 years) of Lithocarpus-a widespread tropical tree genus endemic to the Asian tropics. Loci recovery from historical and recently collected samples was not affected by sample age and preservation history of the study material, underscoring the reliability and flexibility of the MIGseq approach. Phylogenomic inference and biogeographic reconstruction across insular Asia, highlights repeated migration and diversification patterns between continental regions and islands. Results indicate that co-occurring insular species at the extremity of the distribution range are not monophyletic, raising the possibility of multiple independent dispersals along the outer edge of Wallacea. This suggests that dispersal of large seeded tree genera throughout Malesia and across Wallacea may have been less affected by large geographic distances and the presence of marine barriers than generally assumed. We demonstrate the utility of MIGseq in museomic studies using non-model taxa, presenting the first range-wide genomic assessment of Lithocarpus and tropical Fagaceae as a proof-of-concept. Our study shows the potential for developing innovative genomic approaches to improve the capture of novel evolutionary signals using valuable natural history collections of hyperdiverse taxa.

Nuclear phylogenomic analyses of asterids conflict with plastome trees and support novel relationships among major lineages

PREMISE

Discordance between nuclear and organellar phylogenies (cytonuclear discordance) is a well-documented phenomenon at shallow evolutionary levels but has been poorly investigated at deep levels of plant phylogeny. Determining the extent of cytonuclear discordance across major plant lineages is essential not only for elucidating evolutionary processes, but also for evaluating the currently used framework of plant phylogeny, which is largely based on the plastid genome.

METHODS

We present a phylogenomic examination of a major angiosperm clade (Asteridae) based on sequence data from the nuclear, plastid, and mitochondrial genomes as a means of evaluating currently accepted relationships inferred from the plastome and exploring potential sources of genomic conflict in this group.

RESULTS

We recovered at least five instances of well-supported cytonuclear discordance concerning the placements of major asterid lineages (i.e., Ericales, Oncothecaceae, Aquifoliales, Cassinopsis, and Icacinaceae). We attribute this conflict to a combination of incomplete lineage sorting and hybridization, the latter supported in part by previously inferred whole-genome duplications.

CONCLUSIONS

Our results challenge several long-standing hypotheses of asterid relationships and have implications for morphological character evolution and for the importance of ancient whole-genome duplications in early asterid evolution. These findings also highlight the value of reevaluating broad-scale angiosperm and green-plant phylogeny with nuclear genomic data.

Unraveling the Phylogenomic Relationships of the Most Diverse African Palm Genus Raphia (Calamoideae, Arecaceae)

Palms are conspicuous floristic elements across the tropics. In continental Africa, even though there are less than 70 documented species, they are omnipresent across the tropical landscape. The genus Raphia has 20 accepted species in Africa and one species endemic to the Neotropics. It is the most economically important genus of African palms with most of its species producing food and construction material. Raphia is divided into five sections based on inflorescence morphology. Nevertheless, the taxonomy of Raphia is problematic with no intra-generic phylogenetic study available. We present a phylogenetic study of the genus using a targeted exon capture approach sequencing of 56 individuals representing 18 out of the 21 species. Our results recovered five well supported clades within the genus. Three sections correspond to those based on inflorescence morphology. R. regalis is strongly supported as sister to all other Raphia species and is placed into a newly described section: Erectae. Overall, morphological based identifications agreed well with our phylogenetic analyses, with 12 species recovered as monophyletic based on our sampling. Species delimitation analyses recovered 17 or 23 species depending on the confidence level used. Species delimitation is especially problematic in the Raphiate and Temulentae sections. In addition, our clustering analysis using SNP data suggested that individual clusters matched geographic distribution. The Neotropical species R. taedigera is supported as a distinct species, rejecting the hypothesis of a recent introduction into South America. Our analyses support the hypothesis that the Raphia individuals from Madagascar are potentially a distinct species different from the widely distributed R. farinifera. In conclusion, our results support the infra generic classification of Raphia based on inflorescence morphology, which is shown to be phylogenetically useful. Classification and species delimitation within sections remains problematic even with our phylogenomic approach. Certain widely distributed species could potentially contain cryptic species. More in-depth studies should be undertaken using morphometrics, increased sampling, and more variable markers. Our study provides a robust phylogenomic framework that enables further investigation on the biogeographic history, morphological evolution, and other eco-evolutionary aspects of this charismatic, socially, and economically important palm genus.

A Guide to Carrying Out a Phylogenomic Target Sequence Capture Project

High-throughput DNA sequencing techniques enable time- and cost-effective sequencing of large portions of the genome. Instead of sequencing and annotating whole genomes, many phylogenetic studies focus sequencing effort on large sets of pre-selected loci, which further reduces costs and bioinformatic challenges while increasing coverage. One common approach that enriches loci before sequencing is often referred to as target sequence capture. This technique has been shown to be applicable to phylogenetic studies of greatly varying evolutionary depth. Moreover, it has proven to produce powerful, large multi-locus DNA sequence datasets suitable for phylogenetic analyses. However, target capture requires careful considerations, which may greatly affect the success of experiments. Here we provide a simple flowchart for designing phylogenomic target capture experiments. We discuss necessary decisions from the identification of target loci to the final bioinformatic processing of sequence data. We outline challenges and solutions related to the taxonomic scope, sample quality, and available genomic resources of target capture projects. We hope this review will serve as a useful roadmap for designing and carrying out successful phylogenetic target capture studies.

An empirical assessment of a single family-wide hybrid capture locus set at multiple evolutionary timescales in Asteraceae

PREMISE

Hybrid capture with high-throughput sequencing (Hyb-Seq) is a powerful tool for evolutionary studies. The applicability of an Asteraceae family-specific Hyb-Seq probe set and the outcomes of different phylogenetic analyses are investigated here.

METHODS

Hyb-Seq data from 112 Asteraceae samples were organized into groups at different taxonomic levels (tribe, genus, and species). For each group, data sets of non-paralogous loci were built and proportions of parsimony informative characters estimated. The impacts of analyzing alternative data sets, removing long branches, and type of analysis on tree resolution and inferred topologies were investigated in tribe Cichorieae.

RESULTS

Alignments of the Asteraceae family-wide Hyb-Seq locus set were parsimony informative at all taxonomic levels. Levels of resolution and topologies inferred at shallower nodes differed depending on the locus data set and the type of analysis, and were affected by the presence of long branches.

DISCUSSION

The approach used to build a Hyb-Seq locus data set influenced resolution and topologies inferred in phylogenetic analyses. Removal of long branches improved the reliability of topological inferences in maximum likelihood analyses. The Astereaceae Hyb-Seq probe set is applicable at multiple taxonomic depths, which demonstrates that probe sets do not necessarily need to be lineage-specific.

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.

Tackling Rapid Radiations With Targeted Sequencing

In phylogenetic studies across angiosperms, at various taxonomic levels, polytomies have persisted despite efforts to resolve them by increasing sampling of taxa and loci. The large amount of genomic data now available and statistical tools to analyze them provide unprecedented power for phylogenetic inference. Targeted sequencing has emerged as a strong tool for estimating species trees in the face of rapid radiations, lineage sorting, and introgression. Evolutionary relationships in Cyperaceae have been studied mostly using Sanger sequencing until recently. Despite ample taxon sampling, relationships in many genera remain poorly understood, hampered by diversification rates that outpace mutation rates in the loci used. The C4 Cyperus clade of the genus Cyperus has been particularly difficult to resolve. Previous studies based on a limited set of markers resolved relationships among Cyperus species using the C3 photosynthetic pathway, but not among C4 Cyperus clade taxa. We test the ability of two targeted sequencing kits to resolve relationships in the C4 Cyperus clade, the universal Angiosperms-353 kit and a Cyperaceae-specific kit. Sequences of the targeted loci were recovered from data generated with both kits and used to investigate overlap in data between kits and relative efficiency of the general and custom approaches. The power to resolve shallow-level relationships was tested using a summary species tree method and a concatenated maximum likelihood approach. High resolution and support are obtained using both approaches, but high levels of missing data disproportionately impact the latter. Targeted sequencing provides new insights into the evolution of morphology in the C4 Cyperus clade, demonstrating for example that the former segregate genus Alinula is polyphyletic despite its seeming morphological integrity. An unexpected result is that the Cyperus margaritaceus-Cyperus niveus complex comprises a clade separate from and sister to the core C4 Cyperus clade. Our results demonstrate that data generated with a family-specific kit do not necessarily have more power than those obtained with a universal kit, but that data generated with different targeted sequencing kits can often be merged for downstream analyses. Moreover, our study contributes to the growing consensus that targeted sequencing data are a powerful tool in resolving rapid radiations.