Genetic Variation May Have Promoted the Successful Colonization of the Invasive Gall Midge, Obolodiplosis robiniae, in China

Invasive species often cause serious economic and ecological damage. Despite decades of extensive impacts of invasives on bio-diversity and agroforestry, the mechanisms underlying the genetic adaptation and rapid evolution of invading populations remain poorly understood. The black locust gall midge, Obolodiplosis robiniae, a highly invasive species that originated in North America, spread widely throughout Asia and Europe in the past decade. Here, we used 11 microsatellite DNA markers to analyze the genetic variation of 22 O. robiniae populations in China (the introduced region) and two additional US populations (the native region). A relatively high level of genetic diversity was detected among the introduced populations, even though they exhibited lower diversity than the native US populations. Evidence for genetic differentiation among the introduced Chinese populations was also found based on the high Fst value compared to the relatively low among the native US populations. Phylogenetic trees, structure graphical output, and principal coordinate analysis plots suggested that the Chinese O. robiniae populations (separated by up to 2,540 km) cluster into two main groups independent of geographical distance. Genetic variation has been observed to increase rapidly during adaptation to a new environment, possibly contributing to population establishment and spread. Our results provide insights into the genetic mechanisms underlying successful invasion, and identify factors that have contributed to colonization by an economically important pest species in China. In addition, the findings improve our understanding of the role that genetic structure plays during invasion by O. robiniae.

Introduction

Obolodiplosis robiniae (Haldeman, 1847) (Diptera: Cecidomyiidae) is a North American species of gall midge that has recently been extensively introduced throughout Asia and Europe (CABI/EPPO, 2011) and is continuously expanding its range (Cierjacks et al., 2013; Stalazs, 2014; Badmin, 2016; Kostro-Ambroziak and Mieczkowska, 2017). It is specifically associated with host plants from the genus Robinia (Fabaceae) (Stalazs, 2014). Its main host is Robinia pseudoacacia although it is occasionally found on R. pseudoacacia cv. ‘Frisia’ (Badmin, 2016). The gall midge causes leaf rolling and premature leaf shedding, resulting in the deterioration of the host and increased susceptibility to other pests, including wood borers such as longhorn beetles (Yang et al., 2006).

Obolodiplosis robiniae was first recorded in China (Qinhuangdao City, Hebei Province) in 2004 (Yang et al., 2006) and has since spread extensively. Its primary host (R. pseudoacacia) has been planted extensively across China, and O. robiniae is now found in most of these areas (Shang et al., 2015a). Chinese O. robiniae populations may produce between four and six generations per year (Wang, 2009; Mu et al., 2010; Shao et al., 2010; Liu, 2014), which is significantly higher than the rate in regions beyond China. For example, O. robiniae produces three to four generations per year in Italy (Duso et al., 2005) and Serbia (Mihajlovic et al., 2008), and a maximum of three generations per year in Korea (Lee et al., 2009).

Genetic diversity and population structure are important factors affecting the colonization of invasive species (Amouroux et al., 2013; Horst and Lau, 2015; Zhao et al., 2015). Invasive species often exhibit low genetic diversity during founding events, as new habitats are typically colonized by only a few individuals, representing a small proportion of the allelic diversity present in the source population (Nei et al., 1975; Tsutsui et al., 2003). However, when the founding individuals originate from multiple source populations, the genetic diversity of the founder population can be relatively high (Davis, 2009). This can contribute to invasion success by facilitating local adaptation to new environments and increasing new trait diversity (Facon et al., 2006). Besides, the invaders can rapidly evolve in isolation from other individuals of the same species when they were introduced into the new environments (Lee, 2002; Launey et al., 2010).

DNA-based molecular markers have been extensively used to examine the genetic diversity and population structure of a wide range of species. Microsatellite DNA markers (simple sequence repeats, SSRs) are suitable for routine genetic diversity analyses (Varshney et al., 2007; Kong et al., 2014; Kim et al., 2017), as they are ubiquitous among eukaryotes (Sharma et al., 2007), co-dominantly inherited, and highly polymorphic (Zong et al., 2015). Moreover, microsatellite analysis can yield valid results and improved phylogenetic trees compared to analyses involving other molecular markers (Schemerhorn et al., 2015). Due to their feasibility and practicality, microsatellite markers have been widely used in population genetics and ecological studies of various insects (Bonizzoni et al., 2000; Mezghani-Khemakhem et al., 2012; Anjos et al., 2016; Retamal et al., 2016; Duan et al., 2017; Kim et al., 2017; Simonato et al., 2019), including several invasive gall midge species (Bentur et al., 2011; Amouroux et al., 2013).

Previously, Shang et al. (2015b) investigated the genetic variation among Chinese O. robiniae populations using a partial mitochondrial DNA cytochrome c oxidase subunit I (COI) sequence marker. However, only 10 individuals exhibiting haplotypic variation and a mere four haplotypes were detected in 560 O. robiniae samples. Thus, the genetic mechanisms behind successful invasion, the genetic structure in the process of colonization, and the phylogenetic relationships among the Chinese O. robiniae populations remain poorly understood.

Accordingly, to gain further insight into the genetic structure of the Chinese O. robiniae populations and ascertain how the species has spread widely in new regions, we used 11 microsatellite markers to analyze the genetic structure of 22 Chinese O. robiniae populations. Two native populations from the United States (US) were also assessed based on the same loci for comparison with the Chinese populations, in order to explain how genetic diversity is altered during the invasion process.

Materials and Methods

Sample Collection

