A roadmap of phylogenomic methods for studying polyploid plant genera.

Phylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole-genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid-rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long-read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

Travel Tales of a Worldwide Weed: Genomic Signatures of Plantago major L. Reveal Distinct Genotypic Groups With Links to Colonial Trade Routes.

Retracing pathways of historical species introductions is fundamental to understanding the factors involved in the successful colonization and spread, centuries after a species' establishment in an introduced range. Numerous plants have been introduced to regions outside their native ranges both intentionally and accidentally by European voyagers and early colonists making transoceanic journeys; however, records are scarce to document this. We use genotyping-by-sequencing and genotype-likelihood methods on the selfing, global weed, Plantago major, collected from 50 populations worldwide to investigate how patterns of genomic diversity are distributed among populations of this global weed. Although genomic differentiation among populations is found to be low, we identify six unique genotype groups showing very little sign of admixture and low degree of outcrossing among them. We show that genotype groups are latitudinally restricted, and that more than one successful genotype colonized and spread into the introduced ranges. With the exception of New Zealand, only one genotype group is present in the Southern Hemisphere. Three of the most prevalent genotypes present in the native Eurasian range gave rise to introduced populations in the Americas, Africa, Australia, and New Zealand, which could lend support to the hypothesis that P. major was unknowlingly dispersed by early European colonists. Dispersal of multiple successful genotypes is a likely reason for success. Genomic signatures and phylogeographic methods can provide new perspectives on the drivers behind the historic introductions and the successful colonization of introduced species, contributing to our understanding of the role of genomic variation for successful establishment of introduced taxa.

Pioneering polyploids: the impact of whole-genome duplication on biome shifting in New Zealand Coprosma (Rubiaceae) and Veronica (Plantaginaceae).

The role of whole-genome duplication (WGD) in facilitating shifts into novel biomes remains unknown. Focusing on two diverse woody plant groups in New Zealand, Coprosma (Rubiaceae) and Veronica (Plantaginaceae), we investigate how biome occupancy varies with ploidy level, and test the hypothesis that WGD increases the rate of biome shifting. Ploidy levels and biome occupancy (forest, open and alpine) were determined for indigenous species in both clades. The distribution of low-ploidy (Coprosma: 2x, Veronica: 6x) versus high-ploidy (Coprosma: 4-10x, Veronica: 12-18x) species across biomes was tested statistically. Estimation of the phylogenetic history of biome occupancy and WGD was performed using time-calibrated phylogenies and the R package BioGeoBEARS. Trait-dependent dispersal models were implemented to determine support for an increased rate of biome shifting among high-ploidy lineages. We find support for a greater than random portion of high-ploidy species occupying multiple biomes. We also find strong support for high-ploidy lineages showing a three- to eightfold increase in the rate of biome shifts. These results suggest that WGD promotes ecological expansion into new biomes.

Using target sequence capture to improve the phylogenetic resolution of a rapid radiation in New Zealand Veronica.

PREMISE: Recent, rapid radiations present a challenge for phylogenetic reconstruction. Fast successive speciation events typically lead to low sequence divergence and poorly resolved relationships with standard phylogenetic markers. Target sequence capture of many independent nuclear loci has the potential to improve phylogenetic resolution for rapid radiations. METHODS: Here we applied target sequence capture with 353 protein-coding genes (Angiosperms353 bait kit) to Veronica sect. Hebe (common name hebe) to determine its utility for improving the phylogenetic resolution of rapid radiations. Veronica section Hebe originated 5-10 million years ago in New Zealand, forming a monophyletic radiation of ca 130 extant species. RESULTS: We obtained approximately 150 kbp of 353 protein-coding exons and an additional 200 kbp of flanking noncoding sequences for each of 77 hebe and two outgroup species. When comparing coding, noncoding, and combined data sets, we found that the latter provided the best overall phylogenetic resolution. While some deep nodes in the radiation remained unresolved, our phylogeny provided broad and often improved support for subclades identified by both morphology and standard markers in previous studies. Gene-tree discordance was nonetheless widespread, indicating that additional methods are needed to disentangle fully the history of the radiation. CONCLUSIONS: Phylogenomic target capture data sets both increase phylogenetic signal and deliver new insights into the complex evolutionary history of rapid radiations as compared with traditional markers. Improving methods to resolve remaining discordance among loci from target sequence capture is now important to facilitate the further study of rapid radiations.

Polyploidy on Islands: Its Emergence and Importance for Diversification.

Whole genome duplication or polyploidy is widespread among floras globally, but traditionally has been thought to have played a minor role in the evolution of island biodiversity, based on the low proportion of polyploid taxa present. We investigate five island systems (Juan Fernández, Galápagos, Canary Islands, Hawaiian Islands, and New Zealand) to test whether polyploidy (i) enhances or hinders diversification on islands and (ii) is an intrinsic feature of a lineage or an attribute that emerges in island environments. These island systems are diverse in their origins, geographic and latitudinal distributions, levels of plant species endemism (37% in the Galapagos to 88% in the Hawaiian Islands), and ploidy levels, and taken together are representative of islands more generally. We compiled data for vascular plants and summarized information for each genus on each island system, including the total number of species (native and endemic), generic endemicity, chromosome numbers, genome size, and ploidy levels. Dated phylogenies were used to infer lineage age, number of colonization events, and change in ploidy level relative to the non-island sister lineage. Using phylogenetic path analysis, we then tested how the diversification of endemic lineages varied with the direct and indirect effects of polyploidy (presence of polyploidy, time on island, polyploidization near colonization, colonizer pool size) and other lineage traits not associated with polyploidy (time on island, colonizer pool size, repeat colonization). Diploid and tetraploid were the most common ploidy levels across all islands, with the highest ploidy levels (>8x) recorded for the Canary Islands (12x) and New Zealand (20x). Overall, we found that endemic diversification of our focal island floras was shaped by polyploidy in many cases and certainly others still to be detected considering the lack of data in many lineages. Polyploid speciation on the islands was enhanced by a larger source of potential congeneric colonists and a change in ploidy level compared to overseas sister taxa.

The application of high-throughput sequencing for taxonomy: The case of Plantago subg. Plantago (Plantaginaceae).

Plantago is a cosmopolitan genus including over 250 species, concentrated in temperate and high-elevation tropical regions. The taxonomy of Plantago is very difficult, mainly because of its reduced morphology, which features relatively few characters for species classification. Consequently, the infrageneric classification of the genus remains controversial and inadequate. In this study we applied high-throughput plastid genome skimming to provide powerful phylogenetic resolution to clarify the relationships within subg. Plantago, which is the largest, most broadly distributed and poorest understood subgenus of Plantago. Ninety-four samples covering ~56% of all species and representing all sections of subg. Plantago as well as an outgroup were successfully sequenced. The resulting phylogenetic topology was used, complemented by field and herbarium studies, to revise the sectional classification of subg. Plantago and present a complete listing of the accepted species in the subgenus. Our phylogenetic results were also tested for their usefulness in clarifying the taxonomic placement of some taxonomically complicated species in the subgenus. We conclude that a combination of morphological studies and state-of-the art high-throughput DNA data provide a useful toolbox for resolving outstanding taxonomic puzzles exemplified by the genus Plantago.

Bayesian inference of phylogenetic networks from bi-allelic genetic markers.

Phylogenetic networks are rooted, directed, acyclic graphs that model reticulate evolutionary histories. Recently, statistical methods were devised for inferring such networks from either gene tree estimates or the sequence alignments of multiple unlinked loci. Bi-allelic markers, most notably single nucleotide polymorphisms (SNPs) and amplified fragment length polymorphisms (AFLPs), provide a powerful source of genome-wide data. In a recent paper, a method called SNAPP was introduced for statistical inference of species trees from unlinked bi-allelic markers. The generative process assumed by the method combined both a model of evolution for the bi-allelic markers, as well as the multispecies coalescent. A novel component of the method was a polynomial-time algorithm for exact computation of the likelihood of a fixed species tree via integration over all possible gene trees for a given marker. Here we report on a method for Bayesian inference of phylogenetic networks from bi-allelic markers. Our method significantly extends the algorithm for exact computation of phylogenetic network likelihood via integration over all possible gene trees. Unlike the case of species trees, the algorithm is no longer polynomial-time on all instances of phylogenetic networks. Furthermore, the method utilizes a reversible-jump MCMC technique to sample the posterior of phylogenetic networks given bi-allelic marker data. Our method has a very good performance in terms of accuracy and robustness as we demonstrate on simulated data, as well as a data set of multiple New Zealand species of the plant genus Ourisia (Plantaginaceae). We implemented the method in the publicly available, open-source PhyloNet software package.

