- Open Access
Rapid evolution of pyrethroid resistance prevalence in Anopheles gambiae populations from the cities of Douala and Yaoundé (Cameroon)
Malaria Journalvolume 14, Article number: 155 (2015)
The adaptation of malaria vectors to urban areas is becoming a serious challenge for malaria control. The study presents the evolution of pyrethroid resistance in mosquito populations from the cities of Douala and Yaoundé between 2010 and 2013.
Susceptibility tests to permethrin and deltamethrin were carried out with two- to four-day old unfed Anopheles gambiae sensu lato adults raised from larvae collected from the field. Mosquitoes resistant to permethrin and deltamethrin and control were screened to detect the presence of the kdr alleles using the TaqMan assays. Mosquitoes belonging to the An. gambiae complex were subjected to PCR assays designed for species and molecular forms identifications. The genomic region containing the upstream of intron-1 of the voltage-gated sodium channel was sequenced and compared between mosquitoes originating from different breeding habitats.
Anopheles gambiae s.l. specimens collected from the city of Douala were all Anopheles coluzzii. In Yaoundé, both An. gambiae and An. coluzzii were recorded. A rapid decrease of mosquito mortality to permethrin and deltamethrin was recorded between 2010 and 2013 in the two cities. The mortality rate varied from 80.3 to 22.3% and 94.4 to 59.7% for permethrin and deltamethrin, respectively. Both kdr alleles L1014F and L1014S were recorded. The frequency of kdr alleles increased rapidly over the study period, varying from 44 to 88.9% in Yaoundé and from 68 to 81% in Douala. The sequencing of a 1,228 bp region of intro-1 of the voltage-gated sodium channel revealed the presence of five different haplotypes. A high number of these haplotypes were recorded in An. coluzzii samples. No evidence for a recent selective sweep on intron-1 sequence within samples originating from different breeding habitat was detected using Fu’s and Tajima Fs statistics.
The present study supports rapid evolution of pyrethroid resistance in vector populations from the cities of Douala and Yaoundé and calls for immediate action to fight against the increasing prevalence of pyrethroid-resistant mosquitoes.
In the absence of vaccine, malaria prevention relies almost exclusively on vector control and pyrethroids are the most commonly used insecticides for vector control . Their rapid knockdown effect, high remanence and low toxicity for humans at operational doses make them ideal compounds for both bed-net impregnation and indoor residual spraying (IRS). Pyrethroids have as target the insect nervous system and act by inducing a long-term disruption of the sodium permeability which alters the properties of the sodium channel and causes the insect death . However, rapid evolution of pyrethroid resistance resulting from increased use of these compounds in both public health and agriculture is affecting the efficiency of control measures . During the past decades, the use of pyrethroid-based interventions, IRS and/or long-lasting insecticidal nets (LLINs), dramatically increased across the continent and contributed in the reduction of malaria morbidity and mortality [4,5]. There is concern that the expansion of pyrethroid resistance could jeopardize current momentum of reduction and elimination of malaria . Two main physiological mechanisms confer pyrethroid resistance, these are: target site insensitivity and metabolic resistance . However target site resistance remains the most studied because of easier means of assessment [8-10]. Mutations of kdr gene, both leucine to phenylalanine substitution (L1014F) and leucine to serine substitution (L1014S) are the two mutations inducing target site resistance in insects [11,12]. These two mutations are now largely distributed across the continent, yet differences in the prevalence of these resistance alleles have been reported between Anopheles gambiae s.l. species and between sites [13-15]. The genotyping of the kdr locus and partial sequencing of the upstream intron-1 of the voltage-gated sodium channel gene of mosquitoes collected across sub-Saharan Africa also provided evidence for the high genetic diversity of this gene . Although the importance of kdr mutations as a stand-alone mechanism conferring pyrethroid resistance is still subject to debate , the current evolution of pyrethroid resistance and the diversity of resistance mechanisms represent a great challenge that requires further attention.
In Cameroon, despite the introduction of pyrethroid-treated nets (ITNs) in the late 1990s, a fast evolution of resistance was recorded across the country a few years after implementation [17,18]. The increased use of pesticides in agriculture, the expansion of urban agriculture and pollution in urban settings, were also reported to influence the prevalence and distribution of insecticide resistance across the country [19,20]. In the cities of Douala and Yaoundé, An. gambiae larvae emerging from cultivated and organically polluted sites were found to be 325 and 193 times, respectively, more tolerant to deltamethrin than An. gambiae Kisumu larvae . Despite the fact that pyrethroid resistance is dynamic, the influence of local selective pressure induced either by agricultural practices, urban pollution or even increase use of ITNs on the evolution of resistant alleles remain poorly understood. Analysis conducted with mosquitoes originating from different type of habitats in the city of Yaoundé suggested a different resistance pattern in mosquito populations according to their breeding habitat origins . The present study assessed the evolution of insecticide resistance from 2010 to 2013 in the cities of Douala and Yaoundé and the variability of intron-1 of the voltage-gated sodium channel gene in mosquitoes originating from different larval habitats.
The study took place in the cities of Yaoundé (3°51’N11°30’E) and Douala (3°48’N10°08’E), the two major urban cities in Cameroon. Douala and Yaoundé have a population approaching three million inhabitants each . These cities are situated within the Congo-Guinean phytogeographic zone characterized by a typical equatorial climate, with two rainy seasons extending from March to June and from September to November. Douala is situated near the Atlantic coast 1 m above sea level and receives over 3,500 mm of rainfall annually, whereas the annual average rainfall in Yaoundé is 1,700 mm. Yaoundé is located inland 250 km east of Douala. The city is situated 800 m above sea level and is surrounded by many hills. The study was conducted under the ethical clearance N° 216/CNE/SE/09 delivered by the Cameroon National Ethics Committee Ref N° IORG0006538-IRB00007847-FWA00016054.
Insecticide susceptibility assays
Bioassays were carried out using WHO test kits for adult mosquitoes . The following diagnostic concentrations of insecticides were tested: 0.75% permethrin and 0.05% deltamethrin.
Susceptibility tests were carried out with two to four days old unfed An. gambiae s.l. adults raised from larvae collected in breeding sites. Batches of 20 to 25 mosquitoes per tube were exposed to impregnated papers for one hour. The number of mosquito knockdown was recorded every 5 minutes during exposure. After exposure, mosquitoes were supplied with glucose solution as food, and mortality was recorded 24 hours post-exposure. Tests with untreated papers were systematically run as controls. Mortality rate in tested samples was corrected using Abbot formula  when the mortality rate of control was between 5 and 20%. WHO criteria [24,26] were used to evaluate the resistance and susceptibility status of the tested mosquito population. The resistance status was indicated by a mortality rate below 90% while mortality rates greater than 97% were indicative of susceptibility, and mortality rates between 90 and 97% suggested increased tolerance, but resistance should be confirmed.
To screen for the presence of the kdr alleles (L1014F and L1014S) conferring resistance to DDT and pyrethroids, DNA extracted from individuals exposed to insecticide was tested using the hot ligation oligonucleotide assay (HOLA) of Lynd et al.  and/or the TaqMan assay .
Anopheline larvae were identified morphologically using the Gillies and Coetzee keys . Mosquitoes belonging to the An. gambiae complex were subjected to PCR assays designed for species and molecular forms identifications . Genomic DNA used for molecular analysis was extracted from dessicated adult mosquitoes according to either Livak  or Cornel  protocols.
Intron sequence determination
Mosquito used for sequencing included specimen originating from organically polluted, unpolluted and agricultural cultivated sites. The genomic region containing the upstream of intron-1 was PCR amplified and sequenced as presented in Weill et al. . PCR products were then purified using the QIAquick purification Kit (Qiagen, Hilden, Germany). The sample was afterward sequenced on both strands. Sequences from each specimen were analysed manually to detect polymorphic sites using the software BIOEDIT V: 126.96.36.199. Sequence differences in multiple alignments were analysed using Clustal O.
The software DnaSP 5.10  was used to estimate the genetic polymorphism of intron-1 in samples originating from different breeding habitats. Parameters measured were respectively: the number of haplotypes, haplotype diversity, the number of segregating sites, Fu’s Fs statistics, Tajima statistics.
For bioassay analyses, comparison of percentages was performed using the Chi square test. Estimates of mortality rates and the 95% confidence interval were determined using the software MedCalc V188.8.131.52.
Evolution of Anopheles gambiae s.l. susceptibility to permethrin and deltamethrin in different sites across Cameroon
A map describing the evolution of An. gambiae s.l. susceptibility to permethrin and deltamethrin from 2003 onwards after the first countrywide free distribution of ITNs is presented (Figure 1). Data presented here have been extracted from the present study and from different published reports [17,19,33-35]. For a better understanding of the pattern of An. gambiae s.l. susceptibility it is important to keep in mind that, the northern half of the country situated in dry savanna and humid savanna areas have as dominant malaria vector Anopheles arabiensis whereas, An. gambiae and An. coluzzii are the main vectors in the southern half of the country (highlands, coastal and equatorial forest regions).
Annual evolution of Anopheles gambiae susceptibility to pyrethroids in Yaoundé and Douala between 2010 and 2013
To assess the evolution of mosquito susceptibility to permethrin and deltamethrin, susceptibility assays were conducted at least twice each year at different seasons from 2010 to 2013. Data from 2010 have already been published , but are included here for comparison. Apart of mosquitoes collected in cultivated sites, mosquitoes collected in either polluted or non-polluted sites in both Douala and Yaoundé, showed high susceptibility to permethrin at the beginning of the study with a mortality rate ranging from 75 to 85% in 2010. Four years after, less than 40% of An. gambiae s.l. specimens in both cities were detected still susceptible to permethrin, this showing the total establishment of permethrin resistance in these species populations (Figure 2). Deltamethrin mortality decreased significantly between the beginning and the end of the study with mortality rates varying from 95% in 2010 to about 60% in 2013 in both cities (Figure 2). When quantifying the decrease in vector susceptibility between 2010 and 2013, it appeared that An. gambiae s.l. was four to eight times and three to eleven times less susceptible to permethrin and deltamethrin, respectively, in 2013 compared to 2010 (Table 1).
Evolution of kdr alleles prevalence in mosquito populations
Both kdr alleles L1014F and L1014S were recorded in Douala and Yaoundé. The frequency of kdr alleles increased rapidly over the four-year period in the city of Yaoundé, varying from 44% in 2010 to over 88.9% in 2013. In Douala, a similar trend was recorded with the allele frequency varying from 68 to 81% during the same period (Figure 3). Of the 1,449 mosquitoes processed during the study period, eight were found carrying the L1014S allele (seven in Douala and one in Yaoundé) and 1,200 (82.8%) were found with the L1014F allele.
A total of 1,688 samples were identified to species level. Of the 871 mosquitoes processed in Yaoundé, 741 (85%) were An. coluzzii and 130 (15%) An. gambiae. In Douala, of the 817 mosquitoes processed 814 (99.6%) were An. coluzzii and three specimens were An. gambiae. No change in mosquito composition was recorded during the study period in the two cities.
Genotyping of intron-1 of voltage sodium channel in mosquito populations emerging from different type of breeding habitats
A total of 43 samples (2 N = 86) collected in 2011 were sequenced to assess the diversity of the intron-1 of the voltage-gated sodium channel in mosquitoes originating from different types of habitats. Due to the small sample size, samples from Douala and Yaoundé were merged and analysed as a unique population. The sample analysed included 15 specimens from unpolluted sites, nine from polluted sites and 19 from cultivated sites. Samples from unpolluted and polluted sites were all An. coluzzii. Samples from cultivated sites were constituted of eleven An. coluzzii and eight An. gambiae. The L1014F allele was detected in 12 of the 43 samples, the L1014S in one of the 43 samples, and 1014 L in 30 samples.
The sequencing of a 1,228 bp region of intron-1 revealed up to four polymorphic sites resulting in five different haplotypes (Table 2). Haplotypes MS1 and MS3 were found to be the most prevalent. Haplotype MS1 was found associated to both resistant alleles (L1014F and L1014S) and the susceptible allele (L1014L). Combinations of kdr alleles and intron-1 haplotypes in different breeding habitats are presented in Table 2. The highest number of haplotypes was recorded in An. coluzzii samples. Estimates of genetic variability in intron-1 showed a haplotype diversity varying from 0.43 for An. gambiae to 0.81 for An. coluzzii samples from unpolluted sites. No departure from neutrality was recorded using the Fs statistics after pooling mosquitoes according to their breeding habitats. The highest Fs value (−0.614) was recorded in samples originating from unpolluted breeding sites and this was not significant (P >0.15).
A fast evolution of pyrethroid resistance prevalence in vector populations of the cities of Douala and Yaoundé was recorded between 2010 and 2013. The situation was in accordance with ongoing reports across the country supporting rapid expansion of pyrethroid and DDT resistance in vector populations [19,20,36,37]. Because pyrethroid resistance was determined after mosquito exposition to reference diagnostic insecticide doses, it was not clear whether the growing prevalence of pyrethroid resistance confers mosquito-elevated capacities to withstand higher insecticide doses or longer exposure time. Hence, data recorded so far, which are the highest ever reported in the country, are consistent with an increased selective pressure on mosquito populations. The rapid drop in mosquito susceptibility to both permethrin and deltamethrin from 2012 suggests the influence of selection induced by the latest massive distribution of ITNs to the population. Since 2003, two large-scale, free distribution campaigns of ITNs were conducted in Cameroon. The latest, conducted in 2011, permitted the distribution of over eight million LLINs and increased substantially bed-net ownership across the country . It is now estimated that over 60% of the population across the country own nets . These initiatives and other factors such as the widespread use of pesticides in agriculture, particularly in urban agriculture, and the frequent use of repellent or insecticide spray might have contributed to the expansion of pyrethroid resistance across the country [19,20]. Although an important number of malaria outbreaks were reported across the country in 2013 , one could not confirm whether it was a consequence of the current levels of pyrethroid resistance in mosquito populations. Yet clinical and entomological studies conducted so far support no reduction in the intensity of malaria transmission in the cities of Yaoundé and Douala following implementation of control measures [40-42]. Because of the large number of confounding factors influencing malaria transmission, the real epidemiological impact of pyrethroid resistance on malaria transmission is still not easy to quantify, justifying the need for more operational studies [6,43].
Mosquitoes were found to be more susceptible to deltamethrin than to permethrin and this probably suggests different resistance mechanisms implicated. High kdr allele frequencies were recorded in vector populations of the two cities. Both kdr alleles L1014F and L1014S were detected, with, however, the L1014F highly prevalent. Although the two alleles were reported for the first time in Cameroon in 2006 , the L1014S allele frequency has not changed much since then. These findings suggest a probable limited role for L1014S allele in pyrethroid resistance in mosquito populations in Cameroon. Nevertheless the presence of L1014F and L1014S alleles in both An. gambiae and An. coluzzii is in agreement with the existence of introgression events between the two species [13,31].
Analysis of sequence variability of intron-1 voltage-gated sodium channel indicated a high variability of this sequence. Of the five haplotypes recorded, MS1 and MS3 were found to be the most prevalent. These haplotypes were found in previous studies . MS1 was predominant in An. coluzzii samples originating from polluted and unpolluted sites and in An. gambiae, whereas MS3 was present only in An. coluzzii samples. The high prevalence of MS3 haplotype in An. coluzzii and its absence in An. gambiae contrasted with previous findings . Moreover, An. gambiae and An. coluzzii samples from cultivated sites presented different haplotypes, supporting the very limited gene flow between the two species even in an area of sympatry and is consistent with earlier reports . No evidence for a recent selective sweep on intron-1 sequence within samples originating from different breeding habitat was detected using Fu’s and Tajima Fs statistics. The highest negative value obtained with the Fu’s Fs statistic was recorded in unpolluted sites but the P-value was not significant. The diversity of intron-1 haplotypes in different types of habitats was consistent with the influence of different evolutionary forces (mutations, gene flow and selection). It is possible that the influence of evolutionary forces in shaping the genetic structure of mosquito populations might be counterbalanced by the ecological fitness of new genetic variants.
It is likely that the current evolution of pyrethroid resistance is maintained by both kdr-based and metabolic-based resistance mechanisms. Studies conducted so far suggested the implication of metabolic resistance in mosquitoes resistant to DDT, permethrin and deltamethrin in different parts of the country [36,37]. Several cytochrome P450 genes, such as CYP6P3, CYP6M2 and glutathion S transferase, were found to be over-expressed in mosquitoes resistant to DDT and pyrethroids . Due to the large number of factors which could select for pyrethroid resistance  and the complexity of resistance mechanisms, it becomes important to take into consideration all these factors to design efficient vector control interventions.
The present study supports rapid evolution of pyrethroid resistance prevalence in vector populations from the cities of Douala and Yaoundé. These changes need further attention in order to ensure the success of future control interventions. Although recent studies suggested the emergence of carbamate resistance  this compound alongside organophosphates could be used as alternative or associated to the use of pyrethroids for vector control. Other control strategies such as the use of larval source management might also be used as complement to treated nets to improve control of pyrethroid resistant mosquitoes particularly in urban settings.
WHO. World Malaria Report 2012. Geneva, Switzerland: World Health Organization; 2012. p. 259.
Bloomquist J. Ion channels as targets for insecticides. Annu Rev Entomol. 1996;41:163–90.
Ranson H, Abdalla H, Badolo A, Guelbeogo W, Kerah-Hinzoumbe C, Yangalbe-Kalnone E, et al. Insecticide resistance in Anopheles gambiae: data from the first year of a multi-country study highlight the extent of the problem. Malar J. 2009;8:299.
Lengeler C. Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004;2:CD000363.
WHO Malaria Policy Advisory Committee and Secretariat. Malaria policy advisory committee to the WHO: conclusions and recommendations of september 2012 meeting. Malar J. 2012;11:424.
Ranson H, Nguessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27:91–8.
Hemingway J, Hawkes N, Mc Carroll L, Ranson H. The molecular basis of insecticides resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34:653–65.
Lynd A, Ranson H, McCall P, Randle N, Black W, Walker E, et al. A simplified high-throughput method for pyrethroid knock-down resistance (kdr) detection in Anopheles gambiae. Malar J. 2005;4:16.
Pinto J, Lynd A, Elissa N, Donnelly M, Costa C, Gentile G, et al. Co-occurrence of East and West African kdr mutations suggests high levels of resistance to pyrethroid insecticides in Anopheles gambiae from Libreville, Gabon. Med Vet Entomol. 2006;20:27–32.
Moreno M, Vicentez J, Cano J, Berzosa P, de Lucio A, Nzambo S, et al. Knockdown resistance mutations (kdr) and insecticide susceptibility to DDT and pyrethroids in Anopheles gambiae from Equatorial Guinea. Trop Med Int Health. 2008;13:430–3.
Martinez-Torrez D, Chandre F, Williamson M, Darriet F, Berge J, Devonshire A, et al. Molecular characterisation of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.
Ranson H, Jensen B, Vulule J, Wang X, Hemingway J, Collins F. 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 Biol. 2000;9:491–7.
Pinto J, Lynd A, Vicente J, Santolamazza F, Randle N, Gentile G, et al. Multiple origins of knockdown resistance mutations in the afrotropical mosquito vector Anopheles gambiae. PLoS ONE. 2007;2:e1243.
Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, Caccone A, et al. Distribution of knock-down resistance mutations in Anopheles gambiae molecular forms in west and west-central Africa. Malar J. 2008;7:74.
Djegbe I, Boussari O, Sidick A, Martin T, Ranson H, Chandre F, et al. Dynamics of insecticide resistance in malaria vectors in Benin: first evidence of the presence of L1014S kdr mutation in Anopheles gambiae from West Africa. Malar J. 2011;10:261.
Brooke B. Kdr: can a single mutation produce and entire insecticide resistance phenotype? Trans R Soc Trop Med Hyg. 2008;103:524–5.
Etang J, Manga L, Chandre F, Guillet P, Fondjo E, Mimpfoundi R, et al. Insecticide susceptibility status of Anopheles gambiae s.l. (Diptera: Culicidae) in the Republic of Cameroon. J Med Entomol. 2003;40:491–7.
Etang J, Fondjo E, Chandre F, Brengues C, Nwane P, Chouaubou M, et al. First report of knockdown mutations in the malaria vector Anopheles gambiae from Cameroon. Am J Trop Med Hyg. 2006;74:795–7.
Chouaibou M, Etang J, Brevault T, Nwane P, Hinzoumbe C, Mimpfoundi R, et al. Dynamics of insecticide resistance in the malaria vector Anopheles gambiae s.l. from an area of extensive cotton cultivation in northern Cameroon. Trop Med Int Health. 2008;13:476–86.
Nwane P, Etang J, Chouaibou M, Toto J, Kerah-Hinzoumbe C, Mimpfoundi R, et al. Trends in DDT and pyrethroid resistance in Anopheles gambiae s.s. populations from urban and agro-industrial settings in southern Cameroon. BMC Infect Dis. 2009;9:163.
Tene Fossog B, Kopya E, Ndo C, Menze-Djantio B, Costantini C, Njiokou F, et al. Water quality and Anopheles gambiae larval tolerance to pyrethroids in the cities of Douala and Yaounde (Cameroon). J Trop Med. 2012;2012:1–10.
Tene Fossog B, Poupardin R, Costantini C, Awono-Ambene H, Wondji C, Ranson H, et al. Resistance to DDT in an urban setting: Common mechanisms implicated in both M and S forms of Anopheles gambiae in the city of Yaoundé Cameroon. PLoS One. 2013;8:e61408.
BUCREP. Troisième recensement général de la population et de l'habitat. Rapport de présentation des résultats définitifs. 2010. http://www.statistics-cameroon.org/downloads/Rapport_de_presentation_3_RGPH.pdf. Accessed 22 April 2010.
WHO. Test procedures for insecticide resistance monitoring in malaria vectors. Bio-efficacy and persistence of insecticides on treated surfaces. WHO/MAL/98,12 Report of the WHO Informal Consultation. Geneva: World Health Organization; 1998.
Abbott WS. A method of computing the effectiveness of an insecticide. J Eco Entomol. 1925;18:265–7.
WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organization; 2013. p. 1–30.
Bass C, Nikou D, Donnelly M, Williamson M, Ranson H, Ball A, et al. Detection of knockdown resistance (kdr) mutations in Anopheles gambiae: A comparison of two new high-throughput assays with existing methods. Malar J. 2007;6:111.
Gillies M, Coetzee M. A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical region). Pub South Afr Inst Med Res. 1987;55:143.
Fanello C, Petrarca V, della Torre A, Santolamazza F, Dolo G, Coulibaly M, et al. The pyrethroid knock-down resistance gene in the Anopheles gambiae complex in Mali and further indication of incipient speciation within Anopheles gambiae s.s. Insect Mol Biol. 2003;12:241–5.
Livak K. Organization and mapping of a sequence on the Drosophila melanogaster X and Y chromosomes that is transcribed during spermatogenesis. Genetics. 1984;107:611–34.
Cornel A, Collins F. PCR of the ribosomal DNA intergenic spacer regions as a methods for identifying mosquitoes in the Anopheles gambiae complex. Methods Mol Biol. 1996;50:321–32.
Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, et al. The kdr mutation occurs in the Mopti form of Anopheles gambiae s.s. through introgression. Insect Mol Biol. 2000;9:451–5.
Rozas J, Sanchez-Delbarrio J, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–7.
Ndjemai H, Patchoke S, Atangana J, Etang J, Simard F, Bilong Bilong C, et al. The distribution of insecticide resistance in Anopheles gambiae s.l populations from Cameroon: an update. Trans R Soc Trop Med Hyg. 2008;103:1127–38.
Antonio-Nkondjio C, Atangana J, Ndo C, Awono-Ambene P, Fondjo E, Fontenille D, et al. Malaria transmission and rice cultivation in Lagdo, northern Cameroon. Trans R Soc Trop Med Hyg. 2008;102:352–9.
Antonio-Nkondjio C, Fossog B, Ndo C, Djantio B, Togouet S, Awono-Ambene P, et al. Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaounde (Cameroon): influence of urban agriculture and pollution. Malar J. 2011;10:154.
Nwane P, Etang J, Chouaibou M, Toto J, Koffi A, Mimpfoundi R, et al. Multiple insecticide resistance mechanisms in Anopheles gambiae s.l. populations from Cameroon, Central Africa. Parasit Vectors. 2013;6:41.
Muller P, Chouaibou M, Pignatelli P, Etang J, Walker ED, Donnelly MJ. 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. 2008;17:1145–55.
UNICEF. Monthly Humanitarian Situation Report UNICEF Cameroon. 2013. http://www.unicef.org/appeals/files/UNICEF_Cameroon_Monthly_Humanitarian_SitRep_Nov_2013.pdf. Accessed 25th november 2013.
Antonio-Nkondjio C, Defo-Talom B, Tagne-Fotso R, Tene-Fossog B, Ndo C, Lehman L, et al. High mosquito burden and malaria transmission in a district of the city of Douala, Cameroon. BMC Infect Dis. 2012;12:275.
Zeukeng F, Tchinda VHM, Bigoga JD, Seumen CHT, Ndzi ES, Abonweh G, et al. Co-infections of malaria and geohelminthiasis in two rural communities of Nkassomo and Vian in the Mfou health district, Cameroon. PLoS Negl Trop Dis. 2014;8:e3236.
Achidi E, Apinjoh T, Anchang-Kimbi J, Mugri R, Ngwai A, Yafi C. Severe and uncomplicated falciparum malaria in children from three regions and three ethnic groups in Cameroon: prospective study. Malar J. 2012;11:215.
Etang J, Vicente J, Nwane P, Chouaibou M, Morlais I, Do Rosario V, et al. Polymorphism of intron-1 in the voltage-gated sodium channel gene of Anopheles gambiae s.s. populations from Cameroon with emphasis on insecticide knockdown resistance mutations. Mol Ecol. 2009;18:3076–86.
Wondji C, Simard F, Fontenille D. Evidence for genetic differentiation between the molecular forms M and S within the Forest chromosomal form of Anopheles gambiae in an area of sympatry. Insect Mol Biol. 2002;11:11–9.
Antonio-Nkondjio C, Youmsi-Goupeyou M, Kopya E, Tene Fossog B, Njiokou F, Costantini C, et al. Exposure to disinfectants (soap and hydrogen peroxide) increases tolerance to permethrin in Anopheles gambiae populations from the city of Yaoundé Cameroon. Malar J. 2014;13:296.
This work was supported by a Wellcome Trust Intermediate Fellowship in Public Health and Tropical Medicine (WTO86423MA) to CAN.
The authors declare that they have no competing interests.
CAN conceived and designed the study protocol; EK, CN, PAA, YP, BMD, CAN, and CSW participated in field and/or laboratory analyses;. CSW, CN, BTF, TT, and PAA critically revised the manuscript; CAN interpreted, analysed data and wrote the paper. All the authors read and approved the final version.