Translator Disclaimer
1 June 2014 Analysis of Genetic Diversity and Structure of Two Clades of Aphelinus mali (Hymenoptera: Aphelinidae) in China
Author Affiliations +

Our prior research revealed that there are 2 mitochondrial clades of Aphelinus mali (Haldeman) (Hymenoptera: Aphelinidae) in China, which are known as SD clade and LN clade. To further reveal their genetic characteristics and to determine the degrees of hybridization and gene flow between the 2 clades of A. mali in China, we analyzed the genetic diversities and genetic structures of 16 populations from 6 provinces (Shandong, Liaoning, Hebei, Shanxi, Xinjiang, and Yunnan) using 8 micro satellite loci. Our results showed that among the pure populations in the SD and LN clades, the greatest genetic diversities were found in the Qingdao, Shandong (QD) population and in the Dalian, Liaoning (DL) population. QD was the first population of the SD clade to be established, and DL was the first population of the LN clade to be established. In addition, genetic diversity was not substantially lower - and in some cases it was greater - in mixed-clade populations than in QD and DL. Individuals within each mitochondrial clade could not be differentiated based on microsatellite loci. Our data confirmed that the QD and DL populations, which were the first to be established in China, have served as bridgeheads for the other SD and LN populations in China. The results demonstrated that the hybridization or gene flow has occurred between the 2 mitochondrial clades.

Aphelinus mali (Haldeman) (Hymenoptera: aphid, Eriosoma lanigerum (Hausmann) (HemipAphelinidae), a parasitoid of the woolly apple tera: Aphididae) originated from North America (Beers 2012; Lavandero and Tylianakis 2013). Aphelinus mali was introduced into 51 countries during the 20th century to control E. lanigerum, and this parasitoid successfully established populations in 42 of these countries (Zhou et al. 2010).

Aphelinus mali was introduced as an endoparasitoid of E. lanigerum in China during the 1940s-1960s. The first introduction involved the importation of A. mali from Japan to Dalian, Liaoning in 1942 (Long et al. 1960). It was reintroduced from the former Soviet Union into Qingdao, Shandong during 1953–1955 (Long et al. 1960). Since then, A. mali has become established in most regions of China including Shandong, Liaoning, Shanxi, Yunnan, Hebei, Henan, and Xinjiang (Zhou et al. 2010; Zhang et al. 2014). Our field survey showed that all of the mitochondrial cytochrome c oxidase subunit I (mt COI) haplotypes grouped into 2 main clades, i.e., SD and LN. The SD clade was mainly distributed in Shandong and neighboring regions while the LN clade was mainly distributed in Liaoning and neighboring regions (Zhang et al. 2014). However, the genetic diversity and genetic structure of the A. mali in China based on nuclear markers have not been analyzed. Such investigations may provide knowledge of the hybridization and gene flow between the 2 clades of A. mali, which may be helpful in the control of the E. lanigerum in China. Currently, microsatellite markers have been extensively used in the analysis of genetic diversity and genetic structure (e.g., Ascunce et al. 2011; Chu et al. 2013; Dickey et al. 2013; Mwacharo et al. 2013).

Table 1.

Summary population genetics statistics calculated based on 8 microsatellites of 16 Aphelinus MALI POPULATIONS IN CHINA.


In this study, we analyzed the genetic diversity and genetic structure of 16 populations from 6 provinces (Shandong, Liaoning, Hebei, Shanxi, Xinjiang, and Yunnan) using 8 nuclear microsatellite loci. Based on these data, we wanted to reveal the genetic characteristics of the introduced populations of A. mali in China, and to determine the hybridization and gene flow within the 2 clades of A. mali in regions of China where they have been recently introduced.

Supporting material can be found online in Florida Entomologist 97(2) (June, 2014) at

Materials and Methods

Sample Collection and Species Identification

A total of 16 E. lanigerum populations were sampled from apple trees in 6 provinces in China during 2007, 2008, and 2012. They were placed in Petri dishes at room temperature until A. mali eclosed. Adult A. mali that emerged were examined and identified unambiguously with the aid of a microscope and were then stored in 95% ethanol at -20 °C until DNA extraction. The sampling location of each population is listed in Table 1.

