- Research article
- Open Access
Mapping a Quantitative Trait Locus (QTL) conferring pyrethroid resistance in the African malaria vector Anopheles funestus
© Wondji et al; licensee BioMed Central Ltd. 2007
- Received: 31 August 2006
- Accepted: 29 January 2007
- Published: 29 January 2007
Pyrethroid resistance in Anopheles funestus populations has led to an increase in malaria transmission in southern Africa. Resistance has been attributed to elevated activities of cytochrome P450s but the molecular basis underlying this metabolic resistance is unknown. Microsatellite and SNP markers were used to construct a linkage map and to detect a quantitative trait locus (QTL) associated with pyrethroid resistance in the FUMOZ-R strain of An. funestus from Mozambique.
By genotyping 349 F2 individuals from 11 independent families, a single major QTL, rp1, at the telomeric end of chromosome 2R was identified. The rp1 QTL appears to present a major effect since it accounts for more than 60% of the variance in susceptibility to permethrin. This QTL has a strong additive genetic effect with respect to susceptibility. Candidate genes associated with pyrethroid resistance in other species were physically mapped to An. funestus polytene chromosomes. This showed that rp1 is genetically linked to a cluster of CYP6 cytochrome P450 genes located on division 9 of chromosome 2R and confirmed earlier reports that pyrethroid resistance in this strain is not associated with target site mutations (knockdown resistance).
We hypothesize that one or more of these CYP6 P450s clustered on chromosome 2R confers pyrethroid resistance in the FUMOZ-R strain of An. funestus.
- Quantitative Trait Locus
- Quantitative Trait Locus Mapping
- Lambda Cyhalothrin
- Pyrethroid Resistance
- Quantitative Trait Locus Position
Anopheles funestus is widely distributed in Africa from south of the Sahara desert to northern South Africa . The highly anthropophilic and endophilic behaviors of this vector make it an efficient vector of malaria, and in many places, parasite infection rates of An. funestus exceeds those of An. gambiae. In South Africa, indoor house spraying with the insecticide DDT had successfully eradicated malaria transmission by An. funestus. However, after DDT was replaced by pyrethroids in 1996, resistance to the latter insecticides developed rapidly, and there was a dramatic increase in malaria transmission in Kwazulu/Natal Province, accompanied by the finding of An. funestus resting inside insecticide sprayed houses . Subsequent surveys of southern Mozambican An. funestus also showed high levels of resistance to a variety of pyrethroids (permethrin, deltamethrin and lambda-cyhalothrin) and the carbamate propuxur , suggesting that the pyrethroid resistant An. funestus in South Africa had arrived from neighbouring Mozambique where An. funestus has remained abundant  and transmits malaria perennially .
The two major causes of insecticide resistance are alterations in the target sites and increases in the rate of insecticide metabolism. Three enzyme families, the cytochrome P450s, the esterases and glutathione S-transferases (GSTs), are primarily responsible for metabolic resistance . Both field populations and a laboratory-selected colony of An. funestus from Mozambique show elevated levels of cytochrome P450s in pyrethroid resistance populations . Cytochrome P450-based pyrethroid resistance also occurs in East and West African An. gambiae [7, 8] although in West Africa the P450 mechanism is at a much lower frequency than the altered target site [knockdown resistance (kdr)-type] resistance mechanism . The appearance of a P450-based pyrethroid resistance in southern African An. funestus had an immediate and significant effect on malaria transmission in South Africa prompting a switch back to DDT in South Africa and the introduction of the carbamate, bendiocarb in Mozambique [2, 3]. However, resistance to these alternative insecticides has already developed in An. arabiensis in South Africa  and hence effective resistance management strategies are imperative in this region. Identification and characterization of the specific P450 genes involved in pyrethroid metabolism in resistant insects is an important step towards this goal.
The cytochrome P450 enzyme family is very extensive in insects [10, 11]. Only a small subset of these enzymes is likely to be involved in insecticide metabolism. Partial sequences of approximately thirty P450 genes from three families (CYP 4. 6 and 9) have already been identified in An. funestus , and here we report the physical mapping of a subset of these genes to the An. funestus polytene chromosomes. We also describe the first application of quantitative trait locus (QTL) mapping in this malaria vector. QTL mapping techniques have been successfully used in other mosquito vectors to examine complex phenotypes including vectorial competence [13, 14] and insecticide resistance [15, 16]. The QTL approach has never been used in An. funestus mainly because of the difficulties associated with colonization of this species and the lack of available genetic markers. However, the recent establishment of two laboratory colonies of An. funestus  and the construction of a preliminary genetic map  have paved the way for the current study. In this paper we report the genetic and physical map position of a QTL associated with pyrethroid resistance in An. funestus.
Markers with significant association with pyrethroid resistance for each family
No significant marker
No significant marker
For family 11, statistically significant associations were found between genotypes at three loci (7P6P4, 3P6P3 and AFND6) and the resistance trait. These loci showed a similar pattern of association with pyrethroid resistance as in family 4. A high mortality rate was observed in "aa" homozygous F2 individuals (100, 60 and 100% respectively for 7P6P4, 3P6P3 and AFND6) while mortality was low in "bb" individuals (19 and 39%, respectively for 7P6P4 and AF6).
For family 6, only two markers (3BU82 and 6BU40) were statistically associated with pyrethroid resistance. Again, the mortality rate was higher in F2 "aa" homozygotes (65 and 67% respectively for 3BU82 and 6BU40) than in F2 "bb" homozygotes (25 and 35% respectively for 3BU82 and 6BU40).
We next used interval mapping (IM), composite-interval mapping and multiple-interval mapping to predict the location of this QTL. Three families with a sample size of more than 40 F2 progeny (families 4, 6 and 11) were independently analyzed for the presence of QTL associated with pyrethroid resistance.
Multiple interval mapping estimates of QTL position and associated genetic, environmental, phenotypic, additive and dominance effects associated with pyrethroid resistance in An. funestus.
Genetic distance (cM)
A, -5.589 D, -0.660
A, 2.1649 D, 3.6980
A, 1.0312 D, 3.2745
A, -4.825 D, +2.762
A, -2.4513 D, -0.7513
In families 6 and 10 (data not shown in Table 2 for family 10 because of the low sample size), the markers 7P6P4 and 3P6P3 were not informative and only low LOD scores were observed at markers at the end of 2R chromosome. However when using multiple-interval mapping, we detected the same QTL as that observed in families 4 and 11. The LOD score for this QTL in family 6 and 10 was 1.9 and 1.6 respectively with composite-interval mapping, which is considerably smaller than that observed for families 4 and 11. This is probably due to the small sample size (family 10) and the lack of informative markers in this region of the genome. In families 4 and 11, as in family 6 (Table 2), additive effects accounted for the majority of the genetic variance.
Physical mapping of rp1 QTL
Chromosomal locations of P450, Acetylcholinesterase (AChE) and knockdown resistance (kdr) genes in An. funestus and their orthologs in An. gambiae.
Location in An. funestus
Location in An. gambiae
Five CYP6 P450s physically mapped to chromosome arm 2R, within the region of rp1 QTL (Figure 6). These included CYP6P4 and CYP6P3 from which two SNP loci, (7P6P4 and 3P6P3) were identified. In agreement with the physical mapping data, we found that these SNPs were closely linked to rp1 in all families for which these markers were informative.
The kdr locus is located on chromosome 3R, division 36, and is not linked to pyrethroid resistance in the FUMOZ-R strain of An. funestus. Again this agrees with the chromosomal location determined for this gene in A. gambiae (3R in An. funestus is equivalent to 2L in A. gambiae). No amino acid substitutions were found in the amplified 1342 bp fragment of the voltage gated sodium channel gene in the FOMUZ-R strain compared to the field samples of An. funestus and the classical kdr leucine substitution at position 1014, associated with pyrethroid resistance in many insect species, was not found in An. funestus.
The fragment of the AChE gene sequenced encompasses the common ace-1 mutation site associated with insensitive AChE in An. gambiae . We did not detect this mutation in the An. funestus AChE.
To detect QTL associated with resistance to pyrethroids in An. funestus, we constructed a linkage map of this species using combined data from 11 families generated from reciprocal crosses between a susceptible and permethrin resistant strain. The average resolution of 2.7 cM/marker is similar to the resolution of the published An. gambiae microsatellite map . The present An. funestus genetic map uniformly covers the entire genome of An. funestus with a similar marker density for chromosomes 2 and 3. This genetic map represents a significant improvement from the previous published map of .
The rp1 QTL that affects An. funestus susceptibility to pyrethroids was consistently identified in both reciprocal crosses and in the majority of the families (7 out of 11), demonstrating the importance of this QTL. Generally the LOD scores associated with the rp1 QTL were well above the thresholds determined by permutation analysis using both interval mapping and composite-interval mapping. The rp1 QTL was most closely associated with markers 7P6P4 and 3P6P3 on chromosome 2R regardless of the analytical method employed (IM, CIM and MIM). Such correlation between different methods and consistency between different families represent a strong indication that the QTL identified in this study reflects the detection of a true locus involved in resistance to pyrethroids in An. funestus. The percentage of phenotypic variance explained by rp1 differed significantly between families. The different genetic backgrounds of these families could explain the differences observed. The variance explained by rp1 in family 11 (63.4%) could also represent a better reflection of the percentage of variance explained by this QTL than that observed in family 4 (27.9%), as the two closest markers to rp1, 7P6P4 and 3P6P3, were fully informative in family 11 but only semi informative in family 4, hence, only half of the genetic information at these loci was available in family 4. The rp1 QTL appears to present a major effect since it explains more than 60% of the variance for pyrethroid resistance. In general, our results suggest that resistance to pyrethroids is a quantitative trait under the control of at least one QTL.
The family sizes used in this study were generally low because it is difficult to colonize An. funestus. The low sample sizes may have prevented the identification of additional QTLs of minor effect. Small sample sizes are more adequate to detect QTL with large effect such as rp1 but have limited power to detect smaller QTLs [22–24]. This may explain why only one confirmed QTL has been identified in this study. A denser marker map might also have a greater power to detect minor QTL and to separate any linked genes located within rp1. The numbers of actual QTL have been underestimated in studies involving QTL of large effect .
In family 4, two additional QTLs other than rp1 were identified on chromosome 2 by multiple-interval mapping. However when the genotypes of markers linked with these QTLs were tested for association with pyrethroid resistance using the χ2 goodness-of fit test, no association was found. Results from interval mapping and composite-interval mapping did not indicate the presence of QTL around the positions indicated by these MIM results. We believe that these additional QTLs are spurious and need to be confirmed in larger families with a greater number of genetic markers.
The region of the genome where the rp1 QTL was identified contains a cluster of cytochrome P450 genes belonging to the CYP6 family. Indeed, the two SNP markers (7P6P4 and 3P6P3) most closely linked to rp1 are located within two CYP6 genes. Elevated cytochrome P450 activity has been implicated as the major mechanism conferring pyrethroid resistance in An. funestus . The molecular basis of cytochrome P450-mediated insecticide resistance may involve cis or trans regulation of the transcription of a P450 gene . It is also possible, although uncommon, that resistance may be caused by a change in the P450 amino acid sequence increasing the affinity of a P450 for pyrethroids . In An. gambiae, a QTL has been mapped in the vicinity of a P450 gene (CYP6Z1) shown to be over-expressed in a resistant strain . To study the molecular mechanism of the pyrethroid resistance in the FUMOZ-R strain, and to see if additional minor pyrethroid resistance QTLs exist in the population, we are establishing advanced intercross lines and identifying additional genetic markers to better define the interval in which the rp1 QTL is mapped.
The present study reports the first QTL mapping study in An. funestus, one of the most important malaria vectors. This study suggests that metabolic resistance mechanism is playing a significant role in pyrethroid resistance in An. funestus since we did not found any evidence either from direct sequencing or QTL mapping, that mutations in the target site are involved in pyrethroid resistance in An. funestus. The rp1 QTL identified here represents a first step toward a fine mapping of genes involved in this resistance trait and by combining this approach with functional characterization of P450 genes on chromosome 2R we will elucidate the molecular basis of pyrethroid resistance in this malaria vector.
Mosquito strains and bioassays
The two strains of An. funestus used in this study were colonized by the Vector Control Reference Unit of the National Institute for Communicable Diseases in South Africa . The FUMOZ strain originated from southern Mozambique in 2000 and was initially a heterozygous permethrin resistant population. After selection of the parental strain with 0.1% lambda cyhalothrin , a highly resistant strain called FUMOZ-R was generated. The FANG colony which is completely susceptible to all pyrethroids was colonized from Calueque in southern Angola in 2002.
Reciprocal crosses using virgin FANG and FUMOZ-R females with males from the alternative strain were set-up. Blood fed, mated females were left to oviposit singly and the F1 were crossed to generate F2 progeny, producing iso-female lines. Bioassays were carried out on 1–3-day-old adults using the standard WHO bioassay procedure . Mosquitoes were exposed to 0.75% permethrin for 1 hour and mortality was recorded 24 h post-exposure. Dead mosquitoes were considered as susceptible to permethrin and separated from the living progeny which were considered as resistant. Parental females (P1), F1 progeny as well as F2 progeny from each family were collected and kept in silica gel Eppendorf tubes for later DNA extraction and genotype determination.
Genotyping of molecular markers
Scoring of resistance phenotype in F2 progeny
No of progeny
No of survivors genotyped
No. of dead genotyped
♀Fang × Fumoz♂
♀Fumoz × Fang♂
Single Nucleotide Polymorphisms (SNPs) genotyping
Fifteen SNP markers were previously identified in An. funestus . These were scored in the P1 and F1 parents in each of the 11 families and informative SNPs were scored either by the heated ligation oligonucleotide assay (HOLA)  or by single-base pair extension reaction using terminator dyes and CEQ8000 software or using the pyrosequencing method, as described below. An additional 11 informative SNP markers were identified by direct sequencing of PCR products from cytochrome P450 genes or from genes that had been previously mapped to An. funestus polytene chromosomes .
Details on SNP markers scored by pyrosequencing method
Forward primer 5'-3'
Reverse primer Biotine 5'-3'
Product size (bp)
Sequencing primer 5'-3'
Single-stranded biotinylated PCR products were prepared for the pyrosequencing reaction using a Vacuum Prep Tool (Biotage AB). The biotinylated PCR products were immobilized onto Streptavidin Sepharose high performance beads. For a single sample, 3 μl of Streptavidin Sepharose™ HP (Amersham) were added to 37 μl binding buffer (10 mM Tris-HCl pH 7.6, 2 M NaCl, 1 mM EDTA, 0.1% Tween 20) and mixed with 25 μl PCR product and 15 μl water on a mechanical shaker for 5 min at room temperature in a 96-well plate. The beads containing the immobilised templates were captured onto the filter probes after applying the vacuum and then washed with 70% ethanol for 5 sec, denaturizing solution (0.2 M NaOH) for 5 sec and washing buffer (10 mM Tris-Acetate pH 7.6) for 5 sec. The vacuum was switched off and the beads were released into a PSQ 96 well plate containing 38.4 μl annealing buffer (20 mM Tris-Acetate, 2 mM MgAc2 pH 7.6) and 1.6 μl of the corresponding sequencing primer (10 μM). Annealing was achieved by heating the samples to 80°C for 3 min and then allowing them to cool to room temperature. Pyrosequencing reactions were performed according to the manufacturer's instructions using the PSQ 96 SNP Reagent Kit (Biotage AB) and the genotype was determined using SNP Software (Biotage AB).
The JoinMap 2.0 package  was used to build a genetic linkage map for each individual family and for the combined genotype data from all families. Genotype data for each marker were tested for conformity to Mendelian ratio with a χ2 goodness-of-fit analysis using the JMSLA procedure. Loci were separated into linkage groups using the JMGRP1 and JMSPL procedures with minimum and maximum LOD thresholds of 0.0 and 6.0 respectively and LOD increments of 0.1. The JMREC program estimated pairwise cM distances between all pairs of informative loci in each linkage group. The JMMAP program estimated the maximum likelihood map using the Kosambi distances. DrawMap 1.1  software was then used to plot the genetic map.
Associations between genotypes at each locus and the resistance phenotype (dead or alive) were assessed using a contingency χ2 analysis. The null hypothesis was that susceptibility to pyrethroid is equal in each genotype class. The marginal probabilities were the frequencies of each genotype at a locus and the mortality and survival rates after pyrethroid exposure. For loci with a significant χ2 we analyzed the inheritance of the alleles at these loci. The a priori hypothesis was that higher mortality rate would occur among F2 individuals with one or both alleles inherited from the susceptible parent while lower mortality rate will be observed among individuals with alleles inherited from the resistant parent.
The JoinMap linkage map and the genotype/phenotype data were entered into Windows QTL Cartographer 2.5 . Interval mapping , composite-interval mapping  and multiple-interval mapping  procedures were implemented for each family. Interval mapping uses two observable flanking markers to construct an interval within which to search for QTL. The optimum LOD thresholds were estimated by permutation of trait and marker data 1000 times with a walking speed of 2 cM. Composite-interval mapping (CIM) tests whether an interval between two markers contains a QTL while simultaneously controlling for the effect of proximal QTLs located outside the interval. CIM was performed using the standard model with a control marker number of n p = 5 and a window size of w s = 10 cM. We used 1000 permutation analysis to determine the optimum significance threshold of the LOD. Multiple-interval mapping (MIM) analyzes multiple marker intervals simultaneously to fit multiple putative QTLs. An initial MIM model was estimated by forward and backward marker selection with a probability of a partial R2 set to 0.01. We then optimized QTL positions, searched for new QTLs, tested for existing QTL before saving the QTL model. The MIM model summary procedure estimated additive and dominance effects, the QTL likelihood ratio, the confidence interval around QTL positions, and partitioned the variance explained by QTL.
Cloning of partial sequences of sodium channel (kdr) and acetylcholinesterase (AChE) genes in An. funestus
We cloned partial fragments of the two major target sites of insecticides used in malaria control, the acetylcholinesterase gene (AChE) which is the target site of carbamate and organophosphate insecticides and the voltage gated sodium channel gene (kdr) which is the target site of pyrethroids and DDT. Primers designed from a conserved region of the S6 region of the voltage-gated sodium channel gene AAK2 (5'-TGC GGA GAA TGG ATY GAA TC-3') and Ask2 (CTT AGC CTT GCT TTT GTC AAA) were used to amplify a 200 bp fragment of this gene from An. funestus genomic DNA. The PCR conditions were 96°C for 1 min followed by 30 cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 30 s with a final 10 min extension at 72°C. The PCR product was ligated into pGEM-T easy vector and sequenced with M13 forward and reverse primers. From this 200 bp product, the Adf2 primer (5'-TGC AAA ATA GAG TCA TTG GTG AA-3') was designed and used with the Dip3 primer (5'-ATC ATC TTC ATC TTT GC-3')  to amplify a larger fragment of 1342 bp. This fragment was amplified and sequenced from individuals from the resistant FUMOZ-R strain and from field collected populations from Kenya, Mozambique and Malawi. A similar protocol was used to amplify a 1487 bp partial fragment of the acetylcholinesterase gene (AChE) in An. funestus using primers ANace-F1 (5'-CCG GGG GCG ACT ATG TGG AAC-3') and ANace-R3 (5'-GTT GCT GTT CGG GTT GTC CG-3'). The partial sodium channel and acetylcholinesterase of An. funestus sequences have been deposited in GenBank, accession numbers DQ534436 and DQ534435 respectively for kdr and AChE.
Polytene chromosome in situ hybridization
A subset of 12 partial P450 Cytochrome P450 genes were chosen for physical mapping as increased activity of the P450 enzymes has been demonstrated in the FUMOZ-R resistant strain  and because this is an important pyrethroid resistance mechanism in other insects [27, 37]. Half-gravid An. funestus females were collected from Kela village in Chikwawa district in southern Malawi. Ovaries were removed and preserved in Carnoy's fixative solution (3 ethanol: 1 glacial acetic acid by volume) and stored at -20° until use. The identity of specimens used in this study was confirmed by a rDNA diagnostic PCR which distinguishes four members of the An. funestus group . Probes were prepared with the GIBCO BRL in situ hybridization and detection system, using the manufacturer's recommended protocol. Using fluorescence in situ hybridization (FISH), An. funestus clones were mapped on the five arms of the polytene chromosome. After hybridization, chromosomes were washed in 0.2× SSC and counterstained with YOYO-1 and mounted in DABCO antifade solution (Sigma) . Fluorescent signals were detected using Zeiss LSM Pascal scanning confocal microscope.
This work was supported by the National Institute of Health (NIH) grant IU01 A1058271-01. The authors are grateful to Pie Muller and Sonia Casimiro for providing field samples of An. funestus.
- Gillies MT, De Meillon B: The anophelinae of Africa South of the Sahara. 1968, Johannesburg , The South African Institute for Medical Research, 2nd Edition:Google Scholar
- Hargreaves K, Koekemoer LL, Brooke BD, Hunt RH, Mthembu J, Coetzee M: Anopheles funestus resistant to pyrethroid insecticides in South Africa. Med Vet Entomol. 2000, 14 (2): 181-189. 10.1046/j.1365-2915.2000.00234.x.PubMedView ArticleGoogle Scholar
- Brooke BD, Kloke G, Hunt RH, Koekemoer LL, Temu EA, Taylor ME, Small G, Hemingway J, Coetzee M: Bioassay and biochemical analyses of insecticide resistance in southern African Anopheles funestus (Diptera: Culicidae). Bull Entomol Res. 2001, 91 (4): 265-272.PubMedView ArticleGoogle Scholar
- Charlwood JD, Smith T, Lyimo E, Kitua AY, Masanja H, Booth M, Alonso PL, Tanner M: Incidence of Plasmodium falciparum infection in infants in relation to exposure to sporozoite-infected anophelines. Am J Trop Med Hyg. 1998, 59 (2): 243-251.PubMedGoogle Scholar
- Mendis C, Jacobsen JL, Gamage-Mendis A, Bule E, Dgedge M, Thompson R, Cuamba N, Barreto J, Begtrup K, Sinden RE, Hogh B: Anopheles arabiensis and An. funestus are equally important vectors of malaria in Matola coastal suburb of Maputo, southern Mozambique. Med Vet Entomol. 2000, 14 (2): 171-180. 10.1046/j.1365-2915.2000.00228.x.PubMedView ArticleGoogle Scholar
- Hemingway J, Ranson H: Insecticide resistance in insect vectors of human disease. Annual revue of entomology. 2000, 45: 369-389.View ArticleGoogle Scholar
- Etang J, Manga L, Chandre F, Guillet P, Fondjo E, Mimpfoundi R, Toto JC, Fontenille D: Insecticide susceptibility status of Anopheles gambiae s.l. (Diptera: Culicidae) in the Republic of Cameroon. J Med Entomol. 2003, 40 (4): 491-497.PubMedView ArticleGoogle Scholar
- Vulule JM, Beach RF, Atieli FK, Mount DL, Roberts JM, Mwangi RW: Long-term use of permethrin-impregnated nets does not increase Anopheles gambiae permethrin tolerance. Med Vet Entomol. 1996, 10 (1): 71-79.PubMedView ArticleGoogle Scholar
- 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 (4): 417-422. 10.1111/j.1365-2915.2003.00460.x.PubMedView ArticleGoogle Scholar
- Ranson H, Nikou D, Hutchinson M, Wang X, Roth CW, Hemingway J, Collins FH: Molecular analysis of multiple cytochrome P450 genes from the malaria vector, Anopheles gambiae. Insect Mol Biol. 2002, 11 (5): 409-418. 10.1046/j.1365-2583.2002.00350.x.PubMedView ArticleGoogle Scholar
- Tijet N, Helvig C, Feyereisen R: The cytochrome P450 gene superfamily in Drosophila melanogaster: annotation, intron-exon organization and phylogeny. Gene. 2001, 262 (1-2): 189-198. 10.1016/S0378-1119(00)00533-3.PubMedView ArticleGoogle Scholar
- Amenya DA, Koekemoer LL, Vaughan A, Morgan JC, Brooke BD, Hunt RH, Ranson H, Hemingway J, Coetzee M: Isolation and sequence analysis of P450 genes from a pyrethroid resistant colony of the major malaria vector Anopheles funestus. DNA Seq. 2005, 16 (6): 437-445.PubMedView ArticleGoogle Scholar
- Bosio CF, Fulton RE, Salasek ML, Beaty BJ, Black WC: Quantitative trait loci that control vector competence for dengue-2 virus in the mosquito Aedes aegypti. Genetics. 2000, 156 (2): 687-698.PubMed CentralPubMedGoogle Scholar
- Zheng L, Cornel AJ, Wang R, Erfle H, Voss H, Ansorge W, Kafatos FC, Collins FH: Quantitative trait loci for refractoriness of Anopheles gambiae to Plasmodium cynomolgi B. Science. 1997, 276 (5311): 425-428. 10.1126/science.276.5311.425.PubMedView ArticleGoogle Scholar
- Ranson H, Jensen B, Wang X, Prapanthadara L, Hemingway J, Collins FH: Genetic mapping of two loci affecting DDT resistance in the malaria vector Anopheles gambiae. Insect Mol Biol. 2000, 9 (5): 499-507. 10.1046/j.1365-2583.2000.00214.x.PubMedView ArticleGoogle Scholar
- 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 (4): 379-386. 10.1111/j.0962-1075.2004.00495.x.PubMedView ArticleGoogle Scholar
- Hunt RH, Brooke BD, Pillay C, Koekemoer LL, Coetzee M: Laboratory selection for and characteristics of pyrethroid resistance in the malaria vector Anopheles funestus. Med Vet Entomol. 2005, 19 (3): 271-275. 10.1111/j.1365-2915.2005.00574.x.PubMedView ArticleGoogle Scholar
- Wondji CS, Hunt RH, Pignatelli P, Steen K, Coetzee M, Besansky N, Lobo N, Collins FH, Hemingway J, Ranson H: An integrated genetic and physical map for the malaria vector Anopheles funestus. Genetics. 2005, 171 (4): 1779-1787. 10.1534/genetics.105.044800.PubMed CentralPubMedView ArticleGoogle Scholar
- Sharakhov IV, Serazin AC, Grushko OG, Dana A, Lobo N, Hillenmeyer ME, Westerman R, Romero-Severson J, Costantini C, Sagnon N, Collins FH, Besansky NJ: Inversions and gene order shuffling in Anopheles gambiae and A. funestus. Science. 2002, 298 (5591): 182-185. 10.1126/science.1076803.PubMedView ArticleGoogle Scholar
- 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. 10.1111/j.1365-2583.2004.00452.x.PubMedView ArticleGoogle Scholar
- Zheng L, Benedict MQ, Cornel AJ, Collins FH, Kafatos FC: An integrated genetic map of the African human malaria vector mosquito, Anopheles gambiae. Genetics. 1996, 143 (2): 941-952.PubMed CentralPubMedGoogle Scholar
- Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121 (1): 185-199.PubMed CentralPubMedGoogle Scholar
- Zeng ZB: Theoretical basis for separation of multiple linked gene effects in mapping quantitative trait loci. Proc Natl Acad Sci U S A. 1993, 90 (23): 10972-10976. 10.1073/pnas.90.23.10972.PubMed CentralPubMedView ArticleGoogle Scholar
- Zeng ZB: Precision mapping of quantitative trait loci. Genetics. 1994, 136 (4): 1457-1468.PubMed CentralPubMedGoogle Scholar
- Gorman MJ, Severson DW, Cornel AJ, Collins FH, Paskewitz SM: Mapping a quantitative trait locus involved in melanotic encapsulation of foreign bodies in the malaria vector, Anopheles gambiae. Genetics. 1997, 146 (3): 965-971.PubMed CentralPubMedGoogle Scholar
- Liu N, Scott JG: Increased transcription of CYP6D1 causes cytochrome P450-mediated insecticide resistance in house fly. Insect Biochem Mol Biol. 1998, 28 (8): 531-535. 10.1016/S0965-1748(98)00039-3.PubMedView ArticleGoogle Scholar
- Scott JG: Cytochromes P450 and insecticide resistance. Insect Biochem Mol Biol. 1999, 29 (9): 757-777. 10.1016/S0965-1748(99)00038-7.PubMedView ArticleGoogle Scholar
- 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. 10.1016/S0378-1119(03)00763-7.PubMedView ArticleGoogle Scholar
- WHO: Test procedures for insecticide resistance montoring in malaria vectors, bio-efficacy and persistence of insecticides on treated surfaces. Edited by: WHO/CDS/CPC/MAL/98.12 D. 1998, Geneva, Switzerland , World Health OrganizationGoogle Scholar
- Collins FH, Mendez MA, Rasmussen MO, Mehaffey PC, Besansky NJ, Finnerty V: A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex. Am J Trop Med Hyg. 1987, 37 (1): 37-41.PubMedGoogle Scholar
- Lynd A, Ranson H, McCall PJ, Randle NP, Black WC, Walker ED, Donnelly MJ: A simplified high-throughput method for pyrethroid knock-down resistance (kdr) detection in Anopheles gambiae. Malar J. 2005, 4 (1): 16-10.1186/1475-2875-4-16.PubMed CentralPubMedView ArticleGoogle Scholar
- Stam P, Van Ooijen JW: JoinMap (tm) version 2.0: Software for the calculation of genetic linkage maps. 1995, Wageningen, The Netherlands , CPRO-DLOGoogle Scholar
- van Ooijen JW: DrawMap: a computer program for drawing genetic linkage maps. J Hered. 1994, 85 (1): 66-PubMedGoogle Scholar
- Wang S, Basten CJ, Zeng ZB: Windows QTL Cartographer 2.5. 2005, Raleigh, NC. , Department of Statistics, North Carolina State UniversityGoogle Scholar
- Zeng ZB, Kao CH, Basten CJ: Estimating the genetic architecture of quantitative traits. Genet Res. 1999, 74 (3): 279-289. 10.1017/S0016672399004255.PubMedView ArticleGoogle Scholar
- 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. 10.1046/j.1365-2583.1998.72062.x.PubMedView ArticleGoogle Scholar
- Maitra S, Dombrowski SM, Waters LC, Ganguly R: Three second chromosome-linked clustered Cyp6 genes show differential constitutive and barbital-induced expression in DDT-resistant and susceptible strains of Drosophila melanogaster. Gene. 1996, 180 (1-2): 165-171. 10.1016/S0378-1119(96)00446-5.PubMedView ArticleGoogle Scholar
- Koekemoer LL, Lochouarn L, Hunt RH, Coetzee M: Single-strand conformation polymorphism analysis for identification of four members of the Anopheles funestus (Diptera: Culicidae) group. J Med Entomol. 1999, 36 (2): 125-130.PubMedView ArticleGoogle Scholar
- Sharakhov I, Braginets O, Grushko O, Cohuet A, Guelbeogo WM, Boccolini D, Weill M, Costantini C, Sagnon N, Fontenille D, Yan G, Besansky NJ: A microsatellite map of the African human malaria vector Anopheles funestus. J Hered. 2004, 95 (1): 29-34. 10.1093/jhered/esh011.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.