Arnaud Estoup

Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis

Homoplasy has recently attracted the attention of population geneticists, as a consequence of the popularity of highly variable stepwise mutating markers such as microsatellites. Microsatellite alleles generally refer to DNA fragments of different size (electromorphs). Electromorphs are identical in state (i.e. have identical size), but are not necessarily identical by descent due to convergent mutation(s). Homoplasy occurring at microsatellites is thus referred to as size homoplasy. Using new analytical developments and computer simulations, we first evaluate the effect of the mutation rate, the mutation model, the effective population size and the time of divergence between populations on size homoplasy at the within and between population levels. We then review the few experimental studies that used various molecular techniques to detect size homoplasious events at some microsatellite loci. The relationship between this molecularly accessible size homoplasy size and the actual amount of size homoplasy is not trivial, the former being considerably influenced by the molecular structure of microsatellite core sequences. In a third section, we show that homoplasy at microsatellite electromorphs does not represent a significant problem for many types of population genetics analyses realized by molecular ecologists, the large amount of variability at microsatellite loci often compensating for their homoplasious evolution. The situations where size homoplasy may be more problematic involve high mutation rates and large population sizes together with strong allele size constraints.

Reconstructing routes of invasion using genetic data: why, how and so what?

Detailed knowledge about the geographical pathways followed by propagules from their source to the invading populations—referred to here as routes of invasion—provides information about the history of the invasion process and the origin and genetic composition of the invading populations. The reconstruction of invasion routes is required for defining and testing different hypotheses concerning the environmental and evolutionary factors responsible for biological invasions. In practical terms, it facilitates the design of strategies for controlling or preventing invasions. Most of our knowledge about the introduction routes of invasive species is derived from historical and observational data, which are often sparse, incomplete and, sometimes, misleading. In this context, population genetics has proved a useful approach for reconstructing routes of introduction, highlighting the complexity and the often counterintuitive nature of the true story. This approach has proved particularly useful since the recent development of new model‐based methods, such as approximate Bayesian computation, making it possible to make quantitative inferences in the complex evolutionary scenarios typically encountered in invasive species. In this review, we summarize some of the fundamental aspects of routes of invasion, explain why the reconstruction of these routes is useful for addressing both practical and theoretical questions, and comment on the various reconstruction methods available. Finally, we consider the main insights obtained to date from studies of invasion routes.

(CT) n and (GT) n microsatellites: a new class of genetic markers for Salmo trutta L. (brown trout)

Thirteen (GT)n and four (CT)n microsatellite loci (n = 10 or more and n = 20 or more, respectively) have been isolated from a partial genomic library of brown trout and sequenced. On average, a (GT)n repeat sequence occurs approximately every 23 kb and a (CT)n repeat sequence every 76 kb in brown trout genome. Primers for DNA amplifications using the polymerase chain reaction (PCR) were synthesized for three single locus microsatellites. Mendelian inheritance of the observed polymorphisms was confirmed in full-sib families. Four brown trout populations (10 unrelated individuals per population) were screened for polymorphism with these three microsatellite loci. The total number of alleles detected in the four populations is five at one locus, six at the other two microsatellite loci and is three, on average, per population. Heterozygosities range from 0.18 to 0.74. The largest differences in allelic frequencies occurred between the Mediterranean and the Atlantic populations: this result is congruent with previous allozymic data. The gene-centromere distances of the three microsatellite markers were determined on gynogenetic lines: post-reduction rates range from 0.17 to 0.60. For all the three microsatellite loci, the primers designed from brown trout sequences can be used in another closely related species of salmonid, the rainbow trout (Oncorhynchus mykiss). This last aspect supports the view that microsatellite markers may have wide application in genetic studies in salmonid species and fishes in general.

Bridgehead Effect in the Worldwide Invasion of the Biocontrol Harlequin Ladybird

Recent studies of the routes of worldwide introductions of alien organisms suggest that many widespread invasions could have stemmed not from the native range, but from a particularly successful invasive population, which serves as the source of colonists for remote new territories. We call here this phenomenon the invasive bridgehead effect. Evaluating the likelihood of such a scenario is heuristically challenging. We solved this problem by using approximate Bayesian computation methods to quantitatively compare complex invasion scenarios based on the analysis of population genetics (microsatellite variation) and historical (first observation dates) data. We applied this approach to the Harlequin ladybird Harmonia axyridis (HA), a coccinellid native to Asia that was repeatedly introduced as a biocontrol agent without becoming established for decades. We show that the recent burst of worldwide invasions of HA followed a bridgehead scenario, in which an invasive population in eastern North America acted as the source of the colonists that invaded the European, South American and African continents, with some admixture with a biocontrol strain in Europe. This demonstration of a mechanism of invasion via a bridgehead has important implications both for invasion theory (i.e., a single evolutionary shift in the bridgehead population versus multiple changes in case of introduced populations becoming invasive independently) and for ongoing efforts to manage invasions by alien organisms (i.e., heightened vigilance against invasive bridgeheads).

Clonal reproduction by males and females in the little fire ant

Sexual reproduction can lead to major conflicts between sexes and within genomes. Here we report an extreme case of such conflicts in the little fire ant Wasmannia auropunctata. We found that sterile workers are produced by normal sexual reproduction, whereas daughter queens are invariably clonally produced. Because males usually develop from unfertilized maternal eggs in ants and other haplodiploid species, they normally achieve direct fitness only through diploid female offspring. Hence, although the clonal production of queens increases the queens relatedness to reproductive daughters, it potentially reduces male reproductive success to zero. In an apparent response to this conflict between sexes, genetic analyses reveal that males reproduce clonally, most likely by eliminating the maternal half of the genome in diploid eggs. As a result, all sons have nuclear genomes identical to those of their father. The obligate clonal production of males and queens from individuals of the same sex effectively results in a complete separation of the male and female gene pools. These findings show that the haplodiploid sex-determination system provides grounds for the evolution of extraordinary genetic systems and new types of sexual conflict.

