Next Article in Journal
Phenology of the Diamondback Moth (Plutella xylostella) in the UK and Provision of Decision Support for Brassica Growers
Next Article in Special Issue
Phylogeny of the Subtribe Monoctonina (Hymenoptera, Braconidae, Aphidiinae)
Previous Article in Journal
Identification and Expression Patterns of Opsin Genes in a Forest Insect, Dendrolimus punctatus
Previous Article in Special Issue
Microsatellite Marker Discovery in the Stingless Bee Uruçu-Amarela (Melipona rufiventris Group, Hymenoptera, Meliponini) for Population Genetic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Divergence of Two Sitobion avenae Biotypes on Barley and Wheat in China

1
State Key Laboratory of Crop Stress Biology for Arid Areas, Northwest A&F University, Yangling 712100, Shaanxi, China
2
College of Plant Protection, Northwest A&F University, Yangling 712100, Shaanxi, China
3
Department of Foreign Languages, Northwest A&F University, Yangling 712100, Shaanxi, China
*
Author to whom correspondence should be addressed.
Both authors contributed equally to this work.
Insects 2020, 11(2), 117; https://doi.org/10.3390/insects11020117
Submission received: 6 January 2020 / Revised: 7 February 2020 / Accepted: 7 February 2020 / Published: 11 February 2020
(This article belongs to the Special Issue Population Genetics of Insects)

Abstract

:
Host plant affinity and geographic distance can play critical roles in the genetic divergence of insect herbivores and evolution of insect biotypes, but their relative importance in the divergence of insect populations is still poorly understood. We used microsatellite markers to test the effects of host plant species and geographic distance on divergence of two biotypes of the English grain aphid, Sitobion avenae (Fabricius). We found that clones of S. avenae from western provinces (i.e., Xinjiang, Gansu, Qinghai and Shaanxi) had significantly higher genetic diversity than those from eastern provinces (i.e., Anhui, Henan, Hubei, Zhejiang and Jiangsu), suggesting their differentiation between both areas. Based on genetic diversity and distance estimates, biotype 1 clones of eastern provinces showed high genetic divergence from those of western provinces in many cases. Western clones of S. avenae also showed higher genetic divergence among themselves than eastern clones. The Mantel test identified a significant isolation-by-distance (IBD) effect among different geographic populations of S. avenae, providing additional evidence for a critical role of geography in the genetic structure of both S. avenae biotypes. Genetic differentiation (i.e., FST) between the two biotypes was low in all provinces except Shaanxi. Surprisingly, in our analyses of molecular variance, non-significant genetic differentiation between both biotypes or between barley and wheat clones of S. avenae was identified, showing little contribution of host-plant associated differentiation to the divergence of both biotypes in this aphid. Thus, it is highly likely that the divergence of the two S. avenae biotypes involved more geographic isolation and selection of some form than host plant affinity. Our study can provide insights into understanding of genetic structure of insect populations and the divergence of insect biotypes.

1. Introduction

Many insect species are serious pests on various horticultural plants and agricultural crops. In order to develop ecologically based management programs, it is of great importance to understand the genetic diversity and population structure of insect pests. Genetic diversity and divergence among populations are essential in enabling insect species to respond rapidly to evolutionary challenges, thus having significant impacts on the adaptability of insects under constantly changing environmental conditions [1,2,3]. During the long term evolutionary process, insect species could have experienced genetic differentiation and diversification from combined effects of mutation, selection, gene flow, and genetic drift [4,5,6,7]. In this process, gene flow among populations could have been prevented by some environmental barriers, and local selection could be in turn accelerated, thus promoting changes in genetic structure of populations [6,7,8]. Demographic histories involving events of inbreeding and migration can also affect genetic structure of pest insect populations [9,10]. Many studies have shown that different geographic populations can show patterns of IBE (isolation by environment) or IBD (isolation by distance) (IBD) [11,12,13]. In addition, host plants are another significant factor influencing genetic structure of insect populations, since host-plant associated differentiation appears to be a common phenomenon for insect populations [7,14,15,16]. Indeed, various biotypes (usu. host plant associated populations) have been constantly discovered in many pest insects like the pea aphid (Acyrthosiphon pisum), Hessian fly (Mayetiola destructor), black currant leaf midge (Dasineura tetensi), and soybean aphid (Aphis glycines) [17,18,19,20].
Different insect biotypes can be distinguished by their characteristic response patterns (i.e., differential life-history traits or fitness) on different plants or different varieties of the same plant [21,22]. In particular, aphid species are prone to evolve different biotypes on variable plants because of their unique characteristics like common occurrence of local adaptation and phenotypic plasticity [23,24,25], induction of significant plant responses with saliva [26,27,28], and harboring various primary and secondary endosymbionts (e.g., Buchnera, Regiella and Hamiltonella) [29,30,31]. About half of all the insect species with known biotypes are aphids [32]. Indeed, the evolution of biotypes has been found to occur in at least 17 species of aphids, such as the pea aphid (A. pisum), greenbug (Schizaphis graminum), and Russian wheat aphid (Diuraphis noxia) [17,33,34,35]. Since the late 1990s, there have been a great number of studies that use DNA markers to examine geographic or host-plant associated differentiation of various aphids [2,36,37,38,39]. But the relative importance of geography and host plant use on the biotype divergence of aphid populations, as well as genetic relationships among different aphid biotypes, is still not well understood [22,36,37,38].
Here, we use the English grain aphid, Sitobion avenae (Fabricius) (Hemiptera: Aphididae), as a model to address this issue. This aphid can use and survive on various cereal crops and numerous wild grasses [15,40,41]. Multiple biotypes of this aphid have been discovered based on their unique response patterns on different barley and wheat varieties [22]. Biotypes 1 and 2 of S. avenae were found to be the most common in China [22]. Microsatellite markers have been widely used in many studies on population genetics of aphids, involving genotypic diversity [2,42], host range expansion [43], and host plant specialization [14,44], and identification of migration routes [2,45]. In this study, we use microsatellites to interpret the population genetics of the two predominant S. avenae biotypes, which can be critical for relevant crop breeding and integrated management programs. Specifically, our objectives are: (1) to examine genetic diversity and structure of the two S. avenae biotypes; (2) to explore genetic differentiation and relationship between both biotypes of S. avenae; and (3) to assess the importance of geography and host plant use in explaining differentiation of S. avenae biotypes.

2. Materials and Methods

2.1. Insect Sampling

Apterous adults of S. avenae were collected on wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) from nine provinces of China (Table 1). We randomly selected S. avenae samples from more than five wheat or barley fields (cultivars and pest control measures unknown) at the same location. In order to minimize the probability of collecting identical clones, apterous aphid individuals were taken at a distance of >10 m [7,46]. A minimum of 50 (up to 137) aphid clones were collected in each province during April–July 2016.

2.2. Aphid Genotyping

