Research Article
Print
Research Article
Relevant genetic differentiation among Brazilian populations of Anastrepha fraterculus (Diptera, Tephritidae)
expand article infoMosè Manni, Katia Manuela Lima§, Carmela Rosalba Guglielmino, Silvia Beatriz Lanzavecchia|, Marianela Juri|, Teresa Vera, Jorge Cladera#, Francesca Scolari, Ludvik Gomulski, Mariangela Bonnizoni, Giuliano Gasperi, Janisete G Silva¤, Anna Rodolfa Malacrida
‡ University of Pavia, Pavia, Italy
§ Universidade Estadual de Santa Cruz, Ilhéus, Brazil
| Instituto Nacional de Tecnología Agropecuaria (INTA), Castelar, Argentina
¶ Universidad Nacional de Tucumán (UNT), San Miguel de Tucumán, Argentina
# Instituto Nacional de Tecnología Agropecuaria, Castelar, Argentina
¤ Departamento de Ciências Biológicas, Universidade Estadual de Santa Cruz, Bahia, ilheus, Brazil
Open Access

Abstract

We used a population genetic approach to detect the presence of genetic diversity among six populations of A. fraterculus across Brazil. To this aim, we used Simple Sequence Repeat (SSR) markers, which may capture the presence of differentiative processes across the genome in distinct populations. Spatial analyses of molecular variance were used to identify groups of populations that are both genetically and geographically homogeneous while also being maximally differentiated from each other. The spatial analysis of genetic diversity indicates that the levels of diversity among the six populations vary significantly on an eco-geographical basis. Particularly, altitude seems to represent a differentiating adaptation, as the main genetic differentiation is detected between the two populations present at higher altitudes and the other four populations at sea level. The data, together with the outcomes from different cluster analyses, identify a genetic diversity pattern that overlaps with the distribution of the known morphotypes in the Brazilian area.

Keywords

Anastrepha fraterculus, microsatellites, population genetic differentiation, morphotypes

Introduction

The South American fruit fly Anastrepha fraterculus Wiedemann (Diptera: Tephritidae) belongs to the fraterculus group, which comprises a total of 34 formally described species (Norrbom et al. 2012) that can be distinguished only by minor morphological characters. The major diagnostic character for species within this group is the aculeus apex, which shows a large degree of intraspecific variation (Araujo and Zucchi 2006) as a result of genetic and environmental factors (Aluja 1994, Smith-Caldas et al. 2001).

The nominal species A. fraterculus is widely distributed from the Rio Grande Valley in northern Mexico to central Argentina, infesting over 100 hosts (Norrbom 2004), thus being a species of major economic importance in Brazil and other countries in South America (Steck 1991, Steck 1999, Zucchi 2008). Its distribution in South America was thought to be in two broad and unconnected bands along the east coast and the western and northern edges of the continent with a hiatus comprising the Amazon basin (Steck 1999, Malavasi et al. 2000). Yet, recent collections have recorded this species in the Brazilian Amazon in the states of Amapá, Pará, Tocantins, and Maranhão infesting 10 different hosts (Zucchi et al. 2011).

A. fraterculus has long been reported to show extensive morphological variation along its geographic distribution (Lima 1934, Stone 1942, Steck 1991, 1999). In his taxonomic review of the genus Anastrepha, Stone (1942) stated “it is possible that it will eventually be found to represent a complex of species rather than a single one”. Since then, a good deal of research has documented and concluded that the nominal species A. fraterculus actually comprises an unresolved complex of cryptic species. Evidence comes from studies on morphological variation (Steck 1999, Hernández-Ortíz et al. 2012, for review), multivariate morphometric analyses (Hernández-Ortíz et al. 2004, 2012), differences in host use (Steck 1999, for review), behaviour (Yamada and Selivon 2001, Selivon et al. 2005, Vera et al. 2006, Vaníčková et al. 2015), the presence and degree of reproductive isolation (Devescovi et al. 2014 and references therein), and genetic analyses (Silva and Barr 2008, for review). However, the actual number of species within the A. fraterculus complex and their distribution is yet to be elucidated.

Genetic studies performed on A. fraterculus populations so far have revealed the following putative biological entities based on geography: an Andean lineage (Steck 1991), a Mexican species (Steck 1991, Smith-Caldas et al. 2001, Barr et al. 2005 and the morphometric studies of Hernández-Ortíz et al. 2004, 2012), a Guatemalan lineage (Smith-Caldas et al. 2001), a second Venezuelan species (Smith-Caldas et al. 2001), a Peruvian lineage (Steck 1991), and three Brazilian species (Morgante et al. 1980, Smith-Caldas et al. 2001, Selivon et al. 2005, and Ruiz et al. 2007a, 2007b) (Silva and Barr 2008, for review). In addition, morphometric studies revealed that Mexican populations formed a single morphotype, which was distinct from South American populations (Hernández-Ortíz et al. 2004). Hernández-Ortíz et al. (2004) also identified an Andean morphotype (one population from Colombia) and a Brazilian morphotype comprising two populations from Brazil (states of São Paulo and Santa Catarina) and a population from Argentina (Tucumán). A study using isozymes, karyotypes, morphometry, and crossings on populations of A. fraterculus from Brazil recognized two species, Anastrepha sp.1 aff. fraterculus and Anastrepha sp.2 aff. fraterculus (Selivon et al. 2005). More recently, a multivariate morphometric analysis comprising 32 A. fraterculus populations identified seven distinct morphotypes: a Mexican morphotype, a Venezuelan morphotype, an Andean morphotype, a Peruvian morphotype, and three Brazilian morphotypes (Brazilian 1, Brazilian 2, and Brazilian 3) (Hernández-Ortíz et al. 2012). Although Brazilian populations of the nominal species A. fraterculus most likely comprise at least three morphotypes, published studies are based mostly on samples from the south-eastern region, while very few populations from other regions of the country have been examined.

