A target enrichment method for gathering phylogenetic information from hundreds of loci: An example from the Compositae

A Review of Probe-Based Enrichment Methods to Inform Plant Virus Diagnostics

Modern diagnostic techniques based on DNA sequence similarity are currently the gold standard for the detection of existing and emerging pathogens. Whilst individual assays are inexpensive to use, assay development is costly and carries risks of not being sensitive or specific enough to capture an increasingly diverse range of targets. Sequencing can provide the entire nucleic acid content of a sample and may be used to identify all pathogens present in the sample when the depth of coverage is sufficient. Targeted enrichment techniques have been used to increase sequence coverage and improve the sensitivity of detection within virus samples, specifically, to capture sequences for a range of different viruses or increase the number of reads from low-titre virus infections. Vertebrate viruses have been well characterised using in-solution hybridisation capture to target diverse virus families. The use of probes for genotyping and strain identification has been limited in plants, and uncertainty around sensitivity is an impediment to the development of a large-scale virus panel to use within regulatory settings and diagnostic pipelines. This review aims to compare significant studies that have used targeted enrichment of viruses to identify approaches to probe design and potential for use in plant virus detection and characterisation.

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

Development and evolution of the Asteraceae capitulum

Asteraceae represent one of the largest and most diverse families of plants. The evolutionary success of this family has largely been contributed to their unique inflorescences, capitula that mimic solitary flowers but are typically aggregates of multiple florets. Here, we summarize the recent molecular and genetic level studies that have promoted our understanding of the development and evolution of capitula. We focus on new results on patterning of the enlarged meristem resulting in the iconic phyllotactic arrangement of florets in Fibonacci numbers of spirals. We also summarize the current understanding of the genetic networks regulating the characteristic reproductive traits in the family such as floral dimorphism and differentiation of highly specialized floral organs. So far, developmental studies in Asteraceae are still limited to a very narrow selection of model species. Along with the recent advancements in genomics and phylogenomics, Asteraceae and its relatives provide an outstanding model clade for extended evo-devo studies to exploit the morphological diversity and the underlying molecular networks and to translate this knowledge to the breeding of the key crops in the family.

Phylogenomics and divergence pattern of Polygonatum (Asparagaceae: Polygonateae) in the north temperate region

Polygonatum is the largest genus of tribe Polygonateae (Asparagaceae) and is widely distributed in the temperate Northern Hemisphere, especially well diversified in southwestern China to northeastern Asia. Phylogenetic relationships of many species are still controversial. Hence it is necessary to clarify their phylogenetic relationships and infer possible reticulate relationships for the genus. In this study, genome-wide data of 43 species from Polygonatum and its closely related taxa were obtained by Hyb-Seq sequencing. The phylogenetic trees constructed from genome-wide nuclear and chloroplast sequences strongly supported the monophyly of Polygonatum with division into three major clades. A high level of incongruence was detected between nuclear and chloroplast trees as well as among gene trees within the genus, but all occurred within each major clade. However, introgression tests and reticulate evolution analyses revealed low level of gene flow and weak introgression events in the genus, suggesting hybridization and introgression were not dominant during the evolutionary diversification of Polygonatum in the Northern Hemisphere. This study provides important insights into reconstructing evolutionary relationships and speciation pattern of taxa from the north temperate flora.

Compositae-ParaLoss-1272: A complementary sunflower-specific probe set reduces paralogs in phylogenomic analyses of complex systems

Premise

A family-specific probe set for sunflowers, Compositae-1061, enables family-wide phylogenomic studies and investigations at lower taxonomic levels, but may lack resolution at genus to species levels, especially in groups complicated by polyploidy and hybridization.

Methods

We developed a Hyb-Seq probe set, Compositae-ParaLoss-1272, that targets orthologous loci in Asteraceae. We tested its efficiency across the family by simulating target enrichment sequencing in silico. Additionally, we tested its effectiveness at lower taxonomic levels in the historically complex genus Packera. We performed Hyb-Seq with Compositae-ParaLoss-1272 for 19 Packera taxa that were previously studied using Compositae-1061. The resulting sequences from each probe set, plus a combination of both, were used to generate phylogenies, compare topologies, and assess node support.

Results

We report that Compositae-ParaLoss-1272 captured loci across all tested Asteraceae members, had less gene tree discordance, and retained longer loci than Compositae-1061. Most notably, Compositae-ParaLoss-1272 recovered substantially fewer paralogous sequences than Compositae-1061, with only ~5% of the recovered loci reporting as paralogous, compared to ~59% with Compositae-1061.

Discussion

Given the complexity of plant evolutionary histories, assigning orthology for phylogenomic analyses will continue to be challenging. However, we anticipate Compositae-ParaLoss-1272 will provide improved resolution and utility for studies of complex groups and lower taxonomic levels in the sunflower family.

Impact of the climatic changes in the Pliocene-Pleistocene transition on Irano-Turanian species. The radiation of genus Jurinea (Compositae)

The Irano-Turanian region is one of the world's richest floristic regions and the centre of diversity for numerous xerophytic plant lineages. However, we still have limited knowledge on the timing of evolution and biogeographic history of its flora, and potential drivers of diversification remain underexplored. To fill this knowledge gap, we focus on the Eurasian genus Jurinea (ca. 200 species), one of the largest plant radiations that diversified in the region. We applied a macroevolutionary integrative approach to explicitly test diversification hypotheses and investigate the relative roles of geography vs. ecology and niche conservatism vs. niche lability in speciation processes. To do so, we gathered a sample comprising 77% of total genus richness and obtained data about (1) its phylogenetic history, recovering 502 nuclear loci sequences; (2) growth forms; (3) ecological niche, compiling data of 21 variables for more than 2500 occurrences; and (4) paleoclimatic conditions, to estimate climatic stability. Our results revealed that climate was a key factor in the evolutionary dynamics of Jurinea. The main diversification and biogeographic events that occurred during past climate changes, which led to colder and drier conditions, are the following: (1) the origin of the genus (10.7 Ma); (2) long-distance dispersals from the Iranian Plateau to adjacent regions (∼7-4 Ma); and (3) the diversification shift during Pliocene-Pleistocene Transition (ca. 3 Ma), when net diversification rate almost doubled. Our results supported the pre-adaptation hypothesis, i.e., the evolutionary success of Jurinea was linked to the retention of the ancestral niche adapted to aridity. Interestingly, the paleoclimatic analyses revealed that in the Iranian Plateau long-term climatic stability favoured old-lineage persistence, resulting in current high species richness of semi-arid and cold adapted clades; whereas moderate climate oscillations stimulated allopatric diversification in the lineages distributed in the Circumboreal region. In contrast, growth form lability and high niche disparity among closely related species in the Central Asian clade suggest adaptive radiation to mountain habitats. In sum, the radiation of Jurinea is the result of both adaptive and non-adaptive processes influenced by climatic, orogenic and ecological factors.

