Relationship between Knockdown Resistance, Metabolic Detoxification and Organismal Resistance to Pyrethroids in Anopheles sinensis

Abstract

Anopheles sinensis is the most important vector of malaria in Southeast Asia, including China. Currently, the most effective measure to prevent malaria transmission relies on vector control through the use of insecticides, primarily pyrethroids. Extensive use of insecticides poses strong selection pressure on mosquito populations for resistance. Resistance to insecticides can arise due to mutations in the insecticide target site (target site resistance), which in the case of pyrethroids is the para-type sodium channel gene, and/or the catabolism of the insecticide by detoxification enzymes before it reaches its target (metabolic detoxification resistance). In this study, we examined deltamethrin resistance in An. sinensis from China and investigated the relative importance of target site versus metabolic detoxification mechanisms in resistance. A high frequency (>85%) of nonsynonymous mutations in the para gene was found in populations from central China, but not in populations from southern China. Metabolic detoxification as measured by the activity of monooxygenases and glutathione S-transferases (GSTs) was detected in populations from both central and southern China. Monooxygenase activity levels were significantly higher in the resistant than the susceptible mosquitoes, independently of their geographic origin. Stepwise multiple regression analyses in mosquito populations from central China found that both knockdown resistance (kdr) mutations and monooxygenase activity were significantly associated with deltamethrin resistance, with monooxygenase activity playing a stronger role. These results demonstrate the importance of metabolic detoxification in pyrethroid resistance in An. sinensis, and suggest that different mechanisms of resistance could evolve in geographically different populations.

About the Authors:

Daibin Zhong

Affiliation: Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America

Xuelian Chang

Affiliations Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America, Department of Pathogen Biology, Nanjing Medical University, Nanjing, China, Department of Pathogen Biology, Bengbu Medical College, Anhui, China

Guofa Zhou

Affiliation: Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America

Zhengbo He

Affiliation: Institute of Entomology and Molecular Biology, College of Life Sciences, Chongqing Normal University, Chongqing, People's Republic of China

Fengyang Fu

Affiliation: Institute of Entomology and Molecular Biology, College of Life Sciences, Chongqing Normal University, Chongqing, People's Republic of China

Zhentian Yan

Affiliation: Institute of Entomology and Molecular Biology, College of Life Sciences, Chongqing Normal University, Chongqing, People's Republic of China

Guoding Zhu

Affiliations Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America, Division of Malaria Control, Jiangsu Institute of Parasitic Diseases, Wuxi, China

Tielong Xu

Affiliation: National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention, Shanghai, People's Republic of China

Mariangela Bonizzoni

Affiliation: Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America

Mei-Hui Wang

Affiliation: Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America

Liwang Cui

Affiliation: Department of Entomology, the Pennsylvania State University, University Park, Pennsylvania, United States of America

Bin Zheng

* E-mail: [email protected] (BZ); [email protected] (BC); [email protected] (GY)

Affiliation: National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention, Shanghai, People's Republic of China

Bin Chen

* E-mail: [email protected] (BZ); [email protected] (BC); [email protected] (GY)

Affiliation: Institute of Entomology and Molecular Biology, College of Life Sciences, Chongqing Normal University, Chongqing, People's Republic of China

Guiyun Yan

* E-mail: [email protected] (BZ); [email protected] (BC); [email protected] (GY)

Affiliation: Program in Public Health, College of Health Sciences, University of California Irvine, Irvine, California, United States of America

Introduction

The mosquito Anopheles sinensis is the major malaria vector in China and other Southeast Asian countries [1]–[6]. The World Health Organization (WHO) defines vector control as one of the four basic and most effective measures to prevent malaria transmission. Vector control relies primarily on the use of insecticides through insecticide-impregnated bednets and indoor residual spray [7]. Chemical insecticides have been extensively used since the 1940s to control mosquito vectors, and historically, four major categories of insecticides have been sequentially applied: organochlorines, organophosphates, carbamates and pyrethroids [8]–[10]. In the past decade, pyrethroids have become the preferred choice among the currently WHO approved compounds because of their low toxicity to humans, high efficacy against mosquito vectors and short residual action [11].

Synthetic pyrethroids, particularly permethrin and deltamethrin, are widely used for malaria vector control worldwide [12], [13], where the use of impregnated bednets has been shown to be very effective in reducing malaria transmission [14], [15]. Pyrethroids are also used for agricultural pest control. In China, extensive use of insecticides for agricultural pest and public health disease vector control has posed intensive selection pressure on mosquitoes [12], [16]. High levels of resistance to deltamethrin have been reported in An. sinensis populations from China and Korea [12], [17].

It is crucial to understand the mechanism of insecticide resistance in order to establish reliable resistance diagnosis methods and aid in resistance management. Pyrethroid resistance in malaria vectors has been mostly studied in the major African malaria vector, An. gambiae [18]–[23]. Two major resistance mechanisms are known: a) mutations in the para-type sodium channel gene, the target site of pyrethroids, causing a change in affinity between the insecticide and its binding site that reduces sensitivity to the insecticide; and b) metabolic detoxification of pyrethroids before they reach their target site by detoxification enzymes (the multigene family of P450 monooxygenases, esterases, and glutathione-S-transferases) [24].