GENETIC ANALYSIS OF COMPLEX DEMOGRAPHIC SCENARIOS: SPATIALLY EXPANDING POPULATIONS OF THE CANE TOAD,BUFO MARINUS

Abstract Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection‐regression algorithm, to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration‐foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker, demographic, and historical features. The method presented here is the only robust, model‐based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Invasion of Europe by the western corn rootworm, Diabrotica virgifera virgifera: multiple transatlantic introductions with various reductions of genetic diversity

The early stages of invasion involve demographic bottlenecks that may result in lower genetic variation in introduced populations as compared to source population/s. Low genetic variability may decrease the adaptive potential of such populations in their new environments. Previous population genetic studies of invasive species have reported varying levels of losses of genetic variability in comparisons of source and invasive populations. However, intraspecific comparisons are required to assess more thoroughly the repeatability of genetic consequences of colonization events. Descriptions of invasive species for which multiple introductions from a single source population have been demonstrated may be particularly informative. The western corn rootworm (WCR), Diabrotica virgifera virgifera, native to North America and invasive in Europe, offers us an opportunity to analyse multiple introduction events within a single species. We investigated within‐ and between‐population variation at eight microsatellite markers in WCR in North America and Europe to investigate the routes by which WCR was introduced into Europe, and to assess the effect of introduction events on genetic variation. We detected five independent introduction events from the northern USA into Europe. The diversity loss following these introductions differed considerably between events, suggesting substantial variation in introduction, foundation and/or establishment conditions. Genetic variability at evolutionarily neutral loci does not seem to underlie the invasive success of WCR in Europe. We also showed that the introduction of WCR into Europe resulted in the redistribution of genetic variance from the intra‐ to the interpopulational level contrary to most examples of multiple introductions.

Ecological genetics of invasive alien species

There is growing realisation that integrating genetics and ecology is critical in the context of biological invasions, since the two are explicitly linked. So far, the focus of ecological genetics of invasive alien species (IAS) has been on determining the sources and routes of invasions, and the genetic make-up of founding populations, which is critical for defining and testing ecological and evolutionary hypotheses. However an ecological genetics approach can be extended to investigate questions about invasion success and impacts on native, recipient species. Here, we discuss recent progress in the field, provide overviews of recent methodological advances, and highlight areas that we believe are of particular interest for future research. First, we discuss the main insights from studies that have inferred source populations and invasion routes using molecular genetic data, with particular focus on the role of genetic diversity, adaptation and admixture in invasion success. Second, we consider how genetic tools can lead to a better understanding of patterns of dispersal, which is critical to predicting the spread of invasive species, and how studying invasions can shed light on the evolution of dispersal. Finally, we explore the potential for combining molecular genetic data and ecological network modelling to investigate community interactions such as those between predator and prey, and host and parasite. We conclude that invasions are excellent model systems for understanding the role of natural selection in shaping phenotypes and that an ecological genetics approach offers great potential for addressing fundamental questions in invasion biology.

Mating structure and nestmate relatedness in a communal bee, Andrena jacobi (Hymenoptera, Andrenidae), using microsatellites

Complex eusocial insect societies are generally matrifilial, suggesting kin selection has been of importance in their development. For simpler social systems, factors favouring their existence, in particular kin selection, have rarely been studied. Communal nesting is one of these simple social organizations, and is found in a diversity of insect species. To examine whether kin selection may play a role in the evolution and maintenance of communality, we estimated genetic relatedness of nestmate females of the facultatively communal bee, Andrena jacobi. Microsatellite loci were developed for this species and used to analyse individuals from two populations. Loci were variable, they were in heterozygote deficit and showed positive inbreeding coefficients. This may arise from nonrandom mating; previous observations (Paxton & Tengö 1996) indicate that a large proportion of females mate intranidally with nestmate males in their natal nests before first emerging. Nestmate relatedness was low, no different from zero for all loci in one population and for three of four loci in the other population. The large number of nestmates sharing a common nest (up to 594) may explain the low relatedness estimates, although relatedness was also independent of the number of females sharing a nest. Lack of inclusive fitness payoffs could constrain social evolution in this communal species.

The effect of RAD allele dropout on the estimation of genetic variation within and between populations

Inexpensive short‐read sequencing technologies applied to reduced representation genomes is revolutionizing genetic research, especially population genetics analysis, by allowing the genotyping of massive numbers of single‐nucleotide polymorphisms (SNP) for large numbers of individuals and populations. Restriction site–associated DNA (RAD) sequencing is a recent technique based on the characterization of genomic regions flanking restriction sites. One of its potential drawbacks is the presence of polymorphism within the restriction site, which makes it impossible to observe the associated SNP allele (i.e. allele dropout, ADO). To investigate the effect of ADO on genetic variation estimated from RAD markers, we first mathematically derived measures of the effect of ADO on allele frequencies as a function of different parameters within a single population. We then used RAD data sets simulated using a coalescence model to investigate the magnitude of biases induced by ADO on the estimation of expected heterozygosity and FST under a simple demographic model of divergence between two populations. We found that ADO tends to overestimate genetic variation both within and between populations. Assuming a mutation rate per nucleotide between 10−9 and 10−8, this bias remained low for most studied combinations of divergence time and effective population size, except for large effective population sizes. Averaging FST values over multiple SNPs, for example, by sliding window analysis, did not correct ADO biases. We briefly discuss possible solutions to filter the most problematic cases of ADO using read coverage to detect markers with a large excess of null alleles.