We collected the gall midge larvae and pupae contained within rolled leaves of host trees growing in 22 cities across China (Figure 1). Generally, the rolled leaves were randomly picked from different trees; however, when infestation was low, individual trees were singled out for sample collection. Following collection, the rolled leaves were immediately transported to the laboratory in 60 cm × 40 cm plastic bags, in which they were maintained until adult emergence. From the samples collected at each location, 20 larger adults were selected, placed into a 1.5-mL centrifuge tube, and stored at −20°C for subsequent DNA extraction. Additionally, eight O. robiniae adults from two regions of the United States were obtained from the Quarantine Lab at the Institute of Forest Ecology, Environment, and Protection in the Chinese Academy of Forestry. Details of the sample collection and population codes are listed in Table 1.

FIGURE 1

Figure 1. Geographic locations of the 22 Chinese Obolodiplosis robiniae populations sampled. Population codes are listed in Table 1. (The source map was downloaded from the website http://www.webmap.cn/commres.do?method=dataDownload).

TABLE 1

Table 1. Location of Obolodiplosis robiniae populations and the sample size used in this study.

DNA Extraction and Microsatellite Analyses

Genomic DNA was extracted from the entire O. robiniae body following the instructions described by Zhou et al. (2007) and stored at −20°C until needed. The 14 microsatellite loci (W3, W5, W8, W29, W31, W33, W35, W41, W46, W82, W83, W116, W126, and W132) developed by Yao et al. (2015) were initially selected to analyze the genotypes of 20 individuals per collection site (the exception being Zhengzhou, for which 16 individuals were analyzed) (Table 1). For each sample, we attempted to amplify all 14 loci; however, after two attempts, we were unable to amplify five loci (W29, W35, W41, W46, and W116) for many individuals; thus, these loci were not used in subsequent analyses. However, we assessed the applicability of two additional loci (W6 and W107; GenBank numbers: KP260520 and KP260530) that were not characterized by Yao et al. (2015), and we detected sufficient polymorphism among the analyzed samples. Hence, a total of 11 loci were used to genotype 444 O. robiniae individuals.

Microsatellite amplifications were performed in a 15 μL reaction volume containing 1 μL genomic DNA (10 ng), 1 μL of each primer (5 μmol/L), 7.5 μL 2X Taq PCR Master Mix (TIANGEN, Beijing, China), and 4.5 μL ddH₂O. The forward primer of each primer pair was labeled with a fluorescent dye (HEX, ROX, FAM, or TMARA; Sangon Biotech, Shanghai, China). The microsatellite cycling protocol was: 5 min at 95°C (initial denaturation step); followed by 30 cycles of 94°C for 30 s, 53°C (W3, W5, W6, and W8) or 56°C (the remaining loci) for 45 s, 72°C for 45 s, extension at 72°C for 10 min, and finally maintained at 16°C. PCR products were examined using a DNA analyzer (Applied Biosystems, Waltham, CA, United States) and the results were analyzed using Genotyping was carried out using a 3730xl automated DNA sequencer (Applied Biosystems, Waltham, CA, United States). Alleles were scored using GeneMarker software version 2.2. (Softgenetics LLC, State College, PA, United States).

Data Analyses

Genetic diversity was estimated by basic statistical analyses including number of alleles (Na), effective number of alleles (Ne), Shannon’s information index (I), observed heterozygosity (Ho), expected heterozygosity (He), and Nei’s (1973) expected heterozygosity (Nei), which were calculated using GenePop software version 4.3 (Rousset, 2008); genotype number (GN), gene diversity (GD), and polymorphism information content (PIC) were calculated using PowerMarker software version V3.25 (Liu and Muse, 2005). F-statistics and gene flow for each locus across populations were performed using PopGene software version 1.32 (Yeh et al., 2018). Deviations from the Hardy–Weinberg equilibrium (HWE) based on the Markov chain algorithm (10,000 steps) and linkage disequilibrium (LD) (10,000 permutations) were also examined by GenePop software. The genetic relationships between populations were assessed using a neighbor-joining dendrogram generated by PowerMarker and Molecular Evolutionary Genetics Analysis across Computing Platforms (MEGA X) (Kumar et al., 2018).

Populations differentiation was assessed by pairwise Fst values (based on 999 permutations) and gene flow (Nm) through AMOVA (analysis of molecular variance) approach which were performed using the Arlequin program version 3.5 (Excoffier and Lischer, 2010). The analyses can estimate variance and partitioning of the within- and among-population. Genetic structure analysis was performed with 100,000 Markov Chain Monte Carlo repetitions after a burn-in period of 200,000 interactions for each group number (K) using STRUCTURE software version 2.3.4 (Pritchard et al., 2000). The number of subpopulations (K) was assumed to be from 1 to 22, without admixture and with correlated allele frequencies. To determine the most likely number of subpopulations, the optimum K-value was obtained by calculating the ΔK value (Evanno et al., 2005).

In addition, the Mantel test was conducted using the GenALEx 6.5 program (Peakall and Smouse, 2012) to determine correlations between Nei’s genetic distance [was calculated using GenePop software based on Nei (1978)] and both geographical distance (km) and altitude (m). Significance was assessed by conducting 999 permutations. Moreover, a principal coordinate analysis (PCoA) was conducted using the same software.

Results

Microsatellite Polymorphism and Diversity

In this study, locus polymorphism and diversity were determined based on 22 Chinese populations using 11 microsatellite markers, with each population consisting of 20 individual samples (except Zhengzhou, with 16 samples) (Table 1). Amplifying these microsatellite markers loci led to 436 polymorphic bands (Supplementary Table S1), representing 202 genotypes, ranging from 3 to 54 per primer pair. As shown in Table 2, there were 72 alleles among the 22 populations; the number of alleles (Na) observed per locus varied from 2 (W83 and W107) to 14 (W3 and W5), with a mean of 6.5 per locus. The effective numbers of alleles (Ne) varied from 1.4087 (W83) to 7.9971 (W3), with an average of 3.7255 per locus. The gene diversity index (GD) per locus ranged from 0.2901 (W83) to 0.8755 (W3), with an average of 0.6511, indicating that a high level of information was provided by the 11 microsatellite markers. Shannon’s information index (I) ranged from 0.4767 (W83) to 2.2573 (W3), with a mean of 1.3188. The polymorphism information content (PIC) for the microsatellite loci ranged from 0.2494 (W83) to 0.8627 (W3), with an average of 0.6031. The observed heterozygosity (Ho) ranged from 0.2271 (W107) to 0.7477 (W3) and expected heterozygosity (He) ranged from 0.2905 (W83) to 0.8760 (W3). For each locus, both Na and Ho values markedly changed among populations, whereas, He value minorly altered (Supplementary Figure S1). In addition, when we evaluated the within-sample HWE deviations for each locus across all 22 Chinese populations using the Markov chain algorithm (10,000 steps), we detected significant deviation from the expected value (p < 0.05) in 68 of 242 tests (28.10%) (Supplementary Table S2). Examination of genotypic LD between all pairs of alleles across all 22 populations, based on a permutation procedure (10,000 permutations), revealed a significant LD (p < 0.05) in 251 of 1210 tests (20.74%) from 11 loci in the 22 Chinese populations.

TABLE 2

Table 2. Polymorphism of microsatellite loci across Chinese O. robiniae populations.

Population Genetic Diversity

The genetic diversity of 22 Chinese populations and two US populations was assessed. Six indices of genetic diversity (Na, Ne, I, Ho, He, and Nei) were evaluated. As shown in Table 3, for each Chinese population across all loci, the means of the above indices except Na were moderately or considerably lower than those of the native US populations, although the sample size of the Chinese populations (40) was markedly higher than that of the US populations (8). Regarding the three most important indices, I, He, and Nei, the lowest Chinese values (I = 0.7847, He = 0.4279, Nei = 0.4172) occurred in the DY population, and the highest (I = 1.1103, He = 0.6162, Nei = 0.6008) in the TS population. Increased values of these three indices occurred in the SY (I = 1.058, He = 0.6157, Nei = 0.6003), YT (I = 1.0761, He = 0.5825, Nei = 0.5677), and DD (I = 1.016, He = 0.5895, Nei = 0.5748) populations. The inbreeding coefficient (Fis) ranged from −0.0295 (TY) to 0.2011 (CD), with significant various observed for each locus among populations. Both US populations exhibited high He values of 0.6544 (US_f) and 0.6606 (US_g).

TABLE 3

Table 3. Genetic diversity of the O. robiniae populations across 11 microsatellite loci.

Genetic Differentiation in the Chinese Populations

The Fst per locus ranged from 0.1357 to 0.3770, with an average of 0.1994, and the Fis per locus ranged from −0.0037 (W3) to 0.3364 (W82) with an average of 0.0873 alleles per locus across populations (Table 4). Gene flow (Nm) ranged from 0.3093 at W31 to 1.5921 at W33 and averaged 1.0036. Meanwhile, the pairwise Fst (p < 0.001) values (Supplementary Table S3) between populations ranged from 0.022 (HF and WH) to 0.377 (BJ and GY), with an average value of 0.183. A total of 135 of the 213 Fst values (63.38%) were >0.15, while 48 (22.54%) were >0.25, which suggests that significant genetic differentiation exists among the sampling sites and there is some restriction in gene flow between them (Table 4). The most noticeable genetic differentiation occurred between BJ and GY (Fst = 0.377), followed by GY and QH (Fst = 0.372), then DY and QH, DY and GY (Fst = 0.361 for both). According to the coefficient of genetic differentiation (Fst = 0.1830, p < 0.001), genetic variation within populations (81.66%) was substantially higher than that among populations (18.34%) (p < 0.001) (Table 5). Gene flow (Nm) ranged from 0.414 (BJ and GY) to 11.05 (HF and WH) (Supplementary Table S4), with an average of 1.113 (Table 5). All investigated loci contributed to the population differentiation (p < 0.001 for each individual locus). Regarding the native US populations, they had a relatively low Fst value (0.085, p < 0.03), a high Nm value (2.685), and variation within populations of 91%, while variation among populations was 9% (p < 0.02). This further indicates the existence of extensive genetic differentiation among the introduced populations.

TABLE 4

Table 4. Summary of F statistics and gene flow for each locus.

TABLE 5

Table 5. Population genetic variance revealed by 11 microsatellite loci through AMOVA analysis.

Genetic Relationships and Population Structure Analysis

A dendrogram depicting the genetic relationships among the 22 Chinese populations was constructed based on the microsatellite data (Figure 2). The populations were divided into two main clusters, and each cluster was further separated into several sub-clusters. Group I contained populations from Northeast China (Jilin Province) to Southwest China (Guizhou Province), including Liaoning (YK, DL), Jilin (CC), Shanxi (TY), Shaanxi (YA), Shandong (DY, QD), Henan (ZZ), Sichuan (CD), Hubei (WH), Anhui (HF), and Guizhou (GY) provinces. Group II contained populations from North China (with the exception of the NJ population), including Beijing (BJ), Liaoning (SY, DD), Hebei (QH), Shandong (YT, TA), Shaanxi (XA), Gansu (TS), Ningxia (YC), and Jiangsu (NJ) provinces. Besides, some subdivided populations were clustered according to their spatial distribution, such as GY, HF, and WH located in the south of southern China, which clustered together in group I, whereas TA, DD, BJ, and QH located around Bohai Bay clustered together in group II.

FIGURE 2

Figure 2. Unrooted neighbor-joining dendrogram of the 24 O. robiniae populations based on Nei’s distance using the allele frequencies of 11 microsatellite loci. Green represents the subpopulations of group 1, red represents the subpopulations of group 2, and blue represents the US populations.

The 436 Chinese O. robiniae samples were further assessed for population stratification using STRUCTURE software. Microsatellite data were analyzed with possible cluster numbers (K-values) ranging from 1 to 22. ΔK was clearly maximized when K = 2 (ΔK = 4298.3743), indicating the occurrence of two distinct groups among the 22 populations (Figure 3), which validates the dendrogram-based grouping, and the second clade was grouped according to approximate geographical area. These results suggested different degrees of introgression in the populations, detected as differences in allelic frequencies among the populations. In addition, greater structuring (K = 3) revealed that QD, GY, and TY in group I had certain structural similarities to YT, SY, NJ, YC, and XA in group II.

FIGURE 3

Figure 3. Graphical output of the STRUCTURE analysis representing hierarchical data analyses to determine the number of genetic subpopulations (K) of O. robiniae. Each individual is represented by a single vertical bar.

In addition, PCoA based on the marker genotypes also revealed two distinct clusters of the Chinese populations (Figure 4), which were partly related to their geographical regions (group II contained populations from North China). The Mantel test revealed non-significant negative correlations between Nei’s genetic distance (Supplementary Table S5) and geographical distance (km) (r = −0.02, P = 0.456; Supplementary Figure S2), and between Nei’s genetic distance and altitude (r = −0.026, P = 0.419; Supplementary Figure S3), which indicates that genetic differentiation in the 22 Chinese populations may not be caused by geographical isolation.

FIGURE 4

Figure 4. Principal coordinate analysis (PCoA) of 436 Chinese O. robiniae individuals showing two distinct clusters of populations. Each population is represented by a diamond. Coord.1 (39.51%) and Coord.2 (17.63%) refer to the first and second principal component, respectively.

Discussion

In the present study, we detected a high degree of polymorphism among the assessed microsatellite loci. We also identified a relatively high level of genetic diversity among Chinese O. robiniae populations across all loci, with the average expected heterozygosity (He) and Nei’s (1973) expected heterozygosity (Nei) being 0.5346 and 0.5192, respectively. The highest He was 0.6606, which occurred in the US_g population despite the fact that it consisted of only three O. robiniae individuals. Nonetheless, He, Nei (Shang et al., 2016), gene diversity index (GD), and polymorphism information content (PIC) (Ren et al., 2014) are minimally influenced by sample size. Our results indicate significant differences in genetic diversity within the Chinese populations, as well as between the native and invasive populations. Shang et al. (2015b) detected a relatively low level of genetic diversity among Chinese O. robiniae populations using a COI marker. This discrepancy suggests that COI markers may be less suitable than microsatellite DNA markers for population analyses of a new invasive species, such as Chinese O. robiniae, as is the case with another invasive cecidomyiid, Procontarinia mangiferae (Amouroux et al., 2013).

The relatively high level of genetic diversity has likely contributed to the spread of O. robiniae across China in the past decade to the extent that this species is now established in most regions where its host exists (Shang et al., 2015a). O. robiniae has a short history in China, with initial detection occurring in 2004 (Yang et al., 2006), whereas its host was introduced over a century ago (Xu and Yang, 2006). Moreover, for the 22 Chinese populations analyzed, our results show a relatively high level of genetic differentiation (Fst = 0.1830) among populations from sites separated by distances of up to 2,540 km. Hence, the relatively high genetic diversity likely caused by significant differentiation among populations has been conducive to the rapid colonization and establishment of O. robiniae in China, and it is an important factor contributing to the successful invasion of O. robiniae.

The observed population differentiation likely resulted from rapid evolution during adaptation to the new environment. Invasive species may evolve rapidly in response to selection pressures driven by novel habitats (Sakai et al., 2001; Lee, 2002; Ochocki and Miller, 2017), and such rapid genetic adaptation might be important for invasive species (Shine et al., 2011) in order to increase fitness and invasion success (Suarez and Tsutsui, 2008; Roux and Wieczorek, 2009). Increasing the success of both their initial establishment and subsequent range expansion is a particularly effective strategy for introduced populations, as was shown for several invasive species during their colonization processes (Simberloff et al., 2013; Ochocki and Miller, 2017; Wu et al., 2019). Hence, differentiation among Chinese populations was likely accelerated by the rapid evolution of adaptations to the new environments, promoting successful invasion by O. robiniae.

Furthermore, high levels of genetic diversity in an invasive species might be caused by multiple introductions or large founding populations. It is thought that multiple introductions are associated with increased diversity because they supply increased variation and new genetic communities (Dlugosch and Parker, 2008). Multiple introductions are considered to produce invasive populations that are much more genetically diverse than a single source population (Sakai et al., 2001). As such, the successful establishment and invasion of many invasive species have been attributed to multiple introductions (Meixner et al., 2002; Kolbe et al., 2004; Cheng et al., 2008; Zalewski et al., 2010; Michaelides et al., 2018). For O. robiniae, in light of the short history in China, its high genetic diversity and colonization success might be also related to multiple independent invasive events.

Many studies have shown that the genetic structures of invasive species are well developed in their new ranges (Zeisset and Beebee, 2003; Herborg et al., 2007; Rollins et al., 2009; Zalewski et al., 2010). Indeed, the genetic diversity of invasive species is often higher than that of native populations (Marrs et al., 2008). In light of this, our neighbor-joining dendrogram and Bayesian STRUCTURE (K = 2) analyses indicated that the Chinese O. robiniae populations are divided into two independent clusters, although this division appears to be unrelated to geographical distribution. Meanwhile, gene flow was found among some Chinese populations, with QD, GY, and TY in group I having highly similar structures to YT, SY, NJ, YC, and XA in group II. In addition, some subgroups (GY, HF, and WH; TA, DD, BJ, and QH) were clustered according to their geographical distribution, which likely represents different routes of spread in the new environment. Furthermore, our results revealed that each group of the Chinese O. robiniae populations exhibited differences in geographical distribution and genetic distance, suggesting that the two groups do in fact represent two different sources. This implies the introduced populations likely experienced two independent invasive events, which initially shaped the genetic structure of the Chinese O. robiniae populations.

On the other hand, the fact that the division of the two groups was unrelated to geographical distribution (particularly for group I, which contained numerous genetically similar populations located in different geographical regions) suggests that human activity is likely another important contributing factor. Transport, business trips, and long-distance vacations have recently increased not only in frequency but also in distance. High levels of human-mediated dispersion can increase the genetic diversity of an invasive population, thereby substantially modifying the genetic structure and potential management units (Perkins et al., 2013); these factors decrease the success of control measures for O. robiniae populations.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

Y-XY and W-XZ: conceptualization, funding acquisition, project administration, and supervision. R-ZL, W-XH, and JY: data curation. Y-XY, X-PS, and JY: formal analysis, investigation, and resources. Y-XY and X-PS: methodology. Y-XY and R-ZL: software. Y-XY, R-ZL, and W-XH: validation. Y-XY, X-PS, R-ZL, and W-XH: visualization. Y-XY, X-PS, R-ZL, W-XH, JY, and W-XZ: writing – original draft and review and editing.

Funding

This work was funded by the Fundamental Research Funds of the Chinese Academy of Forestry (CAFYBB2017SZ003) and the National Key R&D Program of China (2016YFC1201200).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Dr. Yan-Fei Zeng (Research Institute, Chinese Academy of Forestry) for discussions regarding the data. We also sincerely thank the other members of the Plant Quarantine Lab and Plant Pathology Lab of the Chinese Academy of Forestry for their helpful comments.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00387/full#supplementary-material

References

Amouroux, P., Normand, F., Nibouche, S., and Delatte, H. (2013). Invasive mango blossom gall midge, Procontarinia mangiferae (Felt) (Diptera: Cecidomyiidae) in Reunion Island: ecological plasticity, permanent and structured populations. Biol. Invasions 15, 1677–1693. doi: 10.1017/S0007485314000480

PubMed Abstract | CrossRef Full Text | Google Scholar

Anjos, L. O., Peixoto, R. F. Jr., Chanquinie, D. M., Pinto, L. R., Creste, S. A., Dinardo-Miranda, L. L., et al. (2016). Microsatellite loci and genetic structure of artificial populations of Cotesia flavipes (Hymenoptera, Braconidae). Genet. Mol. Res. 15:gmr.15048851. doi: 10.4238/gmr15048851

PubMed Abstract | CrossRef Full Text | Google Scholar

Badmin, J. (2016). Obolodiplosis robiniae (Haldeman) (Diptera: Cecidomyiidae) on False–acacia cv Frisia. Br. J. Entomol. Nat. Hist. 29:245.

Google Scholar

Bentur, J. S., Sinha, D. K., Padmavathy, R. C., Muthulakshmi, M., and Nagaraju, J. (2011). Isolation and characterization of microsatellite loci in the Asian rice gall midge (Orseolia oryzae) (Diptera: Cecidomyiidae). Int. J. Mol. Sci. 12, 755–772. doi: 10.3390/ijms12010755

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonizzoni, M., Malacrida, A. R., Guglielmino, C. R., Gomulski, L. M., Gasperi, G., and Zheng, L. (2000). Microsatellite polymorphism in the Mediterranean fruit fly, Ceratitis capitata. Insect Mol. Biol. 9, 251–261. doi: 10.1046/j.1365-2583.2000.00184.x

PubMed Abstract | CrossRef Full Text | Google Scholar

CABI/EPPO (2011). Obolodiplosis robiniae. Distribution Maps of Plant Pests. Wallingford: CABI.

Google Scholar

Cheng, X. Y., Cheng, F. X., Xu, R. M., and Xie, B. Y. (2008). Genetic variation in the invasive process of Bursaphelenchus xylophilus (Aphelenchida: Aphelenchoididae) and its possible spread routes in China. Heredity 100, 356–365. doi: 10.1038/sj.hdy.6801082

PubMed Abstract | CrossRef Full Text | Google Scholar

Cierjacks, A., Kowarik, I., Joshi, J., Hempel, S., Ristow, M., Lippe, M., et al. (2013). Biological flora of the British Isles: Robinia pseudoacacia. J. Ecol. 101, 1623–1640.

Google Scholar

Davis, M. A. (2009). Invasion Biology. New York, NY: Oxford University Press.

Google Scholar

Dlugosch, K. M., and Parker, I. M. (2008). Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol. Ecol. 17, 431–449. doi: 10.1111/j.1365-294x.2007.03538.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, X., Wang, K., Su, S., Tian, R., Li, Y., and Chen, M. (2017). De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae). PLoS One 12:e0172513. doi: 10.1371/journal.pone.0172513

