Skip to main content

Agroecosystems shape population genetic structure of the greenhouse whitefly in Northern and Southern Europe



To predict further invasions of pests it is important to understand what factors contribute to the genetic structure of their populations. Cosmopolitan pest species are ideal for studying how different agroecosystems affect population genetic structure within a species at different climatic extremes. We undertook the first population genetic study of the greenhouse whitefly (Trialeurodes vaporariorum), a cosmopolitan invasive herbivore, and examined the genetic structure of this species in Northern and Southern Europe. In Finland, cold temperatures limit whiteflies to greenhouses and prevent them from overwintering in nature, and in Greece, milder temperatures allow whiteflies to inhabit both fields and greenhouses year round, providing a greater potential for connectivity among populations. Using nine microsatellite markers, we genotyped 1274 T. vaporariorum females collected from 18 greenhouses in Finland and eight greenhouses as well as eight fields in Greece.


Populations from Finland were less diverse than those from Greece, suggesting that Greek populations are larger and subjected to fewer bottlenecks. Moreover, there was significant population genetic structure in both countries that was explained by different factors. Habitat (field vs. greenhouse) together with longitude explained genetic structure in Greece, whereas in Finland, genetic structure was explained by host plant species. Furthermore, there was no temporal genetic structure among populations in Finland, suggesting that year-round populations are able to persist in greenhouses.


Taken together our results show that greenhouse agroecosystems can limit gene flow among populations in both climate zones. Fragmented populations in greenhouses could allow for efficient pest management. However, pest persistence in both climate zones, coupled with increasing opportunities for naturalization in temperate latitudes due to climate change, highlight challenges for the management of cosmopolitan pests in Northern and Southern Europe.


The dispersal of phytophagous insect pests can be enhanced by worldwide trade and human movement [1],[2]. In addition, climate change facilitates movement of various taxa polewards [3],[4]. Following introduction to new habitats, the establishment of insect pest populations can be favored by benign climates, as well as by monocultures in agroecosystems, i.e. agricultural fields and greenhouses [5]-[7]. Low genetic diversity of insect pest populations in newly occupied habitats suggests that even a single successful founder event is enough to establish populations [8],[9]. However, further spread of introduced pests into natural ecosystems depends on the environment surrounding the initial introduction and on the origin of the introduced species [10],[11]. For example, pests of tropical origin may be more likely to establish themselves in the Mediterranean than in the boreal climate zone [12],[13]. At southern latitudes there are suitable climatic conditions and year-round availability of host plants, in both natural ecosystems and densely aggregated agroecosystems. In contrast, at northern latitudes, natural habitats are only seasonally available and greenhouses are often sparsely distributed. Thus, the extent of establishment and spread of invasive pests in the North might be more dependent on the distribution of agroecosystems, particularly greenhouses, than it is in the South.

Enclosed greenhouse environments are designed to reduce evaporation, pest entry [14] and loss of expensive biological pest control agents [15], to ensure efficient crop maintenance. Because greenhouses are relatively closed environments, pest populations in greenhouses might be generally more affected by insecticide applications and host plant changes than are pest populations in fields. These crop management practices can lead to reductions in population size and selection for resistant genotypes in the pests leading to increased homozygosity within and differentiation between pest populations [16]. Thus, populations of insects inhabiting greenhouses might show more genetic differentiation than those in fields. Indeed, populations of phytophagous pests inhabiting greenhouses often show population genetic structure, e.g. Tetranychus urticae Koch [17],[18], although dispersal and gene flow can also be restricted among pest populations inhabiting fields, e.g. Leptinotarsa decemlineata Say [8].

The greenhouse whitefly (Trialeurodes vaporariorum Westwood) is an invasive pest which was brought to Europe (UK) on Orchidaceae from Mexico in 1856 [19]. Soon after introduction, it spread to the European continent, and in 1920 it was recorded in greenhouses in Finland [20]. In Finland, it spread through transportation on plant seedlings above the Arctic Circle to Rovaniemi, bringing considerable damage to tomato and cucumber crops, as well as to ornamental plants [21],[22]. T. vaporariorum was reported from the Mediterranean region only later, in 1963 [19]. The species was noticed in Greece (Crete) only in 1978, when it began to cause pest management problems due to its resistance to insecticides [23],[24].

T. vaporariorum currently has an almost cosmopolitan distribution [25],[26]. Its success can be attributed to the worldwide distribution of greenhouse habitats, polyphagy [27], its tolerance of higher or lower temperatures than its biological control agents [17],[26], and its haplodiploid mode of reproduction [28]. However, the absence of an overwintering resting stage [29] potentially limits its spread to natural ecosystems. Since the development of T. vaporariorum ceases at 8.3°C [30], year-round populations might persist at southern latitudes, where some host plants are available during winter, but are not likely to persist at northern latitudes, where crop cultivation in fields is seasonal and wild host plants decay during winter.

To date, population genetic structure has been analyzed in only a few whitefly species, and little is known about population genetic structure in T. vaporariorum. The related Bemisia tabaci species complex, particularly Mediterranean B. tabaci (Med), is characterized by high genetic diversity and differentiation of populations, as indicated by both mitochondrial and microsatellite markers [31],[32] (except in recently introduced populations in Taiwan and France [33],[34]). Populations of B. tabaci (Med) in Greece separated by just a few kilometers show population genetic structure, possibly due to separate founder events or an older population history in this country [35]. Unlike the B. tabaci species complex, T. vaporariorum populations have low genetic diversity in mitochondrial genes [36],[37]. Recent findings indicate that sequences of three mitochondrial genes and composition of endosymbiont communities from populations sampled from different continents show little variation (Kapantaidaki et al., unpubl.). Analysis of a few nuclear genes (allozymes) in T. vaporariorum populations from greenhouses in South Korea revealed their subdivision possibly due to restricted gene flow by natural geographic barriers [38]. However, studies of population genetic structure in T. vaporariorum with other, more polymorphic genetic markers, allowing description of more recent evolutionary processes, have not been performed until now. Recent findings of variation in phenotypic responses, particularly diverse responses to insecticide treatments among geographically close populations, suggest differentiation and low gene flow among invaded greenhouses [39].

To understand how insect pests respond to different environmental conditions, such as climate, habitat, and crop management practices, and the role of agroecosystems in shaping population genetic structure, comparative studies of population genetic diversity of pests in different climate zones are necessary. In this study we present the first extensive genetic data on population structure of the greenhouse whitefly. We compare the genetic structure of T. vaporariorum populations in Finland and in Greece, representing boreal and Mediterranean climate zones, and evaluate the influence of host plants and agricultural practices on the spatial and temporal population genetic structure of this invasive species. We hypothesize that T. vaporariorum populations in Northern Europe are more likely to be genetically differentiated than populations in Southern Europe, because this species is expected to be restricted to greenhouses in the North.



In Finland we sampled commercial greenhouses that operate year-round and produce primarily tomato and cucumber crops of various cultivars (Table 1). Samples were collected from greenhouses belonging to different growers. In total 18 greenhouses were sampled in spring of 2010–2012. Ten of these were sampled twice: in 2010 and in 2011. Sampling was concentrated in Ostrobothnia (16 greenhouses) but also included two distant locations in other parts of the country (Figure 1-I). Ostrobothnia was the focus of our study because in this area we could find multiple greenhouses with different management practices in terms of host plant species, their cultivars (Table 1) and the origin of seedlings. Two to five greenhouses belonging to different growers were sampled within Närpes, Töjby and Pjelax villages. The minimum and maximum distances between these villages were 9 and 32 km, respectively, measured as straight line distance between coordinates. The distances between greenhouses within villages ranged from 1.1 to 3.7 km in Närpes, 0.4 km in Töjby and from 0.28 to 0.9 km in Pjelax.