Previous genetic studies using DNA sequencing of mitochondrial genes from different populations suggested that both the fraterculus group and the A. fraterculus complex have a recent evolutionary history, and thus molecular markers with a higher power of resolution were required to help understand the specific/subspecific differentiation within and between populations of this group (McPheron et al. 1999, Smith-Caldas et al. 2001). Microsatellites are a class of highly polymorphic molecular markers widely distributed in the genome of eukaryotes that can be useful to clarify such patterns of gene flow and to identify the spatial locations of genetic discontinuities (population boundaries) in studies of species complexes (Chambers and MacAvoy 2000, Barbará et al. 2007, Aketarawong et al. 2014). Within the Tephritidae, microsatellites have been successfully developed and applied for several Bactrocera species (Shearman et al. 2006; Aketarawong et al. 2006, 2007, Augustinos et al., 2008, Virgilio et al. 2010; Drew et al. 2011), for Rhagoletis cerasi L. (Augustinos et al. 2011), for a few Ceratitis species (Bonizzoni et al. 2001, 2004, Meixner et al. 2002, Baliraine et al. 2003, 2004, Silva et al. 2003, Delatte et al. 2013, Virgilio et al. 2013), for Anastrepha suspensa (Loew) (Fritz and Schable 2004, Boykin et al. 2010), and for Anastrepha obliqua (Macquart) (Islam et al. 2011). Microsatellites have only recently been isolated in A. fraterculus (Lanzavecchia et al. 2014) and have proven useful for the analysis of population dynamics and differentiation across the distribution range of this polymorphic species.

From an applied perspective, the correct identification of populations and species is an important step in the implementation of biologically-based control methods such as the Sterile Insect Technique (SIT) against this fruit fly complex (Silva and Barr 2008), which represents a serious constraint for fruit production in South America and a hindrance to the export of fresh fruit from regions where it occurs. Currently, control measures for this pest species rely solely on the use of insecticide cover or bait sprays. Therefore, there is demand for the development of the SIT as it would benefit South American countries such as Argentina, Brazil, and Peru (Cladera et al. 2014). However, the application and efficiency of species-specific control methods such as the SIT are critically dependent on the correct identification of the target pest populations and the understanding of the spatial distribution of the pest species, thus the correct delimitation of species within the A. fraterculus complex is paramount.

This paper is centred on the assessment of genetic diversity among A. fraterculus populations from distinct geographic regions across Brazil, most likely belonging to at least three distinct morphotypes (Silva et al. unpubl. data). For this purpose we used highly informative SSR markers, which may capture the presence of eventual differentiative processes across the genome in different populations.

Methods

Six populations from three regions across Brazil were sampled from 2007 to 2013 (Table 1, Figure 1). In the Northeastern region, the populations from Monte Alegre (State of Rio Grande do Norte) and from Una and Porto Seguro (State of Bahia) were sampled. In the Southeastern region, samples from São Mateus (State of Espírito Santo) and from Campos do Jordão (State of São Paulo) were examined. In the Southern region, the population from Vacaria (State of Rio Grande do Sul) was sampled. For each locality, flies emerging from fruits collected from different trees were considered. Adult females were identified as A. fraterculus by Dr. Elton L. Araujo (UFERSA), Dr. Keiko Uramoto (USP), Dr. Miguel Francisco Souza Filho (Instituto Biológico) and Dr. Roberto A. Zucchi (USP) using the aculeus shape following Zucchi (2000) (Table 1). Voucher specimens were deposited at the insect collection of the Escola Superior de Agricultura “Luiz de Queiroz”, USP, Piracicaba, SP, and at the Universidade Estadual de Santa Cruz, Ilhéus, BA, Brazil. According to the classification of Hernández-Ortíz et al. (2012), the flies collected in Una (BA) can be classified as Brazilian morphotype 3, while those from Vacaria (RS) and Campos do Jordão (SP) as Brazilian 1 (Silva et al. unpublished data). For the samples from Monte Alegre, Porto Seguro and São Mateus, no clear-cut information is available to assign them to a specific morphotype.

Figure 1.

Map of the collected samples.

Figure 2.