PubMed Abstract | CrossRef Full Text | Google Scholar

Duso, C., Fontana, P., and Tirello, P. (2005). Spread of the gall midge Obolodiplosis robiniae (Haldeman) injurious to black locust in Italy and Europe. Inform. Fitopatol. 55, 30–33.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Excoffier, L., and Lischer, H. E. L. (2010). Arlequin suite version 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Res. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Facon, B., Genton, B. J., Shykoff, J., Jarne, P., Estoup, A., and David, P. A. (2006). General eco-evolutionary framework for understanding bioinvasions. Trends Ecol. Evol. 21, 130–135. doi: 10.1016/j.tree.2005.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Haldeman, S. S. (1847). Description of several new and interesting animals. Am. J. Agr. Sci. 6, 191–194.

Google Scholar

Herborg, L. M., Weetman, D., Van Oosterhout, C., and Hanfling, B. (2007). Genetic population structure and contemporary dispersal patterns of a recent European invader, the Chinese mitten crab, Eriocheir sinensis. Mol. Ecol. 16, 231–242. doi: 10.1111/j.1365-294x.2006.03133.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Horst, C. P., and Lau, J. A. (2015). Genetic variation in invasive species response to direct and indirect species interactions. Biol. Invasions 17, 651–659. doi: 10.1111/j.1469-185X.2010.00123.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. R., Kim, K. Y., Jeong, J. S., Kim, M. J., Kim, K. H., Choi, K. H., et al. (2017). Population genetic characterization of the Japanese oak silkmoth, Antheraea yamamai (Lepidoptera: Saturniidae), using novel microsatellite markers and mitochondrial DNA gene sequences. Genet. Mol. Res. 16:gmr16029608. doi: 10.4238/gmr16029608

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolbe, J. J., Glor, R. E., Schettino, L. R. G., Lara, A. C., Larson, A., and Losos, J. B. (2004). Genetic variation increases during biological invasion by a Cuban lizard. Nature 431, 177–181. doi: 10.1038/nature02807

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, L., Bai, J., and Li, Q. (2014). Comparative assessment of genomic SSR, EST–SSR and EST–SNP markers for evaluation of the genetic diversity of wild and cultured Pacific oyster, Crassostrea gigas Thunberg. Aquaculture 420-42, S85–S91.