Six microsatellite loci (i.e., Sm10, Sm17, S17b, S5.L, S4Σ, and Sm12) were used to genotype all collected individuals of S. avenae (Supplementary Table S1) [7,47,48]. For each microsatellite locus, three primers (synthesized by Sangon Biotech Co., Ltd., Shanghai, China) were utilized, including a sequence-specific reverse primer, a sequence-specific forward primer with M13 (−21) tail added at its 5′ end, and a M13 (−21) primer fluorescent-labeled with FAM (i.e., 6-carboxy-fluorescine) [49]. The 25 μL-volume PCR reactions contained 12.5 μL 2 × Taq Master Mix (TaKaRa Biomedical Technology Co. Ltd., Beijing, China), 1.5 μL template DNA (15–25 ng/μL), 1 μL forward primer (10 μM), 2 μL reverse primer (10 μM), 1.5 μL M13 (−21) primer (10 μM), and 6.5 μL ddH2O. The PCR amplifications were conducted using a CFX-96 thermocycler (Bio-Rad Laboratories Inc., Hercules, CA, USA) with a temperature program as follows: 94 °C for 2 min, followed by 30 cycles (94 °C for 30 s, 30 s at the annealing temperature of each primer, 72 °C for 30 s), then 8 cycles (94 °C for 30 s, 53 °C for 45 s, and 72 °C for 45 s), then the last step at 72 °C for 10 min. The length of all PCR products was determined with the automated DNA sequencer ABI3730XL (Applied Biosystems, Foster City, CA, USA). Genotypes of S. avenae were distinguished with GENCLONE 2.0 [50]. Detailed information on these genotypes was included in Supplementary Materials (i.e., Table S2), including geography, host plant and allele sizes.

2.3. Biotype Identification

The identification of biotypes 1 and 2 in these provinces was conducted as described in [22]. Briefly, dozens of selected barley and wheat varieties were screened for resistance against S. avenae by using life-history bioassays. All collected S. avenae genotypes were then tested for their response profiles on resistant barley/wheat varieties identified in the above-mentioned screening life-history tests and three susceptible controls (i.e., Aikang 58, Mingxian 169, and Xinong 979). Based on these tests, it was found that five wheat/barley (i.e., wheat: Zhong 4 wumang, and 186.TM12-34; barley: Dulihuang, Zaoshu No.3, and Xiyin No.2) varieties could be used to distinguish between S. avenae biotypes. Biotype 1 was non-virulent on the five wheat/barley varieties mentioned above, whereas biotype 2 was characterized by its virulence on the barley (Hordeum vulgare L.) cultivar Zaoshu No.3.

2.4. Data Analysis

The software GenAlEx version 6.5 was utilized to analyze the following parameters: number of alleles (Na), number of effective alleles (Ne), observed heterozygosities (Ho), expected heterozygosities (He), and Shannon’s information index (I) [51]. We used the software FSTAT Version 2.9.3.2 [52] to determine and compare the estimates of gene diversity (HS) and allelic richness (AR) for biotypes 1 and 2. Using this software, Hardy–Weinberg equilibrium (HWE) was also tested for S. avenae populations. In all the above-mentioned analyses, we used data sets with only a single representative of each multi-locus genotype (MLG) per population following [2,47].
The software ARLEQUIN (version 3.5.1.2) was used to conduct the analysis of molecular variance (AMOVA) [53]. The software ARLEQUIN was also used to assess the pairwise fixation index (FST) and its significance with 9999 bootstraps. Effects of isolation by distance (IBD) were evaluated (9999 permutations) with the Mantel test in the program GenAlEx (version 6.5) [51], for which matrices for both genetic FST/(1-FST) and geographic distances between different populations of S. avenae were used. Based on Nei’s standard genetic distance (Ds), we used MEGA 5 to establish a neighbor-joining (NJ) tree of S. avenae clones [54]. Ds between clones of S. avenae was calculated by using the microsatellite analyzer (MSA) version 3.15 [55]. In order to visualize genetic variation among S. avenae clones, a principal coordinate analysis (PCoA) was conducted with GenAlEx (version 6.5) [51]. The genetic structure of S. avenae biotypes was evaluated with the Bayesian clustering program in STRUCTURE version 2.3.3 [56]. In this analysis, the number of genetic clusters (K) was set from 1 to 13; 20 replicates were performed for each K; in each replicate, a burn-in period of 10,000 iterations and 100,000 Markov chain Monte Carlo (MCMC) iterations were conducted. The most probable value of genetic clusters (K) was identified with the Evanno method in STRUCTURE HARVESTER [57,58]. The genetic assignment software BayesAss (version 1.3) was used to examine the direction and magnitude (m, the proportion of migrant individuals in each population) of gene flow between S. aveane clones [45]. Based on the posterior probability distribution, the proportion of immigrants from one population to another can be determined with a fully Bayesian MCMC (Markov chain Monte Carlo) resampling method [2,45].

3. Results

3.1. Geographic Distribution and Genetic Diversity of Biotypes 1 and 2

Cones of biotype 1 were distributed over all provinces sampled in this study, whereas those of biotype 2 occurred only in western provinces (i.e., Shaanxi, Qinghai, Gansu, and Xinjiang) (Figure 1). The numbers of genotypes identified for biotype 1 in each province ranged from eight to 49 (Table 2), and those for biotype 2 ranged from five to six in each province, indicating higher genotypic diversity for biotype1 compared with biotype 2.
Mean numbers of alleles (Na) ranged from 1.6 for biotype 2 collected from Gansu (GS2) to 11.33 for biotype 1 of Qinghai (QH1) (Table 2), whereas the number of effective alleles (Ne) varied from 0.9 for GS2 to 4.8 for biotype 1 of Shaanxi (SX1). Gene diversity (Hs) was lowest (i.e., 0.651) for biotype 1 of Anhui (AH1), and highest (i.e., 0.801) for SX1. Allelic richness (AR) varied from 3.74 for biotype 1 of Xinjiang (XJ1) to 5.26 for biotype 1 of Jiangsu (JS1). The Shannon’s information index (I) ranged from 1.01 for GS2 to 1.78 for QH1. The values of Ho and He fell in the range of 0.06 for GS2 to 0.77 for biotype 2 of Shaanxi (SX2), and 0.06 for GS2 to 0.75 for QH1, respectively. Based on all six microsatellite loci, significant deviation from HWE (HWE—p < 0.05) was found for all geographic populations but Anhui and Henan (Table S3).
Biotype 1 showed a significantly higher I (1.51 ± 0.06) than biotype 2 (1.01 ± 0.17) (Figure 2B; t = 2.542, p = 0.027), indicating higher genetic diversity for biotype 1 compared to biotype 2. A significantly higher I was also detected for barley clones (1.48 ± 0.11) compared with wheat clones (1.20 ± 0.06) (Figure 2F; t = 2.402, p = 0.036). Aphid clones of western provinces (0.76 ± 0.01) showed a higher gene diversity (Hs) than those of eastern provinces (0.69 ± 0.02) (Figure 2C; t = 2.977, p = 0.024). This result might be in part attributed to larger sample size in western provinces than eastern provinces. Biotype 1 on both plants (i.e., wheat and barley), and biotype 2 on barley showed significant deviation from HWE (Table S4). Biotype 2 on wheat did not show significant departure from HWE.

3.2. Genetic Differentiation between S. avenae Biotypes

Based on pairwise FST values, high genetic divergence (FST > 0.15) between the two biotypes were detected in some cases, such as between AH1 (biotype 1 of Anhui) and GS2 (biotype 2 of Gansu), between HB1 (biotype 1 of Hubei) and GS2, between HN1 (biotype 1 of Henan) and GS2, between ZJ1 (biotype 1 of Zhejiang) and GS2, between GS1 (biotype 1 of Gansu) and SX2 (biotypes 2 of Shaanxi), and between XJ1 (biotype 1 of Xinjiang) and SX2 (Table 3). Indices of genetic differentiation (i.e., FST) between the two biotypes in each province ranged from low to moderate (FST = 0.001–0.078). Moderate differentiation (0.05 < FST ≤ 0.15) between the two biotypes was found only in Shaanxi (FST = 0.078).
A neighbor-joining tree based on Nei’s standard genetic distance (Ds) was created to show the genetic relationships of both S. avenae biotypes from nine locations (Figure 3). Clones of S. avenae formed four clades. One major clade consisted of biotype 1 of eastern provinces (i.e., Zhejiang, Hubei, Henan, Anhui, and Jiangsu). Another major clade included both biotypes from Xinjiang (XJ1 and XJ2) and Gansu (GS1 and GS2), as well as biotype 1 of Qinghai. Biotype 1 of Shaanxi and biotype 2 of Qinghai formed a third clade. Biotype 2 of Shaanxi alone formed a separate clade.
Principal coordinate analyses (PCoA) showed that the first two axes accounted for 76.76% of the total variation (Figure 4; 61.82% and 14.94% for PC1 and PC2, respectively). The PCoA plot showed biotype 1 from eastern provinces of Zhejiang, Jiangsu, Anhui, Henan, and Hubei clustered in the lower left quadrant. Biotype 2 of Shaanxi (SX2) alone fell in the upper right quadrant. All other S. avenae clones fell in the lower right quadrant, and they included biotype 1 of Shaanxi (SX1), and both biotypes of Xinjiang (XJ1 and XJ2), Gansu (GS1 and GS2), and Qinghai (QH1 and QH2).
Based on the Bayesian analysis, the best supported model was K = 2 (Figure 5A), indicating that all clones of S. avenae in this study could be grouped into two genetic clusters. The majority of S. avenae clones for biotype 1 of Henan (98.0%), Hubei (96.2%), Anhui (95%), Jiangsu (96.6%) and Zhejiang (98.5%) belonged to cluster 1 (Figure 5B). On the contrary, most clones belonged to cluster 2 for western provinces with the proportion of these clones varying from 60.8% for biotype 2 of Shaanxi to 98.8% for biotype 1 of Xinjiang. For these western provinces, the genetic structure of biotype 1 clones seemed to be different from that of biotype 2, especially for Xinjiang (biotype 1: 1.2% cluster 1 and 98.8% cluster 2; biotype 2: 22.8% cluster 1 and 77.2% cluster 2) and Shaanxi (biotype 1: 24.3% cluster 1 and 75.7% cluster 2; biotype 2: 39.2% cluster 1 and 60.8% cluster 2).
A significantly positive relationship was identified between genetic and geographic distances (Mantel test, r = 0.696, p < 0.001), indicating a significant effect of IBD (isolation-by-distance) among geographic populations of S. avenae (Figure 6).
AMOVA analyses showed that only 0.26% of the overall molecular variation was explained by variation between the two biotypes (Table 4; p = 0.822), showing no significant genetic differentiation between biotypes 1 and 2 of S. avenae. Another AMOVA analysis showed that 9.70% of the overall molecular variation could be explained by geographic location (p < 0.001), suggesting significant variation between S. avenae clones from eastern (i.e., Anhui, Hubei, Henan, Jiangsu, and Zhejiang) and western (i.e., Gansu, Qinghai, Shaanxi, and Xinjiang) provinces. AMOVA analyses of the host plant effect showed that it could contribute to only 0.48% (p = 0.950) of the overall molecular variation, suggesting that the overall variation between wheat and barley clones of S. avenae was non-significant in our sampling areas of this study.

3.3. Gene Flow

Substantial gene flow (m) was found between the two biotypes (i.e., biotypes 1 and 2), as well as between eastern and western provinces (Table 5). In western provinces, the level of gene flow was found to be high (m = 0.2227) from biotype 1 to biotype 2, but it was non-significant in the reverse direction (i.e., from biotype 2 to biotype 1) (m = 0.0233). The migration of S. avenae toward eastern provinces was very low from both biotypes (biotype 1: m = 0.0053; biotype 2: m = 0.0487) of western provinces. Low rates of migration were also found from biotype 1 of eastern provinces to biotype 1 (m = 0.0258) and biotype 2 (m = 0.0823) of western provinces.
BayesAss analysis was also used to estimate gene flow between the two biotypes from barley and wheat (Table 6). Significant gene flow was identified between biotype 1 of barley and biotype 1 of wheat in both directions (from wheat to barley: 0.1784; from barley to wheat: 0.1205). Gene flow from biotype 1 of barley to biotype 2 of wheat was also significant (m = 0.1322), but it was negligible in the opposite direction (m = 0.0148). The level of gene flow from biotype 1 of wheat to biotype 2 of barley was relatively high (m = 0.2291). On wheat, significant gene flow was found from biotype 1 to biotype 2 (m = 0.1317), but not for the opposite direction. On barley, the level of gene flow between biotypes 1 and 2 was low in both directions (from biotype 1 to biotype 2: 0.0521; from biotype 2 to biotype 1: 0.0146).

4. Discussion

4.1. Geographic Divergence among S. avenae Clones

Geographic structure and isolation of populations can play critical roles in the diversification of herbivorous insect populations. Despite that, a few studies have focused on the genetic differentiation of S. avenae clones in different areas of China [2,6,59], the relationship between geographic divergence and biotype development in S. avenae is still not well understood. In this study, in addition to differences in pairwise FST values, clones of S. avenae from western provinces had significantly higher gene diversity than those from eastern provinces, indicating differentiation between both areas. In addition, based on other genetic parameters like number of effective alleles (Ne), allelic richness (R), and Shannon’s information index (I), clones of western provinces for biotype 1 were also shown to have much higher genetic diversity than those of eastern provinces for the same biotype. It is generally believed that genetic differentiation between populations is low to negligible when FST ≤ 0.05, moderate when 0.05 < FST < 0.15, and high when FST ≥ 0.15 [2,59,60]. In this study, S. avenae clones from eastern provinces, which were all found to be biotype 1, showed little genetic divergence with pairwise FST values ranging from 0.010 to 0.054. However, some clones from western provinces showed high levels of genetic differentiation among them. For example, the pairwise FST value between XJ1 (biotype 1 of Xinjiang) and SX2 (biotype 2 of Shaanxi) was 0.221, and it was 0.186 between SX2 and GS2 (biotype 2 of Gansu). Clones from western provinces showed more scattered distribution in the PCoA plot than those from eastern provinces, suggesting higher levels of genetic divergence for clones from western provinces. Biotype 1 clones of both Xinjiang and Gansu showed high genetic differentiation from those of all five eastern provinces. A similar pattern was also found between GS2 and biotype 1 clones of eastern provinces. In addition to estimates of genetic diversity and distance (i.e., pairwise FST), analyses of genetic structure of S. avenae clones and AMOVA (analyses of molecular variance) also indicated significant genetic divergence between eastern and western populations of this study. Significant geographic divergence of S. avenae populations could be attributed to limited long-distance dispersal from eastern to western provinces. Indeed, movements of biotype 1 of S. avenae were found to be very limited between eastern and western provinces in both directions. In China, the Qinling Mountains are clearly an important genetic barrier between populations to the east or west of the mountains [6]. The Mantel test in this study showed a significant positive relationship between genetic and geographic distances (r = 0.696, p < 0.001), indicating a significant isolation-by-distance (IBD) effect among different geographic populations of S. avenae. Our results indicate that geographic distance can play a significant role in the genetic structure of both S. avenae biotypes.
High genetic differentiation between S. avenae clones from eastern and western provinces can be closely related to geographic differences in their life cycles. In Europe, a majority of holocyclic clones of this aphid were discovered at locations with harsh winters (e.g., Romania), but more clones of obligate parthenogenesis were found in places with mild winters (e.g., southern France) [47,61]. According to the study of [62], overwintering with eggs and sexual reproduction of S. avenae in China could occur in western provinces, but probably not in all the eastern provinces of this study except Jiangsu. Thus, one possible explanation for high genetic differentiation between S. avenae clones from eastern and western provinces in China was that western provinces had lower winter temperatures, which could lead to more occurrences of sexual reproduction (recombination of alleles) for S. avenae in western provinces than in eastern provinces [63,64,65,66,67]. In this study, biotype 2 occurred in all four western provinces, instead of five eastern provinces. Thus, biotype 2 could be more likely to evolve from sexual lineages than asexual lineages. Differences in proportions of asexual and sexual lineages of S. avenae in both areas might contribute significantly to the abovementioned divergence and biotype development of this aphid. Sexual and asexual clones of S. avenae can be distinguished by rearing them under short day conditions in the lab [68]. Further studies in this respect can help to clarify the effects of reproductive mode on population differentiation of this aphid. Other geographic factors, such as selection from local natural enemies, adaptation to different environments (i.e., isolation by environment), and local crop domestication and agriculture, could also influence the divergence of S. avenae clones that we did not include in this study [69,70,71]. Further studies are needed to determine the effects of these specific geographic factors (e.g., relative proportions of sexual and asexual lineages of S. avenae in different areas) on the biotype development of S. avenae.