Two-dimensional plot of Principal Coordinate Analysis (PCoA) based on similarity matrix derived from Anastrepha fraterculus microsatellites data. The Morphotype classification (Hernández-Ortiz et al. 2012), relative to each sample is also reported.

Field collected samples of Anastrepha fraterculus Brazilian populations.

States Sample site Morphotype* Host Coordinate Elevation
Rio Grande do Norte (RN) Monte Alegre ? Guava 6.0678W; 35.3322S 51.816m
Bahia (BA) Una 3 Guava 15.2933W; 39.0753S 27.737m
Bahia (BA) Porto Seguro ? Guava 16.4497W; 39.0647S 48.768m
Espírito Santo (ES) São Mateus ? Araçá 18.7161W; 39.8589S 35.966m
São Paulo (SP) Campos do Jordão 1 Raspberry 22.7394W; 45.5914S 1627.9m
Rio Grande do Sul (RS) Vacaria 1 Guava 28.5122W; 50.9339S 970.79m

Microsatellite analysis

A total of 171 A. fraterculus individuals collected from the above mentioned populations were assessed for their SSR variability. DNA was extracted from three legs of each single fly using the “DNeasy Blood & Tissue” kit (Qiagen, Valencia, CA) following the standard DNeasy protocol. DNA samples were screened using the following ten microsatellite loci: AfD4, AfD105, AfA7, AfA112, AfA115, AfA120, AfA122, AfA117, AfA10, and AfC103 (Lanzavecchia et al. 2014). Allele scoring was performed using an automated ABI PRISM 310 Genetic Analyser (Applied Biosystem) following Aketarawong et al. (2006).

Data analysis

The mean number of alleles (na) and mean null allele frequency (An) (non-amplifying alleles due to changes in the primer binding regions), expected and observed heterozygosity were estimated using GENEPOP version 4.0.7 (Raymond and Rousset 1995) for each population. Deviation from the Hardy-Weinberg equilibrium and linkage disequilibrium, together with their critical levels after the sequential Bonferroni test (Rice 1989), were tested using GENEPOP version 4.0.7 (Raymond and Rousset 1995). The allelic Polymorphic Information Content (PIC) was derived using CERVUS (Kalinowski et al. 2007).

Microsatellite Analyzer (MSA) (Dieringer and Schlötterer 2003) was applied to estimate the pairwise Fst values among populations (Weir and Cockerham 1984). The statistical significance of each Fst value was assessed by comparing the observed values with the values obtained in 10,000 matrix permutations. Spatial analyses of molecular variance were investigated using SAMOVA 2.0 (Dupanloup et al. 2002). This approach identifies groups of populations that are genetically homogeneous and maximally differentiated from each other without the constraint of being geographically close. The method requires the a priori definition of the number of groups (K) of populations that exist and generates F-statistics (FSC, FST and FCT) using an AMOVA approach. Different numbers of groups (K) were tested, and a simulated annealing procedure permitted the identification of the composition of each of the K groups that maximizes the FCT index (proportion of total genetic variance due to differences between groups). The program was run for two to five groups (K = 2 to K = 5) each time with the simulated annealing process repeated 100 times, starting each time with a different partition of the population samples into the K groups. The analysis of molecular variance (AMOVA) was carried out using ARLEQUIN software version 3.11 (Excoffier et al. 2006). Principal Coordinate Analysis (PCoA) in the program GenAIEx 6.5 (Peakall and Smouse 2012) was applied to identify the relationships among populations on the basis of their allele frequencies.

Results

SSR variability

The variability estimates describing the suitability of the ten SSR loci (AfD4, AfD105, AfA7, AfA112, AfA115, AfA120, AfA122, AfA117, AfA10, and AfC103) for detecting the presence of differentiation among the A. fraterculus populations are shown in Table 2. The number of alleles per locus across populations ranges from 7 to 18 with a mean of 12.9, and the mean frequency of null alleles across the loci is generally low (0.04–0.08). Moreover, the Polymorphic Information Content (PIC) estimate for each locus ranges from 0.44 (A10a) to 0.88 (A120a), and the across loci average is 0.74, suggesting that this set of loci is informative for population analyses.

Microsatellite variability detected across the six Brazilian Anastrepha fraterculus populations.

Locus na Min–Max PIC
D4a 7 2–5 0.56
D105a 15 5–11 0.72
A7a 13 7–11 0.83
A112a 18 8–12 0.82
A115a 14 7–12 0.80
A120a 15 7–13 0.88
A122a 12 5–9 0.74
A117a 11 6–9 0.77
A10a 13 2–11 0.44
C103a 11 7–9 0.80
Mean 12.9 5.6–10.2 0.74

Tests for Hardy-Weinberg equilibrium (HWE) using Fisher’s exact test with the sequential Bonferroni correction (Rice 1989) revealed that the populations conformed to Hardy–Weinberg equilibrium (HWE) at most loci. The very few observed locus/populations combinations that were not in HWE were not concentrated at any locus or in any population. Significant linkage disequilibrium was not detected between genotypes at the ten loci. As no evidence of linkage disequilibrium between loci was assessed, these 10 loci can be considered to be independent.