Google Scholar

Kostro-Ambroziak, A., and Mieczkowska, A. (2017). The first record of the black locust gall midge Obolodiplosis robiniae (Haldeman, 1847) (Diptera: Cecidomyiidae) from northeastern Poland. Wiadomosci Entomol. 36:755.

Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096

PubMed Abstract | CrossRef Full Text | Google Scholar

Launey, S., Brunet, G., Guyomard, R., and Davaine, P. (2010). Role of introduction history and landscape in the range expansion of brown trout (Salmo trutta L.) in the Kerguelen Islands. J. Hered. 10, 270–283. doi: 10.1093/jhered/esp130

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C. E. (2002). Evolutionary genetics of invasive species. Trends Ecol. Evol. 17, 386–391. doi: 10.1016/s0169-5347(02)02554-5

CrossRef Full Text | Google Scholar

Lee, J. S., Jung, Y. M., Choi, K. S., Kim, I. K., Kwon, Y. D., Jeon, M. J., et al. (2009). Seasonal fluctuation and distribution of Obolodiplosis robiniae (Diptera: Cecidomyiidae) within crown of Robinia pseudoacacia (Fabaceae). Korean J. App. Entomol. 48, 447–451. doi: 10.5656/ksae.2009.48.4.447

CrossRef Full Text | Google Scholar

