Current Perspectives on Plague Vector Control in Madagascar: Susceptibility Status of Xenopsylla cheopis to 12 Insecticides.

Plague is a rodent disease transmissible to humans by infected flea bites, and Madagascar is one of the countries with the highest plague incidence in the world. This study reports the susceptibility of the main plague vector Xenopsylla cheopis to 12 different insecticides belonging to 4 insecticide families (carbamates, organophosphates, pyrethroids and organochlorines). Eight populations from different geographical regions of Madagascar previously resistant to deltamethrin were tested with a World Health Organization standard bioassay. Insecticide susceptibility varied amongst populations, but all of them were resistant to six insecticides belonging to pyrethroid and carbamate insecticides (alphacypermethrin, lambdacyhalothrin, etofenprox, deltamethrin, bendiocarb and propoxur). Only one insecticide (dieldrin) was an efficient pulicide for all flea populations. Cross resistances were suspected. This study proposes at least three alternative insecticides (malathion, fenitrothion and cyfluthrin) to replace deltamethrin during plague epidemic responses, but the most efficient insecticide may be different for each population studied. We highlight the importance of continuous insecticide susceptibility surveillance in the areas of high plague risk in Madagascar.

Introduction

Arthropod-borne diseases are a major concern worldwide. Every year more than 1 billion cases and over 1 million deaths from vector-borne diseases are estimated [1]. Most of these vectors are bloodsucking arthropods (e.g., mosquitoes, flies, fleas, ticks, lice) living in direct contact with humans or harbored by livestock or commensal animals [2-5]Arthropod-borne diseases such as malaria, chikungunya, dengue fever, Lyme disease, West Nile fever, Rift Valley fever, and plague erupt and cause substantial mortality in humans and livestock [1,5-8]. The risk posed by these diseases can be significantly reduced by the use of insecticides during a public health emergency; insecticide-based intervention can prevent an outbreak or can limit the expansion of the disease [3,5].

In the 1940s, the use of synthetic insecticides led to great improvement in the battle against disease vectors [3,9]. Consequently, the intensive use of insecticides caused selection pressure on insect populations, which developed mechanisms to survive insecticide treatments. All classes of insecticides are currently involved, and the list of pests associated with agriculture and health has been continually increasing [10]. The lack of new insecticidal compounds and the misuse or overuse of insecticides have been identified as reasons of the emergence of insecticide resistance in pests [11]. Hence, an efficient vector control policy must take in account the possibility of insecticide resistance, which can lead to a failure of the control strategy [12]. The early detection and monitoring of insecticide resistance in a vector population may positively impact intervention strategies. Until now, the main defense against resistance is close surveillance of the susceptibility of vector populations [4].

Plague, a rodent disease transmitted to human by infected flea bites, remains an important health problem in Madagascar [13]. The flea, Xenopsylla cheopis is the main plague vector, parasitizing black rats Rattus rattus that live in urban and rural housing [14]. According to the World Health Organization (WHO), the most rapid and effective method for controlling fleas is to apply an appropriate insecticide formulated as a dust or low-volume spray [15,16]. Insecticide dusting in households is the strategy adopted by the National Plague Control Program in Madagascar to control vectors and to limit the expansion of plague epidemics [17,18]. In 1947 the use of DDT (dichlorodiphenyltrichloroethane) insecticide to control fleas in Madagascar gave new hope for combatting plague. Since 1965, resistance of X. cheopis to DDT has been developing in Madagascar [19]. Later, X. cheopis populations were reported to be resistant to different families of insecticide from the early 1980s to 2000 [20-26]. More recently, amongst 32 populations of X. cheopis, only two populations were susceptible to deltamethrin, which is currently the preferred insecticide in Madagascar for flea control [27]. Hence, it is crucial to find insecticide alternatives to deltamethrin. Here we report the results of 12 different insecticide bioassays performed on 8 populations of X. cheopis previously found to be resistant to deltamethrin.

Materials and Methods

Flea populations