LoCoLotive: In silico mining for low-copy nuclear loci based on target capture probe sets and arbitrary reference genomes

Premise

Universal target enrichment probe kits are used to circumvent the individual identification of loci suitable for phylogenetic studies in a given taxon. Under certain circumstances, however, target capture can be inefficient and costly, and lower numbers of marker loci may be sufficient. We therefore propose a computational pipeline that enables the easy identification of a subset of promising candidate loci for a taxon of interest.

Methods and Results

Target sequences used for probe design are filtered based on an assembled reference genome, resulting in presumably intron-containing single-copy loci as present in the reference taxon. The applicability of the proposed approach is demonstrated based on two probe kits (universal and family-specific) in combination with several publicly available reference genomes.

Conclusions

Guided by commercial probe kits, LoCoLotive enables fast and cost-efficient marker mining. Its accuracy mainly depends on the quality of the reference genome and its relatedness to the taxa under study.

A target enrichment probe set for resolving phylogenetic relationships in the coffee family, Rubiaceae

Premise

Rubiaceae is among the most species-rich plant families, as well as one of the most morphologically and geographically diverse. Currently available phylogenies have mostly relied on few genomic and plastid loci, as opposed to large-scale genomic data. Target enrichment provides the ability to generate sequence data for hundreds to thousands of phylogenetically informative, single-copy loci, which often leads to improved phylogenetic resolution at both shallow and deep taxonomic scales; however, a publicly accessible Rubiaceae-specific probe set that allows for comparable phylogenetic inference across clades is lacking.

Methods

Here, we use publicly accessible genomic resources to identify putatively single-copy nuclear loci for target enrichment in two Rubiaceae groups: tribe Hillieae (Cinchonoideae) and tribal complex Palicoureeae+Psychotrieae (Rubioideae). We sequenced 2270 exonic regions corresponding to 1059 loci in our target clades and generated in silico target enrichment sequences for other Rubiaceae taxa using our designed probe set. To test the utility of our probe set for phylogenetic inference across Rubiaceae, we performed a coalescent-aware phylogenetic analysis using a subset of 27 Rubiaceae taxa from 10 different tribes and three subfamilies, and one outgroup in Apocynaceae.

Results

We recovered an average of 75% and 84% of targeted exons and loci, respectively, per Rubiaceae sample. Probes designed using genomic resources from a particular subfamily were most efficient at targeting sequences from taxa in that subfamily. The number of paralogs recovered during assembly varied for each clade. Phylogenetic inference of Rubiaceae with our target regions resolves relationships at various scales. Relationships are largely consistent with previous studies of relationships in the family with high support (≥0.98 local posterior probability) at nearly all nodes and evidence of gene tree discordance.

Discussion

Our probe set, which we call Rubiaceae2270x, was effective for targeting loci in species across and even outside of Rubiaceae. This probe set will facilitate phylogenomic studies in Rubiaceae and advance systematics and macroevolutionary studies in the family.

Next-generation fungal identification using target enrichment and Nanopore sequencing

BACKGROUND

Rapid and accurate pathogen identification is required for disease management. Compared to sequencing entire genomes, targeted sequencing may be used to direct sequencing resources to genes of interest for microbe identification and mitigate the low resolution that single-locus molecular identification provides. This work describes a broad-spectrum fungal identification tool developed to focus high-throughput Nanopore sequencing on genes commonly employed for disease diagnostics and phylogenetic inference.

RESULTS

Orthologs of targeted genes were extracted from 386 reference genomes of fungal species spanning six phyla to identify homologous regions that were used to design the baits used for enrichment. To reduce the cost of producing probes without diminishing the phylogenetic power, DNA sequences were first clustered, and then consensus sequences within each cluster were identified to produce 26,000 probes that targeted 114 genes. To test the efficacy of our probes, we applied the technique to three species representing Ascomycota and Basidiomycota fungi. The efficiency of enrichment, quantified as mean target coverage over the mean genome-wide coverage, ranged from 200 to 300. Furthermore, enrichment of long reads increased the depth of coverage across the targeted genes and into non-coding flanking sequence. The assemblies generated from enriched samples provided well-resolved phylogenetic trees for taxonomic assignment and molecular identification.

CONCLUSIONS

Our work provides data to support the utility of targeted Nanopore sequencing for fungal identification and provides a platform that may be extended for use with other phytopathogens.

hybpiper-nf and paragone-nf: Containerization and additional options for target capture assembly and paralog resolution

Premise

The HybPiper pipeline has become one of the most widely used tools for the assembly of target capture data for phylogenomic analysis. After the production of locus sequences and before phylogenetic analysis, the identification of paralogs is a critical step for ensuring the accurate inference of evolutionary relationships. Algorithmic approaches using gene tree topologies for the inference of ortholog groups are computationally efficient and broadly applicable to non-model organisms, especially in the absence of a known species tree.

Methods and Results

We containerized and expanded the functionality of both HybPiper and a pipeline for the inference of ortholog groups, providing novel options for the treatment of target capture sequence data, and allowing seamless use of the outputs of the former as inputs for the latter. The Singularity container presented here includes all dependencies, and the corresponding pipelines (hybpiper-nf and paragone-nf, respectively) are implemented via two Nextflow scripts for easier deployment and to vastly reduce the number of commands required for their use.

