Next Article in Journal
Multimodal Information Processing and Associative Learning in the Insect Brain
Previous Article in Journal
FKBP39 Controls the Larval Stage JH Activity and Development in Drosophila melanogaster
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Increasing Frequency of G275E Mutation in the Nicotinic Acetylcholine Receptor α6 Subunit Conferring Spinetoram Resistance in Invading Populations of Western Flower Thrips in China

1
Beijing Key Laboratory for Forest Pest Control, College of Forestry, Beijing Forestry University, Beijing 100083, China
2
Institute of Plant Protection, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
3
School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, VIC 3052, Australia
4
Department of Chemistry and Bioscience, Aalborg University, 9220 Aalborg, Denmark
*
Authors to whom correspondence should be addressed.
Insects 2022, 13(4), 331; https://doi.org/10.3390/insects13040331
Submission received: 27 February 2022 / Revised: 19 March 2022 / Accepted: 25 March 2022 / Published: 28 March 2022
(This article belongs to the Section Insect Molecular Biology and Genomics)

Abstract

:

Simple Summary

The western flower thrips (WFT) Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae) is an important invasive pest in agriculture and forestry. It has developed resistance to a frequently used pesticide spinetoram world widely, including the invading area of China. However, the mechanism of resistance to spinetoram is unclear in China. In this study, we found the presence of the G275E mutation in the nicotinic acetylcholine receptor Foα6 in the early invading populations, which has now increased to a high frequency in China. There was a correlation between the frequency of the G275E mutation and resistance to spinetoram as characterized by median lethal concentration. Our results showed that G275E mutation is one of the mechanisms conferring spinetoram resistance in invading populations in China, as in many other countries. Our study highlights the rapid spread of the G275E mutation in China in the 2009–2021 period.

Abstract

The western flower thrips Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae) is an important invasive pest worldwide. Field-evolved resistance to the pesticide spinetoram is an increasing problem in the chemical control of this pest. Here, we examined changes in the frequency of a genetic mutation associated with spinetoram resistance, the G275E mutation in the acetylcholine receptor Foα6, in 62 field populations collected from 2009 to 2021 across areas of China invaded by this pest. We found a low frequency of the G275E mutation in populations collected at the early invasion stage, in contrast to a high frequency in native USA populations. However, the frequency of the G275E mutation has increased to a high level in recently collected populations, with the mutation becoming fixed in some populations. There was a correlation between the frequency of the G275E mutation and resistance to spinetoram as characterized by median lethal concentration, although two populations were outliers. These results showed that G275E mutation is one of the mechanisms conferring spinetoram resistance in many invading populations in China. Ongoing dispersal of the WFT may have facilitated a rapid increase in the G275E mutation across China. Our study highlights the rapid evolution of pesticide resistance in an invasive species and points to a useful marker for molecular diagnostics of spinetoram resistance.

1. Introduction

The western flower thrips (WFT), Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae), is a polyphagous and widespread invasive pest of agriculture and horticulture [1]. It can severely damage more than 500 host crops by direct feeding and oviposition as well as through transmitting plant viruses in both greenhouse and field situations [2]. The WFT originated from North America and has spread extensively through international trade since the late 1970s [3]. WFT has now successfully invaded many countries on all continents except Antarctica [3]. In China, it was first intercepted in Kunming, Yunnan province in 2000, although it was first reported as an invader in pepper cultivated in the Beijing area in 2003 [4]. Since then, WFT has rapidly colonized all suitable areas of China and become a common pest of vegetable crops and flowers [5].
Following rapid invasion of WFT, chemical control became the most reliable measure to control outbreaks of this species in the absence of natural predators and effective physical control measures. The insecticides used to control WFT since the 1960s include organochlorines, pyrethroids, carbamates, neonicotinoids, avermectins and spinosyns [2,6]. However, as WFT rapidly proliferated and control efforts increased, the inevitable problem of insecticide resistance has emerged. WFT populations from California greenhouses had probably already developed resistance to some classes of insecticides before becoming invasive [7]. Invasive populations may therefore have carried resistance alleles, and WFT in China initially showed a low level of susceptibility to some major insecticides [8]. Currently, some WFT populations have developed resistance to spinosad, spinetoram, abamectin, emamectin-benzoate, thiamethoxam and β-cypermethrin [9]. In particular, there have been a number of cases of rapid and strong resistance developing to spinetoram, with high median lethal concentration (LC50) values and a resistance ratio of >16,000 compared to a susceptible strain [10].
Spinetoram belongs to a group of insecticides known as spinosyns (IRAC class 5) and acts through activation of the nicotinic acetylcholine receptors (nAChRs) in the nervous system of insect pest species [11,12]. Spinosyns are extensively used to control a range of agricultural pests, such as thrips [13,14], Helicoverpa armigera (Hübner) (Lepidoptera: Heliothinae) [15] and Plutella xylostella (Linnaeus) (Lepidoptera: Plutellidae) [16]. Resistance mechanisms are mainly due to a single major gene mutation, which involves the subunit α6 of nAChR in Drosophila melanogaster (Meigen) (Diptera: Drosophilidae) [17]. Glycine at amino acid position 275 of the α6 subunit was mutated to glutamic acid in the resistant strains of both WFT [18,19,20] and the melon thrips, Thrips palmi (Karny) (Thysanoptera: Thripidae) [21]. Moreover, Hiruta et al. [22] found G275 mutated to valine in a spinosad-resistant strain of the flower thrips, Frankliniella intonsa (Trybom) (Thysanoptera: Thripidae). In Australia, a G275 mutation (G275E) has also been associated with field resistance to spinetoram in WFT [19]. However, three susceptible strains and two resistant strains of WFT from China and the United States showed no differences in the cDNA of α6 [23], which indicates that this subunit is not the only mutation involved in spinosad resistance.
Here, we genotyped the G275E mutation of nAChR α6 (Foα6) in sixty field populations from China and two populations from the United States collected during 2009–2021 by Kompetitive Allele-Specific PCR (KASP), which is a quick and efficient molecular diagnostic method [24,25]. We also compared our genotyping results with a conventional bioassay for spinetoram resistance in seven field populations to determine the contribution of the G275E mutation to spinetoram resistance. Our study highlights the rapid spread of the G275E mutation in China in the 2009–2021 period, with potential applications to ongoing monitoring of resistance to spinetoram in WFT populations.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

