Legacy of supervolcanic eruptions on population genetic structure of brown kiwi.

Supervolcanoes are volcanoes capable of mega-colossal eruptions that emit more than 1,000 km3 of ash and other particles.1 The earth's most recent mega-colossal eruption was the Oruanui eruption of the Taupo supervolcano 25,580 years before present (YBP) on the central North Island of New Zealand.2 This eruption blanketed major swaths of the North Island in thick layers of ash and igneous rock,2,3 devastating habitats and likely causing widespread population extinctions.4-7 An additional devastating super-colossal eruption (>100 km3) of the Taupo supervolcano occurred approximately 1,690 YBP.8 The impacts of such massive but ephemeral natural disasters on contemporary population genetic structure remain underexplored. Here, we combined data for 4,951 SNPs with spatially explicit demographic and coalescent models within an approximate Bayesian computation framework to test the drivers of genetic structure in brown kiwi (Apteryx mantelli). Our results strongly support the importance of eruptions of the Taupo supervolcano in restructuring pre-existing geographic patterns of population differentiation and genetic diversity. Range shifts due to climatic oscillations-a frequent explanation for genetic structure9-are insufficient to fully explain the empirical data. Meanwhile, recent range contraction and fragmentation due to historically documented anthropogenic habitat alteration adds no explanatory power to our models. Our results support a major role for cycles of destruction and post-volcanic recolonization in restructuring the population genomic landscape of brown kiwi and highlight how ancient and ephemeral mega-disasters may leave a lasting legacy on patterns of intraspecific genetic variation.

Demographic decline and lineage-specific adaptations characterize New Zealand kiwi.

Small and fragmented populations may become rapidly differentiated due to genetic drift, making it difficult to distinguish whether neutral genetic structure is a signature of recent demographic events, or of long-term evolutionary processes that could have allowed populations to adaptively diverge. We sequenced 52 whole genomes to examine Holocene demographic history and patterns of adaptation in kiwi (Apteryx), and recovered 11 strongly differentiated genetic clusters corresponding to previously recognized lineages. Demographic models suggest that all 11 lineages experienced dramatic population crashes relative to early- or mid-Holocene levels. Small population size is associated with low genetic diversity and elevated genetic differentiation (FST), suggesting that population declines have strengthened genetic structure and led to the loss of genetic diversity. However, population size is not correlated with inbreeding rates. Eight lineages show signatures of lineage-specific selective sweeps (284 sweeps total) that are unlikely to have been caused by demographic stochasticity. Overall, these results suggest that despite strong genetic drift associated with recent bottlenecks, most kiwi lineages possess unique adaptations and should be recognized as separate adaptive units in conservation contexts. Our work highlights how whole-genome datasets can address longstanding uncertainty about the evolutionary and conservation significance of small and fragmented populations of threatened species.

Geographic contact drives increased reproductive isolation in two cryptic Empidonax flycatchers.

Geographic contact between sister lineages often occurs near the final stages of speciation, but its role in speciation's completion remains debated. Reproductive isolation may be essentially complete prior to secondary contact. Alternatively, costly interactions between partially reproductively isolated species - such as maladaptive hybridization or competition for resources - may select for divergence, increasing reproductive isolation and driving speciation toward completion. Here, we use coalescent demographic modelling and whole-genome data sets to show that a period of contact and elevated hybridization between sympatric eastern North American populations of two cryptic bird species preceded a major increase in reproductive isolation between these populations within the last 10,000 years. In contrast, substantial introgression continues to the present in a western contact zone where geographic overlap is much narrower and probably of more recent origin. In the sympatric eastern region where reproductive isolation has increased, it is not accompanied by character displacement in key morphometric traits, plumage coloration, or ecological traits. While the precise trait and underlying mechanism driving increased reproductive isolation remains unknown, we discuss several possibilities and outline avenues for future research. Overall, our results highlight how demographic models can reveal the geographic context in which reproductive isolation was completed, and demonstrate how contact can accelerate the final stages of speciation.

Small spaces, big impacts: contributions of micro-environmental variation to population persistence under climate change.

Individuals within natural populations can experience very different abiotic and biotic conditions across small spatial scales owing to microtopography and other micro-environmental gradients. Ecological and evolutionary studies often ignore the effects of micro-environment on plant population and community dynamics. Here, we explore the extent to which fine-grained variation in abiotic and biotic conditions contributes to within-population variation in trait expression and genetic diversity in natural plant populations. Furthermore, we consider whether benign microhabitats could buffer local populations of some plant species from abiotic stresses imposed by rapid anthropogenic climate change. If microrefugia sustain local populations and communities in the short term, other eco-evolutionary processes, such as gene flow and adaptation, could enhance population stability in the longer term. We caution, however, that local populations may still decline in size as they contract into rare microhabitats and microrefugia. We encourage future research that explicitly examines the role of the micro-environment in maintaining genetic variation within local populations, favouring the evolution of phenotypic plasticity at local scales and enhancing population persistence under global change.

Climate change shifts natural selection and the adaptive potential of the perennial forb Boechera stricta in the Rocky Mountains.

Heritable genetic variation is necessary for populations to evolve in response to anthropogenic climate change. However, antagonistic genetic correlations among traits may constrain the rate of adaptation, even if substantial genetic variation exists. We examine potential genetic responses to selection by comparing multivariate genetic variance-covariances of traits and fitness (multivariate Robertson-Price identities) across different environments in a reciprocal transplant experiment of the forb Boechera stricta in the Rocky Mountains. By transplanting populations into four common gardens arrayed along an elevational gradient, and exposing populations to control and snow removal treatments, we simulated future and current climates and snowmelt regimes. Genetic variation in flowering and germination phenology declined in plants moved downslope to warmer, drier sites, suggesting that these traits may have a limited ability to evolve under future climates. Simulated climate change via snow removal altered the strength of selection on flowering traits, but we found little evidence that genetic correlations among traits are likely to affect the rate of adaptation to climate change. Overall, our results suggest that climate change may alter the evolutionary potential of B. stricta, but reduced expression of genetic variation may be a larger impediment to adaptation than constraints imposed by antagonistic genetic correlations.

Genomic evidence of survival near ice sheet margins for some, but not all, North American trees.

Temperate species experienced dramatic range reductions during the Last Glacial Maximum, yet refugial populations from which modern populations are descended have never been precisely located. Climate-based models identify only broad areas of potential habitat, traditional phylogeographic studies provide poor spatial resolution, and pollen records for temperate forest communities are difficult to interpret and do not provide species-level taxonomic resolution. Here we harness signals of range expansion from large genomic datasets, using a simulation-based framework to infer the precise latitude and longitude of glacial refugia in two widespread, codistributed hickories (Carya spp.) and to quantify uncertainty in these estimates. We show that one species likely expanded from close to ice sheet margins near the site of a previously described macrofossil for the genus, highlighting support for the controversial notion of northern microrefugia. In contrast, the expansion origin inferred for the second species is compatible with classic hypotheses of distant displacement into southern refugia. Our statistically rigorous, powerful approach demonstrates how refugia can be located from genomic data with high precision and accuracy, addressing fundamental questions about long-term responses to changing climates and providing statistical insight into longstanding questions that have previously been addressed primarily qualitatively.

Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species.

High-throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double-digest restriction-associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra- and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of FST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high-throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.

Tests of species-specific models reveal the importance of drought in postglacial range shifts of a Mediterranean-climate tree: insights from integrative distributional, demographic and coalescent modelling and ABC model selection.

Past climate change has caused shifts in species distributions and undoubtedly impacted patterns of genetic variation, but the biological processes mediating responses to climate change, and their genetic signatures, are often poorly understood. We test six species-specific biologically informed hypotheses about such processes in canyon live oak (Quercus chrysolepis) from the California Floristic Province. These hypotheses encompass the potential roles of climatic niche, niche multidimensionality, physiological trade-offs in functional traits, and local-scale factors (microsites and local adaptation within ecoregions) in structuring genetic variation. Specifically, we use ecological niche models (ENMs) to construct temporally dynamic landscapes where the processes invoked by each hypothesis are reflected by differences in local habitat suitabilities. These landscapes are used to simulate expected patterns of genetic variation under each model and evaluate the fit of empirical data from 13 microsatellite loci genotyped in 226 individuals from across the species range. Using approximate Bayesian computation (ABC), we obtain very strong support for two statistically indistinguishable models: a trade-off model in which growth rate and drought tolerance drive habitat suitability and genetic structure, and a model based on the climatic niche estimated from a generic ENM, in which the variables found to make the most important contribution to the ENM have strong conceptual links to drought stress. The two most probable models for explaining the patterns of genetic variation thus share a common component, highlighting the potential importance of seasonal drought in driving historical range shifts in a temperate tree from a Mediterranean climate where summer drought is common.

Time to get moving: assisted gene flow of forest trees.

Geographic variation in trees has been investigated since the mid-18th century. Similar patterns of clinal variation have been observed along latitudinal and elevational gradients in common garden experiments for many temperate and boreal species. These studies convinced forest managers that a 'local is best' seed source policy was usually safest for reforestation. In recent decades, experimental design, phenotyping methods, climatic data and statistical analyses have improved greatly and refined but not radically changed knowledge of clines. The maintenance of local adaptation despite high gene flow suggests selection for local adaptation to climate is strong. Concerns over maladaptation resulting from climate change have motivated many new genecological and population genomics studies; however, few jurisdictions have implemented assisted gene flow (AGF), the translocation of pre-adapted individuals to facilitate adaptation of planted forests to climate change. Here, we provide evidence that temperate tree species show clines along climatic gradients sufficiently similar for average patterns or climate models to guide AGF in the absence of species-specific knowledge. Composite provenancing of multiple seed sources can be used to increase diversity and buffer against future climate uncertainty. New knowledge will continue to refine and improve AGF as climates warm further.

The effects of read length, quality and quantity on microsatellite discovery and primer development: from Illumina to PacBio.

The advent of next-generation sequencing (NGS) technologies has transformed the way microsatellites are isolated for ecological and evolutionary investigations. Recent attempts to employ NGS for microsatellite discovery have used the 454, Illumina, and Ion Torrent platforms, but other methods including single-molecule real-time DNA sequencing (Pacific Biosciences or PacBio) remain viable alternatives. We outline a workflow from sequence quality control to microsatellite marker validation in three plant species using PacBio circular consensus sequencing (CCS). We then evaluate the performance of PacBio CCS in comparison with other NGS platforms for microsatellite isolation, through simulations that focus on variations in read length, read quantity and sequencing error rate. Although quality control of CCS reads reduced microsatellite yield by around 50%, hundreds of microsatellite loci that are expected to have improved conversion efficiency to functional markers were retrieved for each species. The simulations quantitatively validate the advantages of long reads and emphasize the detrimental effects of sequencing errors on NGS-enabled microsatellite development. In view of the continuing improvement in read length on NGS platforms, sequence quality and the corresponding strategies of quality control will become the primary factors to consider for effective microsatellite isolation. Among current options, PacBio CCS may be optimal for rapid, small-scale microsatellite development due to its flexibility in scaling sequencing effort, while platforms such as Illumina MiSeq will provide cost-efficient solutions for multispecies microsatellite projects.

Low genetic diversity, moderate local adaptation, and phylogeographic insights in Cornus nuttallii (Cornaceae).

PREMISE OF THE STUDY: Genetic knowledge is completely lacking for Pacific dogwood (Cornus nuttallii), a western North American temperate tree that is pollinated and dispersed by biological vectors. We investigated how history, geography, and climate have affected population genetic structure, local adaptation, and the phylogeography of this species. METHODS: We examined patterns and levels of diversity in nuclear microsatellites (SSRs) and cpDNA haplotypes in populations from across the species range. We compared these results to population differentiation and genetic clines in phenotypic traits in a common garden. KEY RESULTS: Genetic diversity was low for both nuclear SSRs and cpDNA. There was a lack of population structure (F(ST) = 0.090) in the coastal portion of the species range, with estimates of population genetic diversity in microsatellite markers decreasing with latitude from California to British Columbia. A disjunct interior population in Idaho 450 km from the coastal range had the lowest diversity but the highest divergence of all populations studied. Only a single nucleotide polymorphism was discovered after sequencing 5547 base pairs in seven noncoding regions of cpDNA. Both cpDNA haplotypes were widely distributed throughout the species range. Quantitative variation among populations was moderate (0.11 ≤ Q(ST) ≤ 0.63), and weak but significant adaptive clines were found between quantitative traits and population climatic variables (0.09 ≤ R(2) ≤ 0.34). CONCLUSIONS: Cornus nuttallii likely faced a population bottleneck in a single southern refugium during the Last Glacial Maximum. Despite low genetic diversity, it is weakly to moderately locally adapted.