Susceptibility status of the wild-caught Phlebotomus argentipes (Diptera: Psychodidae: Phlebotominae), the sand fly vector of visceral leishmaniasis, to different insecticides in Nepal

Abstract

Background

Visceral leishmaniasis (VL) is targeted for elimination as a public health problem in Nepal by 2023. For nearly three decades, the core vector control intervention in Nepal has been indoor residual spraying (IRS) with pyrethroids. Considering the long-term use of pyrethroids and the possible development of resistance in the vector Phlebotomus argentipes sand flies, we monitored the susceptibility status of their field populations to the insecticides of different classes, in villages with and without IRS activities in recent years.

Methodology/Principal findings

Sand flies were collected from villages with and without IRS in five VL endemic districts from August 2019 to November 2020. The WHO susceptibility test procedure was adopted using filter papers impregnated at the discriminating concentrations of insecticides of the following classes: pyrethroids (alpha-cypermethrin 0.05%, deltamethrin 0.05%, and lambda-cyhalothrin 0.05%), carbamates (bendiocarb 0.1%) and organophosphates (malathion 5%). Pyrethroid resistance intensity bioassays with papers impregnated with 5× of the discriminating concentrations, piperonyl butoxide (PBO) synergist-pyrethroid bioassays, and DDT cross-resistance bioassays were also performed. In the IRS villages, the vector sand flies were resistant (mortality rate <90%) to alpha-cypermethrin and possibly resistant (mortality rate 90–97%) to deltamethrin and lambda-cyhalothrin, while susceptibility to these insecticides was variable in the non-IRS villages. The vector was fully susceptible to bendiocarb and malathion in all villages. A delayed knockdown time (KDT50) with pyrethroids was observed in all villages. The pyrethroid resistance intensity was low, and the susceptibility improved at 5× of the discriminating concentrations. Enhanced pyrethroid susceptibility after pre-exposure to PBO and the DDT-pyrethroid cross-resistance were evident.

Conclusions/Significance

Our investigation showed that P. argentipes sand flies have emerged with pyrethroid resistance, suggesting the need to switch to alternative classes of insecticides such as organophosphates for IRS. We strongly recommend the regular and systematic monitoring of insecticide resistance in sand flies to optimize the efficiency of vector control interventions to sustain VL elimination efforts in Nepal.

Introduction

Indoor residual spraying (IRS) is the core vector control intervention for the elimination of visceral leishmaniasis (VL) in South-East Asia, where about 200 million people are at risk of the disease [1,2]. VL, also called kala-azar in the Indian subcontinent, is a vector-borne neglected tropical disease which is fatal if not treated timely, and mostly affects the lower socio-economic classes of the communities [3]. The causative agent of the disease is the protozoan parasite Leishmania donovani Laveran & Mesnil, 1903 (Kinetoplastea: Trypanosomatida) and is transmitted between human hosts by female sand flies of the species Phlebotomus argentipes Annandale & Brunetti, 1908 (Diptera: Psychodidae: Phlebotominae) [4]. IRS has been shown to reduce the P. argentipes density in Nepal [5–7] and is a key intervention in the country’s VL elimination strategy. Since 1992, mostly synthetic pyrethroids i.e. alpha-cypermethrin, deltamethrin, or lambda-cyhalothrin have been sprayed in houses [8], after the use of dichloro-diphenyl-trichloro-ethane (DDT) was banned in the country due to its health and environmental hazards [9]. The choice of synthetic pyrethroids over other classes of insecticides for IRS was based on the combination of low cost and high efficacy of these insecticides [2]. However, there are concerns that the use of a single class of insecticides for a long period might have potentially selected resistance in sand fly vector and could undermine the effectiveness of vector control in Nepal.

At present, little information is available on the current susceptibility status of P. argentipes to the pyrethroid insecticides used for IRS in Nepal. Previously reported DDT resistance in P. argentipes populations [10] provided a clue for possible cross-resistance to pyrethroids due to the involvement of knockdown resistance (kdr) mutations in a common para-type voltage-gated sodium channel (vgsc) target site gene in nerve cells [11]. Early signs of pyrethroid resistance were observed in the Eastern part of Nepal and also in some Indian areas bordering the country [11,12].

Systematic monitoring of insecticide resistance is considered a key component of VL surveillance. The choice of an appropriate insecticide needs to be dynamic and adapted to the changing insecticide resistance status of the vectors if IRS interventions are to be successful [13]. Besides monitoring resistance to the currently used pyrethroids, it is therefore equally essential to understand the underlying resistance mechanisms and assess the susceptibility status against insecticides of alternative classes with unrelated modes of action, in order to select the most efficacious insecticides and allow for informed decision-making.