Population variability and differentiation

An estimate of variability distribution in and among the six tested populations (AMOVA) indicates that 90% of the variation occurs within populations while only about 10% of total variation is detected among populations. Indeed as shown in Table 3, the intrapopulation genetic variability is similar across the six samples. As a second step, the simulated annealing approach based on the SAMOVA algorithm was applied to identify the presence of genetically homogeneous groups across the considered Brazilian populations. For this, the spatial analyses of molecular variance was performed without constraints for geographic composition of the groups. As observed in Table 4 no great differences in the FCT values were observed when we increased the group number (K): four of the five simulated groupings (2, 4, and 5) produced non-significant FCT values. This implies that the molecular variance due to differences between populations within each group is weak, confirming the AMOVA data. Only in one case did we observe an FCT estimate which statistically maximized the differences between groups. This is the grouping configuration which splits the six populations into three groups: i) Una, Porto Seguro and São Mateus, ii) Monte Alegre, iii) Campos do Jordão and Vacaria, (FCT = 0.182, P = 0.015). These data indicate that the mountain populations of Campos do Jordão (state of São Paulo) and Vacaria (state of Rio Grande do Sul) are genetically homogeneous and differentiated from the group of the three coastal populations (Una, Porto Seguro and São Mateus), and from the other more distant coastal population of Monte Alegre. The pairwise FST estimates (Table 5) confirm the presence of differentiation between the group of the two mountain populations and the coastal populations which in turn share a certain degree of genetic relatedness. Principal Coordinate Analysis (PCoA) was performed to better clarify the genetic relations among the six populations. The first two axes explain a relatively high amount of the genetic variation (88%). The first axis (77.91%) separates Campos do Jordão and Vacaria from the other populations. The second axis (10.39%) mainly differentiates São Mateus from the group of Una, Porto Seguro, and Monte Alegre, but also Campos do Jordão from Vacaria. It is interesting that there is a certain correspondence of the genetic grouping with the known morphotype classification (Hernández-Ortiz et al. 2012). The molecular variance represented by the first axis separates the populations on the basis of both geographical distance and altitude. On this basis, we attempted to disentangle the effect of geographical distance and altitude on the genetic differentiation of Campos de Jordão and Vacaria. The plots in Figure 3 show the correlations between pairwise FST values with the differences in altitude (Figure 3A) and the geographic distance (Figure 3B), respectively. The correlation results clearly indicate that the difference in altitude has a greater impact (r2 = 0.76 vs r2 = 0.28) on genetic differentiation of mountain Vacaria and Campos do Jordão, than on the remaining samples belonging to the coastal plain populations.

Figure 3.

Correlation of FST values with the geographic distances (upper plot) and altitude differences (bottom plot) among the 6 Brazilian samples of Anastrepha fraterculus.

Genetic variability of wild populations of Anastrepha fraterculus from different geographical regions in Brazil estimated using 10 SSRs.

na He Ho FIS
Una-BA 7,7 0,63 0,54 0,14
Porto Seguro-BA 8,2 0,70 0,67 0,02
Monte Alegre-RN 7,5 0,68 0,57 0,20
São Mateus-ES 6,3 0,66 0,66 0,08
Campos do Jordão-SP 8,3 0,71 0,61 0,13
Vacaria-RS 8,5 0,72 0,62 0,12

Spatial Analysis of Molecular Variance (SAMOVA) for different population partitions.

Number of groups (K) FCT P Population partition
2 0.195 0.062 (Una, Porto Seguro, Monte Alegre, São Mateus), (Campos do Jordão, Vacaria)
3 0.182 0.015 (Una, Porto Seguro, São Mateus), (Monte Alegre), (Campos do Jordão, Vacaria)
4 0.190 0.050 (Una, Porto Seguro, São Mateus), (Monte Alegre), (Campos do Jordão), (Vacaria)
5 0.196 0.068 (Una, Porto Seguro), (São Mateus), (Monte Alegre), (Campos do Jordão), (Vacaria)

Pairwise-FST values among 6 population samples of Anastrepha fraterculus as derived from Microsatellite Analyser (Dieringer and Schlotterer 2003).

Una-BA Porto Seguro-BA Monte Alegre-RN São Mateus-ES Campos do Jordão-SP Vacaria-RS
Una-BA -
Porto Seguro-BA 0.015 ns -
Monte Alegre-RN 0.020 0.012 ns -
São Mateus-ES 0.031 0.031 0.042 -
Campos do Jordão-SP 0.163 0.127 0.130 0.160 -
Vacaria-RS 0.165 0.126 0.132 0.175 0.038 -

Discussion

With this paper we initiated, in the polymorphic A. fraterculus complex, an analysis of the underpinning genetic architecture and its interaction with correlated ecological, biological and morphological traits. In this context, we have used microsatellite markers to perform a genetic analysis of populations from a complex ecological area such as Brazil.

