Molecular characterization of the malaria vector Anopheles barbirostris van der Wulp in Sri Lanka
© Gajapathy et al.; licensee BioMed Central Ltd. 2014
Received: 10 April 2014
Accepted: 2 July 2014
Published: 29 July 2014
Anopheles barbirostris is a vector of malaria in Sri Lanka. The taxon exists as a species complex in the Southeast Asian region. Previous studies using molecular markers suggest that there are more than 4 distinct clades within the An. barbirostris complex in Southeast Asia. The present study characterizes Sri Lankan An. barbirostris using mtDNA cytochrome oxidase subunit I (COI) and ribosomal RNA internal transcribed spacer 2 (ITS2) gene sequences.
DNA was extracted from morphologically identified An. barbirostris specimens from Sri Lanka, the COI and ITS2 regions amplified and their sequences analysed by comparison with other GenBank entries. Maximum likelihood trees suggested that Sri Lankan An. barbirostris constitute a different molecular type most closely related to clade I.
Considering the uncorrected p distances between the clade I and Sri Lankan specimens it is fair to assume that the specimens collected from widely separated locations in Sri Lanka with morphology characteristic of An. barbirostris s.l. form a new molecular type with close resemblance to An. barbirostris s.s from Indonesia and Thailand.
Anopheles barbirostris van der Wulp is a vector of malaria in Sri Lanka , India and Southeast Asia, and a vector of Brugian filariasis in Southeast Asia . Anopheles barbirostris and 12 related species form the taxonomically complex Barbirostris Group of malaria vectors in this region ( and references therein). An. barbirostris s.l. is a species complex with up to five sibling species whose individual distributions are not well established [4–7]. Members of this complex are typically zoophillic, develop in still or slow moving fresh water such as in swamps and rice fields, and are relatively tolerant of organic pollution ( and references therein). Larvae were recently also found in brackish water (salinity of 4–15 parts per thousand) in domestic wells and mangrove swamps in the Jaffna district of Northern Sri Lanka . Morphological differentiation of anopheline sibling species is difficult and has previously led to significant misidentification with implications for vector control, e.g. mosquitoes long identified morphologically as Anopheles subpictus sibling species B in Sri Lanka were recently shown to be Anopheles sundaicus s.l. by analysis of ribosomal RNA gene sequences .
Other members of the Barbirostris Group, e.g. An. reidi and An. barbrumbosus, have been identified in Sri Lanka  but the five other species of the Barbirostris Subgroup, viz. An. campestris, An. donaldi, An. franciscoi, An. hodgkini and An. pollicaris have not been reported in Sri Lanka. The present study was designed to characterize morphologically identified An. barbirostris s.l. in Sri Lanka using DNA sequences of the cytochrome oxidase subunit I (COI) of mtDNA and internal transcribed spacer 2 (ITS2) of ribosomal RNA in relation to previous DNA sequence data from An. barbirostris s.l. collected elsewhere in Southeast Asia .
DNA extraction and PCR amplification of COI and ITS2 regions
Among 27 morphologically identified An. barbirostris s.l., DNA was extracted from two, one and two specimens collected from Hambantota, Nainadeevu and Puttalam respectively using previously reported methods . The COI region was PCR amplified with forward 5′ TTG ATT TTT TGG TCA TCC AGA AGT 3′ and reverse 5′ TAG AGC TTA AAT TCA TTG CAC TAA TC 3′primers . Each PCR (50 μL total volume) reaction mixture contained 1x Thermo Scientific Ready Mix (1.5 mM Mg Cl2, 0.625U Taq DNA polymerase, 0.2 mM dNTPs and 1x Taq buffer), 2-4 μl of template DNA and 10pmoles of each primer. The thermal profile for PCR was 94°C for 5 min, followed by 35 cycles of 95°C for 40 s, 50°C for1 min and72°C for1 min. A final extension temperature of 72°C was set for 10 min.
The ITS2 region was amplified using forward 5′ ATC ACT CGG CTC ATG GAT CG 3′ and reverse 5′ ATG CTT AAA TTT AGG GGG TAG TC 3′ primers . Fifty μl reaction volumes were prepared with 1x reaction buffer, 2.5 mM MgCl2; 0.2 mM dNTPs, 20pmoles of each primer, 1.0 U of Taq polymerase and 5 μl of template DNA. The PCR thermal profile was 94°C for 5 min, followed by 35 cycles of 94°C for 1 min, 55°C for 2 min, 72°C for 2 min and a final extension at 72°C for 10 min. The amplified products were separated by agarose gel electrophoresis, stained with ethidium bromide and sequenced at Macrogen, Amsterdam, Netherlands.
Sequences were edited in Finch TV (Geospiza, Seattle, USA) and aligned with Clustal W (MEGA 5.1) . A phylogenetic tree of the sequences was created using the maximum likelihood approach (Phyml 3.0) . A HKY substitution model was selected with a gamma distribution based on the lower Bayesian criterion index (MEGA 5.1) . Genbank deposited sequences from Thailand (16 for COI and 28 for ITS2) and Indonesia (08 samples for COI and 06 samples for ITS2) were also used in the tree construction [5, 6].
Anopheles barbirostris specimens were identified based on standard morphological characteristics . The presence of a median tuft was observed in all flies. Median scales in the abdominal sternum were also observed. Maxillary palps were dark-scaled and the wings had a narrow apical fringe spot at the vein R4+5.
Anopheles barbirostris were previously classified into five distinct clades by Paredes-Esquivel et al. , based on COI and ITS2 data. Clade IV was identified as a species with intermediate morphology between An. barbirostris and An. campestris. Later, clade I was identified unequivocally as An. barbirostris s.s. and Clade II was described as a new species named An. vanderwulpi. Saeung et al.  recognized three sibling species in the An. barbirostris complex based on karyotyping and COI, COII and ITS2 sequences. The sibling species designated A1 and A2 by Saeung et al.  correspond to clades III and IV, respectively, of Paredes-Esquivel et al. . Differential susceptibility of the members of the Barbirostris Complex to infection with malaria parasites has been established in laboratory reared colonies in Thailand .
The disparity between mitochondrial (COI) and nuclear gene (ITS2) trees in the present analysis may be explained in part by different evolutionary rates. Mitochondrial DNA typically mutates at a higher rate than nuclear sequences  such as the ITS2, and thus the latter may be more appropriate for analysing relationships between species.
Considering the uncorrected p distances of 4.5% (COI) and 5.1% (ITS2) between the clade I and Sri Lankan specimens it is fair to assume that the specimens collected from widely separated locations in Sri Lanka with morphology characteristic of An. barbirostris s.l. form a new molecular type with close resemblance to the clade I (An. barbirostris s.s) of Paredes-Esquivel et al. .
The authors acknowledge the Director, National Anti-malaria Campaign, Sri Lanka for permission to collect mosquitoes and members of the Entomology Team of the Anti-Malaria Campaign of the district of Hambantota for field assistance. The study was partly supported by grant NSF/RG/BT/2011/04 from the National Science Foundation of Sri Lanka.
- Amerasinghe PH, Amerasinghe FP, Kondarsen FFP, Fonseka KT, Wirtz RA: Malaria vectors in a traditional dry zone village in Sri Lanka. Am J Trop Med Hyg. 1999, 60 (3): 421-429.PubMedGoogle Scholar
- Lien JC, Kawengian BA, Partono F, Lami B, Cross JH: A brief survey of the mosquitoes of South Sulawesi, Indonesia, with special reference to the identity of Anopheles barbirostris (Diptera: Culicidae) from the Margolembo area. J Med Entomol. 1977, 13: 719-727.View ArticlePubMedGoogle Scholar
- Sinka ME, Bangs MJ, Manguin S, Chareonviriyaphap T, Patil AP, Temperley WH, Gething PW, Elyazarn IRF, Kabaria CW, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2011, 4: 8.View ArticleGoogle Scholar
- Saeung A, Otsuka Y, Baimai V, Somboon P, Pitaswat B, Tuetun B, Junkum A, Takaoka H, Choochote W: Cytogenetic and molecular evidence for two species in the Anopheles barbirostris complex (Diptera: Culicidae) in Thailand. Parasitol Res. 2007, 101: 1337-1344.View ArticlePubMedGoogle Scholar
- Saeung A, Baimai V, Otsuka Y, Rattanarithikul R, Somboon P, Junkum A, Tuetun B, Takaoka H, Choochote W: Molecular and cytogenetic evidence of three sibling species of the Anopheles barbirostris Form A (Diptera: Culicidae) in Thailand. Parasitol Res. 2008, 102: 499-507.View ArticlePubMedGoogle Scholar
- Paredes-Esquivel C, Donnelly MJ, Harbach RE, Townson H: A molecular phylogeny of mosquitoes in the Anopheles barbirostris subgroup reveals cryptic species: implications for identification of disease vectors. Mol Phylogenet Evol. 2009, 50 (1): 141-151.View ArticlePubMedGoogle Scholar
- Jariyapan N, Baimai V, Poovorawan Y, Roytrakul S, Saeung A, Thongsahuan S, Suwannamit S, Otsuka Y, Choochote W: Analysis of female salaivary gland proteins of the Anopheles barbirostris complex (Diptera: Culicidae) in Thailand. Parasitol Res. 2010, 107: 509-516.View ArticlePubMedGoogle Scholar
- Jude PJ, Tharmatha T, Gopika S, Meena S, Kannathasan S, Raveendran S, Ramasamy R, Surendran SN: Salinity-tolerant larvae of mosquito vectors in the tropical coast of Jaffna. Sri Lanka and the effect of salinity on the toxicity of Bacillus thuringiensis to Aedes aegypti larvae. Parasit Vectors. 2012, 5: 269.PubMed CentralView ArticlePubMedGoogle Scholar
- Surendran SN, Singh OP, Jude PJ, Ramasamy R: Genetic evidence for malaria vectors of the Anopheles sundaicus complex in Sri Lanka with morphological characteristics attributed to Anopheles subpictus species B. Malar J. 2010, 9: 343.PubMed CentralPubMedGoogle Scholar
- Harrison BA, Reinert JF, Sirivanakarn S, Huang YM, Peyton EL, de Meillon B: Distributional and biological notes on mosquitoes from Sri Lanka (Ceylon) (Diptera: Culicidae). Mosq Syst. 1974, 6: 142-162.Google Scholar
- Reid JA: The Anopheles barbirostris group (Diptera, Culicidae). Bull Entomol Res. 1962, 53 (1): 1-57.View ArticleGoogle Scholar
- Amerasinghe FP: A guide to the identification of the anopheline mosquitoes (Diptera: Culicidae) of Sri Lanka. Cey J Sci (Biol Sci). 1992, 22: 1-13.Google Scholar
- Gajapathy K, Peiris LBS, Goodacre SL, Silva A, Jude PJ, Surendran SN: Molecular identification of potential leishmaniasis vector species within the Phlebotomus (Euphlebotomus) argentipes species complex in Sri Lanka. Parasit Vectors. 2013, 6: 302.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony method. Mol Biol Evol. 2011, 28: 2731-2739.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2008, 52: 696-704.View ArticleGoogle Scholar
- Townson H, Dyer N, McAlister E, Satoto TBT, Boewono DT, Bangs MJ, Harbach RE: Systematics of Anopheles barbirostris van der Wulp and a sibling species of the Barbirostris Complex (Diptera: Culicidae) in eastern Java, Indonesia. Syst Entomol. 2013, 38: 180-191.View ArticleGoogle Scholar
- Thongsahuan S, Baimai V, Junkum A, Saeung A, Min GS, Joshi D, Park MH, Somboon P, Suwonkerd W, Tippawangkosol P, Jariyapan N, Choochote W: Susceptibility of Anopheles campestris-like and Anopheles barbirostris species complexes to Plasmodium falciparum and Plasmodium vivax in Thailand. Mem Inst Oswaldo Cruz. 2011, 106: 105-112.View ArticlePubMedGoogle Scholar
- Brown WM, George M, Wilson AC: Rapid evolution of animal mitochondrial DNA. Genetics. 1979, 76: 1967-1971.Google 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.