- Research article
- Open Access
Genetic variability in Italian populations of Drosophila suzukii
© The Author(s). 2017
- Received: 12 April 2017
- Accepted: 4 October 2017
- Published: 3 November 2017
Drosophila suzukii is a highly destructive pest species, causing substantial economic losses in soft fruit production. To better understand migration patterns, gene flow and adaptation in invaded regions, we studied the genetic structure of D. suzukii collected across Italy, where it was first observed in 2008. In particular, we analysed 15 previously characterised Simple Sequence Repeat (SSR) markers to estimate genetic differentiation across the genome of 278 flies collected from nine populations.
The nine populations showed high allelic diversity, mainly due to very high heterozygosity. The high Polymorphism Information Content (PIC) index values (ranging from 0.68 to 0.84) indicated good discrimination power for the markers. Negative fixation index (F IS) values in seven of the populations indicated a low level of inbreeding, as suggested by the high number of alleles. STRUCTURE, Principal Coordinate and Neighbour Joining analysis also revealed that the Sicilian population was fairly divergent compared to other Italian populations. Moreover, migration was present across all populations, with the exception of the Sicilian one, confirming its isolation relative to the mainland.
This is the first study characterising the genetic structure of the invasive species D. suzukii in Italy. Our analysis showed extensive genetic homogeneity among D. suzukii collected in Italy. The relatively isolated Sicilian population suggests a largely human-mediated migration pattern, while the warm climate in this region allows the production of soft fruit, and the associated D. suzukii reproductive season occurring much earlier than on the rest of the peninsula.
- Spotted wing drosophila
- SSR markers
- Population structure
- Gene flow
- Human trade
The spotted wing drosophila (SWD), Drosophila suzukii Matsumura (Diptera: Drosophilidae), is a pest species which has spread from its original range in Asia to a number of western countries in the past decade, including the Mediterranean basin , Europe, and USA [1–3]. The history of the geographical spread and infestation of D. suzukii is still under investigation: it is known that in 1939 this species was first recorded in Japan (Kanzawa 1939), while in the 1980s it was collected on the island of Hawaii . Europe and the Americas were colonised much later, possibly during the last 9 years [2, 3, 5, 6]. First adults of D. suzukii were caught contemporaneously in the region of Catalonia, Spain  and in Tuscany, Central Italy, in 2008 . In 2009 D. suzukii individuals were found on both wild hosts (Vaccinium, Fragaria and Rubus spp.) and several species of cultivated berries in Trento Province, North Eastern Italy, where also the first economically important damage by this species in Europe was reported . During the following years, D. suzukii has been spreading rapidly across Europe, with documented infestations ranging from Mediterranean regions (i.e. Greece, Turkey) to northern latitudes (i.e. Sweden, Poland, UK) (EPPO Global Database, Drosophila suzukii – DROSSU, 2017). In Italy, after the first detection, infestations were reported from the regions of Bolzano, Piedmont, Liguria, Campania and Veneto in 2010, from Lombardy, Emilia Romagna, Marche, Aosta Valley, Marche, Calabria and Sicily in 2011 , Sardinia in 2012 , Apulia in 2013 , Umbria in 2014 , and Latium in 2015 (Antonini G., present paper). Invasion dynamics can be studied using molecular markers that can discriminate and characterise the genetic relationships between source and derived populations, migration flows and population expansion patterns [12–14]. In particular, Single Nucleotide Polymorphism (SNP) and Simple Sequence Repeat (SSR) markers have played an increasingly significant role in the study of genetic differentiation across species populations . Thanks to their great discrimination power and high reproducibility and variability, SSRs represent one of the most robust and informative molecular markers available for genotyping individuals . For instance, their use in Drosophila species was pivotal in studying intra-population genetic variation and evolution [14, 17–20].
In relation to D. suzukii, SSRs have been exploited to study genetic aspects of the colonisation process in the USA and Europe. Jeffrey and colleagues based their research on the use of six X-linked genes and suggested that the invasions of the USA and Europe are two independent events . Bahder et al. in particular analysed samples of D. suzukii populations collected in California and Washington and determined that while D. suzukii in the former region had high levels of genetic variation, the latter was highly monomorphic . Furthermore, Fraimout’s group investigated Hawaiian and Spanish populations by exploiting microsatellite markers, finding a significant level of genetic differentiation . Although both studies exploited two different sets of microsatellites and tested different populations, the authors were led to similar conclusions: they demonstrated the presence of a specific differentiation process among ancestral and derived populations and suggested that for D. suzukii a genetic analysis approach is valuable not only to better understanding of the evolutionary history of the species, but also to manage its great potential for invasiveness. Different studies on the invasiveness of species, including D. suzukii, have demonstrated the relationship between their spread and human trade [24–26]. For this reason, it is very important to consider the correlation between gene flow analysis and the sale of soft fruit all over the country. Taking into account this aspect, we used a population genetic approach to characterise genetic diversity among D. suzukii individuals collected in different regions of Italy. In order to perform this work a set of 15 microsatellites validated by Fraimout and colleagues  were employed. The current research is the first study that provides new insights on the trend of genetic diversity in Italian populations of D. suzukii.
D. Suzukii collection, identification and DNA extraction
Each pair of primers was used for PCR amplification in 25 μL final volume, containing 1X GoTaq G2 Master Mix, 0.5 μL of each primer, 10.5 μL of distilled deionized water and 1 μL of genomic DNA. The PCR program was set with an initial period of denaturation at 94 °C (30 s) followed by 32 cycles of additional denaturation at 94 °C (30 s), an annealing phase at 57 °C (1 min 30 s), an elongation phase at 72 °C (1 min), and ending with another extension phase at 72 °C (30 min). PCR products were checked using electrophoresis on 1.5% agarose gel, stained with ethidium bromide and visualised under UV light. Each amplicon was then diluted 1:10 in distilled water and 1 μL of this dilution was added to 12.5 μL of a mixture of deionised Formamide (Sigma-Aldrich) and GeneScan-500 ROX size standard (Life Tech, Waltham, MA USA). Prior to denaturation for 4 min at 94 °C, capillary electrophoresis was carried out in an ABI PRISM 310 Genetic Analyzer (Life Tech) and the fragments were sized with GeneMapper v.4.0 software in binning mode. If no sample amplification was obtained after two PCR attempts, the locus was classified as missing data.
Microsatellite allele data were processed with Tandem program v.1.08 . GenAIEx software v.6.41  was run to study the genetic variability between populations using the following statistics: mean number of alleles (N a), effective number of alleles (N e), expected heterozygosity (H E), observed heterozygosity (H O), number of private alleles (N p), frequency of private alleles (A p) and inbreeding coefficient (F IS). Allelic richness was calculated using FSTAT v.2.9.3 software . Deviation from Hardy-Weinberg equilibrium after the Bonferroni multiple correction test and allelic Polymorphic Information Content (PIC) were tested using CERVUS software v.3.0 . N e and H E were chosen as the basic genetic variability and estimated for each population. N e was analysed with ANOVA using origin as a factor. N e was taken from the formula 1/(1-H E) and then tested with the non-parametric Tukey test . N e was used instead of (N a), considering that it is less sensitive to rare alleles and sample size. H E was taken from the formula H E = 1- (Σq i 2), where q i represents the frequency of the i th allele in the population. H E was converted into 1/H E and then tested with the non-parametric Kruskal Wallis test. All statistical analyses were performed using R software v.3.3.2. The significance level was set below 0.001 (P < 0.001) to minimise sources of uncertainty.
To evaluate the genetic structure of populations, we relied on multiple approaches: Principal Coordinate Analysis (PCoA), Neighbour Joining Tree, AMOVA analysis, measurement of the index of differentiation (F ST) and use of a non-spatial Bayesian algorithm. These approaches were chosen in order to obtain a broad view of the genetic structure of this invasive species in Italy. PCoA, obtained with GenAIEx software, was used to display genetic divergence across D. suzukii in a multidimensional space, considering frequency data. Unrooted Neighbour Joining Tree based on Nei’s genetic distance constructed using DARwin software was complementary to PCoA analysis . AMOVA analysis obtained using the Arlequin v.3.5  program was performed to estimate variability distribution within and between the tested groups. The level of genetic differentiation in populations was detected using the F ST values obtained with Microsatellite Analyzer (MSA) v.4.05 software . The program allows comparison of each observed F ST value with that obtained in 10,000 matrix permutations in order to define the statistical significance of each F ST.
The Bayesian method was implemented with STRUCTURE software v.2.3.3 [36, 37]. This program was employed in order to obtain clusters of individual genotypes. The analysis was run using the admixture hypothesis, which is based on correlated allele frequencies, in which each sample contains a portion of the genome of each ancestral population. This, correlated to the allele frequency model, allows calculation of the log likelihood for the data, L(K). Not knowing the origin and the degree of isolation of the studied populations a priori, this model is considered to be the most appropriate in these situations . Prior probability, i.e. the probability that an individual belongs to any K reference populations, is defined as l/K. The K value was fixed from 1 to 10 with 20 replicates of each K to test the convergence of the Markov chain. A total of 1000,000 simulations per run and 500,000 Markov Chain Monte Carlo MCMC repetitions were fixed. Once the results were obtained, they were scored with STRUCTURE HARVESTER software to detect the number of K groups that best fit the dataset according to the Evanno test [38, 39]. GENECLASS v.2.0  was run to estimate the probability of each individual in a population belonging only to that population, the probability of it being an immigrant from each of the other populations, and the probability of it being a migrant to the other populations. BOTTLENECK v.1.2.02  was run in order to evaluate whether demographic events such as population contraction or expansion took place in each population.
Heterozygosity excess, which is associated with a population expansion, was tested with the two-phase mutation model (TPM) using Wilcoxon signed-rank test, which according to Piry et al. is the most appropriate and powerful test when dealing with less than twenty loci . Parameters were set as 20% multiple-step mutations and 80% single-step mutations with 1000 iterations. In order to verify the effect of isolation by distance, and therefore to find possible correlation between genetic and geographical distances, the ISOLDE option in GENEPOP software was run.
Summary of the genetic variability at 15 microsatellite loci
Level of genetic diversity across 9 populations of D. suzukii
6.40 ± 1.29
3.65 ± 0.93
0.68 ± 0.18
0.70 ± 0.07
7.60 ± 1.50
3.83 ± 1.01
0.66 ± 0.16
0.72 ± 0.08
7.87 ± 1.40
4.37 ± 0.80
0.82 ± 0.10
0.76 ± 0.05
7.80 ± 1.37
4.18 ± 0.87
0.83 ± 0.13
0.74 ± 0.06
5.60 ± 0.98
3.25 ± 0.58
0.87 ± 0.11
0.68 ± 0.05
8.87 ± 3.44
5.22 ± 2.53
0.81 ± 0.13
0.77 ± 0.07
7.47 ± 1.12
4.68 ± 0.98
0.89 ± 0.09
0.77 ± 0.04
7.40 ± 1.50
4.41 ± 0.77
0.87 ± 0.10
0.76 ± 0.03
7.53 ± 1.30
4.14 ± 0.50
0.83 ± 0.16
0.75 ± 0.02
Genetic population structure and gene flow
Analysis of molecular variance test (AMOVA)
Source of variation
Pairwise Fst among the nine populations of D. suzukii
Individual assignment analysis obtained using GENECLASS
Results from the individual population analysis using the TPM model did not support an expansion scenario. On the contrary, population bottleneck, defined by significant heterozygote deficiency was present in the Trentino2 population (P = 0.004). The nine populations showed no significant correlation when comparing genetic and geographical distances [R 2 = 0.014, P = 0.568, F ST/(1-F ST) = 0.049 + Ln (geographical distance) = 0.001].
The introduction of invasive species to new environments poses threats to biodiversity, agriculture, public health and ecosystem integrity [44–47]. For this reason, considerable attention is paid to the rapid spread of alien species [46, 48]. Genetic characteristics deeply affect the capacity for expansion . Therefore, in order to mitigate their impact and define management strategies it is imperative to study these fundamental characteristics. Currently techniques such as genomics [50–52], transcriptomics [53, 54], and metagenomics [55, 56] allow us to investigate these basic traits.
This research investigated the genetic structure of D. suzukii collected in different areas of Italy. In particular, the aim of the analyses was to understand the gene flow of this species in a newly colonised environment. Our findings help to better understand the dynamics and complexity of this invasive species in Italy. The nine populations studied show a high level of genetic variation. The high number of alleles per locus detected clearly demonstrated the discriminatory power of these markers. Taking into consideration N e, H E and H O, it is evident that the level of genetic differentiation is similar in D. suzukii collected across Italy, even in the locations at the greatest distance from the likely spreading centre of the species in Italy . The high level of heterozygosity could be explained by good adaptation to new ranges due to a favourable environment, their reproductive power, and the absence or limited presence of natural competitors and predators [57, 58]. Bahder et al. found that populations from Washington were much less polymorphic than those in California, suggesting a recent strong population bottleneck associated to the recent invasion of the former . Washington has a much cooler climate than California, similar to the contrast between Trentino and the rest of Italy. However, we did not observe such a contrast in heterozygosity, probably due to the highly favourable habitat found in Trentino coupled with a high migration rate with the rest of the Italian populations. Heterozygosity deficiency was detected for one out of the nine investigated populations. In contrast, the remaining eight groups showed heterozygosity excess. Sicily in particular had the lowest heterozygosity value (−0.28). Negative results indicate random mating, therefore a lack of inbreeding among the collected individuals. In contrast, Trentino1 and Trentino2 had a positive F IS value, indicating inbreeding. The Apulian population showed the greatest number of private alleles (25). This could be the consequence of a steady introduction of new alleles due to migration, possibly associated with human-mediated transport [25, 59].
Genetic structure analysis
Moderate genetic differentiation between most of the groups was in evidence for the nine populations, while the Sicilian population was the most differentiated from the others. This is supported by the NJ tree, PCoA data and structural analysis. At the same time, low differentiation between the other populations may be due to gene flow, which can homogenize gene frequency across populations. Data concerning the reduction or expansion of the studied populations indicated that Trentino2 was the only group having indices of genetic bottleneck.
Human transportation is the most probable explanation for the extensive spread of D. suzukii. [25, 59]. When an alien species is introduced into an environment outside its native range, expansion can be identified not only by analysing genetic diversity indices, but also by analysing the genetic flow between populations, which is a direct proof of rapid distribution. [60, 61]. In particular, in the last 40 years, the risk of biotic invaders has increased significantly because of levels of international trade not seen before . This situation facilitates genetic flow between groups located in different locations, and may well apply to the results of our study. For instance, the observation that the level of heterozygosity (F st) does not clearly decline (increase) from the hypothetical source population (Livorno, Tuscany) and that there is a high migration rate among localities, suggests that D. suzukii moves extensively across most of the Italian peninsula. Most of the Sicilian production of vegetables and fruit, including high D. suzukii susceptible hosts, is frequently exported to central and northern Italy. While this could suggest a high probability of flies being transported between Sicily and the rest of the peninsula, our results indicate that there was no gene flow from Sicily to other regions. This is probably due to the fact that ripe fruits are exported from Sicily mostly during the cold season, when moderate temperatures allow the production of berry fruit in Sicily, but not in the rest of Italy. Therefore, any D. suzukii accidentally moving from Sicily to the rest of Italy would arrive at a time when the local population is made up of a few individuals in winter diapause .
A second interesting piece of information revealed by our results is related to the scenario in Sardinia. To satisfy local demand for berry fruit, this region imports fruit from Italy and northern Europe, Spain, the USA and South America. The flies used in this study were collected in Arborea, a town 13 km away from the port of Oristano, one of the most important commercial ports in Italy. Thus, it is likely that the Sardinian population is made up of immigrants from other regions, as suggested by the low differentiation between this population and those on the mainland.
This research represents the first study investigating the pattern of genetic variability for D. suzukii following its introduction to Italy. Defining the population structure of a species, in particular of an invasive species, it is necessary not only to improve our knowledge of the genetic architecture, but also to apply knowledge. Indeed, understanding the current genetic structure of D. suzukii has significant implications in relation to geographical and economic impact. The evaluation of the genetic status of the D. suzukii populations in newly invaded areas and their expansion or reduction phases during defined periods of the year, may thus provide valuable information for predicting population spread, outbreaks, and improve integrated pest management programmes. Proper genetic management practices for D. suzukii and constant monitoring are therefore critical for maintaining populations under control.
The information obtained can be applied in particular to the management of coastal areas; one important action could be to increase monitoring control with the use of traps and other early warning tools in order to limit either multiple reintroductions of the same species or new introductions of exotic organisms.
The authors are grateful to Prof. Giuseppe Firrao (University of Udine) for his availability and support, to Dr. Umberto Salvagnin for his assistance in figure preparation, and to Drs. Lara Giongo, Alberto Grassi and Elisabetta Leonardelli for their helpful comments. We also thank Vivienne Frankell and Linda Brewer who helped with the English language revision.
Funding for this research was provided by Fondazione Edmund Mach.
Availability of data and materials
The data supporting the results of this article are included in the article and the supplementary information.
GT contributed to all the steps, planning the experimental design, sampling populations, conducting laboratory tests and data analysis, as well as writing the manuscript. SV participated in genetic data analysis and drafting the manuscript. SV and LO contributed to data interpretation. FS contributed to laboratory experiments. GfA and GlA contributed to sampling populations, suggesting professional contacts and conceiving the main idea. AB, NB, GS, AC and LT contributed to sampling populations. All the authors have read and approved the final manuscript.
Ethics approval and consent to participate
No specific permits were required for this project. Drosophila suzukii is an agricultural pest, not a protected species. All the insects analysed were collected in the open field and not from national parks or protected areas.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Asplen MK, Anfora G, Biondi A, Choi DS, Chu D, Daane KM, Gibert P, Gutierrez AP, Hoelmer KA, Hutcison WD, Isaacs R, Jiang ZL, Kàrpàti Z, Kinura MT, Pascual M, Philips CR, Plantamp C, Ponti L, Vétek G, Vogt H, Walton VM, Yu Y, Zappalà L, Desneux N. Invasion biology of spotted wing drosophila (Drosophila suzukii): a global perspective and future priorities. J Pest Sci. 2015;88(3):469–94.View ArticleGoogle Scholar
- Hauser M. A historic account of the invasion of Drosophila suzukii (Matsumura) (Diptera: Drosophilidae) in the continental United States, with remarks on their identification. Pest Manag Sci. 2011;67(11):1352–7.View ArticlePubMedGoogle Scholar
- Cini A, Ioriatti C, Anfora G. A review of the invasion of Drosophila suzukii in Europe and a draft research agenda for integrated pest management. Bull Insectol. 2012;65(1):149–60.Google Scholar
- Orhan A, Aslantaş R, Önder BŞ, Tozlu G. First record of the invasive vinegar fly Drosophila suzukii (Matsumura) (Diptera: Drosophilidae) from eastern Turkey. Turkish J Zool. 2016;40(2):290–3.View ArticleGoogle Scholar
- Lee JC, Bruck DJ, Dreves AJ, Ioriatti C, Vogt H, Baufeld P. In focus: spotted wing drosophila, Drosophila suzukii, across perspectives. Pest Manag Sci. 2011;67(11):1349–51.View ArticlePubMedGoogle Scholar
- Deprá M, Poppe JL, Schmitz HJ, De Toni DC, Valente VLS. The first records of the invasive pest Drosophila suzukii in the south American continent. J Pest Sci. 2014;87(3):379–83.View ArticleGoogle Scholar
- Calabria G, Máca J, Bächli G, Serra L, Pascual M. First records of the potential pest species Drosophila suzukii (Diptera: Drosophilidae) in Europe. J Appl Entomol. 2012;136(1–2):139–47.View ArticleGoogle Scholar
- Grassi A, Palmieri L, Giongo L. Nuovo fitofago per i piccoli frutti in Trentino. Terra Trentina. 2009;55(10):19–23. http://hdl.handle.net/10449/19083Google Scholar
- Marongiu G, Marras PM, Ramassini W. First record of Drosophila suzukii in Sardinia (Italy). Protezione delle Colture. 2013;3:34.Google Scholar
- Baser N, Ouantar M, Broutou O, Lamaj F, Verrastro V, Porcelli F. First finding of Drosophila suzukii (Matsumura) (Diptera : Drosophilidae) in Apulia, Italy, and its population dynamics throughout the year. Fruits. 2015;70(4):225–30.View ArticleGoogle Scholar
- Romani R, Natalini G. Prima segnalazione di Drosophila suzukii in Umbria. L'Informatore Agrario. 2015;71(7):64–5.Google Scholar
- Sunnucks P, De Barro PJ, Lushai G, Maclean N, Hales D. Genetic structure of an aphid studied using microsatellites: cyclic parthenogenesis, differentiated lineages and host specialization. Mol Ecol. 1997;6:1059–73.View ArticlePubMedGoogle Scholar
- Jarne P, Lagoda PJL. Microsatellites, from molecules to populations and back. Trends Ecol Evol. 1996;11(10):424–9.View ArticlePubMedGoogle Scholar
- Schlötterer C. The evolution of molecular markers -just a matter of fashion? Nat Rev Genet. 2004;5(1):63–9.View ArticlePubMedGoogle Scholar
- Morin PA, Mccarthy M. Highly accurate SNP genotyping from historical and low-quality samples. Mol Ecol Notes. 2007;7(6):937–46.View ArticleGoogle Scholar
- Murphy KA, West JD, Kwok RS, Chiu JC. Accelerating research on spotted wing drosophila management using genomic technologies. J Pest Sci. 2016;89(3):631–41.View ArticleGoogle Scholar
- Hutter CM, Schug MD, Aquadro CF. Microsatellite variation in Drosophila melanogaster and Drosophila simulans: a reciprocal test of the ascertainment bias hypothesis. Mol Biol Evol. 1998;15(12):1620–36.View ArticlePubMedGoogle Scholar
- Amos W, Hutter CM, Schug MD, Aquadro CF. Directional evolution of size coupled with ascertainment bias for variation in drosophila microsatellites. Mol Biol Evol. 2003;20(4):660–2.View ArticlePubMedGoogle Scholar
- Kauer MO, Dieringer D, Schlötterer C. A microsatellite variability screen for positive selection associated with the “out of Africa” habitat expansion of Drosophila melanogaster. Genetics. 2003;165(3):1137–48.PubMedPubMed CentralGoogle Scholar
- Abdulazeez R, Auta J, Shehu DM, IS N, Umar AU. Determination of genetic similarities in natural populations of Drosophila Melanogaster in savanna zone of Nigeria using microsatellite markers. Mol Biol. 2016;5:177.Google Scholar
- Adrion JR, Kousathanas A, Pascual M, Burrack HJ, Haddad NM, Bergland AO, Machado H, Sackton TB, Schlenke TA, Wegmann MWD, Singh ND. Drosophila suzukii: the genetic footprint of a recent, worldwide invasion. Mol Biol Evol. 2014;31(12):3148–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Bahder BW, Bahder LD, Hamby KA, Walsh DB, Zalom FG. Microsatellite variation of two Pacific coast Drosophila suzukii (Diptera: Drosophilidae) populations. Environ Entomol. 2015;44(5):1449–53.View ArticlePubMedGoogle Scholar
- Fraimout A, Loiseau A, Price DK, Xuereb A, Martin J-F, Vitalis R, Fellous S, Debat V, Estoup A. New set of microsatellite markers for the spotted-wing Drosophila suzukii (Diptera: Drosophilidae): a promising molecular tool for inferring the invasion history of this major insect pest. Eur J Entomol. 2015;112(4):855–9.Google Scholar
- Hulme PE. Trade, transport and trouble: managing invasive species pathways in an era of globalization. J Appl Ecol. 2009;46(1):10–8.View ArticleGoogle Scholar
- Cini A, Anfora G, Escudero-Colomar LA, Grassi A, Santosuosso U, Seljak G, Papini A. Tracking the invasion of the alien fruit pest Drosophila suzukii in Europe. J Pest Sci. 2014;87(4):559–66.View ArticleGoogle Scholar
- Manni M, Guglielmino CR, Scolari F, Vega-Rua A, Failloux A-B, Somboon P, Lisa A, Savini G, Bonizzoni M, Gomulski LM, Malacrida AR, Gasperi G. Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector, Aedes albopictus. PLoS Negl Trop Dis. 2017;11(1):e0005332.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossi-Stacconi MV, Kaur R, Mazzoni V, Ometto L, Grassi A, Gottardello A, Rota Stabelli O, Anfora G. Multiple lines of evidence for reproductive winter diapause in the invasive pest Drosophila suzukii: useful clues for control strategies. J Pest Sci. 2016;89(3):689–700.View ArticleGoogle Scholar
- Voorrips RE. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93(1):77–8.View ArticlePubMedGoogle Scholar
- Matschiner M, Salzburger W. TANDEM: integrating automated allele binning into genetics and genomics workflows. Bioinformatics. 2009;25(15):1982–3.View ArticlePubMedGoogle Scholar
- Peakall R, Smouse PE. GENALEX 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.View ArticleGoogle Scholar
- Goudet J. FSTAT (version 1.2): a computer program to calculate F-statistics. J Hered. 1995;86(6):485–6.View ArticleGoogle Scholar
- Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16(5):1099–106.View ArticlePubMedGoogle Scholar
- Aketarawong N, Isasawin S, Sojikul P, Thanaphum S. Gene flow and genetic structure of Bactrocera carambolae (Diptera, Tephritidae) among geographical differences and sister species, B. dorsalis, inferred from microsatellite DNA data. Zookeys. 2015;2015(540):239–72.View ArticleGoogle Scholar
- Pickersgill B, Botany SE. Genetic diversity of cultivated tropical plants book review editor. Daniel F Austin. 2006;60(1):91–2.Google Scholar
- Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005;1:47–50.Google Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.PubMedPubMed CentralGoogle Scholar
- Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164(4):1567–87.PubMedPubMed CentralGoogle 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.View ArticlePubMedGoogle Scholar
- Earl DA, Von Holdt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.View ArticleGoogle Scholar
- Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A. GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered. 2004;95(6):536–9.View ArticlePubMedGoogle Scholar
- Piry S, Luikart G, Cornuet JM. BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999;90(4):502–3.View ArticleGoogle Scholar
- Chapuis MP, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24(3):621–31.View ArticlePubMedGoogle Scholar
- Morin PA, Leduc RG, Archer FI, Martien KK, Huebinger R, Bickham JW, Taylor BL. Significant deviations from hardy-Weinberg equilibrium caused by low levels of microsatellite genotyping errors. Mol Ecol Resour. 2009;9(2):498–504.View ArticlePubMedGoogle Scholar
- Mooney HA, Cleland EE. The evolutionary impact of invasive species. Proc Natl Acad Sci U S A. 2001;98(10):5446–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Harris RJ, Barker G. Relative risk of invasive ants (hymenoptera: Formicidae) establishing in New Zealand. New Zeal J Zool. 2007;34(3):161–78.View ArticleGoogle Scholar
- Pyšek P, Richardson DM. Invasive species, environmental change and management, and health. Annu Rev Environ Resour. 2010;35(1):25–55.View ArticleGoogle Scholar
- Gasperi G, Bellini R, Malacrida AR, Crisanti A, Dottori M, Aksoy S. A new threat looming over the Mediterranean Basin: emergence of viral diseases transmitted by Aedes albopictus mosquitoes. PLoS Negl Trop Dis. 2012;6(9):e1836.View ArticlePubMedPubMed CentralGoogle Scholar
- Allendorf FW, Lundquist LL. Population biology, evolution, and control of invasive species. Conserv Biol. 2003;17(1):24–30.View ArticleGoogle Scholar
- Lee CEE. Evolutionary genetics of invasive species. Trends Ecol Evol. 2002;17(8):386–91.View ArticleGoogle Scholar
- Ometto L, Cestaro A, Ramasamy S, Grassi A, Revadi S, Siozios S, Moretto M, Fontana P, Varotto C, Pisani D, Dekker T, Wrobel N, Viola R, Pertot I, Cavalieri D, Blazter M, Anfora G, Rota SO. Linking genomics and ecology to investigate the complex evolution of an invasive drosophila pest. Genome Biol Evol. 2013;5(4):745–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Crava CM, Ramasamy S, Ometto L, Anfora G, Rota-Stabelli O. Evolutionary insights into taste perception of the invasive Pest Drosophila suzukii. G3 (Bethesda). 2016;6:4185–96.Google Scholar
- Ramasamy S, Ometto L, Crava MC, Revadi S, Kaur R, Horner DS, Pisani D, Dekker T, Anfora G, Rota SO. The evolution of olfactory gene families in drosophila and the genomic basis of chemical-ecological adaptation in Drosophila suzukii. Genome Biol Evol. 2016;8(8):2297–311.View ArticlePubMedPubMed CentralGoogle Scholar
- Kučerová L, Kubrak OI, Bengtsson JM, Strnad H, Nylin S, Theopold U, Nassel DR. Slowed aging during reproductive dormancy is reflected in genome-wide transcriptome changes in Drosophila melanogaster. BMC Genomics. 2016;17(1):50.View ArticlePubMedPubMed CentralGoogle Scholar
- Giampetruzzi A, Morelli M, Saponari M, Loconsole G, Chiumenti M, Boscia D, Savino VN, Martelli GP, Saldarelli P. Transcriptome profiling of two olive cultivars in response to infection by the CoDiRO strain of Xylella fastidiosa subsp. pauca. BMC Genomics. 2016;17(1):475.Google Scholar
- Champion CJ, Xu J. The impact of metagenomic interplay on the mosquito redox homeostasis. Free Radic Biol Med. 2017;105:79–85.View ArticlePubMedGoogle Scholar
- Sardi S, Somasekar S, Naccache SN, Bandeira AC, Tauro LB, Campos GS, Chiu CY. Co-infections from Zika and Chikungunya virus in Bahia, Brazil identified by metagenomic next-generation sequencing. J Clin Microbiol. 2016;54(9):2348–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller B, Anfora G, Buffington M, Daane KM, Dalton DT, Hoelmer KM, Rossi Stacconi MV, Grassi A, Ioriatti C, Loni A, Miller JC, Quantar M, Wang W, Wiman NG, Walton VM. Seasonal occurrence of resident parasitoids associated with Drosophila suzukii in two small fruit production regions of Italy and the USA. Biol Con. 2015;68(2):255–63.Google Scholar
- Mazzetto F, Marchetti E, Amiresmaeili N, Sacco D, Francati S, Jucker C, Dindo L, Lupi D, Tavella L. Drosophila parasitoids in northern Italy and their potential to attack the exotic pest Drosophila suzukii. J Pest Sci. 2016;89(3):837–50.View ArticleGoogle Scholar
- Lengyel GD, Orosz S, Kiss B, Lupták R, Kárpáti Z. New records and present status of the invasive spotted wing drosophila, Drosophila suzukii (Matsumura, 1931) (Diptera) in Hungary. Acta Zool Acad Sci Hungaricae. 2015;61(1):73–80.View ArticleGoogle Scholar
- Benedict MQ, Levine RS, Hawley WA, Lounibos LP. Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus. Vector Borne Zoonotic Dis. 2007;7(1):76–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Goergen G, Vayssières J-F, Gnanvossou D, Tindo M. Bactrocera invadens (Diptera: Tephritidae), a new invasive fruit fly Pest for the Afrotropical region: host plant range and distribution in west and Central Africa. Environ Entomol. 2011;40(4):844–54.View ArticlePubMedGoogle Scholar
- Pimentel D, Zuniga R, Morrison D. Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecol Econ. 2005;52(3):273–88.View ArticleGoogle Scholar