A previous study reported the susceptibility of 32 populations of fleas to deltamethrin [27]. We chose to study eight populations from different geographical regions of Madagascar (Fig 1). Chosen populations were resistant fleas with mortality rates when exposed to deltamethrin of 2.5% to 65% [27]Flea populations (X. cheopis) were collected from the field and reared in an insectarium [27]. Briefly, rodents were trapped alive, fleas were combed into a large container, and fleas were reared in insectarium at 22–27°C and 75–80% relative humidity until having the sufficient number to perform bioassays. [27]. Fleas used in bioassays were subsequent generation of those collected in field.

Fig 1

Map showing location of the eight sites where fleas were collected.

Source: Institut Pasteur de Madagascar, OCHA.

WHO bioassay

Bioassays were conducted on fleas populations according to the WHO protocol [28]. Ten adult fleas per tube were exposed to insecticide-impregnated paper (1.5 x 6 cm; Vector Control Research Unit, Penang, Malaysia) for specified times and at predetermined insecticide concentrations (Table 1). Each test was replicated at least four times for a total of 40 fleas per insecticide and per population. Negative controls were performed with paper only impregnated with the carrier of each insecticide family. LT50 (the time by which 50% of fleas were knocked down) were estimated for each insecticide during the diagnostic time. At the end of the exposure time and for all bioassays, the impregnated papers were removed and replaced by non-impregnated papers. Final mortality was recorded 24 hours after the beginning of exposure time. Susceptibility status was established according to the WHO guidelines for insecticide susceptibility test. Mortality rates of 98 to 100% indicated susceptibility, 80 to 98% tolerance or suspected resistance, and less than 80% resistance [29]. The test was not validated, and the data not included, if the negative control mortality rate was over 20%. The mortality rate was corrected with the Abbott formula [30] when control values were between 5% and 20%.

Table 1

Insecticides used in the bioassays with their concentration and the diagnostic exposure times.

Insecticide	Family	Concentration (%)	Diagnostic time (hours)
Alphacypermethrin	Pyrethroid	0.025	8
Lambdacyhalothrin	Pyrethroid	0.05	8
Etofenprox	Pyrethroid	0.5	8
Permethrin	Pyrethroid	0.75	8
Cyfluthrin	Pyrethroid	0.15	8
Deltamethrin	Pyrethroid	0.05	8
Bendiocarb	Carbamates	0.1	5
Propoxur	Carbamates	0.1	5
Malathion	Organophosohate	5	5
Fenitrothion	Organophosphate	1	5
DDT	Organochlorine	4	6
Dieldrin	Organochlorine	4	6

Statistical analysis

Analysis of Variance (ANOVA) and Tukey’s b test were used to compare mortality rates. Mean LT50 and the standard errors for each flea population and for each insecticide were estimated with a binomial generalized linear model (glm) analysis. This glm including a probit function is a fitted model giving a prediction and a standard error at each response probability (p.model function with the package MASS). High mortality may not occur with some insecticides for some populations and therefore the LT50 would not be estimated (NE) Correlations between the mortality rates were calculated with Pearson tests (packages: corplot, Hmisc and ggplot2 to generate figures). Statistical analyses were done with R software (RStudio) [31].

Results

Mortality rate

The mortality rate was different among insecticides and populations: mean mortality was significantly different according to populations (F value = 195.34, p < 0.0001) and insecticide (F value = 36.22, p < 0.0001). A strong correlation between insecticides and populations (F value = 9.10, p < 0.0001) was observed. Nonetheless, all populations were at least somewhat resistant to the six insecticides alphacypermethrin, lambdacyhalothrin, etofenprox, deltamethrin, bendiocarb and propoxur, with mortality rate ranging from 0 to 79% (Fig 2). Dieldrin was the only insecticide with 100% mortality rate for all flea populations (Fig 2). The resistance to DDT was substantial for most populations, with mortality rates varying between 5 and 26.4%, with the exception of one tolerant population, Andranomanalina, which had 90% mortality. Apart from the dieldrin, the highest mortality rates were observed for malathion, fenitrothion, cyfluthrin and permethrin (Fig 2). For these four insecticides, the susceptibility profiles were very different for each population (Fig 3). Almost the same resistance profile was observed for the organophoshates: two populations (Ambohimiandra and Ambohipananina) were susceptible with 100% mortality rate and fleas from Tsararano Ambony and Amparaky were both resistant to malathion and fenitrothion. Mortality induced by cyfluthrin ranged from 67.5 to 100% with one susceptible population (Ambohipananina) and four populations were tolerant. The mortality rate of the eight populations with permethrin varied between 50 and 95%, with two resistant populations.

