High genetic diversity in the offshore island populations of the tephritid fruit fly Bactrocera dorsalis
© The Author(s) 2016
Received: 20 July 2016
Accepted: 6 October 2016
Published: 13 October 2016
Geographic isolation is an important factor that limit species dispersal and thereby affects genetic diversity. Because islands are often small and surrounded by a natural water barrier to dispersal, they generally form discrete isolated habitats. Therefore, islands may play a key role in the distribution of the genetic diversity of insects, including flies.
To characterize the genetic structure of island populations of Bactrocera dorsalis, we analyzed a dataset containing both microsatellite and mtDNA loci of B. dorsalis samples collected from six offshore islands in Southern China. The microsatellite data revealed a high level of genetic diversity among these six island populations based on observed heterozygosity (Ho), expected heterozygosity (HE), Nei’s standard genetic distance (D), genetic identity (I) and the percentage of polymorphic loci (PIC). These island populations had low F ST values (F ST = 0.04161), and only 4.16 % of the total genetic variation in the species was found on these islands, as determined by an analysis of molecular variance. Based on the mtDNA COI data, high nucleotide diversity (0.9655) and haplotype diversity (0.00680) were observed in all six island populations. F-statistics showed that the six island populations exhibited low or medium levels of genetic differentiation among some island populations. To investigate the population differentiation between the sampled locations, a factorial correspondence analysis and both the unweighted pair-group method with arithmetic mean and Bayesian clustering methods were used to analyze the microsatellite data. The results showed that Hebao Island, Weizhou Island and Dong’ao Island were grouped together in one clade. Another clade consisted of Shangchuan Island and Naozhou Island, and a final, separate clade contained only the Wailingding Island population. Phylogenetic analysis of the mtDNA COI sequences revealed that the populations on each of these six islands were closely related to different populations on mainland China.
Our study suggests that these island populations have high genetic diversity, experience frequent gene flow and exhibit low or medium levels of genetic differentiation among some island populations. Therefore, the geographic isolation of the six islands does not appear to be a major dispersal barrier to B. dorsalis. Such knowledge is helpful for a better understanding of evolutionary processes of the species of island populations.
Genetic diversity is a critical component of biodiversity and affects the survival and evolution of species [1, 2]. Geographic isolation is an important factor that limits species dispersal and affects the genetic diversity of species [3, 4]. By definition, islands are smaller than continents, surrounded by water, and therefore form discrete habitats isolated from other terrestrial habitats [5, 6]. Generally, the water surrounding islands acts as a geographic barrier to dispersal that limits gene flow both between island populations and between island and mainland populations. Consequently, the genetic diversity of insular species tends to be more complex due to multiple factors that include the natural dispersal distance of flight-capable species and dispersal that is mediated by human activity. In such special cases, island populations have lower genetic diversity than do mainland populations due to founder effects/bottlenecks and continued isolation from the mainland [3, 7–10]. To maintain a population’s fitness and adapt to an island environment, a population may lose genetic diversity . Previous bottlenecks or continued isolation from the mainland can also decrease genetic diversity. Moreover, the loss of genetic diversity in island populations can often be caused by founder effects, breeding rates, and dispersal ability, among other factors. However, the island populations of some species, i.e., mammals and birds, have higher levels of genetic diversity than do mainland populations [8, 11, 12].
Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), the oriental fruit fly, is a quarantined pest that is found worldwide, is highly fecund and highly adaptable. This fly species causes serious economic losses to fruit production in tropical and subtropical areas [13–16]. The original source of B. dorsalis populations has been reported to be from tropical or subtropical Asia . To date, B. dorsalis has spread over Asia and to many regions around the Pacific, including Hawaii . Bactrocera dorsalis was first recorded in Taiwan in 1911 but has since spread to many other provinces in China [19, 20]. Wan et al.  examined the genetic diversity and genetic differentiation of B. dorsalis in the Chongqing region in China. Their results indicated that the height of these mountains in this region was insufficient to prevent the long-distance dispersal of B. dorsalis and suggested that there would be a high frequency of gene flow among fly populations. However, Li et al.  demonstrated that geographic isolation from mountains and canyons distributed across China, Vietnam and Thailand slowed the dispersal of B. dorsalis and has resulted in genetic differentiation between these regions. To date, no studies have explored the genetic diversity of B. dorsalis on islands. In this study, we hypothesized that six island populations of B. dorsalis in South China would have lower levels of genetic diversity due to their isolation from the mainland. Furthermore, the isolation of islands may also result in high levels of genetic differentiation between different island populations. In this study, we generated and analyzed a dataset of both microsatellite and mtDNA loci from B. dorsalis samples collected from six offshore island populations in South China to estimate the genetic divergence and dispersal ability of these flies. Our results may also suggest possible strategies for the control of this species based on their dispersal patterns.
DNA from individual specimens was extracted using the TIANamp Genomic DNA Kit (Qiagen, Hilden, Germany). In total, DNA was extracted from 20 flies from each island. Eight polymorphic microsatellite loci (MS4, MS6, MS12A, 4.3A, 6.8A, 4.6A, Ccmic32 and Bo-D48) were analyzed, all of which were previously developed [14, 21–23] (see Additional file 1: Table S1). Primers were synthesized by Sangon Biotech Co. Ltd (Shanghai, China). PCR reactions were performed in 25 µL volumes containing 2 µL of DNA extract (40 ng/µL), 2.5 µL of 10× PCR buffer (containing Mg2+), 2 µL of a dNTP mixture (10 mmol/L), 1.5 µL of each primer (10 mmol/L), 0.2 µL of Taq DNA polymerase (5 U/µL) and 15.3 µL of ddH2O. Cycling conditions for amplification were as follows: 3 min at 94 °C followed by 35 cycles of 30 s at 94 °C, 30 s at 48 °C, 45 s at 72 °C, and an extension of 5 min at 72 °C.
The PCR-amplified products were separated using a DNF-900 High Sensitivity Large Fragment Analysis Kit in a fragment analyzer (Advanced Analytical Technologies, Inc., USA). Alleles were scored using PROSize 2.0 (Advanced Analytical Technologies, Inc., USA).
A 759-bp fragment of the cytochrome oxidase subunit I (COI) gene was amplified using primers (COI-F: CAACATTTATTTTGATTTTTTGG; COI-R: TCCATTGCACTAATCTGCCATATTA; synthesized by Sangon Biotech Co., Ltd (Shanghai, China)) and following the method described by Tang et al. [24, 25]. The PCR reaction volumes were the same as the microsatellite DNA amplification volumes, except the annealing temperature duration was 45 s.
Deviations from Hardy–Weinberg proportions and linkage disequilibrium were calculated using Genepop version 5.0 (http://www.wbiomed.curtin.edu.au/Genepop/Genepop_op1.html) . We estimated null allele frequencies for each locus using the expectation maximization (EM) algorithm method in the FREENA software (http://www.montepllier.inra.fr/URLB) . Population genetic parameters included the observed number of alleles (NA), the effective number of alleles (N E ), Shannon’s information index (I), observed heterozygosity (H o ), expected heterozygosity (H E ), the percentage of polymorphic loci (P), Nei’s genetic diversity (Nei’s), gene flow (N em ), Nei’s original measures of genetic identity (I) and genetic distance (D) for paired populations. These parameters were then analyzed using Popgene Version 1.3.1 . The population genetic structure was analyzed with STRUCTURE v.5.0 , which uses a Bayesian Markov Chain Monte Carlo (MCMC) method to analyze the numbers of K genetic clusters. We ran K values from 1 to 6 using an admixture model and an allele frequencies correlated model. The first 100,000 repetitions were discarded as burn-in, and 10,000 MCMC repetitions were then run. Multiple runs are used to check consistency between runs. Then, the results after running STRUCTURE were then uploaded to the Structure Harvester website (http://www.taylor0.biology.ucla.edu/structureHarvester/#). Structure Harvester was used to calculate six iterations per K-value and obtained the best K-value using the ΔK method of Evanno. . Nei’s genetic distance (D) for the six island populations was reconstructed for the population phylogenetic tree using the unweighted pair-group method with arithmetic mean (UPGMA) in the tools for population genetic analyses (TFPGA) program . We used a multidimensional factorial correspondence analysis (FCA) within the GENETIX program to identify clusters of individuals with similar genotypes based on allele frequencies and genotype .
Isolation by distance was measured using the TFPGA software. Matrix correlations using pairwise genetic distance versus geographic distance were estimated and significance was determined using the Mantel test.
Analysis of molecular variance (AMOVA) was performed with ARLEQUIN ver. 3.5 to partition the genetic variance within and among populations based on the numbers of different alleles (F ST ) and the sum of squared size differences .
Published COI sequences of B. dorsalis from mainland China populations and mainland populations from other countries (Cambodia, Mandalay, Laos and America) were downloaded from NCBI [16, 34] to characterize the genetic relationships between these mainland populations and the studied six island populations in South China. We cut and joined the amplification sequences using one pair of primers with DNAStar and manually corrected any obvious errors. Next, we aligned all sequences using MEGA6 software . A population phylogenetic tree was reconstructed using the UPGMA method in MEGA6 based on genetic distances. Haplotype diversity (Hd), the average number of differences (k), and nucleotide diversity (π) were calculated for each population using DnaSP 5.0 [20, 36].
The fixation index (pairwise F ST ) and gene flow (N em ) were estimated using ARLEQUIN ver. 3.5 . Generally, gene exchange leading to low genetic differentiation between populations occurs when N em > 4 . The significance level was assessed using 1000 permutations. The index was interpreted as follows: a low degree of genetic differentiation (0 ≤ F ST < 0.05); a medium degree of genetic differentiation (0.05 ≤ F ST < 0.15); a high degree of genetic differentiation (0.15 ≤ F ST ≤ 0.25); and a very high degree of genetic differentiation (F ST > 0.25) [38, 39].
Hardy–Weinberg and genetic diversity
All loci and samples significantly deviated from Hardy–Weinberg proportions (P < 0.01). The average values of F IS for each locus in each population ranged from 0.0365 to 0.8353, which indicated that there were deficiencies in heterozygotes. Over all populations and loci, the locus MS12A for the DAD population had the highest null allele frequency (N a = 0.41419). The null allele frequency of the other loci and populations ranged from 0 to 0.37780. The presence of null alleles was the primary contributor to the deficiencies in heterozygotes and the departure from the Hardy–Weinberg equilibrium. There was no significant linkage disequilibrium between pairs of loci for all populations (P > 0.05). Therefore, the eight loci were inherited independently.
Indices of genetic diversity of Bactrocera dorsalis populations from six islands based on microsatellite data
5.7500 ± 2.1876
3.7076 ± 1.8766
1.3450 ± 0.4904
0.2970 ± 0.2371
0.6740 ± 0.1850
0.6571 ± 0.1804
8.0000 ± 2.3905
4.4046 ± 1.8152
1.6184 ± 0.4118
0.4313 ± 0.2359
0.7532 ± 0.1131
0.7344 ± 0.1103
7.1250 ± 2.1671
3.4046 ± 1.6079
1.3886 ± 0.4488
0.3355 ± 0.2442
0.2442 ± 0.1769
0.1769 ± 0.1725
7.6250 ± 2.8754
3.8322 ± 1.7547
1.4875 ± 0.5302
0.3500 ± 0.3218
0.6697 ± 0.1983
0.6869 ± 0.1933
6.8750 ± 2.7999
3.4946 ± 1.5458
1.4143 ± 0.4487
0.2812 ± 0.2103
0.6796 ± 0.1494
0.6627 ± 0.1457
7.8750 ± 2.7999
4.4092 ± 1.9838
1.6249 ± 0.4329
0.4016 ± 0.2906
0.7546 ± 0.1072
0.7356 ± 0.1046
7.3214 ± 0.8192
3.8755 ± 0.4386
1.4798 ± 0.1192
0.3494 ± 0.0583
0.6292 ± 0.1926
0.6089 ± 0.2143
Population genetic structure
AMOVA results based on microsatellite genotypes
Source of variation
Percentage of variation
FST = 0.04161**
Among individuals within populations
FIS = 0.50715**
FIT = 0.52765**
FCT = 0.04462**
Among populations within groups
FSC = 0.00863
Among individuals within populations
FIS = 0.50631**
FIT = 0.53241**
Nei’s original measures of genetic identity (above diagonal) and genetic distance (below diagonal)
Haplotype diversity (Hd), average number of differences (k) and nucleotide diversity (π) for each population
H1(1), H2(1), H3(5), H4(2), H5(5), H6(1), H7(1), H8(1), H9(1), H10(1), H11(1)
0.889 ± 0.049
H3(3), H5(1), H10(1), H13(1), H19(4), H23(2), H57(1), H105(1), H142(1), H143(1), H144(1), H145(1), H146(1), H147(1)
0.947 ± 0.034
H3(2), H5(4), H8(1), H13(1), H19(1), H25(1), H74(1), H77(1), H102(1), H103(1), H104(1), H105(1), H106(1), H107(1), H108(1), H109(1)
0.963 ± 0.033
H3(2), H5(2), H8(1), H11(1), H13(1), H19(2), H37(1), H77(2), H105(1), H133(1), H134(1), H135(1), H136(1),H137(1), H138(1), H139(1)
0.979 ± 0.021
H3(2), H5(2), H44(1), H62(1), H77(1), H107(1), H129(1), H140(1), H142(1), H148(1), H149(1), H150(1), H151(1), H152(1), H153(2), H154(1), H155(1)
0.984 ± 0.020
H5(2), H8(1), H25(1), H52(1), H77(1), H119(1), H120(1), H121(1), H122(1), H123(1), H124(1), H125(1), H126(1), H127(1), H128(1), H129(1), H130(1), H131(1), H132(1)
0.995 ± 0.018
0.9665 ± 0.008
F ST values (below diagonal) and N em (above diagonal) between the six island populations of Bactrocera dorsalis
Null alleles, geographic factors and likely founder effects among populations caused heterozygote deficiencies and emerged as significant departures from the Hardy–Weinberg equilibrium . In this study, the frequency of null alleles in an isolated microsatellite of Bactrocera dorsalis was high and exceeded 0.4 only among locus 6.4. We checked the heterozygotes and found significant heterozygote deficiencies, except for locus MS12A, in several island populations.
Population genetic diversity is a product of evolutionary change over several generations . Thus, population genetic diversity reflects the ability of populations to adapt to local environments . The percentage of polymorphism, the level of heterozygosity, the number of alleles, nucleotide diversity and haplotype diversity were used to estimate population genetic diversity. Our research showed that populations of the oriental fruit fly from islands in South China had high levels of genetic diversity and were closely related to populations from mainland China based on COI data analysis. Using the mtDNA COI from five Yunnan province populations of B. dorsalis, Shi et al.  showed that the average haplotype diversity and nucleotide diversity of the five populations were 0.9786 and 0.9038 %, respectively. Thus, the five Yunnan populations had higher levels of genetic diversity than did the island populations examined in this study. Analysis of the genetic diversity of B. dorsalis from China, Laos, and Thailand with microsatellite markers showed that the average Nei’s genetic diversity and Shannon’s information index were 0.6464 and 0.7870, respectively. The average percentage of polymorphic loci in all populations was 94.45 %. This finding suggests that mainland South China has high levels of population genetic diversity .
However, many studies have shown that the gene flow of mammals, birds and other species from islands is blocked or reduced due to isolation, which results in lower levels of genetic diversity of island populations compared to mainland populations [8, 14, 45]. Isolated islands have discrete boundaries that are generally thought to reduce migration between populations on islands . Long-term bottlenecks, founder effects, genetic drift and inbreeding may all reduce genetic diversity [46, 47]. Mitochondrial and microsatellite data have been previously used to estimate the genetic structure of Apis mellifera populations from the Canary Islands . The results showed that there was a lower level of genetic variation based on the average number of alleles and heterozygosity in populations on the Canary Islands than mainland populations in Iberia and Morocco . Boessenkool et al. reported that natural populations of robins on the Breaksea and Nukuwaiata Islands have lower levels of genetic diversity than larger mainland populations. In addition, some alleles have been lost, which is a result of bottlenecks and isolation from the mainland . Jensen et al.  reported that the genetic diversity of house sparrows along a coastal latitudinal gradient from middle to Northern Norway tended to decrease in island populations compared to mainland populations. However, in our study, the results showed that the genetic diversity of B. dorsalis from island populations was high. The geographic isolation imposed by the sea may not hinder gene exchange. There are many reasons for high genetic diversity among island populations, including the dispersal ability of this species and both environmental and anthropogenic activities that cause frequent contact with mainland populations. The six islands have been isolated for a long period of time, and most of the islands originated from the accumulation of volcanic matter. As the islands developed, commercial activities on the islands increased, which has incidentally promoted the spread of species throughout the islands. B. dorsalis can use most commercial plants and fruits as host plants, and the flies can disperse between the mainland and islands by relying on the wind. According to the phylogenetic tree analysis, there were two main routes of invasion from inland China. One route was from Southeast China, and the other route was from Southwest China . However, whether B. dorsalis regularly makes migratory flights from the nearby mainland to these islands is unclear.
For COI data, the value of N em among populations all exceeded 4 (N em > 4), which showed that among populations there was a high levels of gene flow and a low and medium level of differentiation among some paired island populations. These results were consistent with the results from our microsatellite data that shows a low Nei’s standard genetic distance (D) and a high genetic identity (I). Pairwise F ST values indicated that there was a medium degree of population differentiation between Weizhou Island and Naozhou Island. The microsatellite data indicated that population genetic variation was mostly partitioned among populations. We divided the six island populations into three groups, which were supported by the FCA. Our results were also similar to the levels of population differentiation between mainland China populations. Wan et al.  analyzed the genetic diversity of B. dorsalis from six populations in Qiongqing using eight microsatellite loci. The results showed that the Qiongqing populations had low levels of population differentiation. Yao et al.  also examined the genetic relationships among populations from Fujian, Hainan, Guangdong, Yunnan and Sichuan provinces. Except for the Fujian population, all populations had low levels of population differentiation because of the proximity between provinces. Using mtDNA from 25 populations on the China mainland, flies from mainland China had lower levels of genetic divergence than those from Thailand, Japanese and American populations . Li et al.  reported that the average Nei’s standard genetic distance was 0.8049 and 0.9397, for South China and Southeast Asia populations, respectively. Li et al.  also reported that F ST was 0.25 between South China and Southeast Asia populations, which suggested that these populations had genetic divergence due to geographic isolation. Our study suggested that island B. dorsalis populations had high levels of genetic diversity and a low or medium level of differentiation among some paired island populations and the genetic distance of pairing populations had no correlation with geographic isolation (R2 = 0.0027, P = 0.855). There are several possible explanations for this pattern. First, Tephritid fruit flies are capable fliers that can fly more than 25 km under windless conditions [52, 53]. The distance between each island and the mainland ranged from 8.2 to 38.4 km. The distance between islands ranged from 32.4 to 519.7 km. Hebao Island, Wailingding Island and Dong’ao Island are located in Zhuhai province and are closer in proximity than are Weizhou Island and Naozhou Island. Thus, Hebao Island populations were barely differentiated from the Wailingding Island, Shangchuan Island and Naozhou Island populations. However, there was a medium degree of genetic differentiation between Dong’ao Island and Weizhou Island populations and between Dong’ao Island and Naozhou Island populations. Bactrocera dorsalis can migrate to and from islands on its own or with the aid of typhoons . The isolation of islands by water limits gene flow. However, in our study, the islands did not appear to restrict gene flow: the population genetic diversity and differentiation on the islands were similar to values observed for mainland populations. Second, based on the results of our study, commercial activities between the island and the mainland have affected the gene pool and diversity of species. These six islands all have suitable host plants for B. dorsalis, such as banana, pawpaw and others. Many tropical fruits, such as guava, carambola, and mango, are imported and exported from the mainland to the islands. Bactrocera dorsalis can lay eggs on these fruits, and the infested fruits can subsequently be transported to other areas. Therefore, gene flow is likely frequent between islands. The low levels of genetic diversity are also likely caused by the short distance between the islands and the mainland.
This study showed that island isolation may not significantly influence the genetic diversity of tephritid fruit flies. We have found that offshore island populations of B. dorsalis have relatively high levels of genetic diversity, whereas populations from offshore islands exhibited low genetic differentiation.
YJX conceived and designed the experiments. CYY performed the experiments. YJX, LZ, CYY and CYZ analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
We wish to thank G. Q. Wang, C. Cui, and Y. Q. Huang from South China Agricultural University for their help in collecting the material used in this study.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset supporting the results of this article is available at http://dx.doi.org/10.5061/dryad.666sn.
This study was supported by the Science Foundation for Excellent Youth Scholars of Guangdong Province (No. Yq2013031). The funders played no role in the design of the study, data collection and analysis or the decision to publish and prepare the manuscript.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Booy G, Hendriks R, Smulders M, Groenendael JV, Vosman B. Genetic diversity and the survival of populations. Plant Biol. 2000;2(4):379–95.View ArticleGoogle Scholar
- Ellstrand NC, Elam DR. Population genetic consequences of small population size: implications for plant conservation. Annu Rev Ecol Syst. 1993;24:217–42.Google Scholar
- Jensen H, Moe R, Hagen IJ, Holand AM, Kekkonen J, Tufto J, Sæther BE. Genetic variation and structure of house sparrow populations: is there an island effect? Mol Ecol. 2013;22(7):1792–805.PubMedView ArticleGoogle Scholar
- Petit RJ, Excoffier L. Gene flow and species delimitation. Trends Ecol Evol. 2009;24(7):386–93.PubMedView ArticleGoogle Scholar
- Gillespie RG, Roderick GK. Arthropods on islands: colonization, speciation, and conservation. Annu Rev Entomol. 2002;47(1):595–632.PubMedView ArticleGoogle Scholar
- Martinkova N, Barnett R, Cucchi T, Struchen R, Pascal M, Pascal M, Fischer MC, Higham T, Brace S, Ho SY, et al. Divergent evolutionary processes associated with colonization of offshore islands. Mol Ecol. 2013;22(20):5205–20.PubMedPubMed CentralView ArticleGoogle Scholar
- De la Rúa P, Galián J, Serrano J, Moritz R. Genetic structure and distinctness of Apis mellifera L. populations from the Canary Islands. Mol Ecol. 2001;10(7):1733–42.PubMedView ArticleGoogle Scholar
- Frankham R. Do island populations have less genetic variation than mainland populations? Heredity. 1997;78(3):311–27.PubMedView ArticleGoogle Scholar
- Miller HC, Allendorf F, Daugherty CH. Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.). Mol Ecol. 2010;19(18):3894–908.PubMedView ArticleGoogle Scholar
- Wolf AT, Harrison SP. Effects of habitat size and patch isolation on reproductive success of the serpentine morning glory. Conserv Biol. 2001;15(1):111–21.Google Scholar
- Boessenkool S, Taylor SS, Tepolt CK, Komdeur J, Jamieson IG. Large mainland populations of South Island robins retain greater genetic diversity than offshore island refuges. Conserv Genet. 2006;8(3):705–14.View ArticleGoogle Scholar
- Brekke P, Bennett PM, Santure AW, Ewen JG. High genetic diversity in the remnant island population of hihi and the genetic consequences of re-introduction. Mol Ecol. 2011;20(1):29–45.PubMedView ArticleGoogle Scholar
- Clarke A, Armstrong K, Carmichael A, Milne J, Raghu S, Roderick G, Yeates D. Invasive phytophagous pests arising through a recent tropical evolutionary radiation: the Bactrocera dorsalis complex of fruit flies. Annu Rev Entomol. 2005;50:293–319.PubMedView ArticleGoogle Scholar
- Dai S, Lin C, Chang C. Polymorphic microsatellite DNA markers from the oriental fruit fly Bactrocera dorsalis (Hendel). Mol Ecol Notes. 2004;4(4):629–31.View ArticleGoogle Scholar
- Huang S, Han R. Advance in the research on the quarantine pest Bactrocera dorsalis. Chin Bull Entomol. 2005;5:000.Google Scholar
- Li Y, Wu Y, Chen H, Wu J, Li Z. Population structure and colonization of Bactrocera dorsalis (Diptera: Tephritidae) in China, inferred from mtDNA COI sequences. J Appl Entomol. 2012;136(4):241–51.View ArticleGoogle Scholar
- Aketarawong N, Bonizzoni M, Thanaphum S, Gomulski L, Gasperi G, Malacrida A, Gugliemino C. Inferences on the population structure and colonization process of the invasive oriental fruit fly, Bactrocera dorsalis (Hendel). Mol Ecol. 2007;16(17):3522–32.PubMedView ArticleGoogle Scholar
- Drew RA, Hancock DL. The Bactrocera dorsalis complex of fruit flies (Diptera: Tephritidae: Dacinae) in Asia. Bull Entomol Res Suppl Ser. 1994;2:1–68.View ArticleGoogle Scholar
- Shi W, Kerdelhué C, Ye H. Genetic structure and inferences on potential source areas for Bactrocera dorsalis (Hendel) based on mitochondrial and microsatellite markers. PLoS One. 2012;7(5):e37083.PubMedPubMed CentralView ArticleGoogle Scholar
- Wan X, Liu Y, Zhang B. Invasion history of the oriental fruit fly, Bactrocera dorsalis, in the Pacific-Asia region: two main invasion routes. PLoS One. 2012;7(5):e36176.PubMedPubMed CentralView ArticleGoogle Scholar
- Augustinos A, Stratikopoulos E, Zacharopoulou A, Mathiopoulos K. Polymorphic microsatellite markers in the olive fly, Bactrocera oleae. Mol Ecol Notes. 2002;2(3):278–80.View ArticleGoogle Scholar
- Baliraine F, Bonizzoni M, Osir E, Lux S, Mulaa F, Zheng L, Gomulski L, Gasperi G, Malacrida A. Comparative analysis of microsatellite loci in four fruit fly species of the genus Ceratitis (Diptera: Tephritidae). Bull Entomol Res. 2003;93(01):1–10.PubMedView ArticleGoogle Scholar
- Wang Y, Yu H, Raphael K, Gilchrist A. Genetic delineation of sibling species of the pest fruit fly Bactocera (Diptera: Tephritidae) using microsatellites. Bull Entomol Res. 2003;93(04):351–60.PubMedView ArticleGoogle Scholar
- Mun J, Bohonak AJ, Roderick GK. Population structure of the pumpkin fruit fly Bactrocera depressa (Tephritidae) in Korea and Japan: pliocene allopatry or recent invasion? Mol Ecol. 2003;12(11):2941–51.PubMedView ArticleGoogle Scholar
- Tang K. Molecular indentification of Bactrocera (Diptera: Tephritidae) larvae by using microsatellite markers and population genetic structure of Bactrocera dorsalis in China base on COI data Master. Guilin: Guangxi Normal University; 2006.Google Scholar
- Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86(3):248–9.Google Scholar
- Chapuis M-P, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24(3):621–31.PubMedView ArticleGoogle Scholar
- Yeh F, Yang R, Boyle T. POPGENE version 1.31 quick user guide. Canada: University of Alberta and Center for International Forestry Research; 1999.Google Scholar
- Vorsino AE, Wieczorek AM, Wright MG, Messing RH. Genetic analysis of an introduced biological control agent reveals temporal and geographic change, with little evidence of a host mediated shift. Biol Control. 2014;77:41–50.View ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.PubMedView ArticleGoogle Scholar
- Miller MP. Tools for population genetic analyses (TFPGA) 1.3: a windows program for the analysis of allozyme and molecular population genetic data. 1997;4:157. (Computer software distributed by author).Google Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire génome, populations, interactions, CNRS UMR. 1996;5000:1996–2004.Google Scholar
- Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010;10(3):564–7.View ArticleGoogle Scholar
- Wu G, Li Y, Li Z, Liu L. The genetic relation of the population of Bactrocera dorsalis (Diptera: Tephritidae) in Shanghai District based on mtDNA CO I gene. J China Agric Univ. 2012;2:94–101.Google Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.Google Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.PubMedView ArticleGoogle Scholar
- Wright S. Variability within and among natural populations. Evolution and the genetics of populations: a treatise in four volumes. Vol 4. Chicago: University of Chicago Press; 1978.Google Scholar
- Chakraborty R, Danker-Hopfe H. 7 Analysis of population structure: a comparative study of different estimators of wright’s fixation indices. Handbook Stat. 1991;8:203–54.View ArticleGoogle Scholar
- Rousset F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997;145(4):1219–28.PubMedPubMed CentralGoogle Scholar
- Wen YFUK, Han WJ, Ueno S, Xie WD, Xu GB, Tsumura Y. Null alleles in microsatellite markers. Biodivers Sci. 2013;21(1):117–26.View ArticleGoogle Scholar
- Rieseberg LH, Soltis D. Phylogenetic consequences of cytoplasmic gene flow in plants: evolutionary trends in Plants; 1991.Google Scholar
- Pogson GH, Mesa KA, Boutilier RG. Genetic population structure and gene flow in the Atlantic cod Gadus morhua: a comparison of allozyme and nuclear RFLP loci. Genetics. 1995;139(1):375–85.PubMedPubMed CentralGoogle Scholar
- Shi W, Ye H. Genetic differentiation in five geographic populations of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) in Yunnan Province. Acta Entomol Sin. 2003;47(3):384–8.Google Scholar
- Li W, Yang L, Tang K, Zeng L, Liang G. Microsatellite polymorphism of Bactrocera dorsalis (Hendel) populations in China. Acta Entomol Sin. 2007;50(12):1255–62.Google Scholar
- Francisco-Ortega J, Santos-Guerra A, Kim S-C, Crawford DJ. Plant genetic diversity in the Canary Islands: a conservation perspective. Am J Bot. 2000;87(7):909–19.PubMedView ArticleGoogle Scholar
- Green EC, Tremetsberger K, Jiménez A, Gómez-González S, Stuessy TF, Baeza CM, López PG. Genetic diversity of pioneer populations: the case of Nassauvia argentea (Asteraceae: Mutisieae) on Volcán Lonquimay, Chile. Plant Syst Evol. 2012;298(1):109–19.View ArticleGoogle Scholar
- Nei M, Maruyama T, Chakraborty R. The bottleneck effect and genetic variability in populations. Evolution. 1975;29(1):1–10.Google Scholar
- Franck P, Garnery L, Solignac M, Cornuet J-M. The origin of west European subspecies of honeybees (Apis mellifera): new insights from microsatellite and mitochondrial data. Evolution. 1998;52(4):1119–34.Google Scholar
- Boessenkool S, Taylor SS, Tepolt CK, Komdeur J, Jamieson IG. Large mainland populations of South Island robins retain greater genetic diversity than offshore island refuges. Conserv Genet. 2007;8(3):705–14.View ArticleGoogle Scholar
- Wang M, Xu L, Zhang R, Zhang G, Yu D. Population genetic differentiation analysis of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae) based on mtDNA COI gene. Acta Entomol Sin. 2014;57(12):1424–38.Google Scholar
- Yao TSHJ, Li HJ, Ran C, Liu HQ, Lei HD, Xiao T. Microsatellite markers for genetic polymorphism in five geographic population of the Bactrocera dorsalis (Hendel). Chin Agric Sci Bull. 2010;26(8):234–6.Google Scholar
- Fletcher B. The biology of dacine fruit flies. Annu Rev Entomol. 1987;32:115–44.View ArticleGoogle Scholar
- Koyama J, Kakinohana H, Miyatake T. Eradication of the melon fly, Bactrocera cucurbitae, in Japan: importance of behavior, ecology, genetics, and evolution. Ann Rev Entomol. 2004;49(1):331–49.View ArticleGoogle Scholar