Strong intraspecific variation in genetic diversity and genetic differentiation in Daphnia magna: the effects of population turnover and population size

Theory predicts that genetic diversity and genetic differentiation may strongly vary among populations of the same species depending on population turnover and local population sizes. Yet, despite the importance of these predictions for evolutionary and conservation issues, empirical studies comparing high‐turnover and low‐turnover populations of the same species are scarce. In this study, we used Daphnia magna, a freshwater crustacean, as a model organism for such a comparison. In the southern/central part of its range, D. magna inhabits medium‐sized, stable ponds, whereas in the north, it occurs in small rock pools with strong population turnover. We found that these northern populations have a significantly lower genetic diversity and higher genetic differentiation compared to the southern/central populations. Total genetic diversity across populations was only about half and average within‐population diversity only about a third of that in southern/central populations. Moreover, an average southern population contains more genetic diversity than the whole metapopulation system in the north. We based our analyses both on silent sites and microsatellites. The similarity of our results despite the contrasting mutation rates of these markers suggests that the differences are caused by contemporary rather than by historical processes. Our findings show that variation in population turnover and population size may have a major impact on the genetic diversity and differentiation of populations, and hence may lead to differences in evolutionary processes like local adaptation, hybrid vigour and breeding system evolution in different parts of a species range.     

genodive version 3.0: Easy‐to‐use software for the analysis of genetic data of diploids and polyploids

genodive version 3.0 is a user‐friendly program for the analysis of population genetic data. This version presents a major update from the previous version and now offers a wide spectrum of different types of analyses. genodive has an intuitive graphical user interface that allows direct manipulation of the data through transformation, imputation of missing data, and exclusion and inclusion of individuals, population and/or loci. Furthermore, genodive seamlessly supports 15 different file formats for importing or exporting data from or to other programs. One major feature of genodive is that it supports both diploid and polyploid data, up to octaploidy (2n = 8x) for some analyses, but up to hexadecaploidy (2n = 16x) for other analyses. The different types of analyses offered by genodive include multiple statistics for estimating population differentiation (φ_ST, F_ST, F?_ST, G_ST, G?_ST, G??_ST, D_est, R_ST, ρ), analysis of molecular variance‐based K‐means clustering, Hardy–Weinberg equilibrium, hybrid index, population assignment, clone assignment, Mantel test, Spatial Autocorrelation, 23 ways of calculating genetic distances, and both principal components and principal coordinates analyses. A unique feature of genodive is that it can also open data sets with nongenetic variables, for example environmental data or geographical coordinates that can be included in the analysis. In addition, genodive makes it possible to run several external programs (lfmm , structure , instruct and vegan ) directly from its own user interface, avoiding the need for data reformatting and use of the command line. genodive is available for computers running Mac OS X 10.7 or higher and can be downloaded freely from: http://www.patrickmeirmans.com/software .     

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R_ST and F_ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad‐scale patterns of genetic structure and closely aligned with the R_ST network. The F_ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.     

Friends and family: A software program for identification of unrelated individuals from molecular marker data

The identification of related and unrelated individuals from molecular marker data is often difficult, particularly when no pedigree information is available and the data set is large. High levels of relatedness or inbreeding can influence genotype frequencies and thus genetic marker evaluation, as well as the accurate inference of hidden genetic structure. Identification of related and unrelated individuals is also important in breeding programmes, to inform decisions about breeding pairs and translocations. We present Friends and Family, a Windows executable program with a graphical user interface that identifies unrelated individuals from a pairwise relatedness matrix or table generated in programs such as coancestry and genalex . Friends and Family outputs a list of samples that are all unrelated to each other, based on a user‐defined relatedness cut‐off value. This unrelated data set can be used in downstream analyses, such as marker evaluation or inference of genetic structure. The results can be compared to that of the full data set to determine the effect related individuals have on the analyses. We demonstrate one of the applications of the program: how the removal of related individuals altered the Hardy–Weinberg equilibrium test outcome for microsatellite markers in an empirical data set. Friends and Family can be obtained from https://github.com/DeondeJager/Friends-and-Family .     

Spatial scale affects landscape genetic analysis of a wetland grasshopper