Fig 2

Box-and-whisker plot of mortality rate after 24 hours for each in secticide for all populations study sites.

Black Diamond-shaped points inside the boxes are mean values. Horizontal bars in boxes are the 50th percentiles (medians), and the bottom and the top of the box represent the 25th and the 75th percentiles, respectively. The two limits of vertical lines above and at the bottom of the box are the wiskers and represent the maximum and the minimum values of the data. Points outside the limit of vertical line are “outlier”, which are values outside 95% the confidence interval.

Fig 3

Mean of mortality rate (24 hours) per insecticide for each population.

X-axis indicates populations and Y-axis indicates mortality rate in percent. Error bars represents standard errors. Diagrams color codes: in red are resistant populations, in yellow, tolerant and in green, susceptible according to WHO. Letters code (a, b, and c) above and on the side of each bar plot indicate significant difference between the mortality for each population according to the Tukey’s b test.

Lethal Time 50 and mortality rate

The curve profile, obtained during exposure time for each insecticide and for each station (Fig 4), and values of LT50 (S1 Table) were in concordance with the results obtained with the average mortality observed after 24 hours. Highly resistant population to insecticide had LT50 values longer than durations of exposure time. For Etofenprox, six tested populations had estimated LT50 > 500 minutes whereas exposure time was 480 minutes. These six populations had mortality below 30% after 24 hours. For DDT, no tested population reached LT50 until the exposure time (LT50 > 360 minutes), except for the tolerant population of Andranomanalina with a LT50 = 142 ± 8.80 minutes (Fig 4). Flea populations susceptible to cyfluthrin (Ambohipananina) had a LT50 equal to 21 ± 1.88 minutes and the tolerant ones had LT50 between 36 ± 4.58 and 128 ± 10.26 minutes. Even though the LT50 value of the population most resistant to cyfluthrin (Tsararano Ambony) was seven times higher than the value for the most susceptible population, the possible emergence of resistant individuals in tolerant populations could be suspected.

Fig 4

Mean mortality rates of flea populations per insecticides through the exposure time.

Similarly for malathion, susceptible population had an LT50 < 104 ± 5.14 minutes. LT50 was reached after 152 ± 6.81 to 285 ± 6.97 minutes for tolerant populations. Resistant populations’ LT50 values (Amparaky and Tsararano Ambony) exceeded the exposure time (300 minutes): LT50 were 370 ± 31.03 minutes and 464 ± 128.71 minutes, corresponding to 57.7% and 10% mean mortality after 24 hours, respectively.

Fig 4 illustrates the heterogeneity of response to insecticides amongst different flea populations. For example, in dieldrin trials, although 100% mortality rate was observed after 24 hours for every population, one population (Ambohimiandra) did not reach its LT50 value until the end of the exposure time (360 minutes). The LT50 value (425 ± 38.51 minutes), was 4 times higher than the minimal value obtained for dieldrin (99±6.84 minutes).

With propoxur, 100% of exposed fleas were knocked down before the exposure time was elapsed in two populations (Ambohimiandra and Ambohipnanina). These two populations had the shortest values of LT50 and lowest mortality rate values to propoxur after 24 hours (highly resistant populations).

Mortality rate correlations between insecticides

Positive correlations were observed between deltamethrin, etofenprox, cyfluthrin, malathion, fenitrothion and propoxur (Fig 5), suggesting possible insecticide cross-resistance mechanisms in fleas. A strong negative correlation was observed between permethrin and etofenprox (r = -0.74, p<0.05). Significant correlations (p<0.05) were observed between fenitrothion and propoxur (r = 0.76, p = 0.03), propoxur and cyfluthrin (r = 0.77, p = 0.02), malathion and propoxur (r = 0.82, p = 0.01), and fenitrothion and cyfluthrin (r = 0.82, p = 0.01).

Fig 5

Schematic illustration of correlation between insecticide mortality rates (mean of 24h mortality per insecticide for all stations).