WFT collections were obtained between 2009 and 2021, coinciding with the period of spread of this pest in China. Adults were collected from 60 locations across their distribution in China with samples from an additional two locations from the US (collected in 2010) also included (Figure 1, Table S1). Chinese samples of WFT were divided into three collection periods coinciding with invasion stages (Table S1). Samples were obtained by collecting adults throughout an area of at least 1000 m2 and then taken to the laboratory in sealed bags. To reduce the likelihood of the collected thrips being close relatives, each individual was collected two meters from another individual. Female adults were used for population genetic analysis because of male-haploid from thrips’ parthenogenetic reproduction. All specimens were kept in absolute ethanol at −80 °C before DNA extraction.
Ten to eighty-seven female adults of WFT from each population were analyzed (Table S1), with a total of 1508 individuals processed. Genomic DNA was extracted from individual WFTs using a high-throughput method; an individual was placed into a well of a 96-well plate containing 80 μL lysis buffer, 1 μL proteinase K [26] and a steel ball (2 mm in diameter). Then, specimens were ground up with a high-throughput tissue grinder (MiniG 1600, SPEX SamplePrep, Metuchen, NJ, USA) for 10 min at 1500 rpm. The grounded specimens were incubated at 65 °C for 30 min followed by 95 °C for 10 min. The supernatant extract containing genomic DNA was kept for subsequent usage after centrifugation for 1 min with a speed of 4000 rpm and stored at −20 °C.

2.2. Bioassay

A Spinetoram (Exalt 6 g·L−1 SC, recommended application rate in 1.5–3.0 μL·m−2, Dow AgroSciences Company, Midland, MI, USA) preparation was used in this study, with the chemical applied through a modified spray tower [27] to test the susceptibility of seven field populations of WFT (Table S1). The recommended concentration of spinetoram for field spray is 10 mg/L. We serially diluted spinetoram into five concentrations ranging from 0.0625–50 mg/L with water containing 0.1% Triton X-100 (Beijing Solar BioScience and Technology Limited Company, China). Cucumber leaves which had never been exposed to pesticides were cut into a circle and pressed into plastic cups (6 cm in diameter and 4 cm in height) with a layer of 1% agar (0.2 cm in height) to seal the leaves and keep them fresh. About 25 female adults were then placed into a cup with an aspirator. The thrips in a cup were anesthetized with carbon dioxide and then cups were transferred to the spray tower. We applied 2 mL of an insecticide solution (or control) to the thrips in the cup. Then, the cups were quickly covered with a tissue and a lid with a circular hole (1.5 cm in diameter). Three replicate cups were tested per concentration and the 0.1% Triton X-100 control treatment. The treated thrips were kept under 25 °C, 60–70% relative humidity and a photoperiod of 16 h:8 h (L:D). Mortality was recorded 48 h after treatment. Thrips were scored as dead when touched with a brush tip and unable to move.

2.3. KASP Primer Design

Two allele-specific forward primers and one common primer were designed by LGC Genomics (UK) based on the DNA segment of the nAChR α6 subunit gene around the G275E mutation [20]. Each forward primer has one of fluorophore (FAM: GAAGGTGACCAAGTTCATGCT, HEX: GAAGGTCGGAGTCAACGGATT) attached as a tail and the sequences are as follows: susceptible allele forward primer, 5′-GAACATGATGCAGTTGAAGTAAGTTC-3′; resistant allele forward primer, 5′-ATGAACATGATGCAGTTGAAGTAAGTTT-3′; and common primer, CGGCCACTGCCTGCGTCTGTTT.

2.4. KASP Assay of G275E Mutation

