Insecticide susceptibility of the sand fly leishmaniasis vector Phlebotomus argentipes in Sri Lanka
Parasites & Vectors volume 13, Article number: 246 (2020)
Leishmania donovani-induced and sand fly-transmitted leishmaniasis is a growing health problem in Sri Lanka. Limited knowledge on biological and behavioral characteristics of probable vector Phlebotomus argentipes hinders disease control. Here, insecticide susceptibility patterns of P. argentipes were investigated with exploration of probable underlying resistance mechanisms.
Adult sand flies were collected using standard cattle baited net traps and CDC light traps from selected sites in four districts. Adult F1 progeny of P. argentipes were exposed to different concentrations of DDT, malathion, deltamethrin and propoxur using WHO susceptibility bioassay kits. Post-1-h knockdown and post-24-h mortality were recorded and analyzed. Metabolic enzyme activity and the sensitivity of the acetylcholinesterase target-site were determined by biochemical assays using wild-caught flies. Extracted fly DNA samples were tested for the presence of knockdown-resistance (kdr) type mutations.
The LC100 values for DDT, malathion, propoxur and deltamethrin were 0.8–1.5%, 0.9–2.0%, 0.017–0.03% and 0.007% respectively. Insecticide-susceptibility levels were higher than the discriminating dosages established for Aedes mosquitoes, except for malathion. The lowest susceptibility levels (except for deltamethrin) were detected in the Mamadala population, whereas the highest levels were detected in the Mirigama population. The percentage of knocked-down sand flies was < 75% at any tested concentration, including those, which exhibited 100% mortality after 24 h. Elevated activity levels of glutathione S-transferase (3%, 7%, 12.5% and 14%) and esterase (2%, 5%, 5.5% and 6.5%) were detected in flies that originated from Mirigama, Pannala, Thalawa and Mamadala respectively, while monooxygenase quantities remained below the cut-off level. Ten to 34.5% of flies were heterozygous for acetylcholinesterases target-site insensitivity, associated with organophosphate and carbamate resistance. Pyrethroid-resistance-associated L1014F kdr-type mutation in the voltage gated sodium channel gene was detected in 30/53 flies.
Populations of P. argentipes in Sri Lanka are largely susceptible to common insecticides, except for malathion (used extensively in the past for malaria control). Their insecticide susceptibility appears negatively associated with past malaria endemicity of the study sites, with signs of early insecticide tolerance. Presence of insecticide target site insensitivity in a notable proportion of flies and enhanced insecticide metabolizing enzyme activities imply potential future challenges for leishmaniasis control, with a call for urgent proactive measures for its containment.
Phlebotominae sand flies are the vectors of leishmaniasis, a disease which extends to over 98 countries across the world . Cutaneous leishmaniasis (CL), muco-cutaneous leishmaniasis (MCL) and visceral leishmaniasis (VL) or “kala-azar” are the three main clinical forms of the disease. The causative agent is an obligate intracellular flagellated protozoan belonging to the genus Leishmania. Phlebotomus argentipes (Diptera: Psycodidae) is the known vector of Leishmania donovani, the causative agent of leishmaniasis in Sri Lanka [2, 3]. The same causative species of parasite exists in India, Nepal and Bangladesh, with VL as the recognized disease . Leishmaniasis is a rapidly growing health threat in Sri Lanka  with no national control programme yet in place to contain the situation. Among the three types of the disease, VL is the most virulent form and accounts for the second largest number of deaths due to parasitic causes . The majority of L. donovani infections in Sri Lanka manifest as CL with only a few MCL and VL cases . Although it was an exclusively imported disease prior to 1990s , it is now a widely prevalent disease, with case numbers increasing since 2001 and has been considered ‘notifiable’ in the health sector since 2008 [6, 8].
Control of leishmaniasis is achieved through the interruption of the cycle of transmission. Widely used methods include early infection detection and treatment of patients and control of the vector and its reservoir hosts. Vector control, although often used as a strategy for leishmaniasis control, is restricted due to the difficulties of locating the larval habitats of sand flies . Therefore, vector control mostly relies on the control of adult flies often through the use of chemical insecticides. Four groups of chemical insecticides, i.e. organochlorines, organophosphates, carbamates and pyrethroids are generally used for vector control . However, insecticide use is advocated only after careful studies on biological properties of the vector, including its resting behavior and insecticide susceptibility patterns. Haphazard or over-use of such chemicals often leads to development of tolerance against those agents.
Insecticides continue to be used in Sri Lanka for mosquito control, with extensive use in the past for anti-malaria activities (pre-elimination phase), which was in the form of regular indoor residual spraying programmes that continued until 2012 . Therefore, sand flies may have been regularly exposed to insecticides, at least in malaria endemic areas in the dry zone, which covers two-thirds of the country . Insecticide resistance in sand flies has been reported in several geographical regions around the world, particularly in malaria endemic areas, including the neighboring country India . In India, DDT was replaced with pyrethroids for indoor residual spraying due to its apparent ineffectiveness in the control of VL, especially in Bihar [14, 15]. Resistance to deltamethrin has also been reported in Pondicherry, India . DDT usage for indoor residual spraying programmes for the control of visceral leishmaniasis vectors in Bangladesh was banned in 1998  and in Nepal 1995  due to its hazardous effects. Thereafter, synthetic pyrethroids were successfully introduced, in Nepal and Bangladesh in 1992 and 2012, respectively . It is believed to have contributed to the achievement of VL-elimination in Bangladesh by 2017, reaching the World Health Organization’s (WHO) set target of less than one case per 10,000 population  with the sand flies largely remaining susceptible to the insecticides in use .
Insects can resist toxic insecticidal compounds through several physiological mechanisms, including interference with enzyme systems via production of detoxifying enzymes and/or effects on target site(s) of the insecticide, rendering them insensitive. The neurotransmitter acetylcholine is broken down by acetylcholinesterase (AChE), which is the target site for the insecticides organophosphates and carbamates. Therefore, mutations in the AChE gene can bring about resistance to these insecticides [10, 22].
Sodium channel regulatory protein is the target site of pyrethroids and organochlorines (DDT and its analogues) and resistance to these groups of insecticides can develop due to the mutations in voltage gated sodium channel (VGSC) genes leading to knockdown resistance (kdr) . If this resistance occurs at higher levels, it is known as super-kdr [24, 25]. To date, more than 37 resistance-associated ‛kdr’-type mutations or combination of mutations have been detected in pyrethroid and DDT-resistant insect populations [26, 27]. Recently two kdr mutations at the codon 1014 including both L1014F and L1014S have been discovered in phlebotominae sand flies in India, which are located in the same codon regions as described in mosquitoes .
The present study describes for the first time, to our knowledge, the insecticide susceptibility status of P. argentipes in Sri Lanka with investigations into the probable underlying metabolic and genetic mechanisms that may confer insecticide resistance.
Collection and rearing of sand flies
Sand fly adults of both sexes were collected in four administrative districts of Sri Lanka: Thalawa (8°14′11.468″N, 80°21′2.782″E) in Anuradhapura district (North-Central Province), Pannala (7°19′43.608″N, 80°1′26.3316″E) in Kurunegala district (North-Western Province), Mamadala (6°07′16.80″N, 81°07′12.60″E) in Hambantota district (Southern Province), and Mirigama (7°13′30.72″N, 80°7′40.439″E) in Gampaha district (Western Province) (Fig. 1). All collections were made between November 2015 and November 2017. The country is arbitrarily divided into dry, intermediate and wet zones based on the annual rainfall. Thalawa and Mamadala are located within the dry zone, Pannala in the intermediate zone and Mirigama in the wet zone (Fig. 1). Study sites were selected based on leishmaniasis case burden , vector prevalence  and to represent different climatic zones . At each site, the standard cattle-baited net traps  and Center for Disease Control (CDC) light traps were used for sand fly collections. Traps were set at 18:00 h each evening and flies were collected each morning at 6:00 h.
Sand flies were transferred by aspirators from trap nets into plaster of Paris-lined plastic pots covered with a piece of net. A hole in the net (for insertion of the aspirator) was covered with a cotton swab soaked with 30% sucrose solution to feed sand flies. The samples were then transferred to the insectary of the Department of Parasitology, Faculty of Medicine, University of Colombo, Sri Lanka, for rearing. The methods of Volf & Volfova  were used thereafter, for isolation of newly emerged flies (F1 progeny) for insecticide bioassay tests. Identification was carried out to the species level by examination under a light microscope at a magnification of 400× using standard taxonomic keys [33, 34].
Insecticides, chemicals and equipment
Chemicals were purchased from Sigma-Aldrich (Poole, UK) unless otherwise stated. DDT (an organochlorine), malathion (an organophosphate), deltamethrin (a pyrethroid) and propoxur (a carbamate) (97–99% pure) were obtained from ChemService (Maidstone, UK). Biochemical assay plates were read using an EPOCH 2 absorbance microtitre plate reader with Gen 5.03 software (Molecular Devices, Bio-Tek, Vermont, USA) Protein assay kit (Bio-Rad, Hemel Hempstead, UK) was used for the detection of protein content in fly homogenates.
Insecticide-impregnated papers and insecticide bioassays
Insecticide-impregnated papers were prepared at the Department of Zoology, Faculty of Science, University of Peradeniya, Sri Lanka, as previously described , using the standard WHO method . In brief, the technical grade insecticides at a pre-determined range of concentrations were mixed with acetone (as the spreading agent). DDT (0.4–4.0%), malathion (0.5–5.0%) and propoxur (0.01–0.10%) solutions were dissolved in olive oil. Deltamethrin (0.005–0.050%) solutions were prepared in Dow-Corning 556 silicone fluid. Rectangles of Whatman No. 1 filter paper (12 × 15 cm) were used for insecticide impregnation. Insecticide/oil solutions (0.7 ml) were mixed with an equal volume of acetone and the mixture was applied uniformly on the filter paper. Acetone enabled an even distribution of the solution on the paper. Insecticide-impregnated papers were left overnight at room temperature to allow complete evaporation of acetone. Papers were then foil wrapped and stored at − 20 °C until use. Each paper was used for repeat experiments, but only up to 5 times .
Insecticide bioassays were performed using standard WHO procedures  to investigate the susceptibility pattern of P. argentipes for DDT, malathion, deltamethrin and propoxur insecticides. In each test, 20 female sand flies from the F1 population were placed in a holding tube lined with an insecticide-free paper for 0.5 h at 26 ± 1 °C to acclimatize. Then they were transferred to the testing tube with insecticide-impregnated filter paper. The experimental tubes were allowed to stand upright for 1 h; insects were then transferred to holding tubes for a further 24 h recovery period with 30% sugar at 26 ± 1 °C. Five replicates were used for each concentration of insecticide. The percentage of flies knocked-down after one hour of insecticide exposure and the percentage mortality after 24 h recovery period were calculated. Papers with oil and acetone were used as the negative controls for each experiment. Data were considered only when the mortality in the control tubes were less than 20%, and the results were validated with the control samples using Abbott’s formula: [(% test mortality − % control mortality) × 100]/[100 − % control mortality)] .
WHO-prescribed biochemical tests were performed, with slight modifications . Specific activity levels of esterase, glutathione S-transferase and the amounts of monooxygenases were determined in two-hundred wild-caught adult sand flies from each study site . Each fly was separately homogenized in 100 µl of ice-cold distilled water. Two replicates of 10 µl of a crude homogenate obtained after centrifugation at 13,000× rpm for 2 min were used for each of the above standard assays.
Another 200 flies from each study site were homogenized individually in 70 µl of ice-cold distilled water to assess the sensitivity of the acetylcholinesterase (AChE). Two replicates of 25 μl of crude homogenate were used .
Detection of kdr mutations
DNA was extracted using DNeasy® Blood and Tissue kit (Qiagen, Hilden, Germany) from individual flies with a total of 10–15 P. argentipes flies per site. Homogenates were incubated overnight after adding proteinase K at 56 °C until tissues were completely lysed.
PCR amplification of domain II of the VGSC gene was carried out using primers Vssc8F (5′-AAT GTG GGA TTG CAT GCT GG-3′) and Vssc1bR (5′-CGT ATC ATT GTC TGC AGT TGG T-3′)  in a 50 µl PCR reaction which contained 2 µl of each primer, 10 µl of 5× Colourless GoTaq® Flexi Reaction Buffer (Promega, Wisconsin, USA), 4 µl of 25 mM MgCl2, 1 µl of 10 mM of each dNTP and 5 units of GoTaq® Flexi DNA polymerase (Promega). PCR cycling conditions included an initial denaturation step of 5 min at 95 °C followed by 30 cycles each of 96 °C for 30 s, 56 °C for 30 s, and 72 °C for 30 s, and a final extension step at 72 °C for 5 min. PCR products were electrophoresed on a 1% agarose gel stained with ethidium bromide. Amplified PCR products were sent to Macrogen, Korea for purification and DNA sequencing using an Applied Biosystems 3730 DNA Analyzer. The sequences (GenBank: MN685050-MN685102) were analyzed using BioEdit software v.7 (http://www.mbio.ncsu.edu/BioEdit/bioedit.html) to determine the possible mutations of the VGSC gene that confer resistance at the pyrethroid target site sodium channel regulatory proteins (kdr-type resistance). Amplified domain II VGSC gene sequences were further analyzed and compared with the NCBI database VGSC domain II sequence of Musca domestica (GenBank: AAB47604).
The concentration versus mortality and knockdown relationships were determined for DDT, malathion, deltamethrin and propoxur using bar charts in Sigma Plot (version 10). Wilcoxonʼs signed-rank test was used to test whether the results are significantly different among the four populations for each insecticide. Boxplots were prepared using Minitab (version 15) for each insecticide to test significant differences using ANOVA of elevated enzymes and to determine the remaining enzyme activities of flies in each location. Allelic frequencies with L1014F mutation of flies in each location were calculated.
The concentration versus mortality and knockdown relationships obtained from the bioassay experiments are presented in Fig. 2. Although the sand flies were tested over a wide range of insecticide concentrations, the mortality percentages moved from 0% to 100% within a very narrow range (Fig. 2; Table 1).
Sand flies from all four field sites survived (LC0) the exposure to concentrations of 0.4% DDT, 0.5% malathion, 0.01% propoxur and 0.005% deltamethrin (Table 1). However, there was 100% mortality (LC100) of sand fly populations from all sites at 0.8% DDT, 0.9% malathion and 0.017% propoxur, except for flies from Mamadala. The Mamadala population reached 100% mortality at higher concentrations of insecticides, i.e. 1.5% DDT, 2.0% malathion and 0.03% propoxur. The lowest concentration which resulted in a mortality of 100% for deltamethrin was 0.007% across all the study sites.
DDT 0.8% and malathion 0.9% killed all flies from Pannala, Thalawa and Mirigama whereas, for the same concentrations of the insecticides, only 75% and 77% mortality, respectively, were observed for Mamadala flies. Similarly, propoxur 0.017% resulted in a mortality of 100% in Pannala, Thalawa and Mirigama populations while Mamadala had a mortality of only 75% for the same propoxur concentration. No mortality was observed for deltamethrin 0.005% in the fly populations except for the Mirigama population, where a 30% mortality was observed for 0.005% deltamethrin. However, 0.007% deltamethrin resulted in a mortality of 100% in all four populations.
Results of the statistical analysis showed that there was no significant difference in mortality between flies originating from the four populations for each of the insecticides (DDT, malathion, propoxur and deltamethrin) with no significant difference in mortality rates between specific concentrations of tested insecticides (P > 0.05; Fig. 2). Details of statistical results are provided in Additional file 1: Table S1.
Flies that were exposed to the tested insecticides for 1 h exhibited zero knockdown at concentrations below 0.8% (DDT), 0.9% (malathion), 0.017% (propoxur) and 0.005% (deltamethrin). The maximum knockdown rate did not exceed 75%, even at concentrations which resulted in a mortality of 100% after the 24 h recovery period (Fig. 2).
Enzyme activity profiles of glutathione S-transferase (GST) and carboxylesterase, and the quantity profiles of monooxygenases of sand flies from all four study sites are shown in Fig. 3a, b and c, respectively. Already established discriminating values for enhanced activities/amounts of these enzymes that contribute to resistance in anopheline mosquitoes are given as interrupted lines (except for cytochrome P450 monooxygenase, since its discriminating value as given for anophelines is 0.35 equivalent units, which is off the scale). Mean values of carboxylesterase and GST activities and the quantity of monooxygenase of the four sand fly populations are given in Table 2. According to the profiles and mean values, the majority of flies in each population did not show elevated activities/amounts of these detoxifying enzymes. Only 3%, 7%, 12.5% and 14% of Mirigama, Pannala, Thalawa and Mamadala populations, respectively, had GST activities higher than discriminating dosages as given for mosquitoes (0.4 µmol min−1 mg−1). Similarly, only 2%, 5%, 5.5% and 6.5%, respectively, had esterase activities higher than that value (0.25 µmol min−1 mg−1), while none of the flies had monooxygenase quantities higher than the mosquito cut-off value of 0.35 equivalent units of cytochrome P450. However, there was no significant difference in mean values of GST activity and the quantity of monooxygenases among four populations in this study (F(3, 788) = 2.02, P = 0.110 and F(3, 792) = 2.46, P = 0.062, respectively). Although the mean values of esterase activities were not significantly different in Mamadala, Thalawa and Pannala flies (F(2, 594) = 0.03, P = 0.972), the mean esterase activity of the Mirigama population was significantly higher than that of the other three populations (F(3, 792) = 5.00, P = 0.002) (Additional file 2: Text S1). Remaining AChE activity of the propoxur-inhibited crude homogenate fraction was calculated for each fly as a percentage of the activity of the uninhibited fraction. Remaining activity profiles and the average activities for all four populations are shown in Fig. 3d and Table 2. According to the WHO standards, remaining activity of < 30% indicates that the population is homozygous susceptible for insecticide inhibition whereas values of 30–70% indicate heterogeneity, and values > 70% indicate homozygous resistance status . After inhibition, less than 30% remaining activity of the AChE target site was shown by 90%, 72%, 68.5% and 65.5% of Mirigama, Pannala, Thalawa and Mamadala sand fly populations, respectively (Fig. 3d). The remaining activities of the AChE target site were comparable in the flies originating from Mamadala, Thalawa and Pannala (F(2, 594) = 0.12, P = 0.890) and significantly higher when compared with flies originating from Mirigama (F(3, 792) = 14.75, P < 0.0001) (Table 2, Additional file 2: Text S1).
Detection of kdr mutations
The L1014 mutation, a known mutation in insects that produces resistance to pyrethroids [26, 27], was discovered in 30 out of 53 sand flies tested. It is a non-synonymous mutation substituting leucine with phenylalanine, or rarely serine at position 1014. Out of 30 mutations detected in the present study, 29 had the L1014F mutation while only one sample had the L1014S mutation (Table 3, Additional file 3: Figure S1). The highest frequency of flies having the L1014F mutation was found in Mamadala (Table 3, Additional file 3: Figure S1).
The study describes, to our knowledge for the first time, early signs of insecticide tolerance as evident through bioassays and biochemical investigations, and the presence of kdr genetic mutation in phlebotomine sand flies in Sri Lanka.
Vector control is an established strategy used to arrest leishmaniasis transmission in endemic areas, especially through regular indoor residual spraying of chemical insecticides. Such strategies are being used in the neighboring countries India, Nepal and Bangladesh, with the aim of eliminating leishmaniasis by the year 2020 . The susceptibility of Sri Lankan sand flies to insecticides was studied using a standard methodology. The results revealed that the Sri Lankan sand fly populations were largely susceptible to the insecticides tested (except perhaps to malathion), when compared against the resistant discriminating dosages for mosquitoes. However, the range of concentrations of insecticides that gave values between LC0 and LC100 in sand flies was very narrow, and therefore, only a limited number of intermediate concentrations could be tested. Furthermore, regression lines for log probit mortality curves could not be obtained due to insufficient data points between 0 and 100%, which precluded the calculation of LC50 values.
Previous studies have also used discriminating concentrations or fixed dosages of insecticides to evaluate the susceptibility status of sand fly populations using mosquitoes as the standard [41,42,43]. For anopheline mosquitoes, the discriminating dosages have been recorded as 4%, 5%, 0.1% and 0.05% for DDT, malathion, propoxur and deltamethrin, respectively . For aedine mosquitoes, the values are 4% DDT, 0.8% malathion, 0.1% propoxur  and 0.03% deltamethrin . For the Sri Lankan sand fly populations tested in the present study, the lowest concentrations of insecticides that resulted in a mortality of 100% were 0.8–1.5 %, 0.9–2.0%, 0.017–0.03% and 0.007%, respectively, for the above insecticides. Therefore, the Sri Lankan populations of P. argentipes can be considered as susceptible to the tested insecticides, when compared against the discriminating dosages specified for either anopheline or aedine mosquitoes. The only exception was malathion for which the LC100 values were 1.5% for Mamadala, which was higher than the discriminating dosages given for aedine mosquitoes, i.e. 0.8% (while it was 0.8% in the other three sand fly populations tested). Similarly, sand fly populations from neighboring countries such as in Bihar State, India, and the border villages of Nepal have shown high malathion-resistance (and DDT-resistance) with ability to survive even at 5% malathion (and 4% DDT) . Moreover, Phlebotomus papatasi from Surogia village in Khartoum State, Sudan, and P. argentipes from Baizalpur, Chandi and Khusroopur villages in India, have also shown such high insecticide resistance .
Except for deltamethrin, the LC100 values were relatively higher in the Mamadala sand fly population when compared to the other three populations studied. Mamadala is located within the dry zone, in the district of Hambantota (Southern Province), where high malaria transmission prevailed in the past  with intensive insecticide usage for mosquito control [48, 49]. Based on the country’s annual malaria incidence reports [50,51,52,53,54], the highest number of malaria cases have been reported from the Hambantota district, consistently over a 5-year period prior to 2012 (when malaria was eliminated) . DDT was used in Sri Lanka for both agricultural insect pest control and malaria vector control, until 1975/1977 when its usage was banned due to environmental concerns and the appearance of widespread DDT resistance among malaria vectors [49, 55]. Organophosphates, mainly malathion, replaced DDT in mosquito control programmes and carbamates were introduced to the agricultural sector [49, 55]. Pyrethroids were first introduced to both health and agricultural sectors in mid-1990s [49, 55]. Such history of prolonged exposure to insecticides may be responsible for the observations made in the present study, especially the malathion resistance observed in the tested sand fly populations. It is noteworthy, that the calculations and cut-offs used in this study were based on the current recommendations and insecticide discriminating dosages described by the WHO for mosquitoes. This is due to the non-availability of either specific recommendations or resistance discriminating values for sand flies. It is tempting to speculate that the actual cut-offs for resistance discrimination in sand flies may be even lower than those given for mosquitoes, considering the smaller body mass and the softness of the cuticle of these tiny insects. However, the toxicity profile of an organism is considered species-specific and is known to be determined by a complex interplay of factors and therefore, no solid data exist to support such speculations.
It was noted that the knockdown effect, recorded as percent knockdown flies after one hour of insecticide exposure, increased with increasing concentrations of an insecticide. However, this was only up to a certain point of concentration of the insecticide, and never exceeded 75% at any of the tested concentrations, which included those that resulted in a mortality of 100% after a 24-hour recovery period. For Pannala and Thalawa populations, certain tested concentrations of DDT, malathion and propoxur showed zero knockdown after one hour of exposure, but a mortality of 100% after 24 hours. This apparent discrepancy between the knockdown concentrations and the LC values may be an early indication of the development of insecticide tolerance in these populations, as described previously .
According to the WHO guidelines, remaining activity of < 30% of propoxur-inhibited AChE target site is considered as the cut-off level for its sensitivity . Our study indicates 10–34.5% of the tested sand fly populations were heterozygous insensitive at the AChE target site that is associated with organophosphate and carbamate resistance. A previous study in Delft Island in Sri Lanka (a small island located close to the northern border of the country) has also shown a similar high prevalence of reduced sensitivity of the AChE target site in sand flies . Using the same biochemical assays Surendran et al.  have shown that 14% of the Delft island population of sand flies had > 80% residual activity of AChEs after propoxur inhibition.
Results of the biochemical assays indicated that GST and esterase activities of the populations were poorly elevated. All four populations showed no elevation of monooxygenases. Sand flies from Delft island of Sri Lanka have shown a slight elevation of esterase and GST but no elevation of monooxygenases . However, the impact of elevated esterases and insensitive AChE on the resistance status of sand flies was not assessed in the study of Surendran et al. . Based on the present observations, it might be assumed that the development of tolerance to all four groups of insecticides is ‘in progress’ in Sri Lankan sand fly populations. This might be attributed to the prolonged and widespread use of insecticides in both agricultural and health sector insect pest/vector control programmes [49, 55, 57]. If the study sites are lined up based on the malaria case burden in the relevant district during the 2007–2012 period, the order would be as: Mirigama < Pannala < Thalawa < Mamadala [11, 50,51,52,53,54]. This order is reversed apparently when the trend of insecticide susceptibility of sand flies in each study site (as demonstrated in this study) is considered. However, statistically there was no significant difference in mortality among the populations studied here for the tested insecticides.
Association of the L1014F kdr-type mutation in pyrethroid resistance was initially discovered in house flies and cockroaches [24, 25, 58,59,60]. Later, the same mutation was recorded from several other pyrethroid-resistant insects . In Sri Lanka, the presence of the L1014F kdr-type mutation has been previously reported the the mosquitoes Culex quinquefasciatus  and Anopheles subpictus  and the bed bug Cimex hemipterus . However, no previous studies have described the presence of kdr type mutations in sand flies in Sri Lanka, although L1014F and L1014S kdr-type mutations have been recorded in the sand fly population in India . It is interesting to note that the L1014F/L1014S mutation was present in a notable proportion of sand flies tested in the present study, with the highest mutation frequency evident in the Mamadala samples, presumably associated with the history of high pyrethroid spraying (N. Subasinghe, Regional Malaria Officer, Hambantota, Sri Lanka, personal communication).
Sri Lankan sand fly populations from the wild were in most part susceptible to the insecticides tested (except perhaps for malathion), when judged against the resistant discriminating dosages of mosquitoes. However, the less than optimal knockdown rates, the presence of marginally higher concentrations of insecticide detoxifying enzymes with the presence of the AChE target site insensitivity, relatively lower levels of insecticide susceptibility (except for deltamethrin) in flies originated from a location with a history of heavy insecticide usage and the presence of a known genetic mutation associated with insecticide resistance in a notable proportion of flies, could be viewed as a warning signal for the development of tolerance against insecticides.
Populations of P. argentipes in Sri Lanka seem to be susceptible to a broad range of insecticides. However, the early signs of insecticide tolerance that were evident in a country where strategic initiatives towards vector control of leishmaniasis are not yet in place, imply substantial challenges for future leishmaniasis control. Therefore, these findings, demand a regionally-coordinated strategic plan to address the apparent threat of insecticide tolerance in local sand fly populations. Such initiatives may include regulations imposed on the use of insecticides, particularly in the agricultural sector. The findings will also pave the way for more extensive investigations to aid a vector control strategy that may need to be adopted within a future national leishmaniasis control programme.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The datasets used and/or analyzed during the present study are available from the corresponding author upon reasonable request.
Center for Disease Control
- kdr :
- LC100 :
the lowest dose, which kills 100% of the tested population after 24 hours
- LC0 :
the dosage, which is not lethal to any flies in the population
- VGSC :
voltage gated sodium channel
World Health Organization
WHO. Leishmaniasis. 2019. Geneva: World Health Organization; 2019. https://www.who.int/gho/neglected_diseases/leishmaniasis/en/. Accessed 30 Aug 2019.
Leishmaniasis Desjeux P. Public health aspects and control. Clin Dermatol. 1996;14:417–23.
Karunaweera ND, Pratlong F, Siriwardane HVYD, Ihalamulla RL, Dedet JP. Sri Lankan cutaneous leishmaniasis is caused by Leishmania donovani zymodeme MON-37. Trans R Soc Trop Med Hyg. 2003;97:380–1.
Picado A, Dash AP, Bhattacharya S, Boelaert M. Vector control interventions for visceral leishmaniasis elimination initiative in South Asia, 2005–2010. Indian J Med Res. 2012;136:22–31.
Weekly Epidemiological reports: epidemiology unit. http://www.epid.gov.lk/. Accessed 30 Aug 2019.
Karunaweera ND. Leishmania donovani causing cutaneous leishmaniasis in Sri Lanka: a wolf in sheep’s clothing? Trends Parasitol. 2009;25:458–63.
Naotunna TD, Rajakulendran S, Abeywickreme W, Kulasiri CD, Perera J, Premaratne UN, et al. Cutaneous leishmaniasis in Sri Lanka. An imported disease linked to the Middle East and African employment boom. Trop Geogr Med. 1990;42:72–4.
Epidemiology unit. Weekly epidemiological report. 2019. http://www.epid.gov.lk/web/images/pdf/wer/2019/vol_46_no_01-english.pdf. Accessed 30 Aug 2019.
Alexander B, Maroli M. Control of phlebotomine sandflies. Med Vet Entomol. 2003;17:1–18.
Karunaratne SHPP, De Silva WAPP, Weeraratne TC, Surendran SN. Insecticide resistance in mosquitoes: development, mechanisms and monitoring. Ceylon J Sci. 2018;47:299–309.
Annual report 2016—Anti malaria campaign. 2016. http://www.malariacampaign.gov.lk/index.php/en/resources/annual-reports/78-annual-report-2016. Accessed 1 Apr 2018.
Surendran SN, Karunaratne SHPP, Adamsn Z, Hemingway J, Hawkes NJ. Molecular and biochemical characterization of a sand fly population from Sri Lanka: evidence for insecticide resistance due to altered esterases and insensitive acetylcholinesterase. Bull Entomol Res. 2005;95:371–80.
Bhattacharya SK, Sur D, Sinha PK, Karbwang J. Elimination of leishmaniasis (kala-azar) from the Indian subcontinent is technically feasible and operationally achievable. Indian J Med Res. 2006;123:195–6.
Dinesh DS, Das P, Picado A, Davies C, Speybroeck N, Ostyn B, et al. Long-lasting insecticidal nets fail at household level to reduce abundance of sand fly vector Phlebotomus argentipes in treated houses in Bihar (India). Trop Med Int Health. 2008;13:953–8.
Thakur CP. A new strategy for elimination of kala-azar from rural Bihar. Indian J Med Res. 2007;126:447–51.
Amalraj DD, Sivagnaname N, Srinivasan R. Susceptibility of Phlebotomus argentipes and P. papatasi (Diptera: Psychodidae) to insecticides. J Commun Dis. 1999;31:177–80.
Rahman M. Insecticide substitutes for DDT to control mosquitoes may be causes of several diseases. Environ Sci Pollut Res Int. 2013;20:2064–9.
Joshi AB, Bhatt LR, Regmi S, Ashford R. An assessment of the effectiveness of insecticide spray in the control of visceral leishmaniasis in Nepal. J Nepal Health Res Council. 2003;1:1–6.
Epidemiology and Disease Control Division (EDCD). The annual internal assessment of malaria and kala-azar control activities 1999. Department of Health Services. Epidemiology and Disease Control Division and Vector Borne Disease Research and Training Centre. Ministry of Health, HMG/N, Nepal; 2000. https://www.worldcat.org/title/annual-internal-assessment-of-malaria-and-kala-azar-control-activities/oclc/607726283. Accessed 30 Aug 2019.
Government of the People’s Republic of Bangladesh. Ministry of Health and Family Welfare. Health Bulletin 2018. Management information system. Directorate General of Health Services, Mohakhali, Dhaka 1212; 2018. https://dghs.gov.bd/images/docs/Publicaations/HB%202018%20final.pdf. Accessed 16 Dec 2019.
Chowdhury R, Das ML, Chowdhury V, Roy L, Faria S, Priyanka J, et al. Susceptibility of field-collected Phlebotomus argentipes (Diptera: Psychodidae) sand flies from Bangladesh and Nepal to different insecticides. Parasites Vectors. 2018;11:336.
Hemingway J, Hawkes NJ, Mccarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34:653–65.
Soderlund DM, Bloomquist JR. Molecular mechanisms of insecticide resistance. In: Roush RT, Tabashnik BE, editors. Pesticide resistance in arthropods. New York: Chapman and Hall; 1990. p. 58–96.
Ingles PJ, Adams PM, Knipple DC, Soderlund DM. Characterization of voltage-sensitive sodium channel gene coding sequences from insecticide-susceptible and knockdown-resistant house fly strains. Insect Biochem Mol Biol. 1996;26:319–26.
Williamson M, Martinez-Torres D, Hick C, Devonshire A. Identification of mutations in the housefly para-type sodium channel gene associated with knockdown resistance (kdr) to pyrethroid insecticides. Mol Gen Genet. 1996;252:51–60.
Rinkevich FD, Yuzhe D, Dong K. Diversity and convergence of sodium channel mutations involved in resistance to pyrethroids. Pestic Biochem Physiol. 2013;106:93–100.
Amelia-yap ZH, Chen CD, Sofian-azirun M, Low VL. Pyrethroid resistance in the dengue vector Aedes aegypti in Southeast Asia: present situation and prospects for management. Parasites Vectors. 2018;11:332.
Gomes B, Purkait B, Deb RM, Rama A, Singh RP, Foster GM, et al. Knockdown resistance mutations predict DDT resistance and pyrethroid tolerance in the visceral leishmaniasis vector Phlebotomus argentipes. PLoS Negl Trop Dis. 2017;11:e0005504.
Senanayake SASC, Abeyewicreme W, Dotson EM, Karunaweera ND. Characteristics of phlebotomine sandflies in selected areas of Sri Lanka. Southeast Asian J Trop Med Public Health. 2015;46:994–1004.
Climatic zones in Sri Lanka: climatic change Secretariat Sri Lanka. Ministry of Mahaweli Development and Environment. http://www.climatechange.lk. Accessed 30 Aug 2019.
WHO. Part II tutor’s guide. In: Entomological field techniques for malaria control. Geneva: World Health Organization; 1992. https://www.who.int/malaria/publications/atoz/9241544392/en/. Accessed 30 Aug 2019.
Volf P, Volfova V. Establishment and maintenance of sand fly colonies. J Vector Ecol. 2011;36(Suppl. 1):S1–9.
Lewis DJ. The phlebotomine sand flies (Diptera: Psychodidae) of oriental region. Bull Br Mus (Nat Hist) Entomol. 1978;37:217–343.
Ilango K. A taxonomic reassessment of the Phlebotomus argentipes species complex (Diptera: Psychodidae: Phlebotominae). J Med Entomol. 2010;47:1–15.
Bandara KMUJ, Karunaratne SHPP. Mechanisms of acaricide resistance in the cattle tick Rhipicephalus (Boophilus) microplus in Sri Lanka. Pestic Biochem Physiol. 2017;139:68–72.
WHO. Insecticide resistance and vector control 1963. Geneva: World Health Organization; 1963. https://apps.who.int/iris/bitstream/handle/10665/40568/WHO_TRS_265_%28part1%29.pdf?sequence=1&isAllowed=y. Accessed 1 Apr 2018.
WHO. Guidelines for efficacy testing of insecticides for indoor and outdoor ground-applied space spray applications 2009. Geneva: World Health Organization; 2009. http://apps.who.int/iris/bitstream/10665/70070/1/WHO_HTM_NTD_WHOPES_2009.2_eng.pdf?ua=1. Accessed 1 Apr 2018.
Abbott WS. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18:265–7.
WHO. Techniques to detect insecticide resistance mechanisms (field and laboratory manual) 1998. Geneva: World Health Organization; 1998. http://whqlibdoc.who.int/hq/1998/WHO_CDS_CPC_MAL_98.6.pdf?ua=1. Accessed 1 Apr 2018.
Rijal S, Sundar S, Mondal D, Das P, Alvar J, Boelaert M. Eliminating visceral leishmaniasis in South Asia: the road ahead. BMJ. 2019;364:k5224.
Dinesh DS, Das ML, Picado A, Roy L, Rijal S, Singh SP, et al. Insecticide susceptibility of Phlebotomus argentipes in visceral leishmaniasis endemic districts in India and Nepal. PLoS Negl Trop Dis. 2010;4:e859.
Faraj C, Ouahabi S, Adlaoui EB, Elkohli ME, Lakraa L, Rhazi ME, et al. Insecticide susceptibility status of Phlebotomus (Paraphlebotomus) sergenti and Phlebotomus (Phlebotomus) papatasi in endemic foci of cutaneous leishmaniasis in Morocco. Parasites Vectors. 2012;5:51.
Hassan MM, Widaa SO, Osman OM, Numiary MSM, Ibrahim MA, Abushama HM. Insecticide resistance in the sand fly, Phlebotomus papatasi from Khartoum State, Sudan. Parasites Vectors. 2012;5:46.
WHO. Test procedures for insecticide resistance, monitoring in malaria vectors, bio-efficacy and persistence of insecticides on treated surfaces. Report of the WHO informal consultation 1998. Geneva: World Health Organization; 1998. https://apps.who.int/iris/bitstream/handle/10665/64879/WHO_CDS_CPC_MAL_98.12.pdf?sequence=1&isAllowed=y. Accessed 1 Apr 2018.
WHO. Guidelines for testing mosquito adulticides for indoor residual spraying and treatment of mosquito nets 2006. Geneva: World Health Organization; 2006. https://apps.who.int/iris/bitstream/handle/10665/69296/WHO_CDS_NTD_WHOPES_GCDPP_2006.3_eng.pdf. Accessed 1 Apr 2018.
WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes 2016. Geneva: World Health Organization; 2016. https://apps.who.int/iris/rest/bitstreams/1063682/retrieve. Accessed 1 Apr 2019.
Karunaweera ND, Galappaththy GN, Wirth DF. On the road to eliminate malaria in Sri Lanka: lessons from history, challenges, gaps in knowledge and research needs. Malar J. 2014;13:59.
Wickramasinghe MB. Malaria and its control in Sri Lanka. Ceylon Med J. 1981;26:107–15.
Konradsen F, Amerasinghe FP, Hoek W Van der, Amerasighe PH. Malaria in Sri Lanka: current knowledge on transmission and control. International Water management Institute: Colombo, Sri Lanka; 2000. http://publications.iwmi.org/pdf/H027692.pdf. Accessed 30 Aug 2019.
Annual reports. Anti malaria campaign. 2007. http://www.malariacampaign.gov.lk/images/Publication%20Repository/Annual%20Reports/Annual%20Report%202007.pdf. Accessed 1 Apr 2018.
Annual reports. Anti malaria campaign. 2008. http://www.malariacampaign.gov.lk/images/Publication%20Repository/Annual%20Reports/Annual%20Report%202008.pdf. Accessed 1 Apr 2018.
Annual reports. Anti malaria campaign. 2009. http://www.malariacampaign.gov.lk/images/Publication%20Repository/Annual%20Reports/Annual%20Report%202009.pdf. Accessed 1 Apr 2018.
Annual reports. Anti malaria campaign. 2011. http://www.malariacampaign.gov.lk/images/Publication%20Repository/Annual%20Reports/Annual%20Report%202011.pdf. Accessed 1 Apr 2018.
Annual reports. Anti malaria campaign. 2012. http://www.malariacampaign.gov.lk/images/Publication%20Repository/Annual%20Reports/Annual%20Report%202012.pdf. Accessed 1 Apr 2018.
WHO. Malaria free Sri Lanka 2017. Geneva: World Health Organization; 2017. http://apps.searo.who.int/PDS_DOCS/B5395.pdf. Accessed 1 Apr 2018.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Bergé JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.
Senewiratne B, Thambilillai S. Pattern of poisoning in a developing agricultural country. Br J Prev Soc Med. 1974;28:32–6.
Williamson MS, Denholm I, Bell CA, Devonshire AL. Knockdown resistance (kdr) to DDT and pyrethroid insecticides maps to a sodium channel gene locus in the housefly (Musca domestica). Mol Gen Genet. 1993;240:17–22.
Knipple DC, Doyle KE, Marsella-Herrick PA, Soderlund DM. Tight genetic linkage between the kdr insecticide resistance trait and a voltage-sensitive sodium channel gene in the house fly. Proc Natl Acad Sci USA. 1994;91:2483–7.
Miyazaki M, Ohyama K, Dunlap D, Matsumura F. Cloning and sequencing of the para -type sodium channel gene from susceptible and kdr-resistant German cockroaches (Blattella germanica) and house fly (Musca domestica). Mol Gen Genet. 1996;252:61–8.
Wondji CS, De Silva WAPP, Hemingway J, Ranson H, Parakrama SHPP. Characterization of knockdown resistance in DDT- and pyrethroid-resistant Culex quinquefasciatus populations from Sri Lanka. Trop Med Int Health. 2008;13:548–55.
Karunaratne SHPP, Hawkes NJ, Perera MDB, Ranson H, Hemingway J. Mutated sodium channel genes and elevated monooxygenases are found in pyrethroid resistant populations of Sri Lankan malaria vectors. Pestic Biochem Physiol. 2007;88:108–13.
Punchiheva R, De Silva WAPP, Weeraratne TC, Karunaratne SHPP. Insecticide resistance mechanisms with novel ‘kdr’ type gene mutations in the tropical bed bug Cimex hemipterus. Parasit Vectors. 2019;12:310.
Perera MDB, Hemingway J, Karunaratne SHPP. Multiple insecticide resistance mechanisms involving metabolic changes and insensitive target sites selected in anopheline vectors of malaria in Sri Lanka. Malar J. 2008;7:168.
We acknowledge Ms Janadaree Bandara, Universirty of Peradeniya, for providing the insecticide-impregnated papers for the bioassay tests, Dr Thilini Weeraratne and Ms Sumudu Samarasinghe for assistance with molecular analysis, Mr M. P. Ariyapala and Mr. Sunil Shantha for field assistance and rearing of sand flies, Mr Wasantha Senadeera for providing a map of study locations, Mr Siraj for assistance with formatting figures, Ministry of Health officials, especially Regional Malaria Officers and Entomological officers in Hambantota, Anuradhapura and Kurunegala and Mirigama for field assistance, the Head and staff of the Department of Parasitology, Faculty of Medicine, University of Colombo for logistical help.
Research reported in this publication was supported by the National Institute of Allergy and Infectious Diseases of the National Institutes of Health, USA, under award number U01AI136033. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Statistical analysis to test for significant differences of sand fly mortality between the studied populations to tested insecticides using the Wilcoxon signed-rank test.
Statistical analysis to test for significant differences using ANOVA of elevated enzymes and remaining enzyme activities in each location.
Sequence alignment of the VGSC gene in P. argentipes with the sequence for the voltage-sensitive sodium channel of Musca domestica (GenBank: AAB47604). Arrow indicates the position of the L1014F mutation. Heterozygotes samples were indicated with an asterisk (*).
About this article
Cite this article
Pathirage, D.R.K., Karunaratne, S.H.P.P., Senanayake, S.C. et al. Insecticide susceptibility of the sand fly leishmaniasis vector Phlebotomus argentipes in Sri Lanka. Parasites Vectors 13, 246 (2020). https://doi.org/10.1186/s13071-020-04117-y