Up-regulation of apoptotic- and cell survival-related gene pathways following exposures of western corn rootworm to B. thuringiensis crystalline pesticidal proteins in transgenic maize roots

BACKGROUND

Resistance of pest insect species to insecticides, including B. thuringiensis (Bt) pesticidal proteins expressed by transgenic plants, is a threat to global food security. Despite the western corn rootworm, Diabrotica virgifera virgifera, being a major pest of maize and having populations showing increasing levels of resistance to hybrids expressing Bt pesticidal proteins, the cell mechanisms leading to mortality are not fully understood.

RESULTS

Twenty unique RNA-seq libraries from the Bt susceptible D. v. virgifera inbred line Ped12, representing all growth stages and a range of different adult and larval exposures, were assembled into a reference transcriptome. Ten-day exposures of Ped12 larvae to transgenic Bt Cry3Bb1 and Gpp34/Tpp35Ab1 maize roots showed significant differential expression of 1055 and 1374 transcripts, respectively, compared to cohorts on non-Bt maize. Among these, 696 were differentially expressed in both Cry3Bb1 and Gpp34/Tpp35Ab1 maize exposures. Differentially-expressed transcripts encoded protein domains putatively involved in detoxification, metabolism, binding, and transport, were, in part, shared among transcripts that changed significantly following exposures to the entomopathogens Heterorhabditis bacteriophora and Metarhizium anisopliae. Differentially expressed transcripts in common between Bt and entomopathogen treatments encode proteins in general stress response pathways, including putative Bt binding receptors from the ATP binding cassette transporter superfamily. Putative caspases, pro- and anti-apoptotic factors, as well as endoplasmic reticulum (ER) stress-response factors were identified among transcripts uniquely up-regulated following exposure to either Bt protein.

CONCLUSIONS

Our study suggests that the up-regulation of genes involved in ER stress management and apoptotic progression may be important in determining cell fate following exposure of susceptible D. v. virgifera larvae to Bt maize roots. This study provides novel insights into insect response to Bt intoxication, and a possible framework for future investigations of resistance mechanisms.

Biases of STRUCTURE software when exploring introduction routes of invasive species

Population genetic methods are widely used to retrace the introduction routes of invasive species. The unsupervised Bayesian clustering algorithm implemented in STRUCTURE is amongst the most frequently used of these methods, but its ability to provide reliable information about introduction routes has never been assessed. We simulated microsatellite datasets to evaluate the extent to which the results provided by STRUCTURE were misleading for the inference of introduction routes. We focused on an invasion scenario involving one native and two independently introduced populations, because it is the sole scenario that can be rejected when obtaining a particular clustering with a STRUCTURE analysis at K = 2 (two clusters). Results were classified as "misleading" or "non-misleading". We investigated the influence of effective size, bottleneck severity and number of loci on the type and frequency of misleading results. We showed that misleading STRUCTURE results were obtained for 10% of all simulated datasets. Our results highlighted two categories of misleading output. The first occurs when the native population has a low level of diversity. In this case, the two introduced populations may be very similar, despite their independent introduction histories. The second category results from convergence issues in STRUCTURE for K = 2, with strong bottleneck severity and/or large numbers of loci resulting in high levels of differentiation between the three populations. Overall, the risk of being misled by STRUCTURE in the context of introduction routes inferences is moderate, but it is important to remain cautious when low genetic diversity or genuine multimodality between runs are involved.

EVOLUTION OF RESISTANCE IN CULEX PIPIENS: ALLELE REPLACEMENT AND CHANGING ENVIRONMENT

Fixation of adaptive mutations in populations is often constrained by pleiotropic fitness costs. The evolutionary pathways that compensate such fitness disadvantages are either the occurrence of modifier genes or replacement of the adaptive allele by less costly ones. In this context, 23 years of evolution of insecticide resistance genes in the mosquito Culex pipiens from southern France are analyzed. The aim of this study is to answer the following points. Is there a fitness cost associated with these resistance genes in natural populations? Does evolution proceed through allele replacement or through selection of modifiers? And finally, how do environmental changes affect the evolution of resistance genes? Samples from the same transect, crossing the boundary between an insecticide-treated and a nontreated area, are analyzed. Clinal analyses indicate a variable fitness cost among the resistance genes and show that allele replacement has been the primary mechanism of resistance evolution in this area. It is also shown that replacement was probably due to environmental changes corresponding to modification in pesticide-treatment intensity.

APPEARANCE AND SWEEP OF A GENE DUPLICATION: ADAPTIVE RESPONSE AND POTENTIAL FOR NEW FUNCTIONS IN THE MOSQUITO CULEX PIPIENS