Most landscape genetic studies assess the impact of landscape elements on species' dispersal and gene flow. Many of these studies perform their analysis on all possible population pairs in a study area and do not explicitly consider the effects of spatial scale and population network topology on their results. Here, we examined the effects of spatial scale and population network topology on the outcome of a landscape genetic analysis. Additionally, we tested whether the relevant spatial scale of landscape genetic analysis could be defined by population network topology or by isolation‐by‐distance (IBD) patterns. A data set of the wetland grasshopper Stethophyma grossum, collected in a fragmented agricultural landscape, was used to analyse population network topology, IBD patterns and dispersal habitats, using least‐cost transect analysis. Landscape genetic analyses neglecting spatial scale and population network topology resulted in models with low fits, with which a most likely dispersal habitat could not be identified. In contrast, analyses considering spatial scale and population network topology resulted in high model fits by restricting landscape genetic analysis to smaller scales (0–3 km) and neighbouring populations, as represented by a Gabriel graph. These models also successfully identified a likely dispersal habitat of S. grossum. The above results suggest that spatial scale and potentially population network topology should be more explicitly considered in future landscape genetic analyses.     

Landscape genetic analysis of the tropical freshwater fish Mogurnda mogurnda (Eleotridae) in a monsoonal river basin: importance of hydrographic factors and population history

1. We performed spatial genetic analyses, incorporating landscape genetic methods using microsatellite data and phylogeographic analyses using mtDNA data, to identify the principal factors that determine population heterogeneity of the tropical freshwater fish, Mogurnda mogurnda, in the Daly River, northern Australia. We tested the individual and interactive effects of several environmental variables on spatial genetic patterns, including metrics relating to connectivity (i.e. stream distance, maximum stream gradient and elevation), habitat size (i.e. mean annual discharge) and a categorical variable relating to population history, as determined by mtDNA phylogeographic analyses. The Daly River is geomorphologically and hydrologically complex, and M. mogurnda has life history traits that limit its dispersal potential at river basin scales. Thus, we predicted that variables relating to connectivity would be the most important landscape factors driving population structure of the species. 2. Tree‐based phylogeographic analyses indicated four divergent mtDNA lineages within M. mogurnda in the Daly River, although three of the lineages were sympatric in various combinations and did not correspond with microsatellite groups identified by assignment tests. The allopatric mtDNA lineage detected in the uppermost part of the catchment was also identified as being highly differentiated by the microsatellite data, strongly suggesting that it may be a cryptic species. This site was therefore excluded from subsequent landscape genetic analyses. 3. Analyses of Molecular Variance indicated that M. mogurnda has a hierarchical population structure in the Daly River, thus supporting theoretical expectations that hierarchically arranged river habitats in dendritic systems impose hierarchal population structures on lotic species. 4. All landscape genetic analyses rejected stream distance, and supported stream gradient, as the major determinant of spatial genetic variation in M. mogurnda in the Daly River. Support for elevation as a determinant of spatial genetic patterns differed among the landscape genetic methods. Several of the landscape genetic methods also indicate that population history, including secondary contact between divergent and formerly allopatric genetic lineages, has a strong influence on spatial genetic patterns within M. mogurnda in the Daly River. 5. This study demonstrates the need to consider multiple environmental factors, especially factors relating to connectivity, and their interactions in spatial genetic analysis, rather than just geographic distance. Importantly, it demonstrates the need to account for population history and evolutionary divergences in landscape genetic analyses.     

Consequences of marine barriers for genetic diversity of the coral‐specialist yellowbar angelfish from the Northwestern Indian Ocean