The DNA samples were assayed in a 384-well plate. The total volume of KASP reaction mixtures was 3 µL, containing 1.5 µL 2 × KASP master mixture (Standard Rox, LGC Genomics, Hoddesdon, Herts., UK), 12 µM of each allele-specific forward primer, 30 µM of common primer, at least 1 ng genomic DNA and 1.5 µL ddH2O.
The KASP reaction was performed on a Hydrocycler-16 machine (LGC Genomics, Hoddesdon, Herts., UK) with the following program: 15 min at 94 °C, followed by 10 touchdown cycles (20 s at 94 °C; 25 s at touchdown at 68 °C initially and decreasing by 0.6 °C per cycle) and then followed by 26 cycles of amplification (20 s at 94 °C; 60 s at 62 °C). Then, the fluorescence signal was read on a Pherastar scanning machine (BMG Latech, Offenburg, Germany). We used SNP viewer 2 (LGC Genomics, UK) to check the KASP assay and output the results.

2.5. Data Analysis

For the bioassay data, the LC50 (median lethal concentration) was computed for each WFT population tested. Parameters for the regression model and LC50 values were estimated in DPS software [28]. The resistance ratio (RR) was calculated by dividing the LC50 of a population by the lowest LC50 obtained from the WFT populations. Differences between populations were considered significant when the 95% confidence intervals (95% CIs) for the LC50 values did not overlap.
Using the KASP results, tests for deviation from the Hardy-Weinberg equilibrium (HWE) at the mutation G275E were performed based on the chi-square statistic with the R package genetics version 1.3.8 after applying a Bonferroni correction (Baltzegar et al., 2021). A one-way analysis of variance (ANOVA) was performed to test for differences in the frequencies of the resistance allele in the three invasion stages using IBM SPSS Statistics version 20 (IBM, Armonk, NY, USA), with pairwise comparisons undertaken based on Tukey’s honestly significant difference (HSD) test. A linear regression was undertaken to investigate whether the frequency of the G275E mutation could be predicted by the (log transformed) LC50 of that population using GraphPad Prism version 9.0 (GraphPad Software, San Diego, CA, USA). We also considered these associations based on the frequency of the resistant homozygote and summed resistant and susceptible homozygotes (to reflect possible different dominance relationships between the resistant allele over the susceptible allele).

3. Results

3.1. Resistance of Seven Field Populations of WFT to Spinetoram

We performed bioassays on seven field populations. The LC50 values of these populations ranged from 0.24 mg·L−1 (BJDJ, 2017) to 52.9 mg·L−1 (BJDJ, 2020) (Table 1, Figure 1). The LC50 value of BJDJ was significantly lower than that of other populations. There were two other populations (BJPG and BJDS) with LC50 values lower than 1.0 mg·L−1. The BJHD population had a higher level of resistance (LC50 = 8.7 mg·L−1, RR = 36.25), while the BJDQ (LC50 = 42.32 mg·L−1, RR = 176.33) and BJDJ (LC50 = 52.90 mg·L−1, RR = 220.42) populations had very high levels of resistance.

3.2. Frequency of G275E Mutations in Spatial and Temporal Populations

We genotyped 1508 individuals for the G275E mutation by KASP. The only genotype detected in the USDZ population was the homozygous resistant genotype (RR). Genotype frequencies in the USA population (USJZ) were 64.7% (RR), 17.6% (RS-heterozygote) and 17.6% (SS-susceptible homozygote). In the 2009–2013 samples from China, the highest (16.7%) frequency of the resistant homozygote (RR) was detected in the XJKL population. No RR or RS individuals were detected in five populations (BJYL, SDQD, BJDX, BJHD2011 and BJHD2013). Genotype frequencies in four other populations (XJKL, BJMT, YNKJ and YNHH) deviated from HWE (Figure 2a). In 2014–2018, the resistance homozygous (RR) frequency was highest (81.8%) in the BJDQ2017 population and also over 50% in the BJDQ2018 population. However, the RR genotype was still relatively rare in the other populations, and eight populations only had the sensitive SS homozygote. Genotype frequencies in three populations (BJCH, BJHD2017 and BJDQ2018) deviated from HWE (Figure 2b). By 2020–2021, the frequency of the RR genotype was over 50% in more populations than in the first period (2009–2013) and second period (2014–2018). The SCCD population only had RR homozygotes, while three populations only had SS homozygote individuals. The genotype frequencies in three populations (BJPG, SCDY and BJYZ2021) deviated from HWE (Figure 2c).
The ANOVA results showed that the frequencies of resistance allele were significantly different across periods (F2,57 = 7.138, p = 0.011, Figure 3). There was a significant difference in the resistance allele frequencies between population groups collected from 2014–2018 and from 2020–2021 (p = 0.012, 95% CI = 4.6045–45.1037, Figure 3). The frequencies of resistance allele increased, being significantly higher in the third period (2020–2021, 37.14%) than in the two earlier periods (8.89% in 2009–2013 and 12.29% in 2014–2018).

3.3. Association between G275E Mutation and Resistance of WFT to Spinetoram