Positive correlations are in blue and negative correlations are in red. The intensity of colors and the size of circles are proportional to the correlation coefficient. In the right, the color legend shows the correlation coefficient values with corresponding colors. Numbers followed by star are significant p values (p<0.05). Positive correlations between insecticides are represented surrounded with red rectangle.

Discussion

X. cheopis resistance and insecticide uses in Madagascar

X. cheopis populations tested in this study were found highly resistant to DDT. Seven of eight populations showed final mortality rate less than 30% to DDT. These results reflected past observations of DDT resistance amongst X. cheopis populations from Madagascar, even though this product has not been used for many decades. The main argument raised to explain X.cheopis resistance to DDT worldwide was the extensive use of this insecticide in plague and malaria vector control [32-37]. In Madagascar, DDT was widely used against rat fleas since the 1940s [19]. DDT and pyrethroids were used in Indoor Residual Spraying (IRS) and long lasting insecticide impregnated nets against malaria vectors [38]. Furthermore, malaria vector treatment could have effect on flea vector resistance; in fact it was demonstrated elsewhere that insecticides used in IRS programs reduced flea loads on indoor rodents. [39]. However, in areas where malaria and plague are endemic, IRS treatment could have the potential to put selective pressure on fleas to develop resistance [40].

Dieldrin, an organochlorine insecticide, also saw widespread use in countries where plague occurred. Dieldrin was used in Madagascar during the period of DDT use, and X. cheopis was already described as resistant to dieldrin [21,26]. Insecticide susceptibility tests done in India showed that fleas resistant to DDT often were resistant to Dieldrin and other cyclodien insecticides [25, 29]. Yet, in our study, X. cheopis populations were all susceptible to this compound.

Pyrethrum was shown to have lethal effects on rat fleas before synthetic pyrethroids were used [41]. X. cheopis resistance to pyrethroid compounds (deltamethrin 0.025% and cyfluthrine 0.15%) was previously described in Madagascar [23,24]. Fleas from the Central Highlands of Madagascar were resistant to low concentrations of deltamenthrin (0.025%). Recently, 94% of studied populations were not sensitive to higher concentrations of deltamethrin (0.05%) [23,24,27]. The use of deltamethrin in plague control since the 1990s likely led to the development of resistant flea populations.

We present the first data illustrating resistance of X. cheopis populations to alphacypermethrin, lambdacyhalothrin and etofenprox, which were never used in mass vector control. This may suggest the involvement of cross-resistance mechanisms between these insecticides and those that were extensively used.

Organophoshates were also described as inducing resistance in rat flea populations. In India, X. cheopis was indicated as resistant to malathion and fenitrothion, although these compounds were never used in the study areas [36]. It was suggested that flea resistance to these compounds was associated with resistance to DDT [42]. Even if resistance to organophosphates was already described in some areas of Madagascar, the majority of populations studied presently showed less resistance to these compounds [21,43]. On the other hand, X. cheopis populations were previously described as susceptible to carbamate insecticides [24]; however, our study demonstrated a high resistance to propoxur and bendiocarb.

Hypothesis on resistance mechanisms

Our results suggest resistance to all insecticides except dieldrin, which produced 100% mortality for all population. However, the LT50 values observed in one population (Ambohimiandra) suggest a progressive development of resistance to this compound. But dieldrin was banned in most of country worldwide because of its high toxicity in mammals and its bioaccumulation in the environment [44,45]. The use of dieldrin was suspended in Madagascar since 1993 [46]. However, other insecticide families having the same mode of action as dieldrin (antagonist of GABBA receptor) such as fiproles could be promising [10,47].

Six insecticides (alphacypermethrin, lambdacyhalothrin, etofenprox, deltamethrin, bendiocarb and propoxur) were relatively ineffective for flea control in all populations. Nonetheless, resistance level to the insecticides (permethrin, cyfluthrin malathion and fenitrothion) was very different among populations, suggesting different selection pressures. Hence, in this study, according to WHO thresholds, some insecticides were still efficient in some localities; thus, insecticides that induce resistance according WHO thresholds still may exhibit high performance in the field [12].

