Arctic-alpine distributions--metapopulations on a continental scale?

Relevance of ddRADseq method for species and population delimitation of closely related and widely distributed wolf spiders (Araneae, Lycosidae)

Although species delimitation is often controversial, emerging DNA-based and classical morphology-based methods are rarely compared using large-scale samplings, even less in the case of widely distributed species that have distant, allopatric populations. In the current study, we examined species boundaries within two wolf spider species of the genus Pardosa (Araneae, Lycosidae), P. riparia and P. palustris. Wolf spiders constitute an excellent model for testing the relevance of traditional vs. modern methods in species and population delimitation because several closely related species are distributed over cross-continental geographic ranges. Allopatric populations of the two Pardosa species were sampled across Europe to Far East Russia (latitudinal range > 150°) and several dozen individuals were studied using morphological characters (morphometry of three measures for both sexes, plus five in males only and two in females only), DNA barcoding (COI sequencing) and double-digest restriction site associated DNA sequencing (ddRADseq). The results obtained allow for changing the taxonomic status of two Far East Russian populations to subspecies and ddRADseq proved to be a powerful tool for taxonomic research despite scarce sampling and inherent subjectivity of species delimitation in allopatry. Overall, this study pleads for both multi-criteria and more population-based studies in taxonomy.

Relicts from Glacial Times: The Ground Beetle Pterostichus adstrictus Eschscholtz, 1823 (Coleoptera: Carabidae) in the Austrian Alps

Simple Summary

The extant distribution of many plants and animals is the result of the dynamics of the last ice ages with their recurrent advances and retreats of the northern ice sheet and the glaciers in the mountains. The arctic-alpine distribution is a special case where a species occurs in the subarctic/arctic regions and locally restricted in the alpine mountain regions of central or southeastern Europe. Among the ground beetles, several species display this type of distribution, one of which is Pterostichus adstrictus. In Europe, this ground beetle has been thought to have its southernmost occurrences in Wales and southern Scandinavia. In this study, we provide the first reliable record of P. adstrictus from the Austrian Alps based on morphological determination and comparison to other closely related species as well as molecular genetic data. Furthermore, the seasonal occurrence as well as empirical habitat preferences of P. adstrictus in the Austrian Alps are described.

Abstract

The last ice age considerably influenced distribution patterns of extant species of plants and animals, with some of them now inhabiting disjunct areas in the subarctic/arctic and alpine regions. This arctic-alpine distribution is characteristic for many cold-adapted species with a limited dispersal ability and can be found in many invertebrate taxa, including ground beetles. The ground beetle Pterostichus adstrictus Eschscholtz, 1823 of the subgenus Bothriopterus was previously known to have a holarctic-circumpolar distribution, in Europe reaching its southern borders in Wales and southern Scandinavia. Here, we report the first findings of this species from the Austrian Ötztal Alps, representing also the southernmost edge of its currently known distribution, confirmed by the comparison of morphological characters to other Bothriopterus species and DNA barcoding data. Molecular data revealed a separation of the Austrian and Finish specimens with limited to no gene flow at all. Furthermore, we present the first data on habitat preference and seasonality of P. adstrictus in the Austrian Alps.

Inter and intra island introgression in a wolf spider radiation from the Galápagos, and its implications for parallel evolution

Parallel radiations within island systems are often assumed to follow a simple scenario in which single colonization events are followed by in situ adaptive divergence. However, subsequent gene exchange after the initial colonization and during the divergence process might have important evolutionary impacts on species radiations. Gene exchange among ecologically similar species from different islands may lead to introgression of adaptive genetic variation and influence the parallel divergence process. In this study, we estimate levels of gene exchange within a wolf spider radiation of the genus Hogna Simon, 1885, from the Galápagos, wherein habitat specialization into 'high elevation' and 'coastal dry' species apparently evolved repeatedly on two islands. By using a multilocus approach we show that low levels of inter-island and relatively higher levels of intra island introgression shaped genetic variation in this species complex. Using these estimates, we demonstrate by means of a coalescence simulation that under these inter- and intra-island migration rates parallel evolution most likely evolves by introgression of adaptive alleles among islands, rather than through independent mutations despite the close genetic relationship of species within islands. As species phylogenies within radiations are frequently used to infer the divergence pattern, even relatively low levels of interspecific gene flow should not be neglected when interpreting parallel trait evolution.

Testing mitochondrial marker efficacy for DNA barcoding in spiders: a test case using the dwarf spider genus Oedothorax (Araneae : Linyphiidae : Erigoninae)

The present study focusses on comparatively assessing the efficacy for DNA barcoding of the two most commonly used mitochondrial markers (cox1 and 16S) in a genus of erigonine spiders. In total, 53 specimens representing five species, including four multi-sampled species, were sampled from several European localities. Initial evaluation of species monophyly was performed through parsimony and Bayesian phylogenetic analyses. Efficacy of mitochondrial markers was tested using operational (including distance-, tree-based measures and Barcode Gap) and evolutionary criteria (using the General Mixed Yule-coalescent Model) for species delimitation. We propose that the cox1 marker can potentially overestimate analyses of biodiversity and thus might not be the preferred marker for DNA species identification and delimitation methods in Oedothorax. Instead, our results suggest that the 16S marker appears to be a promising candidate for such endeavour. Evaluating the contribution and suitability of markers to the re-identification of species, measured by their recovery of well established morphological species, is critical for future studies and for reliable results in species identification in spiders.

Nuclear gene phylogeography using PHASE: dealing with unresolved genotypes, lost alleles, and systematic bias in parameter estimation

Background

A widely-used approach for screening nuclear DNA markers is to obtain sequence data and use bioinformatic algorithms to estimate which two alleles are present in heterozygous individuals. It is common practice to omit unresolved genotypes from downstream analyses, but the implications of this have not been investigated. We evaluated the haplotype reconstruction method implemented by PHASE in the context of phylogeographic applications. Empirical sequence datasets from five non-coding nuclear loci with gametic phase ascribed by molecular approaches were coupled with simulated datasets to investigate three key issues: (1) haplotype reconstruction error rates and the nature of inference errors, (2) dataset features and genotypic configurations that drive haplotype reconstruction uncertainty, and (3) impacts of omitting unresolved genotypes on levels of observed phylogenetic diversity and the accuracy of downstream phylogeographic analyses.

Results

We found that PHASE usually had very low false-positives (i.e., a low rate of confidently inferring haplotype pairs that were incorrect). The majority of genotypes that could not be resolved with high confidence included an allele occurring only once in a dataset, and genotypic configurations involving two low-frequency alleles were disproportionately represented in the pool of unresolved genotypes. The standard practice of omitting unresolved genotypes from downstream analyses can lead to considerable reductions in overall phylogenetic diversity that is skewed towards the loss of alleles with larger-than-average pairwise sequence divergences, and in turn, this causes systematic bias in estimates of important population genetic parameters.

Conclusions

A combination of experimental and computational approaches for resolving phase of segregating sites in phylogeographic applications is essential. We outline practical approaches to mitigating potential impacts of computational haplotype reconstruction on phylogeographic inferences. With targeted application of laboratory procedures that enable unambiguous phase determination via physical isolation of alleles from diploid PCR products, relatively little investment of time and effort is needed to overcome the observed biases.

Inference of population history by coupling exploratory and model-driven phylogeographic analyses

Understanding the nature, timing and geographic context of historical events and population processes that shaped the spatial distribution of genetic diversity is critical for addressing questions relating to speciation, selection, and applied conservation management. Cladistic analysis of gene trees has been central to phylogeography, but when coupled with approaches that make use of different components of the information carried by DNA sequences and their frequencies, the strength and resolution of these inferences can be improved. However, assessing concordance of inferences drawn using different analytical methods or genetic datasets, and integrating their outcomes, can be challenging. Here we overview the strengths and limitations of different types of genetic data, analysis methods, and approaches to historical inference. We then turn our attention to the potentially synergistic interactions among widely-used and emerging phylogeographic analyses, and discuss some of the ways that spatial and temporal concordance among inferences can be assessed. We close this review with a brief summary and outlook on future research directions.

How robust are "isolation with migration" analyses to violations of the im model? A simulation study

Methods developed over the past decade have made it possible to estimate molecular demographic parameters such as effective population size, divergence time, and gene flow with unprecedented accuracy and precision. However, they make simplifying assumptions about certain aspects of the species' histories and the nature of the genetic data, and it is not clear how robust they are to violations of these assumptions. Here, we use simulated data sets to examine the effects of a number of violations of the "Isolation with Migration" (IM) model, including intralocus recombination, population structure, gene flow from an unsampled species, linkage among loci, and divergent selection, on demographic parameter estimates made using the program IMA. We also examine the effect of having data that fit a nucleotide substitution model other than the two relatively simple models available in IMA. We find that IMA estimates are generally quite robust to small to moderate violations of the IM model assumptions, comparable with what is often encountered in real-world scenarios. In particular, population structure within species, a condition encountered to some degree in virtually all species, has little effect on parameter estimates even for fairly high levels of structure. Likewise, most parameter estimates are robust to significant levels of recombination when data sets are pared down to apparently nonrecombining blocks, although substantial bias is introduced to several estimates when the entire data set with recombination is included. In contrast, a poor fit to the nucleotide substitution model can result in an increased error rate, in some cases due to a predictable bias and in other cases due to an increase in variance in parameter estimates among data sets simulated under the same conditions.