DNA Extraction and Microsatellite Genotyping

Genomic DNA was extracted from individual female adults as described in Frohlich et al. (1999). The DNA of A. mali individuals and the PCR primers were used to amplify the microsatellite loci Aml3, Aml4, Aml9, Am27, Am34, Am35, Am36, and Am38 (Lavandero & Dominguez 2010). PCR reactions were performed as described by Lavandero & Dominguez (2010). Products were run on an ABI 3730xl DNA analyzer. Allele size was determined by comparing the mobility of the PCR products to that of the GeneScanTM 400HD size standard (Applied Biosystems).

Genetic Diversity Analysis

For the microsatellite data of each population, POPGENE vl.31 was used to calculate the following indices: the observed number of alleles (Na), the effective number of alleles (Ne), the expected heterozygosity (He), the observed heterozygosity (Ho), and Nei's expected heterozygosity (Nei) (Yeh et al. 1997). We determined these indices for both pure and mixed populations; pure populations were those that contained only one haplotype, and mixed populations were those that contained more than one haplotype (Table S1). For mixed populations, we also determined these indices separately for each clade associated with each haplotype. The program FSTAT v2.9.3.2 was used to calculate gene diversity (Hs) and allelic richness (Ar) (Franks et al. 2010).

GENEPOP v3.4 (Weir & Cockerham 1984) was used to test for linkage disequilibrium between all pairs of loci and to test for deviation from Hardy-Weinberg equilibrium (HWE) at each locus for each population; this was done with Markov chain parameters of 10,000 dememorization steps, followed by 1,000 batches of 10,000 iterations per batch (Raymond & Rousset 1995) in combination with Fisher's exact tests. The BOTTLENECK software was used to test for deviation from mutation-drift equilibrium (Cornuet & Luikart 1996). To identify the bottleneck event for each population, we investigate the occurrence of significant heterozygosity excess (which is a signature of a bottleneck) within all of the populations within clades under 3 mutational models (TPM, two-phase model of mutation; IAM, infinite alleles model; and SMM, stepwise mutation model). The TPM model was performed with the default settings of 30% variation from the IAM model and 70% from the SMM model.

A regular linear regression of the relationship between expected heterozygosity (He) and geographical distance from the nearest likely source population (the region of introduction mentioned in historical reports) was tested using Sigmaplot 12.0. The Google Earth ruler tool (Anderson & Congdon 2013) was used to calculate the Euclidean distance.

Analysis of Genetic Structure and Gene Flow

Differentiation approach based on FST analysis was used. GENEPOP v3.4 (Raymond & Rousset 1995) was used to calculate Weir and Cockerham's estimator of the fixation index FST (Weir & Cockerham 1984). Isolation by distance (IBD) based on genetic differentiation and the logarithm of geographic distance was examined by using 10,000 permutations of the Mantel test implemented in IBDWS v3.15 ( (Jensen et al. 2005; Hasselman et al. 2013). Distances between sampling locations were calculated using the Google Earth ruler tool (Anderson & Congdon 2013). Gene flow among regions was approximated as Nm (analogous to M = (1/FST -1)/2) (Slatkin 1993).

To evaluate the assignment of A. mali individuals to clusters, we also analyzed all of the populations and the 6 mixed populations using STRUCTURE software (Pritchard et al. 2000). For all populations analysis, we set K (number of clusters) from 1 to 16 with 6 iterations of this parameter set, using a burn-in of 20 000 followed by 50 000 iterations of Markov Chain Monte Carlo (MCMC). For 6 mixed populations analysis, we set K (number of clusters) from 1 to 6. Other parameters were similar with aforementioned. The most optimal K was estimated by examining the standardized second-order change of ◃K (Evanno et al. 2005). Sigmaplot 12.0 was used to graphically display the results.


Microsatellite-Based Genetic Diversity