Liu, K., and Muse, S. V. (2005). Powermarker: integrated analysis environment for genetic marker data. BMC Bioinformatics 21:2128–2129. doi: 10.1093/bioinformatics/bti282

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. B. (2014). Management methods to control Obolodiplosis robiniae. Gansu Agric. 2, 87–88.

Google Scholar

Marrs, R. A., Sforza, R., and Hufbauer, R. A. (2008). When invasion increases population genetic structure: a study with Centaurea diffusa. Biol. Invasions 10, 561–572. doi: 10.1007/s10530-007-9153-6

CrossRef Full Text | Google Scholar

Meixner, M. D., Mcpheron, B. A., Silva, J. G., Gasparich, G. E., and Sheppard, W. S. (2002). The Mediterranean fruit fly in California: evidence for multiple introductions and persistent populations based on microsatellite and mitochondrial DNA variability. Mol. Ecol. Notes 11, 891–899. doi: 10.1046/j.1365-294x.2002.01488.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mezghani-Khemakhem, M., Bouktila, D., Casse, N., Maaroufi, H., Makni, M., and Makni, H. (2012). Development of new polymorphic microsatellite loci for the barley stem gall midge, Mayetiola hordei (Diptera: Cecidomyiidae) from an enriched library. Int. J. Mol. Sci. 13, 14446–14450. doi: 10.3390/ijms131114446

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaelides, S. N., Goodman, R. M., Crombie, R. I., Kolbe, J. J., and Cowie, R. (2018). Independent introductions and sequential founder events shape genetic differentiation and diversity of the invasive green anole (Aanolis carolinensis) on pacific islands. Divers. Distrib. 24, 666–679. doi: 10.1111/ddi.12704