4.2. Genetic Diversity and Divergence for S. avenae Biotypes

Host plant use and diet breadth are often assumed to be the critical factors influencing population divergence of insect herbivores. This is especially true for different aphid species, which are prone to develop variable host plant associated populations (e.g., biotypes) [17,20,33,34,35]. In our previous study, biotype 1 of the English grain aphid (S. avenae) was non-virulent on barley, whereas biotype 2 was characterized by its high virulence on barley (e.g., Zaoshu No.3) [22], clearly showing the potential effects of host plants on the divergence of S. aveane clones. In the present study, clones of biotype 1 in western provinces had higher genetic diversity than those of biotype 2 in the same area, based on estimates like number of effective alleles (Ne), allelic richness (R), and Shannon’s information index (I). When all samples were combined, S. avenae biotype 1 showed a higher Shannon’s information index (I) than biotype 2 (Figure 2B). In addition, in western provinces, the genetic structure of biotype 1 clones seemed to be different from that of biotype 2, especially for Xinjiang (biotype 1: 1.2% cluster 1 and 98.8% cluster 2; biotype 2: 22.8% cluster 1 and 77.2% cluster 2) and Shaanxi (biotype 1: 24.3% cluster 1 and 75.7% cluster 2; biotype 2: 39.2% cluster 1 and 60.8% cluster 2). These results indicated that there could be some degree of genetic divergence between both biotypes in relation to host plant use.
However, the genetic divergence between the two biotypes (i.e., biotypes 1 and 2) of S. avenae was not significant in our analyses of molecular variance (AMOVA) with samples from nine provinces (Table 4). Genetic differentiation between barley and wheat clones of S. avenae in our sampling areas of nine provinces in China was not statistically significant either. Low genetic differentiation between both biotypes or between barley and wheat clones of S. avenae might be attributed to significant dispersal between these clones. For example, the levels of gene flow were found to be significant between biotype 1 of barley and biotype 1 of wheat in both directions (from wheat to barley: 0.1784; from barley to wheat: 0.1205). Similarly, substantial gene flow was also found between wheat and barley clones in our previous study [7]. Despite this, based on pairwise FST values and genetic structure analyses, significant levels of genetic divergence between barley and wheat clones were identified in Jiangsu and Zhejiang in our previous study [7]. This inconsistency may be due to different geographic scales used in both studies. Our previous study was limited to four eastern provinces (i.e., Jiangsu, Zhejiang, Hubei and Henan). Patterns of genetic divergence can strongly depend on the geographic scale being studied, because geography can strongly influence levels of gene flow in addition to environment and ecological contexts (i.e., host plant distribution) in a specific area [72]. Thus, geographic factors appeared to be more important than host plant use in the divergence of the two S. avenae biotypes. Some studies have proposed several mechanisms that may promote the maintenance of genetic structure of insect geographic populations and evolution of insect biotypes, such as differential composition of secondary endosymbionts, phenotypic plasticity, habitat isolation, and habitat persistence [2,21,24,30,31,73]. Future studies are needed to determine if these mechanisms can explain the development of biotypes in S. avenae.

5. Conclusions

In summary, host plant use and geographic distance are some of the factors that can influence the genetic divergence of insect herbivores and development of insect biotypes. So far, the relative importance of both factors in the divergence of S. avenae clones is still not well understood. In this study, in terms of comparisons on genetic diversity between groups of S. avenae clones (e.g., eastern vs. western areas, and biotype 1 vs. biotype 2), the results of Nei’s gene diversity (Hs) were different from those of Shannon information index (I). The meaning of these differences is still difficult to interpret [59,74,75,76], which warrants further studies in the future. Surprisingly, we found that host-plant associated differentiation seemed to contribute little to the divergence of biotypes in S. avenae, since the plant factor contributed little to the total genetic differentiation of different clones of this aphid. In our previous study limited in eastern provinces, we found evidence for significant differentiation of barley and wheat clones of S. avenae [7]. Therefore, patterns of genetic divergence for this aphid can strongly depend on the geographic scale being studied. In addition, two biotypes of S. avenae showed significant structure and differentiation between eastern and western populations, and biotype 2 only occurred in western provinces. Therefore, it is very likely that the divergence of the two S. avenae biotypes [22] involved more geographic isolation and selection of some form than host plant affinity. The greater influence of geographic distance than host plant use on S. avenae genetic structure suggests that levels of host plant specialization in this aphid are still low, or that host shift is still recent in the evolutionary history, consistent with our previous life-history studies on this aphid [14,15]. Although host plant affinity is frequently identified to be the basis for evolution of biotypes, our data suggest that geographic factors, instead of host plant use, can be more important in biotype development in some insects (especially for those with wide-area distributions). Further studies are needed to identify key specific geographic elements involved in the process, and the underlying mechanisms. Genome scans, as well as candidate-gene studies, are also needed to identify which genes in insects have evolved in response to geographically variable selection.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/2/117/s1.

Author Contributions

Designed research, D.L. and D.W.; performed research, D.W., X.S., and Y.Y.; writing, review and editing, D.W., X.S., D.L. and Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (31572002), and the State Key Laboratory of Integrated Management of Pest Insects and Rodents (Grant No. IPM1909).

Acknowledgments