Several point mutations in the para-type sodium channel gene have been identified in Anopheles mosquitoes in association with pyrethroid resistance. The most common mutation conferring knockdown resistance (kdr) is a mutation at position 1014 causing a change from Leucine to either Phenylalanine (L1014F) or Serine (L1014S) [17], [25]–[31]. In Musca domestica [32] and Haematobia irritans [33] a mutation at position 918 (M918T) of the para gene was identified as being correlated with very high levels of pyrethroid resistance (super kdr phenotype) [32], [34]. Metabolic detoxification has been linked to elevated levels of P450 monooxygenases, esterases, and glutathione-S-transferases (GSTs) [24], [35]–[37].

In the present study we assessed the level of pyrethroid resistance in An. sinensis from five localities in southern and central China and examined the frequency of kdr mutations and the levels of detoxification enzymes in the phenotyped mosquitoes.

Materials and Methods

Study Sites

Anopheles sinensis larvae and pupae were collected from five different geographical sites in China, including Yunnan (Mengla and Yuanyang), Hunan (Liuyang), Hubei (Wuxue) and Jiangsu (Sihong) provinces (Fig. 1). The Yunnan and Hubei sites are malaria endemic. The Hunan site was historically endemic. Jiangsu has been experiencing sporadic P. vivax malaria outbreaks. Since the 1980s in China, deltamethrin impregnated nets (ITNs) and indoor residual spraying (IRS) have been widely used against mosquitoes, as well as other insects [38]. Rice is the major agricultural crop in these study sites. There are two or three rice harvests each year in the Yunnan sites, but only one harvest per year in the Hunan, Hubei and Jiangsu sites. Due to severe insect pest damage to the rice crop, insecticide use for rice pest control has been very intensive, with several rounds of sprays in each growing season. Pyrethroids are commonly used for agricultural pest control in study areas, but other insecticides such as organic phosphates and carbamates are also being used [16].

[Figure omitted. See PDF.]

Figure 1. A map of China showing the distribution of mosquito sampling sites.

Site name: 1, Mengla county (Yunnan Province); 2, Yuanyang county (Yunnan Province); 3, Liuyang county (Hunan Province); 4, Wuxue county (Hubei Province); and 5, Sihong city (Jiangsu Province).

https://doi.org/10.1371/journal.pone.0055475.g001

Mosquito Sampling

In each site, we collected about 1,000 larvae and pupae from more than 100 larval habitats, mostly rice fields and small ponds with aquatic plants, using the standard 350 ml dippers. Mosquito larvae and pupae were reared to adults for bioassay in order to eliminate the influence of mosquito age and blood feeding history on resistance. After identifying mosquitoes to species morphologically, larvae and pupae of An. sinensis were kept in plastic storage bins and transported to a local rearing shelter close to the collection site. Larvae and pupae were reared to adults under local environmental conditions. The emerged female adults were used in deltamethrin resistance bioassay [39], enzyme activity assay and kdr allele genotyping assays. Mosquito sampling and bioassays were conducted in July–August 2010 for the two sites in Yunnan province. Collections from the Hunan, Hubei and Jiangsu study sites were conducted in July–August 2011. A laboratory susceptible strain that has been maintained in the insectary of the Jiangsu Institute of Parasitic Diseases in Wuxi, China, for more than 10 years with no insecticide exposure was used as the reference susceptible strain.

Insecticide Susceptibility Bioassay

Female adult mosquitoes directly reared from field-collected larvae and pupae, 3 days post adult emergence, were individually tested for susceptibility to deltamethrin, using the standard WHO tube bioassay with 0.05% deltamethrin test papers [39]. For each population, about 100–150 female mosquitoes were tested. Paraffin oil-treated papers without insecticide were used as control. The knockdown time of individual mosquitoes was recorded as minutes within one hour. After 1-hr exposure to the diagnostic concentration of deltamethrin (0.05%), mosquitoes were transferred to recovery cups and maintained on 6% sucrose solution for 24 hrs. Mosquito survival status was examined every 3 hrs, and dead mosquitoes were immediately collected and preserved at −20°C prior to metabolic enzyme analysis. Mosquitoes were considered dead if they were motionless, even when they were mechanically stimulated, following the method of González Audino [40]. This is reasonable because motionless mosquitoes after bioassay exposure rarely recovered and survived. One leg of each mosquito was also removed and preserved individually in 95% alcohol for subsequent DNA analysis. The intact tissue of the remaining mosquito body was immediately tested for metabolic enzyme activity [41], [42].

Metabolic Enzyme Activity Assays