When considering all seven populations, the linear regression analysis showed no strong association between (log transformed) resistance to spinetoram and the frequency of the resistance allele (R2 = 0.3099, p = 0.1943, Figure 4a), the frequency of the RR genotype (R2 = 0.2200, p = 0.2883, Figure 4b) or the summed frequency of the RR and RS genotypes (R2 = 0.3315, p = 0.1761, Figure 4c). When we removed the data of BJDJ2020 and BJHD2021 populations, the linear regression analysis did show an association between resistance and the frequency of the resistance allele (R2 = 0.9543, p = 0.0042, Figure 4d), the frequency of the RR genotype (R2 = 0.8966, p = 0.0146, Figure 4e) and the summed frequency of the RR and RS genotypes (R2 = 0.8799, p = 0.0183, Figure 4f).

4. Discussion

In this study, we considered the detection and distribution of genotypes carrying the nAChR α6 G275E mutation in field populations of WFT during 2009–2021. The results will first be discussed in terms of the distribution of resistance genotypes at different invasion stages and following the local registration of spinetoram against WFT. We then discuss the contribution of the G275E mutation to possible mechanisms underlying spinetoram resistance. Finally, we consider the new high-throughput KASP diagnostic assay developed here as an effective way of detecting the G275E mutation. This diagnostic assay can be used to monitor resistance alleles in field populations in support of integrated pest management (IPM) programs that aim to manage resistance.

4.1. The G275E Resistance Mutation of WFT to Spinetoram in China

Low frequencies of RR and RS genotype of the nAChR α6 G275E mutation were detected in all field WFT populations from China. The G275E mutation was at a low level in the first period (2009–2013) and also in the second period (2014–2018) except in one population in Daxing, Beijing. However, the frequency of the resistant allele markedly increased in the third period by a factor of 4 over the first period. Bioassay results show the same pattern when LC50 values are compared across studies; the susceptibility of populations to spinetoram was high before 2014 with an average LC50 value below 0.07 mg·L−1 [29]. From 2015 to 2018, the susceptibility of WFT to spinetoram decreased in some populations, with an LC50 value for northern populations higher than this value, especially for those from Beijing and the Shandong province [8,9,10]. In our study, the LC50 values of populations in the Beijing area were in the range 0.65–52.9 mg·L−1 in 2020–2021, a further change that points to susceptibility decreasing rapidly after the second period (2014–2018).

4.2. Spinetoram Resistance Mechanisms in WFT

In the association analysis linking the LC50 values of seven populations in the Beijing area to the frequency of the G275E mutation, we only found a weak association unless we removed data from two populations. We suspect that mutation G275E of Foα6 represents an important target-site mutation leading to spinetoram resistance as also noted for the tomato leafminer, Tuta absoluta [30], and consistent with other studies on WFT. The presence of two populations that do not fit this pattern points to another important mechanism being involved in spinetoram resistance in WFT, such as the expression of truncated nAChR α6 variants [31]. In previous studies, metabolism-based resistance mechanisms to spinosyns have been ruled out in laboratory strains of WFT [18,32,33], but we suspect that this requires further study. Multiple mechanisms are known to be involved in other pests such as T. palmi, where spinosyn resistance in the field is conferred by both reduced sensitivity of TPα6 and by cytochrome P450-mediated detoxification [21,34]. However, the detailed genetic basis of spinosyn resistance in these cases has not yet been established, and it is intriguing that the outlier populations have been identified from one region where other populations seem to fit the general pattern of increasing resistance with an increasing frequency of the mutation.

4.3. Association of the G275E Mutation with WFT Invasion and Spread

We detected a high frequency of the resistant mutation in early period collections of American WFT populations, in contrast to a low mutation frequency in the early Chinese WFT populations. With the spread of the WFT in China, the resistance of WFT increased, along with the frequency of the resistant mutation. Ongoing WFT dispersal across the time periods may have accelerated the development of resistance. We have previously noted a population structure in WFT consistent with anthropogenic transmission which can span great distances, introducing new alleles that can then further increase under selection [5]. We have also previously noted the spread of resistance in China in other pests including T. palmi and spider mites [35,36].

4.4. A KASP Diagnostic Assay for the G275E Mutation

For low cost and high throughput genotyping of G275E mutation in WFT, we designed a KASP assay that can genotype the G275E mutation in individual insects from field-collected populations [37]. The diagnostic assay differed in sequence only at the mutation site, with one primer complementary to the susceptible wild-type allele and the other primer complementary to the mutant-resistant allele. This assay allowed us to rapidly genotype 62 field populations from past collections to reveal that 44 out of 62 field populations harbored the G275E mutation. The KASP assay can also be used for high-throughput screening of resistance alleles from small samples. Incorporating our molecular diagnostic assay as an adjunct in support of IPM-based resistance management will provide timely and precise information on resistance frequencies and aid successful management [24,38]. This is particularly important given the high level of variability in resistant allele frequency across populations, although we acknowledge that the likely presence of multiple resistance mechanisms means that a low incidence of the R allele does not preclude resistance being present in a population.

5. Conclusions

