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/)

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

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

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

On the origin of the invasive olives (Olea europaea L., Oleaceae).

The olive tree (Olea europaea) has successfully invaded several regions in Australia and Pacific islands. Two olive subspecies (subspp. europaea and cuspidata) were first introduced in these areas during the nineteenth century. In the present study, we determine the origin of invasive olives and investigate the importance of historical effects on the genetic diversity of populations. Four invasive populations from Australia and Hawaii were characterized using eight nuclear DNA microsatellites, plastid DNA markers as well as ITS-1 sequences. Based on these data, their genetic similarity with native populations was investigated, and it was determined that East Australian and Hawaiian populations (subsp. cuspidata) have originated from southern Africa while South Australian populations (subsp. europaea) have mostly derived from western or central Mediterranean cultivars. Invasive populations of subsp. cuspidata showed significant loss of genetic diversity in comparison to a putative source population, and a recent bottleneck was evidenced in Hawaii. Conversely, invasive populations of subsp. europaea did not display significant loss of genetic diversity in comparison to a native Mediterranean population. Different histories of invasion were inferred for these two taxa with multiple cultivars introduced restoring gene diversity for europaea and a single successful founder event and sequential introductions to East Australia and then Hawaii for cuspidata. Furthermore, one hybrid (cuspidata x europaea) was identified in East Australia. The importance of hybridizations in the future evolution of the olive invasiveness remains to be investigated

High genetic diversity in a potentially vulnerable tropical tree species despite extreme habitat loss

10.1371/journal.pone.0082632PLoS ONE812-POLN

Tigers of Sundarbans in India: Is the Population a Separate Conservation Unit?

The Sundarbans tiger inhabits a unique mangrove habitat and are morphologically distinct from the recognized tiger subspecies in terms of skull morphometrics and body size. Thus, there is an urgent need to assess their ecological and genetic distinctiveness and determine if Sundarbans tigers should be defined and managed as separate conservation unit. We utilized nine microsatellites and 3 kb from four mitochondrial DNA (mtDNA) genes to estimate genetic variability, population structure, demographic parameters and visualize historic and contemporary connectivity among tiger populations from Sundarbans and mainland India. We also evaluated the traits that determine exchangeability or adaptive differences among tiger populations. Data from both markers suggest that Sundarbans tiger is not a separate tiger subspecies and should be regarded as Bengal tiger (P. t. tigris) subspecies. Maximum likelihood phylogenetic analyses of the mtDNA data revealed reciprocal monophyly. Genetic differentiation was found stronger for mtDNA than nuclear DNA. Microsatellite markers indicated low genetic variation in Sundarbans tigers (He= 0.58) as compared to other mainland populations, such as northern and Peninsular (Hebetween 0.67- 0.70). Molecular data supports migration between mainland and Sundarbans populations until very recent times. We attribute this reduction in gene flow to accelerated fragmentation and habitat alteration in the landscape over the past few centuries. Demographic analyses suggest that Sundarbans tigers have diverged recently from peninsular tiger population within last 2000 years. Sundarbans tigers are the most divergent group of Bengal tigers, and ecologically non-exchangeable with other tiger populations, and thus should be managed as a separate "evolutionarily significant unit" (ESU) following the adaptive evolutionary conservation (AEC) concept.Wildlife Institute of India, Dehra Dun (India)

Panmixia in a fragmented and unstable environment: the hydrothermal shrimp Rimicaris exoculata disperses extensively along the Mid-Atlantic ridge

Dispersal plays a fundamental role in the evolution and persistence of species, and especially for species inhabiting extreme, ephemeral and highly fragmented habitats as hydrothermal vents. The Mid-Atlantic Ridge endemic shrimp species Rimicaris exoculata was studied using microsatellite markers to infer connectivity along the 7100-Km range encompassing the sampled sites. Astonishingly, no genetic differentiation was found between individuals from the different geographic origins, supporting a scenario of widespread large-scale dispersal despite the habitat distance and fragmentation. We hypothesize that delayed metamorphosis associated to temperature differences or even active directed migration dependent on physical and/or chemical stimuli could explain these results and warrant further studies on adaptation and dispersal mechanisms

Genetic variability of the grey wolf Canis lupus in the Caucasus in comparison with Europe and the Middle East: distinct or intermediary population?

Despite continuous historical distribution of the grey wolf (Canis lupus) throughout Eurasia, the species displays considerable morphological differentiation that resulted in delimitation of a number of subspecies. However, these morphological discontinuities are not always consistent with patterns of genetic differentiation. Here we assess genetic distinctiveness of grey wolves from the Caucasus (a region at the border between Europe and West Asia) that have been classified as a distinct subspecies C. l. cubanensis. We analysed their genetic variability based on mtDNA control region, microsatellite loci and genome-wide SNP genotypes (obtained for a subset of the samples), and found similar or higher levels of genetic diversity at all these types of loci as compared with other Eurasian populations. Although we found no evidence for a recent genetic bottleneck, genome-wide linkage disequilibrium patterns suggest a long-term demographic decline in the Caucasian population – a trend consistent with other Eurasian populations. Caucasian wolves share mtDNA haplotypes with both Eastern European and West Asian wolves, suggesting past or ongoing gene flow. Microsatellite data also suggest gene flow between the Caucasus and Eastern Europe. We found evidence for moderate admixture between the Caucasian wolves and domestic dogs, at a level comparable with other Eurasian populations. Taken together, our results show that Caucasian wolves are not genetically isolated from other Eurasian populations, share with them the same demographic trends, and are affected by similar conservation problems