Individual females were homogenized in a 1.5-ml tube with 200 µl of phosphate KPO4 buffer (0.25 M, pH 7.2) and then diluted by adding 1,300 µl of phosphate buffer. The tube was mixed and centrifuged, and the supernatant was used to test the activity of GSTs and monooxygenases [25], [43], [44]. GST and monooxygenase assays were carried out in duplicate. A total of 300 µl reduced glutathione solution and 300 µl cDNB solution was added to 300 µl of mosquito supernatant. The absorbance was measured immediately, and then every minute for 5 minutes at 340 nm, using 0.25 M KPO4 buffer as the negative control. The change in absorbance was converted to enzyme activity using the CDNB extinction coefficient of 0.0096 uM−1cm−1 for this reaction, and the estimated protein content of the supernatant was measured by the Bradford method [45]. For the monooxygenase activity assay, a total of 8 µl of the 7-ethoxycoumarine (7-EC) solution was added to 400 µl of mosquito supernatant and samples were incubated at 30°C for 4 hours. The reaction was stopped by the addition of 560 µl of glycine buffer (pH 10.4, 0.1 mM) and sample absorbance was measured at 450 nm using 0.25 M KPO4 buffer as a negative control. The OD values were converted into concentration by using a standard regression based on a serial dilution of 7-hydroxycoumarin and its relevant OD values. The value of monooxygenase content was calculated for each mosquito.

Mosquito DNA Extraction

One leg of each mosquito was used for DNA extraction with the Fast Tissue-to-PCR Kit (Fermentas). Briefly, each mosquito leg was placed at the bottom of a 500 µl Eppendorf tube. A total of 50 µl of tissue lyses solution and 5 µl of protein K solution were added and incubated at 55°C for 20 min, followed by 10 minutes at 95°C. After the incubations, 50 µl of neutralization solution was added and mixed by vortexing. The neutralized tissue was then centrifuged at 14,000 rpm for 10 min. Extracted DNA was stored at 4°C or used immediately for PCR.

Molecular Identification and Detection of kdr Mutation

Molecular identifications of An. sinensis species were conducted by using species-specific primers and amplification of the ITS2 and 28S rDNA regions (D1 and D2) [46]. To determine point mutations of the para-type sodium gene at positions 1014 and 918, we first amplified a 325 bp fragment of the para-type sodium gene including position 1014, and then a 450 bp fragment including position 918 (Table 1). PCR primers were designed based on the An. sinensis sequences of the DIIS6 region of the para-type sodium gene (Genbank acc. no. DQ334052). PCR products were directly sequenced in both forward and reverse directions using the same primers. Direct sequencing was performed by Genewiz Inc (South Plainfield, NJ) on 150 field-derived mosquitoes. Based on data obtained from these individuals, allele specific amplification (AS-PCR) primers were designed to detect four mutations (TTG, TTT, TGT and TTC). To increase allele specific amplification, each primer was designed to incorporate a single base mismatch at third base position from the 3' end terminus. The DNA sequences were used as a gold standard to determine the sensitivity and specificity of the AS-PCR method.

[Figure omitted. See PDF.]

Table 1. Primer sequence of multiplex allele-specific PCR for the detection of mutations in codon 1014 of para-type sodium channel gene and super kdr mutation in codon 918 in Anopheles sinensis.

https://doi.org/10.1371/journal.pone.0055475.t001

Statistical Analysis

WHO (1998) criteria were used to evaluate the resistance or susceptibility of the tested mosquito populations (<80% mortality: resistance, 80–98% mortality: probable resistance; >98% mortality: susceptibility) [39]. Univariate analysis of variance (ANOVA) was conducted using arcsin transformation of mosquito mortality rate to determine among-population differences in mosquito mortality rate in the insecticide susceptibility bioassay. The kdr allele frequency was calculated in each population and statistical differences among populations were examined using the t-test. ANOVA also was used to examine whether monooxygenase and GST activity varied among populations.

Stepwise multivariate logistic regression and linear regression analyses were conducted to determine the role of kdr mutations and the level of detoxification enzymes on deltamethrin resistance. Stepwise regression is a procedure for automatically selecting independent variables in a stepwise manner to maximize the squared multiple correlations coefficient (R2) (for linear regression) or the likelihood (for logistic regression) of the dependent variables against a set of independent variables. In the present study, the dependent variables were the mosquito knockdown time and survival status (dead or survived within 24-hr recovery period) for stepwise multiple linear regression or the stepwise logistic regression, respectively. In all tests, independent variables were the three kdr genotypes and the activity levels of monooxygenases and GSTs. A total of three kdr genotypes were considered: 1) homozygous wild type (TTG/TTG), 2) heterozygous mutations (TTG/TTT, TTG/TTC or TTT/TTC), and 3) homozygous mutations (TTT/TTT, TTC/TTC or TTT/TTC).

Ethics Statement

No specific permits were required for the described field studies. For mosquito collection in rice paddies, oral consent was obtained from field owners in each location. These locations were not protected land, and the field studies did not involve endangered or protected species.

Results

Deltamethrin Susceptibility Bioassay

All five tested populations were resistant to deltamethrin, with mortality rates ranging from 9.8% to 38.8% (Fig. 2A). ANOVA analysis found significant among-population variation in mortality rates (F4, 15 = 97.2, P < 0.0001). In particular, the two populations from Yunnan province in southern China (Mengla and Yuanyang, 29.6%–38.8%) had significantly higher mortality rates than populations from central China (Liuyang, Sihong, and Wuxue, 9.8%–15.0%) (Fig. 2B). The An. sinensis susceptible laboratory strain showed 100% mortality, confirming the quality of resistance test papers.

[Figure omitted. See PDF.]

Figure 2. Mortality rates of Anopheles sinensis in standard WHO tube bioassays.