The different responses of populations to each insecticide reflect also the mode of action of insecticide molecules and the mechanism developed by insects to overcome toxic effects. Pyrethroids and DDT belong to a group of neurotoxic chemicals and share a similar mode of action distinct from other classes of insecticide. The studies on kdr mutation demonstrated the same mode of action of pyrethroids and DDT, which is the reduced target-site sensitivity of sodium channel [48]. Thus, the mechanism of resistance may not be specific to a particular insecticide family or group but the molecule structure of each insecticide can play important role.

For instance, the negative correlation between permethrin and etofenprox may involve the different effect induced by a Type I pyrethroid (permethrin) and a pseudo-pyrethroid (nonester pyrethroid) [49]. In addition, different levels of pyrethroid resistances were observed amongst populations. All studied populations were resistant to etofenprox alphacypermethrin, lambdacyhalothrin and deltamethrin; yet cyfluthrin and permethrin were effective in some localities. In a study of cross resistance amongst pyrethroids, cross resistance between 19 pyrethroid insecticides was assessed in bollworm moth, Helicoverpa armigera [50]. Cross resistance between pyrethroids seemed due to enhanced oxidative metabolism induced by pyrethroid with the same structure. The modification or replacement of any compound (aromatic compound) in the molecule structure could modify the susceptibility of the population [50].

Moreover, DDT and dieldrin belong to the oragnochlorine family, but their structures are very different, conferring different mode of action. The first attempt to elucidate cross resistance between chlorinated insecticides in X. cheopis was performed in 1974 [27]; a DDT-selected population was found to be resistant also to insecticides structurally related to DDT, and exhibited variable resistance to cyclodiene insecticides (such as dieldrin, endrin). But biochemical assays did not show significant difference between susceptible and resistant population. [32]. Other mechanisms like Rdl mutation can confer resistance to cyclodiene like dieldrin [47].

Furthermore, the correlations between deltamethrin, etofenprox, propoxur, and between cyfluthrin, malathion and fenitrothion, may be explained by the same mechanism of resistance developed by X. cheopis for these insecticides. The absence of references on this topic in Siphonaptera illustrates a need for further research on insecticide resistance mechanisms in fleas.

Perspectives on improving fight against fleas vector of plague

Efforts can be undertaken before each epidemic season in order to control the proliferation of vectors and host, such as live rat mass trapping, promotion of rat proofing in houses and environment sanitation [15]. One available method could be the use of insecticide bait box, combining insecticide and delayed toxicity rodenticide [51-55]. The objective is to kill fleas on rodents before the action of the rodenticide. Instead of insecticide dusting in household, the quantity of insecticide is reduced because the insecticide bait box is more focused on fleas harbored by rodent with the host acting as a vehicle for the insecticide, carrying it to its nest. In the same line of thought, the feasibility of “insecticide delivery tubes” in reducing flea loads was studied on commensal rodents, capitalizing on the knowledge of their behavior [56]. Similarly, using rodent bait containing systemic insecticide could be a new avenue for combating or at least, reducing fleas load on rodents in plague endemic area during inter-epidemic season [57]. Besides, novel approaches to fighting vector limiting the use of chemical insecticides should be explored in order to avoid insecticide resistance [58,59]. For instance, research must be undertaken in the way to better understand the interaction between the vector, the pathogen and the insect microbiome. The strategy is based on introduction of microorganism which may affect the insect lifespan or the ability to transmit the pathogen [59-61]. Furthermore, research on bio pesticide is already ongoing with the use of entomopathogen fungi to reduce the survival of flea larvae [62].

Recommendations for public health concern

The main finding in this study is that X. cheopis populations developed resistance mechanisms to the insecticide families most used in vector control. The description of phenotypic resistance to insecticides is important to help practitioners choose the most efficient strategy in pest management. Hence, in a public health context, insecticide susceptibility status of fleas in each plague risk area may be monitored periodically to conduct more focused and adapted flea control. However information available on the mechanism of resistance and cross-resistance about X. cheopis is scarce or nonexistent. Research must be done to understand the mechanisms conferring resistance to insecticides in plague vectors such X. cheopis.

Table of the values of LT50 in minutes and standard error values for each insecticide and station.

(XLSX)

Click here for additional data file.