Fluidigm2PURC: Automated processing and haplotype inference for double-barcoded PCR amplicons

PURC Provides Improved Sequence Inference for Polyploid Phylogenetics and Other Manifestations of the Multiple-Copy Problem

Inferring the true biological sequences from amplicon mixtures remains a difficult bioinformatics problem. The traditional approach is to cluster sequencing reads by similarity thresholds and treat the consensus sequence of each cluster as an "operational taxonomic unit" (OTU). Recently, this approach has been improved by model-based methods that correct PCR and sequencing errors in order to infer "amplicon sequence variants" (ASVs). To date, ASV approaches have been used primarily in metagenomics, but they are also useful for determining homeologs in polyploid organisms. To facilitate the usage of ASV methods among polyploidy researchers, we incorporated ASV inference alongside OTU clustering in PURC v2.0, a major update to PURC (Pipeline for Untangling Reticulate Complexes). In addition, PURC v2.0 features faster demultiplexing than the original version and has been updated to be compatible with Python 3. In this chapter we present results indicating that using the ASV approach is more likely to infer the correct biological sequences in comparison to the earlier OTU-based PURC and describe how to prepare sequencing data, run PURC v2.0 under several different modes, and interpret the output.

Adaptive responses drive the success of polyploid yellowcresses (Rorippa, Brassicaceae) in the Hengduan Mountains, a temperate biodiversity hotspot

Polyploids contribute substantially to plant evolution and biodiversity; however, the mechanisms by which they succeed are still unclear. According to the polyploid adaptation hypothesis, successful polyploids spread by repeated adaptive responses to new environments. Here, we tested this hypothesis using two tetraploid yellowcresses (Rorippa), the endemic Rorippa elata and the widespread Rorippa palustris, in the temperate biodiversity hotspot of the Hengduan Mountains. Speciation modes were resolved by phylogenetic modeling using 12 low-copy nuclear loci. Phylogeographical patterns were then examined using haplotypes phased from four plastid and ITS markers, coupled with historical niche reconstruction by ecological niche modeling. We inferred the time of hybrid origins for both species as the mid-Pleistocene, with shared glacial refugia within the southern Hengduan Mountains. Phylogeographic and ecological niche reconstruction indicated recurrent northward colonization by both species after speciation, possibly tracking denuded habitats created by glacial retreat during interglacial periods. Common garden experiment involving perennial R. elata conducted over two years revealed significant changes in fitness-related traits across source latitudes or altitudes, including latitudinal increases in survival rate and compactness of plant architecture, suggesting gradual adaptation during range expansion. These findings support the polyploid adaptation hypothesis and suggest that the spread of polyploids was aided by adaptive responses to environmental changes during the Pleistocene. Our results thus provide insight into the evolutionary success of polyploids in high-altitude environments.

Phylogeny, classification, and character evolution of tribe Citharexyleae (Verbenaceae)

PREMISE

As a family of Neotropical origin and primarily Neotropical distribution, the Verbenaceae are a good but understudied system with which to understand Neotropical evolution. Tribe Citharexyleae comprises three genera: Baillonia, Citharexylum-one of the largest genera in Verbenaceae-and Rehdera. A molecular phylogenetic approach was taken to resolve intergeneric relationships in Citharexyleae and infrageneric relationships in Citharexylum. The phylogeny is used to elucidate character evolution in a widespread, morphologically diverse Neotropical genus.

METHODS

Seven plastid regions, two nuclear ribosomal spacers, and six low-copy nuclear loci were analyzed for 64 species of Citharexyleae. Phylogenetic analyses were conducted using maximum likelihood, Bayesian inference, and multispecies coalescent approaches. Habit, presence or absence of thorns, inflorescence architecture, flower color, fruit color, and geography were examined to identify diagnostic character states for clades within Citharexylum.

RESULTS

Rehdera is resolved as sister to Citharexylum, and Baillonia nested within Citharexylum. Two species, C. oleinum and C. tetramerum, are not closely related to tribe Citharexyleae, but may be related to members of tribe Duranteae instead. Seven clades within Citharexylum are inferred, each characterized by a combination of geography, fruit color and/or maturation, and inflorescence architecture. There is evidence of correlated evolution between habit, axillary inflorescences, and flower number per inflorescence. Shrubs with reduced inflorescences have evolved repeatedly.