A: after exposure to 0.05% deltamethrin test papers and 24 hr recovery period. B: pairwise Tukey–Kramer HSD test for statistical significance in mortality rates. The vertical bar stands for standard error of the mean.

https://doi.org/10.1371/journal.pone.0055475.g002

Kdr Gene Sequencing and Allele-specific PCR (AS-PCR) for kdr Genotyping

A 325 bp fragment of the para-type sodium channel gene including position 1014 and a 450 bp fragment including position 918 were sequenced from 30 An. sinensis individuals in each of the five populations. The wild-type kdr codon sequence spanning position 1014 was TTG. Three types of kdr mutations were detected: two (TTT and TTC) leading to a change from Leucine to Phenylalanine and one (TGT) leading to a Leucine to Cysteine substitution. A total of ten genotypes were identified in the five populations (Fig. 3). No mutations spanning position 918 were detected. Based on these sequence data, an AS-PCR protocol for the identification of kdr mutations was developed (Fig. 4). The AS-PCR method showed 98.9% sensitivity and 87% specificity with respect to sequencing results. The AS-PCR method was used for kdr genotyping of 456 additional individuals.

[Figure omitted. See PDF.]

Figure 3. Examples of nucleotide sequence chromatograms of kdr genotypes detected in Anopheles sinensis from China.

The position at codon 1014 of the para-type sodium channel gene is indicated by a rectangle box. A: four types of homozygote genotypes detected; and B: six types of heterozygote genotypes detected (K = G/T; Y = T/C; S = G/C).

https://doi.org/10.1371/journal.pone.0055475.g003

[Figure omitted. See PDF.]

Figure 4. Allele-specific PCR (AS-PCR) for kdr genotyping in Anopheles sinensis.

A: Schematic representation of the primer location and predicted size of the PCR product. Arrow indicated the location of PCR primers in which sequences are reported in Table 1. Primer pair TTG-F and kdr-R amplifies a 141 bp fragment for the wildtype allele (for codon TTG). Primer pair TGT-F/kdr-R yields a 158 bp fragment for L1014C allele (codon TTC). Similarly, primer pair kdr-F/TTT-R and kdr-F/TTC-R leads to amplification of a 221 bp and 241 bp fragment diagnostic to the L1014F allele (for codon TTT and TTC), respectively. B: An example of AS-PCR gel electrophoresis. Lane M: 100 bp DNA ladder; lanes 1–4: homozygote for the codons TTG, TGT, TTT, and TTC, respectively; lanes 5–10: six heterozygous individuals (TTG/TTT, TGT/TTT, TTG/TGT, TTT/TTC, TGT/TTC, and TTG/TTC).

https://doi.org/10.1371/journal.pone.0055475.g004

Distribution of kdr allele Frequencies in An. sinensis Populations

