Identification of sibling species status of Anopheles culicifacies breeding in polluted water bodies in Trincomalee district of Sri Lanka
© Gunathilaka and Karunaraj; licensee BioMed Central. 2015
Received: 26 February 2015
Accepted: 6 May 2015
Published: 22 May 2015
Anopheles culicifacies s.l., the major vector of malaria in Sri Lanka, is known to breed in clean and clear water. However, recent findings have confirmed breeding from waste water bodies in urban and semi-urban areas. No study has been conducted to identify whether it is vector or non-vector siblings. The objective of the study was to identify the sibling species status of An. culicifacies s.l.
Anopheles culicifacies s.l. adult samples (reared from larvae) were obtained from the Padavisiripura Entomological team attached to Tropical and Environmental Diseases and Health Associates (TEDHA) Malaria Elimination Programme in Trincomalee District. The collected mosquito specimens were processed for the extraction of genomic DNA individually. The PCR amplifications were carried out using different primer combinations for differentiating species A from D, species B from C, species B from E, and species B, C, and E from each other. The results obtained from the gel electrophoresis were compared with the marker and band sizes of 359 bp, 248 bp, 95 + 248 bp, 166 + 359 bp and 178 + 248 bp were used to compare the sibling species A, B, C, D and E, respectively.
The molecular biological identification of the field caught An. culicifacies s.l. samples observed that only 13.34 % (4/30) was represented sibling species B. About 86.66 % (26/30) of the samples were An. culicifacies sibling species E. This study also provided evidence that An. culicifacies E was able to breed in a wide range of breeding habitats. This is the first time that An. culicifacies E breeding in waste water was confirmed by a molecular method. Malaria control programmes in most parts of the country focus on rural communities as a result of bio-ecology of mosquitoes. Therefore, unusual breeding habitats, such as waste water collections, may mislead the current vector controlling programmes.
These results reconfirm that An. culicifacies s.l. has adapted to breed in a wide range of water bodies, including waste water collections. Since the majority of mosquitoes sampled belong to sibling species E, this may adversely affect the current malaria elimination programme and new strategies should be adopted to control malaria vectors breeding in these unusual breeding habitats in Sri Lanka.
Anopheles culicifacies sensu lato (s.l.) has a wide distribution in India, extending to Ethiopia, Yemen, Iran, Afghanistan and Pakistan in the west and Bangladesh, Myanmar, Thailand, Laos and Vietnam in the east. It is also found in Nepal and Southern China to the north and extends to Sri Lanka in the south [1, 2]. Anopheles culicifacies s.l. is the major vector of malaria in India and contributes for transmission of 60–70 % of the 2–3 million malaria cases reported every year . It is also the best known and most important malaria vector in Sri Lanka . It has thus far been recognized as a complex of five sibling species that are provisionally designated as species A, B, C, D, and E  and the distribution of these species has been mapped .
Out of all five, only sibling B and E has detected from Sri Lanka . Sibling species complex B of An. culicifacies was considered responsible for malaria transmission in Sri Lanka . However recently, the sibling species E has been recognized as the responsible vector for malaria transmission . The sibling species exhibit distinct biological characters, host-feeding preference, biting activity, and susceptibility to commonly used insecticides in public health programmes .
Malaria control programmes in most parts of the country traditionally focus on rural communities, where the disease was prevalent. The bio-ecology of Anopheles breeding habitats in urban areas received very little attention . Recent findings identified a conducive breeding of An. culicifacies in waste water containing drains in urban/semi-urban areas in Trincomalee District of Sri Lanka observing the highest larval density for An. culicifacies .
Since, current control activities target common types of breeding habitats having clear, sunlit, fresh waters such as sand pools, rock pools of drying rivers, margins of slow moving water and irrigation channels [10, 11], unusual breeding habitats are likely to be missed.
Therefore, identification of sibling species status is very important since the ability to transmit the malaria parasite is different from one sibling species to another . Hence, accurate detection of sibling species status of malaria vectors may facilitate in developing cost effective vector controlling strategies under the current malaria elimination programme embarked in 2009 as a result of the significant reduction in malaria cases in the country.
Hence, the objective of the present study was to explore the sibling species status of An. culicifacies breeding in waste water habitats in order to adopt proper cost-effective vector control interventions through accurate vector identification.
Malaria was considered as a major health problem for people living in rural agricultural communities in Sri Lanka and it was endemic in the Dry Zone of the country. The abundance of malaria vector, mosquitoes has not been studied in some parts of the Sri Lanka, especially in the Northern and Eastern Provinces over the past 30 years in view of the security situation. According to the recent publications, there are differences in the bio-ecology of anopheline mosquitoes . Trincomalee district is an area where it showed some diversified breeding habitats for An. culicifacies s.l. The district of Trincomalee was selected for this study in order to determine the molecular identity of the main malaria vector, An. culicifacies s.l. breeding in unusual habitats.
Selecting the sampling sites
Collection of Anopheles mosquito samples from the field
Anopheline mosquitoes from the identified waste water habitats were collected for a period of 3 months from March, 2014 to May, 2014 using standard dipping method as described by WHO .
Mosquito larval specimens were collected using larval ladles [25 mm diameter, 250 ml capacity] with extensible handle. The number of dips was determined by considering the size of the breeding site. The number of dips and positive dips for Anopheles were recorded. Among the positive dips number of different larval stages were recorded (I-II and III-IV). The collected larvae were transferred to the field laboratory established at the Padavisiripura entomology office.
Identification of field caught mosquitoes
Above collected larvae were reared in net covered plastic basins (30 cm diameter) filled with one litre of well water and feed with fish food dust. These larvae were reared until the adults emerged. Emerged adults were identified to the species level by using achromatic magnifying lens (× 10) and the standard morphological taxonomic keys prepared for Sri Lankan anophelines [13, 14].
Processing mosquito samples
Anopheles culicifacies s.l. mosquitoes were separated from the adult mosquito collection obtained from larval rearing and individual specimens of An. culicifacies s.l. were separated into two parts, the head–thoracic region and the abdominal region. Both body parts were kept separately in two micro-centrifuge tubes. The processed samples were transferred to the Biotechnology Unit of Industrial Technology Institute, Colombo. The specimens were stored at −20 °C freezer until the DNA extraction.
Extraction of genomic DNA
The genomic DNA was extracted from thirty specimens from both head-thoracic and abdominal regions of individual mosquito separately by phenol-chloroform method .
Sequences of species specific primers used for multiplex PCR protocol
Name of the primer
Sequence (5’ to 3’)
Optimization of the multiplex PCR
Trial PCR amplifications were performed using different primer combinations for differentiating species A from D, species B from C, species B from E, and species B, C, and E from each other. The PCR conditions were optimized with respect to parameters such as MgCl2 concentration, primer concentration, Taq DNA polymerase concentration, and annealing temperature. After achieving a set of reaction conditions from which products of the correct size were generated from the corresponding primer pair and DNA combinations, an effort was made to minimize the number of PCR assays needed to differentiate all five members of the complex.
Agarose gel electrophoresis
Agarose powder (Promega) was used for the preparation of gels. The 1.5 % agarose gels which stained with ethidium bromide (0.5 μg/ml) were prepared. Amplified PCR products were loaded with Promega lambda marker and gel electrophoresis was carried out at 150 mV for 40 min to get a good separation in the amplified products. Following the gel electrophoresis, the migrated DNA was visualized and documented using gel documentation apparatus and Image Lab 3.0 software protocol (GelDoc ™XR+).
Validation of results
The results obtained from the gel electrophoresis were compared with the marker and band sizes of 359 bp, 248 bp, 95 + 248 bp, 166 + 359 bp and 178 + 248 bp were used to compare the sibling species A, B, C, D and E respectively.
In order to confirm the results in the experiment, about 50 μl of amplified PCR products using species-specific primers and universal primers were sent to Macrogen, Korea (10 F, 254 beotkott-ro, Guemcheon-qu, Seoul, Republic of Korea).
Optimization of multiplex PCR assay for AD- PCR
The amplifications were conducted for species specific primers using 25.0 μl of solution containing 1.0 μl each of ADF (10 μM), ADR (10 μM) and DR (10 μM), 7.5 μl of Master mix [10× PCR buffer + MgCl2 (15 mM) + dNTPs- all four types (2 mM)], 1.0 μl of Taq polymerase (2 U/μl), 11.5 μl PCR water and 2.0 μl of 40× DNA template. The cycling conditions were hot start at 95 °C for 5 min followed by 40 cycles each of denaturation at 95 °C for 40 s, annealing for 48 °C for 1 min, and extension at 68 °C for 1 min, followed by the final extension at 68 °C for 10 min. The expected fragment sizes in the multiplex AD-PCR assay are follows: species A, 359 bp; species D, 166 + 359 bp.
Optimization of multiplex PCR assay for BCE-PCR
Trial PCR amplifications were conducted using 25.0 μl of solution containing 1.0 μl each of BCEF (10 μM), BCR (10 μM), ER (10 μM) and CR (10 μM), 7.5 μl of Master mix [10× PCR buffer + MgCl2 (15 mM) + dNTPs- all four types (2 mM)], 1.0 μl of Taq polymerase (2 U/μl), 10.5 μl PCR water and 2.0 μl of 40× DNA template. The cycling conditions were hot start at 95 °C for 5 min followed by 40 cycles each of denaturation at 95 °C for 40 s, annealing for 48 °C for 1 min, and extension at 68 °C for 1 min, followed by the final extension at 68 °C for 10 min. The expected fragment sizes in the multiplex BCE-PCR assay are follows: species B, 248 bp; species C, 95 + 248 bp; species E, 178 + 248 bp.
Validation of sibling species status of the field caught samples
Accurate identification of anopheline mosquitoes is necessary for planning effective vector control strategies and for a better understanding of their potential role in malaria transmission . Species identification based on morphological features is still practice by the field entomological teams. It must be emphasized still the use of morphological keys is an important and indispensable method for field identification of mosquitoes.
However, when comes to species complexes, they cannot be identified by morphological features, since they are morphologically similar. The importance of accurate identification of sibling species is, the ability to transmit the disease may vary from one sibling species to another. Therefore, it may adversely affect vector control efforts, if the vector is not correctly identified [17–20]. Hence, cytogenetical and molecular biological techniques are mandatory for a reliable identification of sibling species in a species complex.
A previous study conducted by Goswami et al. had used a multiplex BCE-PCR assay , while in the present study, the PCR amplifications were conducted for species-specific primers. Molecular assays targeting internal transcribed spacer 2 (its2) and co-II genes are most commonly used in the identification of sibling species. However, the its2 method has some limitations, since it can only differentiate sibling species into two sets as AD and BCE . The co-II region does not have such drawbacks since this region is a conserved area. Therefore, this method can be easily used to differentiate the sibling species status undoubtedly.
Investigations on larval breeding sites suggests that species E is able to exploit a wide range of breeding habitats with different limnological characteristics . There are some evidence in Sri Lanka to prove that the ability of anopheline mosquitoes breeding in waste water habitats with low dissolved oxygen level . However, it was not characterized the sibling species status of the An. culicifacies. This is the first time the molecular status of An. culicifacies s.l. breeding in waste water habitats was characterized. According to the results, the majority of the species belongs to sibling species E, which is considered as the vector sibling species of An. culicifacies. This study also evidences that the ability of An. culicifacies E to breed in a wide range of breeding habitats. In addition, Anopheles subpictus, Anopheles barbirostris, Anopheles peditaeniatus, Anopheles nigerrimus and Culex tritaeniorhynchus species were also detected along with the An. culicifacies s.l. in these habitats.
There may be similar breeding habitats from other parts of the country, which allow breeding of anopheline mosquitoes. However, malaria control programmes in most part of the countries focus on rural communities, as a result of larval bio-ecology of anophelines. Therefore, unusual breeding habitats, such as waste water collections, may mislead the current vector control programmes. This investigation adds new information to the existing knowledge on malaria molecular entomology, as well as supporting more effective vector control programmes through accurate confirmation of vector species.
The present study reveals, for the first time in Sri Lanka, the ability of An. culicifacies sibling species B and E breeding in drains containing waste water. It is suggested that health authorities should adopt new vector control interventions to eliminate these breeding habitats. This requires entomological surveillance in urban areas to detect potential vector breeding habitats.
This research work was supported by the Global Fund for Aids, Tuberculosis and Malaria (GFATM-Round 8) through TEDHA malaria elimination programme.
- Harrison BA. The Myzomyia series of Anopheles (Cellia) in Thailand, with emphasis on intra-specific variations (Diptera: Culicidae). Contr Am Entomol Inst (Ann Arbor). 1980;17:1–195.Google Scholar
- Van Bortel W, Sochanta T, Harbach RE, Socheat D, Roelants P, Backeljau T, et al. Presence of Anopheles culicifacies B in Cambodia established by the PCR-RFLP assay developed for the identification of Anopheles minimus species A and C and four related species. Med Vet Entomol. 2002;16:329–34.View ArticlePubMedGoogle Scholar
- Manonmani AM, Sadanandane C, Sahu SS, Mathivanan A, Jambulingam P. rDNA-ITS2-PCR assay for grouping the cryptic species of Anopheles culicifacies complex (Diptera: Culicidae). Acta Trop. 2007;104:72–7.View ArticlePubMedGoogle Scholar
- de Silva BGDNK. Anopheles culicifacies complex: Geographical distribution of sibling species and existing methods for their identification. Vidyodaya J Sci. 2009;14:1–28.Google Scholar
- Kar I, Subbarao SK, Eapen A, Ravindran J, Satyanarayana TS, Raghavendra K, et al. Evidence for a new malaria vector species, species E, within the Anopheles culicifacies complex (Diptera: Culicidae). J Med Entomol. 1999;36:595–600.View ArticlePubMedGoogle Scholar
- Subbarao SK. Anopheline species complexes in South-East Asia. New Delhi: India, World Health Organization Regional Office for South-East Asia; 1998.Google Scholar
- Subbarao SK, Sharma VP. Anopheline species complexes and malaria control. Indian J Med Res. 1997;106:164–73.PubMedGoogle Scholar
- Surendran SN, Abhayawardane TA, De Silva BG, Ramasamy R, Ramasamy MS. Anopheles culicifacies Y-chromosomal dimorphism indicates sibling species (B and E) with different malaria vectors potential in Sri Lanka. Med Vet Entomol. 2000;14:437–40.View ArticlePubMedGoogle Scholar
- Gunathilaka N, Fernando T, Hapugoda M, Wickremasinghe R, Wijeyerathne P, Abeyewickreme W. Anopheles culicifacies breeding in polluted water bodies in Trincomalee district of Sri Lanka. Malar J. 2013;12:285.View ArticlePubMed CentralPubMedGoogle Scholar
- Abhayawardana TA: Breeding preference of Anopheles culicifacies in Galketiyagama area in Kurunegala District. Proceedings of the Annual Session of Sri Lanka Association of Advancement Sciences, Colombo, Sri Lanka. 1995: 53–54Google Scholar
- Amarasinghe FP, Munasingha NB. A pre-development mosquito survey of the Mahaweli development project area, Sri Lanka: immatures. J Med Entomol. 1988;25:286–94.View ArticleGoogle Scholar
- WHO. Entomological field techniques for malaria control. Part 1 learner’s guide. Geneva: World Health Organization; 1992. p. 1–77.Google Scholar
- Amerasinghe FP. A guide to the identification of anopheline mosquitoes (Diptera; Culicidae) of Sri Lanka.1 adult females. Ceylon J Sci (Bio Sci). 1990;2:13–29.Google Scholar
- Gunathilaka N. Distribution of major and potential malaria vectors in Mannr and Trincomalee districts and systematics of anophelines in Sri Lanka. Sri Lanka: Ph.D thesis, University of Kelaniya; 2014.Google Scholar
- Coen E, Strachan T, Dover G. Dynamics of concerted evolution of ribosomal DNA and histone gene families in the melanogaster species subgroup of Drosophila. J Mol Biol. 1982;158:17–35.View ArticlePubMedGoogle Scholar
- Goswami G, Singh OP, Nanda N, Raghavendra K, Gakhar SK, Subbarao SK. Identification of all members of the Anopheles culicifacies complex using allele-specific polymerase chain reaction assays. Am J Trop Med Hyg. 2006;75:454–60.PubMedGoogle Scholar
- Gunathilaka N, Fernando T, Hapugoda MD, Wijeyerathne P, Wickramasinghe AR, Abeyewickreme W. Increasing trend of Anopheles culicifacies (Giles) and other potential malaria vector mosquitoes in Trincomalee District, Sri Lanka. Proceedings of the Annual Session of Sri Lanka Association of Advancement Sciences, Colombo, Sri Lanka. 2013;69:17.Google Scholar
- Gunathilaka N, Ranathunga B, Hapugoda MD, Abeyewickreme W, Wickremasinghe AR. Morphological identification and dynamics of malaria vectors in Trincomalee District of Sri Lanka. Proceedings of the Annual Session of Sri Lanka Association of Advancement Sciences, Colombo, Sri Lanka. 2014;70:8.Google Scholar
- Premaratne R, Ortega L, Janakan N, Mendis KN. Malaria elimination in Sri Lanka: what it would take to reach the goal. South- East Asian J Publ Health. 2014;3:85–9.Google Scholar
- Gunathilaka N, Hapugoda M, Abeyewickreme W, Wickremasinghe R: Entomological investigations on malaria vectors in some war-torn areas in the Trincomalee District of Sri Lanka after settlement of 30-year of civil disturbance. Malar Res Treat 2015, 2015:367635 http://0-dx.doi.org.brum.beds.ac.uk/10.1155/2015/367635.
- Goswami G, Raghavendra K, Nanda N, Gakhar SK, Subbarao SK. PCR-RFLP of mitocondrial cytochrome oxidase subunit II and ITS2 of ribosomal DNA: markers for the identification of members of the Anopheles culicifacies complex (Diptera: Culicidae). Acta Trop. 2005;95:92–9.View ArticlePubMedGoogle Scholar
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.