We are indebted to Ying-Ting Zhai (Northwest A&F University) for her laboratory and field assistance.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Xie, P.; Zhao, G.; Niu, J.G.; Wang, J.; Zhou, Q.; Guo, Y.; Ma, X.F. Comprehensive analysis of population genetics of Phoxinus phoxinus ujmonensis in the Irtysh River: Abiotic and biotic factors. Ecol. Evol. 2019, 9, 7997–8012. [Google Scholar] [CrossRef] [PubMed]
  2. He, Y.; Liu, D.; Dai, P.; Wang, D.; Shi, X. Genetic differentiation and structure of Sitobion avenae (Hemiptera: Aphididae) populations from moist, semiarid and arid areas in Northwestern China. J. Econ. Entomol. 2018, 111, 603–611. [Google Scholar] [CrossRef] [PubMed]
  3. Cobben, M.M.; Verboom, J.; Opdam, P.F.; Hoekstra, R.F.; Jochem, R.; Smulders, M.J. Landscape prerequisites for the survival of a modelled metapopulation and its neutral genetic diversity are affected by climate change. Landsc. Ecol. 2012, 27, 227–237. [Google Scholar] [CrossRef]
  4. Via, S.; Conte, G.; Mason-Foley, C.; Mills, K. Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation-with-gene-flow. Mol. Ecol. 2012, 21, 5546–5560. [Google Scholar] [CrossRef] [PubMed]
  5. Lang, M.; Murat, S.; Clark, A.G.; Gouppil, G.; Blais, C.; Matzkin, L.M.; Guittard, É.; Yoshiyama-Yanagawa, T.; Kataoka, H.; Niwa, R. Mutations in the neverland gene turned Drosophila pachea into an obligate specialist species. Science 2012, 337, 1658–1661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Huang, X.; Liu, D.; Wang, D.; Shi, X.; Simon, J.-C. Molecular and quantitative genetic differentiation in Sitobion avenae populations from both sides of the Qinling Mountains. PLoS ONE 2015, 10, e0122343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Wang, D.; Liu, D.; Zhai, Y.; Zhang, R.; Shi, X. Clonal diversity and genetic differentiation of Sitobion avenae (Hemiptera: Aphididae) from wheat and barley in China. J. Econ. Entomol. 2019, 112, 1217–1226. [Google Scholar] [CrossRef]
  8. Nosil, P. Speciation with gene flow could be common. Mol. Ecol. 2008, 17, 2103–2106. [Google Scholar] [CrossRef]
  9. De Barba, M.; Waits, L.; Garton, E.; Genovesi, P.; Randi, E.; Mustoni, A.; Groff, C. The power of genetic monitoring for studying demography, ecology and genetics of a reintroduced brown bear population. Mol. Ecol. 2010, 19, 3938–3951. [Google Scholar] [CrossRef]
  10. McKinnon, J.S.; Rundle, H.D. Speciation in nature: The threespine stickleback model systems. Trends Ecol. Evol. 2002, 17, 480–488. [Google Scholar] [CrossRef]
  11. Collin, H.; Fumagalli, L. Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae). Mol. Ecol. 2011, 20, 4490–4502. [Google Scholar] [CrossRef] [PubMed]
  12. Orsini, L.; Vanoverbeke, J.; Swillen, I.; Mergeay, J.; De Meester, L. Drivers of population genetic differentiation in the wild: Isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol. Ecol. 2013, 22, 5983–5999. [Google Scholar] [CrossRef] [PubMed]
  13. Ruiz-Gonzalez, A.; Cushman, S.A.; Madeira, M.J.; Randi, E.; Gómez-Moliner, B.J. Isolation by distance, resistance and/or clusters? Lessons learned from a forest-dwelling carnivore inhabiting a heterogeneous landscape. Mol. Ecol. 2015, 24, 5110–5129. [Google Scholar] [CrossRef] [PubMed]
  14. Gao, S.X.; Liu, D.G.; Chen, H.; Meng, X.X. Fitness traits and underlying genetic variation related to host plant specialization in the aphid Sitobion avenae. Insect Sci. 2014, 21, 352–362. [Google Scholar] [CrossRef] [PubMed]
  15. Gao, S.; Liu, D. Differential performance of Sitobion avenae (Hemiptera: Aphididae) clones from wheat and barley with implications for its management through alternative cultural practices. J. Econ. Entomol. 2013, 106, 1294–1301. [Google Scholar] [CrossRef]
  16. Simon, J.-C.; Mahéo, F.; Mieuzet, L.; Buchard, C.; Gauthier, J.-P.; Maurice, D.; Bonhomme, J.; Outreman, Y.; Hullé, M. Life on the Edge: Ecological Genetics of a High Arctic Insect Species and Its Circumpolar Counterpart. Insects 2019, 10, 427. [Google Scholar] [CrossRef] [Green Version]
  17. Cartier, J.J. Recognition of three biotypes of the pea aphid from southern Quebec. J. Econ. Entomol. 1959, 52, 293–294. [Google Scholar] [CrossRef]
  18. Ratcliffe, R.H.; Cambron, S.E.; Flanders, K.L.; Bosque-Perez, N.A.; Clement, S.L.; Ohm, H.W. Biotype composition of Hessian fly (Diptera: Cecidomyiidae) populations from the southeastern, midwestern, and northwestern United States and virulence to resistance genes in wheat. J. Econ. Entomol. 2000, 94, 1319–1328. [Google Scholar] [CrossRef]
  19. Hellqvist, S. Biotypes of Dasineura tetensi, differing in ability to gall and develop on black currant genotypes. Entomol. Exp. Appl. 2001, 98, 85–94. [Google Scholar] [CrossRef]
  20. Kim, K.-S.; Hill, C.B.; Hartman, G.L.; Mian, M.A.R.; Diers, B.W. Discovery of Soybean Aphid Biotypes. Crop Sci. 2008, 48, 923–928. [Google Scholar] [CrossRef]
  21. Bansal, R.; Michel, A. Molecular Adaptations of Aphid Biotypes in Overcoming Host-Plant Resistance. In Short Views on Insect Genomics and Proteomics; Springer: Berlin/Heidelberg, Germany, 2015; pp. 75–93. [Google Scholar]
  22. Wang, D.; Zhai, Y.; Liu, D.; Zhang, N.; Li, C.; Shi, X. Identification and genetic differentiation of Sitobion avenae (Hemiptera: Aphididae) biotypes in China. J. Econ. Entomol. 2019. [Google Scholar] [CrossRef] [PubMed]
  23. Huang, X.; Liu, D.; Gao, S.; Chen, H. Differential performance of Sitobion avenae populations from both sides of the Qinling Mountains under common garden conditions. Environ. Entomol. 2013, 42, 1174–1183. [Google Scholar] [CrossRef] [PubMed]
  24. Dai, X.; Gao, S.; Liu, D. Genetic basis and selection for life-history trait plasticity on alternative host plants for the cereal aphid Sitobion avenae. PLoS ONE 2014, 9, e106179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Moran, N.A. The evolution of aphid life cycles. Annu. Rev. Entomol. 1992, 37, 321–348. [Google Scholar] [CrossRef]
  26. Dreyer, D.; Campbell, B. Association of the degree of methylation of intercellular pectin with plant resistance to aphids and with induction of aphid biotypes. Experientia 1984, 40, 224–226. [Google Scholar] [CrossRef]
  27. Lapitan, N.L.; Li, Y.-C.; Peng, J.; Botha, A.-M. Fractionated extracts of Russian wheat aphid eliciting defense responses in wheat. J. Econ. Entomol. 2007, 100, 990–999. [Google Scholar] [CrossRef]
  28. Mutti, N.S.; Louis, J.; Pappan, L.K.; Pappan, K.; Begum, K.; Chen, M.-S.; Park, Y.; Dittmer, N.; Marshall, J.; Reese, J.C. A protein from the salivary glands of the pea aphid, Acyrthosiphon pisum, is essential in feeding on a host plant. Proc. Nat. Acad. Sci. USA 2008, 105, 9965–9969. [Google Scholar] [CrossRef] [Green Version]
  29. Moran, N.A.; Wernegreen, J.J. Lifestyle evolution in symbiotic bacteria: Insights from genomics. Trends Ecol. Evol. 2000, 15, 321–326. [Google Scholar] [CrossRef]
  30. Oliver, K.M.; Degnan, P.H.; Burke, G.R.; Moran, N.A. Facultative symbionts in aphids and the horizontal transfer of ecologically important traits. Annu. Rev. Entomol. 2010, 55, 247–266. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, D.; Shi, X.; Dai, P.; Liu, D.; Dai, X.; Shang, Z.; Ge, Z.; Meng, X. Comparison of fitness traits and their plasticity on multiple plants for Sitobion avenae infected and cured of a secondary endosymbiont. Sci. Rep. 2016, 6, 23177. [Google Scholar] [CrossRef] [Green Version]
  32. Saxena, R.; Barrion, A. Biotypes of insect pests of agricultural crops. Int. J. Trop. Insect Sci. 1987, 8, 453–458. [Google Scholar] [CrossRef]
  33. Wood, E., Jr. Biological studies of a new greenbug biotype. J. Econ. Entomol. 1961, 54, 1171–1173. [Google Scholar] [CrossRef]
  34. Puterka, G.; Burd, J.; Burton, R. Biotypic variation in a worldwide collection of Russian wheat aphid (Homoptera: Aphididae). J. Econ. Entomol. 1992, 85, 1497–1506. [Google Scholar] [CrossRef]
  35. Jankielsohn, A. Distribution and diversity of Russian wheat aphid (Hemiptera: Aphididae) biotypes in South Africa and Lesotho. J. Econ. Entomol. 2011, 104, 1736–1741. [Google Scholar] [CrossRef]
  36. Sunnucks, P.; De Barro, P.; 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–1073. [Google Scholar] [CrossRef]
  37. Shufran, K.; Burd, J.; Anstead, J.; Lushai, G. Mitochondrial DNA sequence divergence among greenbug (Homoptera: Aphididae) biotypes: Evidence for host-adapted races. Insect Mol. Biol. 2000, 9, 179–184. [Google Scholar] [CrossRef]
  38. Vorwerk, S.; Forneck, A. Reproductive mode of grape phylloxera (Daktulosphaira vitifoliae, Homoptera: Phylloxeridae) in Europe: Molecular evidence for predominantly asexual populations and a lack of gene flow between them. Genome 2006, 49, 678–687. [Google Scholar] [CrossRef]
  39. Duan, X.; Peng, X.; Qiao, X.; Chen, M. Life cycle and population genetics of bird cherry-oat aphids Rhopalosiphum padi in China: An important pest on wheat crops. J. Pest Sci. 2017, 90, 103–116. [Google Scholar] [CrossRef]
  40. Blackman, R.L.; Eastop, V.F. Aphids on the World’s Crops: An Identification and Information Guide; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2000. [Google Scholar]
  41. Yang, Y.; Kloos, S.; Mora-Ramírez, I.; Romeis, J.; Brunner, S.; Li, Y.; Meissle, M. Transgenic Winter Wheat Expressing the Sucrose Transporter HvSUT1 from Barley does not Affect Aphid Performance. Insects 2019, 10, 388. [Google Scholar] [CrossRef] [Green Version]
  42. Kanbe, T.; AKIMOTO, S.I. Allelic and genotypic diversity in long-term asexual populations of the pea aphid, Acyrthosiphon pisum in comparison with sexual populations. Mol. Ecol. 2009, 18, 801–816. [Google Scholar] [CrossRef]
  43. Peccoud, J.; Figueroa, C.; Silva, A.; Ramirez, C.; Mieuzet, L.; Bonhomme, J.; Stoeckel, S.; Plantegenest, M.; Simon, J.C. Host range expansion of an introduced insect pest through multiple colonizations of specialized clones. Mol. Ecol. 2008, 17, 4608–4618. [Google Scholar] [CrossRef]
  44. Carletto, J.; Lombaert, E.; Chavigny, P.; Brévault, T.; Lapchin, L.; Vanlerberghe-Masutti, F. Ecological specialization of the aphid Aphis gossypii Glover on cultivated host plants. Mol. Ecol. 2009, 18, 2198–2212. [Google Scholar] [CrossRef]
  45. Wilson, G.A.; Rannala, B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics 2003, 163, 1177–1191. [Google Scholar]
  46. Ahmed, S.S.; Liu, D.; Simon, J.-C. Impact of water-deficit stress on tritrophic interactions in a wheat-aphid-parasitoid system. PLoS ONE 2017, 12, e0186599. [Google Scholar] [CrossRef] [Green Version]
  47. Simon, J.C.; Baumann, S.; Sunnucks, P.; Hebert, P.D.; Pierre, J.S.; Le Gallic, J.F.; Dedryver, C.A. Reproductive mode and population genetic structure of the cereal aphid Sitobion avenae studied using phenotypic and microsatellite markers. Mol. Ecol. 1999, 8, 531–545. [Google Scholar] [CrossRef]
  48. Wilson, A.C.; Massonnet, B.; Simon, J.C.; Prunier-Leterme, N.; Dolatti, L.; Llewellyn, K.S.; Figueroa, C.C.; Ramirez, C.C.; Blackman, R.L.; Estoup, A. Cross-species amplification of microsatellite loci in aphids: Assessment and application. Mol. Ecol. Notes 2004, 4, 104–109. [Google Scholar] [CrossRef]
  49. Schuelke, M. An economic method for the fluorescent labeling of PCR fragments. Nat. Biotechnol. 2000, 18, 233. [Google Scholar] [CrossRef]
  50. Arnaud-Haond, S.; Belkhir, K. GENCLONE: A computer program to analyse genotypic data, test for clonality and describe spatial clonal organization. Mol. Ecol. Notes 2007, 7, 15–17. [Google Scholar] [CrossRef]
  51. Peakall, R.; Smouse, P. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and researchdan update. Bioinformatics 2012, 28, 2537e2539. [Google Scholar] [CrossRef] [Green Version]
  52. Goudet, J. FSTAT (version 1.2): A computer program to calculate F-statistics. J. Hered. 1995, 86, 485–486. [Google Scholar] [CrossRef]
  53. Weir, B.S.; Cockerham, C.C. Estimating F-statistics for the analysis of population structure. Evolution 1984, 38, 1358–1370. [Google Scholar]
  54. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  55. Dieringer, D.; Schlötterer, C. Microsatellite analyser (MSA): A platform independent analysis tool for large microsatellite data sets. Mol. Ecol. Notes 2003, 3, 167–169. [Google Scholar] [CrossRef]
  56. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar]
  57. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  58. Earl, D.A. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  59. Xin, J.-J.; Shang, Q.-L.; Desneux, N.; Gao, X.-W. Genetic diversity of Sitobion avenae (Homoptera: Aphididae) populations from different geographic regions in China. PLoS ONE 2014, 9, e109349. [Google Scholar] [CrossRef] [Green Version]
  60. Wright, S. Variability within and Among Natural Populations Evolution and the Genetics of Populations: A Treatise in Four Volumes; The University of Chicago Press: Chicago, IL, USA, 1988. [Google Scholar]
  61. Papura, D.; Simon, J.; Halkett, F.; Delmotte, F.; Le Gallic, J.; Dedryver, C. Predominance of sexual reproduction in Romanian populations of the aphid Sitobion avenae inferred from phenotypic and genetic structure. Heredity 2003, 90, 397. [Google Scholar] [CrossRef] [Green Version]
  62. Yang, Y.; Lu, Z.; Han, J.; Chen, X. Overwintering eggs of Sitobion avenae in Yangzhou area. Entomol. Knowl. 1994, 31, 207–209. [Google Scholar]
  63. Dedryver, C.-A.; Hullé, M.; Le Gallic, J.-F.; Caillaud, M.C.; Simon, J.-C. Coexistence in space and time of sexual and asexual populations of the cereal aphid Sitobion avenae. Oecologia 2001, 128, 379–388. [Google Scholar] [CrossRef]
  64. Delmotte, F.; Leterme, N.; Gauthier, J.P.; Rispe, C.; Simon, J.C. Genetic architecture of sexual and asexual populations of the aphid Rhopalosiphum padi based on allozyme and microsatellite markers. Mol. Ecol. 2002, 11, 711–723. [Google Scholar] [CrossRef] [PubMed]
  65. Balloux, F.; Lehmann, L.; de Meeûs, T. The population genetics of clonal and partially clonal diploids. Genetics 2003, 164, 1635–1644. [Google Scholar] [PubMed]
  66. Bengtsson, B.O. Genetic variation in organisms with sexual and asexual reproduction. J. Evol. Biol. 2003, 16, 189–199. [Google Scholar] [CrossRef] [PubMed]
  67. Islam, M.S.; Roush, T.L.; Walker, M.A.; Granett, J.; Lin, H. Reproductive mode and fine-scale population genetic structure of grape phylloxera (Daktulosphaira vitifoliae) in a viticultural area in California. BMC Genetics 2013, 14, 123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Dedryver, C.A.; Le Gallic, J.F.; Mahéo, F.; Simon, J.C.; Dedryver, F. The genetics of obligate parthenogenesis in an aphid species and its consequences for the maintenance of alternative reproductive modes. Heredity 2013, 110, 39–45. [Google Scholar] [CrossRef] [Green Version]
  69. Heath, J.J.; Abbot, P.; Stireman, J.O., III. Adaptive Divergence in a Defense Symbiosis Driven from the Top Down. Am. Nat. 2018, 192, E21–E36. [Google Scholar] [CrossRef]
  70. Wang, I.J.; Bradburd, G.S. Isolation by environment. Mol. Ecol. 2014, 23, 5649–5662. [Google Scholar] [CrossRef]
  71. Ramirez-Romero, R.; Garibay-Benítez, D.; Vargas-Ponce, O.; Joyce, A.; Bernal, J.S. Do assortative mating and immigrant inviability help maintain population genetic structuring of an herbivore on a crop and a wild relative? Insect Sci. 2019, 26, 283–296. [Google Scholar] [CrossRef]
  72. Simon, J.-C.; Rispe, C.; Sunnucks, P. Ecology and evolution of sex in aphids. Trends Ecol. Evol. 2002, 17, 34–39. [Google Scholar] [CrossRef]
  73. Vidal, M.C.; Quinn, T.W.; Stireman, J.O.; Tinghitella, R.M.; Murphy, S.M. Geography is more important than host plant use for the population genetic structure of a generalist insect herbivore. Mol. Ecol. 2019, 28, 4317–4334. [Google Scholar] [CrossRef]
  74. Bhattacharyya, P.; van Staden, J. Molecular insights into genetic diversity and population dynamics of five medicinal Eulophia species: A threatened orchid taxa of Africa. Physiol. Mol. Biol. Plants 2018, 24, 631–641. [Google Scholar] [CrossRef] [PubMed]
  75. Tiwari, V.; Meena, B.; Nair, K.N.; Upreti, D.K.; Tamta, S.; Rana, T.S. Assessment of genetic diversity and population structure of Bergenia stracheyi (Saxifragaceae) in the western Himalaya (India). Biochem. Syst. Ecol. 2017, 70, 205–210. [Google Scholar] [CrossRef]
  76. Hennink, S.; Zeven, A.C. The interpretation of Nei and Shannon-Weaver within population variation indices. Euphytica 1991, 51, 235–240. [Google Scholar] [CrossRef]