Ocean circulation, geological history, geographic distance, and seascape heterogeneity play an important role in phylogeography of coral‐dependent fishes. Here, we investigate potential genetic population structure within the yellowbar angelfish (Pomacanthus maculosus) across the Northwestern Indian Ocean (NIO). We then discuss our results with respect to the above abiotic features in order to understand the contemporary distribution of genetic diversity of the species. To do so, restriction site‐associated DNA sequencing (RAD‐seq) was utilized to carry out population genetic analyses on P. maculosus sampled throughout the species’ distributional range. First, genetic data were correlated to geographic and environmental distances, and tested for isolation‐by‐distance and isolation‐by‐environment, respectively, by applying the Mantel test. Secondly, we used distance‐based and model‐based methods for clustering genetic data. Our results suggest the presence of two putative barriers to dispersal; one off the southern coast of the Arabian Peninsula and the other off northern Somalia, which together create three genetic subdivisions of P. maculosus within the NIO. Around the Arabian Peninsula, one genetic cluster was associated with the Red Sea and the adjacent Gulf of Aden in the west, and another cluster was associated with the Arabian Gulf and the Sea of Oman in the east. Individuals sampled in Kenya represented a third genetic cluster. The geographic locations of genetic discontinuities observed between genetic subdivisions coincide with the presence of substantial upwelling systems, as well as habitat discontinuity. Our findings shed light on the origin and maintenance of genetic patterns in a common coral reef fish inhabiting the NIO, and reinforce the hypothesis that the evolution of marine fish species in this region has likely been shaped by multiple vicariance events.     

Morphological and genetic distinctiveness of metallicolous and non-metallicolous populations of Armeria maritima s.l. (Plumbaginaceae) in Poland

Patterns of morphological, genetic and epigenetic variation (DNA methylation pattern) were investigated in metallicolous (M) and non‐metallicolous (NM) populations of Armeria maritima. A morphological study was carried out using plants from six natural populations grown in a greenhouse. Morphological variation was assessed using seven traits. On the basis of this study, three representative populations were selected for molecular analyses using metAFLP to study sequence‐ and methylation‐based DNA variation. Only one morphological trait (length of outer involucral bracts) was common to both metallicolous populations studied; however, the level of variation was sufficient to differentiate between M and NM populations. Molecular analyses showed the existence of naturally occurring epigenetic variation in A. maritima populations, as well as structuring into distinct between and within population components. We show that patterns of population genetic structure differed depending on the information used in the study. Analysis of sequence‐based information data demonstrates the presence of three well‐defined and genetically differentiated populations. Methylation‐based data show that two major groups of individuals are present, corresponding to the division into M and NM populations. These results were confirmed using different analytical approaches, which suggest that the DNA methylation pattern is similar in both M populations. We hypothesise that epigenetic processes may be involved in microevolution leading to development of M populations in A. maritima.     

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.     

Accounting for missing data in the estimation of contemporary genetic effective population size (Ne)

Theoretical models are often applied to population genetic data sets without fully considering the effect of missing data. Researchers can deal with missing data by removing individuals that have failed to yield genotypes and/or by removing loci that have failed to yield allelic determinations, but despite their best efforts, most data sets still contain some missing data. As a consequence, realized sample size differs among loci, and this poses a problem for unbiased methods that must explicitly account for random sampling error. One commonly used solution for the calculation of contemporary effective population size (N_e) is to calculate the effective sample size as an unweighted mean or harmonic mean across loci. This is not ideal because it fails to account for the fact that loci with different numbers of alleles have different information content. Here we consider this problem for genetic estimators of contemporary effective population size (N_e). To evaluate bias and precision of several statistical approaches for dealing with missing data, we simulated populations with known N_e and various degrees of missing data. Across all scenarios, one method of correcting for missing data (fixed‐inverse variance‐weighted harmonic mean) consistently performed the best for both single‐sample and two‐sample (temporal) methods of estimating N_e and outperformed some methods currently in widespread use. The approach adopted here may be a starting point to adjust other population genetics methods that include per‐locus sample size components.     

Scale‐dependent post‐establishment spread and genetic diversity in an invading mollusc in South America