Conclusions

The hybpiper-nf and paragone-nf pipelines are easily installed and provide a user-friendly experience and robust results to the phylogenetic community. They are used by the Australian Angiosperm Tree of Life project. The pipelines are available at https://github.com/chrisjackson-pellicle/hybpiper-nf and https://github.com/chrisjackson-pellicle/paragone-nf.

Nano-Strainer: A workflow for the identification of single-copy nuclear loci for plant systematic studies, using target capture kits and Oxford Nanopore long reads

In modern plant systematics, target enrichment enables simultaneous analysis of hundreds of genes. However, when dealing with reticulate or polyploidization histories, few markers may suffice, but often are required to be single-copy, a condition that is not necessarily met with commercial capture kits. Also, large genome sizes can render target capture ineffective, so that amplicon sequencing would be preferable; however, knowledge about suitable loci is often missing. Here, we present a comprehensive workflow for the identification of putative single-copy nuclear markers in a genus of interest, by mining a small dataset from target capture using a few representative taxa. The proposed pipeline assesses sequence variability contained in the data from targeted loci and assigns reads to their respective genes, via a combined BLAST/clustering procedure. Cluster consensus sequences are then examined based on four pre-defined criteria presumably indicative for absence of paralogy. This is done by calculating four specialized indices; loci are ranked according to their performance in these indices, and top-scoring loci are considered putatively single- or low copy. The approach can be applied to any probe set. As it relies on long reads, the present contribution also provides template workflows for processing Nanopore-based target capture data. Obtained markers are further tested and then entered into amplicon sequencing. For the detection of possibly remaining paralogy in these data, which might occur in groups with rampant paralogy, we also employ the long-read assembly tool canu. In diploid representatives of the young Compositae genus Leucanthemum, characterized by high levels of polyploidy, our approach resulted in successful amplification of 13 loci. Modifications to remove traces of paralogy were made in seven of these. A species tree from the markers correctly reproduced main relationships in the genus, however, at low resolution. The presented workflow has the potential to valuably support phylogenetic research, for example in polyploid plant groups.

Repeatedly Northwards and Upwards: Southern African Grasslands Fuel the Colonization of the African Sky Islands in Helichrysum (Compositae)

The Afromontane and Afroalpine areas constitute some of the main biodiversity hotspots of Africa. They are particularly rich in plant endemics, but the biogeographic origins and evolutionary processes leading to this outstanding diversity are poorly understood. We performed phylogenomic and biogeographic analyses of one of the most species-rich plant genera in these mountains, Helichrysum (Compositae-Gnaphalieae). Most previous studies have focused on Afroalpine elements of Eurasian origin, and the southern African origin of Helichrysum provides an interesting counterexample. We obtained a comprehensive nuclear dataset from 304 species (≈50% of the genus) using target-enrichment with the Compositae1061 probe set. Summary-coalescent and concatenation approaches combined with paralog recovery yielded congruent, well-resolved phylogenies. Ancestral range estimations revealed that Helichrysum originated in arid southern Africa, whereas the southern African grasslands were the source of most lineages that dispersed within and outside Africa. Colonization of the tropical Afromontane and Afroalpine areas occurred repeatedly throughout the Miocene-Pliocene. This timing coincides with mountain uplift and the onset of glacial cycles, which together may have facilitated both speciation and intermountain gene flow, contributing to the evolution of the Afroalpine flora.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

A target capture approach for phylogenomic analyses at multiple evolutionary timescales in rosewoods (Dalbergia spp.) and the legume family (Fabaceae)

Understanding the genetic changes associated with the evolution of biological diversity is of fundamental interest to molecular ecologists. The assessment of genetic variation at hundreds or thousands of unlinked genetic loci forms a sound basis to address questions ranging from micro- to macroevolutionary timescales, and is now possible thanks to advances in sequencing technology. Major difficulties are associated with (i) the lack of genomic resources for many taxa, especially from tropical biodiversity hotspots; (ii) scaling the numbers of individuals analysed and loci sequenced; and (iii) building tools for reproducible bioinformatic analyses of such data sets. To address these challenges, we developed target capture probes for genomic studies of the highly diverse, pantropically distributed and economically significant rosewoods (Dalbergia spp.), explored the performance of an overlapping probe set for target capture across the legume family (Fabaceae), and built the general purpose bioinformatic pipeline CaptureAl. Phylogenomic analyses of Malagasy Dalbergia species yielded highly resolved and well supported hypotheses of evolutionary relationships. Population genomic analyses identified differences between closely related species and revealed the existence of a potentially new species, suggesting that the diversity of Malagasy Dalbergia species has been underestimated. Analyses at the family level corroborated previous findings by the recovery of monophyletic subfamilies and many well-known clades, as well as high levels of gene tree discordance, especially near the root of the family. The new genomic and bioinformatic resources, including the Fabaceae1005 and Dalbergia2396 probe sets, will hopefully advance systematics and ecological genetics research in legumes, and promote conservation of the highly diverse and endangered Dalbergia rosewoods.

The relationship between transposable elements and ecological niches in the Greater Cape Floristic Region: A study on the genus Pteronia (Asteraceae)

Non-coding repetitive DNA (repeatome) is an active part of the nuclear genome, involved in its structure, evolution and function. It is dominated by transposable elements (TEs) and satellite DNA and is prone to the most rapid changes over time. The TEs activity presumably causes the global genome reorganization and may play an adaptive or regulatory role in response to environmental challenges. This assumption is applied here for the first time to plants from the Cape Floristic hotspot to determine whether changes in repetitive DNA are related to responses to a harsh, but extremely species-rich environment. The genus Pteronia (Asteraceae) serves as a suitable model group because it shows considerable variation in genome size at the diploid level and has high and nearly equal levels of endemism in the two main Cape biomes, Fynbos and Succulent Karoo. First, we constructed a phylogeny based on multiple low-copy genes that served as a phylogenetic framework for detecting quantitative and qualitative changes in the repeatome. Second, we performed a comparative analysis of the environments of two groups of Pteronia differing in their TEs bursts. Our results suggest that the environmental transition from the Succulent Karoo to the Fynbos is accompanied by TEs burst, which is likely also driving phylogenetic divergence. We thus hypothesize that analysis of rapidly evolving repeatome could serve as an important proxy for determining the molecular basis of lineage divergence in rapidly radiating groups.

