High frequency of kdr L1014F is associated with pyrethroid resistance in Anopheles coluzzii in Sudan savannah of northern Nigeria
BMC Infectious Diseases volume 14, Article number: 441 (2014)
Malaria burden is high in Nigeria, yet information on the major mosquito vectors is lacking especially in the Sudan savannah region of the country. In order to facilitate the design of future insecticide-based control interventions in the region, this study has established the resistance profile of An. gambiae s.l. populations in two northern Nigeria locations and assessed the contribution of target site resistance mutations.
Larval collection was conducted in two localities in Sudan savannah (Bunkure and Auyo) of northern Nigeria between 2009 and 2011, from which resulting adult, female mosquitoes were used for insecticides bioassays with deltamethrin, lambda-cyhalothrin, DDT and malathion. The mosquitoes were identified to species level and molecular forms and then genotyped for the presence of L1014F-kdr, L1014S-kdr and ace-1 R mutations.
WHO bioassays revealed that An. gambiae s.l. from both localities were highly resistant to lambda-cyhalothrin and DDT, but only moderately resistant to deltamethrin. Full susceptibility was observed to malathion. An. gambiae, M form (now An. coluzzii), was predominant over An. arabiensis in Auyo and was more resistant to lambda-cyhalothrin than An. arabiensis. No ‘S’ form (An. gambiae s.s.) was detected. A high frequency of 1014 F mutation (80.1%) was found in An. coluzzii in contrast to An. arabiensis (13.5%). The presence of the 1014 F kdr allele was significantly associated with resistance to lambda-cyhalothrin in An. coluzzii (OR = 9.85; P < 0.001) but not in An. arabiensis. The L1014S-kdr mutation was detected in a single An. arabiensis mosquito while no ace-1 R mutation was found in any of the mosquitoes analysed.
The predominance of An. coluzzii and its resistance profile to main insecticides described in this study can guide the implementation of appropriate vector control interventions in this region of Nigeria where such information was previously lacking.
Malaria is one of the most widespread infectious diseases of our time, taking the lives of more than 650,000 people a year, most of them in sub-Saharan Africa and under the age of 5 ; 1300 young lives are estimated to be lost to malaria every day . Malaria is holoendemic in Nigeria, accounting for 25% of infant mortality and 30% of childhood mortality  and transmission of malaria is geographically specific . Malaria accounts for 60% of outpatient visits and 30% of hospitalizations among children under five years of age in Nigeria. In 2010, Nigeria with a population of about 150 million and reporting more deaths due to malaria than any country in the world, became the seventeenth President’s Malaria Initiative country . Nigeria and Democratic Republic of Congo account for more than 40% of the total global estimated malaria deaths .
Information about distribution of Anopheles vectors and pattern of resistance to the major insecticides from Nigeria especially northern region is grossly lacking. A recent analysis of data on the distribution of Anopheles mosquitoes across Nigeria from 1900–2010 reported the predominance of mosquitoes from the An. gambiae s.l. complex (65.2%) followed by those from the An. funestus group (17.3%) . However, most of the reports came from studies in southern Nigerian with only few studies from the north. Presently, vector control strategies involve the use of pyrethroid-based long-lasting insecticidal treated nets (LLINs) and indoor residual spray (IRS) although bendiocarb and dichlorodiphenyltrichloroethane (DDT) are also used in some areas for IRS . This widespread use of a single class of insecticide could result in development of more insecticide resistance in the mosquito vectors and lead to a major public health problem given the limited availability of alternative insecticides .
As of May 2011 a total of 35.6 million LLINs had been distributed across 22 states of Nigeria, with a balance of 27.3 million to complete the remaining 15 states. Results from the malaria indicator survey conducted from October to December 2010, showed significant increases in LLIN ownership and use as compared to the 2008 Nigeria’s demographic health survey . The average percentage for ownership of one LLIN in 2008 was 8% and increased to 42% in 2010, with rural ownership higher (45%) than urban (33%). This together with the World Bank supported IRS and insecticide treated nets (ITN) programs in several states in Nigeria  may add to the selective pressure on malaria vectors to develop more resistance against pyrethroids.
The two major causes of insecticide resistance in malaria vectors are alterations in the target sites and an increase in the rate of insecticide metabolism . The former is characterized by knock-down resistance (kdr) to pyrethroids and DDT caused by mutation in the voltage-gated sodium channel (VGSC) and well characterized  with the L1014F-kdr mutation predominantly found in An. gambiae from West Africa and the L1014S-kdr mutation predominantly found in East Africa . Another variant of target site mechanism of insecticide resistance in An. gambiae is a single amino acid substitution of glycine to serine at position 119 in the catalytic domain of the acetylcholinesterase (AChE) gene which confers resistance to both organophosphates and carbamates .
Pyrethroid knockdown resistance in An. gambiae s.s. from Nigeria was first reported in 2001  and most of the recent studies on dominant vector species (DVS) composition and resistance profiles were carried out mostly in the south-western  and north-central Guinea Savannah  of Nigeria. Unfortunately, apart from the Garki Project carried out long time ago in 1960–1970  very little has been done in the recent decades on analysis of the DVS compositions and insecticide resistance profiles of the major malaria vector from Sudan and Sahelian savannah of northern Nigeria . There is thus an overwhelming need for a continued entomological survey of the malaria vector species and dissecting their resistance profile in this and other regions of Nigeria.
Here we characterised two populations of An. gambiae s.l. from northern Nigeria, established their species compositions, resistance profiles to insecticides, particularly pyrethroids and DDT and genotyped the main resistance markers.
Auyo (12°21′N, 9°59′E) located northeast of Dutse, the capital of Jigawa State (Figure 1) is a town known for its history of irrigation activities in which rice and other vegetables are produced. Jigawa state is situated within the Sudan savannah vegetation zone, but there are traces of Guinea savannah in the southern part of the state . Larval collection was conducted in Auyo twice: in rainy season, June 2009 (for most bioassay tests) and during the harmattan of February 2011 (for malathion bioassay only).
Bunkure (11°42′N, 8°33′E), is a town located approximately 88.5 km east of Kano city. Kano State is located within the Sudan savannah zone of West Africa about 840 kilometers from the edge of the Sahara desert . The vegetation of Kano State is semi-arid savannah sandwiched by the Sahel savannah in the north and the Guinea savannah in the south. The state has the largest irrigation projects in Nigeria, with six irrigation projects and more than twenty earth dams. Alongside the importance of these dams to the development of agriculture and provision of food comes the health implication of providing suitable breeding sites for vectors of diseases. Rice paddies in particular have been established and increase the risk of malaria by providing suitable sites for vector development. The locality of Bunkure has a large expanse of irrigable lands where rice is the major crop grown, although other crops such as green vegetables and tomatoes are grown for subsistence. Larval collections were done at the height of the rainy season in August 2010.
Both in Bunkure and Auyo farmers use petrol pumps and a myriad of pesticides to protect crops. These pesticides include organophosphates, carbamates, pyrethroids and organochlorine (including DDT and dieldrin). Reliable data on the impacts of agricultural pesticides in the wetland environments of northern Nigeria is lacking .
Mosquito collection and rearing
Mosquito larvae were collected from several rice paddies within the vicinity of Bunkure town, Kano and twice in Auyo town, Jigawa using classical dipping method as described by . The third and fourth instar larvae identified as belonging to the Anopheles gambiae complex using morphological keys  were transported to the insectary at Aminu Kano Teaching Hospital, Kano. The larvae were maintained under standard insectary condition (25-28°C and ~70-80% humidity, with a 12 hr day/night cycle)  and supplied with baby fish food daily. The adults that emerged were fed 10% sucrose solution and randomly mixed for subsequent experiments.
Mosquito species identification
Mosquitoes were identified as belonging to the An. gambiae s.l. complex using the morphological keys of . Genomic DNA was extracted from mosquitoes obtained from Auyo, Jigawa state and which survived exposure to lambda-cyhalothrin, using the LIVAK method . The species identity of the An. gambiae s.l. mosquitoes and the molecular forms of all the An. gambiae s.s. were determined according to the SINE PCR method .
Insecticide susceptibility bioassays
Insecticide susceptibility assays were carried out on 2–3 days old adult female An. gambiae mosquitoes following the protocol outlined by WHO . Three to fifteen replicates of around 20–25 mosquitoes per tube were exposed to insecticide impregnated papers for 1 hr and then transferred to a holding tube and supplied with 10% sucrose. More replicates were used for lambda-cyhalothrin (twelve and fifteen replicates respectively for Bunkure and Auyo) with the aim of covering a large number of breeding sites in these locations. 20–25 mosquitoes were also used as a control with no exposure to any insecticide. Knockdown was recorded at intervals of 10 minutes up to an hour post-exposure to insecticides. Mortality was recorded 24 hours post-exposure to 0.05% deltamethrin, 0.05% lambda-cyhalothrin, 4% DDT and 5% malathion. Susceptible individuals are defined as individuals that did not survive a discriminating dose of the particular insecticide used. Mosquitoes from Bunkure were not tested with DDT and malathion due to these insecticide papers not being available when mosquito collection was carried out in this location.
Genotyping of target site mutation
Presence of kdr and the ace-1 R mutations were investigated in An. coluzzii and An. arabiensis mosquitoes collected from Auyo and which survived exposure to lambda-cyhalothrin, using the pyrosequencing method as previously described [28–30]. Briefly, three sequence specific primers (Table 1) each for kdr and ace-1 mutations designed using the software provided by Pyrosequencing AB (http://www.pyrosequencing.com) were used for amplification and genotyping of both mutations. Target DNA fragments for the kdr and ace-1 were first PCR-amplified in a reaction containing 10 pmol each of forward and biotinylated reverse primer, 1X HotStarTaq buffer, 0.2 mM dNTPs, 1.5 mM MgCl2, 1U HotStarTaq (Qiagen) and 10 ng genomic DNA. The amplification was carried out using the following conditions: 1 cycle at 95°C for 5 min; 50 cycles at 94°C for 20s, 57°C for 30s and extension at 72°C for 20s; and finally 1 cycle at 72°C for 5 min. Pyrosequencing was carried out according to the manufacturer’s instructions using the PSQ 96 SNP Reagent kit (Qiagen) and the genotype determined using the SNP software (Biotage AB).
The correlation between the L1014F kdr genotypes and lambda-cyhalothrin resistance phenotypes was assessed by estimating the odds ratios (OR) and the statistical significance based on the Fisher exact probability test.
Insecticide susceptibility bioassays
Insecticides bioassay conducted using the mixed F0 adults from Auyo, revealed resistance to DDT with only 12% knockdown after one hour exposure (Figure 2A) and mortality of 44.6% after 24 hours (Figure 3A). These mosquitoes were also resistant to lambda-cyhalothrin and deltamethrin (representative type II pyrethroids) with slow knockdown of about 30% after an hour of exposure and mortality of 52.8% and 78.4% respectively 24 hours post-exposure. However, full susceptibility was obtained with malathion (an organophosphate) with 100% knockdown within the first 40 minutes after exposure and 100% mortality after 24 hours.Mosquitoes from Bunkure were also more resistant to lambda-cyhalothrin with knockdown of only 32.6% (Figure 2B) in the first hour of exposure and mortality of 52.1% after 24 hours (Figure 3B). However, this population was found not to be resistant to deltamethrin, for 71.7% of the mosquitoes were knocked down in the first hour of exposure and 97.2% died 24 hours after exposure.
Species and molecular forms identification
A total of one hundred and twenty one (121) females from Auyo that survived exposure to lambda-cyhalothrin and twenty six (26) dead after exposure were analysed using SINE PCR for species identification (Table 2). Approximately 77% of the mosquitoes were An. coluzzii (formerly M form), 22% were An. arabiensis and a single hybrid mosquito (<1.0%) was also found. No An. gambiae s.s. (formerly S form) were detected. An. coluzzii was predominant among the mosquitoes resistant to lambda-cyhalothrin (86.8%) while An. arabiensis was predominant in the susceptible mosquitoes (69.2%). Overall, 92.9% (105/113) of An. coluzzii were resistant to lambda-cyhalothrin whereas only 45.4% (15/33) An. arabiensis were resistant. Molecular analyses were not conducted with the rest of the mosquitoes used in this research for logistic reasons.
Genotyping of target-site mutations
Knockdown resistance (kdr) mutations
A total of one hundred and fifteen (115) mosquitoes from Auyo which were exposed to lambda-cyhalothrin and identified to molecular level were successfully pyrosequenced simultaneously for both West-kdr (L1014F) and East-kdr (L1014S) mutations. This included ninety-five mosquitoes alive after exposure (81 An. coluzzii, 13 An. arabiensis and one hybrid) and twenty dead (7 An. coluzzii and 13 An. arabiensis). The pyrosequencer unambiguously scored the homozygote resistant (T/T), heterozygote resistant (A/T) and homozygote susceptible (A/A) genotypes for the L1014F-kdr as well as the L1014S-kdr genotypes.
Out of 81 genotyped An. coluzzii individuals alive after exposure, 74.1% (60) were homozygote (T/T) for L1014F-kdr, 19.7% (16) were heterozygotes, while only 6.2% (5) individuals were homozygote susceptible (A/A) (Figure 4A). Despite a low sample size of only seven An. coluzzii among the dead mosquitoes, a significant difference in the genotype distribution could be observed compared with the resistant mosquitoes. Consequently, a significant correlation was observed between the 1014 F allele and the resistance to lambda-cyhalothrin [OR = 9.4 (4.7-48.4); P < 0.001] (Table 2).
In contrast to An. coluzzii, most An. arabiensis resistant individuals (69.2%) were homozygote for the susceptible allele (A/A) while 23.1% (3) individuals were heterozygotes (A/T) and only 7.7% (1) homozygote for the resistant allele (T/T) (Figure 4B). Similar frequencies of the respective genotypes were observed in the dead An. arabiensis. Consequently, no significant correlation is observed between the 1014 F allele and the resistance to lambda-cyhalothrin in An. arabiensis [OR = 2.85 (0.5-16.4); P = 0.4] (Table 2).
In total, a frequency of 80.1% was detected in An. coluzzii for the 1014 F kdr resistant allele in contrast with only 13.5% in An. arabiensis. One An. arabiensis individual was found to be heterozygote for the L1014S-kdr mutation and homozygote susceptible (A/A) for L1014F-kdr. Also, the single hybrid (An. arabiensis/An. gambiae) individual genotype turned out to be heterozygote (A/T) for L1014F-kdr mutation.
Acetylcholinesterase (AChE) mutation
All the mosquitoes (100) used for bioassay with malathion were genotyped for G119S ace-1 R mutation. The pyrosequencer unambiguously scored all the samples with wild type G/G combination. No ace-1 R mutation was found in any of the mosquitoes consistent with the full susceptibility observed with malathion (Figure 3A).
This study has investigated the dynamics of species composition in the An. gambiae complex and their susceptibility status to main insecticides in northern Nigeria. Larval collection was conducted from several rice paddies in both rainy and dry season, but larger numbers of An. coluzzii were obtained than An. arabiensis in the collection during the rainy season of 2009 in Auyo, Jigawa state.
This supports previous observations that An. gambiae s.s. could be predominant in the Sudan savannah ecological zone, compared to An. arabiensis that was spread across Sudan, Sahel and northern Guinea savannah ecological zones [31, 32]. However, further collections are needed to fully establish the malaria vector species distribution in this region of Nigeria.
Mosquitoes from both Auyo and Bunkure localities used in this research showed a high resistance to lambda-cyhalothrin. However, only marginal resistance in mosquitoes from both sites was observed against deltamethrin. Previous studies have shown that same mosquito populations can be resistant and also susceptible to different insecticides from the same family. For example, a kdr-free population of An. arabiensis from Chad was shown to be highly resistant to permethrin but fully susceptible to deltamethrin . One cytochrome P450 (CYP6P4) was discovered to be up-regulated in this population and shown to preferentially metabolise permethrin but not deltamethrin (S. Ibrahim, Personal Communication). Marginal resistance observed with deltamethrin indicates that kdr alone might not be sufficient to confer resistance to pyrethroid; it is likely that metabolic resistance is also involved. Djouaka and colleagues  have documented the involvement of cytochrome P450s irrespective of presence or absence of kdr linked with resistance in ‘M’ form of An. gambiae from southern Benin and south-western Nigeria. Also, Awolola and colleagues  have reported the presence of kdr in susceptible species of An. gambiae from south-western Nigeria and the absence of the kdr alleles in resistant mosquitoes of the molecular ‘M’ form and its presence in ‘S’ form . Thus, metabolic resistance might play a major role in the pyrethroid resistance observed in the localities sampled, especially with respect to An. arabiensis which has low frequency of L1014F-kdr genotype.
Pyrethroid resistance observed in these An. coluzzii populations is similar to previously reported cases of resistance to permethrin and deltamethrin in north-central and south-western Nigeria [17, 34]. In addition, the DDT resistance observed in mosquitoes from Auyo is in keeping with previous studies elsewhere in Nigeria  and suggests a similar mechanism of resistance against pyrethroid and organochlorine insecticides.
Insecticide resistance was first reported in the Sudan savannah of western provinces of Sokoto, northern Nigeria back in 1958  and specifically, resistance to dieldrin was first reported in An. gambiae from the same ecological zone in 1959 . However, to our knowledge this is the first study that investigates insecticide resistance in An. gambiae s.l. from Sudan savannah of northern Nigeria after several decades, particularly against the pyrethroid insecticides. It is also the first of its kind to report the presence of East form of kdr (L1014S mutation) in An. arabiensis from Nigeria though the mutation has been reported in low frequency in An. arabiensis, and in the ‘M’ and ‘S’ molecular forms of An. gambiae s.s. from Burkina-Faso , An. arabiensis from Benin , as well as An. arabiensis from Sudan . The low frequency of the East-kdr in these regions indicates that its selection is very recent or it has just recently migrated into these regions. The high frequency of West-kdr (L1014F) mutation observed in the An. coluzzii from Auyo is in contrast to observations from south-western Nigeria  where the kdr is high in frequency mainly in ‘S’ form and low in frequency in ‘M’ form of An. gambiae. The absence of ‘S’ form of Anopheles gambiae from the field collections conducted in Auyo could be explained by the observation previously made that the ‘M’ form (Mopti chromosomal form) is predominant in this type of ecological setting with a more permanent breeding site from the irrigation system .
The pyrethroid and DDT resistance observed in the An. coluzzii from this region suggests that this population is under local selective pressure. Indeed, increased usage of insecticides for agricultural purposes and/or widespread use of LLINs in the region could explain the high frequency of kdr mutation as recently observed in another ‘M’ form population of An. gambiae in a similar agricultural setting within Burkina Faso .
The few An. arabiensis genotyped exhibited low frequency of L104F-kdr. Unfortunately, in addition to the low number of An. arabiensis obtained in collection from Auyo in 2009, we don’t know the resistance level of An. arabiensis separately, as the bioassay results are for the two species combined. It is possible that An. arabiensis is as resistant as An. coluzzii but with less involvement of kdr as seen in Chad where a kdr-free resistant An. arabiensis population was reported with resistance driven mainly by P450s. Larger collection and profiling of resistance pattern of both An. gambiae s.s. and An. arabiensis from these sites separately can help to further elucidate this. There is also an overwhelming need to explore the involvement of metabolic resistance in Anopheline mosquitoes in this region using synergist assay and molecular tools.
The mosquitoes obtained from Auyo were fully susceptible to malathion and the susceptibility pattern was confirmed by the absence of the ace-1 R mutation in all the mosquitoes analysed. Susceptibility to malathion has already been established in several other An. gambiae populations across Africa such as Burkina-Faso  and Cameroon  and also in An. funestus another major malaria vector [43, 44]. Malathion could therefore be used as an alternative to pyrethroids and DDT in IRS campaigns in Sudan savannah of northern Nigeria.
Control of malaria vector cannot be achieved without containment of resistance to insecticides and implementation of robust resistance management strategy hinges on the knowledge of the distribution and composition of the major malaria vectors and their resistance profiles. In this research we have established the temporal and spatial distribution of the two major malaria vectors An. coluzzii and An. arabiensis from Sudan savannah of northern Nigeria. We also described the pyrethroid, DDT and malathion resistance profiles in these mosquitoes as well as the high frequency of L1014F-kdr mutation in the An. coluzzii and presence of L1014S-kdr in low frequency in An. arabiensis. The ace-1 R mutation that confers resistance to carbamates and organophosphates is absent in the mosquitoes. This type of information could guide the Nigerian Malaria Control Programme in the choice of appropriate resistance management strategy to implement vis-à-vis different ecological zones in Nigeria.
World Health Organization: World Malaria Report. WHO Global Malaria Programme. 2011, Geneva, Switzerland: WHO
World Health Organization: World Malaria Report. WHO Global Malaria Programme. 2013, Geneva, Switzerland: WHO
Onwujekwe O, Uguru N, Etiaba E, Chikezie I, Uzochukwu B, Adjagba A: The economic burden of malaria on households and the health system in Enugu State southeast Nigeria. PLoS One. 2013, 8 (11): e78362-
Gallup JL, Sachs JD: The economic burden of malaria. Am J Trop Med Hyg. 2001, 64 (1–2 Suppl): 85-96.
NIGERIA FY 2011 Malaria Operational Plan. [http://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy11/nigeria_mop-fy11.pdf?sfvrsn=6]
World Health Organization: World Malaria Report. WHO Global Malaria Programme. 2012, Geneva, Switzerland: WHO
Okorie PN, McKenzie FE, Ademowo OG, Bockarie M, Kelly-Hope L: Nigeria Anopheles vector database: an overview of 100 years’ research. PLoS One. 2011, 6 (12): e28347-
Ranson H, N’Guessan 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 (2): 91-98.
Prevention and management of insecticide resistance in vectors of public health importance. [http://www.irac-online.org/documents/moa-brochure/?ext=pdf]
Seventh annual report to congress. [http://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy13/nigeria_mop_fy13.pdf?sfvrsn=8]
Hemingway J, Ranson H: Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000, 45: 371-391.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, Guillet P, Pasteur N, Pauron D: Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998, 7 (2): 179-184.
Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH: 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 (5): 491-497.
Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, Raymond M: The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004, 13 (1): 1-7.
Awolola TS, Brooke BD, Hunt RH, Coetze M: Resistance of the malaria vector Anopheles gambiae s.s. to pyrethroid insecticides, in south-western Nigeria. Ann Trop Med Parasitol. 2002, 96 (8): 849-852.
Oduola AO, Idowu ET, Oyebola MK, Adeogun AO, Olojede JB, Otubanjo OA, Awolola TS: Evidence of carbamate resistance in urban populations of Anopheles gambiae s.s. mosquitoes resistant to DDT and deltamethrin insecticides in Lagos, South-Western Nigeria. Parasites Vectors. 2012, 5: 116-
Adeogun A, Olojede J, Oduola A, Awolola T: Village-scale evaluation of PermaNet 3.0: an enhanced efficacy combination long-lasting insecticidal net against resistant populations of Anopheles gambiae ss. Malar Chemother Control Elimination. 2012, 1: 9-
Molineaux L, Gramiccia G: The Garki Project: Research On The Epidemiology And Control Of Malaria In The Sudan Savanna Of West Africa. 1980, Geneva, Switzerland: WHO
Jigawa State. A New World. [http://www.jigawastate.gov.ng/contentpage.php?id=82]
Ado-Kurawa I: Efforts of the Shekarau Administration in Harnessing Resources for Social and Economic Development. 2006, Jaji: Research and Documentation Department
Kimmage K, Adams WM: Small-scale farmer-managed irrigation in northern Nigeria. Geoforum. 1990, 21 (4): 435-443.
Robert V, Le Goff G, Ariey F, Duchemin JB: A possible alternative method for collecting mosquito larvae in rice fields. Malar J. 2002, 1: 4-
Gillies M, Coetzee M: A supplement to anophelinae of Africa south of Sahara (Afro-tropical region). The South African Institute for Medical Research. 1987, 55: 96-110.
Das S, Garver L, Dimopoulos G: Protocol for mosquito rearing (An. gambiae). J Vis Exp. 2007, 221-doi:10.3791/221., Issue 5
Livak KJ: Organization and mapping of a sequence on the Drosophila melanogaster X and Y chromosomes that is transcribed during spermatogenesis. Genetics. 1984, 107 (4): 611-634.
Santolamazza F, Mancini E, Simard F, Qi Y, Tu Z, Della Torre A: Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar J. 2008, 7: 163-
World Health Organization: Test Procedures for Insecticide Resistance Monitoring in Malaria Vectors, Bio-efficacy and Persistence of Insecticides on Treated Surfaces: Report of the WHO Informal Consultation. 1998, Geneva, Switzerland: WHO, 28-30.
Wondji CS, Morgan J, Coetzee M, Hunt RH, Steen K, Black WC, Hemingway J, Ranson H: Mapping a quantitative trait locus (QTL) conferring pyrethroid resistance in the African malaria vector Anopheles funestus. BMC Genomics. 2007, 8: 34-
Witzig C, Parry M, Morgan JC, Irving H, Steven A, Cuamba N, Kerah-Hinzoumbe C, Ranson H, Wondji CS: Genetic mapping identifies a major locus spanning P450 clusters associated with pyrethroid resistance in kdr-free Anopheles arabiensis from Chad. Heredity. 2013, 110 (4): 389-397.
Kwiatkowska RM, Platt N, Poupardin R, Irving H, Dabire RK, Mitchell S, Jones CM, Diabate A, Ranson H, Wondji CS: Dissecting the mechanisms responsible for the multiple insecticide resistance phenotype in Anopheles gambiae s.s., M form, from Vallee du Kou, Burkina Faso. Gene. 2013, 519 (1): 98-106.
Onyabe DY, Conn JE: The distribution of two major malaria vectors, Anopheles gambiae and Anopheles arabiensis, in Nigeria. Mem Inst Oswaldo Cruz. 2001, 96 (8): 1081-1084.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 73 (5): 483-497.
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 Genomics. 2008, 9: 538-
Awolola TS, Oduola AO, Oyewole IO, Obansa JB, Amajoh CN, Koekemoer LL, Coetzee M: Dynamics of knockdown pyrethroid insecticide resistance alleles in a field population of Anopheles gambiae s.s. in southwestern Nigeria. J Vector Borne Dis. 2007, 44 (3): 181-188.
Awolola TS, Oduola OA, Strode C, Koekemoer LL, Brooke B, Ranson H: Evidence of multiple pyrethroid resistance mechanisms in the malaria vector Anopheles gambiae sensu stricto from Nigeria. Trans R Soc Trop Med Hyg. 2009, 103 (11): 1139-1145.
Armstrong JA, Ramsdale CD, Ramakrishna V: Insecticide resistance in Anopheles gambiae Giles in Western Sokoto, Northern Nigeria. Ann Trop Med Parasitol. 1958, 52 (3): 247-256.
Ramakrishna V, Elliott R: Insecticide resistance in Anopheles gambiae in Sokoto province. Trans R Soc Trop Med Hyg. 1959, 53 (1): 102-109.
Namountougou M, Diabate A, Etang J, Bass C, Sawadogo SP, Gnankinie O, Baldet T, Martin T, Chandre F, Simard F, Dabire RK: First report of the L1014S kdr mutation in wild populations of Anopheles gambiae M and S molecular forms in Burkina Faso (West Africa). Acta Trop. 2013, 125 (2): 123-127.
Djegbe I, Boussari O, Sidick A, Martin T, Ranson H, Chandre F, Akogbeto M, Corbel V: 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-
Himeidan YE, Chen H, Chandre F, Donnelly MJ, Yan G: Short report: permethrin and DDT resistance in the malaria vector Anopheles arabiensis from eastern Sudan. AmJTrop Med Hyg. 2007, 77 (6): 1066-1068.
Coluzzi M: Heterogeneities of the malaria vectorial system in tropical Africa and their significance in malaria epidemiology and control. Bull World Health Organ. 1984, 62 (Suppl): 107-113.
Antonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, Awono-Ambene P, Costantini C, Wondji CS, Ranson H: 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-
Djouaka R, Irving H, Tukur Z, Wondji CS: Exploring Mechanisms of Multiple Insecticide Resistance in a Population of the Malaria Vector Anopheles funestus in Benin. PLoS One. 2011, 6 (11): e27760-
Morgan JC, Irving H, Okedi LM, Steven A, Wondji CS: Pyrethroid resistance in an Anopheles funestus population from Uganda. PLoS One. 2010, 5 (7): e11872-
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/441/prepub
This work was supported by a Wellcome Trust RCD Fellowship (083515/Z/07/Z) to CSW. Molecular work was conducted in The Liverpool School of Tropical Medicine, United Kingdom. We thank the authorities of the Ministry of Health and all the staff of the Malaria Control Unit, Jigawa State for permission and support during this project. Financial support for field work was provided by the World Bank under the Malaria Control Booster Initiative in collaboration with Nigeria’s National Malaria Control Programme. We would also like to thank Mr Ado A.K. and Dr Kabiru A.I. for their cooperation and invaluable support.
The authors declare no competing interest.
YM and CW designed the study. YM, ZT, SS and HI conducted the experiments. SS and CW interpreted the data and wrote the manuscript. YM, HI and ZT reviewed the manuscript. All authors read and approved the final draft of the manuscript.
About this article
Cite this article
Ibrahim, S.S., Manu, Y.A., Tukur, Z. et al. High frequency of kdr L1014F is associated with pyrethroid resistance in Anopheles coluzzii in Sudan savannah of northern Nigeria. BMC Infect Dis 14, 441 (2014). https://doi.org/10.1186/1471-2334-14-441