Aims Our study aimed to characterize the dispersal dynamics and population genetic structure of the introduced golden mussel Limnoperna fortunei throughout its invaded range in South America and to determine how different dispersal methods, that is, human‐mediated dispersal and downstream natural dispersal, contribute to genetic variation among populations. Location Paraná–Uruguay–Río de la Plata watershed in Argentina, Brazil, Paraguay and Uruguay. Methods We performed genetic analyses based on a comprehensive sampling strategy encompassing 22 populations (N = 712) throughout the invaded range in South America, using the mitochondrial cytochrome c oxidase subunit I (COI) gene and eight polymorphic nuclear microsatellites. We employed both population genetics and phylogenetic analyses to clarify the dispersal dynamics and population genetic structure. Results We detected relatively high genetic differentiation between populations (F_ST = ?0.041 to 0.111 for COI, ?0.060 to 0.108 for microsatellites) at both fine and large geographical scales. Bayesian clustering and three‐dimensional factorial correspondence analyses consistently revealed two genetically distinct clusters, highlighting genetic discontinuities in the invaded range. Results of all genetic analyses suggest ship‐mediated ‘jump’ dispersal as the dominant mode of spread of golden mussels in South America, while downstream natural dispersal has had limited effects on contemporary genetic patterns. Main conclusions Our study provides new evidence that post‐establishment dispersal dynamics and genetic patterns vary across geographical scales. While ship‐mediated ‘jump’ dispersal dominates post‐establishment spread of golden mussels in South America, once colonies become established in upstream locations, larvae produced may be advected downstream to infill patchy distributions. Moreover, genetic structuring at fine geographical scales, especially within the same drainages, suggests a further detailed understanding of dynamics of larval dispersal and settlement in different water systems. Knowledge of the mechanisms by which post‐establishment spread occurs can, in some cases, be used to limit dispersal of golden mussels and other introduced species.     

Unexpected cryptic species diversity in the widespread coral Seriatopora hystrix masks spatial‐genetic patterns of connectivity

Mounting evidence of cryptic species in a wide range of taxa highlights the need for careful analyses of population genetic data sets to unravel within‐species diversity from potential interspecies relationships. Here, we use microsatellite loci and hierarchical clustering analysis to investigate cryptic diversity in sympatric and allopatric (separated by 450 km) populations of the widespread coral Seriatopora hystrix on the Great Barrier Reef. Structure analyses delimited unique genetic clusters that were confirmed by phylogenetic and extensive population‐level analyses. Each of four sympatric yet distinct genetic clusters detected within S. hystrix demonstrated greater genetic cohesion across regional scales than between genetic clusters within regions (<10 km). Moreover, the magnitude of genetic differentiation between different clusters (>0.620 G”_ST) was similar to the difference between S. hystrix clusters and the congener S. caliendrum (mean G”_ST 0.720). Multiple lines of evidence, including differences in habitat specificity, mitochondrial identity, Symbiodinium associations and morphology, corroborate the nuclear genetic evidence that these distinct clusters constitute different species. Hierarchical clustering analysis combined with more traditional population genetic methods provides a powerful approach for delimiting species and should be regularly applied to ensure that ecological and evolutionary patterns interpreted for single species are not confounded by the presence of cryptic species.     

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever‐widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F_ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity‐based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration‐specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context‐dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non‐human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high‐throughput technologies in population genetics.     

Utility of internally transcribed spacer region of rDNA (ITS) and β‐tubulin gene sequences to infer genetic diversity and migration patterns of Colletotrichum truncatum infecting Capsicum spp.

Anthracnose is among the most economically important diseases affecting pepper (Capsicum spp.) production in the tropics and subtropics. Of the three species of Colletotrichum implicated as causal agents of pepper anthracnose, C. truncatum is considered to be the most destructive in agro‐ecosystems worldwide. However, the genetic variation and the migration potential of C. truncatum infecting pepper are not known. Five populations were selected for study and a two‐locus (internally transcribed spacer region, ITS1‐5.8S‐ITS2, and β‐tubulin, β‐TUB) sequence data set was generated and used in the analyses. Sequences of the ITS region were less informative than β ‐ tubulin gene sequences based on comparisons of DNA polymorphism indices. Trinidad had the highest genetic diversity and also had the largest effective population size in pairwise comparisons with the other populations. The Trinidad population also demonstrated significant genetic differentiation from the other populations. AMOVA and STRUCTURE analyses both suggested significant genetic variation within populations more so than among populations. A consensus Maximum Likelihood tree based on β‐TUB gene sequences revealed very little intraspecific diversity for all isolates except for Trinidad. Two clades consisting solely of Trinidad isolates may have diverged earlier than the other isolates. There was also evidence of directional migration among the five populations. These findings may have a direct impact on the development of integrated disease management strategies to control C. truncatum infection in pepper.     

