|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
The Anopheles funestus and An. minimus groups belong to the Myzomyia series, subgenus Cellia,1 and are closely related based on various cytogenetic and molecular studies (Garros C and others, unpublished data).2,3 They are probably considered distinct only because of their geographic separation.4 In spite of this genetic proximity, both groups had seldom been studied jointly (Garros C and others, unpublished data).2,3
On the Afrotropical continent, the An. funestus group consists of nine species, namely An. brucei, An. confusus, An. fuscivenosus, An. leesoni, and An. rivulorum, as well as the An. funestus subgroup, which is composed of An. aruni, An. funestus, An. parensis, and An. vaneedeni.57 Although An. leesoni does not belong to the An. funestus group, following the recommendation of Harbach,1 who included it in the An. minimus group, it was included here within the members of the An. funestus group because of its morphologic similarities and sub-Saharan distribution. Beside the anthropophilic and endophilic An. funestus, which is a highly efficient malaria vector, the other species of the group are mainly zoophilic and play little or no role in malaria transmission.5,6,10 Realizing the difficulties of morphologic identification and the need to elucidate the role of individual species in malaria transmission, molecular methods have been introduced for distinguishing members of the An. funestus group.8,9,11,12 A polymerase chain reaction (PCR)-single strand conformation polymorphism (SSCP) assay can identified four species: An. funestus, An. leesoni, An. rivulorum, and An. vaneedeni.8 The inclusion of An. parensis revealed overlapping banding patterns with An. vaneedeni. A sequence length variation in the internal transcribed spacer 2 (ITS2) of ribosomal DNA was identified by its amplification and permitted the identification of An. parensis.11 Very recently, a species-specific PCR assay was developed for differentiating these five members of the group.12
Similarly, in the Oriental region, the An. minimus group is composed of nine species: An. aconitus, An. filipinae, An. flavirostris, An. fluviatilis, An. mangyanus, An. minimus species A and C, belonging to the An. minimus complex,1 An. pampanai, and An. varuna.4 Anopheles minimus s.l. is most commonly associated with rivers and streams throughout Southeast Asia, where it is one of the primary vectors of human malaria. These species are difficult to distinguish at the adult stage due to overlapping characteristics.4 However, associated immature skins can be useful for differentiating the species, except those of the An. minimus complex. These species play different roles in malaria transmission, requiring vector and non-vector species to be differentiated with reliable markers. Molecular methods were recently developed for distinguishing different species of the group. Two PCR-based techniques, an allele-specific amplification (ASA) and an SSCP were elaborated for distinguishing, respectively An. minimus A from C; and both species A and C along with An. aconitus and An. varuna.13 A PCR-restriction fragment length polymorphism (RFLP) method was also designed for the identification of these four species, as well as for An. pampanai, An. culicifacies B, An. jeyporiensis, and the hybrids between An. minimus species A and C.14 Very recently, a single multiplex PCR assay, using sequence characterized amplified region (SCAR) markers derived from individual random amplified polymorphic DNA markers, was developed for an easy and reliable identification of An. minimus A and C, and their hybrids, as well as the three closely related species of the An. minimus group (An. aconitus, An. varuna, and An. pampanai).15,16
Since the An. funestus and An. minimus groups include major malaria vectors on both the African and Asian continents, there is a tremendous need for accurate and precise identification methods for members of these two groups. Moreover, the standardization of molecular identification methods is important for their transfer to field laboratories for faunistic studies, as well as for scientists working on phylogenetic relationships. Therefore, based on a previous study by Koekemoer and others,12 the aim of the present work was to develop a multiplex PCR assay, based on one-step amplification, targeting An. funestus- and An. minimus-specific parts of the gene encoding the ITS2 spacer, for rapid, easy, and reliable identification of the most common species of both groups.
| MATERIALS AND METHODS |
|
|
|---|
|
Amplification by PCR and sequencing of the ITS2 region. The ITS2 primers from the 5.8S and 28S coding regions flanking the variable ITS2 region were used to amplify the genomic DNA.14 The primers were as follows: ITS2A: 5'-TGT GAA CTG CAG GAC ACA T-3'; ITS2B: 5'-TAT GCT TAA ATT CAG GGG GT-3'. The PCR mixture contained 2.5 µL of 10x reaction buffer (Qiagen, Valencia, CA), 200 µM of each dNTP, 64 nmol of primer, 0.5 units of Taq DNA polymerase (Qiagen), and 3 µL of DNA template. Amplification conditions were composed of an initial hot start at 94°C for two minutes, followed by 40 cycles of denaturation at 94°C for 30 seconds, annealing at 50°C for 30 seconds, and extension at 72°C for 40 seconds, and a final extension at 72°C for 10 minutes. One negative control was included. To confirm amplifications, 7 µL of the PCR-product was subjected to electrophoresis on a 1.5% agarose gel stained with ethidium bromide. The remaining PCR product (18 µL) was loaded onto a 1.5% agarose gel and subjected to electrophoresis. The PCR product was excised from the gel and purified using QIAquick spin columns (Qiagen) and products were electrophoresed on a 1% agarose gel to quantify the purification. Double-strand PCR products were directly sequenced by MWG-Biotech (Ebersberg, Germany) with either primer ITS2A or ITS2B. At least two individuals per species were sequenced in both strands. GenBank accession numbers are as follows: An. funestus AY259143, An. parensis AY 259148, An. rivulorum AY255105, An. vaneedeni AY259147, An. leesoni AY255107, An. aconitus AY255106, An. minimus species A AY255108, An. minimus species C AY255109, An. pampanai AY255110, and An. varuna AY255111.
Primer design and allele-specific PCR.
From the sequence of individual specimens from each species, specific primers were designed and their specificity were tested. The ITS2 sequences of all 10 species were aligned using BioEdit.19 Moreover, sequences from previous studies obtained through GenBank were also used in the alignment for improving primer design.12,14 The location for primer sequences was chosen on the criteria of at least three nucleotide differences between species and PCR products easily distinguishable by agarose gel electrophoresis. The oligonucleotide primers were synthesized by MWG-Biotech. Primer names, sequences, and sizes of the PCR products are shown in Table 2
. To check the length of the amplified fragment and assess primer specificity, each primer was tested with samples of different locations of the considered species, as well as the other nine species. Amplification was done in a final volume of 25 µL containing 2.5 µL of 10x buffer (Qiagen), 200 µM of each dNTP, 64 nmol of each primer, 0.5 units of Taq polymerase, and 3 µL of DNA template. After an initial denaturation step at 94°C for two minutes, 30 cycles were programmed as follows: 94°C for 30 seconds, 45°C for 30 seconds, 72°C for 40 seconds, and a final extension at 72°C for five minutes. Amplimers were resolved by electrophoresis on a 3% agarose gel with Tris-borate-EDTA buffer and stained with ethidium bromide.
|
Validation of the assay.
To determine the reliability of the new PCR assay, a large sample of specimens from both groups and from different localities over both continents, Africa and Southeast Asia, were tested (Table 1
). Standards used for the design of the primers were run as positive controls on the electrophoresis gel.
| RESULTS |
|
|
|---|
|
Criteria for primer design was that the PCR products after amplification would be specific for each species and easily visualized on an agarose gel. The oligonucleotide sequence for each primer and the melting temperature (Tm) are shown in Table 2
. Since both groups are closely related and some species exhibited very similar ITS2 sequences, such as the An. minimus complex and An. leesoni, the MIA (An. minimus A) and MIC (An. minimus C) primers had only one base difference from the ITS2 sequence of An. leesoni. In spite of this, the species specificity of all the primers designed was conserved.
The universal forward primer (ITS2A) is located in the conserved 5.8S gene, whereas the species-specific reverse primers are within the ITS2 spacer region. Each of the 10 species-specific primers works in combination with the universal forward primer ITS2A to give a good and reliable amplification signal. Ten species of both the An. funestus and An. minimus groups, five for each group, can be identified by the combination of 11 primers in our multiplex PCR. One of the 11 primers is the forward primer (ITS2A).
Lengths of amplified species-specific products were 90 bp for An. pampanai, 180 bp for An. minimus C, 200 bp for An. aconitus, 235 bp for An. parensis, 260 bp for An. varuna, 280 bp for An. leesoni, 310 bp for An. minimus A, 400 bp for An. rivulorum, 460 bp for An. funestus, and 555 bp for An. vaneedeni (Figure 2
and Table 2
).
|
|
Validation of the multiplex PCR assay.
The multiplex PCR assay was validated on specimens previously identified by molecular methods and/or morphology (Table 1
). The DNA from the standard specimens was subjected to electrophoresis on an agarose gel as positive controls. A total of 213 specimens were tested for the validation of the PCR, including 63 specimens for the An. funestus group and 150 specimens for the An. minimus group. Two specimens of the An. minimus group from Cambodia and all 12 An. minimus C from Thailand (Kanchanaburi Province) were misidentified by morphology due to difficult field conditions for a precise identification. One specimen from Cambodia that was identified as An. minimus A was identified by the multiplex assay as An. pampanai, and second specimen that was identified as An. minimus A was identified by the multiplex assay as An. minimus C. The 12 samples from Kanchanaburi Province, which were all morphologically identified as An. aconitus, were identified by the multiplex assay as An. minimus C. This was confirmed by testing the two specimens from Cambodia and five specimens from Kanchanaburi Province with the RFLP assay of Van Bortel and others.14
| DISCUSSION |
|
|
|---|
Allele-specific strategies for discriminating closely related species involve 1) the use of universal primers able to amplify target sequences in virtually all species, and 2) species-specific primers for the identification of the target species. Our method used a universal forward primer that binds to the 5' end of ITS2 in species in the An. funestus and An. minimus groups, whereas the reverse primers are species-specific and hybridize at different positions along the ITS2. This method has the advantage that any unknown Anopheles species can be amplified by universal primers, and the amplification allows unique identification of the expected species from the An. funestus and/or An. minimus groups with the specific banding pattern. The forward and reverse primers might be so universal that they will potentially amplify the ITS2 sequence of any Diptera. Moreover, the possibility of easily amplifying the entire ITS2 plays the role of a positive control in checking the quality of the amplification. Therefore, in a one-step PCR assay, the anopheline fauna of a specific region can be screened.
The validation of the multiplex PCR showed that two specimens of An. pampanai and An. minimus C from Cambodia were misreferenced in the collection instead of a misidentification by molecular and morphologic methods. The specimens from Thailand originally identified as An. aconitus, but later identified as An. minimus C by molecular methods, illustrates the difficulty in making precise identifications within the An. minimus group using only adult morphology. The presence of An. minimus C in the province of Kanchanaburi Province has already been demonstrated.13,15
No species hybrids between members of the An. funestus group have been reported either from the field or the laboratory. This is not the case for the An. minimus group, in which rare hybrids (less than 1%) between An. minimus species A and C have been reported in some regions of Vietnam.14,15 We tested two hybrids of An. minimus A and C identified by ITS2 restriction.14 The primer MIA, which is specific for An. minimus A, binds to the ITS2 sequence of the hybrids. This was verified after sequencing the ITS2 of the hybrids (Garros C and others, unpublished data). The multiplex PCR assay reached its limit by not identifying hybrids between An. minimus A and C. However, the low number of base differences existing in the ITS2 hybrid sequence, compared with those of both species of the An. minimus complex, makes it difficult to design hybrid primers that will not bind with the two former species. The development of a multiplex PCR that could identify the 10 most common species of both groups, including the hybrids between An. minimus A and An. minimus C, is possible if a more variable locus is used. The intergenic spacer might be a better region. However, male hybrids might be impossible to identify if all ribosomal DNA repeats were located on the X chromosome, as in the case of An. gambiae.26 No information is available on the precise position of ribosomal DNA in An. minimus.
Primer design strategy was targeted at the ribosomal subunit DNA, which characterized by alternate well-conserved regions and variable ones. Molecular identification methods for anopheline species complexes have used mainly the ribosomal DNA locus.12,23,24,2733 This locus has many advantages: it is represented in multiple copies throughout the genome in mosquitoes and leads to high amplification signal. It contains highly variable regions that facilitates the selection of primer binding sites for each species to generate specific amplification products of different lengths. More precisely, both ITS1 and ITS2 regions show relatively high levels of interspecies variations, which allows us to design species-specific diagnostic assays.19,24,32
Many malaria control laboratories in Africa and Asia are already equipped to perform PCR assays; therefore, no additional equipment will be needed to identify members of both An. funestus and An. minimus groups. The multiplex PCR assay is rapid, cheap, sensitive, and easy to use, and can be easily adapted to countries in these regions. Both male and female mosquitoes and any life stage can be identified. Preservation is simple since specimens can be stored desiccated with silica gel or in ethanol. Only small quantities of material (12 legs) are needed for identification, leaving the rest of the body parts for additional analysis, such as sporozoite detection, blood meal analysis, and population genetics or insecticide resistance status.
In conclusion, the multiplex AS-PCR reported in this study represents a rapid and efficient method that is applicable on a routine basis for the identification of members of the An. funestus and An. minimus groups. This one-step PCR method constitutes a very powerful tool in large surveys of anopheline populations and widespread collections from Africa and Southeast Asia, and will assist in the improvement of the current knowledge on species distribution of both groups. Moreover, this method does not rely on skillful interpretation; therefore, no subjective bias is introduced in the identification. This is the first AS-PCR developed for identifying Anopheles species from two different continents. Malaria control programs on both continents could select the primers specific to either the An. funestus or An. minimus fauna present in their study sites. These assays will help improving our knowledge in the Arabic corridor or Indian region where some species distributions are unresolved.4,5 The improvement of the existing assays and their standardization for both continents will aid vector control research. This kind of universal AS-PCR could be applied to other groups of insects or organisms living in sympatry and involved in agricultural or medical problems.
Received July 28, 2003. Accepted for publication September 12, 2003.
Acknowledgments: This study was made possible because of the assistance of a number of people who helped with mosquito sampling on both continents. We are especially grateful to D. Fontenille (Institute of Research for Development, Montpellier, France), S. Laventure (Institut Pasteur de Madagascar, Antananarivo, Madagascar), G. Le Goff (Institute of Research for Development, Antananarivo, Madagascar), J. Mwangi and G. Otsyula (Kenya Medical Research Institute, Nairobi, Kenya), J. La Grange (Department of Health, Nelspruit, Mpumalanga, South Africa), T. Furemele (Department of Health, Tzaneen, Northern Province, South Africa), K. Hargreaves (Department of Health, Jozini, KwaZulu/Natal, South Africa), G. Kloke (Foray Consultants, Maputo, Mozambique), B. Ntomwa (Malaria Control Program, Oshakati, Namibia), T. Byembabazi (Malaria Control Program, Kampala, Uganda), and the entomology team of the National Institute of Malarialogy, Parasitology and Entomology (Hanoi, Vietnam) for assistance with field collections or supply of specimens. We also thank V. Baimai and M. Savanasoojarit (Mahidol University, Bangkok, Thailand) for specimens from Thailand.
Financial support: This work was partially supported by an EC grant, INCO-DC research project ERBIC 18CT970211, and a VIHPAL 2000 grant from the French Ministry of Research to Sylvie Manguin.
Authors addresses: Claire Garros and Sylvie Manguin, Center of Biology and Management of Populations, Campus International de Baillarguet, CS 30016, 34988 Montferrier sur Lez, France, Lizette L. Koekemoer and Maureen Coetzee, Vector Control Reference Unit, National Institute for Communicable Diseases, PO Box 1038, Johannesburg, South Africa, Marc Coosemans, Department of Parasitology, Prince Leopold Institute of Tropical Medicine, Nationalestraat 155, B-2000 Antwerp, Belgium.
Reprint requests: Claire Garros or Sylvie Manguin, Center of Biology and Management of Populations, Campus International de Baillarguet, CS 30016, 34988 Montferrier sur Lez, France, Telephone: 33-4-99-62-33-28/27, Fax: 33-4-99-62-33-45, E-mails: garros{at}mpl.ird.fr and manguin{at}mpl.ird.fr.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
M. A. Zapata, A. V. Cienfuegos, O. I. Quiros, M. L. Quinones, S. Luckhart, and M. M. Correa Discrimination of Seven Anopheles Species from San Pedro de Uraba, Antioquia, Colombia, by Polymerase Chain Reaction-Restriction Fragment Length Polymorphism Analysis of Its Sequences Am J Trop Med Hyg, July 1, 2007; 77(1): 67 - 72. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |