Search

Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole-genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non-monophyly in Primula. Biogeographic analyses support the 'out-of-the-Caucasus' hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome-wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species-rich clades in biodiversity hotspots.

Obligate pollination mutualisms, in which plant and pollinator lineages depend on each other for reproduction, often exhibit high levels of species-specificity. However, cases in which two or more pollinator species share a single host species (host sharing), or two or more host species share a single pollinator species (pollinator sharing), are known to occur in current ecological time. Further, evidence for host switching in evolutionary time is increasingly being recognized in these systems. The degree to which departures from strict specificity differentially affect the potential for hybridization and introgression in the associated host or pollinator is unclear. We addressed this question using genome-wide sequence data from five sympatric Panamanian free-standing fig species (Ficus subgenus Pharmacosycea, section Pharmacosycea) and their six associated fig pollinator wasp species (Tetrapus). Two of the five fig species, F. glabrata and F. maxima, were found to regularly share pollinators. In these species, ongoing hybridization was demonstrated by the detection of several first-generation (F1) hybrid individuals and historical introgression was indicated by phylogenetic network analysis. In contrast, although two of the pollinator species regularly share hosts, all six species were genetically distinct and deeply divergent, with no evidence for either hybridization or introgression. This pattern is consistent with results from other obligate pollination mutualisms, suggesting that, in contrast to their host plants, pollinators appear to be reproductively isolated, even when different species of pollinators mate in shared hosts.

Premise The Ocotea complex contains the greatest diversity of Lauraceae in the Neotropics. However, the traditional taxonomy of the group has relied on only three main floral characters, and previous molecular analyses have used only a few markers and provided limited support for relationships among the major clades. This lack of useful data has hindered the development of a comprehensive classification, as well as studies of character evolution. Methods We used RAD-seq data to infer the phylogenetic relationships of 149 species in the Ocotea complex, generating a reference-based assembly using the Persea americana genome. The results provide the basis for a phylogenetic classification that reflects our current molecular knowledge and for analyses of the evolution of breeding system, stamen number, and number of anther locules. Results We recovered a well-supported tree that demonstrates the paraphyly of Licaria, Aniba, and Ocotea and clarifies the relationships of Umbellularia, Phyllostemonodaphne, and the Old World species. To begin the development of a new classification and to facilitate precise communication, we also provide phylogenetic definitions for seven major clades. Our ancestral reconstructions show multiple origins for the three floral characters that have routinely been used in Lauraceae systematics, suggesting that these be used with caution in the future. Conclusions This study advances our understanding of phylogenetic relationships and character evolution in a taxonomically difficult group using RAD-seq data. Our new phylogenetic names will facilitate unambiguous communication as studies of the Ocotea complex progress.

Genomes are composed of a mosaic of segments inherited from different ancestors, each separated by past recombination events. Consequently, genealogical relationships among multiple genomes vary spatially across different genomic regions. Expectations for the amount of genealogical variation among unlinked (uncorrelated) genomic regions is well described for either a single population (coalescent) or multiple structured populations (multispecies coalescent). However, the expected similarity among genealogies at linked regions of a genome is less well characterized. Recently, an analytical solution was derived for the expected distribution of waiting distances between changes in genealogical trees spatially across a genome for a single population with constant effective population size. Here we describe a generalization of this result, in terms of the expected distribution of waiting distances between changes in genealogical trees and topologies, for multiple structured populations with branch-specific effective population sizes (i.e., under the multispecies coalescent). Our solutions establish an expectation for genetic linkage in multispecies datasets and provide a new likelihood framework for linking demographic models with local ancestry inference across genomes.

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

Orobanchaceae is the largest family of parasitic plants, containing autotrophic and parasitic plants with all degrees of parasitism. This makes it by far the best family for studying the origin and evolution of plant parasitism. Here we provide three high-quality genomes of orobanchaceous plants, the autotrophic Lindenbergia luchunensis and the holoparasitic plants Phelipanche aegyptiaca and Orobanche cumana. Phylogenomic analysis of these three genomes together with those previously published and the transcriptomes of other orobanchaceous species, created a robust phylogenetic framework for Orobanchaceae. We found that an ancient whole-genome duplication (WGD; about 73.48 Mya), which occurred earlier than the origin of Orobanchaceae, might have contributed to the emergence of parasitism. However, no WGD events occurred in any lineage of orobanchaceous parasites except for Striga after divergence from their autotrophic common ancestor, suggesting that, in contrast to previous speculations, WGD is not associated with the emergence of holoparasitism. We detected evident convergent gene loss in all parasites within Orobanchaceae and between Orobanchaceae and dodder Cuscuta australis. The gene families in the orobanchaceous parasites showed a clear pattern of recent gains and expansions. The expanded gene families are enriched in functions related to the development of the haustorium, suggesting that recent gene family expansions may have facilitated the adaptation of orobanchaceous parasites to different hosts. This study illustrates a stepwise pattern in the evolution of parasitism in the orobanchaceous parasites, and will facilitate future studies on parasitism and the control of parasitic plants in agriculture.

The degree of gametophyte dependence on the sporophyte life stage is a major feature that differentiates the life cycles of land plants, yet the evolutionary consequences of this difference remain poorly understood. Most evolutionary models assume organisms are either haploid or diploid for their entire lifespan, which is not appropriate for simulating plant life cycles. Here, we introduceWe implemented evolutionary simulations under three realistic plant life cycle models supported inThe dynamics of single beneficial mutation fixation did not vary dramatically between different models, but the patterns of spatial variation did differ, demonstrating that different life histories and model parameters affect both genetic diversity and linkage disequilibrium. The rate of linkage disequilibrium decay away from selected sites varied depending on model parameters such as cloning and selfing rates, through their impact on effective population sizes.Evolutionary simulations are an exciting, underutilized approach in evolutionary research and education.

Replicated radiations, in which sets of similar forms evolve repeatedly within different regions, can provide powerful insights into parallel evolution and the assembly of functional diversity within communities. Several cases have been described in animals, but in plants we lack well-documented cases of replicated radiation that combine comprehensive phylogenetic and biogeographic analyses, the delimitation of geographic areas within which a set of 'ecomorphs' evolved independently and the identification of potential underlying mechanisms. Here we document the repeated evolution of a set of leaf ecomorphs in a group of neotropical plants. The Oreinotinus lineage within the angiosperm clade Viburnum spread from Mexico to Argentina through disjunct cloud forest environments. In 9 of 11 areas of endemism, species with similar sets of leaf forms evolved in parallel. We reject gene-flow-mediated evolution of similar leaves and show, instead, that species with disparate leaf forms differ in their climatic niches, supporting ecological adaptation as the driver of parallelism. Our identification of a case of replicated radiation in plants sets the stage for comparative analyses of such phenomena across the tree of life.

Two populations evolving in isolation can accumulate genetic differences over time that cause incompatibilities in their hybrid offspring. These “Dobzhansky-Muller incompatibilities” (DMIs) are predicted to accumulate at a rate faster than linear as the number of incompatible gene interactions “snowballs”. Here we show that this snowball prediction is an artifact of two unrealistic modeling assumptions that stem from ignoring demography. We introduce a new alternative “demographic speciation model” in which the rate of DMI accumulation between populations is affected by the efficiency of purifying selection to remove incompatibilities that arise within populations. This model yields new testable predictions for understanding the tempo and mode of speciation based on population demographic parameters. A large-scale empirical analysis of bird and mammal datasets supports a unique prediction of our model: a negative relationship between effective population sizes and speciation rates. Our results challenge views of the snowball theory, and of ecological speciation models rooted in positive selection, showing instead that purifying selection may play a major role in mediating speciation rates.

A key distinction between species tree inference under the multi-species coalescent model (MSC), and the inference of gene trees in sliding windows along a genome, is in the effect of genetic linkage. Whereas the MSC explicitly assumes genealogies to be unlinked, i.e., statistically independent, genealogies located close together on genomes are spatially auto-correlated. Here we use tree sequence simulations with recombination to explore the effects of species tree parameters on spatial patterns of linkage among genealogies. We decompose coalescent time units to demonstrate differential effects of generation time and effective population size on spatial coalescent patterns, and we define a new metric, “phylogenetic linkage,” for measuring the rate of decay of phylogenetic similarity by comparison to distances among unlinked genealogies. Finally, we provide a simple example where accounting for phylogenetic linkage in sliding window analyses improves local gene tree inference.

Background and aims Enlarged sterile flowers on the periphery of inflorescences increase the attractiveness of floral displays, and previous studies have generally demonstrated that these have positive effects on insect visitation and/or reproductive success. However, experiments have not specifically been designed to examine the benefits of sterile flowers under conditions that reflect the early stages in their evolution, i.e. when plants that produce sterile flowers are at low frequency. Methods Over three years, three experiments were performed in natural populations of Viburnum lantanoides, which produces sterile marginal flowers (SMFs). The first experiment established that fruit production in V. lantanoides increases with the receipt of outcross pollen. The second tested the role of SMFs under extant conditions, comparing fruit production in two populations composed entirely of intact plants or entirely of plants with the SMFs removed. The third was designed to mimic the presumed context in which SMFs first evolved; here, SMFs were removed from all but a few plants in a population, and rates of insect visitation and fruit set were compared between plants with intact and denuded SMFs. Key results In comparing whole populations, the presence of SMFs nearly doubled fruit set. Under simulated 'ancestral' conditions within a population, plants with intact SMFs received double the insect visits and produced significantly more fruits than denuded plants. There was no significant effect of the number of inflorescences or fertile flowers on insect visitation or fruit set, indicating that the presence of SMFs accounted for these differences. Conclusions The presence of SMFs significantly increased pollinator attraction and female reproductive success both in contemporary and simulated ancestral contexts, indicating that stabilizing selection is responsible for their maintenance, and directional selection likely drove their evolution when they first appeared. This study demonstrates a novel approach to incorporating historically relevant scenarios into experimental studies of floral evolution.

Major gaps remain in our understanding of the ecology, evolution, biodiversity, biogeography, extinction risk, and adaptive potential of reef building corals. One of the central challenges remains that there are few informative genetic markers for studying boundaries between species, and variation within species. Reduced representation sequencing approaches, such as RADseq (Restriction site Associated DNA sequencing) have great potential for resolving such relationships. However, it is necessary to identify loci in order to make inferences for endosymbiotic organisms such as corals. Here, we examined twenty-one coral holobiont ezRAD libraries from Hawai'i, focusing on P. lobata and P. compressa, two species with contrasting morphology and habitat preference that previous studies have not resolved. We used a combination of de novo assembly and reference mapping approaches to identify and compare loci: we used reference mapping to extract and compare nearly complete mitochondrial genomes, ribosomal arrays, and histone genes. We used de novo clustering and phylogenomic methods to compare the complete holobiont data set with coral and symbiont subsets that map to transcriptomic data. In addition, we used reference assemblies to examine genetic structure from SNPs (Single Nucleotide Polymorphisms). All approaches resolved outgroup taxa but failed to resolve P. lobata and P. compressa as distinct, with mito-nuclear discordance and shared mitochondrial haplotypes within the species complex. The holobiont and 'coral transcriptomic' datasets were highly concordant, revealing stronger genetic structure between sites than between coral morphospecies. These results suggest that either branching morphology is a polymorphic trait, or that these species frequently hybridize. This study provides examples of several approaches to acquire, identify, and compare loci across metagenomic samples such as the coral holobiont while providing insights into the nature of coral variability.

Made available in DSpace on 2020-12-12T02:47:34Z (GMT). No. of bitstreams: 0 Previous issue date: 2020-10-01 Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) The reconstruction of relationships within recently radiated groups is challenging even when massive amounts of sequencing data are available. The use of restriction site-associated DNA sequencing (RAD-Seq) to this end is promising. Here, we assessed the performance of RAD-Seq to infer the species-level phylogeny of the rapidly radiating genus Cereus (Cactaceae). To examine how the amount of genomic data affects resolution in this group, we used datasets and implemented different analyses. We sampled 52 individuals of Cereus, representing 18 of the 25 species currently recognized, plus members of the closely allied genera Cipocereus and Praecereus, and other 11 Cactaceae genera as outgroups. Three scenarios of permissiveness to missing data were carried out in iPyRAD, assembling datasets with 30% (333 loci), 45% (1440 loci), and 70% (6141 loci) of missing data. For each dataset, Maximum Likelihood (ML) trees were generated using two supermatrices, i.e., only SNPs and SNPs plus invariant sites. Accuracy and resolution were improved when the dataset with the highest number of loci was used (6141 loci), despite the high percentage of missing data included (70%). Coalescent trees estimated using SVDQuartets and ASTRAL are similar to those obtained by the ML reconstructions. Overall, we reconstruct a well-supported phylogeny of Cereus, which is resolved as monophyletic and composed of four main clades with high support in their internal relationships. Our findings also provide insights into the impact of missing data for phylogeny reconstruction using RAD loci. Departamento de Biologia Centro de Ciências Humanas e Biológicas Universidade Federal de São Carlos (UFSCar) Programa de pós-graduação em Biologia Comparada Faculdade de Filosofia Ciências e Letras de Ribeirão Preto Universidade de São Paulo (USP) Departamento de Genética e Biologia Evolutiva Universidade de São Paulo (USP) Department of Ecology Evolution and Environmental Biology Columbia University Instituto de Biotecnologia (IBTEC) e Instituto de Biociências de Botucatu (IBB) Universidade Estadual Paulista (UNESP) Singapore Botanic Gardens (National Parks Board), 1 Cluny Road Instituto Tecnológico Vale / Museu Paraense Emilio Goeldi Coord. Botânica Instituto de Biotecnologia (IBTEC) e Instituto de Biociências de Botucatu (IBB) Universidade Estadual Paulista (UNESP) FAPESP: 2014/25227-0 FAPESP: 2018/03428-5

Phylogeny, fossils, biogeography, and biome occupancy provide evidence that reflects the singular evolutionary history of a clade. Despite the connections that bind them together, these lines of evidence are most often studied separately, by first inferring a fossil-dated molecular phylogeny, then mapping on ancestral ranges and biomes inferred from extant species. Here we jointly model the evolution of biogeographic ranges, biome affinities, and molecular sequences, incorporating fossils to estimate a dated phylogeny for all of the 163 extant species of the woody plant cladeViburnum(Adoxaceae) that we currently recognize. Our analyses indicate that while the majorViburnumlineages evolved in the Eocene, the majority of extant species originated since the Miocene.Viburnumradiated first in Asia, in warm, broad-leaved evergreen (lucidophyllous) forests. Within Asia we infer several early shifts into more tropical forests, and multiple shifts into forests that experience prolonged freezing. From Asia we infer two early movements into the New World. These two lineages probably first occupied warm temperate forests and adapted later to spreading cold climates. One of these lineages (Porphyrotinus) occupied cloud forests and moved south through the mountains of the Neotropics. Several other movements into North America took place more recently, facilitated by prior adaptations to freezing in the Old World. We also infer four disjunctions between Asia and Europe: theTinuslineage is the oldest and probably occupied warm forests when it spread, while the other three were more recent and in cold-adapted lineages. These results variously contradict published accounts, especially the view thatViburnumradiated initially in cold forests and, accordingly, maintained vessel elements with scalariform perforations. We explored how the location and biome assignments of fossils affected our inference of ancestral areas and biome states. Our results are sensitive to, but not entirely dependent upon, the inclusion of fossil biome data. We argue that it will be critical to take advantage of all available lines of evidence to decipher events in the distant past, and the joint estimation approach developed here provides cautious hope even when fossil evidence is limited.

The Hengduan Mountains of southwest China is a temperate biodiversity hotspot characterized by high rates of endemism, including approximately 300 endemic species of Pedicularis. Isolation in allopatry has likely contributed to high rates of plant speciation in this region. However, infrastructure development in recent decades has substantially altered connectivity of this landscape -- connecting previously isolated areas by roads and tunnels -- with the potential to threaten endemic biodiversity with genetic homogenization. We are using comparative phylogeographic methods to reconstruct the history of divergence and gene flow among multiple widespread species of Pedicularis to identify major barriers that have contributed to rapid speciation processes. Here, we present preliminary results for one widespread species complex that includes P. cranolopha and P. tricolor. Using RAD-seq data we estimated a rooted phylogeny, population structure, and effective migration surfaces to characterize the history of divergence and migration in this clade.

1AbstractThe fig and pollinator wasp obligate mutualism is diverse (~750 described species), ecologically important, and ancient (~80-90 Ma), providing model systems for generating and testing many questions in evolution and ecology. Once thought to be a prime example of strict one-to-one cospeciation, current thinking suggests that genera of pollinator wasps co-evolve with corresponding subsections of figs, but the degree to which cospeciation or other processes contributes to the association at finer scales is unclear. Here we use genome-wide sequence data from a community of Panamanian strangler figs (Ficus subgenus Urostigma, section Americana) and associated fig wasp pollinators (Pegoscapus spp.) to infer the process of coevolutionary diversification in this obligate mutualism. Using a model-based approach adapted from the study of gene family evolution, our results indicate pervasive and ongoing host switching of pollinator wasps at this fine phylogenetic and regional scale. Although the model estimates a modest amount of cospeciation, simulations reveal this signal to be consistent with levels of co-association expected under a model of free host switching. Our findings provide an outline for testing how ecological and evolutionary processes can be modeled to evaluate the history of association of interacting lineages in a phylogenetic framework.

The past decade has seen a major breakthrough in our ability to easily and inexpensively sequence genome-scale data from diverse lineages. The development of high-throughput sequencing and long-read technologies has ushered in the era of phylogenomics, where hundreds to thousands of nuclear genes and whole organellar genomes are routinely used to reconstruct evolutionary relationships. As a result, understanding which options are best suited for a particular set of questions can be difficult, especially for those just starting in the field. Here, we review the most recent advances in plant phylogenomic methods and make recommendations for project-dependent best practices and considerations. We focus on the costs and benefits of different approaches in regard to the information they provide researchers and the questions they can address. We also highlight unique challenges and opportunities in plant systems, such as polyploidy, reticulate evolution, and the use of herbarium materials, identifying optimal methodologies for each. Finally, we draw attention to lingering challenges in the field of plant phylogenomics, such as reusability of data sets, and look at some up-and-coming technologies that may help propel the field even further.

Introgressive hybridization challenges the concepts we use to define species and infer phylogenetic relationships. Methods for inferring historical introgression from the genomes of extant species, such as ABBA-BABA tests, are widely used, however, their results can be easily misinterpreted. Because these tests are inherently comparative, they are sensitive to the effects of missing data (unsampled species) and nonindependence (hierarchical relationships among species). We demonstrate this using genomic RADseq data sampled from all extant species in the American live oaks (Quercus series Virentes), a group notorious for hybridization. By considering all species and their phylogenetic relationships, we were able to distinguish true hybridizing lineages from those that falsely appear admixed. Six of seven species show evidence of admixture, often with multiple other species, but which is explained by introgression among a few related lineages occurring in close proximity. We identify the Cuban oak as the most admixed lineage and test alternative scenarios for its origin. The live oaks form a continuous ring-like distribution around the Gulf of Mexico, connected in Cuba, across which they could effectively exchange alleles. However, introgression appears highly localized, suggesting that oak species boundaries and their geographic ranges have remained relatively stable over evolutionary time.

Genome-wide assessments allow for fuller characterization of genetic diversity, finer-scale population delineation, and better detection of demographically significant units to guide conservation compared with those based on "traditional" markers. Galapagos giant tortoises (Chelonoidis spp.) have long provided a case study for how evolutionary genetics may be applied to advance species conservation. Ongoing efforts to bolster tortoise populations, which have declined by 90%, have been informed by analyses of mitochondrial DNA sequence and microsatellite genotypic data, but could benefit from genome-wide markers. Taking this next step, we used double-digest restriction-site associated DNA sequencing to collect genotypic data at >26000 single nucleotide polymorphisms (SNPs) for 117 individuals representing all recognized extant Galapagos giant tortoise species. We then quantified genetic diversity, population structure, and compared results to estimates from mitochondrial DNA and microsatellite loci. Our analyses detected 12 genetic lineages concordant with the 11 named species as well as previously described structure within one species, C. becki. Furthermore, the SNPs provided increased resolution, detecting admixture in 4 individuals. SNP-based estimates of diversity and differentiation were significantly correlated with those derived from nuclear microsatellite loci and mitochondrial DNA sequences. The SNP toolkit presented here will serve as a resource for advancing efforts to understand tortoise evolution, species radiations, and aid conservation of the Galapagos tortoise species complex.

Species are the starting point for most studies of ecology and evolution, but the proper circumscription of species can be extremely difficult in morphologically variable lineages, and there are still few convincing examples of molecularly informed species delimitation in plants. Here, we focus on the Viburnum nudum complex, a highly variable clade that is widely distributed in eastern North America. Taxonomic treatments have mostly divided this complex into northern (V. nudum var. cassinoides) and southern (V. nudum var. nudum) entities, but additional names have been proposed. We used multiple lines of evidence, including RADseq, morphological, and geographic data, to test how many independently evolving lineages exist within the V. nudum complex. Genetic clustering and phylogenetic methods revealed three distinct groups-one lineage that is highly divergent, and two others that are recently diverged and morphologically similar. A combination of evidence that includes reciprocal monophyly, lack of introgression, and discrete rather than continuous patterns of variation supports the recognition of all three lineages as separate species. These results identify a surprising case of cryptic diversity in which two broadly sympatric species have consistently been lumped in taxonomic treatments. The clarity of our findings is directly related to the dense sampling and high-quality genetic data in this study. We argue that there is a critical need for carefully sampled and integrative species delimitation studies to clarify species boundaries even in well-known plant lineages. Studies following the model that we have developed here are likely to identify many more cryptic lineages and will fundamentally improve our understanding of plant speciation and patterns of species richness.

Phylogenetic relationships among temperate species of bamboo are difficult to resolve, owing to both the challenge of detecting sufficiently variable markers and their polyploid history. Here, we use restriction site-associated DNA sequencing to identify candidate loci with fixed allelic differences segregating between and within two temperate species of bamboos: Arundinaria faberi and Yushania brevipaniculata. Approximately 27 million paired-end sequencing reads were generated across four samples. From pooled data, we assembled 67685 and 70668 de novo contigs from partial overlap among paired-end reads, with an average length of 240 and 241bp for the two species, respectively, which were used to investigate functional classification of RAD tags in a blastx search. Analysed separately by population, we recovered 29443 putatively orthologous RAD tags shared across the four sampled populations, containing 28023 sequence variants, of which c. 13000 are segregating between species, and c. 3000 segregating between populations within each species. Analyses based on these RAD tags yielded robust phylogenetic inferences, even with data set constructed from surprisingly few loci. This study illustrates the potential for reduced-representation genome data to resolve difficult phylogenetic relationships in temperate bamboos.

Phylogenetic relationships among recently diverged species are often difficult to resolve due to insufficient phylogenetic signal in available markers and/or conflict among gene trees. Here we explore the use of reduced-representation genome sequencing, specifically in the form of restriction-site associated DNA (RAD), for phylogenetic inference and the detection of ancestral hybridization in non-model organisms. As a case study, we investigate Pedicularis section Cyathophora, a systematically recalcitrant clade of flowering plants in the broomrape family (Orobanchaceae). Two methods of phylogenetic inference, maximum likelihood and Bayesian concordance, were applied to data sets that included as many as 40,000 RAD loci. Both methods yielded similar topologies that included two major clades: a “rex-thamnophila” clade, composed of two species and several subspecies with relatively low floral diversity, and geographically widespread distributions at lower elevations, and a “superba” clade, composed of three species characterized by relatively high floral diversity and isolated geographic distributions at higher elevations. Levels of molecular divergence between subspecies in the rex-thamnophila clade are similar to those between species in the superba clade. Using Patterson’s D-statistic test, including a novel extension of the method that enables finer-grained resolution of introgression among multiple candidate taxa by removing the effect of their shared ancestry, we detect significant introgression among nearly all taxa in the rex-thamnophila clade, but not between clades or among taxa within the superba clade. These results suggest an important role for geographic isolation in the emergence of species barriers, by facilitating local adaptation and differentiation in the absence of homogenizing gene flow. [Concordance factors; genotyping-by-sequencing; hybridization; partitioned D-statistic test; Pedicularis; restriction-site associated DNA.]

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment-specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2-5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.

A pervasive hypothesis at the interface of ecology and evolution is that biotic interactions contribute to regional biodiversity by accelerating adaptation and speciation. We investigated this question in the context of closely related, bumble bee-pollinated plants (Pedicularis spp.) in the Hengduan Mountains of south-central China, where they exhibit spectacular levels of richness, endemism, and floral diversity. Because these species co-occur frequently, flower synchronously, and share pollinators during the brief reproductive season, we predict that pollinator-mediated interactions may influence their community assembly and evolutionary diversification. If disparity in floral traits reduces competitive interactions between species, as would happen if floral isolation mitigates reproductive interference caused by heterospecific pollen flow, then species with dissimilar flowers should co-occur more often, yielding greater floral diversity at local scales than expected by chance. Moreover, if such interactions have repeatedly driven character displacement, then floral traits should exhibit homoplasy, the phylogenetic signature of labile evolution. We present evidence supporting these predictions, and find that local species richness is best explained by a model including both floral diversity and phylogenetic distance. Our results suggest that a dynamic mosaic of pollinator-mediated interactions among Pedicularis in the Hengduan region promotes ecological sorting through recurrent selection against reproductive interference, causing rapid species turnover at local scales, and accelerating the rate of floral divergence among species. Together these processes may have contributed to the remarkable accumulation of florally diverse species of Pedicularis endemic to the Hengduan Mountains biodiversity hotspot.

