- Research article
- Open Access
Tracking the global dispersal of a cosmopolitan insect pest, the peach potato aphid
BMC Ecologyvolume 9, Article number: 13 (2009)
Global commerce and human transportation are responsible for the range expansion of various insect pests such as the plant sucking aphids. High resolution DNA markers provide the opportunity to examine the genetic structure of aphid populations, identify aphid genotypes and infer their evolutionary history and routes of expansion which is of value in developing management strategies. One of the most widespread aphid species is the peach-potato aphid Myzus persicae, which is considered as a serious pest on various crops in many parts of the world. The present study examined the genetic variation of this aphid at a world scale and then related this to distribution patterns. In particular, 197 aphid parthenogenetic lineages from around the world were analysed with six microsatellite loci.
Bayesian clustering and admixture analysis split the aphid genotypes into three genetic clusters: European M. persicae persicae, New Zealand M. persicae persicae and Global M. persicae nicotianae. This partition was supported by FST and genetic distance analyses. The results showed two further points, a possible connection between genotypes found in the UK and New Zealand and globalization of nicotianae associated with colonisation of regions where tobacco is not cultivated. In addition, we report the presence of geographically widespread clones and for the first time the presence of a nicotianae genotype in the Old and New World. Lastly, heterozygote deficiency was detected in some sexual and asexual populations.
The study revealed important genetic variation among the aphid populations we examined and this was partitioned according to region and host-plant. Clonal selection and gene flow between sexual and asexual lineages are important factors shaping the genetic structure of the aphid populations. In addition, the results reflected the globalization of two subspecies of M. persicae with successful clones being spread at various scales throughout the world. A subspecies appears to result from direct selection on tobacco plants. This information highlights the ultimate ability of a polyphagous aphid species to generate and maintain ecologically successful gene combinations through clonal propagation and the role of human transportation and global commerce for expanding their range.
Myzus persicae (Sulzer) (Hemiptera: Aphididae) is an exceptional species in many respects. It is extremely polyphagous, highly efficient as a plant-virus vector and one of the most widespread insect pests, as it has been recorded on all continents where crops are grown . The species has a typical aphid annual cycle (cyclical parthenogenesis), i.e., a sexual generation on peach during winter and spring, alternating with many parthenogenetic (all female) generations during spring on peach and on various crop and non-crop annual plant hosts in summer and autumn. The sexual generation may be lost either totally (obligate parthenogenesis) or partially (functional parthenogenesis when a few sexual forms are produced). Genotypes with different reproductive strategies can occur sympatrically in peach growing areas where populations on summer crops consist of new recombinants that have migrated from peach and old clones that survived the previous winter(s) parthenogenetically on winter hosts. Their proportions depend on the availability of peach trees for the sexual phase and the severity of winter which mostly affects the parthenogenetical overwintering [2–4]. The plasticity in the mode of reproduction is a great biological advantage for M. persicae, because as a species, it is able to adapt to different climatic conditions in terms of day length and temperature . Sexual reproduction provides the advantages of cold hardy eggs and new gene combinations in the subsequent generations. Asexual reproduction has the advantage of maintaining successful gene combinations and success in temperate regions where green bridges are available. The effect of reproductive strategy on the population structure of M. persicae has gained much attention recently. Studies conducted in Australia [4, 6] and Europe [7, 8] showed that clonal diversity was greatest in populations capable of a sexual phase compared to parthenogenetic ones. The most extreme case of reduced variability has been documented from Scotland, where the majority of the stable long-term population appears to consist of only three genotypes [9, 10].
Myzus persicae exhibits strong selection with respect to host-plant adaptation on tobacco. Tobacco-feeding populations show consistent morphological differences from those on other crops, regardless of the mode of reproduction or origin [11–13]. The subspecies name, Myzus persicae nicotianae, has been given to the tobacco population  and genetic differences have also been reported between the specialist (nicotianae) and the generalist (persicae) subspecies . The two taxa experience multifarious divergent selection, i.e., selection against cross-host migrants and their subsequent generations, which is crucial for the maintenance of host specialization [16, 17]. In outdoor choice experiments with winged females, it has been shown that the two taxa have evolved an improved host recognition mechanism which is based on chemical cues perceived prior to the initiation of feeding . Gene-flow between the two taxa is reduced due to differences in the mode of reproduction (asexual vs. sexual) or to the existence of prezygotic reproductive isolation mechanism in sexual populations .
Myzus persicae has developed multiple insecticide resistance mechanisms which have spread to many parts of the world due to aphid migration or to human transport activity . The international trade in plants offers considerable potential for the widespread distribution of insect pests  and M. persicae is ideally suited for this. It is associated with many transportable crop hosts and its primary host, the peach tree, has been spread throughout the world. Knowledge of the inter-regional or inter-continental dispersal routes of M. persicae and the between region genetic variation will be helpful in elucidating several aspects of its ecology that offer explanations for the spread and persistence in heterogeneous environments of successful genotypes. This information could also help to predict the evolution and spread of insecticide resistance mechanisms at both regional and global scale, to locate sources of host-plant resistance or biological control agents and to make inferences about the relative fitness and persistence of the pest. To our knowledge no other study has attempted to examine the genetic variation of M. persicae or the spread of certain genotypes at such a global scale. Some studies have focused mainly on insecticide resistance mechanisms. For instance, a phylogenetic analysis of a fragment of the para-sodium channel gene flanking the kdr and super-kdr mutations in samples from various countries and continents suggested multiple independent origins of both mutations . In other recent studies the genetics of M. persicae was examined using microsatellite DNA genotyping analysis, but not at a global scale, mostly in populations of the same or of two neighbouring countries [6–8, 21]. In some cases certain asexual genotypes were widespread and found on herbaceous crops from year to year [8, 21].
In this paper we studied the genetic variation of M. persicae at a global scale by examining worldwide samples. A particular aim was to understand the population structure and to identify the possible routes of global dispersal.
Microsatellite analysis was carried out using four independent runs of 96 samples per run. It was possible to include some of the same samples on all four runs. The accuracy of the capillary system for alleles known to be the same was generally 1 bp or less. The ABI capillary genotyper produced slightly different size calculations for some of the samples that had been analysed previously using older ABI systems. The results were tabulated and used as input for population genetic analysis programmes.
Bayesian clustering analysis
The posterior probabilities (PPD) of the whole data set, which included all the 197 parthenogenetic lineages, were calculated with STRUCTURE software  for K values 1–15 without any prior population information. Following the pointers for choosing K provided in previous studies [22, 23], the best solution for K proved to be 3 in four independent runs. There was a sharp increase of PPD values with K moving from 1 to 3. For K > 3 the gain of information is rather less and exhibits gradually lower values. A plateau appears to be reached at K = 3 (see Additional file 1) and the information brought by the fourth K cluster (and the following) is less important than the information brought by the former three. It seems that splitting the samples in three clusters represents the optimal subdivision of the data and avoids unjustified and less informative oversplitting.
DISTRUCT software  was used to create an admixture clustering plot of the 11 M. persicae predefined samples (Figure 1). Each aphid lineage is represented as a vertical bar partitioned into three coloured segments, the lengths of which are proportional to the estimated membership coefficients of the lineage in each of the three K genetic clusters defined by STRUCTURE. The four independent runs performed using STRUCTURE gave almost identical plots (result not shown). Cluster 1 contained almost all lineages from France and those from pepper and potato from Slovenia (samples FRP and SLO in Figure 1). One French lineage from weed showed a high membership coefficient to Cluster 3. Most of the peach lineages from eastern central Greece, some from peach in northern Greece and from UK crops (other than tobacco, termed here crops), one from peach in Argentina and one from pepper in Spain had high membership coefficients for Cluster 1.
Cluster 2 is characterized by the lineages collected from tobacco (sample WET except the third lineage, sample GRT and the last two lineages in sample FES in Figure 1). Most of the lineages from peach in northern Greece, where tobacco is widely cultivated, showed a high membership coefficient for Cluster 2. Lineages belonging to Cluster 2 were also found on crops in UK, Canada, Turkey and Slovenia, on weeds in Chile and on peach in central eastern Greece and Japan.
Cluster 3 is characterized by the lineages collected from potato in New Zealand and Australia (sample AUO in Figure 1). High membership coefficient to Cluster 3 was observed in various lineages from UK crops as well as in a few lineages from crops in Slovenia, Canada and Japan and weeds in Sri Lanka. Apart from the three samples France (FRP), Greece tobacco (GRT) and New Zealand/Australia (AUO) all the others were not pure and this is also demonstrated in their mean membership coefficients to the three clusters (Table 1).
Allele frequencies, Hardy-Weinberg equilibrium and linkage disequilibrium
A moderate to high allelic diversity was found in the seven populations analysed. Apart from the highly polymorphic locus M49 (10–21 alleles) the number of alleles per locus and population ranged from 3 to 15 (> 7 in 51% of the cases) in the other five loci. The average over all loci and populations was 9.5 alleles/locus (for allele sizes see Additional file 2). Mean allelic richness (R s) ranged among populations from 7.2 to 11.8. Mean observed heterozygosity over all loci ranged from 0.623 to 0.855, with no significant differences among populations (P = 0.542). The UK population had the largest values for every diversity indicator (Table 2). Significant single locus deviation from HW equilibrium was observed in 8 out of the 56 tests. The loci which showed deviations were M35, M40, M49, M63 and myz9. In all cases the deviation was associated with a positive FIS value, i.e., a heterozygote deficiency. The multilocus test showed significant positive FIS values in four out of the seven populations examined (Table 3). Significant linkage disequilibrium was observed in 14 out of the 105 locus pairs examined.
F ST and genetic distance analysis
The pairwise multilocus FST analysis revealed important interpopulation variation and the global FST value was 0.086. Relatively high values were observed in comparisons between most of the European populations and those from Australasia. High values were also observed in most of the comparisons between the tobacco population from Europe (EUT) and those from tobacco-free regions (Europe: FRP, EUO, CGP; New Zealand/Australia: AUO). However, the peach population from a tobacco-growing region in northern Greece (NGP) showed mostly low FST values when compared with the former two sample categories, suggesting similarity to both (Table 4). The UK population showed low or moderate FST values in the comparisons with each of the other regions including Australasia suggesting it has links with all the major populations. The spatial population subdivision was also supported with the single-locus test for allelic differentiation as all pairwise comparisons were significant (Table 4). The NJ tree based on the DAS genetic distance (Figure 2) resulted in two major clusters. The first cluster contained the samples from New Zealand/Australia (AUO) and that from UK (UKO). In the second cluster tobacco lineages (sample EUT) and those from peach from the tobacco-growing region in northern Greece (sample NGP) formed a distinct group separated from the remaining samples.
The present study has revealed a broadly heterogeneous genetic structure of M. persicae at a global scale as evidenced by high allelic differentiation and relatively high FST values between certain populations, and the partitioning of genetic variation by the Bayesian clustering analysis. The observed genetic variation can be attributed to mode of reproduction, host-plant adaptation, differences between regions and dispersal.
Genetic diversity and reproductive mode
A moderate to high genic diversity was found at the intrapopulation level, as well as cases of both single and multi locus deviations from HW equilibrium which was associated with heterozygote deficiency. The populations on peach, which are expected to contain only or mostly cyclical parthenogenetic genotypes, showed HW deviations in one or both of M40 and myz9 loci. In the populations from herbaceous hosts, which probably consist of a mixture of obligate/functional and cyclical parthenogens, some heterozygote deficiency was observed in M35, M49 and M63 loci. Heterozygote deficiency in microsatellite loci appears to be common among aphid species (S. avenae, France , R. padi, France ) including M. persicae (France [7, 21], Australia  and Greece ) and it has been recorded in both sexual and asexual populations. However, some studies have found less deviation from HW equilibrium in sexual populations (M. persicae, Australia  and M. persicae, Greece ). In populations found on secondary hosts, HW deviation might be expected due to the presence of asexual lineages. Asexuality with strong clonal selection is likely to cause deviations from HW equilibrium in polymorphic loci, such as microsatellite markers, via hitch-hiking and evolution in clonal lineages. The direct effects of local clonal propagation were mediated in the current study by removing clone duplicates. In some cases, e.g., hI lineages of Rhopalosiphum padi (L.) in France  and M. persicae in Victoria Australia , a heterozygote excess has been found in asexual lineages. In R. padi this excess is attributed either to ancient loss of sexuality and the consequence of accumulated mutations or to a hybrid origin. In other cases, however, asexual lineages showed heterozygote deficiency and heterozygosity levels close to that of their sexual counterparts (hII R. padi lineages in France  and M. persicae in France ). This has been associated with a recent loss of sexuality and the time has not been sufficient to allow accumulation of mutations in asexual lineages. It has also been suggested that gene flow between sexual and asexual functional parthenogens producing males may be sufficient to prevent differences in heterozygosity accumulating between reproductive modes. Previous studies have discussed reasons why R. padi  and M. persicae  sexual populations show homozygous excess (selection, clonal expansion, Wahlund effect, inbreeding and other population effects). All of these studies concluded that null alleles were not responsible for the effects as these would have been detected during the scoring process. Wahlund effect of sampling from distinct gene pools in the same population may contribute to the homozygote excess at least in some populations that have been examined here. In support of this, the Bayesian analysis showed that some populations contain members of more than one genetic cluster.
Partitioning of genetic variation – host-plant and region
In general, the high FST values obtained in pairwise population comparisons and the estimated overall value (0.086) are among the highest reported in an aphid species using microsatellite markers [25, 26, 28] and in the same order as those reported in previous studies for M. persicae populations from Europe [7, 15] and Australia .
The Bayesian clustering and admixture analysis partitioned the genetic variation into three clusters, European (1), tobacco (2) and Australasian populations (3). Clusters 1 and 3 correspond to the generalist M. persicae persicae while Cluster 2 corresponds to the tobacco-adapted subspecies M. persicae nicotianae. Cluster members are spread over all continents and in most of the countries from which populations have been examined. These results support the hypothesis that the globalization of agriculture will have an immediate impact on the evolution of pest populations. Previous studies have provided more direct evidence of this through the spread of obligate/functional parthenogenetic genotypes (see detailed discussion below). In addition to anthropogenic activity, M. persicae populations will be influenced by natural mating and biological processes according to geographical region (Australasia vs. Southern Europe) and to tobacco adaptation, i.e., nicotianae vs. persicae (tobacco vs. other crops, tobacco or peach in tobacco regions vs. peach in non tobacco regions). In general, the proportion of membership for the tobacco Cluster 2 was greater in tobacco-growing areas. In addition, the genetic distance and the FST analyses supported the separation of the tobacco aphid populations as well as the regional population structure of persicae. It is worth noting that the separation between the peach population from eastern central Greece and the equivalent from northern Greece, as revealed by FST and Bayesian analyses (membership coefficient of CGP to European persicae: 0.55), was not as strong as observed in a previous paper (membership coefficient of the equivalent peach population to persicae cluster: 0.78, ). A possible explanation is that the two peach samples from Greece examined here (NGP and CGP) were a mixture of genotypes of both subspecies at a different ratio according to the region. An influx of the tobacco aphids (Cluster 2) into the peach orchards of eastern central Greece associated with Cluster 1 and the converse in peach orchards in northern Greece could be an explanation.
Our results suggest that certain alleles are associated with the genotypes of M. persicae feeding on tobacco. This can be considered as defining a tobacco adapted aphid 'genome' perhaps encoding a series of important enzyme variants for this specialisation as a result of continuous selection on this plant. This phenomenon is associated with the taxon M. persicae nicotianae and according to our results it is widespread and appears to have moved into countries where tobacco is not cultivated (UK and Slovenia) or the cultivation is limited (Canada). It is not surprising to find nicotianae on other crops, since it is able to colonize and reproduce on various herbaceous hosts within the vast host range of M. persicae s.l. [16, 29]. It is likely that the source of the UK nicotianae genotypes is Europe. Studies using the European suction trap network  have shown that M. persicae s.l. in Europe can migrate over southern England. Evidence that supports the continental origin of UK nicotianae is the identification in this study of a red nicotianae genotype found in the UK and also a tobacco region in southern Greece (and in Bulgaria, Fenton unpublished data). We also noticed that none of the M. persicae genotypes sampled historically in the UK (e.g., C, D, E, I, J, L in Fenton et al. [21, 31], are amongst the tobacco Cluster 2 aphids. This suggests that nicotianae has arrived rather recently in UK. The ability of the tobacco aphid to colonise new territories, even if its optimal host is not present, has interesting evolutionary implications. In addition to marker differences, it also differs physically from persicae. Generally a red nicotianae colour morph predominates in various parts of the world  and red colour populations have been associated with a complete or partial loss of sexuality [3, 32]. The red colour morph might have ecological advantages such as absorption of solar radiation  and lower choice selection by parasitoids . The red form present on tobacco plants in North America may be more resistant to organophosphorus insecticides than the green form . It has also been found that the red form of M. persicae s.l. mostly has an A1,3 autosomal translocation, which is linked to the E4-based resistance mechanism, whereas the same translocation only occasionally appears in the green form [11, 36]. Lastly, parthenogenetic lineages of the red form of M. persicae s.l. have shown better performance on tobacco plants than green ones . It has been hypothesised that adaptation to tobacco arose as a single evolutionary event in sexual populations, probably in East Asia where nicotianae was first reported as a pest [13, 38]. The tobacco-adapted population then established as permanently asexual populations in various regions. In some temperate regions the availability of peach favoured the return to a yearly sexual generation. The Bayesian clustering and admixture analyses in the present study revealed a genetic similarity of the nicotianae genotypes which strengthens the hypothesis that the adaptation to tobacco was a single evolutionary event.
The UK population and a potential link with Australasia
The UK population contained only 23 genotypes. Despite this it was the most diverse by every measurement. Approximately one third of the UK lineages belonged to each of the three Clusters. Given that peach is not openly cultivated in the UK, recent asexual populations appear to mostly develop from successive waves of colonising clones . Surprisingly, unlike the rest of the European populations, elements of the UK population were most like the Australasian population. Taking into account the admixture clustering plot, it seems that the UK is a good candidate as a source of exchange with the gene pool of the Australasian (mainly New Zealand) aphids. The earliest introductions of exogenous aphids to Australia and New Zealand were likely to be associated with settlers from Europe, especially from the UK. In support of this, previous studies revealed many common microsatellite alleles between M. persicae genotypes from Australia  and Europe . Europe does also seem to be the origin of other non-indigenous Australasian aphid pest species such as Elatobium abietinum (Walker) (Hemiptera: Aphididae) . In the present study, parameters were similar or higher when the New Zealand population were compared with European populations. Moderate or high genetic diversity has also been reported in previous New Zealand  and Australian studies [4, 6]. Theoretical and empirical work suggests a general pattern of loss of genetic diversity during colonization [42, 43]; this is because emigrant populations are serially bottlenecked [44, 45]. The substantial genetic variation observed in New Zealand and Australian M. persicae suggests that the species has not been bottlenecked and this could be attributed to sexual reproduction [41, 46] and the time that it has been there. In Australia the species was first recorded in 1910 , although it is believed to have existed there since at least 1893. This period is adequate for the mutation of new microsatellite length alleles in asexual lineages  and for sexual reproduction to give rise to diversified genotypes. In the aphid Schizaphis graminum (Rondani) (Hemiptera: Aphididae), which was introduced in the USA in the 1880s, sexual reproduction was considered as the main reason for the high diversity observed . The UK M. persicae population is believed to lack holocyclic forms and therefore restoration of the full sexual cycle through mutation or though breeding between clones with partial loss of sexual reproduction (functional parthenogenesis)  may have been required. Nevertheless, introduction of sexual clones cannot be excluded. Another factor is possible multiple introductions of M. persicae in New Zealand. Van Toor et al  reported that clones NZ2 and NZ3 appeared to be introduced as they contained many unique alleles when compared to the remaining NZ population. Additional support for the existence of multiple introductions was found when NZ3 was recognised as being a common asexual M. persicae clone found in Scotland (clone D in Fenton et al. ). Exchange of genotypes between New Zealand and Australia could also occur as demonstrated for Sitobion genotypes . The two Australian lineages examined in the current study had 17 of the 21 alleles recorded in common with the New Zealand population suggesting a recent common origin.
The present study revealed some genotypes that were sampled many miles apart in different countries some of which had been identified before, e.g., Clones B (UK and Turkey) and D (UK and New Zealand) [31, 41], and others we identified for the first time in the current study e.g., Clone M in UK and Slovenia; a genotype found in France and Greece, another in UK, Greece and Bulgaria and another in southern Greece and Slovenia. These studies have suggested that widespread clones appear to occur as a result of selection for insecticide resistance in agriculture [9, 31, 41, 52]. We have examined a relatively small number of individuals in the M. persicae population, yet we have detected these clones. This suggests that the number of successful insecticide resistant genotypes is still relatively limited, despite the possibility of resistance genes combining into more and more genotypes in sexual populations every year. In addition to the spread of resistant clones, it has also been found that asexual tobacco aphid lineages have spread between neighbouring countries such as Greece and Italy  and in the current study a widespread nicotianae lineage was found in southern Greece and Slovenia and another in Greece, UK and Bulgaria. We also report here that a distinct tobacco lineage has been found in Greece and Chile. A previous study found only one microsatellite genotype of the tobacco aphid in Chile and it seems likely that the lineage we have identified is the same as that reported by Fuentes-Contreras et al. . These studies suggest an old world origin of southern American nicotianae as the subspecies exhibits genotypic variation in Europe , but none in Chile . During the last decade several studies have revealed that the spread of certain genotypes over distant geographical areas is a common phenomenon among aphid species (e.g., Sitobion , S. avenae [25, 28, 53, 54]) including M. persicae [4, 8, 31, 41]. The rapid spread of the M. persicae s.l. lineages in different countries and continents should be attributed mostly to human transport and commerce. While winged aphids may be transported very rapidly over great distances by low level-jet streams  other studies have found that particular genotypes remain localised . The widespread lineages reported in the current work probably represent asexual genotypes reproducing parthenogenetically all-year-round. This trait enables them to spread because their reproduction will not be altered by temperature, day length or the requirement for peaches to complete their life cycle. These clones might represent 'general-purpose genotypes'  with broad ecological tolerance, which predominate in fluctuating environments through selection, although anthropogenic activities, e.g., insecticide selection pressure might also be involved [9, 52].
The present study is the first attempt to elucidate the pattern of global genetic variation of M. persicae s.l. using high resolution DNA markers. Figure 3 is a representation of some of the processes that influence population structure. The populations of M. persicae in and around peach orchards exhibit considerable genetic diversity and in some cases population parameters are close to those of sexually reproducing insects (Figure 3a). Geographical separation does create discrete gene pools as exhibited by the separation of New Zealand and European sexual populations. Sympatric speciation has generated a tobacco race of M. persicae and this process could clearly have been aided by agriculture with monocultures selecting tobacco adapted genotypes (Figure 3b). In addition to selection by tobacco, selection by insecticides is likely to have played a role in determining population structure. In the early stages of the evolution of resistant genotypes inbreeding could occur, reducing population diversity (Figure 3c). In all of these cases, any well adapted genotype could bypass sexual reproduction to become clone. The most successful of these asexual lines will become superclones spreading naturally as well as being aided by commercial activities (Figure 3). Long lived superclones will accumulate mutations. Examples of these superclones can be found in both nicotianae and persicae lineages. Evidence was reported of movement of the tobacco aphid from Europe to S. America and the spread of some persicae lineages to geographically distant regions (Europe to New Zealand). UK populations consist of asexual clones originating from diverse sources. This information highlights the ultimate ability of polyphagous aphids to adapt to different environmental conditions and the role of commerce in the globalization of insect pests.
The data set consisted of 197 aphid lineages collected from 14 countries in four continents mainly from peach and tobacco but other herbaceous hosts were included (Table 5). To obtain these unique genotypes, thousands of individuals had been sampled and analysed in the various study areas over a ten year period. The samples from peach in northern Greece and Japan were from tobacco growing regions as was the sample from weeds in Chile. All the other samples were collected in non-tobacco growing regions [see 21 (France), 41 (NZ), 10(UK)]. The samples from herbaceous crops from UK, Slovenia and southern Greece were from non peach-growing areas and they should consist mostly of asexual lineages. Most lineages were reared parthenogenetically under laboratory conditions and specimens from each lineage were kept at -80°C or in tubes filled with absolute ethanol until microsatellite analysis. Some of the samples consisted of a single aphid collected directly from the source tree or plant and stored as above.
In some cases the numbers of individuals representing an area were low and in these cases samples collected from the same region or continent and host, as well as data between years, were pooled for some of the analyses (i.e., genetic diversity indices, Hardy-Weinberg equilibrium, linkage disequilibrium, FST and genetic distance analysis). Samples from crops other than tobacco were also combined (termed here 'crops'). Both unique and multicopy genotypes were included in the analysis, but with only one copy of each multicopy genotype per population in order to avoid artificial deviations from the Hardy-Weinberg and linkage equilibria within populations and distorted estimates of allele frequencies .
DNA extraction and microsatellite genotyping
Details on DNA extraction, microsatellite loci amplification, analysis and visualization are present in a previous paper . Six microsatellite loci; M35, M40, M49, M63, M86 and myz9 , were chosen on the basis of their resolution (based on allele numbers of 12, 11, 35, 19, 21 and 18, respectively, giving 2.43 × 1013 possible combinations). Many of the lineages had been analysed in earlier work [e.g., ]. Therefore, to eliminate any doubt over allele size scoring due to technical modifications, the entire collection of genotypes was analysed again for this study using exactly the same equipment and fluorochrome primers.
Bayesian clustering analysis
A Bayesian clustering method  as implemented in the program STRUCTURE version 2.2 was used to infer the number of K unknown genetic populations in which the sampled multilocus genotypes can be split. This model-based Bayesian method also assigns a probability that the individuals belong to a certain population or to more than one population if they are admixed. In this analysis all 197 lineages were used (see Figure 1, Table 5). The data set was analysed using the admixture and uncorrelated allele frequencies models and K values 1–15 without incorporating population information. Four independent runs for each K were conducted with 100,000 iterations after a burn-in period of 20,000 iterations in each run.
Allele frequencies, Hardy-Weinberg equilibrium and linkage disequilibrium
The sample of American and Far East aphids was not included in this analysis due to the low number of the individuals examined. In addition, the samples from Scotland (SCO) and England (ENO) were combined in one sample (UKO) as well as the European samples from tobacco (pooled sample are EUT). Therefore, seven populations were analysed. Allele frequencies, mean number of alleles per locus, observed (HO) and expected (HE) heterozygosity and inbreeding coefficient (FIS) were calculated using GENEPOP version 4.0  (see also http://genepop.curtin.edu.au/). Allelic richness (R s = number of alleles independent of sample size) was also calculated using FSTAT version 22.214.171.124 , see also http://www2.unil.ch/popgen/softwares/fstat.htm). Differences in the average observed heterozygosity over all loci among populations were examined using the STRUC program of GENEPOP v. 3.4, which computes an unbiased estimate of the exact P value of a probability test of homogeneity on R × C contingency tables using a Markov chain method . Deviation from Hardy-Weinberg (HW) equilibrium at each locus was examined separately using the U test  as implemented in GENEPOP version 4.0. A Markov chain (MC) method is used for the unbiased estimation of the exact P value of this test . A multisample score test , which is performed by MC algorithm, was used as a global test across loci. Independence of microsatellite loci was examined with the G log-likelihood based exact test , which uses a simple modification of the MC algorithm described in Raymond & Rousset . The latter two tests were performed using GENEPOP version 4.0.
F ST and genetic distance analysis
Population structure was also assessed by calculating multilocus FST values  for pairwise comparisons of samples using GENEPOP version 4.0. In addition, allelic differentiation between samples was examined using GENEPOP version 4.0. The test statistic is the G log-likelihood based exact test. In these analyses seven populations were used. To further investigate the genetic relationship between populations, a neighbour joining (NJ) tree based on the allele shared distance (DAS)  was constructed using the software POPULATIONS version 1.1.28 (http://bioinformatics.org/~tryphon/populations/). DAS distance counts the number of different alleles between multilocus genotypes. Bootstrap values were calculated by resampling loci, and are presented as percentages over 1000 replications.
Blackman RL, Eastop VF: Aphids on the World's Crops. An Identification and Information Guide. 2000, London: John Wiley and Sons, 2
Blackman RL: Variation in the photoperiodic response within natural populations of Myzus persicae (Sulzer). Bull Entomol Res. 1971, 60: 533-546.
Margaritopoulos JT, Tsitsipis JA, Goudoudaki S, Blackman RL: Life cycle variation of Myzus persicae (Sulzer) (Hemiptera: Aphididae) in Greece. Bull Entomol Res. 2002, 92: 309-320.
Vorburger C, Lancaster M, Sunnucks P: Environmentally related patterns of reproductive modes in the aphid Myzus persicae and the predominance of two 'superclones' in Victoria, Australia. Mol Ecol. 2003, 12: 3493-3504. 10.1046/j.1365-294X.2003.01998.x.
Blackman RL: Life cycle variation of Myzus persicae (Sulz.) (Hom., Aphididae) in different parts of the world, in relation to genotype and environment. Bull Entomol Res. 1974, 63: 595-607.
Wilson ACC, Sunnucks P, Blackman RL, Hales DF: Microsatellite variation in cyclically parthenogenetic populations of Myzus persicae in south-eastern Australia. Heredity. 2002, 88: 258-266. 10.1038/sj.hdy.6800037.
Guillemaud T, Mieuzet L, Simon J-C: Spatial and temporal genetic variability in French populations of the peach-potato aphid, Myzus persicae. Heredity. 2003, 91: 143-152. 10.1038/sj.hdy.6800292.
Blackman RL, Malarky G, Margaritopoulos JT, Tsitsipis JA: Distribution of common genotypes of Myzus persicae (Hemiptera: Aphididae) in Greece, in relation to life cycle and host plant. Bull Entomol Res. 2007, 97: 253-263. 10.1017/S0007485307004907.
Kasprowicz L, Malloch G, Foster S, Pickup J, Zhan J, Fenton B: Clonal turnover of MACE-carrying peach-potato aphids (Myzus persicae (Sulzer), Homoptera: Aphididae) colonizing Scotland. Bull Entomol Res. 2008, 98: 115-124. 10.1017/S0007485307005445.
Kasprowicz L, Malloch G, Pickup J, Fenton B: Spatial and temporal dynamics of Myzus persicae clones in fields and suction traps. Agric For Entomol. 2008, 10: 91-100. 10.1111/j.1461-9563.2008.00365.x.
Blackman RL: Morphological discrimination of a tobacco feeding form from Myzus persicae (Sulzer) (Hemiptera: Aphididae), and a key to New World Myzus (Nectarosiphon) species. Bull Entomol Res. 1987, 77: 713-730.
Margaritopoulos JT, Tsitsipis JA, Zintzaras E, Blackman RL: Host-correlated morphological variation of Myzus persicae (Hemiptera: Aphididae) populations in Greece. Bull Entomol Res. 2000, 90: 233-244. 10.1017/S0007485300000353.
Margaritopoulos JT, Shigehara T, Takada H, Blackman RL: Host-related morphological variation within Myzus persicae group from Japan (Homoptera: Aphididae). Appl Entomol Zool. 2007, 42: 329-335. 10.1303/aez.2007.329.
Eastop VF, Blackman RL: Some new synonyms in Aphididae (Hemiptera: Sternorrhyncha). Zootaxa. 2005, 1089: 1-36.
Margaritopoulos JT, Malarky G, Tsitsipis JA, Blackman RL: Microsatellite DNA and behavioural studies provide evidence of host-mediated speciation in Myzus persicae (Hemiptera: Aphididae). Biol J Linn Soc. 2007, 91: 687-702. 10.1111/j.1095-8312.2007.00828.x.
Nikolakakis NN, Margaritopoulos JT, Tsitsipis JA: Performance of Myzus persicae (Hemiptera: Aphididae) clones on different host-plants and their host preference. Bull Entomol Res. 2003, 93: 235-242. 10.1079/BER2003230.
Margaritopoulos JT, Tsourapas C, Tzortzi M, Kanavaki OM, Tsitsipis JA: Host selection by winged colonisers within the Myzus persicae group: a contribution towards understanding host specialization. Ecol Entomol. 2005, 30: 406-418. 10.1111/j.0307-6946.2005.00700.x.
Devonshire AL, Field LM, Foster SP, Moores GD, Williamson MS, Blackman RL: The evolution of insecticide resistance in the peach-potato aphid, Myzus persicae. Philos Trans R Soc Lond, Ser B. 1998, 353: 1677-1684. 10.1098/rstb.1998.0318.
Frey JE: The analysis of athropod pest movement through trade in ornamental plants. Plant Health and the European Single Market. Brighton Crop Protection Conference Monograph no. 54: April 1993. Edited by: Ebbels DL. 1993, Farnham: British Crop Protection Council, 157-165.
Anstead JA, Williamson MS, Denholm I: Evidence for multiple origins of identical insecticide resistance mutations in the aphid Myzus persicae. Insect Biochem Mol Biol. 2005, 35: 249-256. 10.1016/j.ibmb.2004.12.004.
Fenton B, Malloch G, Navajas M, Hillier J, Birch ANE: Clonal composition of the peach-potato aphid Myzus persicae (Homoptera: Aphididae) in France and Scotland: Comparative analysis with IGS fingerprinting and microsatellite markers. Ann Appl Biol. 2003, 142: 255-267. 10.1111/j.1744-7348.2003.tb00249.x.
Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.
Garnier S, Alibert P, Audiot P, Prieur B, Rasplus J-Y: Isolation by distance and sharp discontinuities in gene frequencies: implications for the phylogeography of an alpine insect species, Carabus solieri. Mol Ecol. 2004, 13: 1883-1897. 10.1111/j.1365-294X.2004.02212.x.
Rosenberg NA: DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004, 4: 137-138. 10.1046/j.1471-8286.2003.00566.x.
Simon J-C, Baumann S, Sunnucks P, Hebert PDN, Pierre J-S, Le Gallic J-F, Dedryver C-A: Reproductive mode and population genetic structure of the cereal aphid Sitobion avenae studied using phenotypic and microsatellite markers. Mol Ecol. 1999, 8: 531-545. 10.1046/j.1365-294x.1999.00583.x.
Delmotte F, Leterme N, Gauthier J-P, Rispe C, Simon J-C: Genetic architecture of sexual and asexual populations of the aphid Rhopalosiphum padi based on allozyme and microsatellite markers. Mol Ecol. 2002, 11: 711-723. 10.1046/j.1365-294X.2002.01478.x.
Delmotte F, Leterme N, Bonhomme J, Rispe C, Simon J-C: Multiple routes to asexuality in an aphid species. Proc R Soc Lond B. 2001, 268: 2291-2299. 10.1098/rspb.2001.1778.
Llewellyn KS, Loxdale HD, Harrington R, Brookes CP, Clark SJ, Sunnucks P: Migration and genetic structure of the grain aphid (Sitobion avenae) in Britain related to climate and clonal fluctuation as revealed using microsatellites. Mol Ecol. 2003, 12: 21-34. 10.1046/j.1365-294X.2003.01703.x.
Semtner PJ, Tilson WM, Dara SK: Performance of the tobacco aphid (Homoptera: Aphididae) on various host plants. J Entomol Sci. 1998, 33: 180-195.
Taylor LR: Synoptic dynamics, migration and the Rothamsted insect survey: Presidential address to the British Ecological Society, December 1984. J Anim Ecol. 1986, 55: 1-38. 10.2307/4690.
Fenton B, Malloch G, Woodford JAT, Foster SP, Anstead J, Denholm I, King L, Pickup J: The attack of the clones: tracking the movement of insecticide-resistant peach-potato aphids Myzus persicae (Hemiptera: Aphididae). Bull Entomol Res. 2005, 95: 483-494. 10.1079/BER2005380.
Fuentes-Contreras E, Figueroa CC, Reyes M, Briones LM, Niemeyer HM: Genetic diversity and insecticide resistance of Myzus persicae (Hemiptera: Aphididae) populations from tobacco in Chile: evidence for the existence of a single predominant clone. Bull Entomol Res. 2004, 94: 11-18. 10.1079/BER2003275.
Markkula M, Rautapää J: The effect of light and temperature on the colour of the English grain aphid, Macrosiphum avenae (F.) (Hom., Aphididae). Ann Entomol Fenn. 1967, 33: 1-13.
Libbrecht R, Gwynn DM, Fellowes MDE: Aphidius ervi p referentially attacks the green morph of the pea aphid, Acyrthosiphon pisum. J Insect Behav. 2007, 20: 25-32. 10.1007/s10905-006-9055-y.
Harlow CD, Southern PS, Lambert ER: Geographic distribution of two colour forms, carboxylesterase activity, and chromosome configuration of the tobacco aphid (Homoptera: Aphididae) in North Carolina. Ann Entomol Soc Am. 1991, 84: 1175-1179.
Blackman RL, Spence JM, Field LM, Devonshire AL: Variation in the chromosomal distribution of amplified esterase (FE4) genes in Greek field populations of Myzus persicae (Sulzer). Heredity. 1999, 82: 180-186. 10.1038/sj.hdy.6884580.
Goundoudaki S, Tsitsipis JA, Margaritopoulos JT, Zarpas KD, Divanidis S: Performance of the tobacco aphid Myzus persicae (Hemiptera: Aphididae) on Oriental and Virginia tobacco varieties. Agric ForEntom. 2003, 5: 285-291. 10.1046/j.1461-9563.2003.00190.x.
De Jong JK: Enkele resultaten von het onderzoek naar de biologie von de tabaksluis Myzus persicae Sulzer. Bulletin van het Deli Proefstation te Medan. 1929, 28: 1-36.
Terradot L, Simon J-C, Leterme N, Bourdin D, Wilson ACC, Gauthier J-P, Robert Y: Molecular characterization of clones of the Myzus persicae complex (Hemiptera: Aphididae) differing in their ability to transmit the potato leafroll luteovirus (PLRV). Bull Entomol Res. 1999, 89: 355-363. 10.1017/S0007485399000498.
Nicol D, Armstrong KF, Wratten SD, Walsh PJ, Straw NA, Cameron CM, Lahmann C, Frampton CM: Genetic diversity of an introduced pest, the green spruce aphid Elatobium abietinum (Hemiptera: Aphididae) in New Zealand and the United Kingdom. Bull Entomol Res. 1998, 88: 537-543.
van Toor RF, Foster SP, Anstead JA, Mitchinson S, Fenton B, Kasprowicz L: Insecticide resistance and genetic composition of Myzus persicae (Hemiptera: Aphididae) on field potatoes in New Zealand. Crop Prot. 2008, 27: 236-247. 10.1016/j.cropro.2007.05.015.
Austerlitz F, Jung-Muller B, Godelle B, Gouyon P-H: Evolution of coalescence times, genetic diversity and structure during colonization. Theor Popul Biol. 1997, 51: 148-164. 10.1006/tpbi.1997.1302.
Hewitt GM: Post-glacial re-colonization of European biota. Biol J Linn Soc. 1999, 68: 87-112. 10.1111/j.1095-8312.1999.tb01160.x.
Ibrahim KM, Nichols RA, Hewitt GM: Spatial patterns of genetic variation generated by different forms of dispersal during range expansion. Heredity. 1996, 77: 282-291. 10.1038/hdy.1996.142.
Fry AJ, Zink RM: Geographic analysis of nucleotide diversity and song sparrow (Aves: Emberizidae) population history. Mol Ecol. 1998, 7: 1303-1313. 10.1046/j.1365-294x.1998.00462.x.
Cottier W: Aphids of New Zealand. N. Z. Department of Scientific and Industrial Research Bulletin. 1953, 106: 210-211.
Zeck EH: The green peach aphid, Myzus persicae (Sulzer). Agric Gaz NSW. 1928, 39: 147-154.
Wilson ACC, Sunnucks P, Hales DF: Heritable genetic variation and potential for adaptive evolution in asexual aphids (Aphidoidea). Biol J Linn Soc. 2003, 79: 115-135. 10.1046/j.1095-8312.2003.00176.x.
Shufran KA, Peters DC, Webster JA: Generation of clonal diversity by sexual reproduction in the greenbug, Schizaphis graminum. Insect Mol Biol. 1997, 6: 203-209. 10.1046/j.1365-2583.1997.00174.x.
Blackman RL: The inheritance of life-cycle differences in Myzus persicae (Sulz.) (Hem., Aphididae). Bull Entomol Res. 1972, 62: 281-294.
Wilson ACC, Sunnucks P, Hales DF: Microevolution, low clonal diversity and genetic affinities of parthenogenetic Sitobion aphids in New Zealand. Mol Ecol. 1999, 8: 1655-1666. 10.1046/j.1365-294x.1999.00751.x.
Zamoum T, Simon J-C, Crochard D, Ballanger Y, Lapchin L, Vanlerberghe-Masutti F, Guillemaud T: Does insecticide resistance alone account for the low genetic variability of asexually reproducing populations of the peach-potato aphid Myzus persicae?. Heredity. 2005, 94: 630-639. 10.1038/sj.hdy.6800673.
Sunnucks P, England PR, Taylor AC, Hales DF: Microsatellite and chromosome evolution of parthenogenetic Sitobion aphids in Australia. Genetics. 1996, 144: 747-756.
Haack L, Simon JC, Gauthier JP, Plantegenest M, Dedryver CA: Evidence for predominant clones in a cyclically parthenogenetic organism provided by combined demographic and genetic analyses. Mol Ecol. 2000, 9: 2055-2066. 10.1046/j.1365-294X.2000.01108.x.
Elton CS: The dispersal of insects to Spitsbergen. Trans Entomol Soc Lond. 1925, 73: 289-299.
Lynch M: Destabilizing hybridization, general-purpose genotypes and geographic parthenogenesis. Q Rev Biol. 1984, 59: 257-290. 10.1086/413902.
Sunnucks P, De Barro PJ, Lushai G, MacLean N, Hales D: Genetic structure of an aphid studied using microsatellites: cyclic parthenogenesis, differentiated lineages and host specialization. Mol Ecol. 1997, 6: 1059-1073. 10.1046/j.1365-294X.1997.00280.x.
Malloch G, Highet F, Kasprowicz L, Pickup J, Neilson R, Fenton B: Microsatellite marker analysis of peach-potato aphids (Myzus persicae, Homoptera: Aphididae) from Scottish suction traps. Bull Entomol Res. 2006, 96: 573-582. 10.1017/BER2006459.
Wilson ACC, Massonnet B, Simon J-C, Prunier-Leterme N, Dolatti L, Llewellyn KS, Figueroa CC, Ramirez CC, Blackman RL, Estoup A, Sunnucks P: Cross-species amplification of microsatellite loci in aphids: assessment and application. Mol Ecol Notes. 2004, 4: 104-109. 10.1046/j.1471-8286.2004.00584.x.
Rousset R: GENEPOP'007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.
Goudet J: FSTAT (version 1.2): a computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.
Raymond M, Rousset F: An exact test for population differentiation. Evolution. 1995, 49: 1280-1283. 10.2307/2410454.
Guo SW, Thompson EA: Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992, 48: 361-372. 10.2307/2532296.
Goudet J, Raymond M, de Meeüs T, Rousset F: Testing differentiation in diploid populations. Genetics. 1996, 144: 1933-1940.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Chakraborty R, Jin L: A unified approach to study hypervariable polymorphisms: statistical considerations of determining relatedness and population distances. DNA Fingerprinting: State of the Science. Edited by: Pena SDJ, Chakraborty R, Epplen JT, Jeffreys AJ. 1993, Basel: Birkhäuser, 153-175.
The authors are grateful to Drs. Claudio Ramìrez, Eduardo Fuentes-Contreras, Tomoko Sigehara, Alex Wilson, Steve Foster, Jon Pickup, Fiona Highet, Mark Stephens, Ron van Toor, Alain Migeon, Maria Navajas and Panagiotis Skouras for help providing aphid material. Philip Smith and Jiasui Zhan made constructive comments on the manuscript. The present study was partially supported by the PENED 2003 project. The project is co-funded from the European Union (European Social Fund) by 80%, from the Hellenic State (Ministry of Development, General Secretariat for Research and Technology) by 20% (public grants) and from Hellenic Private Sector, in the framework of Measure 8.3 of Operational Programme "Competitiveness" (the 3rd Community Support Programme). SCRI is supported by the Scottish Government Rural and Environment Research and Analysis Directorate.
Microsatellite data were obtained and scored by LK, GM and BF. Statistical analyses were performed by JTM. The study was designed by all authors. The paper was written by JTM and BF. The manuscript has been read and approved by all authors.
John T Margaritopoulos, Louise Kasprowicz, Gaynor L Malloch and Brian Fenton contributed equally to this work.