We examined spatial and historical variation of G275E mutation in the invasion populations of the WFT across China. Association analysis showed that G275E mutation is one of the mechanisms conferring spinetoram resistance in invading populations in China, although other mutations can also be important. We found that the G275E mutation present in early invading populations of the WFT increased in recent years. Ongoing dispersal of the WFT may have facilitated a rapid increase in the G275E mutation across China. Our study emphasizes the importance of molecular diagnostics of pesticide resistance in the invading populations to manage biological invasions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects13040331/s1. Table S1: Collection information for Frankliniella occidentalis populations.

Author Contributions

S.-J.W. contributed to the study conception and design; L.-N.S., X.-J.S., L.-J.C., J.-C.C. and L.-J.M. conducted the experiments; L.-N.S. and X.-J.S. analyzed the data and wrote the first draft of the manuscript; A.A.H., S.-J.W. and S.-A.W. critically revised the manuscript and interpreted the data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Promotion and Innovation Program of Beijing Academy of Agricultural and Forestry Sciences (KJCX20220409) and Joint Laboratory of Pest Control Research between China and Australia (Z201100008320013).

Institutional Review Board Statement

Not applicable.

Acknowledgments

We thank Ming-Ming Zhang for collecting the recent samples, Yu-Jie Zhang, Pan Shi, Bin-Shuang Pang for her assistance with KASP bioassays.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kirk, W.D.J. The pest and vector from the West: Frankliniella occidentalis. In Proceedings of the Thrips and Tospoviruses: Proceedings of the 7th International Symposium on Thysanoptera, Reggio Calabria, Italy, 2–7 July 2001; pp. 33–42. [Google Scholar]
  2. Reitz, S.R.; Gao, Y.; Kirk, W.D.J.; Hoddle, M.S.; Leiss, K.A.; Funderburk, J.E. Invasion biology, ecology, and management of western flower thrips. Annu. Rev. Entomol. 2020, 65, 17–37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kirk, W.D.J.; Terry, L.I. The spread of the western flower thrips Frankliniella occidentalis (Pergande). Agric. For. Entomol. 2003, 5, 301–310. [Google Scholar] [CrossRef]
  4. Yang, X.M.; Sun, J.T.; Xue, X.F.; Li, J.B.; Hong, X.Y. Invasion genetics of the western flower thrips in China: Evidence for genetic bottleneck, hybridization and bridgehead effect. PLoS ONE 2012, 7, e34567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cao, L.J.; Wang, Z.H.; Gong, Y.J.; Zhu, L.; Hoffmann, A.A.; Wei, S.J. Low genetic diversity but strong population structure reflects multiple introductions of western flower thrips (Thysanoptera: Thripidae) into China followed by human-mediated spread. Evol. Appl. 2017, 10, 391–401. [Google Scholar] [CrossRef] [Green Version]
  6. Gao, Y.L.; Lei, Z.R.; Reitz, S.R. Western flower thrips resistance to insecticides: Detection, mechanisms and management strategies. Pest Manag. Sci. 2012, 68, 1111–1121. [Google Scholar] [CrossRef]
  7. Immaraju, J.A.; Paine, T.D.; Bethke, J.A.; Robb, K.L.; Newman, J.P. Western flower thrips (Thysanoptera, Thripidae) resistance to insecticides in coastal California greenhouses. J. Econ. Entomol. 1992, 85, 9–14. [Google Scholar] [CrossRef]
  8. Wang, Z.H.; Gong, Y.J.; Jin, G.H.; Li, B.Y.; Chen, J.C.; Kang, Z.J.; Zhu, L.; Gao, Y.L.; Reitz, S.; Wei, S.J. Field-evolved resistance to insecticides in the invasive western flower thrips Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae) in China. Pest Manag. Sci. 2016, 72, 1440–1444. [Google Scholar] [CrossRef]
  9. Zhang, K.; Yuan, J.J.; Wang, J.; Hua, D.K.; Zheng, X.B.; Tao, M.; Zhang, Z.; Wan, Y.R.; Wang, S.Y.; Zhang, Y.J.; et al. Susceptibility levels of field populations of Frankliniella occidentalis (Thysanoptera: Thripidae) to seven insecticides in China. Crop Prot. 2022, 153, 105886. [Google Scholar] [CrossRef]
  10. Gao, Y.F.; Gong, Y.J.; Cao, L.J.; Chen, J.C.; Gao, Y.L.; Mirab-balou, M.; Chen, M.; Hoffmann, A.A.; Wei, S.J. Geographical and interspecific variation in susceptibility of three common thrips species to insecticide, spinetoram. J. Pest Sci. 2021, 94, 93–99. [Google Scholar] [CrossRef]
  11. Geng, C.X.; Watson, G.B.; Sparks, T.C. Nicotinic acetylcholine receptors as spinosyn targets for insect pest management. In Advances in Insect Physiology; Cohen, E., Ed.; Academic Press: Oxford, UK, 2013; Volume 44, pp. 101–210. [Google Scholar]
  12. Sparks, T.C.; Nauen, R. IRAC: Mode of action classification and insecticide resistance management. Pestic. Biochem. Physiol. 2014, 121, 122–128. [Google Scholar] [CrossRef] [Green Version]
  13. Fu, B.L.; Li, Q.; Qiu, H.Y.; Tang, L.D.; Zeng, D.Q.; Liu, K.; Gao, Y.L. Resistance development, stability, cross-resistance potential, biological fitness and biochemical mechanisms of spinetoram resistance in Thrips hawaiiensis (Thysanoptera: Thripidae). Pest Manag. Sci. 2018, 74, 1564–1574. [Google Scholar] [CrossRef] [PubMed]
  14. Gholam, Z.; Sadeghi, A. Management strategies for western flower thrips in vegetable greenhouses in Iran: A review. Plant Prot. Sci. 2016, 52, 87–98. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, D.; Gong, P.Y.; Li, M.; Qiu, X.H.; Wang, K.Y. Sublethal effects of spinosad on survival, growth and reproduction of Helicoverpa armigera (Lepidoptera: Noctuidae). Pest Manag. Sci. 2009, 65, 223–227. [Google Scholar] [CrossRef] [PubMed]
  16. Lima Neto, J.E.; Amaral, M.H.P.; Siqueira, H.A.A.; Barros, R.; Silva, P.A.F. Resistance monitoring of Plutella xylostella (L.) (Lepidoptera: Plutellidae) to risk-reduced insecticides and cross resistance to spinetoram. Phytoparasitica 2016, 44, 631–640. [Google Scholar] [CrossRef]
  17. Watson, G.B.; Chouinard, S.W.; Cook, K.R.; Geng, C.X.; Gifford, J.M.; Gustafson, G.D.; Hasler, J.M.; Larrinua, I.M.; Letherer, T.J.; Mitchell, J.C.; et al. A spinosyn-sensitive Drosophila melanogaster nicotinic acetylcholine receptor identified through chemically induced target site resistance, resistance gene identification, and heterologous expression. Insect Biochem. Mol. Biol. 2010, 40, 376–384. [Google Scholar] [CrossRef] [PubMed]
  18. Bielza, P.; Quinto, V.; Contreras, J.; Torne, M.; Martin, A.; Espinosa, P.J. Resistance to spinosad in the western flower thrips, Frankliniella occidentalis (Pergande), in greenhouses of south-eastern Spain. Pest Manag. Sci. 2007, 63, 682–687. [Google Scholar] [CrossRef]
  19. Chen, Y.; Nguyen, D.T.; Gupta, R.; Herron, G.A. Mutation (G275E) of nAChR subunit Foα6 associated with spinetoram resistance in Australian western flower thrips, Frankliniella occidentalis (Pergande). Mol. Biol. Rep. 2021, 48, 3155–3163. [Google Scholar] [CrossRef]
  20. Puinean, A.M.; Lansdell, S.J.; Collins, T.; Bielza, P.; Millar, N.S. A nicotinic acetylcholine receptor transmembrane point mutation (G275E) associated with resistance to spinosad in Frankliniella occidentalis. J. Neurochem. 2013, 124, 590–601. [Google Scholar] [CrossRef] [Green Version]
  21. Bao, W.X.; Narai, Y.; Nakano, A.; Kaneda, T.; Murai, T.; Sonoda, S. Spinosad resistance of melon thrips, Thrips palmi, is conferred by G275E mutation in α6 subunit of nicotinic acetylcholine receptor and cytochrome P450 detoxification. Pestic. Biochem. Physiol. 2014, 112, 51–55. [Google Scholar] [CrossRef]
  22. Hiruta, E.; Aizawa, M.; Nakano, A.; Sonoda, S. Nicotinic acetylcholine receptor α6 subunit mutation (G275V) found in a spinosad-resistant strain of the flower thrips, Frankliniella intonsa (Thysanoptera: Thripidae). J. Pestic. Sci. 2018, 43, 272–276. [Google Scholar] [CrossRef]
  23. Hou, W.J.; Liu, Q.; Tian, L.X.; Wu, Q.J.; Zhang, Y.J.; Xie, W.; Wang, S.L.; Miguel, K.S.; Funderburk, J.; Scott, J.G. The a6 nicotinic acetylcholine receptor subunit of Frankliniella occidentalis is not involved in resistance to spinosad. Pestic. Biochem. Physiol. 2014, 111, 60–67. [Google Scholar] [CrossRef] [PubMed]
  24. Alvarez-Fernandez, A.; Bernal, M.J.; Fradejas, I.; Martin Ramirez, A.; Md Yusuf, N.A.; Lanza, M.; Hisam, S.; Perez de Ayala, A.; Rubio, J.M. KASP: A genotyping method to rapid identification of resistance in Plasmodium falciparum. Malar. J. 2021, 20, 16. [Google Scholar] [CrossRef] [PubMed]
  25. Semagn, K.; Babu, R.; Hearne, S.; Olsen, M. Single nucleotide polymorphism genotyping using kompetitive allele specific PCR (KASP): Overview of the technology and its application in crop improvement. Mol. Breed. 2014, 33, 1–14. [Google Scholar] [CrossRef]
  26. Mardulyn, P.; Goffredo, M.; Conte, A.; Hendrickx, G.; Meiswinkel, R.; Balenghien, T.; Sghaier, S.; Lohr, Y.; Gilbert, M. Climate change and the spread of vector-borne diseases: Using approximate Bayesian computation to compare invasion scenarios for the bluetongue virus vector Culicoides imicola in Italy. Mol. Ecol. 2013, 22, 2456–2466. [Google Scholar] [CrossRef]
  27. Martin, N.A.; Workman, P.J.; Butler, R.C. Insecticide bioassays for western flower thrips (Frankliniella occidentalis) (Thysanoptera: Thripidae) and greenhouse whitefly (Trialeurodes vaporariorum) (Hemiptera: Aleyrodidae). N. Z. J. Crop Hortic. Sci. 2005, 33, 177–184. [Google Scholar] [CrossRef]
  28. Tang, Q.Y.; Zhang, C.X. Data processing system (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research. Insect Sci. 2013, 20, 254–260. [Google Scholar] [CrossRef]
  29. Wang, Z.H.; Hou, W.J.; Hao, C.Y.; Wu, Q.J.; Xu, B.Y.; Zhang, Y.J. Monitoring the insecticide resistance of field populations of western flower thrips, Frankliniella occidentalis, in the Beijing area. Chin. J. Appl. Entomol. 2011, 48, 542–547. [Google Scholar]
  30. Silva, W.M.; Berger, M.; Bass, C.; Williamson, M.; Moura, D.M.N.; Ribeiro, L.M.S.; Siqueira, H.A.A. Mutation (G275E) of the nicotinic acetylcholine receptor α6 subunit is associated with high levels of resistance to spinosyns in Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae). Pestic. Biochem. Physiol. 2016, 131, 1–8. [Google Scholar] [CrossRef] [Green Version]
  31. Wan, Y.; Yuan, G.; He, B.; Xu, B.; Xie, W.; Wang, S.; Zhang, Y.; Wu, Q.; Zhou, X. Foccalpha6, a truncated nAChR subunit, positively correlates with spinosad resistance in the western flower thrips, Frankliniella occidentalis (Pergande). Insect Biochem. Mol. Biol. 2018, 99, 1–10. [Google Scholar] [CrossRef]
  32. Bielza, P.; Quinto, V.; Fernandez, E.; Gravalos, C.; Contreras, J. Genetics of spinosad resistance in Frankliniella occidentalis (Thysanoptera: Thripidae). J. Econ. Entomol. 2007, 100, 916–920. [Google Scholar] [CrossRef]
  33. Zhang, S.Y.; Kono, S.; Murai, T.; Miyata, T. Mechanisms of resistance to spinosad in the western flower thrip, Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae). Insect Sci. 2008, 15, 125–132. [Google Scholar] [CrossRef]
  34. Shi, P.; Guo, S.K.; Gao, Y.F.; Chen, J.C.; Gong, Y.J.; Tang, M.Q.; Cao, L.J.; Li, H.; Hoffmann, A.A.; Wei, S.J. Association between susceptibility of Thrips palmi to spinetoram and frequency of G275E mutation provides basis for molecular quantification of field-evolved resistance. J. Econ. Entomol. 2021, 114, 339–347. [Google Scholar] [CrossRef] [PubMed]
  35. Shi, P.; Cao, L.J.; Gong, Y.J.; Ma, L.; Song, W.; Chen, J.C.; Hoffmann, A.A.; Wei, S.J. Independently evolved and gene flow-accelerated pesticide resistance in two-spotted spider mites. Ecol. Evol. 2019, 9, 2206–2219. [Google Scholar] [CrossRef] [PubMed]
  36. Shi, P.; Guo, S.K.; Gao, Y.F.; Cao, L.J.; Gong, Y.J.; Chen, J.C.; Yue, L.; Li, H.; Hoffmann, A.A.; Wei, S.J. Variable resistance to spinetoram in populations of Thrips palmi across a small area unconnected to genetic similarity. Evol. Appl. 2020, 13, 2234–2245. [Google Scholar] [CrossRef] [PubMed]
  37. Yang, G.; Chen, S.P.; Chen, L.K.; Sun, K.; Huang, C.H.; Zhou, D.H.; Huang, Y.T.; Wang, J.F.; Liu, Y.Z.; Wang, H.; et al. Development of a core SNP arrays based on the KASP method for molecular breeding of rice. Rice 2019, 12, 21. [Google Scholar] [CrossRef]
  38. Xu, X.Y.; Li, G.Q.; Bai, G.; Bernardo, A.; Carver, B.F.; St Amand, P.; Bian, R. Characterization of an incomplete leaf rust resistance gene on chromosome 1RS and development of KASP markers for Lr47 in wheat. Phytopathology 2021, 111, 649–658. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Sample locations (red triangle) of Frankliniella occidentalis populations collected in (a) 2009–2013, (b) 2014–2018 and (c) 2020–2021. The larger red triangle represents multiple populations collected from Beijing.