Transcriptomic resources and marker validation for diploid and polyploid Veronica (Plantaginaceae) from New Zealand and Europe.

PREMISE OF THE STUDY: Polyploidy may generate novel variation, leading to adaptation and species diversification. An excellent natural system to study polyploid evolution in a comparative framework is Veronica (Plantaginaceae), which comprises several parallel, recently evolved polyploid series. METHODS: Over 105 million Illumina paired-end sequence reads were generated from cDNA libraries of leaf tissue from eight individuals representing three European and four New Zealand species. Forty-eight simple sequence repeat (SSR) and 48 low-copy nuclear (LCN) markers were developed and validated with Fluidigm microfluidic PCR and Illumina MiSeq amplicon sequencing on 48 different individuals each. RESULTS: Individual Trinity assemblies were similar regarding annotated transcripts (13,009-14,271), mean contig length (635-742 bp), N50 value (916-1133 bp), E90N50 value (1099-1308 bp), contigs with positive BLAST hits (42-63%), and gene ontology terms. Analyses of 29,738 single-nucleotide polymorphisms (8746 phylogenetically informative) mined from these transcriptomes plus two outgroups (Picrorhiza kurrooa and Plantago ovata) showed moderate to high bootstrap support for all branches and reticulation among sampled European Veronica. DISCUSSION: The transcriptome sequences themselves, as well as the validated SSR (40/48) and LCN (11/48) markers derived from them, show inter- and intraspecific genetic variation. These resources will be invaluable for future population genetic, phylogenetic, and functional genetic investigations in polyploid Veronica.

Microsatellite markers for the New Zealand endemic Myosotis pygmaea species group (Boraginaceae) amplify across species.

PREMISE OF THE STUDY: Microsatellite loci were developed as polymorphic markers for the New Zealand endemic Myosotis pygmaea species group (Boraginaceae) for use in species delimitation and population and conservation genetic studies. METHODS AND RESULTS: Illumina MiSeq sequencing was performed on genomic DNA from seedlings of M. drucei. From trimmed paired-end sequences >400 bp, 484 microsatellite loci were identified. Twelve of 48 microsatellite loci tested were found to be polymorphic and consistently scorable when screened on 53 individuals from four populations representing the geographic range of M. drucei. They also amplify in all other species in the M. pygmaea species group, i.e., M. antarctica, M. brevis, M. glauca, and M. pygmaea, as well as 18 other Myosotis species. CONCLUSIONS: These 12 polymorphic microsatellite markers establish an important resource for research and conservation of the M. pygmaea species group and potentially other Southern Hemisphere Myosotis.

Cost-effective fluorescent amplified fragment length polymorphism (AFLP) analyses using a three primer system.

The amplified fragment length polymorphism (AFLP) technique is a widely used multi-purpose DNA fingerprinting tool. The ability to size-separate fluorescently labelled AFLP fragments on a capillary electrophoresis instrument has provided a means for high-throughput genome screening, an approach particularly useful in studying the molecular ecology of nonmodel organisms. While the 'per-marker-generated' costs for AFLP are low, fluorescently labelled oligonucleotides remain costly. We present a cost-effective method for fluorescently end-labelling AFLPs that should make this tool more readily accessible for laboratories with limited budgets. Both standard fluorescent AFLPs and the end-labelled alternatives presented here are repeatable and produce similar numbers of fragments when scored using both manual and automated scoring methods. While it is not recommended to combine data using the two approaches, the results of the methods are qualitatively comparable, indicating that AFLP end-labelling is a robust alternative to standard methods of AFLP genotyping. For researchers commencing a new AFLP project, the AFLP end-labelling method outlined here is easily implemented, as it does not require major changes to PCR protocols and can significantly reduce the costs of AFLP studies.