Evolution of a new gene function is a fundamental process of adaptation. Gene duplication followed by divergence due to relaxed selection on redundant copies has been viewed as the predominant mechanism involved in this process. At a macroevolutionary scale, evidence for this scenario came from the analysis of sequences of genes families. However, even if several genetic models have described the different potential microevolutionary scenario for a new function to evolve, little is really known about the initial evolutionary dynamics of such processes. We analyze such early dynamics in natural populations of the mosquito Culex pipiens polymorphic for a duplication at Ace.1, a locus involved in insecticide resistance. The date of occurrence and the selective advantages of the duplication were estimated using frequency data. We propose a scenario where the spread of a duplication is driven, from the very beginning, by selection due to insecticide treatment.

Development of 23 polymorphic microsatellite loci in invasive silver wattle, Acacia dealbata (Fabaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed for silver wattle, Acacia dealbata (Fabaceae), which is both an ornamental and an invasive weed species. It is native to southeastern Australia and invasive in Europe, Africa, Asia, and the Americas.

METHODS AND RESULTS

The pyrosequencing of a microsatellite-enriched genomic DNA library of A. dealbata produced 33,290 sequences and allowed the isolation of 201 loci with a minimum of seven repeats of microsatellite motifs. Amplification tests led to the setup of two multiplex PCR mixes allowing the amplification of 21 loci. The polymorphism of these markers was evaluated on a sample of 32 individuals collected in southeastern Australia. The number of alleles and the expected heterozygosity varied between two and 11, and between 0.11 and 0.88, respectively.

CONCLUSIONS

The level of polymorphism of this set of 23 microsatellites is large enough to provide valuable information on the genetic structure and the invasion history of A. dealbata.

No clear effect of admixture between two European invading outbreaks of Diabrotica virgifera virgifera in natura

In this study, we challenged the hypothesis that admixture may have had a positive impact in the context of the European invasion of the western corn rootworm (WCR), Diabrotica virgifera virgifera, LeConte. This beetle was introduced in Europe from the USA several times since the 1980's. The multiple introductions of this major pest of cultivated corn led to the formation of two major outbreaks in North Western (NW) Italy and in Central and South Eastern (CSE) Europe that eventually merged into a secondary contact zone where insects from both outbreaks interbreed. We collected about 600 insects from this contact zone and genotyped them using 13 microsatellite markers. Three types of information were obtained from the collected individuals: (i) their survival under starvation; (ii) their admixed status, determined through a Bayesian method of genetic clustering and (iii) their mating probability, studied via the detection, isolation and genotyping of sperm in female spermathecae. Twenty six % and 12% of the individuals were assigned to the NW Italy or the CSE Europe parental types, respectively, and 23% and 39% to the F1 and backcross hybrid types, respectively. Globally, our results do not reveal any significant impact of the admixed status on the mating probability and on the choice of mating partners. However the admixed status had a sex- and sampling site-dependent effect on survival in adults under starvation. In addition sex had an effect on survival, with mortality hazard about 3 times larger in males than in females. The consequences of these findings for the evolution of the admixture zone of northern Italy are discussed.

Cannibalism in invasive, native and biocontrol populations of the harlequin ladybird

Background

Cannibalism is widespread in both vertebrates and invertebrates but its extent is variable between and within species. Cannibalism depends on population density and nutritional conditions, and could be beneficial during colonisation of new environments. Empirical studies are needed to determine whether this trait might facilitate invasion of a new area in natural systems. We investigated whether the propensity for cannibalism in H. axyridis differs both between native and invasive populations and between invasive populations from the core and from the front of the invasive area in Western Europe. We also compared the propensity for cannibalism of these natural populations with that of laboratory-reared biocontrol populations. We measured the cannibalism rates of eggs by first instar larvae and adult females at two different individual densities of ladybirds from three types of population (invasive, native and biocontrol), in laboratory-controlled conditions.

Results

Cannibalism was significantly greater in larvae from invasive populations compared to native or biocontrol populations, but there was no difference in cannibalism rates between populations from the core or front of the invaded range. Cannibalism was significantly lower in larvae from biocontrol populations compared to wild (invasive and native) populations. No differences in cannibalism rates of adult females were found between any populations. While high population density significantly increased cannibalism in both larvae and adults, the norm of reaction of cannibalism to individual density did not change significantly during the invasion and/or laboratory rearing processes.

Conclusion

This study is the first to provide evidence for a higher propensity for cannibalism in invasive populations compared to native ones. Our experiments also shed light on the difference in cannibalism evolution with respect to life stages. However, we are still at an early stage in understanding the underlying mechanisms and several different research perspectives are needed to determine whether the higher propensity for cannibalism is a general feature of the invasion process.

Phenotypic trait changes in laboratory--reared colonies of the maize herbivore, Diabrotica virgifera virgifera

The North American and European maize pest Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae) was used to assess whether conditions of the natal field, subsequent laboratory rearing, or genetic population origin affect phenotypic traits of fitness, activity, or morphometrics. Standardized laboratory bioassays with large sample sizes revealed that none of the 16 tested traits, except crawling behaviours, appeared consistently stable across all seven tested colonies. Environmental conditions in the natal field of the F 0 generation affected trait averages of the subsequently reared F 1 generation in laboratory in ca. 47% of cases, and trait variability in 67% of cases. This was apparent for fitness and morphometrics, but less obvious for activity traits. Early generation laboratory rearing affected trait averages in ca. 56% of cases: morphometrics changed; fecundity and egg survival increased from F 1 to F 2. Trait variability increased or decreased in 38% of cases. Laboratory rearing for over more than 190 generations affected the trait averages in 60% of cases, reflected by decreases in flight activity and increases in body size, weight, and fecundity to some extent. It had little effect on trait variability, especially so for morphometric variability. The genetic population origin affected average levels of 55% and variability of 63% of phenotypic traits. A comparison among D. v. virgifera studies might be difficult if they use different populations or laboratory colonies. It is advised to consider possible effects of original field conditions, laboratory rearing, and population genetics when planning comparative studies targeting fitness, activity, or morphometric questions regarding Diabrotica species.