Utilizing paralogs for phylogenetic reconstruction has the potential to increase species tree support and reduce gene tree discordance in target enrichment data

The analysis of target enrichment data in phylogenetics lacks optimization toward using paralogs for phylogenetic reconstruction. We developed a novel approach of detecting paralogs and utilizing them for phylogenetic tree inference, by retrieving both ortho- and paralogous copies and creating orthologous alignments, from which the gene trees are built. We implemented this approach in ParalogWizard and demonstrate its performance in plant groups that underwent a whole genome duplication relatively recently: the subtribe Malinae (family Rosaceae), using Angiosperms353 as well as Malinae481 probes, the genus Oritrophium (family Asteraceae), using Compositae1061 probes, and the genus Amomum (family Zingiberaceae), using Zingiberaceae1180 probes. Discriminating between orthologs and paralogs reduced gene tree discordance and increased the species tree support in the case of the Malinae, but not for Oritrophium and Amomum. This may relate to the difference in the proportion of paralogous loci between the datasets, which was highest for the Malinae. Overall, retrieving paralogs for phylogenetic reconstruction following ParalogWizard has the potential to increase the species tree support and reduce gene tree discordance in target enrichment data, particularly if the proportion of paralogous loci is high.

New Insights Into the Relationships Within Subtribe Scorzonerinae (Cichorieae, Asteraceae) Using Hybrid Capture Phylogenomics (Hyb-Seq)

Subtribe Scorzonerinae (Cichorieae, Asteraceae) contains 12 main lineages and approximately 300 species. Relationships within the subtribe, either at inter- or intrageneric levels, were largely unresolved in phylogenetic studies to date, due to the lack of phylogenetic signal provided by traditional Sanger sequencing markers. In this study, we employed a phylogenomics approach (Hyb-Seq) that targets 1,061 nuclear-conserved ortholog loci designed for Asteraceae and obtained chloroplast coding regions as a by-product of off-target reads. Our objectives were to evaluate the potential of the Hyb-Seq approach in resolving the phylogenetic relationships across the subtribe at deep and shallow nodes, investigate the relationships of major lineages at inter- and intrageneric levels, and examine the impact of the different datasets and approaches on the robustness of phylogenetic inferences. We analyzed three nuclear datasets: exon only, excluding all potentially paralogous loci; exon only, including loci that were only potentially paralogous in 1-3 samples; exon plus intron regions (supercontigs); and the plastome CDS region. Phylogenetic relationships were reconstructed using both multispecies coalescent and concatenation (Maximum Likelihood and Bayesian analyses) approaches. Overall, our phylogenetic reconstructions recovered the same monophyletic major lineages found in previous studies and were successful in fully resolving the backbone phylogeny of the subtribe, while the internal resolution of the lineages was comparatively poor. The backbone topologies were largely congruent among all inferences, but some incongruent relationships were recovered between nuclear and plastome datasets, which are discussed and assumed to represent cases of cytonuclear discordance. Considering the newly resolved phylogenies, a new infrageneric classification of Scorzonera in its revised circumscription is proposed.

The genomic basis of the plant island syndrome in Darwin's giant daisies

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the ‘plant island syndrome’, include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin’s giant daisies.

Many island plant species share a syndrome of characteristic phenotype and life history. Cerca et al. find the genomic basis of the plant island syndrome in one of Darwin’s giant daisies, while separating ancestral genomes in a chromosome-resolved polyploid assembly.

Phylogenomic loci define the generic boundaries of Gochnatieae and improve resolution at the species level in Moquiniastrum (Compositae)

Understanding the evolution of the tribe Gochnatieae (Compositae) has been the subject of considerable effort in the past decade. This is due to the key position of this tribe in the phylogeny of the sunflower family and the corresponding implications for biogeographic and morphological evolution of Compositae. Previous studies have confirmed the monophyly of this tribe as well as most of the genera that belong to it. However, phylogenetic resolution of Gochnatieae at both the genus- and species-level has remained poor. A subset of new phylogenomic loci used in this study has proven effective and has improved phylogenetic resolution in this group. The results of this work demonstrate Gochnatieae is a well-supported clade comprised of nine genera (Anastraphia, Cnicothamnus, Cyclolepis, Gochnatia, Moquiniastrum, Nahuatlea, Pentaphorus, Richterago, Tehuasca). One recently described genus, Vickia, was not included in this study; but its placement in Gochnatieae as a tenth genus in the tribe is well-justified. The monospecific Cyclolepis, which had been circumscribed within the tribe since its inception but was subsequently removed and designated as incertae sedis since 2014, is also shown to belong to Gochnatieae. We confirmed the monophyletic Moquiniastrum with two well-supported subclades. Ancestral area reconstruction analyses show that Gochnatieae originated in Eastern South America about 53 my. Apparently, except for Cyclolepis and Richterago, the ancestors of the other genera of Gochnatieae originated about 44 my from an area that now corresponds to the central Andes. The presence of the genera in the Chaco phytogeographic province, central Chile, and Mexico-United States-Caribbean is a result of dispersal from the central Andes. The ancestral distribution of Moquiniastrum corresponds to a large area comprising Eastern South America and the current central Andes, about 32 my. Ancestral character state reconstruction that included four characters indicates several states associated with complex plant reproductive biology such as gynodioecy, gynomonoecy, and polygamodioecy are derived in Gochnatieae as are heterogamous capitula (in Moquiniastrum and Richterago), dimorphic and subdimorphic corollas (in Cnicothamnus, Moquiniastrum, and Richterago), and the presence of marginal female corollas (in Moquiniastrum and Richterago). Within Moquiniastrum, two subclades (Densicephalum and Polymorphum) exhibit divergent patterns of trait evolution associated with these reproductive characters which suggests this genus can serve as a model to understand the sexual system evolution in plants.

Target enrichment and extensive population sampling help untangle the recent, rapid radiation of Oenothera sect. Calylophus

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking non-coding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic and morphological diversity in the group. We found that the combination of exons and flanking non-coding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation.

Effectiveness of Two Universal Angiosperm Probe Sets Tested In Silico for Caryophyllids Taxa with Emphasis on Cacti Species

In angiosperms, huge advances in massive DNA sequencing technologies have impacted phylogenetic studies. Probe sets have been developed with the purpose of recovering hundreds of orthologous loci of targeted DNA sequences (TDS) across different plant lineages. We tested in silico the effectiveness of two universal probe sets in the whole available genomes of Caryophyllids, emphasizing phylogenetic issues in cacti species. A total of 870 TDS (517 TDS from Angiosperm v.1 and 353 from Angiosperms353) were individually tested in nine cacti species and Amaranthus hypochondriacus (external group) with ≥17 Gbp of available DNA data. The effectiveness was measured by the total number of orthologous loci recovered and their length, the percentage of loci discarded by paralogy, and the proportion of informative sites (PIS) in the alignments. The results showed that, on average, Angiosperms353 was better than Angiosperm v.1 for cacti species, since the former obtained an average of 275.6 loci that represent 123,687 bp, 2.48% of paralogous loci, and 4.32% of PIS in alignments, whereas the latter recovered 148.4 loci (37,683 bp), 10.38% of paralogous loci, and 3.49% of PIS. We recommend the use of predesigned universal probe sets for Caryophyllids, since these recover a high number of orthologous loci that resolve phylogenetic relationships.

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.].

Whence Came These Plants Most Foul? Phylogenomics and Biogeography of Lowiaceae (Zingiberales)

Lowiaceae (order Zingiberales) is a small family of forest herbs in Southeast Asia. All species belong to the genus Orchidantha. They are known for possessing orchid-like flowers that are smelly, apparently mimicking dead animals, feces, or mushrooms. Little is known of the biogeographic patterns or character evolution of the family. We sampled the family extensively, including many recently discovered species, and reconstructed the phylogeny of the family using HybSeq with Lowiaceae-specific RNA baits. Our phylogenetic reconstructions confirm that the family is most closely related to Strelitziaceae, and that species with dark, foul-smelling flowers form a grade in which a clade of species with paler flowers are embedded. The pale-flowered species produce a distinct odor, resembling edible mushrooms. Apart from a single species, the species from Borneo form a clade, and the same is true for Indochinese species. The remaining species form a more widespread clade. A biogeographic analysis shows that the distribution of Lowiaceae can explained by vicariance and gradual dispersal from a shared ancestral range of Borneo and Indochina. There is no evidence of long-distance dispersal, only a later extension in distribution to Peninsular Malaysia which coincides with the presence of a land bridge. Different directions of spread are possible, but none require long-distance dispersal. The results are consistent with the geological history of Southeast Asia. In particular, the relatively early isolation between Indochina and Borneo could be explained by the presence of a sea barrier that developed 10–15 MYA, and the continuous movement of plant species between Borneo and Peninsular Malaysia could be explained by a land bridge that existed until c. 5 MYA. The lack of an extensive land bridge with a suitable habitat may explain the absence of this genus from Sumatra and other Indonesian islands aside from Borneo. The strict reliance on a continuous habitat for the range expansion of Lowiaceae can be explained by their fruits and seeds, which lack obvious adaptations for long-distance dispersal. The inability to disperse to new areas may also explain why the extant species have very restricted distributions.

Analysis of Paralogs in Target Enrichment Data Pinpoints Multiple Ancient Polyploidy Events in Alchemilla s.l. (Rosaceae)

Target enrichment is becoming increasingly popular for phylogenomic studies. Although baits for enrichment are typically designed to target single-copy genes, paralogs are often recovered with increased sequencing depth, sometimes from a significant proportion of loci, especially in groups experiencing whole-genome duplication (WGD) events. Common approaches for processing paralogs in target enrichment data sets include random selection, manual pruning, and mainly, the removal of entire genes that show any evidence of paralogy. These approaches are prone to errors in orthology inference or removing large numbers of genes. By removing entire genes, valuable information that could be used to detect and place WGD events is discarded. Here, we used an automated approach for orthology inference in a target enrichment data set of 68 species of Alchemilla s.l. (Rosaceae), a widely distributed clade of plants primarily from temperate climate regions. Previous molecular phylogenetic studies and chromosome numbers both suggested ancient WGDs in the group. However, both the phylogenetic location and putative parental lineages of these WGD events remain unknown. By taking paralogs into consideration and inferring orthologs from target enrichment data, we identified four nodes in the backbone of Alchemilla s.l. with an elevated proportion of gene duplication. Furthermore, using a gene-tree reconciliation approach, we established the autopolyploid origin of the entire Alchemilla s.l. and the nested allopolyploid origin of four major clades within the group. Here, we showed the utility of automated tree-based orthology inference methods, previously designed for genomic or transcriptomic data sets, to study complex scenarios of polyploidy and reticulate evolution from target enrichment data sets.[Alchemilla; allopolyploidy; autopolyploidy; gene tree discordance; orthology inference; paralogs; Rosaceae; target enrichment; whole genome duplication.].

Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA

BACKGROUND AND AIMS

With the advance of high-throughput sequencing, reduced-representation methods such as target capture sequencing (TCS) emerged as cost-efficient ways of gathering genomic information, particularly from coding regions. As the off-target reads from such sequencing are expected to be similar to genome skimming (GS), we assessed the quality of repeat characterization in plant genomes using these data.

METHODS

Repeat composition obtained from TCS datasets of five Rhynchospora (Cyperaceae) species were compared with GS data from the same taxa. In addition, a FISH probe was designed based on the most abundant satellite found in the TCS dataset of Rhynchospora cephalotes. Finally, repeat-based phylogenies of the five Rhynchospora species were constructed based on the GS and TCS datasets and the topologies were compared with a gene-alignment-based phylogenetic tree.