Table 1 Description of the samples collected in Finland and Greece
Figure 1
figure 1

Maps of sampling locations. Sample codes are listed in Table 1. I - Finland, II - Greece.

In Greece we sampled 16 agroecosystems in different seasons over several years: 2004–2011. These included eight greenhouses and eight fields growing various crop plants (Table 1) which were distributed throughout mainland Greece, the Peloponnese and the Island of Crete (Figure 1-II). Sampling was concentrated in the fields of West Peloponnese because open environments in this region cover the expected range (7–20 km) of the potential dispersal abilities of whiteflies (natural or by wind, as known from the B. tabaci species complex; [40],[41]). In West Peloponnese, the minimum distance between samples ranged from 3.4 km (between WP3 and WP4) to 24.4 km (between WP3 and WP5), and the maximum distance between locations reached 100 km.

During sampling, both genders of whiteflies were collected using a mouth aspirator. The whiteflies were preserved in 90% Ethanol and stored at 4°C until they were sexed and used for genotyping. Since T. vaporariorum is a haplodiploid species, only adult females were chosen for genotyping.

DNA extraction and microsatellite genotyping

Total genomic DNA was extracted from each individual female as described in Tsagkarakou et al. [42]. Nine microsatellite markers (Table 2) out of the 13 characterized in Molecular Ecology Resources Primer Development Consortium et al. [43] were used to genotype 1274 T. vaporariorum females, 800 from Finland and 474 from Greece (Table 3). Four of the microsatellite markers described previously did not amplify consistently, and thus, were excluded from this study. Three multiplex amplification reactions were performed as described in Molecular Ecology Resources Primer Development Consortium et al. [43] with slight modifications. Diluted amplification products were separated on an ABI 3130xl genetic analyzer (Applied Biosystems). Allele sizes were scored against GeneScan™ 500 LIZ standard using GeneMapper® v 4.0 software (both Applied Biosystems) and were confirmed manually.

Table 2 Characteristics of the nine polymorphic microsatellite loci analysed in T. vaporariorum
Table 3 Genetic diversity estimated over the nine microsatellite loci for samples of T. vaporariorum

Data analysis

To analyze genetic distance between samples, both between and within countries, pairwise estimates of F ST were calculated in Arlequin 3.11 [44]. The significance of the genetic distances at the 0.05 level was tested by permuting the individuals or genotypes between the samples 110 times and adjusting P values with strict Bonferroni correction. Since all pairwise F ST between samples from Finland and Greece showed statistically significant differences, and due to the low probability of gene flow between the two distant countries, all further analyses were done for each country separately.

The samples taken in two consecutive years from the same greenhouse in Finland showed no genetic differentiation, except for one sampling location (TJ-2 a and TJ-2 b). Therefore, in most analyses we used a combined dataset, which pooled samples collected from the same greenhouse in 2010 and 2011, (except for TJ-2 a and TJ-2 b, which were considered as separate samples). For other analyses (specified below) we used a separated dataset, in which each sampling effort in Finland was considered a separate sample. Each sampling effort in Greece was considered a separate sample in all analyses, because none of the locations were sampled in consecutive years.

For each sample, mean observed (H O ) and expected heterozygosity (H E ), and mean number of alleles (N A ) per locus were calculated using GenAlEx v. 6.5 [45]. Observed and expected heterozygosities were also calculated for each locus over the total data from each country. Departure from Hardy-Weinberg expectations (HWE) was tested with 1000 permutations using a global test across loci or samples as implemented in GENEPOP v. 4.2 [46]. The test was performed using Fisher’s method, testing hypotheses of heterozygote deficiency and heterozygote excess [47], and producing global p value estimates for each sample over all loci and for each locus over samples from Finland and Greece. Genotypic linkage disequilibrium for each pair of loci in the samples was tested using the log likelihood ratio statistic (G-test) as implemented in GENEPOP v. 4.2. For multiple tests, statistical significance was adjusted using strict Bonferroni corrections [48]. The samples were analyzed for potential scoring errors in all loci using MICRO-CHECKER v. 2.2.3 and the frequency of null alleles (f) was estimated [49].

To investigate the relationship between genetic and geographic distance, isolation by distance was analyzed in GenAlEx v. 6.5 [45]. Genetic distance was defined by pairwise linear F ST, (F ST /(1- F ST )), and geographic distance was defined as pairwise distances generated from geographical coordinates expressed in decimal degrees. The correlation between the two data matrices was assessed using a Mantel test and its significance estimated by P values, the regression coefficient (R2), and the mean correlation coefficient (R XY ) over 999 random permutations of linear F ST values as implemented in GenAlex v.6.5. Isolation by distance was assessed with smaller subsets of the data as well as using the full datasets, to evaluate the influence of scale on the relationships.

Analyses of molecular variance (AMOVA) were performed using Arlequin 3.11 [44] to estimate and compare the percentage of genetic variation explained by different hierarchical groups (i.e. individual, sample, group of samples). Four analyses were constructed to test the following groups: 1) country (Finland vs. Greece), 2) host plant species (cucumber vs. tomato) in samples from Finland only, 3) host plant botanical family (Cucurbitaceae vs. Solanaceae vs. Fabaceae vs. Rosaceae) in samples from Greece only, and 4) habitat (greenhouse vs. field) in samples from Greece only. For analysis 2, samples HR, ML and NL were excluded since the whiteflies in these greenhouses might have been exposed to both cucumber and tomato grown in the same compartment or greenhouse. However, in analysis 3, sample MA 2 was not excluded because several hosts grown in the same compartment belonged to the same family (Solanaceae) (see Table 1).

To assess the level of genetic differentiation between groups defined above for AMOVA, we compared summary statistics calculated for the different groups: F IS (inbreeding coefficient measuring heterozygote deficit within populations), F ST (a measure of population structure and heterozygote deficit among populations), allelic richness (measure of the number of alleles independent of sample size), HE (unbiased expected heterozygosity) and HO (observed heterozygosity). FSTAT v. 2.9.3 [50] was used to calculate the average (over samples and loci; weighted by sample size) of the chosen statistics for each group and for their comparison. Statistical significance was assessed after 1000 permutations. As in the AMOVA groups, some samples were excluded because multiple hosts were grown in the same compartment or greenhouse (see above).

The relationship between environmental variables and genetic structure of the studied populations was estimated using default settings of the software GESTE v. 2.0 [51]. The software gives the highest posterior probability (Pr) to the model explaining genetic structure the best, evaluating environmental variables separately and in combination through a generalized linear model. In this analysis, the separated data set for the samples from Finland was used. Latitude, host plant species, cultivar, crop source and year of sampling were evaluated as explanatory variables of population structure in Finland. Latitude, longitude, four host plant families and habitat (field or greenhouse) were evaluated as explanatory variables of population structure in Greece.