Kdr genotyping was performed on both susceptible (dead) and resistant (alive) mosquitoes recovered from the WHO bioassays for deltamethrin resistance. No kdr mutations were detected in the two populations from southern China (Mengla and Yuanyang) despite high levels of phenotypic resistance. An opposite situation was found in three populations from central China ((Liuyang, Wuxue, and Sihong), where deltamethrin-resistant mosquitoes carried a high frequency (88.5–94.8%) of both the L1014F and the L1014C mutations (Table 2). Kdr mutations were identified also in susceptible mosquitoes from populations in central China, but at a significantly lower frequency than the wild type allele (Paired t-test, t = 7.79, d.f. = 2, P<0.01).

[Figure omitted. See PDF.]

Table 2. Frequencies (in percentage) of kdr alleles in relation to mosquito survival phenotype determined by the deltamethrin susceptibility bioassay in five Anopheles sinensis populations from China.

https://doi.org/10.1371/journal.pone.0055475.t002

Metabolic Enzyme Assay

Monooxygenase and GST activity varied both within and among the five populations tested (Table 3). The average GST activity was 0.23 umoles cDNB/min/mg protein in the laboratory strain and 0.39 umoles cDNB/min/mg protein for the five field populations. The GST activity was significantly higher in the three populations from central China (Liuyang, Wuxue and Sihong) than in the laboratory susceptible strain (F4, 15 = 4.04, P<0.05) (Table 3). The Mengla and Yuanyang population from southern China exhibited similar GST activity to the laboratory strain. The monooxygenase activity was two to three folds higher in the two populations from southern China than in the laboratory population (F4, 15 = 25.2, P<0.0001). The three populations from central China showed comparable monooxygenase activity to the laboratory strain (Table 3).

[Figure omitted. See PDF.]

Table 3. Glutathione s-transferases (GSTs) and monooxygenases activity in Anopheles sinensis female adults in five populations in China when compared with lab susceptible strain.

https://doi.org/10.1371/journal.pone.0055475.t003

Role of kdr Genotypes and Metabolic Detoxification Enzyme in Resistance

We conducted three analyses to determine the relative contribution of the different kdr genotypes and the activity of monooxygenase and GSTs enzymes on deltamethrin resistance. First, we tested whether the frequency of the different kdr genotypes varied between phenotypically -resistant and -susceptible mosquitoes in the three populations from central China (Liuyang, Wuxue and Sihong). The percentage of resistant mosquitoes that were homozygotes for the wild-type kdr genotype was low (1.1%). Total of 44.0% and 54.8% resistant mosquitoes showed kdr genotypes with mutations in a heterozygote or homozygote state, respectively (Fig. 5A). Second, we determined whether the activity of monooxygenases and GSTs were different between phenotypically - resistant and -susceptible mosquitoes. GST activity was not significantly different between resistant and susceptible individuals (Fig. 5B) in all tested populations. Monooxygenase activity was significantly higher in resistant than in susceptible mosquitoes (t-test, P<0.05) (Fig. 5C) in all tested populations.

[Figure omitted. See PDF.]

Figure 5. Relationship between kdr, metabolic detoxifying enzyme activity and pyrethroid resistance in Anopheles sinensis.

Association between kdr genotypes and survival rate of mosquitoes after the standard WHO tube deltamethrin susceptibility bioassay (A), and variation of resistant (survived the bioassay) and susceptible (died in the bioassay) individuals in glutathione s-transferases activity (B) and monooxygenases activity (C). **, P < 0.01.

https://doi.org/10.1371/journal.pone.0055475.g005

Finally, we performed a stepwise multivariate logistic regression and a linear regression analysis which considers kdr genotypes, GSTs and monooxygenase activities simultaneously to minimize the effects of collinearity between the two variables. Stepwise linear regression analyses with knockdown time as the dependent variable identified kdr mutations in either a homozygous or heterozygous state and levels of monooxygenase activity as the significant independent variables (ANOVA F3, 293 = 13.078, P<0.0001) (Table 4). In our model, positive standard coefficients were identified for homozygous kdr mutations and monooxygenase activity, with the latter showing higher standardized coefficient. This indicates a positive effect on resistance, with monooxygenase activity having a stronger impact (Table 4). Consistent with the knockdown time analysis, stepwise logistic regression analysis using survival status as the dependent variable selected homozygous mutations and monooxygenase activity as the significant variables. High odds ratio for monooxygenase activity suggests that mosquitoes with high monooxygenase activity were far more likely to be deltamethrin resistant (Table 5).

[Figure omitted. See PDF.]

Table 4. Summary results of stepwise multivariate regression analysis on knockdown time in Anopheles sinensis.

https://doi.org/10.1371/journal.pone.0055475.t004

[Figure omitted. See PDF.]

Table 5. Stepwise logistic regression for factors significantly associated with deltamethrin resistance in Anopheles sinensis.

https://doi.org/10.1371/journal.pone.0055475.t005

Discussion

This study examined phenotypic resistance to deltamethrin in An. sinensis from China in relation to mutations in the insecticide target site and levels of detoxification enzymes. All five studied populations were highly resistant to deltamethrin, with particularly high levels (>85%) in populations from central China. When target site mutations and levels of detoxification enzymes were tested in phenotypically- resistant and -susceptible mosquitoes, major differences were found between the populations from southern and central China. Interestingly, no kdr mutations were found in the two populations from southern China, but high mutation frequencies (88.5–94.8%) were identified in the three populations from central China. In these populations three kdr genotypes were identified, leading to either a L1014F or a L1014C substitution. We did not detect any L1014S mutation as previously reported in An. sinensis populations from the Republic of Korea and Vietnam [17], [31], [47], [48]. Regarding metabolic detoxification mechanisms of resistance, GST activity was not detected as a significant variable associated with deltamethrin resistance in the multiple regression analyses. On the contrary, monooxygenase enzymes were significantly correlated to phenotypic resistance in all tested populations, with particularly high levels of enzyme activity in populations from southern China. The results show that metabolic detoxification mechanisms are wide-spread in An. sinensis populations from southern and central China, but target site mutations confined to central China suggest that different resistance mechanisms can arise in locally different settings, probably in relation to different selection pressures. Additionally, we noticed higher levels of phenotypic resistance in populations from central China, where both kdr and detoxification mechanisms were identified, suggesting an additive effect of target site mutation and metabolic detoxification on phenotypic resistance.

Yunnan is a typical longitudinal range-gorge region with complex terrain and climate, where great mountains and deep valleys run parallel and connect the Indochinese peninsula and the inner area of western China [49]. The mountain ranges in the longitudinal range-gorge region of Yunnan effectively restrict gene flow in the latitudinal direction. Population genetic studies of An. sinensis revealed the role of the geographical barriers in obstructing gene flow among the mosquito populations [2], [5]. In the mountainous regions of southern China, agricultural crops are diverse and insecticide use has, traditionally, not been intensive. On the other hand, insecticide use has been intense in central China, with several rounds of spray per year [16]. Rice is the most widely planted crop in this region and rice paddies are the major breeding sites of An. sinensis mosquitoes. It will be interesting to test for traces of insecticides in breeding site water as it has been suggested that phenotypic resistance in adults is due to biological influence of the breeding sites on mosquito immature stages [50].

All specimens used in the bioassay were adult mosquitoes reared from larvae collected from the field. Although this was very time consuming and labor intensive, this design allowed using mosquitoes of fixed age (3 days post emergence) and with known physiological status (no blood feeding) for the bioassay, thus minimizing the confounding effects of age or blood-feeding status on resistance phenotype and subsequent metabolic enzyme activity measurement. Because the resistance bioassays were conducted with individual mosquitoes, we could measure the precise knockdown time and resistance status within the 24-hr recovery period. This enabled us to accurately determine the relationship between knockdown time, resistance status, kdr genotypes and monooxygenase and GST activity. Our results indicate that overall knockdown time in bioassays is a good approximation of pyrethroid susceptibility [51], [52]. None of the resistant mosquitoes were knocked down during the 60 min exposure time. This suggests that a complete loss of knockdown is indicative of deltamethrin resistance. On the other hand, susceptible mosquitoes were being knocked-down and showed high levels of variation in the knockdown time. The majority of susceptible mosquitoes were dead within 10 hours after the end of the bioassay. This result is possibly related to the slow killing effect of deltamethrin [53] and suggests phenotypic plasticity, which should be further investigated.

Considering the high pyrethroid resistance in An. sinensis populations in China, resistance management becomes a particularly important issue. Insecticide rotation or synergistic uses of different classes of insecticides to enhance the killing effects are being considered as alternative methods of vector control to reduce the human-vector contact rate and to slow down the spread of resistance genes [54]–[57]. Therefore, resistance to multiple classes of insecticides and the role of selection by insecticides used for both agricultural pest control and public health pest control are also important topics for future research. Such knowledge will help guide rational use of insecticides in both public health and agriculture. Additionally, the identification of widespread metabolic detoxification mechanisms with the confinement of kdr mutations to central China suggests that the development of reliable resistance surveillance tools should focus on monooxygenase metabolic detoxification mechanisms.

Acknowledgments

We thank Ziyong Bai, Zhimin Bai, Lingzhong Kong, Yujiang Xiao, Jihui Duan and Guangquan Huang for assistance with mosquito collection. Two anonymous reviewers provided constructive criticisms and helpful suggestions.

Author Contributions

Conceived and designed the experiments: DZ BC BZ LC GY. Performed the experiments: DZ BC XC ZH FF ZY G. Zhu TX. Analyzed the data: DZ G. Zhou GY. Contributed reagents/materials/analysis tools: BC BZ G. Zhou. Wrote the paper: DZ BC G. Zhu MB MW GY.

Citation: Zhong D, Chang X, Zhou G, He Z, Fu F, Yan Z, et al. (2013) Relationship between Knockdown Resistance, Metabolic Detoxification and Organismal Resistance to Pyrethroids in Anopheles sinensis. PLoS ONE 8(2): e55475. https://doi.org/10.1371/journal.pone.0055475

References

1. Chareonviriyaphap T, Bangs MJ, Ratanatham S (2000) Status of malaria in Thailand. Southeast Asian J Trop Med Public Health 31: 225–237.

2. Jung J, Jung Y, Min G-S, Kim W (2007) Analysis of the population genetic structure of the malaria vector Anopheles sinensis in South Korea based on mitochondrial sequences. Am J Trop Med Hyg 77: 310–315.

3. Jung J, Lee E, Kim W (2006) Isolation and characterization of polymorphic microsatellite markers of Anopheles sinensis, a malaria vector mosquito in the East Asia region. Mol Ecol Notes 6: 1272–1274.

4. Liu X-B, Liu Q-Y, Guo Y-H, Jiang J-Y, Ren D-S, et al. (2011) The abundance and host-seeking behavior of culicine species (Diptera: Culicidae) and Anopheles sinensis in Yongcheng city, people's Republic of China. Parasit Vectors 4: 221.

5. Ma Y, Yang M, Fan Y, Wu J, Ma Y, et al. (2011) Population structure of the malaria vector Anopheles sinensis (Diptera: Culicidae) in China: two gene pools inferred by microsatellites. PLoS ONE 6: e22219.

6. Sinka M, Bangs M, Manguin S, Chareonviriyaphap T, Patil A, et al. (2011) The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic precis. Parasit Vectors 4: 89.

7. WHO (2006) Global Malaria Program: Use of indoor residual spraying for scaling up global malaria control and elimination. WHO position statement. Geneva: World Health Organization. Available: http://www.who.int/malaria/docs/IRS/IRS-position.pdf. Accessed 2009 Apr 06.

8. Kamgang B, Marcombe S, Chandre F, Nchoutpouen E, Nwane P, et al. (2011) Insecticide susceptibility of Aedes aegypti and Aedes albopictus in Central Africa. Parasit Vectors 4: 79.

9. Lima E, Paiva M, de Araujo A, da Silva E, da Silva U, et al. (2011) Insecticide resistance in Aedes aegypti populations from Ceara, Brazil. Parasit Vectors 4: 5.

10. Liu N, Xu Q, Zhu F, Zhang LEE (2006) Pyrethroid resistance in mosquitoes. Insect Sci 13: 159–166.

11. WHO (2005) Malaria control today: current WHO recommendations. Geneva: Working document March, World Health Organization.

12. Cui F, Raymond M, Qiao CL (2006) Insecticide resistance in vector mosquitoes in China. Pest Manag Sci 62: 1013–1022.

13. Overgaard HJ, Sandve SR, Suwonkerd W (2005) Evidence of Anopheline mosquito resistance to agrochemicals in northern Thailand. Southeast Asian J Trop Med Public Health 36 Suppl 4152–157.

14. Sexton JD (1994) Impregnated bed nets for malaria control: biological success and social responsibility. Am J Trop Med Hyg 50(6 Suppl): 72–81.

15. Zuzi L, Mancheng Z, Yuguang W, Binglin Z, Guangyu L, et al. (1989) Trial of deltamethrin impregnated bed nets for the control of malaria transmitted by Anopheles sinensis and Anopheles anthropophagus. Am J Trop Med Hyg 40: 356–359.

16. Cui F, Tan Y, Qiao CL (2007) Filariasis vector in China: insecticide resistance and population structure of mosquito Culex pipiens complex. Pest Manag Sci 63: 453–458.

17. Kim H, Baek JH, Lee W-J, Lee SH (2007) Frequency detection of pyrethroid resistance allele in Anopheles sinensis populations by real-time PCR amplification of specific allele (rtPASA). Pestic Biochem Physiol 87: 54–61.

18. Chandre F, Darrier F, Manga L, Akogbeto M, Faye O, et al.. (1999) Status of pyrethroid resistance in Anopheles gambiae sensu lato. Bulletin of the World Health Organization, 77.

19. Fanello C, Kolaczinski JH, Conway DJ, Carnevale 427 P, Curtis CF (1999) The kdr pyrethroid resistance gene in Anopheles gambiae: tests of non-pyrethroid insecticides and a new detection method for the gene. Parassitologia 41: 323–326.

20. Kawada H, Dida GO, Ohashi K, Komagata O, Kasai S, et al. (2011) Multimodal pyrethroid resistance in malaria vectors, Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus s.s. in Western Kenya. PLoS ONE 6: e22574.

21. Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, et al. (1998) Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Bio 7: 179–184.

22. Ranson H, Abdallah H, Badolo A, Guelbeogo W, Kerah-Hinzoumbe C, et al. (2009) Insecticide resistance in Anopheles gambiae: data from the first year of a multi-country study highlight the extent of the problem. Malaria J 8: 299.

23. Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, et al. (2008) Distribution of knock-down resistance mutations in Anopheles gambiae molecular forms in west and west-central Africa. Malaria J 7: 74.

24. Hemingway J, Ranson H (2000) Insecticide resistance in insect vectors of human disease. Annu Rev Entomol 45: 371–391.

25. Chen H, Githeko AK, Githure JI, Mutunga J, Zhou G, et al. (2008) Monooxygenase levels and knockdown resistance (kdr) allele frequencies in Anopheles gambiae and Anopheles arabiensis in Kenya. Med Vet Entomol 45: 242–250.

26. Enayati AA, Vatandoost H, Ladonni H, Townson H, Hemingway J (2003) Molecular evidence for a kdr-like pyrethroid resistance mechanism in the malaria vector mosquito Anopheles stephensi. Med Vet Entomol 17: 138–144.

27. Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, et al. (2000) Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Bio 9: 491–497.

28. Singh O, Dykes C, Lather M, Agrawal O, Adak T (2011) Knockdown resistance (kdr)-like mutations in the voltage-gated sodium channel of a malaria vector Anopheles stephensi and PCR assays for their detection. Malaria J 10: 59.

29. Syafruddin D, Hidayati A, Asih P, Hawley W, Sukowati S, et al. (2010) Detection of 1014F kdr mutation in four major Anopheline malaria vectors in Indonesia. Malaria J 9: 315.

30. Verhaeghen K, Van Bortel W, Roelants P, Backeljau T, Coosemans M (2006) Detection of the east and west African kdr mutation in Anopheles gambiae and Anopheles arabiensis from Uganda using a new assay based on FRET/Melt curve analysis. Malaria J 5: 16.

31. Verhaeghen K, Van Bortel W, Trung HD, Sochantha T, Keokenchanh K, et al. (2010) Knockdown resistance in Anopheles vagus, An. sinensis, An. paraliae and An. peditaeniatus populations of the Mekong region. Parasit Vectors 3: 59.

32. Williamson MS, Martinez-Torres D, Hick CA, Devonshire AL (1996) Identification of mutations in the housefly para-type sodium channel gene associated with knockdown resistance (kdr) to pyrethroid insecticides. Mol Genet Genomics 252: 51–60.

33. Guerrero FD, Jamroz RC, Kammlah D, Kunz SE (1997) Toxicological and molecular characterization of pyrethroid-resistant horn flies, Haematobia irritans: identification of kdr and super-kdr point mutations. Insect Biochem Mol Biol 27: 745–755.

34. Jamroz RC, Guerrero FD, Kammlah DM, Kunz SE (1998) Role of the kdr and super-kdr sodium channel mutations in pyrethroid resistance: correlation of allelic frequency to resistance level in wild and laboratory populations of horn flies (Haematobia irritans). Insect Biochem Mol Biol 28: 1031–1037.

35. Brooke BD, Kloke G, Hunt RH, Koekemoer LL, Temu EA, et al. (2001) Bioassay and biochemical analyses of insecticide resistance in southern African Anopheles funestus (Diptera: Culicidae). Bull Entomol Res 91: 265–272.

36. David J-P, Strode C, Vontas J, Nikou D, Vaughan A, et al. (2005) The Anopheles gambiae detoxification chip: A highly specific microarray to study metabolic-based insecticide resistance in malaria vectors. Proc Natl Acad Sci USA 102: 4080–4084.

37. Muller P, Warr E, Stevenson BJ, Pignatelli PM, Morgan JC, et al. (2008) Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genetics 4: e1000286.

38. Li ZZ, Zhang MC, Wus YG, Zhong BL, Lin GY, et al. (1989) Trial of deltamethrin impregnated bed nets for the control of malaria transmitted by Anopheles sinensis and Anopheles anthropophagus. Am J Trop Med Hyg 40: 356–359.

39. WHO (1998) Test procedures for insecticide resistance monitoring in malaria vectors, bioefficacy and persistence of insecticides on treated surfaces. Report of an informal consultation, 1981, WHO/CDS/CPC/MAL/9812.

40. Gonzalez Audino P, Vassena C, Barrios S, Zerba E, Picollo MI (2004) Role of enhanced detoxication in a deltamethrin-resistant population of Triatoma infestans (Hemiptera, Reduviidae) from Argentina. Mem Inst Oswaldo Cruz 99: 335–339.

41. Bouvier J-C, Cuany A, Monier C, Brosse V, Sauphanor B (1998) Enzymatic diagnosis of resistance to deltamethrin in diapausing larvae of the codling moth, Cydia pomonella (L.). Arch Insect Biochem Physiol 39: 55–64.

42. Picollo MI, Vassena C, Santo Orihuela P, Barrios S, Zaidemberg M, et al. (2005) High resistance to pyrethroid insecticides associated with ineffective field treatments in Triatoma infestans (Hemiptera: Reduviidae) from Northern Argentina. J Med Entomol 42: 637–642.

43. Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, et al. (2003) Pyrethroid and DDT cross-resistance in Aedes aegypti is correlated with novel mutations in the voltage-gated sodium channel gene. Med Vet Entomol 17: 87–94.

44. Penilla RP, Rodriguez AD, Hemingway J, Torres JL, Arredondo-Jimenez JI, et al. (1998) Resistance management strategies in malaria vector mosquito control. Baseline data for a large-scale field trial against Anopheles albimanus in Mexico. Med Vet Entomol 12: 217–233.

45. Bradford MM (1976) A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72: 248–254.

46. Joshi D, Park MH, Saeung A, Choochote W, Min GS (2010) Multiplex assay to identify Korean vectors of malaria. Mol Ecol Resour 10: 748–750.

47. Kang S, Jung J, Lee S, Hwang H, Kim W (2012) The polymorphism and the geographical distribution of the knockdown resistance (kdr) of Anopheles sinensis in the Republic of Korea. Malaria J 11: 151.

48. Tan W-L, Wang Z-M, Li C-X, Chu H-L, Xu Y, et al. (2011) First report on co-occurrence knockdown resistance mutations and susceptibility to beta-cypermethrin in Anopheles sinensis from Jiangsu Province, China. PLoS ONE 7: e29242.

49. He D, Wu S, Peng H, Yang Z, Ou X, et al. (2005) A study of ecosystem changes in longitudinal range-gorge region and transboundary eco-security in Southwest China. Adv Earth Sci 20: 338–344 (in Chinese)..

50. Djouaka R, Bakare A, Coulibaly O, Akogbeto M, Ranson H, et al. (2008) Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria. BMC Genomics 9: 538.

51. Elissa N, Mouchet J, Riviere F, Meunier JY, Yao K (1993) Resistance of Anopheles gambiae s.s. to pyrethroids in Cote d'Ivoire. Ann Soc Belg Med Trop 73: 291–294.

52. Kang W, Gao B, Jiang H, Wang H, Yu T, et al. (1996) Tests for possible effects of selection by domestic pyrethroids for resistance in culicine and anopheline mosquitoes in Sichuan and Hubei, China. Ann Trop Med Parasitol 89: 677–684.

53. Khambay BPS, Jewess PJ (2010) Pyrethroids. In Insect 528 Control Biological and Synthetic Agents. Edited by Gilbert LI and Gill SS: Oxford, UK: Elsivier. 1–29.

54. Corbel V, Stankiewicz M, Bonnet J, Grolleau F, Hougard JM, et al. (2006) Synergism between insecticides permethrin and propoxur occurs through activation of presynaptic muscarinic negative feedback of acetylcholine release in the insect central nervous system. Neurotoxicology 27: 508–519.

55. Farenhorst M, Knols BGJ, Thomas MB, Howard AFV, Takken W, et al. (2010) Synergy in efficacy of fungal entomopathogens and permethrin against West African insecticide-resistant Anopheles gambiae mosquitoes. PLoS ONE 5: e12081.

56. Mathenge EM, Gimnig JE, Kolczak M, Ombok M, Irungu LW, et al. (2001) Effect of permethrin impregnated nets on exiting behavior, blood feeding success, and time of feeding of malaria mosquitoes (Diptera: Culicidae) in western Kenya. J Med Entomol 38: 531–536.

57. Munhenga G, Masendu H, Brooke B, Hunt R, Koekemoer L (2008) Pyrethroid resistance in the major malaria vector Anopheles arabiensis from Gwave, a malaria-endemic area in Zimbabwe. Malaria J 7: 247.