Figure 1. Sample locations (red triangle) of Frankliniella occidentalis populations collected in (a) 2009–2013, (b) 2014–2018 and (c) 2020–2021. The larger red triangle represents multiple populations collected from Beijing.
Insects 13 00331 g001
Figure 2. Genotype frequencies of G275E mutation in Frankliniella occidentalis populations during (a) 2009–2013, (b) 2014–2018 and (c) 2020–2021. RR, resistant homozygote genotype; RS, heterozygote genotype; SS, sensitive homozygote genotype. The “*” represents deviations from Hardy-Weinberg equilibrium (p < 0.05); NA represents a population consisting only of homozygous genotypes.
Figure 2. Genotype frequencies of G275E mutation in Frankliniella occidentalis populations during (a) 2009–2013, (b) 2014–2018 and (c) 2020–2021. RR, resistant homozygote genotype; RS, heterozygote genotype; SS, sensitive homozygote genotype. The “*” represents deviations from Hardy-Weinberg equilibrium (p < 0.05); NA represents a population consisting only of homozygous genotypes.
Insects 13 00331 g002
Figure 3. Plot of G275E-resistant allele frequencies in three collecting periods. The horizontal black solid line represents the median for each period; the red dot represents the mean. Red error bars indicate the standard error.
Figure 3. Plot of G275E-resistant allele frequencies in three collecting periods. The horizontal black solid line represents the median for each period; the red dot represents the mean. Red error bars indicate the standard error.
Insects 13 00331 g003
Figure 4. Regression between (ln) LC50 and the frequencies of (a) the resistant allele, (b) RR resistant genotype and (c) the sum of homozygous and heterozygote genotypes with the resistant G275E mutation across all populations tested. We have also provided the regressions when the BJDJ2020 and BJHD2021 populations are excluded (df).
Figure 4. Regression between (ln) LC50 and the frequencies of (a) the resistant allele, (b) RR resistant genotype and (c) the sum of homozygous and heterozygote genotypes with the resistant G275E mutation across all populations tested. We have also provided the regressions when the BJDJ2020 and BJHD2021 populations are excluded (df).
Insects 13 00331 g004
Table 1. Regressions for spinetoram responses and frequencies of G275E R/S genotypes in seven Frankliniella occidentalis populations.
Table 1. Regressions for spinetoram responses and frequencies of G275E R/S genotypes in seven Frankliniella occidentalis populations.
PopulationYearToxicity Regression
Equation (y = bx + a)
LC50 (95% CI) (mg·L−1) *SE (of b)χ2p-ValueRRR%RR%RS%
BJDJ20172.76x + 6.730.24 (0.18–0.30) D0.43990.57550.902012.204.5
BJDQ20202.07x + 4.2842.32 (35.82–50.19) A0.19410.47370.9246176.3347.239.116.3
BJPG20202.45x + 5.230.81 (0.62–1.04) BC0.38471.40030.49653.3817.210.313.8
BJDJ20202.45x + 0.7852.90 (43.87–64.28) A0.29549.46550.0237220.4214.3028.6
BJTZ20203.35x + 4.901.07 (0.92–1.30) B0.36270.41640.81204.4615.9031.8
BJDS20212.13x + 5.400.65 (0.47–0.84) C0.22509.44010.05102.714.609.1
BJHD20210.91x + 4.158.70 (3.91–54.323) A0.107819.17110.001836.251.603.3
SE, standard error; RR, resistance ratio. * Differences in LC50 were considered significant when the 95% CIs did not overlap, and these are marked by different letters.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, L.-N.; Shen, X.-J.; Cao, L.-J.; Chen, J.-C.; Ma, L.-J.; Wu, S.-A.; Hoffmann, A.A.; Wei, S.-J. Increasing Frequency of G275E Mutation in the Nicotinic Acetylcholine Receptor α6 Subunit Conferring Spinetoram Resistance in Invading Populations of Western Flower Thrips in China. Insects 2022, 13, 331. https://doi.org/10.3390/insects13040331