Bayesian clustering analysis implemented in STRUCTURE v.2.3.4 [52] was used to infer the number of genetically distinct clusters (K) in each country using a model of no admixture, correlated allele frequencies and including the sampling location as a prior [53]. Initial analyses were performed both with admixture and no admixture models, but the later was selected since visualization of the results was more straightforward and no differences in the most likely number of clusters were observed for the two models. Analysis parameters included a burn-in period of 250,000 followed by 500,000 MCMC iterations. For each dataset, Finland and Greece, we tested K from 2 to 10, with ten replicate analyses per value of K. Subsets of each dataset were analyzed with the same settings. The most likely number of clusters in our samples was determined using the ΔK approach [54] as implemented in Structure Harvester v. 0.56.3 [55]. Results were visualized as bar plots by finding the optimal alignment of the ten replicate analyses of the “best” K in CLUMPP v. 1.1.2 [56] using the Greedy algorithm and 1000 random input orders, and then by creating graphics in Distruct v. 1.1 [57]. For Finland, the combined dataset was used first, and then two subsets of the data were created. In these subsets, the samples were grouped by host plant species and samples HR a, b and ML a, b were separated since they had been collected from different hosts (see Table 1). Samples HR, ML and NL were included in both data subsets since these samples might have been exposed to several hosts grown in the same compartment or greenhouse. For Greece, all samples were first analyzed together, then data subsets were created grouping samples by habitat (field or greenhouse). The definition of the data subsets (by host plant species or habitat) was chosen after considering the results of initial analyses with the full data sets and our analysis with GESTE v. 2.0 [51].


Genetic diversity of microsatellite loci and samples

Significant deviations from HWE through heterozygote deficiencies were detected at one locus (Tvap 3–1) and 4 loci (Tvap-4-2, Tvap-1-4, Tvap-1-5 and Tvap-3-1) in the Finnish and Greek samples, respectively (Table 2). At the sample level, a test of HWE across the nine microsatellite loci indicated significant heterozygote deficiency in two Finnish and three Greek samples (Table 3). There were no cases of significant heterozygote excess.

Three loci showed a null allele frequency > 0.2: Tvap-1-1C (2 samples), Tvap-3-1 (12 samples) and Tvap-1-4 (1 sample). For each of these loci, the frequency of null alleles within a sample varied: f = 0.110-0.258 (Tvap-1-1C), f = 0.164-0.401 (Tvap-3-1), and f = 0.166-0.238 (Tvap-1-4), and the average frequency of null alleles over the samples ranged from 0.119 to 0.250. No cases of large allele drop out were found. Even though null alleles are present, deviations from HWE could be also due to significant homozygosity in populations inhabiting the human-mediated environment (i.e. due to population bottlenecks and inbreeding), rather than due to significant genotyping errors.

Genotypic linkage disequilibrium tested for each pair of loci for each sample revealed a potential association between loci Tvap-1-1 and Tvap 3–1 in sample PJ 4. Since locus Tvap-3-1 was characterized by homozygote excess and had a high frequency of null alleles only in sample PJ 4 (f = 0.401), we suspect that the linkage disequilibrium indicated for this sample does not reflect a true association between the loci. Therefore, data from all nine loci were used in the analyses.

Differences between Finland and Greece

AMOVA indicated significant genetic structure between the two geographic areas (Table 4). The percentage of variation explained by country of origin, Finland vs. Greece, was higher than that among the samples within each country (9.90% and 6.87%, respectively; Table 4), indicating that overall genetic variation might be explained by these groups. T. vaporariorum from the two countries also differed significantly in their global observed and expected heterozygosities (Finland: H O /H E  = 0.350/0.385 vs. Greece: H O /H E  = 0.451/0.496), and in allelic richness (2.498 vs. 3.234 for Finland vs. Greece, respectively) (all P = 0.001). However, F statistics calculated for each country did not differ statistically (Finland/Greece; F IS : 0.091/0.090, P = 0.157; F ST : 0.093/0.055, P = 0.976). Nevertheless, the range of pairwise F ST values between samples within countries was broader for Finland (−0.006 < F ST  < 0.533) than it was for Greece (−0.007 < F ST  < 0.164) (Tables 5A and B).

Table 4 Distribution of the molecular variance between and within four groups of samples of T. vaporariorum
Table 5 Pairwise estimates of F ST between samples in Finland (A) and Greece (B) over the nine microsatellite loci

Population structure in Finland

Seventy nine percent of pairwise F ST comparisons (121 of 154) between samples from Finland showed significant population differentiation (Table 5A). Some populations (HR, KR, PR, TJ 2b and LH) were differentiated from all other samples (Table 5A). However, one of the samples most distant from the Ostrobothnia region (NL) was not significantly different in pairwise F ST from one of the Ostrobothnian samples (OV). For samples collected from different greenhouses at the same location (NR, PJ and TJ), there was no significant genetic structure, except for TJ: TJ1 was not differentiated from TJ 2a, but both of these samples were differentiated from TJ 2b. Moreover, some of these samples were also not differentiated from other samples in a neighboring village (Table 5A).

There was no evidence for isolation by distance when all Finnish samples were analyzed together (R XY  = 0.259, R2 = 0.067, P = 0.190), nor when samples outside Ostrobothnia were excluded from the dataset (R XY  = 0.138, R2 = 0.019, P = 0.211). However, it was moderately strong and significant (R XY  = 0.480, R2 = 0.231, P = 0.001), when only the core samples in Ostrobothnia were analyzed (outlying samples HR, ML and KR). The data was explained best when only NR, PJ and TJ samples were analyzed (R XY  = 0.618, R2 = 0.382, P = 0.001) and at the smallest geographic scale, when including only NR and PJ samples (R XY  = 0.876, R2 = 0.767, P = 0.007).

In AMOVA, the percentage of variation explained by host plant species (cucumber vs. tomato) was low (1.69%) but statistically significant (P = 0.036). The percentage of variation explained by the groups was lower than that among samples within each group (8.59%, Table 4). However, T. vaporariorum collected from cucumber and tomato hosts differed significantly in allelic richness (P = 0.005, 2.25 vs.2.61), heterozygosities (P = 0.002, H O /H E : 0.278/0.300 vs. 0.371/0.413), and genetic differentiation (P = 0.002, F ST : 0.228 vs.0.043) (cucumber vs. tomato respectively). Overall, samples collected from cucumber exhibited less genetic diversity and a higher degree of genetic differentiation than samples collected from tomato.

Furthermore, host plant species was indicated to be the most important factor explaining population genetic structure according to our analysis with GESTE. The model with highest posterior probability (Pr = 0.81) was the model with host plant species. All other models including other environmental factors or their combination had very low probability values. For example, Pr = 0.0446 was determined for the model including crop origin in combination with host plant species and Pr = 0.0443 was determined for the model including crop origin in combination with cultivar, to mention a few.

The Bayesian analysis of population structure indicated that the 18 samples from Finland (Table 1) represent three main genetic clusters (K = 3; Figure 2A). The genetic clusters could be characterized by both host plant species and geographic location. Most samples were clearly assigned to one of the three clusters (assignment higher than 80%). However, for samples PR and NR3, the assignment was mixed between cluster 1 and 2 (PR: K 1  = 0.453 and K 2  = 0.517; NR3: K 1  = 0.769 and K 2  = 0.195). Majority of samples in K 1 were collected from tomato, with two exceptions: LH and NL were collected from cucumber and these samples were geographically distant from the other samples in this cluster. K 2 included four samples collected from cucumber (NR 2, TJ 1, TJ 2a and TJ 2b). Notably, NR 2 was assigned to a different cluster than the samples from other greenhouses in Närpes. K 3 corresponded to a single greenhouse growing cucumber (KR), which, being surrounded by forests was somewhat geographically isolated from other greenhouses in Ostrobothnia.

