|Year : 2020 | Volume
| Issue : 2 | Page : 153-160
Genetic characterization of Aedes aegypti (Diptera: Culicidae) in Sri Lanka based on COI gene
H.A.K.M. Dharmarathne1, O.V.D.S.J. Weerasena1, K.L.N.S Perera2, GH Galhena3
1 Institute of Biochemistry, Molecular Biology and Biotechnology, University of Colombo, Colombo 03, Sri Lanka
2 Genetech, Colombo 08, Sri Lanka
3 Department of Zoology and Environment Sciences, Faculty of Science, University of Colombo. Colombo 03, Sri Lanka
|Date of Submission||01-Nov-2018|
|Date of Acceptance||23-Jul-2019|
|Date of Web Publication||14-Jul-2021|
Dr G H Galhena
Department of Zoology and Environment Sciences, Faculty of Science, University of Colombo, Colombo 00300
Source of Support: None, Conflict of Interest: None
Background & objectives: Aedes aegypti is the most prominent vector for dengue virus worldwide. Accurate identification of the species and understanding its colonization pattern are essential prerequisites in vector control. Thus, the present study was aimed to genetically characterize Ae. aegypti mosquitoes collected from different regions of Sri Lanka based on mitochondrial COI gene.
Methods: Thirty-three Ae. aegypti larval samples were collected from 19 districts. A 735bp region of the mitochondrial COI gene was amplified and analyzed for genetic diversity indices. Phylogenetic trees were constructed using Sri Lankan samples and also including mosquito samples reported from other parts of the world.
Results: High genetic diversity was observed within the samples analysed (gene diversity: 0.949; average number of nucleotide differences: 6.371). There were 20 haplotypes presented within the 19 localities investigated. The phylogenetic tree derived two main clades. However, no distinguishable clustering pattern was observed in the phylogenetic tree except for the districts in the northern corner indicating extensive admixing among different populations. When samples from other countries were included in the phylogenetic tree, Anuradhapura, and Mannar samples were clustered together with samples from India, Venezuela, USA, Portugal and Cambodia while Rathnapura was clustered with Bolivia and France.
Interpretation & conclusion: Our results suggest that Sri Lanka has undergone multiple invasions of Ae. aegypti from various parts of the world over an extensive period. Further, the mosquito control campaigns had not caused a significant effect on the Ae. aegypti populations which is existing in mutation-drift equilibrium.
Keywords: Aedes aegypti; COI gene; genetic characterization; genetic diversity; vector control; population admixture
|How to cite this article:|
Dharmarathne H, Weerasena O, Perera K, Galhena G H. Genetic characterization of Aedes aegypti (Diptera: Culicidae) in Sri Lanka based on COI gene. J Vector Borne Dis 2020;57:153-60
|How to cite this URL:|
Dharmarathne H, Weerasena O, Perera K, Galhena G H. Genetic characterization of Aedes aegypti (Diptera: Culicidae) in Sri Lanka based on COI gene. J Vector Borne Dis [serial online] 2020 [cited 2021 Aug 5];57:153-60. Available from: https://www.jvbd.org/text.asp?2020/57/2/153/310871
| Introduction|| |
The threat delineates by mosquitoes has grown rapidly all over the world during the past few decades. Enhanced tolerance towards environmental stress along with progressive gain of insecticide-resistance is partly responsible for their enhanced vectorial capacities. The adaptive changes exhibited by different isolates/populations of the same mosquito species found in different parts of the world are found to be strikingly different. As a result, successful vector mitigation strategies are likely to rely on genetic characterization of vector populations, which would provide valuable information on their genetic makeup, vectorial capacity, temporal stability and the processes relevant to disease transmission, including movement of the vectors across various regions.
Due to its many desirable characteristics as a molecular marker,, at present, frequently used barcode region for animals is a 5’-segment of the mitochondrial gene Cytochrome Oxidase I (COI) which is referred to as the ‘Universal’ or ‘Folmer’ region. This region is the standard marker chosen by the Barcode of Life Data System (www. boldsystems.org) (BOLD), an online platform for assembling and curating DNA barcoding data from all over the world. While most of mosquito barcoding studies use this region, some studies have used different areas of COI.
Several previous studies have successfully attempted the use of COI gene in genetic characterization of Aedes aegypti, a vector for many arboviral diseases including dengue and chikungunya, at a regional scale. For example, Scarpassa et al in 2008, showed the existence of two mitochondrial lineages separated by eight fixed mutations, distributed among 14 localities in Brazil, based on 852-bp region of COI. The use of both COI and NADH4 together, replicated the same results for Brazilian Aedes aegypti populations indicating the relative robustness of COI region in phylogenetic analysis. Another study conducted with COI region, covering the north, southeast, northeast and central regions of South India, has reported the presence of three genetic lineages within the 31 localities studied, but with limited genetic diversity.
An effort to use COI region to investigate Sri Lankan Aedes aegypti mosquitoes was conducted recently. However, this study has covered only five administrative districts and thus provides only a limited understanding of the phylogenetic relationships of mosquito samples countrywide. Hence, the current study was designed for molecular characterization of the Aedes aegypti mosquitoes collected from all parts of Sri Lanka based on mito-chondrial COI to provide a wider understanding of the phylogenetic relationships of Ae. aegypti mosquitoes that exist throughout Sri Lanka. The genetic relatedness between Sri Lankan mosquitoes and those reported from other parts of the world was also investigated to understand historical migratory routes.
| Material & Methods|| |
The fourth stage larval samples of Aedes aegypti were collected from different regions of Sri Lanka using ovitraps or standard dipping method. Initially, a total of 50 sampling sites were selected from all the 25 administrative districts (2 sites per district) for sample collection with a considerable distance (at least 10 km) between the two sites. However, only 33 sites representing 19 districts were positive for Ae. aegypti larvae and were thus included in the study [Figure 1] and [Table 1]. All samples were examined under a light microscope and morphologically identified to species level using taxonomic keys. Identified larvae were preserved in 95% ethanol. Only one larval sample from each site was used for the study.
|Figure 1: District map of Sri Lanka showing sample collection locations (1-33 denote the collection sites).|
Click here to view
PCR amplification and sequencing
Genomic DNA was extracted from the samples using Wizard® Genomic DNA Purification Kit. A 735bp region of the mitochondrial COI gene was amplified using oli-gonucleotide primers as described by Gupta et al (2016). PCR products were run on 2% agarose and the excised amplicon were purified using Wizard® SV gel and PCR clean up system (Promega, USA). Bidirectional sequencing of amplified fragments was conducted with the Big Dye® Terminator v3.1 cycle sequencing kit and 3500Dx Genetic analyzer (Thermo Fisher Scientific, USA).
The consensus sequence was generated for each sample using BioEdit version 6.0.7 and were searched over the GenBank database using Basic Local Alignment Search Tool (BLAST) against the Ae. aegypti genomes in GenBank. The complete sequences were deposited in GenBank [Table 1].
Local samples were compared with the Ae. aegypti found from other countries [Table 2] using Multiple Sequence Alignment (MSA) based on the sequences available in GenBank. Phylogenetic trees were built using maximum likelihood method (with 1000 bootstraps) with Kimura a cluster containing >50% bootstrap support, was considered significant. All sequences were analyzed with DnaSP version 5.10.01 (Librado and Rozas, 2009) for nucleotide diversity indices (the number of monomorphic sites, polymorphic sites, singleton sites and parsimony informative sites) and for the indices of haplotype diversity (the number of haplotypes, haplotype diversity and the average number of nucleotide differences).
|Table 2: GenBank Accession numbers of previously published Aedes aegypti COI sequences used in the current analysis|
Click here to view
This study was approved by the Ethics Review Committee of Institute of Biology, Sri Lanka (ERC IOBSL 122 04 15)
| Results|| |
The sequence analysis was carried out using a 683bp fragment of the amplified region of COI gene using a total of 33 samples derived from the present study representing 19 districts. All sequences were deposited in GenBank repository (https://www.ncbi.nlm.nih.gov/genbank/) and the accession numbers are given in [Table 1].
The sequences consisted of 661 monomorphic sites (96.78%), 22 polymorphic sites (3.22%), 4 singleton sites (0.59%) and 18 parsimony informative sites (2.63%). This diversity was observed to generate 20 haplotypes as shown in [Table 3]. The haplotype (gene) diversity and average number of nucleotide differences were 0.949 and 6.371 respectively. The Tajima’s D value for Sri Lankan isolates was 0.60607; p >0.05.
|Table 3: Haplotype diversity, frequency and distribution in different districts in Sri Lanka.|
Click here to view
Phylogenetic analysis of Aedes aegypti
The analysis of evolutionary divergence within Sri Lankan samples of Ae. aegypti, resulted in two main clades with the second clade bearing two subclades as shown in [Figure 2]. The first clade includes more samples from the southern part of the country (i.e. Hambantota, Colombo, Gampaha, Monaragala and Rathnapura) while second clade carries more from the northern region (i.e. Mannar, Vavuniya, Ampara, Jaffna, Anuradhapura, Puttalam, Batticaloa). Interestingly, all districts in the northern corner of the country i.e. Jaffna, Vavuniya, Ampara and Mannar have clustered separately. Even though samples from Gampaha, Colombo, Hambantota, Kandy and Batticaloa districts were distributed among both main clades, samples from other districts were found within only one main clade.
|Figure 2: Phylogenetic tree showing the relationship of COI sequence of Sri Lankan Ae. aegypti samples collected from different districts.|
Click here to view
An effort to understand the relationship between Sri Lankan Ae. aegypti samples and samples reported from other countries, were also made by constructing a phylogenetic tree that consisted of 16 other COI sequences from GenBank database [Figure 3]. Since India is the closest country to Sri Lanka, with frequent cargo and passenger transportation throughout history, it is likely that passive transportation of mosquitoes between the two countries has resulted in a higher genetic similarity among the samples derived from the two countries. To analyse this possibility, five COI sequences representing five different areas of India is also included in the phylogenetic tree. As shown, Sri Lankan samples of the present study did not cluster separately when analyzed together with the samples collected from different regions of the world, but showed varying degrees of relationship with them.
|Figure 3: Phylogenetic tree showing the relationship among COI sequence of Sri Lankan samples and samples of Ae. aegypti derived from various parts of the world.|
Click here to view
In the phylogenetic tree derived from this analysis, two main clades were observed where Madagascar sample was separated from all other countries in a separate clade. Anuradhapura, and Mannar samples were clustered together with isolates from India, Venezuela, USA, Portugal and Cambodia. Sample from Rathnapura district that was clustered with Bolivia and France was found in the second clade.
| Discussion|| |
The overall objective of this study was to genetically characterize the Aedes aegypti mosquitoes collected from different regions of Sri Lanka in relation to Ae. aegypti mosquitoes found from various parts of the world based on 683bp fragment of mitochondrial COI gene. Sequences of 33 mosquito samples from 19 administrative regions were analysed for genetic diversity and phylogenetic relationships. Our results suggest that Sri Lanka has undergone multiple invasions of Ae. aegypti from various parts of the world over on extensive period. The mosquito population is existing in mutation-drift equilibrium unrestrained by breeding habitat removal or chemical control.
According to our results 20 haplotypes were present within the 19 localities investigated. Although the number of samples analysed does not allow the detection of the pattern of haplotype distribution, the presence of 20 haplotypes among the 33 mosquito isolates provides evidence for extensive genetic diversity. Additionally, this observation supports the suggestion made previously on the relative high variability of Ae. aegypti COI gene. According to Mousson et al (2005), Ae. aegypti COI gene is more variable than the ND5 (nicotinamide adenine dinucleotide dehydrogenase subunit 5) or Cytb (cytochrome b) genes in Brazilian samples. The high haplotype diversity values (ranging from 0.2 to 0.8) observed among the Brazilian samples supported this suggestion. The value we have observed in the present study is even greater suggesting an extensive genetic diversity of the mosquitoes. This diversity might be indicative of several facts such as a large effective population size of mosquitoes within the island, high gene flow among local populations and favorable environmental factors for mosquito breeding, irrespective of the ongoing mosquito control campaigns. The observed nucleotide diversity (Pi=0.0093) and the average number of nucleotide differences (K=6.371) was also relatively high and higher than what had been reported from Brazil for the same gene region (Pi=0.0065; K=5.546). This implies the presence of relatively large differences between haplotypes suggesting a long-established mosquito population within the country that has not suffered recent bottlenecks. However, it is equally possible that the presence of genetically distant haplotypes in sympatry to produce the same observations.
The Tajima D value calculated is also in favour of this hypothesis. Tajima’s D test distinguishes between DNA sequences evolving randomly (neutrally) from one’s that are evolving under a non-random process. Statistically non-significant Tajima’s D value obtained for the present study accepts the null hypothesis of neutral evolution. This indicates that the mosquito populations within the country is at mutation-drift equilibrium maintaining a stable large population.
When the phylogenetic tree drawn for the Sri Lankan samples from the 19 districts were considered, there was no distinguishable clustering observed except for the three districts of the northern corner [Mannar, Vavunia (one sample each), Jaffna and Ampara]. They clustered together when analysed along with 20 samples around the globe including the five samples from different regions in India. All other isolates which represented various districts in Sri Lanka were randomly distributed over the phylogenetic tree. However, one would expect mosquito samples from geographically closer regions to be more genetically similar than those collected from far apart when considering their dispersal distances. Honório et al (2003), reported Ae. aegypti to have a general dispersal distance around 800m. This suggests of restricted genetic mixing of the mosquitoes collected from several hundred kilometers away from each other. Contradictory results observed in the present study to this assumption suggest mosquitoes were able to surpass their biological dispersal barrier perhaps through the means of passive transportation. It is possible that passenger and cargo transportation through various motor vehicles among different districts would have served as an ideal passage for both the adult, larvae and mosquitoes eggs from one district to another allowing extensive admixing. The limited admixture observed between the mosquitoes of northern most districts and the rest of the Sri Lanka is likely to have been caused by the civil war that lasted for decades until it ended in 2009. The main routes of transportation were re-opened for civil usage since then, and this might have caused some amount of mosquito migrations between northern area and others. The grouping of two northern samples (MNA1 and VVA2) with samples from western and central provinces of the country could have resulted from these recent advances in transportation. This observation agrees with the previous reports of low levels of genetic differentiation among Aedes aegypti samples from different regions of Sri Lanka based on microsatellite markers.
Despite the extensive admixing observed among Sri Lankan samples, two main clades were identified in the phylogenetic tree. This pattern of clustering was consistent when the samples were analysed with samples derived from various parts of the world. Interestingly, Venuzuela, Andhra Pradesh in India, USA, Portugal, Cameroon and Cambodia samples clustered together with the samples derived from the northern parts of Sri Lanka. Vietnam, Tanzania, Mexico, Brazil, Thailand and Thirumala, Andhra Pradesh in India made a separate sub-clade which did not carry a single Sri Lankan sample indicating a subtle but distinct genetic difference between the Sri Lankan samples and samples from these countries. The second main clade of Sri Lankan phylogenetic tree clustered with samples of Bolivia and France shows the close genetic similarity of some of the Sri Lankan samples to these countries. Unfortunately, the phylogenetic tree with global samples does not have sufficient bootstrap values to make strong inferences, probably due to the low sample numbers used in the analysis. Nevertheless, it is highly likely that Sri Lanka received multiple invasions from Ae. aegypti mosquitoes from various parts of the world as indicated in the phylogenetic tree (Vietnam, Tanzania, Mexico, Brazil, Thailand, Thirumala, Andhra Pradesh in India, Bolivia and France). This agrees with the observation made previously by other authors upon the recurrent spreading of mosquitoes over various continents producing genetically similar mosquito samples among geographically faraway regions like South America, Asia and Africa.
The extensive genetic mixing of Ae. aegypti population in Sri Lanka and around the world as evident from our results has important implications on mosquito control and disease transmission. Firstly, high genetic diversity and extensive genetic mixing increases the evolutionary fitness of the mosquitoes and stabilizes their populations. This would be important for mosquitoes for surviving in an environment where numerous environmental stresses are constantly available due to the continuous mosquito control measures. Secondly, such genetic mixing would favour the spread of evolutionary important traits in mosquitoes such as insecticide resistance and salinity tolerance over a wider geographic area enhancing their adaptive fitness. As an invasive species, this would make their status even stronger. Thirdly, as a highly successful vector of many arboviral diseases, extensive dispersal of mosquitoes indicates spreading of diseases over larger geographic area resulting in emergence of mosquito borne diseases in new localities. Thus, the passive long distance migrations of Aedes aegypti indicated in our results highlights the necessity to prevent human-aided mosquito dispersals to ensure successful dengue control.
However, it is important to understand the relativity of these inferences in the context of genetic stability of different marker regions in the genome. In a study where both ND4 and COI genes of Aedes aegypti have been used to measure the nucleotide diversity, COI gene showed significantly lower diversity compared to the ND4 gene. It was postulated that this difference may have arose due to the higher constraints (i.e. balancing selection) acting upon the mutation rate of the COI gene. This indicates that COI gene in Ae. aegypti is more conservative than the ND4 gene. Thus, it is possible that a higher genetic diversity with a somewhat different phylogenetic picture is to be yielded, if ND4 gene has been used in the analysis along with the COI gene.
Despite such limitations, knowledge on vector genetics and the gene flow between mosquito populations residing at various geographically distant localities is beneficial in designing effective vector control strategies. This is the first study that genetically characterized Aedes aegypti samples from all the 19 districts of Sri Lanka, where the species is routinely reported. Though this paper only presents the results of a preliminary investigation conducted using a limited number of samples and with a single genomic region, the results revealed are of significant epidemiological value. As a country which is extensively affected by dengue fever related morbidity and mortality, the need to look for means of reducing passive mosquito dispersal within the country would be of utmost importance to keep the disease transmission under control. The mosquito population with high evolutionary potential points at the necessity to revisit the vector control strategies to yield success. While more elaborate studies are required to delve into the role of control programs in manipulating the genetic composition of the mosquito populations, it is recommended to incorporate the information on other genetic markers such as mitochondrial NADH dehydrogenase subunit 4 (ND4) in such analysis to improve the robustness of the reported observations.
Conflict of interest: None
| Acknowledgements|| |
This work was supported jointly by the Collaborative Research Grant (AP/3/2/204/RG/05) and the Institute of Biochemistry, Molecular Biology and Biotechnology of the University of Colombo, Sri Lanka.
| References|| |
Brown JE, Obas V, Morley V, Powell JR. Phylogeography and spatio-temporal genetic variation of Aedes aegypti
(diptera: culicidae) populations in the Florida keys. J Med Entomol
Lunt DH, Zhang DX, Szymura JM, Hewltt OM. The insect cytochrome oxidase I gene: Evolutionary patterns and conserved primers for phylogenetic studies. Insect Mol Biol
Scarpassa VM, Cardoza TB, Junior RPC. Population genetics and phylogeography of Aedes aegypti
(Diptera: Culicidae) from Brazil. Am J Trop Med Hyg
Batovska J, Blacket MJ, Brown K, Lynch SE. Molecular identification of mosquitoes (Diptera: Culicidae) in southeastern Australia. Ecol Evol
Gloria-Soria A, Ayala D, Bheecarry A, Calderon-Arguedas O, Chadee D, Chiappero M et al
. Global genetic diversity of Aedes aegypti. Mol Ecol
Paduan KDS, Ribolla PEM. Mitochondrial DNA Polymorphism and heteroplasmy in populations of Aedes aegypti
in Brazil. J Med Entomol
Vadivalagan C, Karthika P, Murugan K, Panneerselvam C, Paulpandi M, Madhiyazhagan P, et al
. Genetic deviation in geographically close populations of the dengue vector Aedes aegypti
(Diptera: Culicidae): Influence of environmental barriers in South India. Parasitol Res
Fernando HSD, Hapugoda MD, De Silva BGDNK. A preliminary study to analyze Aedes aegypti
diversity in Sri Lanka using mitochondrial DNA variations. Proceedings of International Conference on Multidisciplinary Approaches 2015; University of Sri Jayewardenepura, Nugegoda.
Rattanarithikul R, Panthusiri P. Illustrated keys to the medically important mosquitoes of Thailand. Southeast Asian J Trop Med Public Health
1994; 25(suppl 1): 66–9.
Gupta K, Dhawan R, Kajla M, Kumar S, Jnanasiddhy B, Naveen K, et al
. Molecular identification of Aedes aegypti
mosquitoes from Pilani region of Rajasthan, India. J Vector Borne Dis
Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics
Mousson L, Dauga C, Garrigues T, Schaffner F, Vazeille M, Failloux A. Phylogeography of Aedes
(L.) and Aedes
(Skuse) (Diptera: Culicidae) based on mitochondrial DNA variations. Genet Res
Honório N, Silva W, Leite P, Gonçalves J, Lounibos L, Lourenço-de-Oliveira R. Dispersal of Aedes aegypti
and Aedes albopictus
(Diptera: Culicidae) in an urban endemic dengue area in the State of Rio de Janeiro, Brazil. Mem Inst Oswaldo Cruz
Nirmani MD, Gangani PD, Perera KLN, Galhena GH. Long distance passive migrations may explain the low genetic differentiation in principal dengue vector (Aedes aegypti)
in Sri Lanka. Proceedings of the 36th Annual Sessions of the Institute of Biology 2016.
Failloux AB, Vazeille M, Rodhain, F. Geographic genetic variation in populations of the dengue virus vector Aedes aegypti. J Mol Evol
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2], [Table 3]