Rapid increase in dispersal during range expansion in the invasive ladybird Harmonia axyridis

The evolutionary trajectories associated with demographic, genetic and spatial disequilibrium have become an issue of growing interest in population biology. Invasive species provide unique opportunities to explore the impact of recent range expansion on life-history traits, making it possible to test for a spatial arrangement of dispersal abilities along the expanding range, in particular. We carried out controlled experiments in laboratory conditions to test the hypothesis of an increase in dispersal capacity with range expansion in Harmonia axyridis, a ladybird that has been invading Europe since 2001. We found a marked increase in the flight speed of the insects from the core to the front of the invasion range in two independent sampling transects. By contrast, we found that two other traits associated with dispersal (endurance and motivation to fly off) did not follow the same spatial gradient. Our results provide a striking illustration of the way in which predictable directional genetic changes may occur rapidly for some traits associated with dispersal during biological invasions. We discuss the consequences of our results for invasion dynamics and the evolutionary outcomes of spatially expanding populations.

Role of a γ-aminobutryic acid (GABA) receptor mutation in the evolution and spread of Diabrotica virgifera virgifera resistance to cyclodiene insecticides

The western corn rootworm, Diabrotica virgifera virgifera, is a damaging pest of cultivated corn that was controlled by applications of cyclodiene insecticides from the late 1940s until resistance evolved ∼10 years later. Range expansion from the western plains into eastern USA coincides with resistance development. An alanine to serine amino acid substitution within the Rdl subunit of the gamma-aminobutyric acid (GABA) receptor confers resistance to cyclodiene insecticides in many species. We found that the non-synonymous single nucleotide polymorphism (SNP) G/T at the GABA receptor cDNA position 838 (G/T(838)) of D. v. virgifera resulted in the alanine to serine change, and the codominant SNP allele T(838) was genetically linked to survival of beetles in aldrin bioassays. A phenotypic gradient of decreasing susceptibility from west to east was correlated with higher frequencies of the resistance-conferring T(838) allele in the eastern-most populations. This pattern exists in opposition to perceived selective pressures since the more eastern and most resistant populations probably experienced reduced exposure. The reasons for the observed distribution are uncertain, but historical records of the range expansion combined with the distribution of susceptible and resistant phenotypes and genotypes provide an opportunity to better understand factors affecting the species' range expansion.

Estimation of the dispersal of a major pest of maize by cline analysis of a temporary contact zone between two invasive outbreaks

Dispersal is a key factor in invasion and in the persistence and evolution of species. Despite the importance of estimates of dispersal distance, dispersal measurement remains a real methodological challenge. In this study, we characterized dispersal by exploiting a specific case of biological invasion, in which multiple introductions in disconnected areas lead to secondary contact between two differentiated expanding outbreaks. By applying cline theory to this ecological setting, we estimated σ, the standard deviation of the parent-offspring distance distribution, of the western corn rootworm, Diabrotica virgifera virgifera, one of the most destructive pests of maize. This species is currently invading Europe, and the two largest invasive outbreaks, in northern Italy and Central Europe, have recently formed a secondary contact zone in northern Italy. We identified vanishing clines at 12 microsatellite loci throughout the contact zone. By analysing both the rate of change of cline slope and the spatial variation of linkage disequilibrium at these markers, we obtained two σ estimates of about 20 km/generation(1/2). Simulations indicated that these estimates were robust to changes in dispersal kernels and differences in population density between the two outbreaks, despite a systematic weak bias. These estimates are consistent with the results of direct methods for measuring dispersal applied to the same species. We conclude that secondary contact resulting from multiple introductions is very useful for the inference of dispersal parameters and should be more widely used in other species.