Figure 2
figure 2

Genetic clusters of T. vaporariorum sampled in Finland delineated by Bayesian analysis implemented in STRUCTURE [[52]] . In each graph individuals are represented by vertical bars broken into K differently colored genetic clusters, with length proportional to assignment in each cluster. Probability of assignment of each individual is indicated from 0 to 1. Individuals are grouped by sample with the sample code listed at the bottom of the bar graph. A: Analysis of all 18 samples, 800 individuals, K = 3. B: Analysis of a subset of 14 samples collected from tomato, 570 individuals, K = 4. C: Analysis of a subset of 11 samples collected from cucumber, 380 individuals, K = 3. The combined dataset was used in the initial run (A). Samples HR, ML and NL were included in analysis of both host plant groups (B and C) since these samples might have been exposed to several hosts grown in the same compartment or greenhouse. Samples HR a, b and ML a, b were separated in analyses B and C as they were exposed to different hosts.

To resolve sub-clustering, further Bayesian analyses with STRUCTURE were conducted for samples collected from tomato and cucumber hosts separately (Figure 2B and C), since host plant species was the major component of genetic structure revealed in our initial analysis and our analysis with GESTE. Some samples were included in both subsequent runs (see Methods). The subset of samples from tomato (K 1 in the initial analysis) was characterized best by four sub-clusters (K = 4) (Figure 2B). Whiteflies collected from PJ formed sub-cluster 1 (K 1.1 ). NR 1 and NR 3 were not clearly resolved and partially formed sub-cluster 2 (K 1.2 ) with NL (NR 1a, b: K 1.2  = 0.584 and K 1.1  = 0.287; NR 3a, b: K 1.2  = 0.542 and K 1.3  = 0.372). PR and OV formed sub-cluster 3 (K 1.3 ), whereas HR and ML formed sub-cluster 4 (K 1.4 ). Assignment was not resolved well for ML b: K 1.4  = 0.543 and K 1.3  = 0.297. Additional analysis of an even smaller subset of samples (PJ 1-PJ 5) collected from greenhouses that were all growing tomato and located within the natural dispersal range expected for the whitefly, revealed complete genetic homogeneity (data not shown) and indicated high gene flow at this location. The subset of samples collected from cucumber (K 2 and K 3 in the initial analysis) was characterized best by three well-resolved sub-clusters (K = 3), all with 85-100% assignment (Figure 2C). HR, ML and NL formed sub-cluster 1 (K 2.1 ), whereas LH, NR 2, TJ 1, TJ 2a and TJ 2b formed sub-cluster 2 (K 2.2 ). KR alone formed sub-cluster 3 (K 3.1 ).

Population structure in Greece

Seventy five per cent of pairwise F ST comparisons (91 out of 120) between samples from Greece showed significant population differentiation (Table 5B). Samples AT, CR 1, CR 3 and MA 2 differed significantly from all other samples (Table 5B). Significant differentiation among samples was absent only in some geographically close sites, in particular among the samples collected from fields in West Peloponnese (except WP 2 vs. WP 7). Some samples collected from greenhouses in West Greece and Peloponnese were not differentiated from the samples collected from fields: WG was not differentiated from WP 5–7 and EP, and WP 6 was not differentiated from EP. One sample from Crete (CR 2) was similar to some samples from West Peloponnese (WP 6 and WP 7), as well as the northern samples (MA 1 and NP), which were similar to each other as well (Table 5B).

There was evidence for isolation by distance in the Greek dataset (R XY  = 0.436, R2 = 0.190, P = 0.011), and this relationship was stronger when the three samples from Crete were removed from the analysis (R XY  = 0.579, R2 = 0.335, P = 0.001). There was also a strong relationship when only samples from Peloponnese were analyzed (R XY  = 678, R2 = 0.460, P = 0.001). In Crete, however, isolation by distance was not detected (R XY  = 0.963, R2 = 0.927, P = 0.186). Isolation by distance in the field samples was much greater than that between greenhouses (R XY  = 0.791/-0.219, R2 = 0.625/0.048, P = 0.005/0.113; fields/greenhouses).

AMOVA indicated that genetic variation was not explained by host plant families, but there was a significant percentage of variation explained by habitat (fields vs. greenhouses), and among the samples within each of these groups (Table 4). Samples from fields (all in Peloponnese) displayed significantly greater genetic diversity and less population differentiation than those from greenhouses (allelic richness: 3.380 vs. 3.139, P = 0.006; Ho/H E : 0.497/0.536 vs. 0.404/0.454, P < 0.001; F ST : 0.012 vs. 0.078, P = 0.003; fields vs. greenhouses, respectively). Comparison of groups of samples differing in host plant family (Table 1) did not indicate differences in F statistics, heterozygosity or allelic richness (all P > 0.05) (data not shown).

Analysis with GESTE revealed that a model including habitat (field vs. greenhouse) had the highest posterior probability (Pr = 0.407) of explaining the genetic structure. The population structure in Greece was also partially explained by longitude (Pr = 0.260) followed by a model which combined the habitat and longitude factors (Pr = 0.174).

The Bayesian analysis of population structure in Greece indicated the presence of three major genetic clusters (Figure 3A). Samples from fields of West Peloponnese (WP 1–7) formed cluster 1 (K 1 ). Cluster 2 (K 2 ) consisted of samples from Crete and from two mainland locations (NP and MA 1), and the sample from West Greece (WG) had a mixed assignment (WG: K 1  = 0.475 and K 2  = 0.397). Samples AT and MA 2 formed cluster 3 (K 3 ). Mixed assignments were indicated for MA 2 (K 3  = 0.599 and K 2  = 0.279) and EP, which grouped either with the other samples from the Peloponnese (K 1 ), or with cluster 3 (K 3 ) (EP: K 3  = 0.515 and K 1  = 0.461).

Figure 3
figure 3

Genetic clusters of T. vaporariorum sampled in Greece assigned by Bayesian analysis implemented in STRUCTURE [[52]] . In each graph individuals are represented by vertical bars broken into K differently colored genetic clusters, with length proportional to assignment in each cluster. Probability of assignment of each individual is indicated from 0 to 1. Individuals are grouped by sample with the sample code listed at the bottom of the bar graph. A: Analysis of all 16 samples, 474 individuals, K = 3. B: Analysis of a subset of 8 samples collected from fields, 239 individuals, K = 3. C: Analysis of a subset of 8 samples collected from greenhouses, 235 individuals, K = 3.