CrossRef Full Text | Google Scholar

Mihajlovic, L., Glavendekic, M. M., Jakovljevic, I., and Marjanovic, S. (2008). Obolodiplosis robiniae (Haldeman) (Diptera: Cecidomyiidae) – a new invasive insect pest on black locust in Serbia. Glas. Sumarskog Fakult. Univ. Beogr. 97, 197–207. doi: 10.2298/gsf0897197m

CrossRef Full Text | Google Scholar

Mu, X. F., Sun, J. S., Lu, W. F., Li, M., Qu, H. X., and Gao, Z. Y. (2010). Bionomics and control of Obolodiplosis robiniae in Beijing. For. Pest Dis. 29, 15–18.

Google Scholar

Nei, M. (1973). Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. 70, 3321–3323. doi: 10.1073/pnas.70.12.3321

PubMed Abstract | CrossRef Full Text | Google Scholar

Nei, M. (1978). Estimation of average heterozygosity and genetic distance from a Small Number of Individuals. Genetics 89, 583–590.

Google Scholar

Nei, M., Maruyama, T., and Chakraborty, R. (1975). The bottleneck effect and genetic variability in populations. Evolution 29, 1–10. doi: 10.1111/j.1558-5646.1975.tb00807.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ochocki, B. M., and Miller, T. E. X. (2017). Rapid evolution of dispersal ability makes biological invasions faster and more variable. Nat. Commun. 8:14315. doi: 10.1038/ncomms14315