Figure 1. Collection of Sitobion avenae biotypes 1 and 2 from nine provinces of China.
Figure 1. Collection of Sitobion avenae biotypes 1 and 2 from nine provinces of China.
Insects 11 00117 g001
Figure 2. Comparisons of estimates for gene diversity (Hs) and Shannon’s information index (I). (A,B): between biotypes 1 and 2; (C,D): between eastern and western provinces; (E,F): between wheat and barley clones; eastern provinces include Anhui, Hubei, Henan, Jiangsu, and Zhejiang; western provinces include Gansu, Qinghai, Shaanxi, and Xinjiang; (*, significant differences at the p < 0.05 level; NS, non-significant).
Figure 2. Comparisons of estimates for gene diversity (Hs) and Shannon’s information index (I). (A,B): between biotypes 1 and 2; (C,D): between eastern and western provinces; (E,F): between wheat and barley clones; eastern provinces include Anhui, Hubei, Henan, Jiangsu, and Zhejiang; western provinces include Gansu, Qinghai, Shaanxi, and Xinjiang; (*, significant differences at the p < 0.05 level; NS, non-significant).
Insects 11 00117 g002
Figure 3. A neighbor-joining dendrogram of Sitobion avenae biotypes from nine provinces (AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively; E and W in brackets indicate eastern and western provinces, respectively).
Figure 3. A neighbor-joining dendrogram of Sitobion avenae biotypes from nine provinces (AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively; E and W in brackets indicate eastern and western provinces, respectively).
Insects 11 00117 g003
Figure 4. A plot of Sitobion avenae biotypes from nine provinces based on the principal coordinate analysis (PCoA) in GenAlEx (AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively).
Figure 4. A plot of Sitobion avenae biotypes from nine provinces based on the principal coordinate analysis (PCoA) in GenAlEx (AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively).
Insects 11 00117 g004
Figure 5. Clustering analyses of clones for two Sitobion avenae biotypes from nine provinces by using STRUCTURE (A: inference of the number of genetic clusters (K); B: the proportion of each cluster for both biotypes at each location; AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively; for more details, see Figure S1).
Figure 5. Clustering analyses of clones for two Sitobion avenae biotypes from nine provinces by using STRUCTURE (A: inference of the number of genetic clusters (K); B: the proportion of each cluster for both biotypes at each location; AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively; for more details, see Figure S1).
Insects 11 00117 g005
Figure 6. A scattered plot of genetic variation FST/(1-FST) vs. geographic distance for Sitobion avenae populations.
Figure 6. A scattered plot of genetic variation FST/(1-FST) vs. geographic distance for Sitobion avenae populations.
Insects 11 00117 g006
Table 1. Collection information for Sitobion avenae samples.
Table 1. Collection information for Sitobion avenae samples.
ProvinceCodeSample SizeHostCoordinatesCollection Date
ZhejiangZJ50wheat, barley120°54′ E; 30°52′ N2016.04
JiangsuJS64wheat, barley120°13′ E; 33°24′ N2016.04
AnhuiAH82wheat, barley116°53′ E; 33°59′ N2016.04
HenanHN83wheat, barley114°02′ E; 33°00′ N2016.04
HubeiHB105wheat, barley112°14′ E; 32°01′ N2016.04
GansuGS75wheat, barley100°46′ E; 38°38′ N2016.06
XinjaingXJ64wheat, barley92°53′ E; 43°36′ N2016.06
QinghaiQH113barley101°44′ E; 36°43′ N2016.07
ShaanxiSX117wheat108°05′ E; 34°17′ N2016.04
Table 2. Genetic diversity indices of two Sitobion avenae biotypes from nine provinces.
Table 2. Genetic diversity indices of two Sitobion avenae biotypes from nine provinces.
GroupsNNaNeHsARIHoHe
Western Provinces
GS1347.8333.4830.7324.0101.4040.4660.645
QH14911.3334.3360.7634.8511.7790.5730.753
SX1107.0004.8410.8015.0911.6390.6170.752
XJ1245.5003.1920.7963.7391.2740.4930.634
Mean29.37.9173.9630.7734.4231.5240.5370.696
GS251.5580.8920.754.0000.2090.0610.057
QH264.1673.3680.7834.0561.3040.5720.698
SX254.0003.1750.7214.0001.1900.7670.653
XJ254.8333.9730.7464.8331.3500.7330.670
Mean5.33.6402.8520.7504.2221.0130.5330.520
Eastern Provinces
AH185.3333.6140.6514.4501.3700.6460.668
HB1218.5004.3010.7614.7331.6550.5870.739
HN1136.1673.1280.6784.1211.3600.7310.654
JS1128.0004.6570.7175.2591.7110.5140.751
ZJ1105.3333.8120.6584.2671.3780.5920.687
Mean12.86.6673.9020.6934.5661.4950.6140.700
Note: AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu, respectively; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang; Na, number of alleles; Ne, number of effective alleles; Hs, gene diversity; AR, allelic richness; I, Shannon’s information index; Ho, observed heterozygosity; He, expected heterozygosity.
Table 3. Pairwise genetic distance estimates (FST) for two Sitobion avenae biotypes from nine provinces.
Table 3. Pairwise genetic distance estimates (FST) for two Sitobion avenae biotypes from nine provinces.
AH1HB1HN1JS1ZJ1GS1QH1SX1XJ1GS2QH2SX2XJ2
AH1
HB10.014
HN10.0010.019
JS10.0030.0130.014
ZJ10.0540.0120.0520.020
GS10.1740.1640.1950.1700.213
QH10.1030.0930.1250.0910.1210.055
SX10.0370.0460.0570.0380.0910.1170.064
XJ10.2110.1910.2440.1970.2410.0290.0950.140
GS20.1520.1510.1880.1290.1900.0310.0510.1060.067
QH20.0710.0610.0900.0590.1120.0760.0010.0120.1230.051
SX20.1160.0920.1170.0820.1300.1750.0880.0780.2210.1860.090
XJ20.1170.0980.1390.1040.1450.0470.0400.0530.0410.0890.0610.120
Note: AH1, biotype 1 of Anhui; HB1, biotype 1 of Hubei; HN1, biotype 1 of Henan; JS1, biotype 1 of Jiangsu; ZJ1, biotype 1 of Zhejiang; GS1 and GS2, biotypes 1 and 2 of Gansu, respectively; QH1 and QH2, biotypes 1 and 2 of Qinghai; SX1 and SX2, biotypes 1 and 2 of Shaanxi; XJ1 and XJ2, biotypes 1 and 2 of Xinjiang, respectively; high levels of genetic divergence (FST > 0.15) highlighted in bold.
Table 4. Analyses of molecular variance (AMOVA) for effects of biotype, geography and plant in the divergence of Sitobion avenae clones.
Table 4. Analyses of molecular variance (AMOVA) for effects of biotype, geography and plant in the divergence of Sitobion avenae clones.
Sources of Variationd. f.Sum of SquaresVariance ComponentsPercentage of Variationp-Value
Biotype effect
Among groups12.510.05 Va0.260.822
Among populations within groups111100.25 Vb11.56<0.001
Within populations391839.322.15 Vc88.18<0.001
Geographic effect
Among groups149.760.24 Va9.70<0.001
Among populations within groups1162.740.13 Vb4.98<0.001
Within populations391839.322.15 Vc85.32<0.001
Plant effect
Among groups14.6160.05 Va0.480.950
Among populations within groups14117.990.28 Vb12.40<0.001
Within populations388829.212.14 Vc87.12<0.001
Note: Clones of S. avenae were divided into two groups in terms of biotype [group 1 (i.e., AH1,GS1, HB1, HN1, JS1, QH1, SX1, XJ1, ZJ1) and group 2 (i.e., GS2, QH2, SX2, XJ2)], two geographic groups [group 1 (i.e., eastern provinces: Anhui, Hubei, Henan, Jiangsu, Zhejiang) and group 2 (i.e., western provinces: Gansu, Qinghai, Shaanxi, Xinjiang)], and two plant-associated groups (group 1, clones on wheat; group 2, clones on barley); significant effects highlighted in bold.
Table 5. Estimates of gene flow between Sitobion avenae biotypes of eastern and western provinces inferred from BayesAss analyses.
Table 5. Estimates of gene flow between Sitobion avenae biotypes of eastern and western provinces inferred from BayesAss analyses.
PopulationE1W1W2
E1 0.0258
(0.0102–0.0463)
0.0823
(0.0339–0.1431)
W10.0053
(0.0001–0.0177)
0.2227
(0.1568–0.2804)
W20.0487
(0.0190–0.0859)
0.0233
(0.0107–0.0407)
Note: 95% confidence intervals presented in parentheses; direction of gene flow is from populations in the left column to those along the top row; E1, biotype 1 of eastern provinces (i.e., Anhui, Hubei, Henan, Jiangsu, and Zhejiang); W1, biotype 1 of western provinces (i.e., Gansu, Qinghai, Shaanxi, and Xinjiang); W2, biotype 2 of western provinces.
Table 6. Estimates of gene flow between Sitobion avenae biotypes collected on barley and wheat inferred from BayesAss analyses.
Table 6. Estimates of gene flow between Sitobion avenae biotypes collected on barley and wheat inferred from BayesAss analyses.
PopulationBA1BA2WH1WH2
BA1 0.0521
(0.0057–0.1265)
0.1205
(0.0833–0.1613)
0.1322
(0.0429–0.2421)
BA20.0246
(0.0058–0.0519)
0.0173
(0.0042–0.0397)
0.0295
(0.0004–0.1066)
WH10.1784
(0.1410–0.2140)
0.2291
(0.1477–0.3003)
0.1317
(0.0426–0.2348)
WH20.0148
(0.0006–0.0372)
0.0128
(0.0001–0.05513)
0.0031
(0.0001–0.0127)
Note: 95% confidence intervals presented in parentheses; direction of gene flow is from populations in the left column to those along the top row; BA1, biotype 1 on barley; BA2, biotype 2 on barley; WH1, biotype 1 on wheat; WH2, biotype 2 on wheat.

Share and Cite

MDPI and ACS Style

Wang, D.; Shi, X.; Liu, D.; Yang, Y.; Shang, Z. Genetic Divergence of Two Sitobion avenae Biotypes on Barley and Wheat in China. Insects 2020, 11, 117. https://doi.org/10.3390/insects11020117

AMA Style

Wang D, Shi X, Liu D, Yang Y, Shang Z. Genetic Divergence of Two Sitobion avenae Biotypes on Barley and Wheat in China. Insects. 2020; 11(2):117. https://doi.org/10.3390/insects11020117

Chicago/Turabian Style

Wang, Da, Xiaoqin Shi, Deguang Liu, Yujing Yang, and Zheming Shang. 2020. "Genetic Divergence of Two Sitobion avenae Biotypes on Barley and Wheat in China" Insects 11, no. 2: 117. https://doi.org/10.3390/insects11020117

APA Style

Wang, D., Shi, X., Liu, D., Yang, Y., & Shang, Z. (2020). Genetic Divergence of Two Sitobion avenae Biotypes on Barley and Wheat in China. Insects, 11(2), 117. https://doi.org/10.3390/insects11020117

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop