- Open Access
Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiae molecular form identification
Malaria Journal volume 10, Article number: 215 (2011)
Anopheles gambiae M and S molecular forms, the major malaria vectors in the Afro-tropical region, are ongoing a process of ecological diversification and adaptive lineage splitting, which is affecting malaria transmission and vector control strategies in West Africa. These two incipient species are defined on the basis of single nucleotide differences in the IGS and ITS regions of multicopy rDNA located on the X-chromosome. A number of PCR and PCR-RFLP approaches based on form-specific SNPs in the IGS region are used for M and S identification. Moreover, a PCR-method to detect the M-specific insertion of a short interspersed transposable element (SINE200) has recently been introduced as an alternative identification approach. However, a large-scale comparative analysis of four widely used PCR or PCR-RFLP genotyping methods for M and S identification was never carried out to evaluate whether they could be used interchangeably, as commonly assumed.
The genotyping of more than 400 A. gambiae specimens from nine African countries, and the sequencing of the IGS-amplicon of 115 of them, highlighted discrepancies among results obtained by the different approaches due to different kinds of biases, which may result in an overestimation of MS putative hybrids, as follows: i) incorrect match of M and S specific primers used in the allele specific-PCR approach; ii) presence of polymorphisms in the recognition sequence of restriction enzymes used in the PCR-RFLP approaches; iii) incomplete cleavage during the restriction reactions; iv) presence of different copy numbers of M and S-specific IGS-arrays in single individuals in areas of secondary contact between the two forms.
The results reveal that the PCR and PCR-RFLP approaches most commonly utilized to identify A. gambiae M and S forms are not fully interchangeable as usually assumed, and highlight limits of the actual definition of the two molecular forms, which might not fully correspond to the two A. gambiae incipient species in their entire geographical range. These limits are discussed and operational suggestions on the choice of the most convenient method for large-scale M- and S-form identification are provided, also taking into consideration technical aspects related to the epidemiological characteristics of different study areas.
The mosquito vector species responsible for most Plasmodium falciparum- malaria transmission in sub-Saharan Africa, Anopheles gambiae sensu stricto (hereafter A. gambiae), is ongoing a process of ecological diversification and adaptive lineage splitting which is changing patterns of malaria transmission and affecting vector control strategies in West Africa [1–4]. Two morphologically indistinguishable incipient species (provisionally named M and S molecular forms) have been described within A. gambiae, based on form-specific single nucleotide polymorphisms (SNPs) on the intergenic spacer (IGS) and internal transcribed spacer (ITS) regions of multicopy ribosomal DNA (rDNA) located on the X-chromosome [5, 6]. S-form is distributed across sub-Saharan Africa and breeds mostly in association with rain-dependent pools and temporary puddles. M-form distribution overlaps with that of S-form in West and Central Africa, but the former form is apparently absent east of the Great Rift Valley; it is able to exploit relatively more permanent breeding sites, often closely associated with human activities, such those created by irrigation, rice cultivation and urbanization [2, 3, 7, 8]. This adaptation allows the M-form to breed throughout the year, thus causing a shift from seasonal to year-round malaria transmission. Importantly, genetic traits conferring resistance to insecticides commonly used against these vectors are differently distributed between the two forms [9, 10].
Genetic divergence between M and S forms has been recently shown to be widespread across the genome [11, 12]. However, the most widely used methods for M and S specimen identification are based on genotyping procedures for the form-specific SNPs in the IGS rDNA region on the centromere of the X-chromosome. These are performed either by PCR using form-specific primers [13, 14] or PCR-RFLP [15–17] (Figure 1). More recently, a PCR-method to detect the M-specific insertion of a SINE200 (short interspersed transposable element), mapping about 1 Mb apart from IGS SNPs in the chromosome-X centromeric region, was also developed . Results from large scale identification of A. gambiae field specimens by means of any of these approaches highlighted the virtual absence of hybrid M/S patterns, thus contributing substantial evidence of M and S reproductive isolation in nature. However, high frequencies of M/S IGS-patterns have been recently reported from The Gambia  and Guinea Bissau . A preliminary comparison of different identification approaches in samples from these westernmost geographical areas highlighted inconsistencies in the results and the occurrence of possible biases due to the routine practice of identifying M and S mosquitoes based on a single assay .
The results of a large-scale comparative analysis of four widely used PCR or PCR-RFLP genotyping methods for M and S identification provide operational recommendations to medical entomologists dealing with M and S identification in the Afrotropical region.
Samples, genotyping and sequencing
Table 1 lists the indoor-resting female samples processed in this study, the dates of collections and the numbers of specimens genotyped. Figure 2 shows the location of the sampling sites. The specimens from Guinea Bissau and The Gambia  were selected based on inconsistent results from IGS  and SINE200 insertion  genotyping.
DNA was extracted from either legs or other parts of the carcasses not including the abdomen, to avoid possible biases associated to the risk of contamination with DNA from sperm harboured in spermathecae.
Anopheles gambiae samples were first identified based on results from PCR-RFLP approach (from now on PCR-RFLP581) recognizing a T/C SNP (T = M-form; C = S-form) at position 581 of IGS rDNA region (from now on IGS581). Subsequently, the following genotyping approaches were applied (Figure 1 shows the position of primers and restriction sites on the IGS amplicon sequence): i) the PCR-RFLP approach (from now on PCR-RFLP690) recognizing a A/T SNP (A = M-form; T = S-form) at position 690 of IGS rDNA region (from now on IGS690); ii) the PCR approach using allele-specific primers designed to detect the IGS581 SNP (from now on AS-PCR ); iii) the PCR approach based on the specific and irreversible single-locus insertion of a SINE200 transposable element in the X-chromosome centromeric region (from now on SINE-PCR ), about 1 Mb from the IGS rDNA region including the IGS form-specific SNPs. A sub-sample of specimens from Angola, Burkina Faso, The Gambia and Guinea Bissau were also identified by a PCR approach utilizing Intentional Mismatch Primers containing single base mismatches at the third nucleotide from their 3' end (from now on IMP-PCR).
An IGS fragment of 367 bp (from now on "IGS-amplicon") was amplified using UN and GA primers by Fanello et al (Figure 1) from selected specimens and sequenced using ABI Big Dye Terminator v.2 chemistry and an ABI Prism 3700 DNA Analyser. Chromatograms were inspected for double peaks by eye. PCR and sequence analyses were carried out in Rome and/or Lisbon. Selected samples were analysed in both laboratories for results validation.
QSVanalyzer software - which allows the extraction of quantitative sequence variant (QSV) information from sequence electropherograms - was applied to estimate the relative proportions of the double peaks (i.e., c opy n umber p roportions: CNP)) observed in electropherograms of IGS amplicon at positions 581  (hereafter CNP581) and 690  (hereafter CNP690) in sequences of the IGS locus from single A. gambiae specimens. The programme analyses each trace and adjusts it in relation to the peak heights of upstream/downstream nucleotides, allowing rapid batchwise analysis of DNA sequence traces for estimation of the relative proportions of two QSVs at a given site. Kruskal-Wallis and multiple comparison tests were carried out by STATISTICA 6.1 (StatSoft, Inc. 2003).
Four-hundred-twenty-four A. gambiae specimens from nine African countries by PCR-RFLP581, PCR-RFLP690, AS-PCR and SINE-PCR were genotyped (Table 1). Most (97%) of the specimens were consistently identified by all approaches in samples from Tanzania to Ghana (N = 250), while the percentage of consistent identifications was lower (46%) in samples from Guinea Bissau and The Gambia - which were selected based on previous inconsistent results from PCR-RFLP581 and SINE-PCR  - and in samples from Burkina Faso (41%). Inconsistent identifications were confirmed at least twice by PCR and PCR-RFLP genotyping carried out in different laboratories. The IGS-amplicon was sequenced in 115 specimens (Angola N = 5; Cameroon N = 1; Nigeria N = 1; Burkina Faso N = 16; Mali N = 1; Guinea Bissau N = 62; The Gambia N = 29). The latter samples plus additional 110 specimens (Tanzania N = 11; Angola N = 21; Cameroon N = 9; Nigeria N = 11; Burkina Faso N = 20; Mali N = 12; Ghana N = 18; Guinea Bissau N = 1; The Gambia N = 7), were also genotyped by IMP-PCR.
The results were as follows:
Tanzania. All specimens were consistently identified as S-form by the four approaches utilized. Eleven identifications were confirmed also by IMP-PCR.
Angola. All individuals from Cabinda (N = 32) were identified as S-form by the four approaches. Twenty-seven out of 32 individuals from Luanda were identified as M-form by all approaches, while five of them showed a MS690 heterozygous pattern. Sequence analysis of these specimens revealed the presence of an A/C polymorphism at position 690 (instead of the expected A/T polymorphism, corresponding to the MseI restriction site), which does not allow the form-specific cleavage of the PCR-amplified band. Ten M- (including the above 5) and 11 S-identifications were confirmed also by IMP-PCR.
Cameroon. Twenty individuals were consistently identified as M-form and 26 as S-form. One single specimen showed a MS690 pattern, but a M-form pattern by the other approaches. IMP-PCR and IGS-sequencing confirmed the MM genotype.
Nigeria. Nine individuals were consistently identified as M-form and 17 as S-form. One single specimen showed a S-pattern by both PCR-RFLPs and MS by AS-PCR. IMP-PCR and IGS-sequencing confirmed the SS genotype.
Burkina Faso. Eleven individuals were consistently identified as M-form and 13 as S-form; 32 specimens (MM581-MM690 = 17 and SS581-SS690 = 15) showed a MS-pattern only by AS-PCR. The PCR-RFLP genotypes were confirmed either by IMP-PCR and/or by sequencing of the IGS-amplicon, which did not reveal any MS-heterozygous pattern. Two MM581-MS690 specimens were found, one showing MS-AS-PCR, the other one showing M-AS-PCR pattern. Both were genotyped as MM by IMP-PCR and IGS-sequencing.
Mali. Three individuals were consistently identified as M-form and 35 as S-form. One single specimen showed a MS-pattern by AS-PCR and a M-form pattern by the other approaches: both IMP-PCR and IGS-sequencing confirmed the MM genotype.
Ghana. All specimens were consistently identified as S-form by the four approaches utilized.
Guinea Bissau. Results from PCR-RFLP581, PCR-RFLP690, AS-PCR and sequencing are shown in Table 2: 61% (41/67) of the specimens showed consistent genotypes with the three approaches (i.e. 5 M, 29 MS, 7 S; Table 2, lines 1, 5, 10). PCR-RFLP581and PCR-RFLP690 patterns were congruent in 73% (49/67, lines 1, 2, 5,7,10) of the specimens: the remaining were either MM581-MS690 (10%; lines 3 and 4) or MS581-SS690 (16%, lines 6, 8 and 9). No opposite MM581-SS690 nor SS581-MM690 identifications were observed. An AS-PCR heterozygous MS pattern was obtained from: i) all specimens genotyped as heterozygotes by at least one PCR-RFLP, with the exception of four MS581-SS690 specimens showing a SS AS-PCR pattern (line 9); ii) 6 out of 11 MM581-MM690 homozygotes (line 2); and iii) 2 out of 9 SS581-SS690 (line 7). A match between the results of SINE-PCR and IGS genotypes (as defined by either consistent results of the different genotyping approaches and/or by direct sequencing of the IGS-amplicon, see below) was observed in 69% of the specimens (46/67: 10 MM, 25 MS and 11 SS). The mismatches were either due to SS (16/21) and MM (1/21) SINE-homozygotes with a heterozygous MSIGS genotype, and to MS SINE-heterozygotes with a homozygous MMIGS (1/21) or SSIGS (3/21) genotype.
The Gambia. Results from PCR-RFLP581, PCR-RFLP690, AS-PCR and sequencing are shown in Table 2: 69% (34/49) of the specimens analysed showed consistent genotypes with the three approaches (i.e. 9 M, 16 MS, 9 S; Table 2, lines 11, 16, 20). PCR-RFLP581and PCR-RFLP690 patterns were congruent in 86% (42/49) of the specimens. The remaining were either MM581-MS690 (8%; lines 13-15) or MS581-SS690 (6%; line 17-18). No MM581-SS690 nor SS581-MM690 genotypes were observed. An AS-PCR heterozygous MS pattern was obtained for: i) all specimens genotyped as heterozygotes by at least one PCR-RFLP, with the exception of one MM581-MS690 specimens showing a MM AS-PCR pattern (line 13); ii) 7 out of 16 MM581-MM690 homozygotes (line 12) and 1 out of 10 SS581-SS690 (line 19). A match between the results of SINE-PCR and IGS genotype (as defined from consistent results of the different genotyping approaches and/or by direct sequencing of the IGS-amplicon, see below) was observed in 72% of the specimens (N = 36: 16 MM, 9 MS and 11 SS). The mismatches were either due to SS (9/13) and MM (1/13) SINE-homozygotes with a heterozygous MSIGS genotype or to MS SINE-heterozygotes with MMIGS(2/13) or SSIGS(1/13) genotypes.
The electrophoregrams of the overall sequenced sample were further scored by QSV analyser  to quantify the proportion of sequences containing C versus T (M-form = T; S-form = C) or A versus T (M-form = A; S-form = T), based on relative peak heights at position 581 (CNP581) and 690 (CNP690), respectively. As expected based on the proximity of the IGS581 and IGS690 SNPs, the CNP scores were strongly correlated (r = 0,97 p << 0.001). The median CNP scores of the two SNPs were significantly different among the 7 IGS-types classified based on both IGS581 and IGS690 PCR-RFLPs (IGS581: KW-H = 116, p << 0.001; IGS690: KW-H = 109, p << 0.001), with specimens identified as MM and SS by both PCR-RFLPs showing median CNP scores near 0.1 and 0.9, respectively (as expected if one allele is fixed) and specimens identified as MS by both PCR-RFLPs showing intermediate scores (Figure 3). The heterozygotes (MS581-MS690) were statistically different from the homozygotes (pairwise comparisons: p < 0.001; Additional file 1). Interestingly, specimens from Guinea Bissau and The Gambia, characterized by inconsistent PCR-RFLP patterns showed intermediate scores between those of M or S homozygotes (MM581-MM690 or SS581-SS690) and MS heterozygotes (MS581-MS690), suggesting that these specimens are characterized by an unequal number of copies of M- and S-arrays (pairwise comparisons: p > 0.05; Additional file 1). CNPs scores of individuals subdivided based on SINE-genotypes revealed that some SINE-XMM and SINE-XSS individuals are characterized by an unequal number of copies of M- and S-arrays . It is relevant to note that the interpretation of single versus double peaks at the two IGS diagnostic sites determined either by eye inspection of chromatograms or by CNP score (i.e. SNP581: homozygous T, CNP < 0.15; homozygous C, CNP > 0.85; SNP690: homozygous A, CNP < 0.18; homozygous T, CNP > 0.85) were consistent in 97% of the cases. IMP-PCR confirmed the results from sequencing, with the exception of 2 S-form individuals from The Gambia genotyped as MS by both IMP-PCR and AS-PCR.
Sequence analysis also showed that all samples were homozygotes (G) at position IGS580, where Favia et al reported an addtional form-specific SNP (M-form = G; S-form = A). Moreover, alignment of IGS amplicon from sequenced individuals from Guinea Bissau (M-form: N = 31; S-form: N = 32), showed 4 IGS-polymorphic sites in addition to IGS581 and IGS690: i) site IGS485 was found heterozygous (C/T) in one M/S specimen; ii) site IGS489 was found heterozygous (G/A) in 2 M and in 1 M/S specimens; iii) site IGS491 was found homozygous (T) in one M specimen and heterozygous (T/C) in 2 M/S specimens. Finally, a A/G polymorphism at site IGS612 was found in linkage with both IGS581 and IGS690 diagnostic sites, as already reported by Oliveira et al.
Since their description 10 years ago , A. gambiae M and S molecular forms have been the focus of extensive field studies aimed to evaluate their macro- and micro-geographic distribution and to analyse their population genetics, their relative role as malaria vectors and their resistance to insecticides used in malaria vector control campaigns. All these studies have exploited different approaches for M and S identification - based on either direct PCR-amplification of form-specific bands by allele-specific primers designed in the same region (AS-PCR, cited in 17 out of 56 papers since 2006, see Additional file 2) or the restriction of form-specific SNPs in the IGS rDNA region (IGS581 and IGS690 PCR-RFLPs: cited in  and , respectively, Additional file 2). In most of these papers only a single identification approach was used, while in four of them both AS-PCR and PCR-RFLP581 methods were applied. Moreover, a few authors (nine papers; Additional file 2) still exploited the original approach developed by Favia et al, which is based on the same restriction site of Fanello et al, but requires the pre-identification of A. gambiae s.s. specimens. Overall, so far, the general attitude has been to consider all these approaches fully interchangeable; however, recent data from the westernmost extreme of M and S range (i.e. The Gambia and Guinea Bissau ) have highlighted that this assumption is not correct. The results here presented confirm this conclusion and allow to pinpoint the bases of the observed inconsistencies among results of the three approaches, as follows.
Incorrect match of M and S specific primers used in the AS-PCR approach
A higher number of MS heterozygous patterns resulted from AS-PCR than from PCR-RFLP genotyping, mostly in specimens from Burkina Faso, Guinea Bissau and The Gambia. This is due to a low specificity of the AS-PCR approach (using form-specific primers differing only for the SNP variant at 3' end) which is affected by the inability of this single 3' mismatch to prevent extension of the non-specific primer by the polymerase . In fact, this low specificity has been recently circumvented using primers containing an additional intentional mismatch at the third nucleotide from the 3' end which increases the power of Taq polymerase to extend from the 'right' primer and to partly optimize the reaction thermodynamics when both primers anneal on the template, thus providing more power to identify MS hybrids (IMP-PCR ). The comparison between the results of the AS-PCR and of IMP-PCR highlighted a higher specificity of the latter. In fact, IMP-PCR produced patterns consistent with those obtained either by sequencing or by the two PCR-RFLP approaches in all specimens tested (N = 146), with the exception of two S-form specimens from Gambia genotyped as MS by IMP-PCR (Table 2, line 18).
Presence of polymorphisms in the recognition sequence of restriction enzymes used in the PCR-RFLP approaches
An A/C heterozygous pattern was observed in the recognition sequence of the enzyme utilised in IGS690 PCR-RFLP (i.e. MseI) in five out of 32 M-form specimens from Angola. This polymorphism did not allow the complete cleavage of the M-specific PCR-amplified band, thus producing a false heterozygous MS690 pattern.
Incomplete cleavage during the restriction reaction
A few specimens from Burkina Faso, Cameroon, Guinea Bissau and The Gambia were incorrectly genotyped as MS by PCR-RFLP (IGS581: N = 8; IGS690: N = 5), due to incomplete digestion of the PCR-amplified fragment during restriction. A second round of PCR-RFLP reactions did not change the observed PCR-pattern and the specimens were confirmed to be homozygous at each site by sequencing.
Presence of different number of copies of M and S-specific IGS-arrays in single individuals
This has been already hypothesized by Caputo et al based on the inconsistent results from PCR-RFLP581 and SINE-PCR on samples from Guinea Bissau and The Gambia, where a secondary contact zone between the two molecular forms has been hypothesized based on the high frequencies of MS putative hybrids reported [19, 20]. The results obtained confirm this hypothesis and highlight the technical bias which emerged when the same samples were identified by PCR-RFLP690. In fact, the restriction enzyme used for the PCR-RFLP581 (i.e. HhaI) recognizes a S-specific restriction site, while the enzyme used for the PCR-RFLP690 (i.e MseI) recognizes a M-specific restriction site. It is possible to hypothesize that the PCR-amplification of individuals characterized by a number of copies of the M-IGS type higher than of S-IGS type exponentially increases this difference, producing a strong M581 band and a weak S581 one. The latter may not be visible on the agarose gel after the restriction step resulting in a MM581/MS690 RFLP pattern. Conversely, individuals characterized by a number of copies of the S-IGS type higher than of M-IGS type are likely to produce a MS581/SS690 RFLP pattern. This hypothesis is further supported by the relative high frequency of MM581/MS690 (9%) and MS581/SS690 (12%) specimens in the sample analysed, and by the absence of SS581/MS690 and MS581/MM690 genotypes. The QSV analysis of IGS sequences confirms that MM581/MS690 or MS581/SS690 individuals have proportions of array copy number intermediate between those of either MM581/MM690 and MS581/MS690 or SS581/SS690 and MS581/MS690 individuals, respectively (Table 2).
The comparison between the results of the IGS-genotyping (including direct sequencing, in case of inconsistencies among the approaches utilized) and of SINE-PCR showed consistent identifications in all samples, with the exception of those from Guinea Bissau and The Gambia. In these populations mismatches were observed, mostly due to SS and MM SINE-homozygotes with a heterozygous MSIGS genotype or, less frequently, to MS SINE-heterozygotes with MMIGS and SSIGS genotypes, while no opposite MM-SINE/SSIGS or SS-SINE/MMIGS were found. As discussed in Caputo et al, the former individuals are likely to represent Fn progenies of inter-form crosses occurring in this "secondary contact zone", where the reproductive isolation mechanisms between M- and S-forms appear to be less effective than in the rest of the molecular forms sympatric distribution range. In fact, discrepancies between results from PCR-RFLP581 and SINE-PCR led to hypothesize that the high frequencies of MS581 patterns found in Guinea Bissau and in The Gambia were due to the presence of both M- and S-arrays in the multi-copy IGS rDNA region of single individuals, suggesting inter-locus recombination . In this scenario, the SINE-PCR genotyping allows to discriminate putative MS hybrids from progenies of Fn-backcrosses (i.e. MM or SS SINE-homozygotes showing both M- and S-specific IGS arrays). In fact, the SINE-PCR genotyping of four MS581 specimens reported in della Torre et al (from Benin, Mali, Guinea and The Gambia) confirmed their putative hybrid origins.
On the other hand, the finding of high frequencies of consistent MS IGS/SINE patterns in larval samples from Burkina Faso, led Riehle et al. to carry out a deeper genetic characterization of these individuals and to hypothesize that they may represent a new A. gambiae "sub-form" highly differentiated from M and S. This "sub-form" seems to be also characterized by a MS SINE-polymorphism in Hardy-Weinberg equilibrium consistent with IGS-patterns, a very unexpected scenario which needs to be taken into consideration when speculating on the origin of this putative "sub-form". In fact, based on their evolutionary dynamics, both IGS and SINE markers are expected to undergo rapid fixation in a randomly mated diverging taxa rather than being found at equilibrium in a taxon separated from M- and S-form.
Overall, the results here presented, as well as those by Riehle et al, do not only highlight limits in the approaches currently applied to discriminate M- and S-forms, but also on the actual definition of the two molecular forms, which might not fully correspond to the two A. gambiae incipient species in their entire geographical range. The M and S molecular forms are, in fact, defined specifically based on SNPs in the IGS region, which were initially used to discriminate between Mopti and Savanna/Bamako chromosomal forms in Mali and Burkina Faso  and, later, to identify two incipient species in other geographical regions, where the correlation with specific karyotypes was more complex [1, 5, 7, 25]. Since their initial description, all data on the genetic, ecological and behavioural divergence of M and S forms were obtained based on the IGS diagnostics, leading to a general acceptance of the IGS-SNPs as form-specific characters possibly linked to genes or genomic regions instrumental to the speciation process. This view was reinforced by the fact that the IGS lies within X-chromosome centromeric region, where most genetic divergence between M- and S-forms is observed [11, 12, 25, 26] and by the consistent almost complete absence of MSIGS genotypes in nature. The finding of different number of copies of M- and S- IGS-arrays in single individuals from the western extreme of the molecular form range  highlighted how the genetic definition of the two A. gambiae incipient species is not fully tenable along their entire range. The recent sequencing of the genome of M- and S-colonies from Mali  and the availability of affordable SNP micro-array platforms , will probably allow in the near future a relatively easy processing of A. gambiae populations from the entire range. Moreover, the likely detection of multiple markers along their M-and S-form genome and their association will possibly allow a more precise definition of the two incipient species, as in the case of the allelic variant of TEP1 immune gene found to be fixed in M samples from Mali and Burkina Faso but absent in sympatric S populations .
The results obtained reveal that the PCR and PCR-RFLP approaches most commonly utilized to identify A. gambiae molecular forms are not fully interchangeable, as usually assumed. Different kinds of technical biases have been highlighted, which may result in an overestimation of MS putative hybrids. This is particularly relevant in settings of realised gene flow between molecular forms, such as the areas at the extreme West African distribution of A. gambiae. Moreover, the IMP-PCR developed by Wilkins et al, and so far applied almost exclusively on laboratory samples, was shown to be more specific than AS-PCR thus encouraging its exploitation in large scale screenings of field A. gambiae samples. However, the risk of biases due to local polymorphisms in the annealing sequences should be always be taken into consideration.
From an operational perspective, it needs to be highlighted that the choice of the most convenient method for large-scale M- and S-form identification, also depends from technical considerations (e.g. laboriousness of the different approaches) and from the sympatric presence of other members of the A. gambiae complex in the study area. In fact, only IGS581 PCR-RFLP allows the simultaneous identification of all species and molecular forms and could be the method of choice whenever the presence of other A. gambiae complex members (i.e. Anopheles melas, Anopheles merus, Anopheles quadriannulatus and/or Anopheles bwambae) cannot be excluded. Alternatively, the IMP-PCR approach could be used after A. gambiae s.s. specimens are identified by the species-specific PCR developed by Scott et al, thus avoiding the risks connected to the restriction step in IGS581 PCR-RFLP. The choice of one or the other approach should also be linked to the relative frequencies of A. gambiae s.s. in the sample (i.e. if this frequency is low, the species-specific PCR + IMP-PCR approach could be more convenient, as only few specimens would require to be identified by IMP-PCR; if high, the IGS581 PCR-RFLP could be a more straightforward approach). It should be noted that the IGS690 PCR-RFLP is more sensitive in cases of degraded DNA samples and that it could simultaneously identify M- and S-forms and A. arabiensis. On the other hand, the use of AS-PCR would require preliminary identification of A. gambiae s.s. specimens and is shown to be subject to more biases than the other approaches. The PCR-RFLP originally developed by Favia et al and still recently utilized by some authors (see Additional file 2) is comparatively less suitable for large-scale studies as it requires previous complex species identification and yet it targets the same SNP as the PCR-RFLP581. The latter method is more efficient since it allows for simultaneous species and molecular form identification on a much smaller amplicon (367 bp compared to 1.3 kb ).
It is also proposed that, due to the straightforward amplification of small DNA fragments (i.e. 249 and 479 bp for S- and M-forms, respectively), SINE-PCR could be conveniently applied to easily identify M- and S-forms (even without preliminary species-specific PCR identification in areas where exclusive sympatry with A. arabiensis is found). However, it is important to keep in mind that the M-form specific SINE insertion is a character linked to the IGS-SNPs defining the M- and S-forms along most of their range, but with a different evolutionary history (i.e. its origin and rapid fixation in M-form).
Finally, it is recommend to apply more than one genotyping approach (and/or sequencing of the IGS-amplicon) when identifying samples from previously unexplored geographic areas within M- and S-form sympatric range and whenever MS hybrid patterns are observed (Additional file 3). In this latter case, in fact, presence of both M- and S-specific IGS arrays in single individuals could lead to a misleading calculation of frequency of hybridization between M and S forms, as shown in populations from Guinea Bissau and The Gambia, where the simultaneous use of SINE-PCR allowed a better understanding of the local situation.
- IGS 581 SNP :
T/C Single Nucleotide Polymorphism at position 581 of IGS rDNA region
- IGS 690 SNP :
A/T Single Nucleotide Polymorphism at position 690 of IGS rDNA region
- PCR-RFLP 581 :
PCR-RFLP recognising IGS581 SNP
- PCR-RFLP 690 :
PCR-RFLP recognising IGS690 SNP
- AS-PCR :
PCR based on Allele-Specific primers designed on IGS581 SNP
- SINE-PCR :
PCR approach based on the specific and irreversible single-locus insertion of a Short Interspersed Transposable Element in the X-chromosome centromeric region
- IMP-PCR :
PCR based on Intentional Mismatch Primers annealing on IGS581 and IGS690 SNPs
- CNP 581 :
Copy Number Proportion of T/C alleles at positions 581 in the IGS amplicon
- CNP 690 :
Copy Number Proportion of A/T alleles at positions 690 in the IGS amplicon
- MM 581 :
SS581, MS581: M- and S-form specific genotypes at IGS581 SNP
- MM 690 :
SS690, MS690: M- and S-form specific genotypes at IGS690 SNP
- MM-SINE :
SS-SINE, MS-SINE: M- and S-form specific SINE genotypes.
della Torre A, Costantini C, Besansky NJ, Caccone A, Petrarca V, Powell JR, Coluzzi M: Speciation within Anopheles gambiae--the glass is half full. Science. 2002, 298: 115-117. 10.1126/science.1078170.
Costantini C, Ayala D, Guelbeogo WM, Pombi M, Some CY, Bassole IH, Ose K, Fotsing JM, Sagnon N, Fontenille D, Besansky NJ, Simard F: Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae. BMC Ecol. 2009, 9: 16-10.1186/1472-6785-9-16.
Simard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, Fotsing JM, Fontenille D, Besansky NJ, Costantini C: Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009, 9: 17-10.1186/1472-6785-9-17.
White BJ, Cheng C, Simard F, Costantini C, Besansky NJ: Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae. Mol Ecol. 2010, 19: 925-939. 10.1111/j.1365-294X.2010.04531.x.
della Torre A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, Petrarca V, Coluzzi M: Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol Biol. 2001, 10: 9-18. 10.1046/j.1365-2583.2001.00235.x.
Gentile G, Slotman M, Ketmaier V, Powell JR, Caccone A: Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 25-32. 10.1046/j.1365-2583.2001.00237.x.
della Torre A, Tu Z, Petrarca V: On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms. Insect Biochem Mol Biol. 2005, 35: 755-769. 10.1016/j.ibmb.2005.02.006.
Lehmann T, Diabate A: The molecular forms of Anopheles gambiae: a phenotypic perspective. Infect Genet Evol. 2008, 8: 737-746. 10.1016/j.meegid.2008.06.003.
Djogbenou L, Dabire R, Diabate A, Kengne P, Akogbeto M, Hougard JM, Chandre F: Identification and geographic distribution of the ACE-1R mutation in the malaria vector Anopheles gambiae in south-western Burkina Faso, West Africa. Am J Trop Med Hyg. 2008, 78: 298-302.
Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, Caccone A, Donnelly MJ, Petrarca V, Simard F, Pinto J, della Torre A: Distribution of knock-down resistance mutations in Anopheles gambiae molecular forms in west and west-central Africa. Malar J. 2008, 7: 74-10.1186/1475-2875-7-74.
Lawniczak MK, Emrich SJ, Holloway AK, Regier AP, Olson M, White B, Redmond S, Fulton L, Appelbaum E, Godfrey J, Farmer C, Chinwalla A, Yang SP, Minx P, Nelson J, Kyung K, Walenz BP, Garcia-Hernandez E, Aguiar M, Viswanathan LD, Rogers YH, Strausberg RL, Saski CA, Lawson D, Collins FH, Kafatos FC, Christophides GK, Clifton SW, Kirkness EF, Besansky NJ: Widespread divergence between incipient Anopheles gambiae species revealed by whole genome sequences. Science. 2011, 330: 512-514.
Neafsey DE, Lawniczak MK, Park DJ, Redmond SN, Coulibaly MB, Traore SF, Sagnon N, Costantini C, Johnson C, Wiegand RC, Collins FH, Lander ES, Wirth DF, Kafatos FC, Besansky NJ, Christophides GK, Muskavitch MA: SNP genotyping defines complex gene-flow boundaries among African malaria vector mosquitoes. Science. 2011, 330: 514-517.
Favia G, Lanfrancotti A, Spanos L, Siden-Kiamos I, Louis C: Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 19-23. 10.1046/j.1365-2583.2001.00236.x.
Wilkins EE, Howell PI, Benedict MQ: IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis. Malar J. 2006, 5: 125-10.1186/1475-2875-5-125.
Favia G, della Torre A, Bagayoko M, Lanfrancotti A, Sagnon N, Toure YT, Coluzzi M: Molecular identification of sympatric chromosomal forms of Anopheles gambiae and further evidence of their reproductive isolation. Insect Mol Biol. 1997, 6: 377-383. 10.1046/j.1365-2583.1997.00189.x.
Fanello C, Santolamazza F, della Torre A: Simultaneous identification of species and molecular forms of the Anopheles gambiae complex by PCR-RFLP. Med Vet Entomol. 2002, 16: 461-464. 10.1046/j.1365-2915.2002.00393.x.
Santolamazza F, della Torre A, Caccone A: Short report: A new polymerase chain reaction-restriction fragment length polymorphism method to identify Anopheles arabiensis from An. gambiae and its two molecular forms from degraded DNA templates or museum samples. Am J Trop Med Hyg. 2004, 70: 604-606.
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-10.1186/1475-2875-7-163.
Caputo B, Nwakanma D, Jawara M, Adiamoh M, Dia I, Konate L, Petrarca V, Conway DJ, della Torre A: Anopheles gambiae complex along The Gambia river, with particular reference to the molecular forms of An. gambiae s.s. Malar J. 2008, 7: 182-10.1186/1475-2875-7-182.
Oliveira E, Salgueiro P, Palsson K, Vicente JL, Arez AP, Jaenson TG, Caccone A, Pinto J: High levels of hybridization between molecular forms of Anopheles gambiae from Guinea Bissau. J Med Entomol. 2008, 45: 1057-1063. 10.1603/0022-2585(2008)45[1057:HLOHBM]2.0.CO;2.
Caputo B, Santolamazza F, Vicente JL, Nwakanma DC, Jawara M, Palsson K, Jaenson T, White BJ, Mancini E, Petrarca V, Conway DJ, Besansky NJ, Pinto J, della Torre A: The "Far-West" of Anopheles gambiae Molecular Forms. PLoS One. 2011, 6: e16415-10.1371/journal.pone.0016415.
Carr IM, Robinson JI, Dimitriou R, Markham AF, Morgan AW, Bonthron DT: Inferring relative proportions of DNA variants from sequencing electropherograms. Bioinformatics. 2009, 25: 3244-3250. 10.1093/bioinformatics/btp583.
Black WCt, Vontas JG: Affordable assays for genotyping single nucleotide polymorphisms in insects. Insect Mol Biol. 2007, 16: 377-387. 10.1111/j.1365-2583.2007.00736.x.
Riehle MM, Guelbeogo WM, Gneme A, Eiglmeier K, Holm I, Bischoff E, Garnier T, Snyder GM, Li X, Markianos K, Sagnon N, Vernick KD: A cryptic subgroup of Anopheles gambiae is highly susceptible to human malaria parasites. Science. 2011, 331: 596-598. 10.1126/science.1196759.
Turner TL, Hahn MW, Nuzhdin SV: Genomic islands of speciation in Anopheles gambiae. PLoS Biol. 2005, 3: e285-10.1371/journal.pbio.0030285.
Stump AD, Fitzpatrick MC, Lobo NF, Traore S, Sagnon N, Costantini C, Collins FH, Besansky NJ: Centromere-proximal differentiation and speciation in Anopheles gambiae. Proc Natl Acad Sci USA. 2005, 102: 15930-15935. 10.1073/pnas.0508161102.
White BJ, Lawniczak MK, Cheng C, Coulibaly MB, Wilson MD, Sagnon N, Costantini C, Simard F, Christophides GK, Besansky NJ: Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium. Proc Natl Acad Sci USA. 2011, 108: 244-249. 10.1073/pnas.1013648108.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Gentile G, Della Torre A, Maegga B, Powell JR, Caccone A: Genetic differentiation in the African malaria vector, Anopheles gambiae s.s., and the problem of taxonomic status. Genetics. 2002, 161: 1561-1578.
Calzetta M, Santolamazza F, Carrara GC, Cani PJ, Fortes F, Di Deco MA, della Torre A, Petrarca V: Distribution and chromosomal characterization of the Anopheles gambiae complex in Angola. Am J Trop Med Hyg. 2008, 78: 169-175.
Kristan M, Fleischmann H, della Torre A, Stich A, Curtis CF: Pyrethroid resistance/susceptibility and differential urban/rural distribution of Anopheles arabiensis and An. gambiae s.s. malaria vectors in Nigeria and Ghana. Med Vet Entomol. 2003, 17: 326-332. 10.1046/j.1365-2915.2003.00449.x.
Anges Y, Thibaud M, Gil P, Fabrice C, Alex A, Luc D, Roch D, Rock A, Michel B, Isabelle G, Martin A: Cotton pest management practices and the selection of pyrethroid resistance in Anopheles gambiae population in Northern Benin. Parasit Vectors. 2011, 4: 60-10.1186/1756-3305-4-60.
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: 181-188.
Awolola TS, Oyewole IO, Amajoh CN, Idowu ET, Ajayi MB, Oduola A, Manafa OU, Ibrahim K, Koekemoer LL, Coetzee M: Distribution of the molecular forms of Anopheles gambiae and pyrethroid knock down resistance gene in Nigeria. Acta Trop. 2005, 95: 204-209. 10.1016/j.actatropica.2005.06.002.
Dery DB, Brown C, Asante KP, Adams M, Dosoo D, Amenga-Etego S, Wilson M, Chandramohan D, Greenwood B, Owusu-Agyei S: Patterns and seasonality of malaria transmission in the forest-savannah transitional zones of Ghana. Malar J. 2010, 9: 314-10.1186/1475-2875-9-314.
Djogbenou L, Pasteur N, Bio-Bangana S, Baldet T, Irish SR, Akogbeto M, Weill M, Chandre F: Malaria vectors in the Republic of Benin: distribution of species and molecular forms of the Anopheles gambiae complex. Acta Trop. 2010, 114: 116-122. 10.1016/j.actatropica.2010.02.001.
Koekemoer LL, Spillings BL, Christian RN, Lo TC, Kaiser ML, Norton RA, Oliver SV, Choi KS, Brooke BD, Hunt RH, Coetzee M: Multiple Insecticide Resistance in Anopheles gambiae (Diptera: Culicidae) from Pointe Noire, Republic of the Congo. Vector Borne Zoonotic Dis. 2011, 0: 1-8.
Tripet F, Dolo G, Lanzaro GC: Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects. Genetics. 2005, 169: 313-324.
Vezenegho SB, Brooke BD, Hunt RH, Coetzee M, Koekemoer LL: Malaria vector composition and insecticide susceptibility status in Guinea Conakry, West Africa. Med Vet Entomol. 2009, 23: 326-334. 10.1111/j.1365-2915.2009.00840.x.
Yadouleton AW, Padonou G, Asidi A, Moiroux N, Bio-Banganna S, Corbel V, N'Guessan R, Gbenou D, Yacoubou I, Gazard K, Akogbeto MC: Insecticide resistance status in Anopheles gambiae in southern Benin. Malar J. 2010, 9: 83-10.1186/1475-2875-9-83.
Boulesteix M, Simard F, Antonio-Nkondjio C, Awono-Ambene HP, Fontenille D, Biemont C: Insertion polymorphism of transposable elements and population structure of Anopheles gambiae M and S molecular forms in Cameroon. Mol Ecol. 2007, 16: 441-452.
Sharp BL, Ridl FC, Govender D, Kuklinski J, Kleinschmidt I: Malaria vector control by indoor residual insecticide spraying on the tropical island of Bioko, Equatorial Guinea. Malar J. 2007, 6: 52-10.1186/1475-2875-6-52.
Tchuinkam T, Simard F, Lele-Defo E, Tene-Fossog B, Tateng-Ngouateu A, Antonio-Nkondjio C, Mpoame M, Toto JC, Njine T, Fontenille , Awono-Ambene HP: Bionomics of Anopheline species and malaria transmission dynamics along an altitudinal transect in Western Cameroon. BMC Infect Dis. 2010, 10: 119-
Tripet F, Wright J, Cornel A, Fofana A, McAbee R, Meneses C, Reimer L, Slotman M, Thiemann T, Dolo G, Traore S, Lanzaro G: Longitudinal survey of knockdown resistance to pyrethroid (kdr) in Mali, West Africa, and evidence of its emergence in the Bamako form of Anopheles gambiae s.s. Am J Trop Med Hyg. 2007, 76: 81-87.
Tanga MC, Ngundu WI: Ecological transition from natural forest to tea plantations: effect on the dynamics of malaria vectors in the highlands of Cameroon. Trans R Soc Trop Med Hyg. 2010, 104: 659-668. 10.1016/j.trstmh.2010.07.009.
Mourou JR, Coffinet T, Jarjaval F, Pradines B, Amalvict R, Rogier C, Kombila M, Pages F: Malaria transmission and insecticide resistance of Anopheles gambiae in Libreville and Port-Gentil, Gabon. Malar J. 2010, 9: 321-10.1186/1475-2875-9-321.
Horton AA, Lee Y, Coulibaly CA, Rashbrook VK, Cornel AJ, Lanzaro GC, Luckhart S: Identification of three single nucleotide polymorphisms in Anopheles gambiae immune signaling genes that are associated with natural Plasmodium falciparum infection. Malar J. 2010, 9: 160-10.1186/1475-2875-9-160.
Aguilar R, Simard F, Kamdem C, Shields T, Glass GE, Garver LS, Dimopoulos G: Genome-wide analysis of transcriptomic divergence between laboratory colony and field Anopheles gambiae mosquitoes of the M and S molecular forms. Insect Mol Biol. 2010, 19: 695-705. 10.1111/j.1365-2583.2010.01031.x.
Ranson H, Abdallah H, Badolo A, Guelbeogo WM, Kerah-Hinzoumbe C, Yangalbe-Kalnone E, Sagnon N, Simard F, Coetzee M: 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-10.1186/1475-2875-8-299.
Nwane P, Etang J, Chouaibou M, Toto JC, Kerah-Hinzoumbe C, Mimpfoundi R, Awono-Ambene HP, Simard F: 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-10.1186/1471-2334-9-163.
Lee Y, Cornel AJ, Meneses CR, Fofana A, Andrianarivo AG, McAbee RD, Fondjo E, Traore SF, Lanzaro GC: Ecological and genetic relationships of the Forest-M form among chromosomal and molecular forms of the malaria vector Anopheles gambiae sensu stricto. Malar J. 2009, 8: 75-10.1186/1475-2875-8-75.
Etang J, Vicente JL, Nwane P, Chouaibou M, Morlais I, Do Rosario VE, Simard F, Awono-Ambene P, Toto JC, Pinto J: 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-3086. 10.1111/j.1365-294X.2009.04256.x.
Diabate A, Dao A, Yaro AS, Adamou A, Gonzalez R, Manoukis NC, Traore SF, Gwadz RW, Lehmann T: Spatial swarm segregation and reproductive isolation between the molecular forms of Anopheles gambiae. Proc Biol Sci. 2009, 276: 4215-4222. 10.1098/rspb.2009.1167.
Ridl FC, Bass C, Torrez M, Govender D, Ramdeen V, Yellot L, Edu AE, Schwabe C, Mohloai P, Maharaj R, Kleinschmidt I: A pre-intervention study of malaria vector abundance in Rio Muni, Equatorial Guinea: their role in malaria transmission and the incidence of insecticide resistance alleles. Malar J. 2008, 7: 194-10.1186/1475-2875-7-194.
Reimer L, Fondjo E, Patchoke S, Diallo B, Lee Y, Ng A, Ndjemai HM, Atangana J, Traore SF, Lanzaro G, Cornel AJ: Relationship between kdr mutation and resistance to pyrethroid and DDT insecticides in natural populations of Anopheles gambiae. J Med Entomol. 2008, 45: 260-266. 10.1603/0022-2585(2008)45[260:RBKMAR]2.0.CO;2.
Diabate A, Dabire RK, Heidenberger K, Crawford J, Lamp WO, Culler LE, Lehmann T: Evidence for divergent selection between the molecular forms of Anopheles gambiae: role of predation. BMC Evol Biol. 2008, 8: 5-10.1186/1471-2148-8-5.
Turner TL, Hahn MW: Locus- and population-specific selection and differentiation between incipient species of Anopheles gambiae. Mol Biol Evol. 2007, 24: 2132-2138. 10.1093/molbev/msm143.
Slotman MA, Tripet F, Cornel AJ, Meneses CR, Lee Y, Reimer LJ, Thiemann TC, Fondjo E, Fofana A, Traore SF, Lanzaro GC: Evidence for subdivision within the M molecular form of Anopheles gambiae. Mol Ecol. 2007, 16: 639-649.
Coulibaly MB, Pombi M, Caputo B, Nwakanma D, Jawara M, Konate L, Dia I, Fofana A, Kern M, Simard F, Conway DJ, Petrarca V, della Torre A, Traore S, Besansky NJ: PCR-based karyotyping of Anopheles gambiae inversion 2Rj identifies the BAMAKO chromosomal form. Malar J. 2007, 6: 133-10.1186/1475-2875-6-133.
Caputo B, Dani FR, Horne GL, N'Fale S, Diabate A, Turillazzi S, Coluzzi M, Costantini C, Priestman AA, Petrarca V, della Torre A: Comparative analysis of epicuticular lipid profiles of sympatric and allopatric field populations of Anopheles gambiae s.s. molecular forms and An. arabiensis from Burkina Faso (West Africa). Insect Biochem Mol Biol. 2007, 37: 389-398. 10.1016/j.ibmb.2007.01.002.
Etang J, Fondjo E, Chandre F, Morlais I, Brengues C, Nwane P, Chouaibou M, Ndjemai H, Simard F: First report of knockdown mutations in the malaria vector Anopheles gambiae from Cameroon. Am J Trop Med Hyg. 2006, 74: 795-797.
Cuamba N, Choi KS, Townson H: Malaria vectors in Angola: distribution of species and molecular forms of the Anopheles gambiae complex, their pyrethroid insecticide knockdown resistance (kdr) status and Plasmodium falciparum sporozoite rates. Malar J. 2006, 5: 2-10.1186/1475-2875-5-2.
Reimer LJ, Tripet F, Slotman M, Spielman A, Fondjo E, Lanzaro GC: An unusual distribution of the kdr gene among populations of Anopheles gambiae on the island of Bioko, Equatorial Guinea. Insect Mol Biol. 2005, 14: 683-688. 10.1111/j.1365-2583.2005.00599.x.
Diabate A, Dabire RK, Kim EH, Dalton R, Millogo N, Baldet T, Simard F, Gimnig JE, Hawley WA, Lehmann T: Larval development of the molecular forms of Anopheles gambiae (Diptera: Culicidae) in different habitats: a transplantation experiment. J Med Entomol. 2005, 42: 548-553. 10.1603/0022-2585(2005)042[0548:LDOTMF]2.0.CO;2.
Kaiser ML, Koekemoer LL, Coetzee M, Hunt RH, Brooke BD: Staggered larval time-to-hatch and insecticide resistance in the major malaria vector Anopheles gambiae S form. Malar J. 2010, 9: 360-10.1186/1475-2875-9-360.
Choi KS, Spillings BL, Coetzee M, Hunt RH, Koekemoer LL: A comparison of DNA sequencing and the hydrolysis probe analysis (TaqMan assay) for knockdown resistance (kdr) mutations in Anopheles gambiae from the Republic of the Congo. Malar J. 2010, 9: 278-10.1186/1475-2875-9-278.
Chandre F, Dabire RK, Hougard JM, Djogbenou LS, Irish SR, Rowland M, N'Guessan R: Field efficacy of pyrethroid treated plastic sheeting (durable lining) in combination with long lasting insecticidal nets against malaria vectors. Parasit Vectors. 2010, 3: 65-10.1186/1756-3305-3-65.
Ahoua Alou LP, Koffi AA, Adja MA, Tia E, Kouassi PK, Kone M, Chandre F: Distribution of ace-1R and resistance to carbamates and organophosphates in Anopheles gambiae s.s. populations from Cote d'Ivoire. Malar J. 2010, 9: 167-10.1186/1475-2875-9-167.
Ndiath MO, Brengues C, Konate L, Sokhna C, Boudin C, Trape JF, Fontenille D: Dynamics of transmission of Plasmodium falciparum by Anopheles arabiensis and the molecular forms M and S of Anopheles gambiae in Dielmo, Senegal. Malar J. 2008, 7: 136-10.1186/1475-2875-7-136.
Moreno M, Vicente JL, Cano J, Berzosa PJ, de Lucio A, Nzambo S, Bobuakasi L, Buatiche JN, Ondo M, Micha F, Do Rosario VE, Pinto J, Benito A: Knockdown resistance mutations (kdr) and insecticide susceptibility to DDT and pyrethroids in Anopheles gambiae from Equatorial Guinea. Trop Med Int Health. 2008, 13: 430-433. 10.1111/j.1365-3156.2008.02010.x.
Dabire KR, Diabate A, Djogbenou L, Ouari A, N'Guessan R, Ouedraogo JB, Hougard JM, Chandre F, Baldet T: Dynamics of multiple insecticide resistance in the malaria vector Anopheles gambiae in a rice growing area in South-Western Burkina Faso. Malar J. 2008, 7: 188-10.1186/1475-2875-7-188.
Cohuet A, Krishnakumar S, Simard F, Morlais I, Koutsos A, Fontenille D, Mindrinos M, Kafatos FC: SNP discovery and molecular evolution in Anopheles gambiae, with special emphasis on innate immune system. BMC Genomics. 2008, 9: 227-10.1186/1471-2164-9-227.
Kleinschmidt I, Torrez M, Schwabe C, Benavente L, Seocharan I, Jituboh D, Nseng G, Sharp B: Factors influencing the effectiveness of malaria control in Bioko Island, equatorial Guinea. Am J Trop Med Hyg. 2007, 76: 1027-1032.
Yaro AS, Dao A, Adamou A, Crawford JE, Traore SF, Toure AM, Gwadz R, Lehmann T: Reproductive output of female Anopheles gambiae (Diptera: Culicidae): comparison of molecular forms. J Med Entomol. 2006, 43: 833-839. 10.1603/0022-2585(2006)43[833:ROOFAG]2.0.CO;2.
Le Goff G, Leong Pock Tsy JM, Robert V: Molecular characterization of the malaria vector Anopheles gambiae s.s. in Madagascar. Med Vet Entomol. 2006, 20: 259-260. 10.1111/j.1365-2915.2006.00628.x.
Diabate A, Chandre F, Rowland M, N'Guessan R, Duchon S, Dabire KR, Hougard JM: The indoor use of plastic sheeting pre-impregnated with insecticide for control of malaria vectors. Trop Med Int Health. 2006, 11: 597-603. 10.1111/j.1365-3156.2006.01605.x.
Dabire RK, Diabate A, Baldet T, Pare-Toe L, Guiguemde RT, Ouedraogo JB, Skovmand O: Personal protection of long lasting insecticide-treated nets in areas of Anopheles gambiae s.s. resistance to pyrethroids. Malar J. 2006, 5: 12-10.1186/1475-2875-5-12.
Tripet F, Thiemann T, Lanzaro GC: Effect of seminal fluids in mating between M and S forms of Anopheles gambiae. J Med Entomol. 2005, 42: 596-603. 10.1603/0022-2585(2005)042[0596:EOSFIM]2.0.CO;2.
Czeher C, Labbo R, Arzika I, Duchemin JB: Evidence of increasing Leu-Phe knockdown resistance mutation in Anopheles gambiae from Niger following a nationwide long-lasting insecticide-treated nets implementation. Malar J. 2008, 7: 189-10.1186/1475-2875-7-189.
Acknowledgements and funding
We are grateful to all scientists and entomology teams who provided samples utilized in this study; we especially thank: K. Adasi, T. Baldet, G. Carrara, P.J. Cani, D.J. Conway, C. Curtis, F. Fortes, M. Jawara, T.G.T. Jaenson, A. Mendjibe, D.C. Nwakanma, K. Palsson and Y.T. Touré. The project was supported by INFRAVEC project (FP7-INFRASTRUCTURES-2008: Grant Agreement Number 228421) and by SAPIENZA University (Progetto Ateneo 2009). FS was partially supported by Compagnia di San Paolo (Torino, Italy) in the context of and C.I.R.M. - Italian Malaria Network. BC was partially supported by I Faculty of Medicine, SAPIENZA University; EM was supported by "Ateneo Federato delle Scienze delle Politiche Pubbliche e Sanitarie" (SAPIENZA University) and C.I.R.M. - Italian Malaria Network.
The authors declare that they have no competing interests.
FS, MC and JLV carried out the molecular genotyping and participated in the sequence alignment. FS, BC and EM performed the statistical analyses. VP and JP participated in the design of the study and contributed to data analyses. AdT conceived of the study, participated in its design and coordination, and drafted the manuscript. FS, BC, MC, EM, VP and JP contributed to the manuscript finalization. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: P values of pairwise comparisons of CNP581and CNP690scores The data show P values of pairwise comparisons of CNP581 and CNP690 scores calculated by QSV analyser in Anopheles gambiae specimens classified by IGS581/IGS690 PCR-RFLPs. (DOC 50 KB)
Additional file 3: Frequency of M and S molecular forms and putative MS hybrids in 10/56 papers on Anopheles gambiae s.s. published since 2006. (DOC 40 KB)
About this article
Cite this article
Santolamazza, F., Caputo, B., Calzetta, M. et al. Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiae molecular form identification. Malar J 10, 215 (2011). https://doi.org/10.1186/1475-2875-10-215
- Internal Transcribe Spacer
- Molecular Form
- Incipient Species
- Afrotropical Region
- Secondary Contact Zone