Since habitat (field vs. greenhouse) was a major component of the genetic structure indicated in our initial analysis and our analysis with GESTE, we analyzed subsets of samples from greenhouses and fields separately in STRUCTURE to investigate possible sub-structuring (Figure 3B and C). Analysis of the subset sampled from fields (K 1 in the initial analysis) revealed the presence of three sub-clusters (K = 3). Majority of samples (WP 3-WP 6) were not well-resolved and assigned to both K 1.1 and K 1.2 , however samples WP 1 and WP 7 had higher assignment probability to K 1.1 (0.80 for WP1, and 0.825 for WP7). Sub-cluster 2 (K 1.2 ) was best represented by WP 2 (with assignment probability of 0.90). Sub-cluster 3 (K 1.3 ) consisted of the sample from East Peloponnese (EP: K 1.3  = 0.754). The subset of samples from greenhouses (K 2 and K 3 in the initial analysis) was characterized best by three well-resolved genetic sub-clusters (K = 3) (Figure 3C), with 71-98% assignment rate. Sub-cluster 1 (K 2.1 ) consisted of two Cretan samples (CR 1 and CR 2) and samples from Peloponnese and the mainland (NP, MA 1 and, in part, MA 2). CR 3, however, formed sub-cluster 2 (K 2.2 ) with WG and MA 2. MA 2 was only partially resolved and grouped both with K 2.1 (0.499) and K 2.2 (0.450). Thus, populations were differentiated even within a relatively short distance on the Island of Crete. AT formed a sub-cluster 3 (K 3.1 ).


Agricultural ecosystems can serve as temporal oases and increase distribution of many beneficial and pest species [5],[58]. The invasive pest greenhouse whitefly (Trialeurodes vaporariorum) has a cosmopolitan distribution [25],[26]. In this study, we examined population genetic structure among samples from its Northern (Finland) and Southern (Greece) European distribution range. There was no evidence of gene flow between the different climate zones. Dispersal of T. vaporariorum appears to be limited, since we found significant spatial population genetic structure among samples in both countries. Samples from Finland were less diverse and showed greater genetic differentiation than samples from Greece, which could be explained by differences in agroecosystems found in the different climate zones. In Greece, habitat (field vs. greenhouse) explained population genetic structure, but in Finland, genetic structure was dictated by host plant species. Related whiteflies in the B. tabaci complex also show population genetic structure [59]. However, population differentiation of B. tabaci (Med) in Tunisia is not related to host plant species nor is it related to type of agroecosystem in Greece [35],[60].

The role of agroecosystems

We hypothesized that greenhouses contribute to population genetic structure in T. vaporariorum by limiting dispersal and gene flow among populations, and that populations in the North are more genetically differentiated than populations in the South. In Finland, our samples of T. vaporariorum were limited to greenhouses (since in this country this species is not able to persist year-round in agricultural fields), whereas in Greece we sampled from both greenhouses and fields. The extent of the area inhabited by homogeneous populations was approximately 1–10 km straight line distance in Finland and approximately 100 km in Greece, indicated in our analyses as samples with non-significant differences in pairwise F ST (Finland: PJ 1–5, NR 1; Greece: WP 1–7), and shared cluster assignment in Bayesian analysis (Finland: PJ 1–5; Greece: WP 1–7). Whiteflies sampled at these spatial scales also showed significant isolation by distance. Although whiteflies are poor fliers, natural dispersal of up to 20 km has been observed previously for the related B. tabaci species [41]. Low pairwise F ST values and significant isolation by distance even between remote populations in Greece suggest that the presence of agricultural fields and wild host plants year-round [61] enables greater connectivity among T. vaporariorum populations. Nevertheless, greenhouse agroecosystems increase genetic structure in both climate zones. Similarly, populations of whiteflies, B. tabaci (Med), and moths, Trichoplusia ni, inhabiting greenhouses in the United States and Canada, respectively, were characterized by higher genetic structure than those inhabiting fields [62],[63].

Although frequent high mortality events occur in agroecosystems in both climate zones (via chemical and biological control methods), the high abundance of less enclosed/hermetic agroecosystems (fields, plastic tunnels and greenhouses with constantly open vents) and suitable climate in the South reduces effectiveness of pest management. On the other hand, isolation of populations in more enclosed agroecosystems, such as those in the North, can create an opportunity for the development of insecticide resistance through natural selection [39]. In open environments resistance genes might be less easily fixed, but they could spread over longer distances, as indicated by the higher connectivity of field populations in our study. Populations of insect pests inhabiting greenhouses are often characterized by higher insecticide resistance than those in the fields [58].

Agroecosystems, particularly greenhouses, could affect genetic structure of pests not only by placing limitations on pest dispersal, but also by limiting their population size. Crop management practices and frequent chemical insecticide exposure can cause population bottlenecks, leading to reduction in within-population genetic variation (increases in homozygosity), as well as increase in between-population genetic differentiation [18]. Higher genetic diversity (allelic richness and heterozygosities) of Greek T. vaporariorum populations indicated a larger gene pool and overall population size and possibly a lower frequency of bottlenecks in the South than in the North. This poses a threat to effective pest management in the Mediterranean region. However, we observed higher than expected homozygosity in samples from both countries (reflected by deviations from HWE; Table 3). Although deviations in HWE could reflect technical problems in genotyping, they could also result from population bottlenecks and inbreeding, and homozygote excess in insect populations inhabiting agroecosystems is not uncommon [59],[64]. Nevertheless, samples deviating from HWE did not occur more frequently in the North than in the South, and in Greece, samples deviating from HWE were collected from fields as well as from greenhouses.

Despite the frequent stochastic events in agroecosystems that can reduce genetic diversity, our results indicated that T. vaporariorum is able to persist over years in the same greenhouses in Finland. Samples collected from the same greenhouses in 2010 and 2011 were not genetically differentiated from each other in all but one case. Prevalence of the pest year-round might eventually allow further spread into natural ecosystems. Pest persistence in agroecosystems can create propagule pressure to natural habitats and favor utilization of wild host plants that surround greenhouses and fields (Ovcarenko et al. unpubl.). Adaptation to crop species could also lead to the development of a preference for particular wild host plants with similar chemistry, as was observed for Tetranychus urticae[15],[65].

The role of host plants

Occupation of agroecosystems and differences in their individual management in the two climate zones could potentially allow T. vaporariorum to specialize and adapt to particular host plant species or their cultivars. Although T. vaporariorum is a polyphagous insect, it is able to develop preference not only for certain plant species [66], but also for particular varieties or cultivars [67]. In Finland, major genetic clusters were characterized by the two common host plants (Figure 2). In Greece, however, host plant taxonomic family did not explain population structure. Absence of host associations in Greece may reflect the frequent alterations of host crops that are dictated by the market. On the other hand, the limited number of samples restricted to particular host plant species may have prevented detecting a possible association.

It is known that T. vaporariorum prefer cucumber hosts to tomato, and have higher fecundity and shorter development time on cucumber as well [28]. Thus, larger and more diverse populations would be expected on cucumber. On the contrary, our study indicated that Finnish samples from tomato were characterized by larger population size, with higher heterozygosity and allelic richness and lower values of F statistics than samples from cucumber. This might also be due to differences in individual management of agroecosystems. Cucumber crops are changed every three to four months, whereas tomato crops are maintained in greenhouses for nine to ten months, leading to more frequent reductions in pest population size in greenhouses growing cucumber than in those growing tomato. The result might also be a sampling effect: there were more greenhouses growing only tomato than those growing only cucumber (nine and five, respectively) (Table 1). However, our sampling reflects the tendency in Finland for more common cultivation of tomato than cucumber [68].