First insights into the genetic diversity of the pinewood nematode in its native area using new polymorphic microsatellite loci

The pinewood nematode, Bursaphelenchus xylophilus, native to North America, is the causative agent of pine wilt disease and among the most important invasive forest pests in the East-Asian countries, such as Japan and China. Since 1999, it has been found in Europe in the Iberian Peninsula, where it also causes significant damage. In a previous study, 94 pairs of microsatellite primers have been identified in silico in the pinewood nematode genome. In the present study, specific PCR amplifications and polymorphism tests to validate these loci were performed and 17 microsatellite loci that were suitable for routine analysis of B. xylophilus genetic diversity were selected. The polymorphism of these markers was evaluated on nematodes from four field origins and one laboratory collection strain, all originate from the native area. The number of alleles and the expected heterozygosity varied between 2 and 11 and between 0.039 and 0.777, respectively. First insights into the population genetic structure of B. xylophilus were obtained using clustering and multivariate methods on the genotypes obtained from the field samples. The results showed that the pinewood nematode genetic diversity is spatially structured at the scale of the pine tree and probably at larger scales. The role of dispersal by the insect vector versus human activities in shaping this structure is discussed.

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.

Microsatellite and mitochondrial data provide evidence for a single major introduction for the Neartic leafhopper Scaphoideus titanus in Europe

Scaphoideus titanus, a leafhopper native to North America and invasive in Europe, is the vector of the Flavescence dorée phytoplasma, the causal agent of the most important form of grapevine yellows in European vineyards. We studied 10 polymorphic microsatellite loci and a 623 bp fragment of the mitochondrial cytochrome oxidase II gene in native S. titanus from north-eastern America and introduced European populations, to elucidate the colonization scenario. Consistent with their recent history, invasive European populations were less genetically diverse than American populations for both types of markers, suggesting a recent bottleneck. Significant isolation by distance was detected between American populations but not between European populations. None of the European mitochondrial haplotypes was found in the American vineyards, from which they are assumed to have originated. The precise source of the invasive S. titanus populations therefore remains unclear. Nevertheless, the high heterozygosity of North-East American populations (which contained 92% of the observed alleles) suggests that this region is part of the native range of S. titanus. Clustering population genetics analyses with microsatellite and mitochondrial data suggested that European populations originated from a single introduction event. Most of the introduced populations clustered with populations from Long Island, the Atlantic Coast winegrowing region in which Vitis aestivalis occurs.

Permanent genetic resources added to Molecular Ecology Resources Database 1 August 2011-30 September 2011

This article documents the addition of 299 microsatellite marker loci and nine pairs of single-nucleotide polymorphism (SNP) EPIC primers to the Molecular Ecology Resources (MER) Database. Loci were developed for the following species: Alosa pseudoharengus, Alosa aestivalis, Aphis spiraecola, Argopecten purpuratus, Coreoleuciscus splendidus, Garra gotyla, Hippodamia convergens, Linnaea borealis, Menippe mercenaria, Menippe adina, Parus major, Pinus densiflora, Portunus trituberculatus, Procontarinia mangiferae, Rhynchophorus ferrugineus, Schizothorax richardsonii, Scophthalmus rhombus, Tetraponera aethiops, Thaumetopoea pityocampa, Tuta absoluta and Ugni molinae. These loci were cross-tested on the following species: Barilius bendelisis, Chiromantes haematocheir, Eriocheir sinensis, Eucalyptus camaldulensis, Eucalyptus cladocalix, Eucalyptus globulus, Garra litaninsis vishwanath, Garra para lissorhynchus, Guindilla trinervis, Hemigrapsus sanguineus, Luma chequen. Guayaba, Myrceugenia colchagüensis, Myrceugenia correifolia, Myrceugenia exsucca, Parasesarma plicatum, Parus major, Portunus pelagicus, Psidium guayaba, Schizothorax richardsonii, Scophthalmus maximus, Tetraponera latifrons, Thaumetopoea bonjeani, Thaumetopoea ispartensis, Thaumetopoea libanotica, Thaumetopoea pinivora, Thaumetopoea pityocampa ena clade, Thaumetopoea solitaria, Thaumetopoea wilkinsoni and Tor putitora. This article also documents the addition of nine EPIC primer pairs for Euphaea decorata, Euphaea formosa, Euphaea ornata and Euphaea yayeyamana.

