1. Introduction
The green peach aphid (GPA) Myzus persicae (Sulzer, 1776) is an enormously polyphagous pest worldwide that has been reported to breed on more than 400 kinds of plants belong to 50 families, including crop, vegetable, fruit, and tobacco. Infestations of GPA reduced plant growth rates and caused reductions in crop quantity and quality [
1]. In addition, the aphid vectors over 100 plant viruses, accounting for 67.7% of the aphid–vector viruses, which can lead to large yield reduction [
2]. As a highly variable species, aphids evolve different strains and biotypes, distinguished by biology, host-plant preference, and insecticide resistance [
3], which favor spread and rapid colonization. The cabbage Brassica oleracea var. capitata Linnaeus is the primary vegetable, having the largest vegetable planting area in China of 0.4-million hectares, which accounts for 25%–30% of Chinese vegetables. GPA is a destructive pest of B. oleracea throughout the country. Effective control is crucial to protect B. oleracea from aphid infection and prevent viral transmission.
Although some other Integrated Pest Management measures, such as natural enemies, cultural control, and mechanical control, have been used to control the GPA, chemical control is still the main way to suppress the pest in practice. Based on IRAC, there are 16 groups of effective insecticides that kill the GPA by acting on nerve and muscle (10 groups), growth and development (four groups) and respiration targets (two groups) targets [
4]. In China, the systemic use of insecticides has played a vital role in controlling this pest since the introduction of organophosphates and carbamates insecticides in the 1950s [
5]. To date, more than 600 pesticide products have been registered to manage the cabbage GPA in China. Among these products, 336 products, accounting for 55.17%, contain neonicotinoids, and 158 products, accounting for 25.94%, contain pyrethroids [
6], which indicates that the two kinds of insecticides are sill the major insecticides used to eliminate the GPA of cabbage, and they have been unitized for several decades. Overreliance on these two kinds of insecticides has led to the occurrence and expansion of insecticide resistance in the GPA and threatens their continued efficacy. Insecticide-resistance monitoring data in China indicates that GPA field populations have high resistance to organophosphate insecticides, such as dimethoate and malathion, carbamate insecticides, such as pirimicarb and methomyl, pyrethroids, such as fenvalerate and cypermethrin, and neonicotinoids, such as imidacloprid [
5,
7,
8].
The evolution and spread of resistance in GPA represent major threats to its sustainable control [
9]. Understanding resistance mechanisms is critical for the development of rational strategies to prolong the lives of current insecticides. To date, seven independent mechanisms of resistance have been found in this species, including alterations in target site, metabolic mechanisms, and physical resistance mechanisms [
9]. The enhanced expression of insecticide detoxifying enzymes is the most common mechanism in pests and often presents the greatest challenge. The overproduction of carboxylesterases (E4 and EF4) causes resistance of M. persicae to organophosphates and carbamates [
10,
11]. The amplification of the cytochrome P450 gene CYP6CY3 is associated with low level resistance to neonicotinoids [
12,
13,
14], whereas the upregulation of glutathione S transferase (GST) and UDP-glycosyltransferases (UGTs) have been shown to affect insecticide resistance [
15,
16]. In addition, the protective antioxidant enzymes, superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), which play important roles in immunity, help maintain a state of equilibrium in insects after insecticide exposure [
17,
18,
19,
20]. Field populations are exposed to different insecticides under variable selection pressures, which favors the evolution of multiple resistance through resistance mechanisms. Understanding these mechanisms will help us estimate the prevalence of resistant phenotypes when bioassays cannot be used to detect resistance [
21,
22].
Brassica species are cultivated year-round in Fujian Province, southeast China and are vulnerable to GPA and other notorious pests, such as Plutella xylostella, Phyllotreta striolata (Fabricius), Spodoptera litura Fabricius, and Pieris rapae. Various insecticides are frequently applied by farmers to control the damage from these pests. Presently, some active ingredients, such as flonicamid, cyantraniliprole, spirotetramat, afidopyropen, sulfoxaflor, and flupyradifurone, having excellent aphicidal activities, have been registered for GPA management in China [
23]. Furthermore, other insecticides, such as broflanilide, have also entered the Chinese market and showed high potential in aphid control, although they have not been registered for aphid management [
24]. These compounds belong to different groups of insecticides and may be alternated with pyrethroids or neonicotinoids in the field. Resistance to some insecticides, such as spirotetramat [
25], flupyradifurone [
26], and sulfoxaflor [
16], has been detected in other countries, but the resistance status of GPA is unavailable in China, including Fujian Province. To warn, delay, or prevent the development of resistance in GPA to these compounds and better guide management practices, resistance monitoring in a GPA field population collected from Fujian Province was conducted. To further investigate the resistance mechanisms, biochemical assays and molecular methods were used to determine the activities and expression levels of detoxifying and protective enzymes, as well as related candidate genes, that might help monitor insecticide resistance.
2. Materials and methods
2.1. Insects
The laboratory susceptible strain of the GPA (FFJ-S) was continuously reared without insecticides exposure in the Fujian Institute of Plant Protection, Fujian Academy of Agricultural Science, Fuzhou, China, since it was collected from Arabidopsis thaliana in Jianxin District of Fuzhou, China in 2008. The resistant field population (SEF-R) was collected from cabbage in the Jin’an District of Fuzhou (26°12′84′′N, 119°33′59′′E) in 2022, where heavy insecticide spraying was carried out to control the aphid and other pests. Both GPA strains were reared on seedlings of Brassica campestris ssp. Chinensis (L.) at 21 ± 2°C and 65% ± 5% relative humidity, under a photoperiod of 16:8 h (light: dark). The SEF-R was reared one generation before being subject to a bioassay.
2.2. Pesticides, synergists, and other chemicals
The insecticides used for bioassays included flupyradifurone (96% purity; Bayer AG, Germany), thiacloprid (97.5% purity; Bayer AG), imidacloprid (97% purity; Bayer AG), esfenvalerate (90% purity; Shandong Huayang Technology Co., Ltd., China), beta-cypermethrin (95% purity; Nanjing Red Sun Co., LTD, China), flonicamid (96% purity; Ishihara Sangyo Kaisha, Ltd., Japan), cyantraniliprole (FMC Corporation, USA), spirotetramat (96% purity; Bayer AG), afidopyropen (92.5% purity; BASF SE, Germany), broflanilide (98% purity; Mitsui Chemical lnc., Japan), thiamethoxam (98% purity; Syngenta Group, Switzerland), Pymetrozine (95% purity; Syngenta Group, Switzerland), and sulfoxaflor (95.9% purity; Corteva Agriscience, USA). Sulfoxaflor
Synergists piperonyl butoxide (PBO; reagent grade), S,S,S-tributyl phosphorotrithioate (DEF; reagent grade), and diethyl maleate (DEM; reagent grade) were purchased from Sinopharm Group (Holding) Co. Ltd, China. Ethylenediaminetetraacetic acid, albumin bovine (BSA), and sodium dodecyl sulfate were purchased from Shanghai Aladdin Bio-Chem Technology Co., LTD (Shanghai, China). Eserine, α-naphthyl acetate, fast blue B salt, p-nitroanisole, n-phenylthiourea, Coomassie Brilliant Blue G250, DL-dithiothreitol, phenylmethylsulfonyl fluoride, Triton X-100, and other chemicals were purchased from Sigma-Aldrich, USA.
2.3. Bioassay methods
To assess the susceptibility of GPA to selected insecticides, the leaf-dip method recommended by IRAC [
4] was used under laboratory conditions. All the tested insecticides were dissolved in 0.1% Triton X-100 to prepare five to seven concentrations. Clean leaf-discs of B. oleracea L. were dipped in insecticide solutions for 10 s. After being air dried at room, the leaf-discs were placed onto 1% agar plate (20-mm depth) in a petri dish (30-mm diameter and 40-mm depth). Then 20 2-d-old apterous aphids (for spirotetramat and pymetrozine test) or apterous adults (for other insecticides) were each transferred onto a treated leaf-disc using a paint brush, and each unit was covered with a close-fitting, ventilated lid. Mortalities were observed 4 d later to determine the impact of spirotetramat and pymetrozine and 3 d later to determine other insecticide effects. Aphids that could not right themselves within 10 s once turned on their backs were considered as dead.
2.4. Resistance recession
To test the stability of resistance, the SEF-R strain was reared for 15 generations (equivalent of approximately 1 year in the field) without insecticide exposure to investigate whether resistance to imidacloprid, beta-cypermethrin, sulfoxaflor, and pymetrozine was stable. GPA of the 1st (G1), 4th (G4), 7th (G7), 10th (G10), 12th (G12), and 15th (G15) generations were subjected to a bioassay. The average rate of change in response per generation (R) was estimated in accordance with Tabashnik [
27], as follows:
where n represents the number of generations. An increase or decrease in resistance is indicated by positive and negative R values, respectively.
2.5. Synergism experiment
The effects of three synergists (PBO, DEF, and DEM) each in combination with a thiacloprid or cypermethrin mixture against GPD were evaluated using the bioassay method. The highest doses of PBO, DEF, and DEM on the susceptible strain that led to zero mortality were 0.08, 0.05, and 0.06 g L−1, respectively, using the bioassay method. Apterous adult aphids were exposed to leaf discs that were treated with PBO, DEF, or DEM combined with the thiacloprid and cypermethrin mixture. Four replicates were performed for each concentration in the bioassay.
2.6. Detoxification enzyme activity assays
The esterase activity was measured using a microplate reader in accordance with Byrne and Devonshire [
28]. Individual aphids were homogenized in 60 μL of ice-cold phosphate buffered saline (PBS) (0.02 M, pH 6.5) and centrifuged at 4°C and 10,000 × g for 10 min. The supernatant was used as the esterase source. After incubating for 5 min at room, 100 μL substrate α-naphthyl acetate (107 mM) and 100 μL of eserine (107 mM) were mixed with the color development reagent (6 mg fast blue RR salt dissolved in 10 mL of 0.02 M PBS, pH 6.5) in the dark. Then, 150 μL of this reaction liquid was added to a 96-well microplate (Corning Life Sciences, USA), followed by the addition of 50 μL of supernatant fluid and 50 μL of PBS per well. The absorbance at 450 nm was measured for 5 min continuously using the kinetic model in a microplate reader (SPECTRA max PLUS384, Molecular Devices, USA). The esterase activity was reported as mOD450/min/aphid. For each population, approximately 100 individuals were measured.
The GST activity was determined in accordance with a published method involving 1-chloro-2,4-dinitrobenzene, with slight modifications [
29]. Briefly, 60 adult aphids were homogenized in 2.0 mL ice-cold PBS (0.04 mol L−1, pH 7.5), and the supernatant solution was used after being centrifuged at 10,000 × g for 10 min at 4°C. Briefly, a 300-μL reaction mixture containing 100 μL diluted enzyme solution, 100 μL 1-chloro-2,4-dinitrobenzene (1.2 mM) substrate solution, and 100 μL glutathione ((6 mM) was prepared, after which the absorbance was measured at 340 nm for 10 min using the kinetic model. The results were determined based on the protein concentration of an enzyme source, and the specific activity was converted from an OD value.
The monooxygenase enzyme (MFO) activity was measured in accordance to Shang’s method [
30]. Briefly, 60 adult aphids were homogenized in 2.0 mL ice-cold PBS (0.04 mol L−1, pH 7.8). The supernatant obtained by centrifuging 10,000 × g for 10 min at 4°C was added to a reaction unit containing NADPH and nitroanisole (0.05 mol L−1 in acetone) as a substrate. Hydrochloric acid (1 mol L−1) was added to terminate the reaction after incubation for 30 min at 37°C. Then, the reaction unit was extracted by a sodium hydroxide (NaOH) and chloroform solution. Finally, the optical density (OD) of the enzyme source was recorded at 400 nm using a microplate reader. The specific activity was obtained using a nitrophenol standard curve and the protein concentration of the enzyme source.
2.7. Determination of protective enzyme activities in M. persicae
The activities of the protective enzymes CAT (catalog no. A007-2), POD (catalog no. A084-1), and SOD (catalog no. A001-3) were determined using commercial assay kits in accordance with the manufacturer’s instructions (Nanjing Jiancheng Bioengineering Institute, China). A spectrophotometric method based on the ultraviolet absorption of peroxide released from the activity of CAT on hydrogen peroxide was used to determine the CAT level, and the nitroblue tetrazolium reduction method and the guaiacol method were used to determine the SOD and POD activities, respectively.
The protein contents of the enzyme solutions were determined using the Bradford method [
31]. Serial dilutions of a BSA solution, and the samples were measured together. The protein contents of the samples were calculated using a standard curve based on the BSA solution. The diluted enzyme solutions (50 μL) were mixed with Coomassie Brilliant Blue (200 μL). After incubating at 25°C for 10 min, the absorbance at 595 nm was measured. All the experiments were repeated three times, and the average values were obtained from the triple-replicated data.
2.8. Target-site resistance (Detection of R81T, L1014F and M918T mutations)
Genetic analyses of GPA from SEF-R were conducted for known targe site resistance. This included the nicotinic acetylcholine receptor R81T mutation and the voltage-gated sodium channel knockdown resistance mutations L1014F (kdr) and M918L (Super-kdr), responsible for neonicotinoids and pyrethroids resistance in GPA, respectively. The real-time TaqMan assays were used for the mutation identification of genome DNA iin a MasterCycler ep RealPlex4 machine (Eppendorf, Hamburg, GER) according to the methods described in other studies [
32,
33].
2.9. Gene expression analysis
The mRNA levels of different enzyme-encoding genes in the FFJ-S and SEF-R strains were measured using RT-qPCR with SYBR® Green Supermix (ThermoFisher, USA) in a qTOWER 2.2 real-time quantitative PCR system (Analytikjena, Germany). Total RNA was extracted as described above and quantified using a ScanDrop 100 spectrophotometer (Analytikjena) in accordance with the manufacturer’s protocols. The RNA concentration was diluted to 0.8 μg μL–1 with diethyl pyrocarbonate-treated H2O, and 0.8 μg of RNA was then reverse transcribed in a 20-μL reaction volume using a TUREscript 1st Stand cDNA Synthesis Kit (Aidlab, Beijing, China) with the actin gene as an internal control (NCBI gene ID: 836110). In total, 15 genes, CYP6CY3 (Gene ID: KF218356), CYP6K1 (Gene ID: XM_022323152), CYP6CY23 (Gene ID: MF155662), CYP380C40 (Gene ID: OM677847), CYP6CY9 (Gene ID: XM_022316140), CYP6CY56 (Gene ID: MT472683), CYP4G15 (Gene ID: XM_022323310), CYP3CY06 (Gene ID: XM_022312311), E4 (Gene ID: X74554), THEM (Gene ID: XM_022326659), GST (Gene ID: 111036096), UGT344P2 (Gene ID: OM677846), POD (Gene ID: 112683716), SOD (Gene ID: 111035379), and CAT (Gene ID: 111041019), were quantified. The primers for these genes were designed using the Primer Express 3.0 software based on the sequences of target genes in NCBI or according other research, and they are listed in Supplementary
Table 1.
Each RT-qPCR was conducted in a 20-μL mixture containing 1 μL of sample cDNA, 1 μL of each primer (200 nM), 6 μL of diethyl pyrocarbonate-treated H2O, and 10 μL of 2×SYBR® Green Supermix. The qPCR cycling parameters were as follows: 95°C for 3 min, followed by 39 cycles of 95°C for 10 s and 58°C for 30 s. A plate reader was used for the analysis. Melting curve generation was performed from 60 to 95°C (+1℃ (+1ti, holding time 4 s). After adding all the components, samples were centrifuged at 6,000 ×g for 1 min to keep all the components in the tube bottom. To check the reproducibility of the assay results, the qPCR for each sample was performed using three technical replicates and three biological replicates. The relative gene expression was calculated automatically using qPCRsoft3.2 software. The comparative 2−ΔΔCT method [
34] was used for the relative quantification calculation.
2.10. In vivo RNAi of CYP6CY3 and bioassays
Fragments of CYP6CY3 were cloned into pGEM-T (Promega, USA) and used as templates for dsRNA-CYP6CY3 synthesis using a T7 RiboMAX™ Express RNAi System (Promega) in accordance with Peng et al. [
35]. The dsRNA-CYP6CY3 (100 ng mL–1 in diet) was added to the artificial diet [
36], and DEPC water and dsGFP (dsRNA of GFP) were also added to the artificial diet as the control. The dsRNA-CYP6CY3-treated artificial diet containing imidacloprid (100 mg L–1), thiacloprid (40 mg L–1), thiamethoxam (30 mg L–1), sulfoxaflor (50 mg L–1), or flupyradifurone (10 mg L–1) was used in the toxicity assessment, and the diet supplemented with DEPC water or dsGFP, but without insecticides, was used as the control. The experiments were performed in triplicate. To determine the efficiency of the dsRNA knockdowns of CYP6CY3, the aphids were collected after they fed on the artificial diet for 48 h, and the samples were used for qRT-PCR analysis. The mortality of the aphids was also calculated for 48 h. Each treatment was replicated in triplicate.
2.11. Data analysis
The enzymatic activities and the gene expression data between SEF-R and FJJ-S were compared with Student’s t-tests, and the mortalities of the GPA that were RNAi plus insecticide-treated versus mortalities of the control were compared with a one-way analysis of variance in combination with a Fisher’s least significant difference multiple comparison tests. These experiments were analyzed using DPS Data Processing System (Hangzhou Ruifeng Technology LTD, China) [
37]. All the data are presented as the means ± standard errors (SEs). Differences were considered significant at P-value < 0.05 or 0.01. Median lethal concentrations (LC50) of each insecticide for adult apterous aphids were obtained by DPS Data Processing System. The LC50 values without overlap among 95% confidence intervals (CIs) were considered significantly different. The Resistance Factor (RF) was calculated as the ratio of the LC50 value of the field population to the LC50 value of the FFJ-S strain. Resistance levels were classified in accordance with the standard reported by World Health Organization [
38]: susceptibility (RF < 3-fold), decreased susceptibility (RF = 3 to 5-fold), low resistance (RF = 5 to 10-fold), moderate resistance (RF = 10 to 40-fold), high resistance (RF = 40 to 160-fold), and extremely high resistance (RR > 160-fold).
4. Discussion
The evolution and spread of insecticide resistance in the GPA are considered major challenges in the management and elimination of the pest and its transmitted virus. To date, the GPA has developed resistance to almost all insecticides used worldwide. An aphid monitoring survey suggested that the GPA in a cabbage field located in southeastern China, Fujian Province, was not susceptible to the 13 tested insecticides categorized to six MoA groups [
4]. Multi-insecticide resistance might drastically reduce the number of insecticides that can be used to control the insect in this district. To delay resistance development and prolong the product’s effectiveness, implementing insecticide resistance management strategies, including rational applications of insecticides, is urgently needed according to monitoring data.
Similar to most other countries, neonicotinoids (Group 4A) are one of the most important chemical classes of insecticides used against the GPA in China [
39]. More than 55% of all registered products in China contain neonicotinoids, which have been a key instrument in GPA management [
6]. However, the GPA have the second-greatest number of cases of neonicotinoid resistance, behind only those of Bemisia tabaci, according to the Arthropod Pesticide Resistance Database (APRD), and this represents a serious threat to the effective control of GPA by neonicotinoids [
40]. The high level of imidacloprid-resistance detected in SEF-R reflected its status in cabbage fields, which is consistent with other monitoring data. Since susceptibility in the GPA was shown to decrease in 2005 [
41], elevated levels of resistance to imidacloprid have been observed in most of the provinces of China, such as Beijing [
42], Hunan [
43], Jiangsu, and Shanghai [
44], in which GPA have been tested. According to published documents, confirmed resistance to imidacloprid is much higher in vegetable aphid populations (RF = 32–350.9-fold) than in tobacco aphid populations (RF = 6.03–22-fold) in China [
8]. Additionally, GPA collected from cabbage field in Longyao, Heibei Province exhibited the highest resistance (2,853-fold) [
44]. We speculated that selection pressures rather than host plants shaped the evolution of imidacloprid resistance because GPA collected from tobacco fields in Europe evolved moderate to high resistance to imidacloprid (RF >15 detected in 86.7% populations in 2015) [
45]. Less insecticide needed to be used against GPA because of the short growing season of tobacco in China, which also indicated that the limited application of imidacloprid delayed resistance development. The present study revealed moderate level resistance in the SEF-R clone to the other two neonicotinoids, thiacloprid and thiamethoxam, which were registered in 2015 to kill aphids on cabbage. Zhang [
46] found that GPA collected from cabbage fields in Guangdong and Hunan Provinces, in southern and southeastern China, respectively, are still susceptible to thiamethoxam. There are few products of thiamethoxam (eight products) and thiacloprid (four products) allowed in cabbage aphid management, which may reduce the number of their spray applications in fields [
6]. We inferred that the moderate resistance to the two nicotinoids in SEF-R was caused by cross-resistance between thiacloprid or thiamethoxam with imidacloprid verified in other studies [
47,
48].
As nAChR competitive modulators, sulfoxaflor (Subgroup 4C) and flupyradifurone (Subgroup 4D) belong to the novel classes sulfoximines and butenolides, respectively. Previous studies showed that the two compounds are effective against some GPA populations exhibiting neonicotinoid resistance [
49,
50], but later researches confirmed the cross-resistance between sulfoxaflor or flupyradifurone with imidacloprid in some field populations [
48,
51]. SEF-R exhibited resistance to the two compounds. Today, sulfoxaflor (Tefuli®) and flupyradifurone (Jinxian®) are only registered by Corteva Agriscience and Bayer AG, respectively, which make us easily know their application history. Growers did not apply the two compounds in the field in the recent years. However, field populations of GPA in other countries have developed high resistance to the two compounds. In total, 25% of the field-collected GPA in Greece had an RF˃119.2 for flupyradifurone [
26], and field populations collected in Western Australia and Spain exhibited 4–23-fold and 0.6-61-fold resistance to sulfoxaflor, respectively [
16] (Mezei et al. 2022). Our study first reported the resistance of GPA to sulfoxaflor and flupyradifurone in China, which might cause by high resistance to imidacloprid. The loss of sensitivity of GPA to sulfoxaflor and flupyradifurone threatened the applications of the two insecticides as alternative tools for integrated pest and resistance management. The resistance of SEF-R to the five Group 4 insecticides makes it difficult select related products. Stability experiments showed that the RFs of imidacloprid and sulfoxaflor dropped to 18.7- and 9-fold, respectively, after 15 generations without selection (
Figure 1A). The study highlighted that the cessation of selective pressure could help restore the susceptibility of GPA to imidacloprid and sulfoxaflor, but they should be prohibited for at least one year. For Group 4 compounds, thiamethoxam and flupyradifurone should be recommended as alternative insecticides to manage GPA in the field, but avoid repeated applications.
Synthetic pyrethroids (Group 4A) form another common insecticide group that have been used for several decades for controlling the GPA in China. Since GPA collected from Beijing and Hebei exhibited extremely high resistance levels to fenvalerate (RF = 675) and deltamethrin (RF = 1,468) in 1990 [
52], pyrethroid resistance became widespread across China by the late 1990s. Surveys from 2011 to 2014 revealed that more than 85.5% of 62 collected field populations had developed a high resistance to beta-cypermethrin [
44] in China. The SEF-R population exhibited extremely high resistance to beta-cypermethrin (RF = 324-fold) and high resistance to esfenvalerate (RF = 79.1-fold). The recommended doses for GPA control ranged from 25 to 90 mg L–1, far below the LC50 values of esfenvalerate (764.4 mg L–1) and beta-cypermethrin (2,371.54 mg L–1) (
Table 1), which led to control failures of GPA after applications of the two pyrethroids at recommended concentrations. After the SEF-R was reared for 15 generations without exposure to any insecticides, the GPA maintained a 91.3-fold resistance to beta-cypermethrin (LC50= 667.97 mg L–1) (
Figure 1B). The two pyrethroids were not advised for the control of the GPA in this field.
IRAC classified pymetrozine [
53] and flonicamid [
54] to Group 9 “Selective homopteran feeding blockers”, with pymetrozine in Subgroup 9B and flonicamid in Subgroup 9C. The two compounds target the chordotonal organ, resulting in an inhibitory effect on aphid feeding, and they offer attractive alternatives to conventional insecticides due to a lack of cross-resistance. Pymetrozine has entered the Chinese market for about two decades and flonicamid has also been used for ten years. According to the Arthropod Pesticide Resistance Database [
55], there have been no reports of field resistance to flonicamid in GPA worldwide. In Greece, RFs to pymetrozine ranged from 1.1 to 2.2) in 2018 [
56], even lower than those found in 2010 (1.0–3.7 and 1.0–6.0 in clones from peach and tobacco, respectively) [
57] although ~10-year passed. Astonishingly, the SEF-R evolved 5.8- and 29.89-fold resistance to flonicamid and pymetrozine, respectively. However, 205.9-fold resistance to pymetrozine has been documented in A. gossypii collected from Cheongju, South Korea [
58], and 1,100-fold resistance to pymetrozine has also been detected in B. tabaci [
59]. These reports suggest, possible resistance development in GAP to the two compounds. Furthermore, resistance to pymetrozine decreased quickly during the first four generation and remained stable in the subsequent nine generations. Owing to the quick increase in susceptibility due to the suspended application of pymetrozine and the relative susceptibility of the SEF-R to flonicamid, they could be used as alternative to other insecticides. The insecticides spirotetramat (Group 23), afidopyropen (Group 9D), and cyantraniliprole (Group 28), were recently used as aphicides against GPA in China. Spirotetramat has a new mode of action and is from a class of chemicals called tetramic acids, which act as Lipid Biosynthesis Inhibitors. Afidopyropen is a derivative of pyripyropene A that has a new mode of action as a chordotonal organ Transient Receptor Potential Vanilloid-type channel modulator, and cyantraniliprole is a member of the bisamides class of insecticides with an excellent aphicidal activity [
23]. Except for the detection of high levels of resistance to spirotetramat in two populations of GPA collected from Queensland (170.7-fold resistance in Alloway171 and 105.6-fold resistance in Osborne171) [
25], there has been little evidence of variation in susceptibility to afidopyropen, or cyantraniliprole. No obvious resistance has been observed to either afidopyropen or cyantraniliprole, and only low resistance to spirotetramat (RF=8.1) was identified in SEF-R. The three insecticides should be applied as alternatives to insecticides provoking high resistance levels.
Broflanilide, a new meta-diamide pesticide that acts on the γ-aminobutyric acid receptor in insect pests through a novel mechanism of action, is highly active against a broad range of pest, including lepidoptera, coleoptera, and homoptera [
60]. Resistance to broflanilide in the GPA has not been documented in previous studies. To date, broflanilide is not permitted for the control of GPA, but it has been registered for the management of Plutella xylostella and Phyllotreta striolata (Fabricius) on cruciferous vegetables in China. We confirmed that growers used broflanilide spray more than 5 times per year. Selection pressures owing to the extensive use of broflanilide to control P. xylostella and P. striolata might lead to the evolvement of resistance in the GPA in the field. In the future, risk assessment of resistance in GPA is necessary prior to the registration of broflanilide use on cabbage.
Resistance monitoring data indicated that SEF-R developed multifold resistance to different insecticides through multiple mechanisms. The study of resistance in GPA has showed the importance of target-site mechanisms in developing high levels of resistance by affecting insecticide binding to receptors [
9]. The R81T mutation in the nAChR gene have been proved responsible for high resistance to neonicotinoids and also found conferred cross-resistance to other insecticide classes function as agonists on the nAChR, such as sulfoxaflor [
39,
49,
51] and butenolides [
48] since it was first detected in France in GPA [
61]. This resistance mutation has been detected in SEF-R, which has also been observed in Chinese populations of GPA [
62]. We speculated the R81T mutation conferred the resistance in SEF-R to all five tested Group 4 compounds, and sulfoxaflor and butenolides are not used as alternative to rotate with neonicotinoids. Otherwise, the SEF-R possessed kdr mutation (L1014F) in sodium channel with 79.35% genotypic frequencies of RS heterozygotes, but super-kdr mutation was not found. Similaryly, Tang [
44] found that the L1014F kdr mutation was detected in 35 of 45 field popoulations with 2.7-100% genotypic frequencies of RS heterozygotes in China. The kdr mutation also meant it difficult to select other pyrethroid compounds, such as deltamethrin, bifenthrin and fenpropathrin to control the aphid in the field.Enhanced detoxification by pesticide detoxification enzymes is an important mechanism for reducing the susceptibility to insecticides. PBO and DEF had useful synergistic effects with imidacloprid and beta-cypermethrin, but DEM had only useful effects with cypermethrin. Synergism of PBO with neonicotinoids or permethrins has been comprehensively documented in field-collected resistant GPA and shows usefulness with imidacloprid [
8,
63] and cyhalothrin [
63], which is consistent with our results. Increases associated with MFO, CarE, and GST were also observed. The results suggested the possible involvement of the metabolic resistance mechanisms of the three detoxification enzymes in SEF-R. This agreed with reports on field resistant populations collected from Hunan and Chongqing Provinces [
46,
60] and Korea [
64].
To further investigate the functions of the MFOs that were related to multi-resistance in SEF-R, we chose eight P450 genes for an expression analysis. Compared with in FFJ-S, seven P450 genes were more than 2-fold upregulated in the SEF-R, with CYP6CY3 (3.63-fold), CYP6CY23 (3.13-fold), CYP380C40 (2.51-fold), CYP6CY9 (2.46-fold), and CYP4G15 (5.24-fold) being significantly overexpressed. The overexpression of CYP6CY3 (10–75-fold) is responsible for neonicotinoid resistance in the GPA, which confers a moderate level of resistance to neonicotinoids [
9,
12]. Additionally, two GPA clones with different copy numbers of CYP6CY3 display low levels (3–20-fold) of resistance to imidacloprid and thiamethoxam, but the two clones survive and reproduce well on Brassica napus seedlings grown from seed treated with commercial levels of neonicotinoids [
14] in a large-scale semi-field trial. The knockdown of CYP6CY3 in SEF-R increased the susceptibility of SEF-R to imidacloprid, thiacloprid, and thiamethoxam, but it did not affect the susceptibility to sulfoxaflor and flupyradifurone, which is consistent with the results of Nakao et al. [
13]. The CYP6CY23 also ungraduated greatly in SEF-R, but it showed no capacity to metabolize nicotine [
65], its function in GPA is needed to study.
The reduced penetration of insecticide is an additional mechanism of resistance in the imidacloprid-resistant 5191A clone of GPA [
9,
12]. Additionally, CYP4G15 was significantly overexpressed in SEF-R. The CYP4 clan in insects is involved in the last step of cuticular hydrocarbon biosynthesis [
5,
66]. Although there was no evidence of the CYP4 clan’s involvement in insecticides resistance in the GPA before our study, the overexpression levels of CYP4EM10 [
67] and CYP4PR1 [
69] in Triatoma infectants and CYP4G19 in Blattella germanica [
70] are associated with pyrethroid resistance, and the increased expression levels of CYP4C68 and CYP4G70 are associated with imidacloprid resistance in Diaphorina citri Kuwayama [
71]. Further studies need to confirm the function of CYP4G15 in the GPA. The overexpression of CYP380C40 (21–76-fold) and UDP-glucuronosyltransferase UGT344P2 (6–33-fold) has been identified in field populations of the GPA in Australia, and this proved partially responsible for sulfoxaflor resistance in the aphid [
16]. The expression levels of CYP380C40 and UGT344P2 were estimated to be 2.51- and 18-fold, respectively, and the overexpression might be responsible for the 32.4-fold increase in resistance to sulfoxaflor in SEF-R compared with in FFJ-S. Although the function of some P450s in GPA has been well identified, co-up-regulation of several P450 genes are often observed and how these genes work together is unclear.
Additionally, we evaluated the differences in the expression levels of two esterase genes E4, and THEM between SEF-R and FFJ-S and found that they all had increased expression levels in SEF-R. The carboxylesterase, E4, mainly confer resistance to carbamates and organophosphates [
11]. Resistance to these two kinds of insecticides in the GPA are widely distributed throughout China, although resistance levels have not been investigated in SEF-R. THEM is a type of thioesterase that catalyzes the cleavage of thioester bonds present in a wide range of glutathione, acyl-carrier proteins, fatty acyl-coenzyme A substrates, and other cellular molecules [
72]. Thioester-containing proteins may be an outcome of different environments and selective pressures that trigger the innate immune response. [
73,74]. The GST gene was not upregulated in SEF-R, which indicated that it was not associated with insecticide resistance. Antioxidant mechanisms were also investigated. A biochemical analysis revealed that the activities of SOD and POD were higher, and the CAT activity was lower, in SEF-R than in FFJ-S. Higher SOD and POD activities have also been observed in another type of field resistance, which is consistent with our results, but the CAT activity also increased in that resistant strain, which is contrary to our result [
64]. Quantitative PCR showed increased expressions of SOD and POD, and a decreased expression of CAT, in SEF-R, consistent with their activity levels. The protective enzymes often work as part of a basic immune response.
Previous resistance monitoring surveys in China put emphasis on conventional insecticides, such as neonicotinoids and pyrethroids. Some new ingredients with aphicidal activities are often not included in monitoring program. SEF-R exhibited resistance to some new introduced or unregistered insecticides, such as broflanilide, which meant careful monitoring of susceptibility to these compounds is necessary. The risk of insecticide resistance including establishment of susceptibility baseline needs to be assessed when a new ingredient will be registered for new target pest insects in China. However, other target pest insects are not considered in the risk assessment. Based on our results, neonicotinoids and pyrethroids are not suitable to control SEF-R, but afidopyropen, and cyantraniliprole are good alternative to the two kinds of insecticdes effectively control SEF-R. Flonicamid and spirotetramat are also recommended in SEF-R management in order to maintain diversity of modes of action and implement resistance management in practice, which is critical in preventing or delaying resistance development. Other control tactics, such as biological control, should integrate with insecticides for enhanced protection from aphid damage and reduced insecticide application.