Associations between particular host plant species and genetic structure of T. vaporariorum populations might alternatively reflect different introduction sources, as well as varying frequencies of repeated introductions [69]. When there is a single introduction source, e.g. infested plant material from a single supplier, even isolated populations experiencing no gene flow will not show genetic differentiation, since all populations will share the same alleles present in the source [70]. One sample in Finland (KR) that formed a distinct cluster in our Bayesian analysis was collected from a site that has been producing its own cucumber cultivar since 2005. Although such conditions could favor adaptation of T. vaporariorum to this cultivar, KR is also somewhat geographically isolated from other agroecosystems and is surrounded by forests, which potentially limits insect dispersal. Therefore, it is difficult to determine if its genetic differentiation is due to adaptation or isolation. Other samples from Finland collected from cucumber (except NL) clustered together (LH, NR 2 and TJ 1–2; Figure 2C). Such genetic structure could reflect either adaptation to the host, or a common vector or origin of T. vaporariorum. Although growing different crops, greenhouses PJ1-5, NR1-3, NL, HR, ML, LH and TJ2 had all obtained seedlings from the same producer, and some similarities between these samples could be seen in our analysis with STRUCTURE (Figure 2A). However, origin of seedlings did not explain local F ST in our analysis with GESTE, whereas the model with host plant species explained the data best. In Greece, populations from greenhouses in Crete (CR 1 and CR 2) clustered with distant locations in the mainland (NP and MA 1; Figure 3A, C), suggesting that human-mediated transfer of whiteflies has occurred between these distant locations. Human-mediated transfer of B. tabaci between northern and southern regions of Greece has been noted previously [32]. Unfortunately, no information on the origin of seedlings in Greece was available for our analysis.

Limitations of the study

Our study compares pest populations living in different climate zones and subjected to different kinds of agroecosystems, so it may not be surprising that we have detected differences in population structure in the North and South. Moreover, our results might be affected by differences in the timing of sampling of T. vaporariorum in Finland and Greece, which matched periods of high insect abundance. In Finland maximum abundance peaks are in spring, whereas in Greece whiteflies are abundant year-round. Sampling T. vaporariorum in Finland in spring, when dispersal is low or non-existent might have facilitated detecting significant population structure. However, persistence of the populations over two years in the same greenhouses indicates that genetic structure is present despite the sampling period. Furthermore, no differences in population structure of B. tabaci in Greece were detected between late and early sampling periods from both field and greenhouse populations [35].