PubMed Abstract | CrossRef Full Text | Google Scholar

Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analyses in Excel. Population genetic software for teaching and research—an update. BMC Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460

PubMed Abstract | CrossRef Full Text | Google Scholar

Perkins, T. A., Phillips, B. L., Baskett, M. L., and Hastings, A. (2013). Evolution of dispersal and life history interact to drive accelerating spread of an invasive species. Ecol. Lett. 16, 1079–1087. doi: 10.1111/ele.12136

PubMed Abstract | CrossRef Full Text | Google Scholar

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

Google Scholar

Ren, X. P., Jiang, H. F., Yan, Z. Y., Chen, Y., Zhou, X. J., Huang, L., et al. (2014). Genetic diversity and population structure of the major peanut (Arachis hypogaea L.) cultivars grown in China by SSR markers. PLoS One 9:e88091. doi: 10.1371/journal.pone.0088091

PubMed Abstract | CrossRef Full Text | Google Scholar

Retamal, R., Zaviezo, T., Malausa, T., Fauvergue, X., Le, G. I., and Toleubayev, K. (2016). Genetic analyses and occurrence of diploid males in field and laboratory populations of Mastrus ridens (Hymenoptera: Ichneumonidae), a parasitoid of the codling moth. Biol. Control 101, 69–77. doi: 10.1016/j.biocontrol.2016.06.009

CrossRef Full Text | Google Scholar

Rollins, L. A., Woolnough, A. P., Wilton, A. N., Sinclair, R., and Sherwin, W. B. (2009). Invasive species can’t cover their tracks: using microsatellites to assist management of starling (Sturnus vulgaris) populations in Western Australia. Mol. Ecol. 18, 1560–1573. doi: 10.1111/j.1365-294x.2009.04132.x

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Roux, L. J., and Wieczorek, A. M. (2009). Molecular systematics and population genetics of biological invasions: towards a better understanding of invasive species management. Ann. Appl. Biol. 154, 1–17. doi: 10.1111/j.1744-7348.2008.00280.x

CrossRef Full Text | Google Scholar

Sakai, A. K., Allendorf, F. W., Holt, J. S., Lodge, D. M., Molofsky, J., With, K. A., et al. (2001). The population biology of invasive species. Ann. Rev. Ecol. Syst. 32, 305–332.

Google Scholar

Schemerhorn, B. J., Crane, Y. M., Cambron, S. E., Crane, C. F., and Shukle, R. H. (2015). Use of microsatellite and SNP markers for biotype characterization in Hessian fly. J. Insect Sci. 15:158. doi: 10.1093/jisesa/iev138

PubMed Abstract | CrossRef Full Text | Google Scholar

Shang, X. P., Yao, Y. X., Huai, W. X., and Zhao, W. X. (2016). Effects of sample size on genetic diversity index for population of Obolodiplosis robiniae in use of microsatellite DNA marker. Plant Quar. 1, 32–35.

Google Scholar

Shang, X. P., Yao, Y. X., and Zhao, W. X. (2015a). Geographic distribution of an invasive insect pest,Obolodiplosis robiniae (Diptera: Cecidomyiidae) in China. For. Pest Dis. 34, 33–36.

Google Scholar

Shang, X. P., Yao, Y. X., Huai, W. X., and Zhao, W. X. (2015b). Population genetic differentiation of the black locust gall midge Obolodiplosis robiniae (Haldeman) (Diptera: Cecidomyiidae): a North American pest invading Asia. Bull. Entomol. Res. 105, 736–742. doi: 10.1017/s000748531500070x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao, X. K., Ma, X. G., Shao, K. F., Lv, J., and Han, G. S. (2010). Occurrence, damage and control of Obolodiplosis robiniae. J. Liaon. For. Sci. Technol. 4, 31–32.

Google Scholar