Among the 6 pure populations in the SD clade (i.e., populations that were 100% Hapl) (Zhang et al. 2014), the genetic diversity was highest in the QD population, which was the first population in that clade that was introduced, and was lowest in the ZT population (Table 1). Genetic diversity was not substantially lower and in some cases was higher in the CZ (Changzhi, Shanxi), YC (Yuncheng, Shanxi), WF (Weifang, Shandong), and YT (Yantai, Shandong) mixed populations of the SD clade than in the QD population (Table 1; Table S1).

Among the 4 pure populations in the LN clade (i.e., populations that were 100% Hap2) (Zhang et al. 2014), the genetic diversity was highest in the DL population, which was the first population in that clade that was introduced (Table 1). As was the case with mixed populations in the SD clade, genetic diversity was not substantially lower and in some cases was higher in the CZ (Changzhi, Shanxi), BD (Baoding, Hebei), WF(Weifang, Shandong), and YT (Yantai, Shandong) mixed populations of the LN clade (Table 1; Table S1). Geographic distance was not correlated with expected heterozygosity (He) for the SD clade (R2 = 0.3825, F = 4.9548, df = 1, P = 0.0567) or the LN clade (R2 = 0.0431, F = 0.3152, df = l,P = 0.5920) (Fig. 1).

Except for the YL population, all populations with small sample sizes significantly deviated from HWE. A strong inbreeding was evident (0.424 < Fis < 1.000), and the deviations were associated with a significant positive Fis value. Null alleles were not a major factor contributing to deviations from HWE because PCR amplifications were successful for 98.1% of each locus across all populations (Table S2). We could not use software to test for null alleles because the software assumes random mating. When tested for deviation from mutation-drift equilibrium in BOTTLENECK under the IAM, 9 of 20 populations had significant heterozygosity excess (Wilcoxon test P < 0.05) (Table S3). However, a significant heterozygosity excess was detected in only 4 populations under the TPM model and in only one population under the SMM model (Table S3).

Population Genetic Structure and Gene Flow When paired populations were considered, all of the 120 F ST values were significant (Table 2). Isolation by distance analysis revealed a significant positive correlation between genetic distance and the logarithm of geographical distance (Fig. 2).

Estimates of gene flow between paired populations (Table 2) revealed that QHD and BD had the highest gene flow (Nm = 39.183) and that gene flow was also relatively high (Nm > 10) for YC and CZ, BD and DL, QHD and DL, TA and HZ, QD and WF, and QD and YT (see population codes in Table 1). While ZT and YL had the lowest gene flow (Nm = 0.141), and both of these populations also had low gene flow with other populations.

Fig. 1.

Linear regression analysis of the relationship between expected heterozygosity (He) and log-transformed Euclidean geographical distance from the nearest source population of Aphelinus mali in China. Data from mixed populations were separated with respect to the nearest likely source population; for example, WF data for Hapl (indicating the SD clade and the QD source population) and Hap 2 (indicating the LN clade and the DL source population) are shown in panel A and B, respectively. (A) The regression with QD as the source population of the SD clade (Y = -9.6526 X + 0.4943, R2 = 0.3825, F = 4.9548, df = 1, P = 0.0567). (B) The regression with DL as the source population of LN clade (Y = -8.2850 X + 0.4421, R2 = 0.0431, F = 0.3152, df = 1, P = 0.5920). Data for He in mixed populations are from Table S1, and data for pure populations are from Table 2. The Google Earth ruler tool was used to calculate the Euclidean geographical distances. Reduced Major Axis regression lines are shown.


According to STRUCTURE software, the genetic structure of all of the 16 populations in China was best described as consisting of 2 clusters (K = 2) (Fig. 3). The result demonstrated that the individuals within mitochondrial clades cannot be differentiated based on nuclear microsatellite loci.


Our study indicates that the nuclear genetic diversity among pure populations was highest in the initially introduced population of each clade, which was the QD population for the SD clade and the DL population for the LN clade; genetic diversity was much lower in the other pure populations (Table 1). The result suggests that the current distribution of A. mali in the introduced range is closely associated with the expansion from the 2 initially established populations, which acted as bridgehead populations, as has been shown for several invasive insects (Miller et al. 2005; Ciosi et al. 2008; Ascunce et al. 2011; Lombaert et al. 2010; Lombaert et al. 2011; Kajita et al. 2012; Yang et al. 2012). The result also confirms that the first two introduced populations, i.e., those in Qingdao, Shandong and Dalian, Liaoning, have served as the bridgeheads for the other populations in China as revealed by mitochondrial DNA markers (Zhang et al. 2014).

Table 2.



Fig. 2.

Relationship between genetic distance and logarithm of geographical distance based on all possible pairs of 16 Chinese populations of Aphelinus mali. Google Earth ruler tool was used to calculate the Euclidean geographical distances. Reduced Major Axis regression line is shown (Y = 0.3568 X - 0.6997, R = 0.7479, F = 149.7562, df = 1, P < 0.0001).


The initial founder effects may result in high genetic diversity near introduction points, but serial founder effects of other introduced populations may cause lower genetic diversity (Kajita et al. 2012). If the dispersal of A. mali is limited by distance, genetic diversity and geographical distances should be negatively correlated. However, the genetic diversity of A. mali was not correlated with geographical distance for either clade in China. That might be explained by anthropogenic influence on the spread of A. mali; the frequent transport of fruit seedlings may have facilitated the spread of the natural enemy of the woolly apple aphid. Another possibility is that the bottleneck effect during the subsequent introductions of A. mali was small or was mitigated by gene flow. Alien species normally experience founder or bottleneck effects, resulting in a loss of genetic diversity (Nei et al. 1975; Templeton 1980; Barton & Charlesworth 1984).

Fig. 3.

Bayesian clustering analysis of 16 Aphelinus mali populations in China using STRUCTURE (here K = 2). Each individual is represented by a vertical bar displaying membership coefficients to each genetic cluster. Population codes (City, Province): CZ, Changzhi, Shanxi; JZ, Jinzhong, Shanxi; YC, Yuncheng, Shanxi; YL, Yili, Xinjiang; DL, Dalian, Liaoning; HLD, Huludao, Liaoning; SJZ, Shijiazhuang, Hebei; BD, Baoding, Hebei; QHD, Qinhuangdao, Hebei; ZT, Zhaotong, Yunnan; HZ, Heze, Shandong; LC, Liaocheng, Shandong; TA, Taian, Shandong; WF, Weifang, Shandong; YT, Yantai, Shandong; QD; and Qingdao, Shandong


In our study, the reduction of genetic diversity within some populations might have been caused by bottlenecks (Table S3). However, most of the populations did not exhibit a heterozygosity excess (especially according to TPM or SMM models), which suggests that the bottleneck effect is a transient feature. High gene flow between populations may also obscure the genetic effects of a bottleneck (Table 2). For example, in the 4 mixed populations of the SD clade, the nuclear genetic diversity had not substantially decreased in YT and YC relative to that in QD and had even increased in WF and CZ (Fig. 1A). A similar phenomenon occurred in LN clade (Fig. 1B). Compared to the nuclear genetic diversity of the bridgehead population (DL) of the LN clade, the nuclear genetic diversity had not substantially decreased in YT and BD and had even increased in WF and CZ. In conclusion, bottleneck effects or gene flow may affect the correlation between genetic diversity and distance from the first established populations.

Hybridization or gene flow occurred between the 2 mitochondrial clades based on nuclear microsatellite loci (Fig. 3). The individuals within mitochondrial clades could not be differentiated based on nuclear microsatellite loci. A positive Fis index significantly different from zero (Table 1) indicates heterozygote deficiency and suggests inbreeding. This finding is consistent with the possibility of high gene flow between populations, which may obscure the genetic effects of a bottleneck. High gene flow between populations is not surprising, given the geographical proximity of the 2 release sites.

High gene flow can limit the genetic distance between nearby populations while low gene flow can increase the genetic distance between distant populations (Scott et al. 2005). Because of the high gene flow among plant species and locations, A. mali populations in Chile exhibit very low genetic differentiation (Lavandero et al. 2011). In the current study, however, FST analyses showed that A. mali populations in China have clearly undergone genetic differentiation in the expansion range. The positive correlation between genetic and geographic distance (Fig. 2) is associated with the low flight capacity of A. mali, which limits its natural dispersal ability. With low or moderate migration, genetic differentiation between regions and a significant isolationby-distance effect (Scott et al. 2005) occurred. Additionally, these results suggest a tendency for A. mali immigration from nearby rather than from distant populations in spite of accidental human transportation, a tendency that is consistent with the data concerning gene flow between populations (Table 2).

Invasive species have become serious problems in many regions because of their widespread and large effects on ecosystems, economies, and societies (Lockwood et al. 2007). Unraveling the dispersal pattern of invasive species is crucial for predicting and preventing additional invasions (Lockwood et al. 2007; Wilson et al. 2009). The study of the dispersal pattern of the invasive species is difficult (Chu et al. 2013), however, because the initial establishment of the invasive species can rarely be detected; the invasive species are typically not detected until a substantial time has passed after the initial introduction. Information about the further spread of “invasive” species can be gained, however, by the study of intentionally introduced biological control agents. Such agents provide an opportunity to test the effects of bridgehead populations and further spread because the timing, duration, and intensity of a demographic bottleneck are often well documented (Sakai et al. 2001; Hufbauer & Roderick 2005; Marsico et al. 2009; Franks et al. 2010).

To our knowledge, this is the first study to use the nuclear microsatellite loci to determine the hybridization and gene flow within 2 clades of A. mali in China. The genetic admixture of populations arising from different bridgeheads may help the biological control agent to adapt to the conditions in the introduced regions. The results obtained with an introduced biological control agent suggest that invasion by a species with limited dispersal ability may be facilitated by the establishment of multiple bridgehead populations.

Supplementary material for this article in Florida Entomologist 97(2) (June, 2014) is online at This material includes the following 3 tables of additional supporting information:


  • Table S2. Allele sizes at eight microsatellite loci SCREENED FOR APHELINUS MALI

  • Table S3. P values for within-population tests for HETEROZYGOSITY EXCESS IN 16 APHELINUS MALI POPULATIONS IN CHINA.


Hong-Xu Zhou and Rui-Ming Zhang contributed equally to this work. We are grateful to Dr. Gabor Lovei (Department of Agroecology, Aarhus University) and Blas Lavandero (Instituto de Biología Vegetal y Biotecnología, Laboratorio de Interacciones Insecto-Planta, Universidad de Talca, Chile) for his detailed comments on this manuscript. This work was supported by the National Natural Science Foundation (31371994), the National Key Basic Research Development Plan Project (2013CB127600), the Science and Technology Development Planning Program of Qingdao (13-1-3-108-nsh), the Special Fund for Agro-scientific Research in the Public Interest (201103026-5-2), and the Taishan Mountain Scholar Constructive Engineering Foundation of Shandong, China.

References Cited


K. L. Anderson , and B. C. Concdon 2013. Population genetics suggest that multiple invasion processes need to be addressed in the management plan of a plant disease vector. Evol. Appl. 4: 660–672. Google Scholar


M. S. Ascunce , C. C. Yang , J. Oakey , L. Calcaterra , and W. J. Wu et al. 2011. Global invasion history of the fire ant Solenopsis invicta. Science 331: 1066–1068. Google Scholar


N. H. Barton , and B. Charlesworth 1984. Genetic revolutions, founder effects and speciation. Annu. Rev. Ecol. Syst. 15: 133–164. Google Scholar


E. H. Beers 2012. Effect of trap color and orientation on the capture of Aphelinus mali (Hymenoptera: Aphelinidae), a parasitoid of woolly apple aphid (Hemiptera: Aphididae). J. Econ. Entomol. 105: 1342–1349. Google Scholar


