Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiaemolecular form identification

Background: 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. Results: 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. Conclusions: 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.