With this study, we aimed to provide an update on the susceptibility status of wild-caught P. argentipes sand flies to pyrethroid insecticides currently used for IRS as well as to insecticides from two alternative classes with unrelated modes of action (i.e., carbamate and organophosphate). Additionally, we assessed the strength of resistance in pyrethroid resistant phenotypes of sand flies with papers impregnated with five times (5×) the discriminating concentrations of pyrethroids, i.e., alpha-cypermethrin, deltamethrin, and lambda-cyhalothrin. We also explored possible underlying metabolic resistance mechanisms by piperonyl butoxide (PBO) synergist-pyrethroid bioassays, and by cross-resistance studies with DDT to ascertain the involvement of kdr mutations in the pyrethroid resistant phenotypes.

Methods

Ethics statement

Ethical clearance for the study was obtained from the Institutional Ethical Review Board of BP Koirala Institute of Health Sciences, Dharan, Nepal (Registration no.- IRC/1729/019). Informed verbal consent was obtained from the household and cattle shed owners where sand flies were collected.

Study sites

The study was carried out in five VL endemic districts in Eastern and Central Nepal, namely Morang, Saptari, Siraha, Mahottari, and Sarlahi (Fig 1). These districts are located at 26°23′–26°55′ N latitude and 85°37′–87°21′ E longitude in the Southern lowland plains of Nepal, which is also known as the ‘Terai’ region. The area has a tropical savanna climate with a mean annual temperature of 19–31°C and a mean annual rainfall of 1600–2500 mm [14,15]. These districts were selected based on (i) a high average annual VL incidence reported in the years 2016–2018 and (ii) accessibility by road. In each district, two villages were selected; one intervention village, which received regular or intermittent IRS during the last three or more years due to the presence of VL cases (subsequently called IRS village), and one non-intervention or control village, which did not receive IRS in the three years before 2019 as no VL case was reported (subsequently called non-IRS village). The selected intervention and control villages were Majhare and Pokhariya in Morang district, Belhi Chapena and Gamhariya Parwaha in Saptari district, Bishanpur and Dahipaudi in Siraha district, Madhepura and Anarban in Mahottari district, and Riterkhor and Lodba in Sarlahi district respectively.

[Figure omitted. See PDF.]

Fig 1. Map of Nepal showing visceral leishmaniasis endemicity and the locations of the study districts.

The map was produced with QGIS (version 3.10) with open access shapefile (New political and administrative boundaries Shapefile of Nepal—Datasets—Open Data Nepal).

https://doi.org/10.1371/journal.pntd.0010304.g001

Sand fly collection

Sand flies resting in human dwellings, cattle sheds, and dwellings shared by both humans and cattle (mixed dwellings) in the selected study villages were collected from August 2019 to November 2020. Collections were done in the early morning between 04:30 and 06:30 hours using a mouth aspirator and a torch. About 5–10 households per village were visited by four insect collectors to collect enough sand flies to perform at least one set of bioassays in a day. A temporary laboratory was set up in each village to avoid long distance transportation, which can stress or damage the captured sand flies.

Bioassays

The WHO bioassay test kit and insecticide-impregnated papers were obtained from Universiti Sains Malaysia, Penang, Malaysia. The susceptibility tests with all insecticides and intensity bioassays with pyrethroids were carried out in two phases following the standard WHO susceptibility test procedure [16] (Fig 2). The usual test procedures for mosquitoes were adopted for sand flies. In Phase I of the study, the baseline susceptibility status of wild-caught female P. argentipes was assessed in all the selected villages with filter papers impregnated at the discriminating concentrations of insecticides of the following classes: pyrethroids (alpha-cypermethrin 0.05%, deltamethrin 0.05%, and lambda-cyhalothrin 0.05%), carbamates (bendiocarb 0.1%) and organophosphates (malathion 5%). Silicone oil and olive oil impregnated papers were used as control papers for pyrethroid and organophosphate/carbamate bioassays respectively.

[Figure omitted. See PDF.]

Fig 2. A graphic representation of susceptibility test sequence and probable outcomes for insecticide resistance monitoring in Phlebotomus argentipes wild populations (adopted from the WHO test procedures [16]).

https://doi.org/10.1371/journal.pntd.0010304.g002

Before exposure, the sand flies were allowed to settle down for one hour in holding tubes lined with insecticide-free clean white papers, each tube holding 20–25 sand flies. They were then transferred from the holding tubes to the exposure tubes lined with insecticide-impregnated papers. Each test for a given insecticide-impregnated paper comprised of five replicates of 20–25 sand flies, while each control test in parallel comprised of one control tube with a similar number of sand flies. In case of insufficient numbers of sand flies collected i.e., <120 sand flies per set, the test was repeated on subsequent days until 100 or more sand flies had been included in the five replicates, in addition to one control group for each batch of the test.

The exposure time was fixed to 60 minutes. Knockdown of the sand flies was recorded at 10, 15, 20, 30, 40, 50, and 60 minutes after the start of the exposure with insecticides. After 1h of the exposure, the sand flies were transferred back to the holding tubes, which were kept standing upright with the screen end up, provided with 10% glucose-soaked cotton wool pads on the fine mesh screen, and placed in an area with diffused illumination for the next 24h recovery period. On hot and dry days, the holding tubes were placed in a container box covered with wet towels to stabilize the ambient temperature below 30°C and the relative humidity at 80% ± 10%. The numbers of dead and alive sand flies from each of the exposure and control holding tubes were counted and recorded on a pre-printed form at the end of the recovery period. The species of the sand flies were identified after completing tests with the help of a stereomicroscope, regional taxonomic keys, and other resources [17–19]. The physiological stages of the sand flies were also recorded. Only female P. argentipes sand flies were considered for mortality assessment. The following criteria were used to classify phenotypic resistance frequency in wild sand flies: fully susceptible (mortality rate ≥98%), possible resistance (mortality rate 90–97%), and confirmed resistance (mortality rate <90%) (Fig 2) [16].

In Phase II of the study, the possibly resistant sand fly populations observed during Phase I were reconfirmed for their resistance status. In such repeat bioassays, the mortality rate of <98% was interpreted as confirmed resistance in line with WHO test procedures [16]. The population with confirmed resistance in both Phases I and II were subjected to three additional tests (Fig 2), as follows:

1. The pyrethroid impregnated papers with 5× discriminating concentrations were used to determine resistance intensity.

2. The synergist-insecticide bioassays were performed using PBO 4% impregnated papers in which one set of tests comprised of at least four types of exposures and a batch of 20–25 sand flies for each exposure group; (i) PBO alone–the sand flies were exposed to PBO 4% impregnated paper for 60 minutes, (ii) PBO followed by insecticide–the sand flies were first exposed to PBO 4% impregnated papers for 60 minutes in tubes and then exposed to papers impregnated with either alpha-cypermethrin 0.05%, deltamethrin 0.05%, or lambda-cyhalothrin 0.05% for another 60 minutes, (iii) insecticide alone–the sand flies were exposed to alpha-cypermethrin 0.05%, deltamethrin 0.05%, or lambda-cyhalothrin 0.05% impregnated papers for 60 minutes, and (iv) control–the sand flies were exposed to silicone oil impregnated paper for 60 minutes. A test was considered invalid in the case of >10% mortality in the ‘PBO alone’ exposure tube. The interpretation of the synergist-insecticide bioassays was done as follows: (a) the effect of PBO could not be reliably assessed if the mean mortality in tests with ‘insecticide alone’ papers was ≥90%, (b) a monooxygenase-based resistance mechanism was fully accountable for expression of resistance in the test population if the mean mortality in ‘insecticide alone’ tests was <90% and the mean mortality in the ‘PBO followed by insecticide’ tests was ≥98%, and (c) a monooxygenase-based resistance mechanism was partially accountable for the expression of resistance if the mean mortality in the ‘insecticide alone’ tests was <90% and the mean mortality in the ‘PBO followed by insecticide’ tests was more than the mean mortality in the ‘insecticide alone’ tests but <98%.

3. Bioassays with DDT 4% papers along with control were conducted to detect possible cross-resistance to pyrethroids. Risella oil impregnated papers were used in control. The procedures, numbers of replicates and interpretation of mortality percentages of resistance intensity and DDT bioassays were the same as those followed in Phase I.

Data management and statistical analysis

The number of sand flies found dead in all the test tubes was counted for each replicate and the mortality rate was expressed as a percentage of the dead among the total number of the exposed sand flies. A similar calculation was done for the control tubes. In case the mortality rate in the control tubes was between 5% to <20%, the mortality rate in the test replicates was corrected using Abbott’s formula [20]. Corrected mortality rates were interpreted to determine susceptibility status as described in the WHO guidelines [16] (Fig 2). The tests were discarded if the mortality rate in the control tube was ≥20%.

R software version 4.1.0 [21] was used to conduct all statistical analyses and plotting bar graphs with their mean and standard errors representing variations between the test replicates. For the main statistical analysis of the sand fly mortality, a Poisson regression model was chosen with the number of dead sand flies as the outcome and the total number of the sand flies included in the tests as the exposure. To account for clustering of observations that were done in the same district, or the same village, or that were part of the same susceptibility test, we opted for a Generalized Linear Mixed Model (GLMM) containing random effects for these three levels. In practice, this yielded boundary estimates for the variances of the random effects, which casted doubt on the stability of the model. We used a Bayesian approach for this analysis to better account for the uncertainty of the estimates using the R package “rstanarm” [22,23] with default, weakly informative priors. Central credible intervals at 95% were reported for each coefficient and these were interpreted based on the weight of the evidence rather than whether they exceeded 1. Additionally, a probit analysis [24–26] was performed on the R package “ecotox” [27] to estimate the time and fiducial confidence limits at which 50% of the sand fly population were expected to be knocked down (KDT50) after exposure to the insecticide-impregnated papers.

Results