Sharma, P. C., Grover, A., and Kahl, G. (2007). Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 25, 490–498. doi: 10.1016/j.tibtech.2007.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Shine, R., Brown, G. P., and Phillips, B. L. (2011). An evolutionary process that assembles phenotypes through space rather than through time. Proc. Natl. Acad. Sci. U.S.A. 108, 5708–5711. doi: 10.1073/pnas.1018989108

PubMed Abstract | CrossRef Full Text | Google Scholar

Simberloff, D., Martin, J. L., Genovesi, P., Maris, V., Wardle, D. A., Aronson, J., et al. (2013). Impacts of biological invasions: what’s what and the way forward. Trends Ecol. Evol. 28, 58–66. doi: 10.1016/j.tree.2012.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Simonato, M., Pilati, M., Magnoux, E., Courtin, C., Sauné, L., Rousselet, J., et al. (2019). A population genetic study of the egg parasitoid Baryscapus servadeii reveals large scale automictic parthenogenesis and almost fixed homozygosity. Biol. Control 139:104097. doi: 10.1016/j.biocontrol.2019.104097

CrossRef Full Text | Google Scholar

Stalazs, A. (2014). New records of some dipterans (Diptera: Cecidomyidae, Tephritidae) in north-eastern Lithuania. Zool. Ecol. 24, 55–57.

Google Scholar

Suarez, A. V., and Tsutsui, N. D. (2008). The evolutionary consequences of biological invasions. Mol. Ecol. 17, 351–360. doi: 10.1111/j.1365-294X.2007.03456.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsutsui, N. D., Suarez, A. V., and Grosberg, R. K. (2003). Genetic diversity, asymmetrical aggression, and recognition in a widespread invasive species. Proc. Natl. Acad. Sci. U.S.A. 100, 1078–1083. doi: 10.1073/pnas.0234412100

PubMed Abstract | CrossRef Full Text | Google Scholar

Varshney, R. K., Chabane, K., Hendre, P. S., Aggarwal, R. K., and Graner, A. (2007). Comparative assessment of EST-SSR, EST-SNP and AFLP markers for evaluation of genetic diversity and conservation of genetic resources using wild, cultivated and elite barleys. Plant Sci. 173, 638–649. doi: 10.1016/j.plantsci.2007.08.010

CrossRef Full Text | Google Scholar

Wang, G. Y. (2009). Primary Study on Biology, Ecology and Chemical Control of Obolodiplosis robiniae. Master’s dissertation, Shangdong Agriculture University, Taian.

Google Scholar

Wu, N. N., Zhang, S. F., Li, X. W., Cao, Y. H., Liu, X. J., Wang, Q. H., et al. (2019). Fall webworm genomes yield insights into rapid adaptation of invasive species. Nat. Ecol. Evol. 3, 105–115. doi: 10.1038/s41559-018-0746-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X. Q., and Yang, M. S. (2006). Review on utilization of Robinia pseudoacacia. J. Hebei For. Sci. Technol. S1:54.

Google Scholar

Yang, Z. Q., Qiao, X. R., Bu, W. J., Yao, Y. X., Xiao, Y., and Han, Y. S. (2006). First discovery of an important invasive insect pest, Obolodiplosis robiniae (Diptera: Cecidomyiidae) in China. Acta Entomol. Sin. 49, 1050–1053.

Google Scholar

Yao, Y. X., Zhao, W. X., and Shang, X. P. (2015). Development of polymorphic microsatellite markers of Obolodiplosis robiniae (Haldeman) (Diptera: Cecidomyiidae), a North American pest invading Asia. J. Insect Sci. 15:127. doi: 10.1093/jisesa/iev104

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeh, F. C., Yang, R. C., Boyle, T. B. J., Ye, Z. H., and Mao, J. X. (2018). POPGENE, the User-Friendly Shareware for Population Genetic Analysis. Available at: http://wanfang.bjast.com.cn/D/ExternalResourceswdyx200204012%5E45.aspx, (accessed November 10, 2018).

Google Scholar

Zalewski, A., Michalska-Parda, A., Bartoszewicz, M., Kozakiewicz, M., and Brzeziński, M. (2010). Multiple introductions determine the genetic structure of an invasive species population: American mink Neovison vison in Poland. Biol. Conserv. 143, 1355–1363. doi: 10.1016/j.biocon.2010.03.009

CrossRef Full Text | Google Scholar

Zeisset, I., and Beebee, T. J. C. (2003). Population genetics of a successful invader: the marsh frog Rana ridibunda in Britain. Mol. Ecol. 12, 639–646. doi: 10.1046/j.1365-294x.2003.01775.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, S. Y., Sun, S. G., Dai, C., Gituru, R. W., Chen, J. M., Wang, Q. F., et al. (2015). Genetic variation and structure in native and invasive Solidago canadensis populations. Weed Res. 55, 163–172.

Google Scholar

Zhou, Z. X., Wan, F. H., Zhang, G. F., and Chen, B. (2007). A rapid method for extraction of genomic DNA of Bemisia tabaci. Plant Protect. 33, 131–133. doi: 10.3390/genes10090632

PubMed Abstract | CrossRef Full Text | Google Scholar

Zong, J. W., Zhao, T. T., Ma, Q. H., Liang, L. S., and Wang, G. X. (2015). Assessment of genetic diversity and population genetic structure of Corylus mandshurica in China using SSR markers. PLoS One 10:e0137528. doi: 10.1371/journal.pone.0137528

PubMed Abstract | CrossRef Full Text | Google Scholar