AMA Style

Sun L-N, Shen X-J, Cao L-J, Chen J-C, Ma L-J, Wu S-A, Hoffmann AA, Wei S-J. Increasing Frequency of G275E Mutation in the Nicotinic Acetylcholine Receptor α6 Subunit Conferring Spinetoram Resistance in Invading Populations of Western Flower Thrips in China. Insects. 2022; 13(4):331. https://doi.org/10.3390/insects13040331

Chicago/Turabian Style

Sun, Li-Na, Xiu-Jing Shen, Li-Jun Cao, Jin-Cui Chen, Li-Jun Ma, San-An Wu, Ary Anthony Hoffmann, and Shu-Jun Wei. 2022. "Increasing Frequency of G275E Mutation in the Nicotinic Acetylcholine Receptor α6 Subunit Conferring Spinetoram Resistance in Invading Populations of Western Flower Thrips in China" Insects 13, no. 4: 331. https://doi.org/10.3390/insects13040331

APA Style

Sun, L. -N., Shen, X. -J., Cao, L. -J., Chen, J. -C., Ma, L. -J., Wu, S. -A., Hoffmann, A. A., & Wei, S. -J. (2022). Increasing Frequency of G275E Mutation in the Nicotinic Acetylcholine Receptor α6 Subunit Conferring Spinetoram Resistance in Invading Populations of Western Flower Thrips in China. Insects, 13(4), 331. https://doi.org/10.3390/insects13040331

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