Why to account for finite sites in population genetic studies and how to do this with Jaatha 2.0

With the advent of next‐generation sequencing technologies, large data sets of several thousand loci from multiple conspecific individuals are available. Such data sets should make it possible to obtain accurate estimates of population genetic parameters, even for complex models of population history. In the analyses of large data sets, it is difficult to consider finite‐sites mutation models (FSMs). Here, we use extensive simulations to demonstrate that the inclusion of FSMs is necessary to avoid severe biases in the estimation of the population mutation rate θ, population divergence times, and migration rates. We present a new version of Jaatha, an efficient composite‐likelihood method for estimating demographic parameters from population genetic data and evaluate the usefulness of Jaatha in two biological examples. For the first application, we infer the speciation process of two wild tomato species, Solanum chilense and Solanum peruvianum. In our second application example, we demonstrate that Jaatha is readily applicable to NGS data by analyzing genome‐wide data from two southern European populations of Arabidopsis thaliana. Jaatha is now freely available as an R package from the Comprehensive R Archive Network (CRAN).     

Reduced genetic variation and strong genetic population structure in the freshwater killifish Valencia letourneuxi (Valenciidae) based on nuclear and mitochondrial markers

The genetic variation of the critically endangered Corfu killifish (Valencia letourneuxi), an endemic freshwater fish species of the western Balkans, was assessed for nine populations sampled in eight water systems in western continental Greece, the Peloponnese and the Ionian Island of Corfu, using mitochondrial and microsatellite markers. The analyses were based on data from three mtDNA regions (D‐loop, COI and 16S rRNA sequences) and 14 microsatellite loci. Samples from the congeneric species Valencia hispanica and the phylogenetically closely related species Aphanius fasciatus were also used in the study as outgroups. Both the mitochondrial and the microsatellite analyses revealed three distinct population groupings associated with the geographical distribution of the populations: one southern group occupying rivers draining to the Patraikos Gulf, the second one including the populations flowing into the Amvrakikos Gulf and the third, more northern group, including the other populations from rivers in Corfu Island and Epirus flowing into the Ionian Sea. Within these groupings there is limited genetic differentiation between populations; in addition, there is reduced intrapopulation genetic variation, evidenced by low heterozygosity values, number of alleles and haplotype diversity. In terms of taxonomic implications and appropriate management actions for conservation, our data suggest that the major population groups should be regarded at least as three distinct conservation units (CUs), with translocation and restocking actions to take place only within the geographical range of the CU concerned. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 334–349.     

Dispersal ability and habitat requirements determine landscape‐level genetic patterns in desert aquatic insects

Species occupying the same geographic range can exhibit remarkably different population structures across the landscape, ranging from highly diversified to panmictic. Given limitations on collecting population‐level data for large numbers of species, ecologists seek to identify proximate organismal traits—such as dispersal ability, habitat preference and life history—that are strong predictors of realized population structure. We examined how dispersal ability and habitat structure affect the regional balance of gene flow and genetic drift within three aquatic insects that represent the range of dispersal abilities and habitat requirements observed in desert stream insect communities. For each species, we tested for linear relationships between genetic distances and geographic distances using Euclidean and landscape‐based metrics of resistance. We found that the moderate‐disperser Mesocapnia arizonensis (Plecoptera: Capniidae) has a strong isolation‐by‐distance pattern, suggesting migration–drift equilibrium. By contrast, population structure in the flightless Abedus herberti (Hemiptera: Belostomatidae) is influenced by genetic drift, while gene flow is the dominant force in the strong‐flying Boreonectes aequinoctialis (Coleoptera: Dytiscidae). The best‐fitting landscape model for M. arizonensis was based on Euclidean distance. Analyses also identified a strong spatial scale‐dependence, where landscape genetic methods only performed well for species that were intermediate in dispersal ability. Our results highlight the fact that when either gene flow or genetic drift dominates in shaping population structure, no detectable relationship between genetic and geographic distances is expected at certain spatial scales. This study provides insight into how gene flow and drift interact at the regional scale for these insects as well as the organisms that share similar habitats and dispersal abilities.     

Source–sink dynamics structure a common montane mammal