KEY RESULTS

All the major repetitive DNA families were identified in TCS, including repeats that showed abundances as low as 0.01 % in the GS data. Rank correlations between GS and TCS repeat abundances were moderately high (r = 0.58-0.85), increasing after filtering out the targeted loci from the raw TCS reads (r = 0.66-0.92). Repeat data obtained by TCS were also reliable in developing a cytogenetic probe of a new variant of the holocentromeric satellite Tyba. Repeat-based phylogenies from TCS data were congruent with those obtained from GS data and the gene-alignment tree.

CONCLUSIONS

Our results show that off-target TCS reads can be recycled to identify repeats for cyto- and phylogenomic investigations. Given the growing availability of TCS reads, driven by global phylogenomic projects, our strategy represents a way to recycle genomic data and contribute to a better characterization of plant biodiversity.

Phylogeny, origin, and dispersal of Dubyaea (Asteraceae) based on Hyb-Seq data

Dubyaea DC. is a small genus of Asteraceae that is almost exclusively endemic to the Pan-Himalayan region. Within Dubyaea, phylogenetic relationships remain poorly understood. Here, our well-supported phylogeny based on Hyb-Seq data shows that all samples of Dubyaea in this study belong to a monophyletic group, which is sister to the clade of Soroseris, Syncalathium, and Nabalus. Dubyaea (s. str.) can be divided into three major clades, which are supported by flower color as well as morphological features of the stems and basal leaves. Based on our phylogenetic results, we performed biogeographic analyses and inferred that Dubyaea arose in the late Miocene in Hengduan Mountains and its eastern areas. Following its evolutionary origin, Dubyaea underwent diversification in situ as well as spread to the Himalayas.

Relative performance of customized and universal probe sets in target enrichment: A case study in subtribe Malinae

PREMISE

Custom probe design for target enrichment in phylogenetics is tedious and often hinders broader phylogenetic synthesis. The universal angiosperm probe set Angiosperms353 may be the solution. Here, we test the relative performance of Angiosperms353 on the Rosaceae subtribe Malinae in comparison with custom probes that we specifically designed for this clade. We then address the impact of bioinformatically altering the performance of Angiosperms353 by replacing the original probe sequences with orthologs extracted from the Malus domestica genome.

METHODS

To evaluate the relative performance of these probe sets, we compared the enrichment efficiency, locus recovery, alignment length, proportion of parsimony-informative sites, proportion of potential paralogs, the topology and support of the resulting species trees, and the gene tree discordance.

RESULTS

Locus recovery was highest for our custom Malinae probe set, and replacing the original Angiosperms353 sequences with a Malus representative improved the locus recovery relative to Angiosperms353. The proportion of parsimony-informative sites was similar between all probe sets, while the gene tree discordance was lower in the case of the custom probes.

DISCUSSION

A custom probe set benefits from data completeness and can be tailored toward the specificities of the project of choice; however, Angiosperms353 was equally as phylogenetically informative as the custom probes. We therefore recommend using both a custom probe set and Angiosperms353 to facilitate large-scale systematic studies, where financially possible.

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.

Phylotranscriptomic insights into Asteraceae diversity, polyploidy, and morphological innovation

Biodiversity is not evenly distributed among related groups, raising questions about the factors contributing to such disparities. The sunflower family (Asteraceae, >26,000 species) is among the largest and most diverse plant families, but its species diversity is concentrated in a few subfamilies, providing an opportunity to study the factors affecting biodiversity. Phylotranscriptomic analyses here of 244 transcriptomes and genomes produced a phylogeny with strong support for the monophyly of Asteraceae and the monophyly of most subfamilies and tribes. This phylogeny provides a reference for detecting changes in diversification rates and possible factors affecting Asteraceae diversity, which include global climate shifts, whole-genome duplications (WGDs), and morphological evolution. The origin of Asteraceae was estimated at ~83 Mya, with most subfamilies having diverged before the Cretaceous-Paleocene boundary. Phylotranscriptomic analyses supported the existence of 41 WGDs in Asteraceae. Changes to herbaceousness and capitulescence with multiple flower-like capitula, often with distinct florets and scaly pappus/receptacular bracts, are associated with multiple upshifts in diversification rate. WGDs might have contributed to the survival of early Asteraceae by providing new genetic materials to support morphological transitions. The resulting competitive advantage for adapting to different niches would have increased biodiversity in Asteraceae.

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.

Lineage-specific vs. universal: A comparison of the Compositae1061 and Angiosperms353 enrichment panels in the sunflower family

PREMISE

Phylogenetic studies in the Compositae are challenging due to the sheer size of the family and the challenges they pose for molecular tools, ranging from the genomic impact of polyploid events to their very conserved plastid genomes. The search for better molecular tools for phylogenetic studies led to the development of the family-specific Compositae1061 probe set, as well as the universal Angiosperms353 probe set designed for all flowering plants. In this study, we evaluate the extent to which data generated using the family-specific kit and those obtained with the universal kit can be merged for downstream analyses.

METHODS

We used comparative methods to verify the presence of shared loci between probe sets. Using two sets of eight samples sequenced with Compositae1061 and Angiosperms353, we ran phylogenetic analyses with and without loci flagged as paralogs, a gene tree discordance analysis, and a complementary phylogenetic analysis mixing samples from both sample sets.

RESULTS

Our results show that the Compositae1061 kit provides an average of 721 loci, with 9-46% of them presenting paralogs, while the Angiosperms353 set yields an average of 287 loci, which are less affected by paralogy. Analyses mixing samples from both sets showed that the presence of 30 shared loci in the probe sets allows the combination of data generated in different ways.

DISCUSSION

Combining data generated using different probe sets opens up the possibility of collaborative efforts and shared data within the synantherological community.

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.

Effects of ethyl acetate extract from Coreopsis tinctoria on learning and memory impairment in d-galactose-induced aging mice and the underlying molecular mechanism