CONCLUSIONS

A subgeneric classification for Citharexylum is proposed. Although suites of associated traits are found, character morphology has been labile throughout Citharexylum's evolutionary history. Morphological diversity may be related to adaptation to differing mesic and xeric habitats.

Unveiling the Patterns of Reticulated Evolutionary Processes with Phylogenomics: Hybridization and Polyploidy in the genus Rosa

Reticulation, caused by hybridization and allopolyploidization, is considered an important and frequent phenomenon in the evolution of numerous plant lineages. Although both processes represent important driving forces of evolution, they are mostly ignored in phylogenetic studies involving a large number of species. Indeed only a scattering of methods exists to recover a comprehensive reticulated evolutionary history for a broad taxon sampling. Among these methods, comparisons of topologies obtained from plastid markers with those from a few nuclear sequences are favored, even though they restrict in-depth studies of hybridization and polyploidization. The genus Rosa encompasses c. 150 species widely distributed throughout the northern hemisphere and represents a challenging taxonomic group in which hybridization and polyploidization are prominent. Our main objective was to develop a general framework that would take patterns of reticulation into account in the study of the phylogenetic relationships among Rosa species. Using amplicon sequencing we targeted allele variation in the nuclear genome as well as haploid sequences in the chloroplast genome. We successfully recovered robust plastid and nuclear phylogenies and performed in-depth tests for several scenarios of hybridization using a maximum pseudo-likelihood approach on taxon subsets. Our diploid-first approach followed by hybrid and polyploid grafting resolved most of the evolutionary relationships among Rosa subgenera, sections, and selected species. Based on these results, we provide new directions for a future revision of the infrageneric classification in Rosa. The stepwise strategy proposed here can be used to reconstruct the phylogenetic relationships of other challenging taxonomic groups with large numbers of hybrid and polyploid taxa.

Polyploid phylogenetics

Polyploidy is a dominant feature of extant plant diversity. However, major research questions, including whether polyploidy is important to long-term evolution or is just 'evolutionary noise', remain unresolved due to difficulties associated with the generation and analysis of data from polyploid lineages. Many of these difficulties have been recently overcome, such that it is now often relatively straightforward to infer the full and often reticulate phylogenetic history of groups with recently formed polyploids. This nascent field of 'polyploid phylogenetics' allows researchers to tackle long-standing questions of polyploid macroevolution, supplies the foundation for mechanistic models of ploidy change, and provides the opportunity to include a more complete and representative sample of plant taxa in our analyses in general.

Examining the utility of DNA barcodes for the identification of tallgrass prairie flora

PREMISE

The tallgrass prairies of North America are one of the most threatened ecosystems in the world, making efficient species identification essential for understanding and managing diversity. Here, we assess DNA barcoding with high-throughput sequencing as a method for rapid plant species identification.

METHODS

Using herbarium collections representing the tallgrass prairie flora of Oak Lake Field Station, South Dakota, USA, we amplified and examined four common nuclear and plastid barcode regions (ITS, matK, psbA-trnH, and rbcL), individually and in combination, to test their success in identifying samples to family, genus, and species levels using BLAST searches of three databases of varying size.

RESULTS

Concatenated barcodes increased performance, although none were significantly different than single-region barcodes. The plastid region psbA-trnH performed significantly more poorly than the others, while barcodes containing ITS performed best. Database size significantly affected identification success at all three taxonomic levels. Confident species-level identification ranged from 8-44% for the global database, 13-56% for the regional database, and 21-80% for the sampled species database, depending on the barcode used.

DISCUSSION

Barcoding was generally successful in identifying tallgrass prairie genera and families, but was of limited use in species-level identifications. Database size was an important factor in successful plant identification. We discuss future directions and considerations for improving the performance of DNA barcoding in tallgrass prairies.

Diet of a rare herbivore based on DNA metabarcoding of feces: Selection, seasonality, and survival

In herbivores, survival and reproduction are influenced by quality and quantity of forage, and hence, diet and foraging behavior are the foundation of an herbivore's life history strategy. Given the importance of diet to most herbivores, it is imperative that we know the species of plants they prefer, especially for herbivorous species that are at risk for extinction. However, it is often difficult to identify the diet of small herbivores because: (a) They are difficult to observe, (b) collecting stomach contents requires sacrificing animals, and (c) microhistology requires accurately identifying taxa from partially digested plant fragments and likely overemphasizes less-digestible taxa. The northern Idaho ground squirrel (Urocitellus brunneus) is federally threatened in the United States under the Endangered Species Act. We used DNA metabarcoding techniques to identify the diet of 188 squirrels at 11 study sites from fecal samples. We identified 42 families, 126 genera, and 120 species of plants in the squirrel's diet. Our use of three gene regions was beneficial because reliance on only one gene region (e.g., only trnL) would have caused us to miss >30% of the taxa in their diet. Northern Idaho ground squirrel diet differed between spring and summer, frequency of many plants in the diet differed from their frequency within their foraging areas (evidence of selective foraging), and several plant genera in their diet were associated with survival. Our results suggest that while these squirrels are generalists (they consume a wide variety of plant species), they are also selective and do not eat plants relative to availability. Consumption of particular genera such as Perideridia may be associated with higher overwinter survival.

Phylogeny of Paullinia L. (Paullinieae: Sapindaceae), a diverse genus of lianas with dynamic fruit evolution

Paullinia L. is a genus of c. 220 mostly Neotropical forest-dwelling lianas that display a wide diversity of fruit morphologies. Paullinia resembles other members of the Paullinieae tribe in being a climber with stipulate compound leaves and paired inflorescence tendrils. However, it is distinct in having capsular fruits with woody, coriaceous, or crustaceous pericarps. While consistent in this basic plan, the pericarps of Paullinia fruits are otherwise highly variable-in some species they are winged, whereas in others they are without wings or covered with spines. With the exception of the water-dispersed indehiscent spiny fruits of some members of Paullinia sect. Castanella, all species are dehiscent, opening their capsules while they are still attached to the branch, to reveal arillate animal-dispersed seeds. Here we present a molecular phylogeny of Paullinia derived from 11 molecular markers, including nine newly-developed single-copy nuclear markers amplified by microfluidics PCR. This is the first broadly sampled molecular phylogeny for the genus. Paullinia is supported as monophyletic and is sister to Cardiospermum L., which together are sister to Serjania Mill + Urvillea Kunth. We apply this novel phylogenetic hypothesis to test previous infrageneric classifications and to infer that unwinged fruits represent the ancestral condition, from which there were repeated evolutionary transitions and reversals. However, because the seeds of both winged and unwinged fruits are dispersed by animals, we conclude that the repeated transitions in fruit morphology may relate to visual display strategies to attract animal dispersers, and do not represent transitions to wind dispersal.

Extensive allopolyploidy in the neotropical genus Lachemilla (Rosaceae) revealed by PCR-based target enrichment of the nuclear ribosomal DNA cistron and plastid phylogenomics

PREMISE OF THE STUDY

Polyploidy has been long recognized as an important force in plant evolution. Previous studies had suggested widespread occurrence of polyploidy and the allopolyploid origin of several species in the diverse neotropical genus Lachemilla (Rosaceae). Nonetheless, this evidence has relied mostly on patterns of cytonuclear discordance, and direct evidence from nuclear allelic markers is still needed.

METHODS

Here we used PCR target enrichment in combination with high throughput sequencing to obtain multiple copies of the nuclear ribosomal (nr) DNA cistron and 45 regions of the plastid genome (cpDNA) from 219 accessions representing 48 species of Lachemilla and to explore the allopolyploid origin of species in this group.

KEY RESULTS

We were able to identify multiple nrDNA ribotypes and establish clear evidence of allopolyploidy in 33 species of Lachemilla, showing that this condition is common and widespread in the genus. Additionally, we found evidence for three autopolyploid species. We also established multiple, independent origins of several allopolyploid species. Finally, based solely on the cpDNA phylogeny, we identified that the monotypic genus Farinopsis is the sister group of Lachemilla and allied genera within subtribe Fragariinae.

CONCLUSIONS

Our study demonstrates the utility of the nuclear ribosomal DNA cistron to detect allopolyploidy when concerted evolution of this region is not complete. Additionally, with a robust chloroplast phylogeny in place, the direction of hybridization events can be established, and multiple, independent origins of allopolyploid species can be identified.