Although the chromosomal location of the considered SSR loci remains unknown, the assessed linkage equilibrium between them suggests they are statistically independent and that their variability patterns might reflect genome-wide patterns across populations. The six considered ecogeographic populations are here represented by highly polymorphic samples, which reflect a high degree of intrapopulation genetic diversity. Indeed only 10% of the total variability (AMOVA) is represented by the differences between the six geographic populations, while greater variability was found within populations.

The spatial analysis of genetic diversity indicates that the levels of diversity among the six populations vary significantly on an eco-geographical basis, as indicated by SAMOVA and PCoA data. More than by geographical distance, the genetic differentiation is influenced by altitude. The multivariate analysis of ten microsatellites depicts a structural pattern, which clearly separates populations on climatic distribution both on latitudinal and altitudinal basis. Particularly, altitude seems to represent a differentiating adaptation, as the main genetic differentiation is that detected between the populations present at higher altitudes (Campos de Jordão and Vacaria) and those populations from sea level. Genetic divergence between populations from low and high altitude areas has already been observed for populations within the A. fraterculus complex using isozymes (Morgante et al. 1980, Steck 1991, Selivon et al. 2005) and mtDNA (Steck and Sheppard 1993, Santos 1999, McPheron et al. 1999, Smith-Caldas et al. 2001, Barr et al. 2005). Steck (1991) concluded that the strong differences in allelic frequency between lowland and Andean populations of A. fraterculus in Venezuela was due to the fact that they actually represent two genetically distinct species albeit morphologically indistinguishable. These allopatric populations can be subject to divergent selection in response to ecological factors that change over large geographical scales such as altitude. Altitude may act as a barrier to gene flow as levels of life history divergence between high- and low-altitude populations can be correlated with levels of post-zygotic reproductive isolation (Orr and Smith 1998).

One interesting observation, which arises from our data, is that the observed structure of Brazilian populations is entangled with the presence of morphotypes. The actual number of these entities and their respective geographic range are questions that remain to be further elucidated. At the moment three different morphotypes are identified in Brazil (Hernández-Ortíz et al. 2012). As it appears clearly, the PCoA analysis depicts a genetic differentiative pattern that overlaps with the distribution of the known morphotypes. Now the open questions are: 1) is the observed population differentiation contributing to the underpinning genetic architecture of the morphotypes associated to these populations? and 2) do morphotypes track environmental variability? In retrospect, a further aim is to clarify the evolutionary relationships between populations, ecotypes, and morphotypes.

Conclusion

The population genetic approach, in addition to improving our knowledge of the underpinning genetic architecture of the A. fraterculus complex, is also important from an applied perspective. The overall level of genetic variability and the presence of differentiation that we detected among the Brazilian populations of A. fraterculus constitute an important contribution for any potential future application of SIT for the control of populations of this fruit fly pest in Brazil.

Acknowledgements

We are grateful to Iara Joachim-Bravo, Adalécio Kovaleski, Adalton Raga, and Miguel Francisco de Souza Jr for providing some of the specimens used in this study and to Roberto A. Zucchi, Keiko Uramoto, and Elton L. Araujo for specimen identification. We thank Jorge Hendrichs for his support for the development of this study and for his comments on an earlier version of this manuscript. Thanks are also due to Carter Robert Miller for reviewing this manuscript. This study had financial support from a FAO/IAEA Research Contract No. 16059 as part of the Coordinated Research Project on “Resolution of cryptic species complexes of tephritid pests to overcome constraints to SIT and international trade” to JGS in which ARM was also involved as a member. Further financial support came from FAO/IAEA Technical Contract No. 16966 to GG. This work is part of the PhD thesis of KLM, which was performed at the Dept. of Biology & Biotechnology, Univ. of Pavia. KLM was a recipient of a doctoral fellowship from Fundação de Amparo à Pesquisa do Estado da Bahia (FAPESB) and a “sandwich fellowship” from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) during her stay at the University of Pavia. JGS is a fellow researcher (grant 305452/2012-6) of the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).

References

  • Aketarawong N, Bonizzoni M, Malacrida AR, Gasperi G, Thanaphum S (2006) Seventeen novel microsatellite markers from an enriched library of the pest species Bactrocera dorsalis sensu stricto. Molecular Ecology Notes 6: 1138–1140. doi: 10.1111/j.1471-8286.2006.01463.x
  • Aketarawong N, Bonizzoni M, Thanaphum S, Gomulski LM, Gasperi G, Malacrida AR, Gugliemino CR (2007) Inferences on the population structure and colonization process of the invasive oriental fruit fly, Bactrocera dorsalis (Hendel). Molecular Ecology 16: 3522–3532. doi: 10.1111/j.1365-294X.2007.03409.x
  • Aketarawong N, Isasawin S, Thanaphum S (2014) Evidence of weak genetic structure and recent gene flow between Bactrocera dorsalis s.s. and B. papayae, across Southern Thailand and West Malaysia, supporting a single target pest for SIT applications. BMC Genetics 15: 70. doi: 10.1186/1471-2156-15-70
  • Araujo EL, Zucchi RA (2006) Medidas do acúleo na caracterização de cinco espécies de Anastrepha do grupo fraterculus (Diptera: Tephritidae). Neotropical Entomology 35: 1–9. doi: 10.1590/S0103-90161998000100017
  • Augustinos AA, Stratikopoulos EE, Drosopoulou E, Kakani EG, Mavragani-Tsipidou P, Zacharopoulou A, Mathiopoulos KD (2008) Isolation and characterization of microsatellite markers from the olive fly, Bactrocera oleae, and their cross-species amplification in the Tephritidae family. BMC Genomics 9: 618. doi: 10.1186/1471-2164-9-618
  • Augustinos AA, Asimakopoulou AK, Papadopoulos NT, Bourtzis K (2011) Cross amplified microsatellites in the European cherry fly, Rhagoletis cerasi: medium polymorphic highly informative markers. Bulletin of Entomological Research 101: 45–52. doi: 10.1017/S0007485310000167
  • Baliraine FN, Bonizzoni M, Osir EO, Lux SA, Mulaa FJ, Zheng L, Gomulski LM, Gasperi G, Malacrida AR (2003) Comparative analysis of microsatellite loci in four fruit fly species of the genus Ceratitis (Diptera: Tephritidae). Bulletin of Entomological Research 93: 1–10. doi: 10.1079/BER2002212
  • Baliraine FN, Bonizzoni M, Guglielmino CR, Osir EO, Lux SA, Mulaa FJ, Gomulski LM, Zheng L, Quilici S, Gasperi G, Malacrida AR (2004) Population genetics of the potentially invasive African fruit fly species, Ceratitis rosa and Ceratitis fasciventris (Diptera: Tephritidae). Molecular Ecology 13: 683–695. doi: 10.1046/j.1365-294X.2004.02105.x
  • Barbará T, Palma-Silva C, Paggi GM, Bered F, Fay MF, Lexer C (2007) Cross-species transfer of nuclear microsatellite markers: potential and limitations. Molecular Ecology 16: 3759–3767. doi: 10.1111/j.1365-294X.2007.03439.x
  • Bonizzoni M, Zheng L, Guglielmino CR, Haymer DS, Gasperi G, Gomulski LM, Malacrida AR (2001) Microsatellite analysis of medfly bioinfestations in California. Molecular Ecology 10: 2515–2524. doi: 10.1046/j.0962-1083.2001.01376.x
  • Bonizzoni M, Guglielmino CR, Smallridge CJ, Gomulski LM, Malacrida AR, Gasperi G (2004) On the origins of medfly invasion and expansion in Australia. Molecular Ecology 13: 3845–3855. doi: 10.1111/j.1365-294X.2004.02371
  • Boykin LM, Shatters RG, Hall DG, Dean D, Beerli P (2010) Genetic Variation of Anastrepha suspensa (Diptera: Tephritidae) in Florida and the Caribbean Using Microsatellite DNA Markers. Molecular Entomology 103: 2214–2222. doi: 10.1603/EC10128
  • Chambers GK, MacAvoy ES (2000) Microsatellites: consensus and controversy. Comparative Biochemistry and Physiology B 126: 455–476. doi: 10.1016/S0305-0491(00)00233-9
  • Cladera JL, Vilardi JC, Juri M, Paulin LE, Giardini MC, Cendra PVG, Segura DF, Lanzavecchia SB (2014) Genetics and biology of Anastrepha fraterculus: research supporting the use of the sterile insect technique (SIT) to control this pest in Argentina. BMC Genetics 15(Suppl 2): S12. doi: 10.1186/1471-2156-15-S2-S12
  • Delatte H, Virgilio M, Simiand C, Quilici S, De Meyer M (2013) Isolation and characterization of sixteen microsatellite markers cross-amplifying in a complex of three African agricultural pests (Ceratitis rosa, C. anonae and C. fasciventris, Diptera: Tephritidae). Conservation Genetics Resources 5: 31–34. doi: 10.1007/s12686-012-9722-6
  • Devescovi F, Abraham S, Roriz AKP, Nolazco N, Castaneda R, Tadeo E, Caceres C, Segura DF, Vera MT, Joachim-Bravo I, Canal N, Rull J (2014) Ongoing speciation within the Anastrepha fraterculus cryptic species complex: the case of the Andean morphotype. Entomologia Experimentalis et Applicata 104: 376–382. doi: 10.1111/eea.12219
  • Dieringer D, Schlötterer C (2003) Microsatelite Analyser (MSA): a platform independent analysis tool for large microsatellite data sets. Molecular Ecology Notes 3: 107–111. doi: 10.1046/j.1471-8286.2003.00351.x
  • Drew RAI, Ma J, Smith S, Hughes JM (2011) The taxonomy and phylogenetic relationships of species in the Bactrocera musae complex of fruit flies (Diptera: Tephritidae: Dacinae) in Papua New Guinea. The Raffles Bulletin of Zoology 59: 145–162. 00b4953cd7f092458a000000-2.pdf
  • Dupanloup I, Schneider S, Excoffier L (2002) A simulated annealing approach to define the genetic structure of populations. Molecular Ecology 11: 2571–2581. doi: 10.1046/j.1365-294x.2002.01650.x
  • Fritz AH, Schable N (2004) Microsatellite loci from the Caribbean Fruit Fly, Anastrepha suspensa (Diptera: Tephritidae). Molecular Ecology Notes 4: 443–445. doi: 10.1111/j.1471-8286.2004.00699.x
  • Hernández-Ortíz V, Gómez-Anaya JA, Sánchez A, McPheron BA, Aluja M (2004) Morphometric analysis of Mexican and South American populations of the Anastrepha fraterculus complex (Diptera: Tephritidae) and recognition of a distinct Mexican morphotype. Bulletin of Entomological Research 94: 487–499. doi: 10.1079/BER2004325
  • Hernández-Ortíz V, Bartolucci AF, Morales-Valles P, Frías D, Selivon D (2012) Cryptic species of the Anastrepha fraterculus complex (Diptera: Tephritidae): A multivariate approach for the recognition of South American morphotypes. Annals of the Entomological Society of America 105: 305–318. doi: 10.1603/AN11123
  • Islam MS, Ruiz-Arce R, McPheron BA (2011) Microsatellite markers for the West Indian fruit fly (Anastrepha obliqua) and cross species amplification in related pest species. Conservation Genetic Resources 3: 549–551. doi: 10.1007/s12686-011-9401-z
  • Kalinowski ST, Taper ML, Marshall TC (2007) Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Molecular Ecology 16: 1099–106. doi: 10.1111/j.1365-294X.2007.03089.x
  • Lanzavecchia SB, Juri M, Bonomi A, Gomulski L, Scannapieco AC, Segura DF, Malacrida A, Cladera JL, Gasperi G (2014) Microsatellite markers from the ‘South American fruit fly’ Anastrepha fraterculus: a valuable tool for population genetic analysis and SIT applications. BMC Genetics 15(Suppl 2): S13. doi: 10.1186/1471-2156-15-S2-S13
  • Malavasi A, Zucchi RA, Sugayama RL (2000) Biogeografia. In: Malavasi A, Zucchi RA (Eds) Moscas-das-frutas de importância econômica no Brasil: conhecimento básico e aplicado. Holos, Ribeirão Preto, 93–98.
  • McPheron BA, Han HY, Silva JG, Norrbom AL (1999) Phylogeny of the genus Anastrepha and Toxotrypana (Trypetinae: Toxotrypanini) based upon 16S rRNA mitochondrial DNA sequences. In: Aluja M, Norrbom AL (Eds) Fruit flies (Tephritidae): Phylogeny and evolution of behavior. CRC Press, Boca Raton, 343–361.
  • Meixner MD, McPheron BA, Silva JG, Gasparich GE, Sheppard WS (2002) The Mediterranean fruit fly in California: evidence for multiple introductions and persistent populations based on microsatellite and mitochondrial DNA variability. Molecular Ecology 11: 891–899. doi: 10.1046/j.1365-294X.2002.01488.x
  • Morgante JS, Malavasi A, Bush GL (1980) Biochemical systematics and evolutionary relationships of Neotropical Anastrepha. Annals of the Entomological Society of America 73: 622–630. doi: 10.1093/aesa/73.6.622
  • Norrbom AL (2004) Host Plant Database for Anastrepha and Toxotrypana (Diptera: Tephritidae: Toxotrypanini), Diptera Data Dissemination Disk. North American Dipterist’s Society, Washington, DC, USA.
  • Norrbom AL, Korytkowski CA, Zucchi RA, Uramoto K, Venable GL, McCormick J, Dallwitz MJ (2012) Anastrepha and Toxotrypana: descriptions, illustrations, and interactive keys. Version: 28th September 2013. http://delta-intkey.com
  • Peakall R, Smouse PE (2012) GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research – an update. Bioinformatics 28: 2537–2539. doi: 10.1093/bioinformatics/bts460
  • Raymond M, Rousset F (1995) GENEPOP v.1.2: population genetics software for exact tests and ecumenicism. Journal of Heredity 86: 248–249.
  • Rice WR (1989) Analyzing tables of statistical tests. Evolution 43: 223–225. doi: 10.2307/2409177
  • Ruiz MF, Milano A, Salvemini M, Eirin-Lopez JM, Perondini ALP, Selivon D, Polito C, Saccone G, Sánchez L (2007a) The gene transformer of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects. PLoS ONE 2: e1239. doi: 10.1371/journal.pone.0001239
  • Ruiz MF, Eirín-López JM, Stefani RN, Perondini ALP, Selivon D, Sánchez L (2007b) The gene doublesex of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects. Development Genes and Evolution 217: 725–731. doi: 10.1007/s00427-007-0178-8
  • Santos P (1999) Variabilidade genética, isolamento reprodutivo e regra de Haldane em espécies de Anastrepha (Diptera: Tephritidae). PhD Dissertation, Universidade de São Paulo, São Paulo, Brazil.
  • Selivon D, Vretos C, Fontes L, Perondini ALP (2006) New variant forms in the Anastrepha fraterculus complex (Diptera: Tephritidae). In: Sugayama RL, Zucchi RA, Ovruski SM, Sivinski J (Eds) Seventh International Symposium on Fruit Flies of Economic Importance, Salvador, Brazil, September 2006. Color Press, Salvador-BA, 253–258.
  • Selivon D, Perondini ALP, Morgante JS (2005) A genetic morphological characterization of two cryptic species of three Anastrepha fraterculus complex (Diptera:Tephritidae). Annals of the Entomological Society of America 98: 367–381. doi: 10.1603/0013-8746(2005)098[0367:AGCOTC]2.0.CO;2
  • Shearman DCA, Gilchrist AS, Crisafulli D, Graham G, Lange C, Frommer MA (2006) Microsatellite markers for the pest fruit fly, Bactrocera papayae (Diptera: Tephritidae) and other Bactrocera species. Molecular Ecology Notes 6: 4–7. doi: 10.1111/j.1471-8286.2006.01024.x
  • Silva JG, Barr NB (2008) Recent advances in molecular systematics of Anastrepha Schiner. In: Sugayama RL, Zucchi RA, Ovruski SM, Sivinski J (Eds) Seventh International Symposium on Fruit Flies of Economic Importance, Salvador, Brazil, September 2006. Color Press, Salvador-BA, 13–18.
  • Silva JG, Meixner MD, McPheron BA, Steck GJ, Sheppard WS (2003) Recent Mediterranean Fruit Fly (Diptera: Tephritidae) Infestations in Florida – A Genetic Perspective. Journal of Economic Entomology 96: 1711–1718. doi: 10.1093/jee/96.6.1711
  • Smith-Caldas MRB, McPheron BA, Silva JG, Zucchi RA (2001) Phylogenetic relationships among species of the fraterculus group (Anastrepha: Diptera: Tephritidae) inferred from DNA sequences of mitochondrial cytochrome oxidase I. Neotropical Entomolology 30: 565–573. doi: 10.1590/S1519-566X2001000400009
  • Steck GJ (1991) Biochemical systematics and population genetic structure of Anastrepha fraterculus and related species (Diptera: Tephritidae). Annals of the Entomological Society of America 84: 10–28. doi: 10.1093/aesa/84.1.10
  • Steck GJ (1999) Taxonomic status of Anastrepha fraterculus. In: Proceedings workshop, The South American fruit fly Anastrepha fraterculus (Wied.): advances in artificial rearing, taxonomic status and biological studies Vienna, Austria; 1999, 13–20, IAEA-TEC. DOC-1064.
  • Steck GJ, Sheppard WS (1993) Mitochondrial DNA variation in Anastrepha fraterculus. In: Aluja M, Liedo P (Eds) Fruit fly biology and management. Springer, New York, 9–14. doi: 10.1007/978-1-4757-2278-9_2
  • Stone A (1942) The fruit flies of the genus Anastrepha. United States Department of Agriculture, Animal and Plant Miscellaneous Publication 439: 1–112.
  • Vaníčková L, Hernández-Ortiz V, Bravo ISJ, Dias V, Roriz AKP, Laumann RA, Mendonça AL, Paranhos BAJ, do Nascimento RR (2015) Current knowledge of the species complex Anastrepha fraterculus (Diptera, Tephritidae) in Brazil. In: De Meyer M, Clarke AR, Vera MT, Hendrichs J (Eds) Resolution of Cryptic Species Complexes of Tephritid Pests to Enhance SIT Application and Facilitate International Trade. ZooKeys 540: 211–237. doi: 10.3897/zookeys.540.9791
  • Vera MT, Cáceres C, Wornoayporn V, Islam A, Robinson AS, De La Vega MH, Hendrichs J, Cayol JP (2006) Mating incompatibility among populations of the South American fruit fly Anastrepha fraterculus (Diptera: Tephritidae). Annals of the Entomological Society of America 99: 387–397. doi: 10.1603/0013-8746(2006)099[0387:MIAPOT]2.0.CO;2
  • Virgilio M, Delatte H, Backeljau T, De Mayer M (2010) Macrogeographic population structuring in the cosmopolitan agricultural pest Bactrocera cucurbitae (Diptera: Tephritidae). Molecular Ecology 19: 2713–2714. doi: 10.1111/j.1365-294X.2010.04662.x
  • Virgilio M, Delatte H, Quilici S, Backeljau T, De Meyer M (2013) Cryptic diversity and gene flow among three African agricultural pests: Ceratitis rosa, Ceratitis fasciventris and Ceratitis anonae (Diptera, Tephritidae). Molecular Ecology 22: 2526–2539. doi: 10.1111/mec.12278
  • Weir BS, Cockerham CC (1984) Estimating F-Statistics for the Analysis of Population Structure. Evolution 38: 1358–1370. doi: 10.2307/2408641
  • Yamada S, Selivon D (2001) Rose, an eye color mutation in a species of the Anastrepha fraterculus complex (Diptera, Tephritidae). Annals of the Entomological Society of America 94: 592–595. doi: 10.1603/0013-8746(2001)094[0592:RAECMI]2.0.CO;2