The aim of this study was to investigate the effects of ethyl acetate extract from Coreopsis tinctoria (EACC) on learning and memory impairment in d-galactose-induced aging mice and the underlying molecular mechanism. The composition of EACC was analyzed by UPLC-MS, and the targets and pathways of EACC to improve learning and memory impairment were predicted and analyzed by the network pharmacology method. A mouse aging model was established by subcutaneous injection of d-galactose in mice, and EACC and piracetam were given to the model mice by gavage to observe their behavioral changes and changes in their SOD and GSH-Px activities in MDA contents in their peripheral blood serum and in the contents of Glu and GABA in their brain tissues. Then the hippocampus of the three mice selected from each of the MOD group and EACC-H group was separated for RT-qPCR assay. The results of the animal experiments showed that EACC could improve the learning and memory impairment of model mice by affecting the level of oxidative stress enzymes in serum and the content of neurotransmitters in the brain tissue. The results of network pharmacology analysis showed that the EACC components corresponded to 74 learning and memory-related targets, of which 13 were enriched in the long-term potentiation pathway. The results of RT-qPCR showed that 12 of the 13 detected targets were consistent with the predicted targets, and 9 of them were located in the NMDA receptor-related pathway of the long-term potentiation process and the pathway played an important regulatory role. It is believed that EACC could improve the learning and memory impairment of d-galactose-induced aging mice by acting on the nine targets Grin1, Grin2a, Camk2a, Camk2b, Kras, Raf1, Mapk1, Mapk3 and Creb to affect the NMDA receptor-related pathway of long-term potentiation.

High-throughput methods for efficiently building massive phylogenies from natural history collections

PREMISE

Large phylogenetic data sets have often been restricted to small numbers of loci from GenBank, and a vetted sampling-to-sequencing phylogenomic protocol scaling to thousands of species is not yet available. Here, we report a high-throughput collections-based approach that empowers researchers to explore more branches of the tree of life with numerous loci.

METHODS

We developed an integrated Specimen-to-Laboratory Information Management System (SLIMS), connecting sampling and wet lab efforts with progress tracking at each stage. Using unique identifiers encoded in QR codes and a taxonomic database, a research team can sample herbarium specimens, efficiently record the sampling event, and capture specimen images. After sampling in herbaria, images are uploaded to a citizen science platform for metadata generation, and tissue samples are moved through a simple, high-throughput, plate-based herbarium DNA extraction and sequencing protocol.

RESULTS

We applied this sampling-to-sequencing workflow to ~15,000 species, producing for the first time a data set with ~50% taxonomic representation of the "nitrogen-fixing clade" of angiosperms.

DISCUSSION

The approach we present is appropriate at any taxonomic scale and is extensible to other collection types. The widespread use of large-scale sampling strategies repositions herbaria as accessible but largely untapped resources for broad taxonomic sampling with thousands of species.

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.

How to Tackle Phylogenetic Discordance in Recent and Rapidly Radiating Groups? Developing a Workflow Using Loricaria (Asteraceae) as an Example

A major challenge in phylogenetics and -genomics is to resolve young rapidly radiating groups. The fast succession of species increases the probability of incomplete lineage sorting (ILS), and different topologies of the gene trees are expected, leading to gene tree discordance, i.e., not all gene trees represent the species tree. Phylogenetic discordance is common in phylogenomic datasets, and apart from ILS, additional sources include hybridization, whole-genome duplication, and methodological artifacts. Despite a high degree of gene tree discordance, species trees are often well supported and the sources of discordance are not further addressed in phylogenomic studies, which can eventually lead to incorrect phylogenetic hypotheses, especially in rapidly radiating groups. We chose the high-Andean Asteraceae genus Loricaria to shed light on the potential sources of phylogenetic discordance and generated a phylogenetic hypothesis. By accounting for paralogy during gene tree inference, we generated a species tree based on hundreds of nuclear loci, using Hyb-Seq, and a plastome phylogeny obtained from off-target reads during target enrichment. We observed a high degree of gene tree discordance, which we found implausible at first sight, because the genus did not show evidence of hybridization in previous studies. We used various phylogenomic analyses (trees and networks) as well as the D-statistics to test for ILS and hybridization, which we developed into a workflow on how to tackle phylogenetic discordance in recent radiations. We found strong evidence for ILS and hybridization within the genus Loricaria. Low genetic differentiation was evident between species located in different Andean cordilleras, which could be indicative of substantial introgression between populations, promoted during Pleistocene glaciations, when alpine habitats shifted creating opportunities for secondary contact and hybridization.

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.

Paralogs and off-target sequences improve phylogenetic resolution in a densely-sampled study of the breadfruit genus (Artocarpus, Moraceae)

We present a 517-gene phylogenetic framework for the breadfruit genus Artocarpus (ca. 70 spp., Moraceae), making use of silica-dried leaves from recent fieldwork and herbarium specimens (some up to 106 years old) to achieve 96% taxon sampling. We explore issues relating to assembly, paralogous loci, partitions, and analysis method to reconstruct a phylogeny that is robust to variation in data and available tools. While codon partitioning did not result in any substantial topological differences, the inclusion of flanking non-coding sequence in analyses significantly increased the resolution of gene trees. We also found that increasing the size of datasets increased convergence between analysis methods but did not reduce gene tree conflict. We optimized the HybPiper targeted-enrichment sequence assembly pipeline for short sequences derived from degraded DNA extracted from museum specimens. While the subgenera of Artocarpus were monophyletic, revision is required at finer scales, particularly with respect to widespread species. We expect our results to provide a basis for further studies in Artocarpus and provide guidelines for future analyses of datasets based on target enrichment data, particularly those using sequences from both fresh and museum material, counseling careful attention to the potential of off-target sequences to improve resolution.

A two-tier bioinformatic pipeline to develop probes for target capture of nuclear loci with applications in Melastomataceae

PREMISE

