Inferring introduction routes of invasive species using approximate Bayesian computation on microsatellite data

Determining the routes of introduction provides not only information about the history of an invasion process, but also information about the origin and construction of the genetic composition of the invading population. It remains difficult, however, to infer introduction routes from molecular data because of a lack of appropriate methods. We evaluate here the use of an approximate Bayesian computation (ABC) method for estimating the probabilities of introduction routes of invasive populations based on microsatellite data. We considered the crucial case of a single source population from which two invasive populations originated either serially from a single introduction event or from two independent introduction events. Using simulated datasets, we found that the method gave correct inferences and was robust to many erroneous beliefs. The method was also more efficient than traditional methods based on raw values of statistics such as assignment likelihood or pairwise F(ST). We illustrate some of the features of our ABC method, using real microsatellite datasets obtained for invasive populations of the western corn rootworm, Diabrotica virgifera virgifera. Most computations were performed with the DIYABC program (http://www1.montpellier.inra.fr/CBGP/diyabc/)

Lineage Divergence and Historical Gene Flow in the Chinese Horseshoe Bat (Rhinolophus sinicus)

PMCID: PMC3581519This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Secondary contact and admixture between independently invading populations of the Western corn rootworm, diabrotica virgifera virgifera in Europe

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a Bayesian analysis of the population structure and in an approximate Bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed

Lack of Detectable HIV-1 Molecular Evolution during Suppressive Antiretroviral Therapy

A better understanding of changes in HIV-1 population genetics with combination antiretroviral therapy (cART) is critical for designing eradication strategies. We therefore analyzed HIV-1 genetic variation and divergence in patients' plasma before cART, during suppression on cART, and after viral rebound. Single-genome sequences of plasma HIV-1 RNA were obtained from HIV-1 infected patients prior to cART (N = 14), during suppression on cART (N = 14) and/or after viral rebound following interruption of cART (N = 5). Intra-patient population diversity was measured by average pairwise difference (APD). Population structure was assessed by phylogenetic analyses and a test for panmixia. Measurements of intra-population diversity revealed no significant loss of overall genetic variation in patients treated for up to 15 years with cART. A test for panmixia, however, showed significant changes in population structure in 2/10 patients after short-term cART (<1 year) and in 7/10 patients after long-term cART (1-15 years). The changes consisted of diverse sets of viral variants prior to cART shifting to populations containing one or more genetically uniform subpopulations during cART. Despite these significant changes in population structure, rebound virus after long-term cART had little divergence from pretherapy virus, implicating long-lived cells infected before cART as the source for rebound virus. The appearance of genetically uniform virus populations and the lack of divergence after prolonged cART and cART interruption provide strong evidence that HIV-1 persists in long-lived cells infected before cART was initiated, that some of these infected cells may be capable of proliferation, and that on-going cycles of viral replication are not evident

Environmental variables, habitat discontinuity and life history shaping the genetic structure of Pomatoschistus marmoratus

Coastal lagoons are semi-isolated ecosystems exposed to wide fluctuations of environmental conditions and showing habitat fragmentation. These features may play an important role in separating species into different populations, even at small spatial scales. In this study, we evaluate the concordance between mitochondrial (previous published data) and nuclear data analyzing the genetic variability of Pomatoschistus marmoratus in five localities, inside and outside the Mar Menor coastal lagoon (SE Spain) using eight microsatellites. High genetic diversity and similar levels of allele richness were observed across all loci and localities, although significant genic and genotypic differentiation was found between populations inside and outside the lagoon. In contrast to the FST values obtained from previous mitochondrial DNA analyses (control region), the microsatellite data exhibited significant differentiation among samples inside the Mar Menor and between lagoonal and marine samples. This pattern was corroborated using Cavalli-Sforza genetic distances. The habitat fragmentation inside the coastal lagoon and among lagoon and marine localities could be acting as a barrier to gene flow and contributing to the observed genetic structure. Our results from generalized additive models point a significant link between extreme lagoonal environmental conditions (mainly maximum salinity) and P. marmoratus genetic composition. Thereby, these environmental features could be also acting on genetic structure of coastal lagoon populations of P. marmoratus favoring their genetic divergence. The mating strategy of P. marmoratus could be also influencing our results obtained from mitochondrial and nuclear DNA. Therefore, a special consideration must be done in the selection of the DNA markers depending on the reproductive strategy of the species

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

BACKGROUND: The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. RESULTS: Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. CONCLUSION: ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results

Genetic characterisation of farmed rainbow trout in Norway: intra- and inter-strain variation reveals potential for identification of escapees

<p>Abstract</p> <p>Background</p> <p>The rainbow trout (<it>Oncorhynchus mykiss</it>) is one of the most important aquaculture species in the world, and Norway is one of the largest producers. The present study was initiated in response to a request from the Norwegian police authority to identify the farm of origin for 35 escaped rainbow trout captured in a fjord. Eleven samples, each consisting of approximately 47 fish, were collected from the three farms operating in the fjord where the escapees were captured. In order to gain a better general understanding of the genetic structure of rainbow trout strains used in Norwegian aquaculture, seven samples (47 fish per sample) were collected from six farms located outside the region where the escapees were captured. All samples, including the escapees, were genotyped with 12 microsatellite loci.</p> <p>Results</p> <p>All samples displayed considerable genetic variability at all loci (mean number of alleles per locus per sample ranged from 5.4–8.6). Variable degrees of genetic differentiation were observed among the samples, with pair-wise <it>F</it>_STvalues ranging from 0–0.127. Self-assignment tests conducted among the samples collected from farms outside the fjord where the escapees were observed gave an overall correct assignment of 82.5%, demonstrating potential for genetic identification of escapees. In the "real life" assignment of the 35 captured escapees, all were excluded from two of the samples included as controls in the analysis, and 26 were excluded from the third control sample. In contrast, only 1 of the escapees was excluded from the 11 pooled samples collected on the 3 farms operating in the fjord.</p> <p>Conclusion</p> <p>Considerable genetic variation exists within and among rainbow trout strains farmed in Norway. Together with modern statistical methods, this will provide commercial operators with a tool to monitor breeding and fish movements, and management authorities with the ability to identify the source of escapees. The data generated in this study were used by the Norwegian police to initiate an investigation of the company operating the three farms in the fjord where escapees were observed.</p

New Insight into the History of Domesticated Apple: Secondary Contribution of the European Wild Apple to the Genome of Cultivated Varieties

The apple is the most common and culturally important fruit crop of temperate areas. The elucidation of its origin and domestication history is therefore of great interest. The wild Central Asian species Malus sieversii has previously been identified as the main contributor to the genome of the cultivated apple (Malus domestica), on the basis of morphological, molecular, and historical evidence. The possible contribution of other wild species present along the Silk Route running from Asia to Western Europe remains a matter of debate, particularly with respect to the contribution of the European wild apple. We used microsatellite markers and an unprecedented large sampling of five Malus species throughout Eurasia (839 accessions from China to Spain) to show that multiple species have contributed to the genetic makeup of domesticated apples. The wild European crabapple M. sylvestris, in particular, was a major secondary contributor. Bidirectional gene flow between the domesticated apple and the European crabapple resulted in the current M. domestica being genetically more closely related to this species than to its Central Asian progenitor, M. sieversii. We found no evidence of a domestication bottleneck or clonal population structure in apples, despite the use of vegetative propagation by grafting. We show that the evolution of domesticated apples occurred over a long time period and involved more than one wild species. Our results support the view that self-incompatibility, a long lifespan, and cultural practices such as selection from open-pollinated seeds have facilitated introgression from wild relatives and the maintenance of genetic variation during domestication. This combination of processes may account for the diversification of several long-lived perennial crops, yielding domestication patterns different from those observed for annual species

Estimating Parameters of Speciation Models Based on Refined Summaries of the Joint Site-Frequency Spectrum

Understanding the processes and conditions under which populations diverge to give rise to distinct species is a central question in evolutionary biology. Since recently diverged populations have high levels of shared polymorphisms, it is challenging to distinguish between recent divergence with no (or very low) inter-population gene flow and older splitting events with subsequent gene flow. Recently published methods to infer speciation parameters under the isolation-migration framework are based on summarizing polymorphism data at multiple loci in two species using the joint site-frequency spectrum (JSFS). We have developed two improvements of these methods based on a more extensive use of the JSFS classes of polymorphisms for species with high intra-locus recombination rates. First, using a likelihood based method, we demonstrate that taking into account low-frequency polymorphisms shared between species significantly improves the joint estimation of the divergence time and gene flow between species. Second, we introduce a local linear regression algorithm that considerably reduces the computational time and allows for the estimation of unequal rates of gene flow between species. We also investigate which summary statistics from the JSFS allow the greatest estimation accuracy for divergence time and migration rates for low (around 10) and high (around 100) numbers of loci. Focusing on cases with low numbers of loci and high intra-locus recombination rates we show that our methods for the estimation of divergence time and migration rates are more precise than existing approaches