First phylogenetic and biogeographical study of the southern bluebells (Wahlenbergia, Campanulaceae).

Wahlenbergia is a largely southern hemisphere genus of at least 260 species; within Campanulaceae only Campanula is larger. This first phylogeny of Wahlenbergia was reconstructed using about 20% of the 260 species in the genus based on the nuclear ribosomal ITS marker and the chloroplast trnL-F marker with samples from South Africa, Europe, Australia and New Zealand. Wahlenbergia was confirmed to be non-monophyletic, though most of the species form a clade. Our tree topology and date estimates indicate that Wahlenbergia diverged in South Africa about 29.6 mya, then dispersed to Australasia about 4.8 mya, thus indicating the radiation of Wahlenbergia occurred relatively recently. Radiations occurred in both of these main centres; there are currently about 170 species in South Africa and 45 species and subspecies in Australasia. New Zealand species comprise two clades, both rooted within the Australasian clade. We thus propose two dispersals from Australia to New Zealand, one leading to a radiation of species with the rhizomatous herbaceous growth form ca. 1.6 mya, and the other leading to a radiation of species with the radicate growth form 0.7 mya. Dispersals from Australia to New Zealand match the expected direction, following the west wind drift and ocean currents. The herbaceous growth form was shown to be ancestral for the genus as a whole, and polyploidy has been a mechanism of the evolution of the genus in Australasia.

Phylogeny of Veronica in the Southern and Northern Hemispheres based on plastid, nuclear ribosomal and nuclear low-copy DNA.

The cosmopolitan and ecologically diverse genus Veronica with approximately 450 species is the largest genus of the newly circumscribed Plantaginaceae. Previous analyses of Veronica DNA sequences were in stark contrast to traditional systematics. However, analyses did not allow many inferences regarding the relationship between major groups identified, hindering further analysis of diversification and evolutionary trends in the genus. To resolve the backbone relationships of Veronica, we added sequences from additional plastid DNA regions to existing data and analyzed matching data sets for 78 taxa and more than 5000 aligned characters from nuclear ribosomal DNA and plastid DNA regions. The results provide the best resolved and supported estimate of relationships among major groups in the Northern (Veronica s. str.) and Southern Hemisphere (hebes). We present new informal names for the five main species groups within the Southern Hemisphere sect. Hebe. Furthermore, in two instances we provide morphological and karyological characters supporting these relationships. Finally, we present the first evidence from nuclear low-copy CYCLOIDEA2-region to compare results from the plastid genome with the nuclear genome.

Species delimitation and phylogeny of a New Zealand plant species radiation.

BACKGROUND: Delimiting species boundaries and reconstructing the evolutionary relationships of late Tertiary and Quaternary species radiations is difficult. One recent approach emphasizes the use of genome-wide molecular markers, such as amplified fragment length polymorphisms (AFLPs) and single nucleotide polymorphisms (SNPs), to identify distinct metapopulation lineages as taxonomic species. Here we investigate the properties of AFLP data, and the usefulness of tree-based and non-tree-based clustering methods to delimit species and reconstruct evolutionary relationships among high-elevation Ourisia species (Plantaginaceae) in the New Zealand archipelago. RESULTS: New Zealand Ourisia are shown to comprise a geologically recent species radiation based on molecular dating analyses of ITS sequences (0.4-1.3 MY). Supernetwork analyses indicate that separate tree-based clustering analyses of four independent AFLP primer combinations and 193 individuals of Ourisia produced similar trees. When combined and analysed using tree building methods, 15 distinct metapopulations could be identified. These clusters corresponded very closely to species and subspecies identified on the basis of diagnostic morphological characters. In contrast, Structure and PCO-MC analyses of the same data identified a maximum of 12 and 8 metapopulations, respectively. All approaches resolved a large-leaved group and a small-leaved group, as well as a lineage of three alpine species within the small-leaved group. We were unable to further resolve relationships within these groups as corrected and uncorrected distances derived from AFLP profiles had limited tree-like properties. CONCLUSION: Ourisia radiated into a range of alpine and subalpine habitats in New Zealand during the Pleistocene, resulting in 13 morphologically and ecologically distinct species, including one reinstated from subspecies rank. Analyses of AFLP identified distinct metapopulations consistent with morphological characters allowing species boundaries to be delimited in Ourisia. Importantly, Structure analyses suggest some degree of admixture with most species, which may also explain why the AFLP data do not exhibit sufficient tree-like properties necessary for reconstructing some species relationships. We discuss this feature and highlight the importance of improving models for phylogenetic analyses of species radiations using AFLP and SNP data.