D. Chu , H. P. Pan , X. C. Li , D. Guo , Y. L. Tao , B. M. Liu , and Y. J. Zhang 2013. Spatial genetic heterogeneity in populations of a newly invasive whitefly in China revealed by a nation-wide field survey. PLoS ONE 8: e26694. Google Scholar


M. Ciosi , N. J. Miller , K. S. Kim , R. Giordano , and A. Estoup et al. 2008. Invasion of Europe by the western corn rootworm, Diabrotica virgifera virgifera: multiple transatlantic introductions with various reductions of genetic diversity. Mol. Ecol. 17: 3614–3627. Google Scholar


J. M. Cornuet , and G. Luikart 1996. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 144: 2001–2014. Google Scholar


A. M. Dickey , L. S. Osborne , and R. G. Shatters , P. M. Hall , and C. L. Mckenzie 2013. Population genetics of invasive Bemisia tabaci (Hemiptera: Aleyrodidae) cryptic species in the United States based on microsatellite markers. J. Econ. Entomol. 106: 1355–1364. Google Scholar


G. Evanno , S. Regnaut , and J. Goudet 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14: 2611–2620. Google Scholar


S. J. Franks , P. D. Pratt , and N. D. Tsutsui 2010. The genetic consequences of a demographic bottleneck in an introduced biological control insect. Conse rv. Genet. 12: 201–211. Google Scholar


D. R. Frohlich , I. Torres-Jerez , and I. D. Bedford et al. 1999. A phylogeographical analysis of the Bemisia tabaci species complex based on mitochondrial DNA markers. Mol. Ecol. 8: 1683–1691. Google Scholar


D. J. Hasselman , D. Ricard , and P. Bentzen 2013. Genetic diversity and differentiation in a wide ranging anadromous fish, American shad (Alosa sapidissima), is correlated with latitude. Mol. Ecol. 22: 1558–1573. Google Scholar


R. A. Hufbauer , and G. K. Roderick 2005. Microevolution in biological control: mechanisms, patterns, and processes. Biol. Control 35: 227–239. Google Scholar


J. L. Jensen , A. J. Bohonak , and S. T. Kelley 2005. Isolation by distance web service. BMC Genetics 6: 13. Google Scholar


Y. Kajita , E. M. O'Neill , Y. B. Zheng , J. J. Obrycki , and D. W. Weisrock 2012. A population genetic signature of human releases in an invasive ladybeetle. Mol. Ecol. 21: 5473–5483. Google Scholar


B. Lavandero , M. And Dominguze 2010. Isolation and characterization of nine microsatellite loci from Aphelinus mali (Hymenoptera: Aphelinidae), a parasitoid of Eriosoma lanigerum (Hemiptera: Aphididae). Insect Sci. 17: 549–552. Google Scholar


B. Lavandero , C. C. Figueroa , P. Franck , and A. Mendze 2011. Estimating Gene Flow between Refuges and Crops: A Case Study of the Biological Control of Eriosoma lanigerum byAphelinus mali in Apple Orchards. PLoS ONE 6: e26694. Google Scholar


B. Lavandero , And J. M. Tylianakis 2013. Genotype matching in a parasitoid-host genotypic food web: an approach for measuring effects of environmental change. Mol. Ecol. 22: 229–238. Google Scholar


J. L. Lockwood , M. F. Hoopes , and M. P. Marchetti 2007. Invasion Ecology. Malden, MA: Blackwell Press. Google Scholar


E. Lombaert , T. Guillemaud , J. M. Cornuet , T. MaLausa , B. And Facon et al. 2010. Bridgehead effect in the worldwide invasion of the biocontrol harlequin ladybird. PLoS ONE 5: e9743. Google Scholar


E. Lombaert , T. Guillemaud , C. E. Thomas , L. J. Lawson Handley , and J. Li et al. 2011. Inferring the origin of populations introduced from a genetically structured native range by approximate Bayesian computation: case study of the invasive ladybird Harmonia axyridis. Mol. Ecol. 20: 4654– 4670. Google Scholar