Of 10,564 female P. argentipes collected from the study villages to perform the bioassays during both phases, 8,009 were used in the test replicates and the remaining in the control. Most of the sand flies were blood-fed (>60%), followed by gravid and unfed physiological stages (S1 Fig).

Phase I bioassays

In the IRS villages, the sand fly populations showed resistance to alpha-cypermethrin 0.05% in four villages and possible resistance in one village, with mortality rates between 69.28% ± 3.82% and 92.77% ± 3.01%. The sand flies from four villages showed possible resistance to deltamethrin 0.05% with mortality rates between 95.44% ± 0.12% to 97.95% ± 1.26% and lambda-cyhalothrin 0.05% with mortality rates between 90.12% ± 4.73% to 95.70% ± 2.16%, while full susceptibility was registered in the other remaining villages. The sand fly populations were fully susceptible to bendiocarb 0.1% and malathion 5% (mortality rate >99%) in all the five IRS villages (Fig 3).

[Figure omitted. See PDF.]

Fig 3.

Mean mortality percentage after 24h of holding in the IRS (A) and the non-IRS villages (B). Error bars represent the standard error of the mean. The green-dashed intercept represents the level of 98% mortality; the red dashed intercept represents the level of 90% mortality.

https://doi.org/10.1371/journal.pntd.0010304.g003

In the non-IRS villages, the sand fly populations were resistant to alpha-cypermethrin 0.05% with mortality rate 83.08% ± 2.66% in one village, possibly resistant in three villages with mortality rates between 90.77% ± 3.07% to 95.36% ± 2.24%, and susceptible in one village. All villages except one had vector populations susceptible to deltamethrin 0.05%. In three villages, the sand fly populations showed possible resistance to lambda-cyhalothrin 0.05% with mortality rates 95.04% ± 0.09% to 96.69% ± 1.72%, whereas in two others they were susceptible. Like the IRS villages, the sand flies in all non-IRS villages were susceptible to bendiocarb 0.1% and malathion 5%, with a 100% mortality rate.

Mortality rates were further assessed in the regression model. Overall, mortality in the sand flies was lower in the IRS villages than in the non-IRS villages (IRR = 0.9). While taking into account physiological stages, blood-fed and gravid vector sand flies had slightly higher mortalities than unfed ones (IRR >1). Mortality rates were higher (IRR >1) for other insecticides viz. deltamethrin, lambda-cyhalothrin, bendiocarb and malathion than for alpha-cypermethrin. The difference between alpha-cypermethrin and other insecticides in terms of mortality was larger in the sprayed villages than in the unsprayed ones (interaction effect IRR >1). There was little variation between the districts and between the villages, which was not already explained by variables in the model (standard deviations of random effects near 0, less than 0.001) (Table 1).

[Figure omitted. See PDF.]

Table 1. Bayesian Poisson regression of the sand fly mortality rate in previously sprayed villages, sand fly subtype, and insecticide used in the susceptibility test, with interaction effects between village spraying and insecticide used, and 95% central credible intervals.

https://doi.org/10.1371/journal.pntd.0010304.t001

Proportion of not knocked down and KDT50.

There was no significant difference in the proportion of sand flies that were not knocked down after 1h exposure to insecticides in the IRS and the non-IRS villages (Fig 4). The estimated time for knockdown of 50% of sand flies (KDT50) upon exposure to pyrethroids at the discriminating concentrations (1×) varied between 31.27 to 38.64 min for the IRS villages and 29.73 to 34.86 min for the non-IRS villages. The KDT50 was four times shorter (8.70–11.27 min) at 5× concentrations of pyrethroids. Higher KDT50 (45.27 min) was observed in tests with DDT 4% (Table 2).

[Figure omitted. See PDF.]

Fig 4. Time until knockdown of P. argentipes with different insecticides tested in Phase I in the IRS and the non-IRS villages.

https://doi.org/10.1371/journal.pntd.0010304.g004

[Figure omitted. See PDF.]

Table 2. Estimated 50% knockdown times of wild populations of P. argentipes sand flies exposed to various insecticide papers impregnated at the 1× and the 5× discriminating concentrations in the IRS and the non-IRS villages.

https://doi.org/10.1371/journal.pntd.0010304.t002

Phase II bioassays

Repeat bioassays.

The bioassays were repeated in the villages where sand flies were found to be possibly resistant to pyrethroids during Phase I. In all these villages, the presence of resistant phenotypes was confirmed, with mortality rates in all repeated bioassays below 98% (Table 3).

[Figure omitted. See PDF.]

Table 3. Assessment of mean mortality percentage of P. argentipes to discriminating concentrations of pyrethroids in bioassays in selected villages with possible resistant vector populations detected in Phase I.

https://doi.org/10.1371/journal.pntd.0010304.t003

Resistance intensity bioassays.

The vector population from Bishanpur with confirmed resistance at 1× concentrations of all three pyrethroids showed 100% mortality at 5× of the discriminating concentrations (Table 4). As the mortality rate at 5× concentrations was 100%, bioassays at 10× concentration were not required.

[Figure omitted. See PDF.]

Table 4. Assessment of mean mortality percentage of P. argentipes to 5× discriminating concentrations of pyrethroids in a selected village with resistant sand fly populations in Phase I bioassays.

https://doi.org/10.1371/journal.pntd.0010304.t004

Synergist bioassays.

The mortality rate of the sand flies that were pre-exposed to PBO 4% was higher than that of those exposed to a given insecticide only. However, the synergistic effect of PBO to alpha-cypermethrin 0.05% could not be interpreted for the sand fly population in Belhichapena and to lambda-cyhalothrin in Bishanpur as the mean mortality rate in the tests with ‘insecticide alone’ was ≥90%. The susceptibility level was partially restored in Bishanpur and Riterkhor populations for alpha-cypermethrin 0.05%. Interestingly, a complete restoration of susceptibility was seen in the Bishanpur population for deltamethrin 0.05% (Fig 5).

[Figure omitted. See PDF.]

Fig 5. Comparative mortality percentage of P. argentipes to pyrethroid insecticides alone or pre-exposed to PBO 4% papers followed by pyrethroid insecticides.

The error bars represent the standard error of the mean. The green-dashed intercept represents the level of 98% mortality; the red dashed intercept represents the level of 90% mortality.

https://doi.org/10.1371/journal.pntd.0010304.g005

Bioassays with DDT 4%.

The mortality rate in the pyrethroid resistant P. argentipes populations exposed to DDT 4% impregnated papers ranged between 55.14% ± 1.49% and 76.15% ± 3.51%. Mortality percentage of pyrethroids in the corresponding villages were included from Phase I observations for the assessment (Fig 6).

[Figure omitted. See PDF.]

Fig 6. DDT-pyrethroid cross-resistance in P. argentipes from the selected IRS and non-IRS villages.

The error bars represent the standard error of the mean. The green-dashed intercept represents the level of 98% mortality; the red dashed intercept represents the level of 90% mortality.

https://doi.org/10.1371/journal.pntd.0010304.g006

Discussion

In our study, we found that pyrethroid resistance in the P. argentipes populations is well established in all the IRS villages, particularly so for alpha-cypermethrin 0.05% which has already been used for more than 15 years now in the IRS campaigns in Nepal. Moreover, pyrethroid resistance was also seen in the non-IRS villages. Development of pyrethroid resistance in the vector populations was also supported by a similar KDT50 in both the IRS and the non-IRS villages. Nevertheless, vector populations in all the study villages were fully susceptible to the selected carbamate and organophosphate insecticides. The resistance intensity bioassays confirmed the intensity of pyrethroid resistance is low, suggesting that a possible operational failure of pyrethroid-based vector control is unlikely in the immediate future, although it signals to implement an insecticide resistance management approach. Incomplete restoration of susceptibility in PBO synergist bioassays was indicative of a partial involvement of monooxygenase-based metabolic resistance mechanism. Resistance to DDT in pyrethroid-resistant populations is indicative of the presence of kdr mutations in the target site gene (vgsc) of the vector sand flies. As bioassays were conducted using wild captured sand flies, blood-fed stages outnumbered gravid and unfed in both categories of the villages (IRS and non-IRS), exposures (test and control) and outcomes (dead and alive) with expected mortalities slightly higher in blood-fed than the other two physiological stages.

Monitoring the susceptibility of wild populations of P. argentipes is essential to sustain the effectiveness of insecticides, as this can guide insecticide resistance management. While sporadic information on the susceptibility status of P. argentipes in Nepal were available over the last 25 years, resistance management has not been pursued with priority within the national VL vector control programme framework. The first known investigation in Nepal, conducted in Thilla village of Dhanusha district, dates back to 1997 and reported that P. argentipes populations were susceptible to lambda-cyhalothrin 0.1% and deltamethrin 0.25% [28]. In subsequent publications until 2001, the P. argentipes populations were reported to be susceptible to DDT 4%, permethrin 0.25%, deltamethrin 0.25%, lambda-cyhalothrin 0.1%, malathion 5% and bendiocarb 0.1% [28,29]. In 2009, bioassays conducted in four villages from Eastern Nepal revealed the presence of DDT resistance in one of the four populations of P. argentipes, as well as possible resistance to deltamethrin 0.05%, thereby providing an early sign of development of selection pressure against pyrethroids [10]. Bioassays in three districts of Eastern Nepal in 2015 and 2016, on the other hand, showed susceptible vector populations at 60 min exposure to alpha-cypermethrin 0.05% and deltamethrin 0.05% [12], which is contrary to the results obtained in our study.

Most of the VL endemic districts in Nepal share close borders with VL endemic districts of Bihar state in India, where DDT-based IRS was replaced with a pyrethroid in 2015 due to growing evidence of DDT resistance in vector sand flies [11]. A similar situation was experienced in the bordering part of Nepal in a previous study [10] as well as in our current study. Since then, several reports have been published reporting phenotypic resistance to pyrethroids [30–33]. Recently, the presence of kdr mutations in vector populations in some parts of India hinted at the emergence of pyrethroid resistance at molecular level [11,30]. Although P. argentipes populations were found susceptible to malathion in several locations within South-East Asia [30,31,34,35], resistant populations have also been detected [35]. In a recent study in three states of India, P. argentipes populations showed full susceptibility to all insecticides except for DDT [36]. In Bangladesh, another VL endemic country neighboring to Nepal, the P. argentipes populations have shown to be susceptible to insecticides in major classes used in IRS till date [12]. Only a few hundred kilometers apart, information on the susceptibility status of P. argentipes in VL endemic areas in India and Bangladesh gives insight into the necessity of such investigations in Nepal.

The delayed knockdown time and reduced mortality detected in our study illustrate the onset of pyrethroid resistance in vector populations. Both the IRS and the non-IRS villages showed similar knockdown times in vector populations, suggesting the spread of pyrethroid resistance in endemic areas irrespective of IRS interventions in recent years. Other causes contributing to the development of resistance in this region might be the widespread use of agricultural pesticides. Unfortunately, no baseline data on KDT for susceptible reference strains of P. argentipes are available in South-East Asia to compare our findings with. The lethal time estimated for 50% mortality (LT50) of the P. argentipes colony population of Rajendra Memorial Research Institute of Medical Sciences, India was <15 min for deltamethrin 0.05% [10], almost half of the time predicted for knockdown in our study. Another study from India [30] reported similar knockdown times as observed in our study (KDT50) for wild-caught P. argentipes against different classes of insecticides. The fact that the highest KDT50 was observed for DDT 4% was not surprising, as resistance in the vector population for this insecticide was already well-established and acknowledged.

One of the major strengths of this study is the baseline information generated on the susceptibility status of the alternative classes of insecticides, carbamate, and organophosphate, which will be valuable for evidence-based decision making in insecticide resistance management of sand flies. In addition, the observed low intensity of resistance to pyrethroids suggests that these insecticides are not completely obsolete and supports the possibility of reversal to susceptibility in vector sand flies provided that alternative insecticides are used for a few years as resistance mitigation measures. Also, for the first time, we attempted PBO synergist-pyrethroid bioassays to explore possible underlying monooxygenase-based metabolic resistance mechanism. The presence of pyrethroid-DDT cross-resistance can be assessed as a proxy indicator of the presence of kdr mutations within the vgsc gene in nerve cells for the expression of resistance [11,30]. In the context of confirmed resistance against pyrethroids, if the VL control programme in Nepal plans to switch to alternative classes of insecticides for IRS, organophosphates could be a choice considering the monooxygenase-based metabolic cross-resistance between pyrethroid and carbamate classes of insecticides [37].

The procedures for conventional bioassays are sophisticated in terms of resources and skills, however, these are WHO approved and provide strong evidence on the focal status of phenotypic susceptibility of a particular vector species to a specific insecticide. One of the limitations of this study is that we had to implement the methodology and insecticide discriminating concentrations meant for testing the susceptibility of mosquitoes due to the absence of WHO standardized procedures and specific insecticide discriminating concentrations for testing the susceptibility of sand flies. Nonetheless, several authors have preceded us in the use of this adopted methodology to assess the susceptibility status of sand flies in the Indian subcontinent [10,12,30,36]. A recent study in Sri Lanka has quantified far lower lethal concentrations of insecticides to kill 100% (LC100) of P. argentipes as compared to discriminating concentrations used to monitor the resistant status of mosquitoes [38]. In that case, the resistance status observed in our study would likely be an underestimation. To overcome this limitation, WHO is undertaking a global multi-centre study for determining insecticide discriminating concentrations of various insecticides for sand flies and the results are expected in 2022 (RS Yadav, personal communication). Another limitation of the study was that we could not standardize the age and physiological stages as we conducted the bioassays with wild captured adult sand flies. The fact that the age distribution and physiological stages of the sand flies caught were representatives of the real-time wild vector population [16], this could be considered a strength of the study as well.

We cannot assume that the susceptibility status of the vector sand flies observed in the five selected VL endemic districts of Eastern and Central Nepal are representative for other endemic districts from the Western part or the hilly districts located in diverse geo-ecological zones and having unlike IRS history. Monitoring of the susceptibility status should therefore be carried out in diverse geo-ecological areas, to allow for area-specific recommendations. Non-endemic districts should not be excluded from this exercise, as sporadic VL outbreaks might warrant effective evidence-based vector control interventions as well.

Conclusion