Optimizing automated AFLP scoring parameters to improve phylogenetic resolution.

The amplified fragment length polymorphism (AFLP) technique is an increasingly popular component of the phylogenetic toolbox, particularly for plant species. Technological advances in capillary electrophoresis now allow very precise estimates of DNA fragment mobility and amplitude, and current AFLP software allows greater control of data scoring and the production of the binary character matrix. However, for AFLP to become a useful modern tool for large data sets, improvements to automated scoring are required. We design a procedure that can be used to optimize AFLP scoring parameters to improve phylogenetic resolution and demonstrate it for two AFLP scoring programs (GeneMapper and GeneMarker). In general, we found that there was a trade-off between getting more characters of lower quality and fewer characters of high quality. Conservative settings that gave the least error did not give the best phylogenetic resolution, as too many useful characters were discarded. For example, in GeneMapper, we found that bin width was a crucial parameter, and that although reducing bin width from 1.0 to 0.5 base pairs increased the error rate, it nevertheless improved resolution due to the increased number of informative characters. For our 30-taxon data sets, moving from default to optimized parameter settings gave between 3 and 11 extra internal edges with >50% bootstrap support, in the best case increasing the number of resolved edges from 14 to 25 out of a possible 27. Nevertheless, improvements to current AFLP software packages are needed to (1) make use of replicate profiles to calibrate the data and perform error calculations and (2) perform tests to optimize scoring parameters in a rigorous and automated way. This is true not only when AFLP data are used for phylogenetics, but also for other applications, including linkage mapping and population genetics.

Phylogeographic patterns in the Australasian genus Chionohebe (Veronica s.l., Plantaginaceae) based on AFLP and chloroplast DNA sequences.

The alpine genus Chionohebe is one of seven genera in the southern hemisphere Hebe complex. The main aims of this study were to infer the evolutionary relationships and assess phylogeographic patterns among the six species of Chionohebe, determine the origin of the two species with trans-Tasman distributions, and test species delimitations and specimen identifications based on morphology. Analyses of AFLP data recovered five major lineages within Chionohebe, some of which corresponded to species and varieties as currently circumscribed. Although the cushion chionohebes were strongly supported as monophyletic, the sole non-cushion species, C. densifolia, was sister to Parahebe trifida, and thus the AFLP data do not support a monophyletic Chionohebe as usually circumscribed. Strong north/south and west/east phylogeographic patterns were found among and within the main AFLP lineages in New Zealand. Analyses of chloroplast DNA (cpDNA) revealed eight haplotypes in Chionohebe, but these did not correspond to current taxonomy or geography due to widespread interspecific haplotype sharing. Based on both AFLP and cpDNA results, the two trans-Tasman species are shown to have originated in New Zealand and dispersed to Australia independently.

Almost forgotten or latest practice? AFLP applications, analyses and advances.

Amplified fragment length polymorphism (AFLP) DNA fingerprinting is a firmly established molecular marker technique, with broad applications in population genetics, shallow phylogenetics, linkage mapping, parentage analyses, and single-locus PCR marker development. Technical advances have presented new opportunities for data analysis, and recent studies have addressed specific areas of the AFLP technique, including comparison to other genotyping methods, assessment of errors, homoplasy, phylogenetic signal and appropriate analysis techniques. Here we provide a synthesis of these areas and explore new directions for the AFLP technique in the genomic era, with the aim of providing a review that will be applicable to all AFLP-based studies.