Assessing the relative role of evolutionary processes on genetic diversity is critical for understanding species response to climatic change. However, many processes, independent of climate, can lead to the same genetic pattern. Because effective population size and gene flow are affected directly by abundance and dispersal, population ecology has the potential to profoundly influence patterns of genetic variation over microevolutionary timescales. Here, we use aDNA data and simulations to explore the influence of population ecology and Holocene climate change on genetic diversity of the Uinta ground squirrel (Spermophilus armatus). We examined phylochronology from three modern and two ancient populations spanning the climate transitions of the last 3000 years. Population genetic analyses based on summary statistics suggest that changes in genetic diversity and structure coincided with the Medieval Warm Period (MWP), c. 1000 years ago. Serial coalescent simulations allowed us to move beyond correlation with climate to statistically compare the likelihoods of alternative population histories given the observed data. The data best fit source–sink models that include large, mid‐elevation populations that exchange many migrants and small populations at the elevational extremes. While the MWP is likely to have reduced genetic diversity, our model‐testing approach revealed that MWP‐driven changes in genetic structure were not better supported for the range of models explored. Our results point to the importance of species ecology in understanding responses to climate, and showcase the use of ancient genetic data and simulation‐based inference for unraveling the relative roles of microevolutionary processes.     

Reconstruction of a windborne insect invasion using a particle dispersal model,historical wind data,and Bayesian analysis of genetic data

Understanding how invasive species establish and spread is vital for developing effective management strategies for invaded areas and identifying new areas where the risk of invasion is highest. We investigated the explanatory power of dispersal histories reconstructed based on local‐scale wind data and a regional‐scale wind‐dispersed particle trajectory model for the invasive seed chalcid wasp Megastigmus schimitscheki (Hymenoptera: Torymidae) in France. The explanatory power was tested by: (1) survival analysis of empirical data on M. schimitscheki presence, absence and year of arrival at 52 stands of the wasp's obligate hosts, Cedrus (true cedar trees); and (2) Approximate Bayesian analysis of M. schimitscheki genetic data using a coalescence model. The Bayesian demographic modeling and traditional population genetic analysis suggested that initial invasion across the range was the result of long‐distance dispersal from the longest established sites. The survival analyses of the windborne expansion patterns derived from a particle dispersal model indicated that there was an informative correlation between the M. schimitscheki presence/absence data from the annual surveys and the scenarios based on regional‐scale wind data. These three very different analyses produced highly congruent results supporting our proposal that wind is the most probable vector for passive long‐distance dispersal of this invasive seed wasp. This result confirms that long‐distance dispersal from introduction areas is a likely driver of secondary expansion of alien invasive species. Based on our results, management programs for this and other windborne invasive species may consider (1) focusing effort at the longest established sites and (2) monitoring outlying populations remains critically important due to their influence on rates of spread. We also suggest that there is a distinct need for new analysis methods that have the capacity to combine empirical spatiotemporal field data, genetic data, and environmental data to investigate dispersal and invasion.     

Low genetic differentiation among reed warbler Acrocephalus scirpaceus populations across Europe

Migratory birds generally have higher dispersal propensity than resident species and are thus expected to show less genetic differentiation. On the other hand, specific migration patterns may promote genetic structure, such as in situations where migratory divides impede random mixing of individuals. Here we investigated population genetic structure and gene flow patterns in a polytypic passerine, the reed warbler Acrocephalus scirpaceus which shows a migratory divide in central Europe. Using ten polymorphic microsatellite loci and extensive sampling we found low but significant overall genetic differentiation (F_ST=0.013, G’_ST=0.078, D=0.063). Hierarchical F‐statistics and barrier analyses showed low but significant genetic differentiation of Iberian populations, and also slight genetic differences across the migratory divide and between subspecies (A. s. scirpaceus and A. s. fuscus). Three individual‐based Bayesian methods, however, inferred a single genetic unit. Our study thus found low levels of genetic differentiation among reed warbler populations but this genetic differentiation was not pronounced enough to detect a clear population structure using the microsatellite data and no prior information on geographic location of the sampled individuals. This result indicates high levels of gene flow and suggests a possibly recent divergence of European populations after a rapid range expansion. Further studies are necessary to assess divergence times and to reveal the evolutionary history of the reed warbler populations.