After three decades of IRS with a single class of insecticides in Nepal, we confirmed the presence of pyrethroid resistance in P. argentipes sand flies. It is strongly recommended that the national vector-borne disease control programme takes this evidence into account when deciding on insecticides to be used for future IRS activities, and switches to an alternative class of insecticides with a different mode of action such as organophosphates or novel products prequalified by WHO. Also, this study endorses the need for regular and systematic monitoring of the susceptibility status of the vector of VL to inform evidence-based decision making for the use of insecticides in the fight against VL.

Supporting information

S1 Fig. Percentage composition of physiological stages of female P. argentipes used in the test (T) and the control (C) replicates.

Inside the doughnut chart: A, alive and D, dead sand flies after 24h recovery period. UF, unfed; BF, blood-fed and G, gravid P. argentipes.

https://doi.org/10.1371/journal.pntd.0010304.s001

(TIF)

S1 Data. Dataset of sand flies collected, insecticide tested and results of bioassays in detail.

https://doi.org/10.1371/journal.pntd.0010304.s002

(XLS)

S1 File. R-scripts for all statistical analysis and figures.

https://doi.org/10.1371/journal.pntd.0010304.s003

(DOC)

Acknowledgments

We would like to thank all the female community health volunteers from the respective study villages for guiding us through the village and facilitating communication with the community people. We are also grateful to the vector control officers from the respective district public health offices. We appreciate members of our field team; Manish Karn, Kailash Majhi, Jibachh Sah, and Durga Hasda for their sincere effort in the collection of adequate sand flies and to assist in the conduction of bioassays. The co-authors UK and RSY are staff members of the World Health Organization. They are responsible for the views expressed in this article, which do not necessarily represent the decisions, policies, or views of the WHO.

Citation: Roy L, Uranw S, Cloots K, Smekens T, Kiran U, Pyakurel UR, et al. (2022) Susceptibility status of the wild-caught Phlebotomus argentipes (Diptera: Psychodidae: Phlebotominae), the sand fly vector of visceral leishmaniasis, to different insecticides in Nepal. PLoS Negl Trop Dis 16(7): e0010304. https://doi.org/10.1371/journal.pntd.0010304

References

1. World Health Organization. Regional Office for South-East A. Regional strategic framework for elimination of kala azar from the South-East Asia Region (2005–2015). New Delhi: WHO Regional Office for South-East Asia; 2005.

2. Epidemiology and Disease Control Division. National Guideline on Kala-azar Elimination Program (Updated) 2019: Department of Health Services, Ministry of Health and Population, Government of Nepal, Teku, Kathmandu. 2019.

3. Boelaert M, Meheus F, Sanchez A, Singh SP, Vanlerberghe V, Picado A, et al. The poorest of the poor: a poverty appraisal of households affected by visceral leishmaniasis in Bihar, India. Trop Med Int Health. 2009;14(6):639–44. pmid:19392741

4. 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. pmid:30670453

5. Banjara MR, Das ML, Gurung CK, Singh VK, Joshi AB, Matlashewski G, et al. Integrating Case Detection of Visceral Leishmaniasis and Other Febrile Illness with Vector Control in the Post-Elimination Phase in Nepal. Am J Trop Med Hyg. 2019;100(1):108–14. pmid:30426921

6. Das ML, Roy L, Rijal S, Paudel IS, Picado A, Kroeger A, et al. Comparative study of kala-azar vector control measures in eastern Nepal. Acta Trop. 2010;113(2):162–6. pmid:19879851

7. Joshi AB, Das ML, Akhter S, Chowdhury R, Mondal D, Kumar V, et al. Chemical and environmental vector control as a contribution to the elimination of visceral leishmaniasis on the Indian subcontinent: cluster randomized controlled trials in Bangladesh, India and Nepal. BMC Med. 2009;7:54. pmid:19804620

8. Epidemiology and Disease Control Division. 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.

9. Joshi AB, Bhatta LR, Regmi S, Ashford RW. An assessment of the effectiveness of insecticide spray in the control of visceral leishmaniasis in Nepal. Journal of Nepal Health Research Council. 2003;1:1–6.

10. 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(10):e859. pmid:21049013

11. 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(4):e0005504. pmid:28414744

12. 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. Parasit Vectors. 2018;11(1):336. pmid:29866195

13. World Health Organization. Management of insecticide resistance in vectors of public health importance. Report of the ninth meeting of the global collaboration for development pesticides for public health. 9–10 September 2014. 2014.

14. DHM. Observed climate trend analysis in the districts and physiographic regions of Nepal (1971–2014). Department of Hydrology and Meteorology, Kathmandu; 2017.

15. Karki R, Talchabhadel R, Aalto J, Baidya S. New climatic classification of Nepal. Theoretical and Applied Climatology. 2015.

16. World Health Organization. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes (2nd edition). Global Malaria Program. Geneva: WHO; 2016. p. 56.

17. Kalra NL, Bang YM. Manuals on Entomology in Visceral Leishmaniasis: World Health Organization, SEARO, New Delhi; 1988.

18. Lewis DJ. The Phlebotominae sandflies (Diptera: Psychodidae) of the oriental region. Bulletin of British Museum (Natural History) of Entomology. 1978;37:217–343.

19. Lewis DJ. A taxonomic review of the genus Phlebotomus (Diptera: Psychodidae). Bulletin of British Museum of Entomology (Natural History). 1982;45:121–209.

20. Abbott WS. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18:265–7.

21. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.

22. Goodrich B, Gabry J, Ali I, Brilleman S. rstanarm: Bayesian applied regression modeling via Stan. R package version 2.21.1. 2020.

23. Brilleman SL, Crowther MJ, Moreno-Betancur M, Buros Novik J, Wolfe R. Joint longitudinal and time-to-event models via Stan. StanCon 2018. 10–12 Jan 2018. Pacific Grove, CA, USA. 2018.

24. Finney DJ. Probit Analysis: Cambridge University Press, Cambridge, England, ISBN: 052108041X; 1971.

25. Wheeler MW, Park RM, Bailer AJ. Comparing median lethal concentration values using confidence interval overlap or ratio tests. Environ Toxicol Chem. 2006;25(5):1441–4. pmid:16704080

26. Robertson JL, Savin NE, Russel RM, Preisler HK. Bioassays with Arthropods. CRC press. ISBN: 9780849323317; 2007.

27. Hlina BL, Birceanu O, Robinson CS, Dhiyebi H, Wilkie MP. Seasonal Variation in the Sensitivity of Invasive Sea Lampreys to the Lampricide TFM: Importance of Energy Reserves and Temperatures. N Am J Fish Manage. 2019.

28. Mittal PK, Wijeyaratne P, Pandey S. Status of insecticide resistance of Malaria, Kala-azar and Japanese Encephalitis vectors in Bangladesh, Bhutan, India and Nepal (BBIN). Office of Health, Infectious Diseases and Nutrition, Bureau for Global Health, US agency for International Development, Washington, DC 20523; 2004. Report No.: 129.

29. Anonymous. The annual internal assessment of malaria and kala-azar control activities 2000. Epidemiology and Disease Control Division, Department of Health Services, Ministry of Health, His Majesty’s Government of Nepal and Vector Borne Disease Research and Trainining Centre, November 2000, Kathmandu, Nepal; 2000.

30. Sardar AA, Saha P, Chatterjee M, Bera DK, Biswas P, Maji D, et al. Insecticide susceptibility status of Phlebotomus argentipes and polymorphisms in voltage-gated sodium channel (vgsc) gene in Kala-azar endemic areas of West Bengal, India. Acta Trop. 2018;185:285–93. pmid:29890155

31. Singh RK, Mittal PK, Dhiman RC. Insecticide susceptibility status of Phlebotomus argentipes, a vector of visceral leishmaniasis in different foci in three states of India. J Vector Borne Dis. 2012;49(4):254–7. pmid:23428526

32. Kumar V, Shankar L, Kesari S, Bhunia GS, Dinesh DS, Mandal R, et al. Insecticide susceptibility of Phlebotomus argentipes & assessment of vector control in two districts of West Bengal, India. Indian J Med Res. 2015;142(2):211–5. pmid:26354219

33. Nandi J, Chaudhuri I, Das C, Mukhopadhyay A. Susceptibility Status of Phlebotomus argentipes, to Synthetic Pyrethroid Insecticides in Kala Azar Endemic Parts of Bihar, India in Relation to Elimination of Visceral Leishmaniasis. Indian Journal of Public Health Research & Development. 2016;7:36. https://doi.org/10.5958/0976-5506.2016.00185.6

34. Dhiman RC, Raghavendra K, Kumar V, Kesari S, Kishore K. Susceptibility status of Phlebotomus argentipes to insecticides in districts Vaishali and Patna (Bihar). J Commun Dis. 2003;35(1):49–51. pmid:15239308

35. Singh R, Kumar P. Susceptibility of the sandfly Phlebotomus argentipes Annandale and Brunetti (Diptera: Psychodidae) to insecticides in endemic areas of visceral leishmaniasis in Bihar, India. Jpn J Infect Dis. 2015;68(1):33–7. pmid:25420660

36. Dinesh DS, Hassan F, Kumar V, Kesari S, Topno RK, Yadav RS. Insecticide susceptibility of Phlebotomus argentipes sandflies, vectors of visceral leishmaniasis in India. Trop Med Int Health. 2021;26(7):823–8. pmid:33733549

37. World Health Organization. Global plan for insecticide resistance management in malaria vectors (GPIRM). World Health Organization, 20 Avenue Appia, 1211 Geneva 27, Switzerland. 2012.

38. Pathirage DRK, Karunaratne S, Senanayake SC, Karunaweera ND. Insecticide susceptibility of the sand fly leishmaniasis vector Phlebotomus argentipes in Sri Lanka. Parasit Vectors. 2020;13(1):246. pmid:32404115