Putatively single-copy nuclear (SCN) loci, which are identified using genomic resources of closely related species, are ideal for phylogenomic inference. However, suitable genomic resources are not available for many clades, including Melastomataceae. We introduce a versatile approach to identify SCN loci for clades with few genomic resources and use it to develop probes for target enrichment in the distantly related Memecylon and Tibouchina (Melastomataceae).

METHODS

We present a two-tiered pipeline. First, we identified putatively SCN loci using MarkerMiner and transcriptomes from distantly related species in Melastomataceae. Published loci and genes of functional significance were then added (384 total loci). Second, using HybPiper, we retrieved 689 homologous template sequences for these loci using genome-skimming data from within the focal clades.

RESULTS

We sequenced 193 loci common to Memecylon and Tibouchina. Probes designed from 56 template sequences successfully targeted sequences in both clades. Probes designed from genome-skimming data within a focal clade were more successful than probes designed from other sources.

DISCUSSION

Our pipeline successfully identified and targeted SCN loci in Memecylon and Tibouchina, enabling phylogenomic studies in both clades and potentially across Melastomataceae. This pipeline could be easily applied to other clades with few genomic resources.

Strategies for reducing per-sample costs in target capture sequencing for phylogenomics and population genomics in plants

The reduced cost of high-throughput sequencing and the development of gene sets with wide phylogenetic applicability has led to the rise of sequence capture methods as a plausible platform for both phylogenomics and population genomics in plants. An important consideration in large targeted sequencing projects is the per-sample cost, which can be inflated when using off-the-shelf kits or reagents not purchased in bulk. Here, we discuss methods to reduce per-sample costs in high-throughput targeted sequencing projects. We review the minimal equipment and consumable requirements for targeted sequencing while comparing several alternatives to reduce bulk costs in DNA extraction, library preparation, target enrichment, and sequencing. We consider how each of the workflow alterations may be affected by DNA quality (e.g., fresh vs. herbarium tissue), genome size, and the phylogenetic scale of the project. We provide a cost calculator for researchers considering targeted sequencing to use when designing projects, and identify challenges for future development of low-cost sequencing in non-model plant systems.

Utility of targeted sequence capture for phylogenomics in rapid, recent angiosperm radiations: Neotropical Burmeistera bellflowers as a case study

Targeted sequence capture is a promising approach for large-scale phylogenomics. However, rapid evolutionary radiations pose significant challenges for phylogenetic inference (e.g. incomplete lineages sorting (ILS), phylogenetic noise), and the ability of targeted nuclear loci to resolve species trees despite such issues remains poorly studied. We test the utility of targeted sequence capture for inferring phylogenetic relationships in rapid, recent angiosperm radiations, focusing on Burmeistera bellflowers (Campanulaceae), which diversified into ~130 species over less than 3 million years. We compared phylogenies estimated from supercontig (exons plus flanking sequences), exon-only, and flanking-only datasets with 506-546 loci (~4.7 million bases) for 46 Burmeistera species/lineages and 10 outgroup taxa. Nuclear loci resolved backbone nodes and many congruent internal relationships with high support in concatenation and coalescent-based species tree analyses, and inferences were largely robust to effects of missing taxa and base composition biases. Nevertheless, species trees were incongruent between datasets, and gene trees exhibited remarkably high levels of conflict (~4-60% congruence, ~40-99% conflict) not simply driven by poor gene tree resolution. Higher gene tree heterogeneity at shorter branches suggests an important role of ILS, as expected for rapid radiations. Phylogenetic informativeness analyses also suggest this incongruence has resulted from low resolving power at short internal branches, consistent with ILS, and homoplasy at deeper nodes, with exons exhibiting much greater risk of incorrect topologies due to homoplasy than other datasets. Our findings suggest that targeted sequence capture is feasible for resolving rapid, recent angiosperm radiations, and that results based on supercontig alignments containing nuclear exons and flanking sequences have higher phylogenetic utility and accuracy than either alone. We use our results to make practical recommendations for future target capture-based studies of Burmeistera and other rapid angiosperm radiations, including that such studies should analyze supercontigs to maximize the phylogenetic information content of loci.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

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.

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.

A fully resolved backbone phylogeny reveals numerous dispersals and explosive diversifications throughout the history of Asteraceae

Flowering plant species represent at least 95% of all vascular plants on Earth, and members of the sunflower family comprise roughly 10% of this diversity. The family is often considered taxonomically difficult primarily because it is enormous in size and cosmopolitan in distribution. Using phylogenomics, we were able to fully resolve the backbone of the sunflower family tree. We provide evidence for a late Cretaceous origin followed by explosive diversifications and dispersals during the middle Eocene—ultimately resulting in the family’s 25,000+ extant species. Our results provide a framework to interpret the spatiotemporal patterns of migration out of South America and the family’s explosive diversifications out of Africa that led to its global evolutionary and ecological success.

The sunflower family, Asteraceae, comprises 10% of all flowering plant species and displays an incredible diversity of form. Asteraceae are clearly monophyletic, yet resolving phylogenetic relationships within the family has proven difficult, hindering our ability to understand its origin and diversification. Recent molecular clock dating has suggested a Cretaceous origin, but the lack of deep sampling of many genes and representative taxa from across the family has impeded the resolution of migration routes and diversifications that led to its global distribution and tremendous diversity. Here we use genomic data from 256 terminals to estimate evolutionary relationships, timing of diversification(s), and biogeographic patterns. Our study places the origin of Asteraceae at ∼83 MYA in the late Cretaceous and reveals that the family underwent a series of explosive radiations during the Eocene which were accompanied by accelerations in diversification rates. The lineages that gave rise to nearly 95% of extant species originated and began diversifying during the middle Eocene, coincident with the ensuing marked cooling during this period. Phylogenetic and biogeographic analyses support a South American origin of the family with subsequent dispersals into North America and then to Asia and Africa, later followed by multiple worldwide dispersals in many directions. The rapid mid-Eocene diversification is aligned with the biogeographic range shift to Africa where many of the modern-day tribes appear to have originated. Our robust phylogeny provides a framework for future studies aimed at understanding the role of the macroevolutionary patterns and processes that generated the enormous species diversity of Asteraceae.