Introgressive hybridization challenges the concepts we use to define species and our ability to infer their evolutionary relationships. Methods for inferring historical introgression from the genomes of extant species are now widely used, however, few guidelines have been articulated for how best to interpret their results. Because these tests are inherently comparative, we show that they are sensitivite to the effects of missing data (unsampled species) and to non-independence (hierarchical relationships among species). We demonstrate this using genomic RAD data sampled from populations across the geographic ranges of all extant species in the American live oaks (QuercusseriesVirentes), a group notorious for hybridization. By considering all species in the clade, and their phylogenetic relationships, we were able to distinguish true hybridizing lineages from those that falsely appear admixed due to phylogenetic structure among hybridizing relatives. Six of seven species show evidence of admixture, often with multiple other species, but which can be explained by hybrid introgression among few related lineages where they occur in close proximity. We identify the Cuban oak as a highly admixed lineage and use an information-theoretic model comparison approach to test alternative scenarios for its origin. Hybrid speciation is a poor fit compared to a model in which a population from Central America colonized Cuba and received subsequent gene flow from Florida. The live oaks form a continuous ring-like distribution around the Gulf of Mexico, connected in Cuba, across which they could effectively exchange alleles. However, introgression appears to remain localized to areas of sympatry, suggesting that oak species boundaries, and their geographic ranges have remained relatively stable over evolutionary time.

The nature and timing of evolution of niche differentiation among closely related species remains an important question in ecology and evolution. The American live oak clade, Virentes, which spans the unglaciated temperate and tropical regions of North America and Mesoamerica, provides an instructive system in which to examine speciation and niche evolution. We generated a fossil-calibrated phylogeny of Virentes using RADseq data to estimate divergence times and used nuclear microsatellites, chloroplast sequences and an intron region of nitrate reductase (NIA-i3) to examine genetic diversity within species, rates of gene flow among species and ancestral population size of disjunct sister species. Transitions in functional and morphological traits associated with ecological and climatic niche axes were examined across the phylogeny. We found the Virentes to be monophyletic with three subclades, including a southwest clade, a southeastern US clade and a Central American/Cuban clade. Despite high leaf morphological variation within species and transpecific chloroplast haplotypes, RADseq and nuclear SSR data showed genetic coherence of species. We estimated a crown date for Virentes of 11 Ma and implicated the formation of the Sea of Cortes in a speciation event ~5 Ma. Tree height at maturity, associated with fire tolerance, differs among the sympatric species, while freezing tolerance appears to have diverged repeatedly across the tropical–temperate divide. Sympatric species thus show evidence of ecological niche differentiation but share climatic niches, while allopatric and parapatric species conserve ecological niches, but diverge in climatic niches. The mode of speciation and/or degree of co-occurrence may thus influence which niche axis plants diverge along.

Determining phylogenetic relationships among very closely related species has remained a challenge for evolutionary biologists due to interlocus phylogenetic discordance and the difficulty of obtaining variable markers. Here, we used a Genotyping-by-Sequencing (GBS) approach to sample a reduced representation genomic data set and infer the phylogeny of seven closely related species in the genus Carex (Cyperaceae). Past attempts to reconstruct phylogenetic relationships among these species produced conflicting and poorly-supported results. We inferred a robust phylogeny based on >3000 GBS loci and >1300 SNPs (with a minimum sequence depth within individuals of 10) using maximum likelihood and Bayesian inference. We also tested for historical introgression using the D-statistic test. We compared these analyses with partitioned RAD analysis, which is designed to identify suboptimal trees reflecting secondary phylogenetic signal that may be obscured by the dominant signal in the data. Phylogenetic analyses yielded fully resolved trees with high support. We found two main clades, one grouping Carex scoparia populations and C. waponahkikensis, and a second clade grouping C. longii, C. vexans, C. suberecta and C. albolutescens. We detected marginally significant signals of introgression between C. scoparia and C. suberecta or C. albolutescens, and we rejected a hybrid origin hypothesis for C. waponahkikensis. Our results demonstrate the power of NGS data sets for resolving some of the most difficult phylogenetic challenges where traditional phylogenetic markers have failed.