Inferring the origin of populations introduced from a genetically structured native range by approximate Bayesian computation: case study of the invasive ladybird Harmonia axyridis

Correct identification of the source population of an invasive species is a prerequisite for testing hypotheses concerning the factors responsible for biological invasions. The native area of invasive species may be large, poorly known and/or genetically structured. Because the actual source population may not have been sampled, studies based on molecular markers may generate incorrect conclusions about the origin of introduced populations. In this study, we characterized the genetic structure of the invasive ladybird Harmonia axyridis in its native area using various population genetic statistics and methods. We found that native area of H. axyridis most probably consisted of two geographically distinct genetic clusters located in eastern and western Asia. We then performed approximate Bayesian computation (ABC) analyses on controlled simulated microsatellite data sets to evaluate (i) the risk of selecting incorrect introduction scenarios, including admixture between sources, when the populations of the native area are genetically structured and sampling is incomplete and (ii) the ability of ABC analysis to minimize such risks by explicitly including unsampled populations in the scenarios compared. Finally, we performed additional ABC analyses on real microsatellite data sets to retrace the origin of biocontrol and invasive populations of H. axyridis, taking into account the possibility that the structured native area may have been incompletely sampled. We found that the invasive population in eastern North America, which has served as the bridgehead for worldwide invasion by H. axyridis, was probably formed by an admixture between the eastern and western native clusters. This admixture may have facilitated adaptation of the bridgehead population.

Weak spatial and temporal population genetic structure in the rosy apple aphid, Dysaphis plantaginea, in French apple orchards

We used eight microsatellite loci and a set of 20 aphid samples to investigate the spatial and temporal genetic structure of rosy apple aphid populations from 13 apple orchards situated in four different regions in France. Genetic variability was very similar between orchard populations and between winged populations collected before sexual reproduction in the fall and populations collected from colonies in the spring. A very small proportion of individuals (∼2%) had identical multilocus genotypes. Genetic differentiation between orchards was low (F(ST)<0.026), with significant differentiation observed only between orchards from different regions, but no isolation by distance was detected. These results are consistent with high levels of genetic mixing in holocyclic Dysaphis plantaginae populations (host alternation through migration and sexual reproduction). These findings concerning the adaptation of the rosy apple aphid have potential consequences for pest management.

Stratified dispersal and increasing genetic variation during the invasion of Central Europe by the western corn rootworm, Diabrotica virgifera virgifera

Invasive species provide opportunities for investigating evolutionary aspects of colonization processes, including initial foundations of populations and geographic expansion. Using microsatellite markers and historical information, we characterized the genetic patterns of the invasion of the western corn rootworm (WCR), a pest of corn crops, in its largest area of expansion in Europe: Central and South-Eastern (CSE) Europe. We found that the invaded area probably corresponds to a single expanding population resulting from a single introduction of WCR and that gene flow is geographically limited within the population. In contrast to what is expected in classical colonization processes, an increase in genetic variation was observed from the center to the edge of the outbreak. Control measures against WCR at the center of the outbreak may have decreased effective population size in this area which could explain this observed pattern of genetic variation. We also found that small remote outbreaks in southern Germany and north-eastern Italy most likely originated from long-distance dispersal events from CSE Europe. We conclude that the large European outbreak is expanding by stratified dispersal, involving both continuous diffusion and discontinuous long-distance dispersal. This latter mode of dispersal may accelerate the expansion of WCR in Europe in the future.

Genome-wide survey and analysis of microsatellites in nematodes, with a focus on the plant-parasitic species Meloidogyne incognita

Background

Microsatellites are the most popular source of molecular markers for studying population genetic variation in eukaryotes. However, few data are currently available about their genomic distribution and abundance across the phylum Nematoda. The recent completion of the genomes of several nematode species, including Meloidogyne incognita, a major agricultural pest worldwide, now opens the way for a comparative survey and analysis of microsatellites in these organisms.

Results

Using MsatFinder, the total numbers of 1-6 bp perfect microsatellites detected in the complete genomes of five nematode species (Brugia malayi, Caenorhabditis elegans, M. hapla, M. incognita, Pristionchus pacificus) ranged from 2,842 to 61,547, and covered from 0.09 to 1.20% of the nematode genomes. Under our search criteria, the most common repeat motifs for each length class varied according to the different nematode species considered, with no obvious relation to the AT-richness of their genomes. Overall, (AT)_n, (AG)_n and (CT)_n were the three most frequent dinucleotide microsatellite motifs found in the five genomes considered. Except for two motifs in P. pacificus, all the most frequent trinucleotide motifs were AT-rich, with (AAT)_n and (ATT)_n being the only common to the five nematode species. A particular attention was paid to the microsatellite content of the plant-parasitic species M. incognita. In this species, a repertoire of 4,880 microsatellite loci was identified, from which 2,183 appeared suitable to design markers for population genetic studies. Interestingly, 1,094 microsatellites were identified in 801 predicted protein-coding regions, 99% of them being trinucleotides. When compared against the InterPro domain database, 497 of these CDS were successfully annotated, and further assigned to Gene Ontology terms.

Conclusions

Contrasted patterns of microsatellite abundance and diversity were characterized in five nematode genomes, even in the case of two closely related Meloidogyne species. 2,245 di- to hexanucleotide loci were identified in the genome of M. incognita, providing adequate material for the future development of a wide range of microsatellite markers in this major plant parasite.

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.

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.

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

Inferring population history with DIY ABC: a user-friendly approach to approximate Bayesian computation

Summary: Genetic data obtained on population samples convey information about their evolutionary history. Inference methods can extract part of this information but they require sophisticated statistical techniques that have been made available to the biologist community (through computer programs) only for simple and standard situations typically involving a small number of samples. We propose here a computer program (DIY ABC) for inference based on approximate Bayesian computation (ABC), in which scenarios can be customized by the user to fit many complex situations involving any number of populations and samples. Such scenarios involve any combination of population divergences, admixtures and population size changes. DIY ABC can be used to compare competing scenarios, estimate parameters for one or more scenarios and compute bias and precision measures for a given scenario and known values of parameters (the current version applies to unlinked microsatellite data). This article describes key methods used in the program and provides its main features. The analysis of one simulated and one real dataset, both with complex evolutionary scenarios, illustrates the main possibilities of DIY ABC.

Availability: The software DIY ABC is freely available at http://www.montpellier.inra.fr/CBGP/diyabc.

Contact:j.cornuet@imperial.ac.uk

Supplementary information: Supplementary data are also available at http://www.montpellier.inra.fr/CBGP/diyabc

A core set of microsatellite markers for western corn rootworm (Coleoptera: Chrysomelidae) population genetics studies

Interest in the ecological and population genetics of the western corn rootworm, Diabrotica virgifera virgifera LeConte, has grown rapidly in the last few years in North America and Europe. This interest is a result of a number of converging issues related to the increasing difficulty in managing this pest and the need to characterize and understand gene flow in the context of insect resistance management. One of the key components needed for successful population genetics studies is the availability of suitable molecular markers. Using a standard group of microsatellite markers enables researchers from different laboratories to directly compare and share their data, reducing duplication of effort and facilitating collaborative work among laboratories. We screened 22 candidate microsatellite loci against five criteria to create a core set of microsatellite markers for D. v. virgifera population genetics studies. The criteria for inclusion were moderate to high polymorphism, unambiguous readability and repeatability, no evidence of null alleles, apparent selective neutrality, and no linkage between loci. Based on our results, we recommend six microsatellite markers to be included as a core set in future population genetics studies of D. v. virgifera along with any other microsatellite or genetic markers. As more microsatellites are developed, those meeting the criteria can be added to the core set. We encourage other groups of researchers with common interests in a particular insect species to develop their own core sets of markers for population genetics applications.

Absence of genetic divergence between western corn rootworms (Coleoptera: Chrysomelidae) resistant and susceptible to control by crop rotation

The western corn rootworm, Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae), is a major pest of corn, Zea mays L., in North America that has recently invaded Europe. A loss of ovipositional fidelity to cornfields has allowed the species to circumvent crop rotation as a means of control in part of its range in the United States. Analyses of variation at eight microsatellite loci provided no evidence for general genetic differentiation between samples of western corn rootworm collected in soybean, Glycine max (L.) Merr., fields and those collected in cornfields both inside and outside the rotation-resistance problem area. This result suggests that few or no barriers to gene flow exist between rotation-resistant and -susceptible rootworm populations. The implications of this result for the management of western corn rootworm in North America and Europe are discussed.

Within-species variability of the response to 20-hydroxyecdysone in peach-potato aphid (Myzus persicae Sulzer)

Phytoecdysteroids have been proposed as new tools for controlling crop pests because of their endocrine disruption and deterrent effects on insects and nematodes. There is increasing evidence of variability between taxa in sensitivity to phytoecdysteroids, but the genetic variability of this sensitivity within species is unknown. However, knowledge about this intraspecies variability is required for predicting evolution of the pest's response to new control methods. We assessed the variability of the response of the aphid Myzus persicae Sulzer, a major agricultural pest, to 20-hydroxyecdysone (20E). We determined the number of nymphs produced by six clones of M. persicae exposed to various concentrations of 20E and the capacity of these clones to detect 20E in choice experiments. High concentrations of 20E significantly decreased the number of nymphs produced for two clones and both increases and decreases in the number of offspring were detected at low concentrations. Two clones significantly avoided food with 20E, while one significantly preferred it, suggesting that 20E does not always act as a deterrent in this species. We conclude that genetic variability in the response to 20E exists in natural populations of M. persicae. The consequences of this finding on the sustainability of control methods using 20E are discussed.

Multiple Transatlantic Introductions of the Western Corn Rootworm

The invasion of Europe by the western corn rootworm, North America's most destructive corn pest, is ongoing and represents a serious threat to European agriculture. Because this pest was initially introduced in Central Europe, it was believed that subsequent outbreaks in Western Europe originated from this area. Using model-based Bayesian analyses of the genetic variability of the western corn rootworm, we demonstrate that this belief is false: There have been at least three independent introductions from North America during the past two decades. This result raises questions about changing circumstances that have enabled a sudden burst of transatlantic introductions.

Does insecticide resistance alone account for the low genetic variability of asexually reproducing populations of the peach-potato aphid Myzus persicae?

The typical life cycle of aphids includes several parthenogenetic generations and a single sexual generation (cyclical parthenogenesis), but some species or populations are totally asexual (obligate parthenogenesis). Genetic variability is generally low in these asexually reproducing populations, that is, few genotypes are spread over large geographic areas. Both genetic drift and natural selection are often invoked to account for this low genetic variability. The peach-potato aphid, Myzus persicae, which encompasses both cyclical and obligate parthenogens, has developed several insecticide resistance mechanisms as a consequence of intense insecticide use since the 1950s. We collected asexually reproducing M. persicae from oilseed rape and examined genetic variability at eight microsatellite loci and three insecticide resistance genes to determine whether their genetic structure was driven by drift and/or selection. We identified only 16 multilocus microsatellite genotypes among 255 individuals. One clone, which combined two insecticide resistance mechanisms, was frequently detected in all populations whatever their location over a large geographical area (the northern half of France). These unexpected findings suggest that drift is not the unique cause of this low variability. Instead, the intensification of both insecticide treatments and oilseed rape cultivation may have favored a few genotypes. Thus, we propose that selective pressures resulting from human activities have considerably modified the genetic structure of M. persicae populations in northern France in a relatively short period of time.

Asymmetry in host and parasitoid diffuse coevolution: when the red queen has to keep a finger in more than one pie

BACKGROUND: Coevolution between pairs of antagonistic species is generally considered an endless "arms race" between attack and defense traits to counteract the adaptive responses of the other species. PRESENTATION OF THE HYPOTHESIS: When more than two species are involved, diffuse coevolution of hosts and parasitoids could be asymmetric because consumers can choose their prey whereas preys do not choose their predator. This asymmetry may lead to differences in the rate of evolution of the antagonistic species in response to selection. The more long-standing the coevolution of a given pair of antagonistic populations, the higher should be the fitness advantage for the consumer. Therefore, the main prediction of the hypothesis is that the consumer trophic level is more likely to win the coevolution race. TESTING THE HYPOTHESIS: We propose testing the asymmetry hypothesis by focusing on the tritrophic system plant/aphid/aphid parasitoid. The analysis of the genetic variability in the virulence of several parasitoid populations and in the defenses of several aphid species or several clones of the same aphid species could be compared. Moreover, the analysis of the neutral population genetic structure of the parasitoid as a function of the aphid host, the plant host and geographic isolation may complement the detection of differences between host and parasitoid trophic specialization. IMPLICATIONS OF THE HYPOTHESIS: Genetic structures induced by the arms race between antagonistic species may be disturbed by asymmetry in coevolution, producing neither rare genotype advantages nor coevolutionary hotspots. Thus this hypothesis profoundly changes our understanding of coevolution and may have important implications in terms of pest management.

Fitness costs of insecticide resistance in natural breeding sites of the mosquito Culex pipiens

Genetic changes conferring adaptation to a new environment may induce a fitness cost in the previous environment. Although this prediction has been verified in laboratory conditions, few studies have tried to document this cost directly in natural populations. Here, we evaluated the pleiotropic effects of insecticide resistance on putative fitness components of the mosquito Culex pipiens. Experiments using different larval densities were performed during the summer in two natural breeding sites. Two loci that possess alleles conferring organophosphate (OP) resistance were considered: ace-1 coding for an acetylcholinesterase (AChE1, the OP target) and Ester, a ''super locus" including two closely linked loci coding for esterases A and B. Resistance ace-1 alleles coding for a modified AChE1 were associated with a longer development time and shorter wing length. The pleiotropic effects of two resistance alleles Ester1 and Ester4 coding for the overproduced esterases A1 and A4-B4, respectively, were more variable. Both A1 and A4-B4 reduced wing length, although only A1 was associated with a longer preimaginal stage. The fluctuating asymmetry (FA) of the wing did not respond to the presence or to the interaction of resistance alleles at the two loci at any of the density levels tested. Conversely, the FA of one wing section decreased when larval density increased. This may be the consequence of selection against less developmentally stable individuals. The results are discussed in relation to the local evolution of insecticide resistance genes.

Genetic study of Coris julis (Osteichthyes, Perciformes, Labridae) evolutionary history and dispersal abilities

Microsatellite markers have been used to study the genetic variability of rainbow wrasse (Coris julis) Mediterranean and Atlantic populations. Differentiation tests failed to reveal any significant genetic differentiation among samples from continental Portugal and the Azores, despite more than 1800 km of geographical separation. Preliminary results tended to indicate a significant genetic differentiation among Atlantic and Mediterranean samples. It also supported the specific status of Cape Verde populations (Coris atlantica). We compare these results with previous mtDNA analyses and propose a biogeographic scenario that could explain our results.

Spatial and temporal genetic variability in French populations of the peach-potato aphid, Myzus persicae

The peach-potato aphid, Myzus persicae (Sulzer), has a complex reproductive mode in which some lineages reproduce by continuous parthenogenesis, whereas others reproduce sexually once a year. The climate is thought to act directly on the reproductive mode, because sexual eggs are the only form that can resist frost in cold regions. Sexual reproduction necessitates an obligatory host alternation that may result in long-distance dispersal. Here, we examined the genetic variability at seven microsatellite loci of populations of M. persicae in France, where both reproductive modes occur. We provide clear genetic evidence that the breeding system affects genotypic variability, as cyclically parthenogenetic aphids are far more variable than their obligately parthenogenetic counterparts. A temporal decrease in genetic variability and a temporal genetic differentiation effect suggest the existence of selective factors that play an important role in shaping the genetic structure of M. persicae populations. Lastly, differences in the population structure between reproductive modes suggest that the migration associated with the change of host during sexual reproduction lowers the level of population differentiation.

Incidence of insecticide resistance alleles in sexually-reproducing populations of the peach-potato aphid Myzus persicae (Hemiptera: Aphididae) from southern France

Intensive chemical treatments have led to the development of a number of insecticide resistance mechanisms in the peach-potato aphid Myzus persicae (Sulzer). Some of these mechanisms are known to be associated with negative pleiotropic effects (resistance costs). Molecular and biochemical methods were used to determine the genotypes or phenotypes associated with four insecticide resistance mechanisms in single aphids from sexually-reproducing populations in southern France. The mechanisms considered were E4 and FE4 carboxylesterase overproduction, modified acetycholinesterase, and kdr and rdl resistance-associated mutations. A new method for determining individual kdr genotypes is presented. Almost all resistant individuals overproduced FE4 carboxylesterase, whereas modified acetylcholinesterase was rare. Both the kdr and rdl resistance mutations were present at high frequencies in French sexually-reproducing populations. The frequencies of insecticide resistance genes were compared before and after sexual reproduction in one peach orchard at Avignon to evaluate the potential impact of selection on the persistence of resistance alleles in the over-wintering phase. The frequencies of the kdr and rdl mutations varied significantly between autumn and spring sampling periods. The frequency of the kdr mutation increased, probably due to pyrethroid treatments at the end of the winter. Conversely, the frequency of the rdl mutation decreased significantly during winter, probably because of a fitness cost associated with this mutation.

Interspecific Utility of Microsatellites in Fish: A Case Study of (CT)(n) and (GT)(n) Markers in the Shanny Lipophrys pholis (Pisces: Blenniidae) and Their Use in Other Blennioidei

We report the development of new microsatellite markers that can be used for population analyses in the shanny Lipophrys pholis. The procedure involved the construction of a microsatellite-enriched genomic bank. Five (GT)(n) and (CT)(n) microsatellites have been characterized, four of which are polymorphic. The analysis of one population allowed us to verify their usefulness as markers in population studies. Moreover, interspecific amplifications have been performed using primers defined in other species to amplify Lipophrys pholis, or using the primers defined in Lipophrys pholis to amplify other species. We use these results to discuss the hypothesis that microsatellites are highly conserved in fish.