Greenhouse agroecosystems contribute to population genetic structure in T. vaporariorum by limiting gene flow among populations. Populations in Finland sampled from greenhouses are less diverse and more genetically differentiated than populations in Greece, collected from both greenhouses and fields. Within Greece, pest populations inhabiting greenhouse agroecosystems were more genetically differentiated than those inhabiting fields, and habitat (field vs. greenhouse) together with longitude explained population genetic structure. In contrast, host plant species (tomato vs. cucumber) explained population genetic structure in Finland. The differing influence of type of agroecosystem and potential host plant adaptation on population genetic structure of the pest in different climate zones highlights challenges for the management of a cosmopolitan invasive pest species.


  1. Hsieh C-H, Wang C-H, Ko C-C: Evidence from molecular markers and population genetic analyses suggests recent invasions of the Western North Pacific region by biotypes B and Q of Bemisia tabaci (Gennadius). Environ Entomol. 2007, 36: 952-961. 10.1603/0046-225X(2007)36[952:EFMMAP]2.0.CO;2.

    Article  PubMed  CAS  Google Scholar 

  2. Margaritopoulos JT, Kasprowicz L, Malloch GL, Fenton B: Tracking the global dispersal of a cosmopolitan insect pest, the peach potato aphid. BMC Ecol. 2009, 9: 13-10.1186/1472-6785-9-13.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bebber DP, Ramotowski MAT, Gurr SJ: Crop pests and pathogens move polewards in a warming world. Nat Clim Chang. 2013, 3: 985-988.

    Article  Google Scholar 

  4. Hickling R, Roy DB, Hill JK, Fox R, Thomas CD: The distributions of a wide range of taxonomic groups are expanding polewards. Glob Chang Biol. 2006, 12: 450-455. 10.1111/j.1365-2486.2006.01116.x.

    Article  Google Scholar 

  5. Hanafi A: Invasive pests and diseases: a challenge to IPM in greenhouse crops. Phytoparasitica. 2005, 33: 423-426. 10.1007/BF02981390.

    Article  Google Scholar 

  6. Van Driesche R, Enkegaard E: Invasive Species as Pests in Greenhouses: Forecasting, Preventing and Remediating Future Invasions. IOBC/WPRS Bull. 2002, 25 (1): 277-280.

    Google Scholar 

  7. Altieri MA: The ecological role of biodiversity in agroecosystems. Agric Ecosyst Environ. 1999, 74: 19-31. 10.1016/S0167-8809(99)00028-6.

    Article  Google Scholar 

  8. Grapputo A, Boman S, Lindström L, Lyytinen A, Mappes J: The voyage of an invasive species across continents: genetic diversity of North American and European Colorado potato beetle populations. Mol Ecol. 2005, 14: 4207-4219. 10.1111/j.1365-294X.2005.02740.x.

    Article  PubMed  CAS  Google Scholar 

  9. Piiroinen S, Lindström L, Lyytinen A, Mappes J, Chen YH, Izzo V, Grapputo A: Pre-invasion history and demography shape the genetic variation in the insecticide resistance-related acetylcholinesterase 2 gene in the invasive Colorado potato beetle. BMC Evol Biol. 2013, 13: 13-10.1186/1471-2148-13-13.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Richardson DM, Pyšek P, Rejmánek M, Barbour MG, Panetta FD, West CJ: Naturalization and invasion of alien plants: concepts and definitions. Divers Distrib. 2000, 6: 93-107. 10.1046/j.1472-4642.2000.00083.x.

    Article  Google Scholar 

  11. Loxdale H, Lushai G: Slaves of the environment: the movement of herbivorous insects in relation to their ecology and genotype. Philos Trans R Soc Lond B Biol Sci. 1999, 354 (January): 1479-1495. 10.1098/rstb.1999.0492.

    Article  PubMed Central  Google Scholar 

  12. Freeman J: Problems of stored products entomology in Britain arising out of the import of tropical products. Ann Appl Biol. 1976, 84: 120-124. 10.1111/j.1744-7348.1976.tb01740.x.

    Article  Google Scholar 

  13. Hill DS: Pests of Crops in Warmer Climates and Their Control. 2008, Springer, Dordrecht, The Netherlands

    Book  Google Scholar 

  14. Muñoz P, Antón A, Nuñez M, Paranjpe A, Ariño J, Castells X, Montero JI, Rieradevall J: Environmental impacts of greenhouse versus open-field tomato production in the Mediterranean region. Acta Hort. 2008, 801: 1591-1596.

    Article  Google Scholar 

  15. Perdikis D, Kapaxidi E, Papadoulis G: Biological control of insect and mite pests in greenhouse Solanaceous crops. Eur J Plant Sci Biotech. 2008, 2: 125-144.

    Google Scholar 

  16. Hoffmann AA, Willi Y: Detecting genetic responses to environmental change. Nat Rev Genet. 2008, 9: 421-432. 10.1038/nrg2339.

    Article  PubMed  CAS  Google Scholar 

  17. Tsagkarakou A, Navajas M: Genetic differentiation in Tetranychus urticae (Acari: Tetranychidae) from greenhouses in France. Exp Appl Acarol. 1999, 23: 365-378. 10.1023/A:1006293627880.

    Article  Google Scholar 

  18. Tsagkarakou A, Navajas M, Papaioannou-Souliotis P, Pasteur N: Gene flow among Tetranychus urticae (Acari: Tetranychidae) populations in Greece. Mol Ecol. 1998, 7: 71-79. 10.1046/j.1365-294x.1998.00305.x.

    Article  CAS  Google Scholar 

  19. Mound LA, Halsey SH: Whitefly of the world. a systematic catalogue of the Aleyrodidae (Homoptera) with host plant and natural enemy data. 1978, John Wiley and Sons, British Museum (Natural History)

    Google Scholar 

  20. Linnaniemi WM: Ästerochiton vaporiorum (Westw.) Suomessa. Medd Soc Fauna Flora Fenn. 1921, 7: 66-68.

    Google Scholar 

  21. Hulden L: The whiteflies and their parasites in Finland. Not Entomol. 1986, 66: 1-40.

    Google Scholar 

  22. Vappula NA: Pests of cultivated plants in Finland. Acta Entomol Fenn. 1965, 19: 163-169.

    Google Scholar 

  23. Michelakis SE: Problems in the application of biological control against Trialeurodes vaporariorum in unheated plastic glasshouses in Crete. Bull OEPP. 1986, 16: 423-427. 10.1111/j.1365-2338.1986.tb00299.x.

    Article  Google Scholar 

  24. Pelekassis C: A catalogue of the more important insects and other animals harmful to the agricultural crops of Greece during the last thirty-year period. Ann Inst Phytopathol Benaki. 1962, 5: 9-75.

    Google Scholar 

  25. Martin JH, Mifsud D, Rapisarda C: The whiteflies (Hemiptera: Aleyrodidae) of Europe and the Mediterranean Basin. Bull Entomol Res. 2000, 90: 407-448.

    PubMed  CAS  Google Scholar 

  26. Van Lenteren JC, Van Roermund H, Sütterlin S: Biological Control of Greenhouse Whitefly (Trialeurodes vaporariorum) with the Parasitoid Encarsia formosa: How Does It Work?. Biol Control. 1996, 6: 1-10. 10.1006/bcon.1996.0001.

    Article  Google Scholar 

  27. Inbar M, Gerling D: Plant-mediated interactions between whiteflies, herbivores, and natural enemies. Annu Rev Entomol. 2008, 53: 431-448. 10.1146/annurev.ento.53.032107.122456.

    Article  PubMed  CAS  Google Scholar 

  28. Van Lenteren JC, Noldus LPJJ: Whitefly Plant Relationships, Behavioral and Ecological Aspects. Whiteflies: Their Bionomics, Pest Status and Management. Edited by: Gerling D. 1990, Intercept Ltd, Andover, UK, 47-49.

    Google Scholar 

  29. Stenseth C: Cold-hardiness in eggs of Greenhouse Whitefly (Trialeurodes vaporiorum). IOBC-WPRS Bull. 1983, 6: 84-86.

    Google Scholar 

  30. Osborne L: Temperature-dependent development of greenhouse whitefly and its parasite, Encarsia formosa. Environ Entomol. 1982, 11: 483-485.

    Article  Google Scholar 

  31. Boykin LM, Bell CD, Evans G, Small I, De Barro PJ: Is agriculture driving the diversification of the Bemisia tabaci species complex (Hemiptera: Sternorrhyncha: Aleyrodidae)?: Dating, diversification and biogeographic evidence revealed. BMC Evol Biol. 2013, 13: 228-10.1186/1471-2148-13-228.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Gauthier N, Clouet C, Perrakis A, Kapantaidaki D, Peterschmitt M, Tsagkarakou A: Genetic structure of Bemisia tabaci Med populations from home range countries inferred by nuclear and cytoplasmic markers: impact on the distribution of the insecticide resistance genes. Pest Manag Sci. 2014, ᅟ: ᅟ-in press

    Google Scholar 

  33. Hsieh C-H, Chiang Y-H, Ko C-C: Population genetic structure of the newly invasive Q biotype of Bemisia tabaci in Taiwan. Entomol Exp Appl. 2011, 138: 263-271. 10.1111/j.1570-7458.2011.01097.x.

    Article  CAS  Google Scholar 

  34. Dalmon A, Halkett F, Granier M, Delatte H, Peterschmitt M: Genetic structure of the invasive pest Bemisia tabaci: evidence of limited but persistent genetic differentiation in glasshouse populations. Heredity. 2008, 100: 316-325. 10.1038/sj.hdy.6801080.

    Article  PubMed  CAS  Google Scholar 

  35. Tsagkarakou A, Mouton L, Kristoffersen JB, Dokianakis E, Grispou M, Bourtzis K: Population genetic structure and secondary endosymbionts of Q Bemisia tabaci (Hemiptera: Aleyrodidae) from Greece. Bull Entomol Res. 2012, 102: 353-365. 10.1017/S0007485311000757.

    Article  PubMed  CAS  Google Scholar 

  36. Prijović M, Skaljac M, Drobnjaković T, Zanić K, Perić P, Marčić D, Puizina J: Genetic variation of the Greenhouse Whitefly, Trialeurodes vaporariorum (Hemiptera: Aleyrodidae), among populations from Serbia and neighboring countries, as inferred from COI sequence variability. Bull Entomol Res. 2014, 1: 1-10.

    Google Scholar 

  37. Roopa HK, Kumar NKK, Asokan R, Rebijith KB, Mahmood R, Verghese A: Phylogenetic Analysis of Trialeurodes Spp. (Hemiptera: Aleyrodidae) from India Based on Differences in Mitochondrial and Nuclear DNA. Fla Entomol. 2012, 95: 1086-1094. 10.1653/024.095.0438.

    Article  Google Scholar 

  38. Shin D, Mo H, Lee S-E, Park J-J, Cho K: Elucidation of the genetic differences in Trialeurodes vaporariorum populations under vegetable greenhouse conditions by using the allozyme approach. Entomol Res. 2013, 43: 271-281. 10.1111/1748-5967.12041.

    Article  CAS  Google Scholar 

  39. Ovcarenko I, Lindström L, Saikkonen K, Vänninen I: Variation in mortality among populations is higher for pymetrozine, than for imidacloprid and spiromesifen in Trialeurodes vaporariorum in greenhouses in Finland. Pest Manag Sci. 2014, ᅟ: ᅟ-in press

    Google Scholar 

  40. Byrne DN, Bellows TS: Whitefly biology. Annu Rev Entomol. 1991, 36: 431-457. 10.1146/annurev.en.36.010191.002243.

    Article  Google Scholar 

  41. Berlinger MJ, Lehmann-Sigura N, Taylor RAJ: Survival of Bemisia tabaci adults under different climatic conditions. Entomol Exp Appl. 1996, 80: 511-519. 10.1111/j.1570-7458.1996.tb00967.x.

    Article  Google Scholar 

  42. Tsagkarakou A, Tsigenopoulos CS, Gorman K, Lagnel J, Bedford ID: Biotype status and genetic polymorphism of the whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) in Greece: mitochondrial DNA and microsatellites. Bull Entomol Res. 2007, 97: 29-40. 10.1017/S000748530700466X.

    Article  PubMed  CAS  Google Scholar 

  43. Aksoy S, Almeida-Val VMF, Azevedo VCR, Baucom R, Bazaga P, Beheregaray LB, Bennetzen JL, Brassaloti RA, Burgess TI, Caccone A, Chang S-M, Ciampi AY, Ciancaleoni S, Clímaco GT, Clouet C, Coimbra MRM, Coutinho LL, Dantas HL, De Vega C, Echodu R, Enyaru J, Figueira A, Filho MAG, Foltz B, Fressigné L, Gadomski M, Gauthier N, Herrera CM, Hyseni C, et al: Permanent genetic resources added to molecular ecology resources database 1 October 2012–30 November 2012. Mol Ecol Resour. 2013, 13: 341-343. 10.1111/1755-0998.12061.

    Article  PubMed  Google Scholar 

  44. Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.

    CAS  PubMed Central  Google Scholar 

  45. Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research--an update. Bioinformatics. 2012, 28: 2537-2539. 10.1093/bioinformatics/bts460.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  46. Rousset F: Genepop’007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Resour. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.

    Article  PubMed  Google Scholar 

  47. Raymond M, Rousset F: An exact test for population differentiation. Evolution. 1995, 49: 1280-1283. 10.2307/2410454.

    Article  Google Scholar 

  48. Rice WR: Analyzing tables of statistical tests. Evolution. 1989, 43: 223-225. 10.2307/2409177.

    Article  Google Scholar 

  49. Oosterhout C, Van Hutchinson WF, Wills DPM, Shipley P: Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.

    Article  Google Scholar 

  50. Goudet J: Fstat 2.9.3: A Program to Estimate and Test Gene Diversities and Fixation Indices (Updated from Goudet 1995). 2002, Switzerland, Lausanne

    Google Scholar 

  51. Foll M, Gaggiotti O: Identifying the environmental factors that determine the genetic structure of populations. Genetics. 2006, 174: 875-891. 10.1534/genetics.106.059451.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.

    PubMed  CAS  PubMed Central  Google Scholar 

  53. Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9: 1322-1332. 10.1111/j.1755-0998.2009.02591.x.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.

    Article  PubMed  CAS  Google Scholar 

  55. Earl DA, VonHoldt BM: STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2011, 4: 359-361. 10.1007/s12686-011-9548-7.

    Article  Google Scholar 

  56. Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806. 10.1093/bioinformatics/btm233.

    Article  PubMed  CAS  Google Scholar 

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

    Article  Google Scholar 

  58. Parrella MP: Arthropod Fauna. Ecosystems of the World. Greenhouse Ecosystems. Edited by: Stanhill G, Enoch HZ. 1999, Elsevier Science Publishers, Amsterdam; New York, 213-250.

    Google Scholar 

  59. Delatte H, David P, Granier M, Lett JM, Goldbach R, Peterschmitt M, Reynaud B: Microsatellites reveal extensive geographical, ecological and genetic contacts between invasive and indigenous whitefly biotypes in an insular environment. Genet Res. 2006, 87: 109-124. 10.1017/S0016672306008135.

    Article  PubMed  CAS  Google Scholar 

  60. Saleh D, Laarif A, Clouet C, Gauthier N: Spatial and host-plant partitioning between coexisting Bemisia tabaci cryptic species in Tunisia. Popul Ecol. 2012, 54: 261-274. 10.1007/s10144-012-0303-z.

    Article  Google Scholar 

  61. Roditakis NE: Host plants of greenhouse whitefly Trialeurodes vaporariorum westwood (Homoptera: Aleyrodidae) in Crete. Attractiveness and impact on whitefly life stages. Agric Ecosyst Environ. 1990, 31: 217-224. 10.1016/0167-8809(90)90221-X.

    Article  Google Scholar 

  62. Franklin MT, Ritland CE, Myers JH: Spatial and temporal changes in genetic structure of greenhouse and field populations of cabbage looper, Trichoplusia ni. Mol Ecol. 2010, 19: 1122-1133. 10.1111/j.1365-294X.2010.04548.x.

    Article  PubMed  Google Scholar 

  63. Dickey AM, Osborne LS, Shatters RG, Hall PM, Mckenzie CL: Population genetics of invasive Bemisia tabaci (Hemiptera: Aleyrodidae) cryptic species in the United States based on microsatellite markers. J Econ Entomol. 2013, 106: 1355-1364. 10.1603/EC12512.

    Article  PubMed  Google Scholar 

  64. Simón B, Cenis JL, De La Rúa P: Distribution patterns of the Q and B biotypes of Bemisia tabaci in the Mediterranean Basin based on microsatellite variation. Entomol Exp Appl. 2007, 124: 327-336. 10.1111/j.1570-7458.2007.00586.x.

    Article  Google Scholar 

  65. Navajas M, Tsagkarakou A, Lagnel J, Perrot-Minnot MJ: Genetic differentiation in Tetranychus urticae (Acari: Tetranychidae): polymorphism, host races or sibling species?. Exp Appl Acarol. 2000, 24: 365-376. 10.1023/A:1006432604611.

    Article  PubMed  CAS  Google Scholar 

  66. Ma R-Y, Kong W-N, Hao L-J: Host preference of Greenhouse Whitefly (Trialeurodes vaporariorum) for several horticultural plants in greenhouse. Entomol Knowl. 2005, 3: 016-

    Google Scholar 

  67. Thomas DC: Host plant adaptation in the Glasshouse Whitefly. J Appl Entomol. 1993, 115: 405-415. 10.1111/j.1439-0418.1993.tb00408.x.

    Article  Google Scholar 

  68. Official Statistics of Finland. ᅟ. ᅟ, ᅟ: ᅟ-[], []

  69. Whitlock M: Selection and Drift in Metapopulations. Ecology, Genetics and Evolution of Metapopulations. Edited by: Gaggiotti O, Hanski I. 2004, Elsevier Academic Press, Amsterdam, 153-173.

    Chapter  Google Scholar 

  70. Sakai A, Allendorf F, Holt J: The population biology of invasive species. Annu Rev Ecol Syst. 2001, 32: 305-332. 10.1146/annurev.ecolsys.32.081501.114037.

    Article  Google Scholar 

Download references


The Finnish Ministry of Education and Culture, Finnish Ministry of Agriculture and Forestry, BIOINT Doctoral Programme at the University of Jyväskylä and the Oiva Kuusisto Foundation provided financial support to I. Ovcarenko, and MTT and the Finnish Ministry of Agriculture and Forestry to I. Vänninen. The Academy of Finland provided financial support to L. Lindström (project 250248) and to K. E. Knott (project 258365). The Centre of Excellence in Biological Interactions Research provided laboratory materials for the research. The EU European Social Fund (ESF) and the Greek State, as well as the General Secretariat General for Research and Technology (SYNERGASIA 09SYN-22-617) provided financial support to D. E. Kapantaidaki and A. Tsagkarakou.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Irina Ovčarenko.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

IV, IO, AT were the authors of original idea. IO, LL and IV took part in sample collection in Finland, whereas DEK and AT collected samples from Greece. NG participated in microsatellite development. DEK extracted DNA of samples from Greece, whereas IO and KEK carried all further sample and data analysis steps of all samples. KEK, NG and AT advised on data analysis methods. IO drafted the manuscript. All authors participated in preparation of the manuscript, read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ovčarenko, I., Kapantaidaki, D.E., Lindström, L. et al. Agroecosystems shape population genetic structure of the greenhouse whitefly in Northern and Southern Europe. BMC Evol Biol 14, 165 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Trialeurodes vaporariorum
  • Pest management
  • Microsatellite markers
  • Climate zone
  • Host adaptation