Restriction-site associated genomic markers are a powerful tool for investigating evolutionary questions at the population level, but are limited in their utility at deeper phylogenetic scales where fewer orthologous loci are typically recovered across disparate taxa. While this limitation stems in part from mutations to restriction recognition sites that disrupt data generation, an alternative source of data loss comes from the failure to identify homology during bioinformatic analyses. Clustering methods that allow for lower similarity thresholds and the inclusion of indel variation will perform better at assembling RADseq loci at the phylogenetic scale.PyRADis a pipeline to assemblede novoRADseq loci with the aim of optimizing coverage across phylogenetic data sets. It utilizes a wrapper around an alignment-clustering algorithm which allows for indel variation within and between samples, as well as for incomplete overlap among reads (e.g., paired-end). Here I comparePyRADwith the programStacksin their performance analyzing a simulated RADseq data set that includes indel variation. Indels disrupt clustering of homologous loci inStacksbut not inPyRAD, such that the latter recovers more shared loci across disparate taxa. I show through re-analysis of an empirical RADseq data set that indels are a common feature of such data, even at shallow phylogenetic scales.PyRADutilizes parallel processing as well as an optional hierarchical clustering method which allow it to rapidly assemble phylogenetic data sets with hundreds of sampled individuals.

Restriction-site associated DNA (RAD) sequencing and related methods rely on the conservation of enzyme recognition sites to isolate homologous DNA fragments for sequencing, with the consequence that mutations disrupting these sites lead to missing information. There is thus a clear expectation for how missing data should be distributed, with fewer loci recovered between more distantly related samples. This observation has led to a related expectation: that RAD-seq data are insufficiently informative for resolving deeper scale phylogenetic relationships. Here we investigate the relationship between missing information among samples at the tips of a tree and information at edges within it. We re-analyze and review the distribution of missing data across ten RAD-seq data sets and carry out simulations to determine expected patterns of missing information. We also present new empirical results for the angiosperm clade Viburnum (Adoxaceae, with a crown age >50 Ma) for which we examine phylogenetic information at different depths in the tree and with varied sequencing effort. The total number of loci, the proportion that are shared, and phylogenetic informativeness varied dramatically across the examined RAD-seq data sets. Insufficient or uneven sequencing coverage accounted for similar proportions of missing data as dropout from mutation-disruption. Simulations reveal that mutation-disruption, which results in phylogenetically distributed missing data, can be distinguished from the more stochastic patterns of missing data caused by low sequencing coverage. In Viburnum, doubling sequencing coverage nearly doubled the number of parsimony informative sites, and increased by >10X the number of loci with data shared across >40 taxa. Our analysis leads to a set of practical recommendations for maximizing phylogenetic information in RAD-seq studies. [hierarchical redundancy; phylogenetic informativeness; quartet informativeness; Restriction-site associated DNA (RAD) sequencing; sequencing coverage; Viburnum.].

Previous phylogenetic studies in oaks (Quercus, Fagaceae) have failed to resolve the backbone topology of the genus with strong support. Here, we utilize next-generation sequencing of restriction-site associated DNA (RAD-Seq) to resolve a framework phylogeny of a predominantly American clade of oaks whose crown age is estimated at 23-33 million years old. Using a recently developed analytical pipeline for RAD-Seq phylogenetics, we created a concatenated matrix of 1.40 E06 aligned nucleotides, constituting 27,727 sequence clusters. RAD-Seq data were readily combined across runs, with no difference in phylogenetic placement between technical replicates, which overlapped by only 43-64% in locus coverage. 17% (4,715) of the loci we analyzed could be mapped with high confidence to one or more expressed sequence tags in NCBI Genbank. A concatenated matrix of the loci that BLAST to at least one EST sequence provides approximately half as many variable or parsimony-informative characters as equal-sized datasets from the non-EST loci. The EST-associated matrix is more complete (fewer missing loci) and has slightly lower homoplasy than non-EST subsampled matrices of the same size, but there is no difference in phylogenetic support or relative attribution of base substitutions to internal versus terminal branches of the phylogeny. We introduce a partitioned RAD visualization method (implemented in the R package RADami; http://cran.r-project.org/web/packages/RADami) to investigate the possibility that suboptimal topologies supported by large numbers of loci--due, for example, to reticulate evolution or lineage sorting--are masked by the globally optimal tree. We find no evidence for strongly-supported alternative topologies in our study, suggesting that the phylogeny we recover is a robust estimate of large-scale phylogenetic patterns in the American oak clade. Our study is one of the first to demonstrate the utility of RAD-Seq data for inferring phylogeny in a 23-33 million year-old clade.