- Open Access
DDT and pyrethroid resistance in Anopheles arabiensis from South Africa
Parasites & Vectorsvolume 6, Article number: 229 (2013)
Pyrethroid resistance has been well documented in Anopheles arabiensis, one of the major African malaria vectors, and the predominant malaria vector in South Africa.
In this study, the genetic basis of pyrethroid resistance in a selected laboratory strain of An. arabiensis from South Africa was investigated using a custom-made microarray, known as the An. gambiae detoxification chip.
A large number of P450 genes were over-transcribed, as well as a suite of redox genes and glutathione S-transferases. The five genes that showed the highest level of gene transcription when compared with an insecticide susceptible strain were: CYP6AG2, CYPZ1, TPX2, CYPZ2 and CYP6P1.
Permethrin resistance in South African An. arabiensis is associated with increased transcription of multiple genes, and a large proportion of these genes were also previously recorded as over-transcribed in another An. arabiensis strain selected for resistance to DDT with cross-resistance to deltamethrin. The deltamethrin resistance developed de novo in the DDT-selected strain and is most likely due to increased transcription of those genes associated with DDT resistance. However, of particular interest was the fact that the strain selected for resistance to pyrethroids did not develop de novo resistance to DDT. These differences are compared and discussed.
Pyrethroids are the most commonly used insecticides for malaria vector control. They are most often used for indoor residual spraying, and are the only class of insecticide approved for treatment of bednets due to their low levels of mammalian toxicity [1–3]. However, pyrethroid resistance in African malaria vectors, including Anopheles arabiensis, is well documented and has been reported in many African countries.
The Global Plan for Insecticide Resistance Management in Malaria Vectors (GPIRM)  provides a generic guideline for managing insecticide resistance where it occurs. Strategies for preventing and managing insecticide resistance are imperative to ensure that the limited number of insecticides available for vector control are protected. As part of this strategy, insecticide resistance should be monitored through routine surveillance and the mechanisms of insecticide resistance should be characterised where possible . Information on insecticide resistance mechanisms in malaria vectors can be obtained using several techniques that vary in sophistication and cost. The more costly methodologies are not feasible in many African countries where the disease burden is the greatest.
Resistance to pyrethroids is mainly due to two mechanisms; enhanced enzyme detoxification or reduced target site sensitivity. Cytochrome oxidases (P450s), glutathione S-transferases (GSTs) and esterases are generally associated with enzymatic detoxification of insecticides, while mutations in the voltage-gated sodium channel (commonly referred to as knockdown resistance or “kdr”) effectively decrease target site sensitivity [4, 5]. The occurrence of kdr mutations is easily detected through the use of a polymerase chain reaction (PCR) assay  but this mechanism is not clearly correlated with all instances of pyrethroid resistance, especially deltamethrin resistance . Identifying metabolic based resistance is more challenging and costly, especially when using microarrays . The use of microarray technology is particularly helpful as it allows for rapid analysis of numerous genes simultaneously, and as such provides substantial genetic information relating to a particular phenotype.
Genes associated with pyrethroid resistance include the esterase, GST and P450 super families [8–13]. Metabolic detoxification of pyrethroids has been demonstrated in An. gambiae[12, 14, 15], An. arabiensis[8, 10, 16] and An. funestus (reviewed by ). Anopheles arabiensis is a major vector of malaria in Africa including South Africa [18, 19]. Historically, this species was considered to be susceptible to all classes of insecticide. In 1996 an extensive survey of insecticide susceptibility in An. arabiensis was conducted in three areas in South Africa (KwaZulu-Natal [KZN], Limpopo Province and the Kruger National Park) . Bioassays revealed that samples of An. arabiensis from all three areas were susceptible to DDT and deltamethrin .
Current malaria vector control in South Africa is based on a mosaic approach using DDT and pyrethroids . These interventions are the mainstay of mosquito vector control in South Africa even though An. arabiensis is both exophilic and endophilic. In 2003, An. arabiensis sampled from two areas in KZN were found to be resistant to DDT, but susceptible to deltamethrin .
Subsequent studies were initiated to investigate the molecular basis for the reported DDT resistance in KZN An. arabiensis. The use of laboratory strains to study resistance mechanisms carries the advantage of excluding factors that can confound the analysis (e.g. exposure to agricultural insecticides, effect of temperature, larval diet, etc). Laboratory based DDT selection pressure was exerted on this population over a period of several years and resulted in the fixation of the L1014F kdr mutation . In addition to target site resistance, metabolic detoxification was also identified and numerous gene transcripts were over-expressed and associated with the resistant phenotype .
Three years after the discovery of DDT resistance, pyrethroid resistance was also reported in wild caught An. arabiensis from KZN . This population was colonised in 2005 and selected for resistance to permethrin in order to characterise the resistance mechanisms involved. Pyrethroid resistance in this population (named KWAG-perm) was primarily associated with increased activity of monooxygenases (P450’s) . Analysis of six specific P450 genes showed increased transcription for three genes in association with the pyrethroid resistant phenotype .
The aim of this study was to evaluate the levels of transcription of detoxification enzymes associated with the South African permethrin resistant strain (KWAG-perm) using a high throughput microarray analysis to supplement the information obtained by Munhenga and Koekemoer . In addition to this, the transcript profiles between the pyrethroid-selected (KWAG-perm, permethrin resistant) and DDT-selected (MBN-DDT, DDT and pyrethroid resistance) strains were compared to identify both strain specific transcripts as well as transcripts shared between these two strains. The implications of these results for malaria control are discussed.
Pyrethroid resistant strain (KWAG-perm)
The An. arabiensis permethrin resistant strain, KWAG-perm, and the equivalent susceptible strain, KWAG-base, were used for this experiment. KWAG-base originates from Mamfene in KwaZulu-Natal (KZN) province, South Africa. This population was characterised as permethrin resistant at the time of colonisation . However, in the absence of selection by permethrin exposure, the laboratory strain reverted back to being susceptible (after a period of one year). When permethrin selection was reintroduced, resistance rapidly returned (within four generations) and this strain was named KWAG-perm (Table 1) . Mosquitoes were reared under standard insectary conditions (26°C ± 2°C, 70-80% relative humidity, and 12:12 light:dark cycle).
DDT resistant strain (MBN-DDT)
Details on this strain can be found in Nardini et al.. Briefly, MBN-base was colonised from Mamfene, KZN, and a sub-strain was selected for resistance to DDT. The selected strain, MBN-DDT, is resistant to DDT (0% mortality) as well as permethrin (4%), deltamethrin (34%), bendiocarb (77.5%) and propoxur (65.3%) through cross-resistance mechanisms (Table 1). The strain is homozygous for the L1014F kdr mutation which is mainly responsible for the observed DDT and permethrin resistance . The deltamethrin resistance is mainly due to enhanced monooxygenase enzymatic detoxification .
RNA extractions and preparation of amplified mRNA
Pyrethroid resistant strain (KWAG-perm)
RNA extractions were carried out according to Nardini et al.. Briefly, female mosquitoes were maintained on a 10% sugar solution and RNA was extracted from 15 mosquitoes three days after emergence, representing one biological repeat. RNA was extracted for three biological repeats using the PicoPure® RNA Isolation Kit (Arcturus) according to supplied methodology (with the inclusion of a DNase treatment, as described in the user manual). RNA was quantified using a NanoDrop and quality was assessed using an Agilent 2100 Bioanalyzer. Extracted RNA was stored at -70°C. In order to obtain larger quantities of mRNA for microarrays, the RiboAmp™ RNA Amplification Kit (Arcturus) was used. Amplified mRNA was quantified using a NanoDrop.
Preparation of microarrays and data analyses
Pyrethroid resistant strain (KWAG-perm)
Three independent biological repeats were used, and for each biological repeat, two technical repeats were performed that included dye swaps in order to compensate for dye bias. Microarrays were prepared as described in Nardini et al..
The arrays were scanned using the Genepix 4000B scanner (Molecular Devices, USA). The photomultiplier tube (PMT) values were adjusted to give a pixel ratio of approximately 1. Spot quality and background intensities were examined and corrected using Genepix Pro 6.0 software (Axon Instruments, USA). Saturated features were excluded from the analysis.
Gene expression data were analysed using Limma version 2.12.0 (Bioconductor)  in R, version 2.8.0 (http://cran.r-project.org/bin/windows/base/old/2.8.0/). Raw intensity values for each spot were calculated, and then background corrected by the method “normexp” with an offset of 50 . The corrected intensity values were transformed to log-ratios and then normalised. Composite Loess was used for within array normalisation. In this method, control spots and features, per sub-array, are used for producing non-linear, best-fit lines . The use of control spots ensures that the resulting best-fit line is not biased by differential expression of genes. Conversely, the use of all genes for normalisation improves stability with respect to the number of spots and, most importantly, provides flexibility in terms of print-tip group trends that might be observed where sub-array Loess curves are used . Print-tip peculiarities were present in some slides, hence the choice of the normalisation method. The “Aquantile” method was used for between array normalisation. Genes with adjusted p-values ≤ 0.05 and fold-changes (FCs) ≥ 2.0 were considered to be statistically significant. These data have been deposited into Vectorbase (https://www.vectorbase.org).
Quantitative real-time PCR (qPCR)
Pyrethroid resistant strain (KWAG-perm)
Real-time PCR was carried out in order to validate the outcome of the microarray experiments. RNA was extracted from 15 three day old An. arabiensis females (represents one biological repeat, and three biological repeats were prepared) that had been supplied with 10% sugar solution. TRI-Reagent® solution (Sigma-Aldrich) was used according to the supplied methodology with a DNase treatment included. Samples were reverse-transcribed into cDNA using the QuantiTect® Reverse Transcription Kit (Qiagen). The cDNA was stored at -70°C until required for qPCR.
Two genes, CYP6AG2 and TPX2, were evaluated using real-time PCR. Beacon Designer™ (Premier Biosoft) or Invitrogen’s free online primer design tool, OligoPerfect™ Designer, were used to design primers. These were based on An. gambiae sequence information (CYP6AG2: GB AY745225; TPX2: TIGR TC48596). The reference gene used for qPCR was based on that of Munhenga and Koekemoer , where studies were conducted in the laboratory on the same mosquito strains. These authors reported that 18S rRNA showed the most stable expression of the six potential reference genes tested. The forward and reverse primer sequences for 18S were the same as those of Munhenga and Koekemoer . Each PCR reaction comprised 12.5 μl IQ™ SYBR super-mix (Bio-Rad), 4 μl primer (primer concentration optimised for each primer set), 1 μl cDNA (100 ng/μl) and nuclease-free water to a final volume of 25 μl. Cycling conditions were as follows: 93°C for 3 minutes, followed by 35 cycles of (94°C/20 seconds, ×°C/25 seconds, 72°C/30 seconds), and a final extension at 72°C for 10 minutes (followed by melt curve analysis). Primer sequences and annealing temperatures are described in Table 2. Data were analysed using the Pfaffl  method. For each gene of interest, the PCR product was sequenced in both directions, so that the presence of the correct product was confirmed (in addition to melt-curve analysis). Where necessary (for small amplicons i.e. < 110 base pairs), samples were cloned and then sequenced.
KWAG-base and KWAG-perm have been evaluated for insecticide resistance on an ongoing basis in the laboratory. The base strain showed very little resistance to permethrin (97.8% mean mortality 24 hours post exposure to 0.75% permethrin), while the resistant strain showed a high frequency of survivors on permethrin (42% mean mortality 24 hours post exposure to 0.75% permethrin) .
Microarray experiments indicated that 29 genes were over-transcribed according to the criteria outlined above. Most of these were P450 genes (55%), followed by redox genes (21%), and GSTs (14%). A group of additional genes, such as cytochrome c, a ribosomal gene (RPS26) and a receptor protein (GPR npy 3) accounted for 10% of the over-transcribed genes. The four genes with the highest transcript abundance were CYP6AG2, CYP6Z1, TPX2 and CYP6Z2 (Table 3), in order of statistical significance (i.e. adjusted p-value).
Two significantly over-transcribed genes, CYP6AG2 and TPX2, were selected for qPCR validation of microarray data. The two other “top” genes (CYP6Z1 and CYP6Z2) were already evaluated with qPCR . CYP6AG2 was the most significantly over-transcribed gene, and TPX2 is of interest because the TPX genes have been found to be important in other studies reporting on insecticide resistance [8, 10]. The genes CYP6AG2 and TPX2 produced fold change values of 2.9 and 4.5 respectively (Figure 1). The FC value of CYP6AG2 and TPX2 were comparable to those obtained by microarray analysis.
Munhenga and Koekemoer  used qPCR to evaluate transcription of specific genes (identified through published literature) that have been implicated in permethrin resistance in An. arabiensis. These included CYP6Z1, CYP6Z2, CYP6Z3, CYP6M2, CYP6P3 and CYP4G16. Of these, five appeared on the list of over-transcribed genes in this microarray study (Table 3) and three of these genes [CYP6Z1 (4.7-fold), CYP6Z2 (1.7-fold) and CYP6M2 (2.2-fold)] were significantly over-transcribed in the qPCR study . The FC values of all six genes analysed by Munhenga and Koekemoer  were compared with the FC values obtained here for the same genes using microarray analysis (Figure 2). The FC values obtained by microarray analyses were compared with those obtained using qPCR for each gene using a t-test. The FC values of only three of the genes were significantly different to the FC values obtained by microarray analysis. These were CYP6Z2 (p = 0.0251), CYP6Z3 (p = 0.0219) and CYP6M2 (p = 0.0092).
Transcript profile comparison between the DDT-selected strain and the pyrethroid-selected strain (this study)
Both MBN-DDT and KWAG-perm originate from Mamfene (KZN). However, they were placed under different insecticide selection pressures (depending on the insecticide resistance phenotypes present at the time of colonisation). The list of over-transcribed probes obtained in the permethrin-resistant strain (KWAG-perm, pyrethroid resistance only) was compared with that obtained in the DDT selected strain (MBN-DDT, pyrethroid and DDT resistance ) (Figure 3). In MBN-DDT, kdr was closely associated with DDT and permethrin resistance, while deltamethrin resistance was largely based on enhanced enzymatic detoxification. Five GST transcripts were unique in the DDT-selected strain (MBN-DDT), while 14 transcripts were unique to the permethrin-selected strain (KWAG-perm). Fifteen gene transcripts were shared between these of which 10 (66%) belonged to the P450-enzyme family.
The use of controlled insecticide selection in the laboratory aids in the study of changes in resistance status and its effect on certain adaptive features . A list of 29 over-transcribed genes according to the criteria used (adjusted p-value ≤ 0.05; FC ≥ 2) was produced for the KWAG-perm strain. Most of the genes that were over-transcribed belong to the cytochrome P450 superfamily, supporting the findings of previous studies [16, 22]. The four most statistically significant P450 genes included CYP6AG2, CYP6Z1, CYP6Z2 and CYP6P1, all of which belong to the CYP6 family and are well-known for their role in insecticide resistance [14, 29]. The most significant gene, CYP6AG2, is not well characterised and it would be valuable to determine the functionality of this enzyme. This gene is located on the chromosome arm 2R, a region that is not typically linked to pyrethroid resistance. Another gene in this family, CYP6AG1, was identified in a wild deltamethrin resistant An. arabiensis population from Cameroon . CYP6Z1, CYP6Z2 and CYP6M2 showed increased transcription in this study as well as in the study by Munhenga and Koekemoer .
Pyrethroids are also known to cause oxidative stress by inducing lipid peroxidation, protein oxidation, and depletion of reduced glutathione . This effectively increases the toxicity of the insecticide. A system of enzymes, including the superoxide dismutases (SODs), catalases, peroxidases, cytochrome C and GSTs, is present and provides defense against these reactive oxygen species (ROS) [31, 32]. If ROS are not metabolised, they damage important compounds such as lipids, proteins, nucleic acids and carbohydrates, and ultimately cause cell death . A number of redox enzymes showed elevated levels of transcription in KWAG-perm, which may be associated with resistance to pyrethroids. These included thioredoxin peroxidases (TPX2 and TPX4), glutathione peroxidases (GPX1 and GPX3), superoxide dismutases (SOD1 and SOD2) and cytochrome C. Superoxide radicals are converted to hydrogen peroxide and oxygen by SODs, and hydrogen peroxide is converted to water and oxygen by catalases; or to water, by peroxidases .
Four GSTs (GSTU1, GSTD1-3, GSTS1-1, GSTS1-2) were also over-transcribed in KWAG-perm. Over-transcription of GSTs has been observed before in pyrethroid resistant An. arabiensis, An. gambiae[7, 9, 34] and in other insects [35–37]. In conjunction with glutathione, GSTs function as antioxidants by limiting peroxidation and by limiting (termination) “free-radical cascades” (see ). In addition, the GSTs are able to bind to pyrethroids and provide protection by sequestration .
In the recent study of Nardini et al., the authors identified genes that were over-transcribed in a DDT and pyrethroid resistant strain (MBN-DDT). These genes were mainly found to be associated with the observed cross resistance to deltamethrin (the use of synergists clarified the roles of the over-transcribed genes, and implicated kdr in DDT and permethrin resistance). The gene transcripts that were identified and which are shared between the permethrin-selected (KWAG-perm) and DDT-selected strains (MBN-DDT), might explain the protection against deltamethrin in the DDT selected strain. In the process of selecting for DDT resistance it appears that there was also selection for increased transcription of genes (n = 15) that coincidentally protect against pyrethroids as well. However, the converse was not true and the pyrethroid selected strain did not develop DDT resistance. This implies that DDT selection pressure can also result in the de novo development of resistance to pyrethroids. Very few reports of DDT resistance without resistance to pyrethroids are available (e.g. ), which supports this hypothesis.
It is possible that the genes that are unique in each instance (KWAG-perm versus MBN-DDT) are more important for resistance to either type I (e.g. GSTD2) or type II pyrethroid (CYPZ1), while those that are common are likely to play a role in resistance to both kinds of pyrethroid (e.g. CYP6M2). Those genes that are shared between the pyrethroids and DDT selected strain might be able to metabolise both classes of insecticides (DDT as well as pyrethroids). This is specifically true for CYP6M2, but to our knowledge similar metabolic studies are not available for the other 10 genes. The role of genes conferring carbamate resistance identified in the DDT-selected strain still needs to be clarified . However, P450s have been shown to provide protection against both pyrethroids and carbamates in An. funestus.
The challenge with detoxification based resistant mechanisms is that it is not possible to develop a single molecular assay for the detection of resistance (as in the case of target-site mediated resistance, for example, kdr detection). Furthermore, in target-site resistance mechanisms, genotype usually correlates with phenotypic outcome. In the case of enzymatic resistance, this is not necessarily true. For example, three GSTs were over-transcribed in the strain here (the activity of GSTs is known to confer DDT resistance in mosquitoes [41, 42]), but no DDT resistance is present in KWAG-perm. In addition, different genes were over-transcribed in the permethrin and deltamethrin resistant strains.
Microarrays are extremely useful for identifying a large number of genes that are associated with a particular resistance phenotype and in particular, transcripts that are unique to a particular phenotype. However, this technology is extremely expensive and not available in most African countries and an alternative to this approach is needed. Due to this problem, the results from the costly microarray approach used in this study were compared against those obtained by a less expensive qPCR study using the same An. arabiensis strain . Selection of genes to be included in this qPCR study were identified through published literature [7, 9, 10, 14, 34]. This approach was used successfully by Munhenga and Koekemoer  who found ≥1.5 FC in transcription in four of the six genes that were selected for analysis. Five of the six genes from Munhenga and Koekemoer  were also identified through this microarray study. However, CYP6AG2 and TPX2 were novel in this study.
These data suggest that relatively inexpensive qPCR studies can successfully be used to identify increased transcription of specific metabolic genes, especially in countries that do not have access to advanced molecular systems eg. microarray analysis. However, many transcripts might be overlooked using qPCR due to the specificity of the technique. Microarray studies will in future expand the number of transcripts available for qPCR screening.
The GPIRM has initiated a “five pillar” approach against insecticide resistance. One of these pillars is to gain knowledge on mechanisms of insecticide resistance and on the effects of current insecticide resistance management strategies . The present study and similar studies, contribute to achieving this goal.
Permethrin-resistant An. arabiensis from South Africa
Linear Models in Microarray Analysis
DDT and pyrethroid resistant An. arabiensis from South Africa
Quantitative real-time PCR
Reactive oxygen species
World Health Organization.
World Health Organization: Global Plan for Insecticide Resistance Management in Malaria Vectors. 2012
Coosemans M, Carnevale P: Malaria vector control: a critical review on chemical methods and insecticides. Ann Soc Belg Med Trop. 1995, 75: 13-31.
Liu N, Xu Q, Zhu F, Zhang L: Pyrethroid resistance in mosquitoes. Insect Sci. 2006, 13: 159-166. 10.1111/j.1744-7917.2006.00078.x.
Hemingway J, Hawkes NJ, McCarroll L, Ranson H: The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004, 34: 653-665. 10.1016/j.ibmb.2004.03.018.
Bass C, Nikou D, Donnelly MJ, Williamson MS, Ranson H, Ball A, Vontas J, Field LM: Detection of knockdown resistance (kdr) mutations in Anopheles gambiae: a comparison of two new high-throughput assays with existing methods. Malaria J. 2007, 6: 111-10.1186/1475-2875-6-111.
Brooke BD, Koekemoer LL: Major effect genes or loose confederations? The development of insecticide resistance in the malaria vector Anopheles gambiae. Parasit Vectors. 2010, 3: 74-10.1186/1756-3305-3-74.
David JP, Strode C, Vontas J, Nikou D, Vaughan A, Pignatelli PM, Louis C, Hemingway J, Ranson H: The Anopheles gambiae detoxification chip: a highly specific microarray to study metabolic-based insecticide resistance in malaria vectors. Proc Natl Acad Sci USA. 2005, 102: 4080-4084. 10.1073/pnas.0409348102.
Nardini L, Christian RN, Coetzer N, Ranson H, Coetzee M, Koekemoer LL: Detoxification enzymes associated with insecticide resistance in laboratory strains of Anopheles arabiensis of different geographic origin. Parasit Vectors. 2012, 5: 113-10.1186/1756-3305-5-113.
Müller P, Donnelly MJ, Ranson H: Transcription profiling of a recently colonised pyrethroid resistant Anopheles gambiae strain from Ghana. BMC Genom. 2007, 8: 36-10.1186/1471-2164-8-36.
Müller P, Chouaïbou M, Pignatelli P, Etang J, Walker ED, Donnelly MJ, Simard F, Ranson H: Pyrethroid tolerance is associated with elevated expression of antioxidants and agricultural practice in Anopheles arabiensis sampled from an area of cotton fields in Northern Cameroon. Mol Ecol. 2007, 17: 1145-1155. 10.1111/j.1365-294X.2007.03617.x.
Vontas J, David JP, Nikou D, Hemingway J, Christophides GK, Louis C, Ranson H: Transcriptional analysis of insecticide resistance in Anopheles stephensi using cross-species microarray hybridization. Insect Mol Biol. 2007, 16: 315-324. 10.1111/j.1365-2583.2007.00728.x.
Müller P, Warr E, Stevenson BJ, Pignatelli PM, Morgan JC, Steven A, Yawson AE, Mitchell SN, Ranson H, Hemingway J: Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genet. 2008, 4: e1000286-10.1371/journal.pgen.1000286.
Christian RN, Strode C, Ranson H, Coetzer N, Coetzee M, Koekemoer LL: Microarray analysis of a pyrethroid resistant African malaria vector, Anopheles funestus, from Southern Africa. Pestic Biochem Physiol. 2011, 99: 140-147. 10.1016/j.pestbp.2010.11.010.
Nikou D, Ranson H, Hemingway J: An adult-specific CYP6 P450 gene is overexpressed in a pyrethroid-resistant strain of the malaria vector, Anopheles gambiae. Gene. 2003, 318: 91-102.
Ranson H, Paton MG, Jensen B, McCarroll L, Vaughan A, Hogan JR, Hemingway J, Collins FH: Genetic mapping of genes conferring permethrin resistance in the malaria vector, Anopheles gambiae. Insect Mol Biol. 2004, 13: 379-386. 10.1111/j.0962-1075.2004.00495.x.
Mouatcho JC, Munhenga G, Hargreaves K, Brooke BD, Coetzee M, Koekemoer LL: Pyrethroid resistance in a major African malaria vector Anopheles arabiensis from Mamfene, northern KwaZulu-Natal, South Africa. S Afr J Sci. 2009, 105: 127-131.
Coetzee M, Koekemoer LL: Molecular systematics and insecticide resistance in the major African malaria vector Anopheles funestus. Annu Rev Entomol. 2013, 58: 393-412. 10.1146/annurev-ento-120811-153628.
Coetzee M, Craig M, le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000, 16: 74-77. 10.1016/S0169-4758(99)01563-X.
Maharaj R, Mthembu DJ, Sharp BL: Impact of DDT re-introduction on malaria transmission in KwaZulu-Natal. S Afr Med J. 2005, 95: 871-874.
Gericke A, Govere JM, Durrheim DN: Insecticide susceptibility in the South African malaria mosquito Anopheles arabiensis (Diptera: Culicidae). S Afr J Sci. 2002, 98: 205-208.
Hargreaves K, Hunt RH, Brooke BD, Mthembu J, Weeto MM, Awolola TS, Coetzee M: Anopheles arabiensis and An. quadriannulatus resistance to DDT in South Africa. Med Vet Entomol. 2003, 17: 417-422. 10.1111/j.1365-2915.2003.00460.x.
Munhenga G, Keoekemoer LL: Differential Expression of cytochrome P450 genes in a laboratory selected Anopheles arabiensis colony. Afr J Biotech. 2011, 10: 12711-12716.
Smyth GK: Limma: linear models for microarry data. 2005, New York: Springer
Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.
van Heerden J, Walford S-A, Shen A, Illing N: A framework for informed normalization of printed microarrays. S Afr J Sci. 2007, 103: 381-390.
Smyth GK, Speed TP: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Martins AJ, Ribeiro CD, Bellinato DF, Peixoto AA, Valle D, Lima JB: Effect of insecticide resistance on development, longevity and reproduction of field or laboratory selected Aedes aegypti populations. PLoS One. 2012, 7: e31889-10.1371/journal.pone.0031889.
Scott JG: Cytochromes P450 and insecticide resistance. Insect Biochem Mol Biol. 1999, 29: 757-777. 10.1016/S0965-1748(99)00038-7.
Vontas JG, Small GJ, Hemingway J: Glutathione S-transferases as antioxidant defence agents confer pyrethroid resistance in Nilaparvata lugens. Biochem J. 2001, 357: 65-72. 10.1042/0264-6021:3570065.
Bauer H, Kanzok SM, Schirmer RH: Thioredoxin-2 but not thioredoxin-1 is a substrate of thioredoxin peroxidase-1 from Drosophila melanogaster: isolation and characterization of a second thioredoxin in D. Melanogaster and evidence for distinct biological functions of Trx-1 and Trx-2. J Biol Chem. 2002, 277: 17457-17463. 10.1074/jbc.M200636200.
James RR, Xu J: Mechanisms by which pesticides affect insect immunity. J Invertebr Pathol. 2012, 109: 175-182. 10.1016/j.jip.2011.12.005.
Fridovich I: The biology of oxygen radicals. Science. 1978, 201: 875-880. 10.1126/science.210504.
Djouaka RF, Bakare AA, Coulibaly ON, Akogbeto MC, Ranson H, Hemingway J, Strode C: 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 Genom. 2008, 9: 538-10.1186/1471-2164-9-538.
Hemingway J, Dunbar SJ, Monro AG, Small GJ: Pyrethroid resistance in German cockroaches (Dictyoptera: Blattelidae): resistance levels and underlying mechanisms. J Econ Entomol. 1993, 86: 1631-1638.
Lagadic L, Cuany A, Bergé JB, Echaubard M: Purification and partial characterization of glutathione S-transferases from insecticide-resistant and lindane-induced susceptible Spodoptera littoralis (Boisd.) larvae. Insect Biochem Mol Biol. 1993, 23: 467-474. 10.1016/0965-1748(93)90054-V.
Kostaropoulos I, Papadopoulos AI, Metaxakis A, Boukouvala E, Papadopoulou-Mourkidou E: Glutathione S-transferase in the defence against pyrethroids in insects. Insect Biochem Mol Biol. 2001, 31: 313-319. 10.1016/S0965-1748(00)00123-5.
Namountougou M, Simard F, Baldet T, Diabate A, Ouedraogo JB, Martin T, Dabire RK: Multiple insecticide resistance in Anopheles gambiae s.l. populations from Burkina Faso, West Africa. PLoS One. 2012, 7: e48412-10.1371/journal.pone.0048412.
Mitchell SN, Stevenson BJ, Muller P, Wilding CS, Egyir-Yawson A, Field SG, Hemingway J, Paine MJ, Ranson H, Donnelly MJ: Identification and validation of a gene causing cross-resistance between insecticide classes in Anopheles gambiae from Ghana. Proc Natl Acad Sci USA. 2012, 109: 6147-6152. 10.1073/pnas.1203452109.
Kloke RG, Nhamahanga E, Hunt RH, Coetzee M: Vectorial status and insecticide resistance of Anopheles funestus from a sugar estate in southern Mozambique. Parasit Vectors. 2011, 4: 16-10.1186/1756-3305-4-16.
LA Prapanthadara HJ, Ketterman AJ: Partial purification and characterisation of glutathione s-transferases involved in DDT resistance from the mosquito Anopheles gambiae. Pest Biochem Physiol. 1993, 47: 119-133. 10.1006/pest.1993.1070.
DE Grant DF, Hammock BD: Glutathione-S-transferase isozymes in Aedes aegypti - purification, characterization, and isozyme-specific regulation. Insect Biochem. 1991, 21: 421-433. 10.1016/0020-1790(91)90009-4.
The authors acknowledge the following funding institutions: National Research Foundation Scarce Skills Scholarships and DAAD (LN), National Research Foundation and the National Health Laboratory Service-Research Trust to LLK, partial funding from the Department of Science and Technology, National Research Foundation Research Chair Initiative grant to Maureen Coetzee. The authors also acknowledge Maureen Coetzee and Basil Brooke for comments on the manuscript. The grant holders acknowledge that opinions, findings and conclusions or recommendations expressed in any publication generated by the NRF supported research are that of the author(s), and that the NRF accepts no liability whatsoever in this regard.
The authors declare that they have no competing interests.
LN performed all the experimental work, analyses and drafted the MS. RC provided technical support for the microarray work and reviewed the MS. NC assisted with the microarray data analysis and reviewed the MS. LK conceived the project, provided funding and reviewed the MS. All authors read and approved the final manuscript.