C. D Long , Y. P. Wang , and P. Z. Tang 1960. Investigations on the biology and utilization of Aphelinus mali Hald., the specific parasite of the woolly apple aphis, Eriosoma Lanigerum Hausm. Acta Entomol. Sinica 10: 1–39. Google Scholar


T. D. Marsico , J. W. Burt , E. K. Espeland , G. W. Gilchrist , and M. A. Jamieson et al. 2009. Underutilized resources for studying the evolution of invasive species during their introduction, establishment, and lag phases. Evol. Appl. 3: 203–219. Google Scholar


J. M. Mwacharo , K. Nomura , and H. Kanada et al. 2013. Reconstructing the origin and dispersal patterns of village chickens across East Africa: insights from autosomal markers. Mol. Ecol. 22: 2683–2697. Google Scholar


N. Miller , A. Estoup , S. Toepfer , D. Bourguet , L. Lapchin et al. 2005. Multiple transatlantic introductions of the western corn rootworm. Science 310: 992. Google Scholar


M. Nei , T. Maruyama , and R. Chakraborty 1975. The bottleneck effect and genetic variability in populations. Evolution 29: 1–10. Google Scholar


J. K. Pritchard , M. Stephens , and P. Donnelly 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945–959. Google Scholar


M. Raymond , and F. Rousset 1995. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J. Hered. 86: 248–249. Google Scholar


A. K. Sakai , F. W. Allendorf , J. S. Holt , D. M. Lodge , and J. Molofsky et al. 2001. The population biology of invasive species. Annu. Rev. Ecol. Syst. 32: 305–332. Google Scholar


K. D. Scott , K. S. Wilkinson , N. Lawrence , C. L. Lange , AND L. J. Scott et al. 2005. Gene-flow between populations of cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) is highly variable between years. B. Entomol. Res. 95: 381–392. Google Scholar


M. Slatkin 1993. Isolation by distance in equilibrium and nonequilibrium populations. Evolution 47: 264– 279. Google Scholar


A. R. Templeton 1980. The theory of speciation via the founder principle. Genetics 94: 11–20. Google Scholar


B. S. Weir , and C. C. Cockerham 1984. Estimating F-statistics for the analysis of population structure. Evolution 38: 1358–1370. Google Scholar


J. R. Wilson , E. E. Dormontt , P. J. Prentis , A. J. Lowe , and D. M. Richardson 2009. Something in the way you move: dispersal pathways affect invasion success. Trends Ecol. Evol. 24: 136–144. Google Scholar


X. M. Yang , J. T. Sun , X. F. Xue , J. B. Li , and X. Y. Hong 2012. Invasion genetics of the western flower thrips in China: evidence for genetic bottleneck, hybridization and bridgehead effect. PLoS ONE 7: e34567. Google Scholar


F. C. Yeh , R. C. Yang , T. Boyle , Z. H. Ye , and J. X. Mao 1997. POPGENE, the user-friendly shareware for population genetic analysis. Molecular Biology and Biotechnology Centre, Univ. Alberta, Edmonton, Canada. Google Scholar


R. M. Zhang , H. X. Zhou , D. Guo , Y. L. Tao , F. H. Wan , Q. Wu , and D. Chu 2014. Two putative bridgehead populations of Aphelinus mali (Hymenoptera: Aphelinidae) in china as revealed by mitochondrial DNA marker. Florida Entomol. 97(2): In press. Google Scholar


H. X. Zhou , J. Y. Guo , F. H. Wan , and Y. Yu 2010. Natural control of Aphelinus mali on Eriosoma lanigerum and the parasitoid’s over wintering protection measures. Acta Phytophylacica Sinica 37: 153–158. Google Scholar
Hong-Xu Zhou, Rui-Ming Zhang, Dong Guo, Yun-Li Tao, Fang-Hao Wan, Qiang Wu, and Dong Chu "Analysis of Genetic Diversity and Structure of Two Clades of Aphelinus mali (Hymenoptera: Aphelinidae) in China," Florida Entomologist 97(2), 699-706, (1 June 2014).
Published: 1 June 2014

Back to Top