Skip to main content

Malaria infection in mosquitoes decreases the personal protection offered by permethrin-treated bednets

Abstract

Background

Insecticides targeting adult mosquitoes are the main way of controlling malaria. They work not only by killing mosquitoes, but also by repelling and irritating them. Indeed their repellent action gives valuable personal protection against biting mosquitoes. In the context of malaria control this personal protection is especially relevant when mosquitoes are infectious, whereas to protect the community we would prefer that the mosquitoes that are not yet infectious are killed (so, not repelled) by the insecticide. As the infectious stage of malaria parasites increases the motivation of mosquitoes to bite, we predicted that it would also change their behavioural response to insecticides.

Results

With two systems, a laboratory isolate of the rodent malaria Plasmodium berghei infecting Anopheles gambiae and several isolates of P. falciparum obtained from schoolchildren in Tanzania that infected Anopheles arabiensis, we found that mosquitoes harbouring the infectious stage (the sporozoites) of the parasite were less repelled by permethrin-treated nets than uninfected ones.

Conclusions

Our results suggest that, at least in the laboratory, malaria infection decreases the personal protection offered by insecticide-treated nets at the stage where the personal protection is most valuable. Further studies must investigate whether these results hold true in the field and whether the less effective personal protection can be balanced by increased community protection.

Background

Indoor residual spraying (IRS) and insecticide-treated bednets (ITNs) are among the most commonly used and cost-effective control measures against malaria [1]. Indeed, these measures significantly reduced malaria transmission, for example, in sub-Saharan Africa [2]. Both protect against malaria by repelling mosquitoes and by killing mosquitoes that rest on the treated surface after having taken a blood meal [3,4,5,6,7]. ITNs give additional protection by creating a physical barrier that reduces human-mosquito contacts.

The repellent action of IRS and ITNs offers effective personal protection [8,9,10,11] by reducing the number of bites (infectious bites in particular) on users; their insecticidal action offers protection to the community, first, by killing mosquitoes before the malaria parasite has completed its development and become transmissible and secondly, by killing infectious mosquitoes [3]. A striking example of community protection is that a 75% coverage of ITNs in a rural community in Tanzania decreased the entomological inoculation rate of the people that did not use an ITN by a factor of 18 [12]. To some extent, the two levels of protection oppose each other: strongly repellent nets, which give strong personal protection, let only few mosquitoes contact the insecticide on the nets, thus killing few mosquitoes and giving weak community protection [3, 13,14,15].

To weaken this trade-off, we would want an ITN to repel mainly infectious mosquitoes, thus protecting users against being infected, but to repel other mosquitoes only weakly, thus maintaining community protection by lowering the number of mosquitoes that become infectious. Our prediction, however, was the opposite: as malaria changes the host-seeking behaviour of mosquitoes, its infectious stages would decrease repellency to give less rather than more personal protection against infection.

Indeed, infection with the infectious stage (sporozoites) increases the number of human-mosquito contacts in several ways. Sporozoite-infected mosquitoes have been found to be more attracted to human odours than uninfected ones [16, 17], though not all experiments show the same pattern [18]. Sporozoites also increase the mosquitoā€™s persistence and biting rate [19,20,21,22] and increase the likelihood that the mosquitoes bite a second time to top up their blood meal [23]. These changes of behaviour probably evolved in order to increase the transmission of the parasite [24]. The general pattern of the behavioural changes - that sporozoite-infected mosquitoes are more motivated to bite than uninfected ones - may well have as a consequence that sporozoite-infected mosquitoes are less repelled by insecticides and therefore more likely to be killed by them.

The aim of this study was to test whether permethrin-treated bednets repel mosquitoes harbouring sporozoites less than uninfected mosquitoes (we chose to work with permethrin because it is one of the most used pyrethroids for insecticide-treated nets [25]). To do so, we performed two experiments in the laboratory that differed in the malaria parasite, the mosquito, and the way we measured the repellency of permethrin. In a first experiment with the rodent malaria Plasmodium berghei infecting Anopheles gambiae, mosquitoes chose between feeding through a permethrin-treated net or through an untreated net. In a second experiment with the human malaria P. falciparum infecting An. arabiensis, mosquitoes had to pass through a treated or untreated net if they wanted to reach the blood source. Note that, while for simplicity we use the term repellency in both experiments, we measure different behaviours. In the first, we measure preference of untreated nets, whereas in the second we measure the motivation to bite despite the presence of the insecticide. Note also that in both experiments, repellency refers to a combination of spatial repellency and contact irritancy [7], for our experimental design does not let us distinguish between the two.

Methods

Plasmodium berghei and Anopheles gambiae

Mosquito and parasite

We used the Kisumu colony of An. gambiae originating from western Kenya [26] and the mouse malaria P. berghei ANKA, modified to express the green fluorescent protein (GFP). Infected mice harbouring gametocytes were obtained from Volker Heussler (Institute of Cell Biology, University of Bern) and used to feed and infect the mosquitoes.

Mosquito rearing and infection

We haphazardly selected larvae from the colony and reared them individually at 26 Ā± 1 Ā°C in 12-well plates with a standard level of fish food (Tetramine) [26]: day of hatching of eggs, 0.04 mg per larva; 1 day after hatching, 0.06 mg; day 2, 0.08 mg; day 3, 0.16 mg; day 4, 0.32 mg; day 5 and later, 0.6 mg. Upon pupation, each mosquito was moved to a 120 ml plastic cup.

Three to four days after emergence, we moved the mosquitoes to an insectary maintained at 19 Ā± 1 Ā°C and 70 Ā± 5 % of humidity for an efficient development of P. berghei in mosquitoes [27]. We haphazardly let females feed for 15 min on an infected or an uninfected mouse that had been anaesthetized. To replicate the infection-treatment, we used two infected and two uninfected mice, and let 100 mosquitoes feed on each one. Twenty-four hours after feeding, engorged mosquitoes were placed into 120 ml individual plastic cups; unfed mosquitoes were discarded. We gave the mosquitoes access to a cotton ball soaked with a 10% sugar solution up to 24 h before the behavioural test.

Measuring repellency

We defined repellency of a permethrin-treated net as the preference of mosquitoes for an untreated net when given the choice of feeding on a blood source through an untreated net or a permethrin-treated one. We measured the preference with a two-way choice apparatus (Fig 1), where mosquitoes could choose to fly towards and attempt to bite an arm held onto a cage containing either a 30 Ɨ 30 cm piece of net with permethrin (270 mg/m2) or netting without permethrin (Trek 1 or Trek, Katadyn France, Grenoble, France). A ventilator guided the odours (at a speed of about 20 cm/s) from the two cages into a central cage, into which the mosquitoes had been placed. The side of the insecticide was alternated among tests so that any side preference was avoided. Moreover, the potential bias for one of the sides of our setup was controlled by testing shortly before our experiment that there was no preference for either side when uninfected mosquitoes were given the choice of two untreated net. In eight replicates (each with 10 uninfected mosquitoes aged 20 to 22 days) a repeated G-test of goodness-of-fit [28] showed strong support for no preference [total-G = 1.48, df = 8, P = 0.99, pooled-G = 0.385, df = 1, P = 0.53, with no heterogeneity among the replicates (G = 1.1, df = 7, P = 0.99); Additional file 1: Table S3].

Fig. 1
figure 1

Two-way choice apparatus to test the effects of malaria infection on insecticide repulsion in the experiment with P. berghei infecting An. gambiae. The central cage (25 Ɨ 25 Ɨ 25 cm) was connected to two test cages (20 Ɨ 20 Ɨ 20 cm) with 15 cm length tubes (diameter 10 cm). The two small cages had a hole covered by the permethrin-treated net or the untreated net on the top. Mosquitoes were placed in the 15 Ɨ 10 Ɨ 5 cm acclimatization chamber for 2 min before their behaviour was tested. The arrows represent the direction of the air flow

Twenty to twenty-two days after the blood meal we tested 90 uninfected mosquitoes and 80 infected ones in groups of 10. For each replicate, we placed the 10 mosquitoes into the central acclimation chamber and KT put his arms as attractants onto the nets of the two small cages. After letting the mosquitoes acclimatize for 2 min we opened the chamber and gave the mosquitoes 8 min to make a choice. We recorded the number of mosquitoes in each cage and put the mosquitoes into a freezer (-20 Ā°C) until further analysis.

Malaria infection

DNA was extracted with DNAzol (MRC inc., Cinncinati, OH, USA) according to the manufacturerā€™s instructions modified by Rider et al. [29]. The head and thorax of mosquitoes were smashed with a pellet in 200 Ī¼l of DNAzol. The solution was incubated for 20 min at 55 Ā°C and then centrifuged at 20,000Ɨ g for 10 min. One hundred and seventy microliters of the supernatant were transferred to a tube together with 3 Ī¼l of PolyAcryl Carrier (MRC inc.) and 100 Ī¼l of cold (-20 Ā°C) 100% ethanol. After centrifuging at 15,000Ɨ g for 8 min we discarded the ethanol. We then washed the DNA by adding 0.8 ml of cold 75% ethanol and centrifuging the tube for 5 min at 15,000Ɨ g. Ethanol was again discarded and the pellet of DNA was dried in a Speedvac for 15 min at 45 Ā°C. Finally, the DNA was eluted in 50 Ī¼l of autoclaved, deionised water.

Sporozoite infection was confirmed by PCR using a modified protocol established by Rider et al. [29]. PCR amplification was done with the T3000 thermocycler [Analytik Jena AG (formerly Biometra), Jena, Germany] with 3 Ī¼l of DNA and 1 Ī¼l of a 10 Ī¼M solution of each primer, forward (5'-ACG ATG ATA TAG ATC AAA T-3') and reverse (5'-TAC CTA AGC TTC TTG CGT A-3'), which amplify a 111 bp sequence of the merozoite surface protein-1 (MSP-1) gene of Plasmodium berghei NK65 and ANKA. DNA was amplified during 40 cycles with denaturation at 95 Ā°C, annealing at 54 Ā°C and extension at 72 Ā°C, each for 45 s. Infection was confirmed from the correct band of DNA after electrophoresis on a 1.8% agarose gel stained with the RedSafeTM Nucleic Acid stain (iNtRON Biotechnology, Seongnam-si, South Korea).

Mosquito body size

Wing length was taken as a measure of mosquito body size. Wings were placed onto a glass slide and scanned. We used the software ImageJ to measure the distance from the distal end of the alula to the tip of the wing (vein R3) without the fringes [30, 31].

Statistical analysis

We analysed the proportion of mosquitoes that responded with a generalized linear mixed model (GLMM) with a binomial distribution where the infection status was set as a nominal factor and the replicate as a random factor (Additional file 1: Table S2). Using only the mosquitoes that responded, we analysed the preference with a GLMM, where the choice of mosquitoes was set as the response variable, the infection status as a nominal factor, the wing length and the day of test (i.e. the mosquitoā€™s age) as covariates, and the replicate as a random factor (Additional file 1: Table S1). We excluded exposed but not infected mosquitoes from the analyses, as there were too few (n = 7) that responded during the test to enable any meaningful statistical analysis or interpretation. The significance of the effects was determined with the ANOVA function (type II) from the ā€œcarā€ package [32]. The analysis was performed with the software R [33] and the Rstudio interface [34]. For the mixed effect models, we used the package ā€œlme4ā€ [35].

Plasmodium falciparum and Anopheles arabiensis

Study site and patients

The study was conducted from the beginning of May to the end of July 2016 near Bagomoyo in the coastal region of Tanzania, where P. falciparum is the main parasite causing malaria. All laboratory-based work was done at the Kingani insectary site of the Bagamoyo Research and Training Unit (a branch of the Ifakara Health Institute). Malaria parasites were obtained from 6 to 14 year-old children from villages around Bagomoyo. A first blood screening was done at the schools of the villages. We used malarial rapid diagnosis test (SD BIOLINE Malaria Ag P.f/Pan, Standard Diagnostic, Gyeonggi-do, South Korea) to determine if children were infected or not. Children carrying both P. falciparum and another Plasmodium species were treated with antimalarial drugs (CoartemĀ©, Novartis, Basel, Switzerland) according to the national guidelines. Blood smears of children infected only with P. falciparum were coloured with 10% Giemsa and checked under a microscope (oil film 100Ɨ magnification) for the presence of gametocytes. Children that harboured only the asexual form of P. falciparum were treated with antimalarial drugs (CoartemĀ©, Novartis) at school. Children that harboured gametocytes were driven together with their parents or their teacher to the laboratory of Kingani, where we took a sample of 3 ml of blood to infect mosquitoes (see below). They were then driven back to their villages and the children received an antimalarial treatment.

Mosquitoes

We used the mosquito An. arabiensis from a partially permethrin-resistant colony established in 2006 in Sakamaganga, Kilombero, Tanzania and maintained in an insectary at the Bagamoyo Research and Training Unit. Mosquitoes were maintained at a temperature of 27 Ā± 2 Ā°C with a relative humidity of 70ā€“80% and 12:12 h of dark:light photocycle. Larvae were reared in groups of 250 mosquitoes in trays containing 1 l of water and given the quantities of Tetramine fish food mentioned above. Adults were maintained in 35 Ɨ 35 Ɨ 35 cm cages and allowed to feed on a 10% glucose solution. Males were removed from the cages every 24 h to avoid fertilization. Twenty-four hours before the infection, female mosquitoes were starved and placed in groups of 50 into paper cups.

Blood preparation and infection of mosquitoes

Mosquitoes were fed on the blood of 5 children containing gametocytes of the human malaria parasite, Plasmodium falciparum, obtained as described above. The blood was prepared according to the procedure described by Bousema et al. [36]. The blood was centrifuged for 2 min at 3000Ɨ g. The serum was removed and replaced with a European AB serum, so that transmission-blocking factors were avoided [37]. Half of the blood of each child was maintained at 37 Ā°C; the other half was heated to 43 Ā°C for 15 min to kill the malaria parasites, preventing infection of mosquitoes without changing other blood characteristics [38, 39]. Mosquitoes were then allowed to feed in the dark for 45 min from micro-feeders containing 300 Ī¼l of infected or heat-killed blood maintained at 37 Ā°C. Twenty-four hours after infection, unfed mosquitoes were removed. The fed ones were placed in groups of 20 to 25 individuals into new cups and maintained in a 2.5 security laboratory at 25ā€“27 Ā°C and 70ā€“80% humidity, with 12:12 h of dark:light photocycle. They were allowed to feed on a piece of cotton soaked in a 10% glucose solution until 24 h before the behavioural tests. The tests were done 14 days after feeding, enabling the malaria parasite to produce sporozoites.

Measuring repellency and 24 h mortality

In this experiment we assessed whether the attraction of a host to a hungry (unfed) mosquito can overcome the repellency by the permethrin. We measured the repellency of the bednets to infectious or uninfected mosquitoes with a modified WHO tunnel (see Fig. 2). Mosquitoes were placed according to infection status and isolate in groups of 20 into one chamber (the release chamber) of the tunnel. KT placed his foot into the chamber on the opposite side (the stimulus chamber) to attract mosquitoes. His foot was protected from being bitten by a cloth with 150 Ī¼m mesh, which enables volatiles but not mosquitoes to pass through. The central chamber was separated from the stimulus chamber by a commercially available permethrin-treated mosquito net (Olyset net, Sumitomo Chemical Co. Ltd., Tokyo, Japan) or by an untreated mosquito net. The nets contained 7 holes (0.8 cm diameter), so that mosquitoes could pass through. The tunnel was covered with black tissue, for An. arabiensis prefers to bite during the night. We left mosquitoes in the tunnel for 20 min and then recorded the number of mosquitoes in each part of the tunnel. They were considered as unresponsive, if they remained in the release chamber; responsive but repelled if they were found in the central chamber; and responsive and not repelled if they were found behind the net in the stimulus chamber. After the test, mosquitoes were placed in paper cups for 24 h with access to 10% glucose solution and their 24 h mortality was recorded. They were then placed into Eppendorf tubes containing silica gel for transport to the University of NeuchĆ¢tel for further analysis.

Fig. 2
figure 2

Modified WHO tunnel to test the repellency of ITNs during the experiment with P. falciparum infecting An. arabiensis. Mosquitoes were placed in the released chamber and allowed to fly in the direction of the stimulus chamber where KT put one of his feet (protected by chiffon with 150 Ī¼m mesh). They had to pass through an opening made with cardboard with a rectangular hole of 8 Ɨ 12 cm to be considered as responsive, and then could pass through a net with holes to approach KTā€™s foot if they were not repelled. The arrows represent an example route for mosquitoes to go to the stimulus chamber

Malaria infection

Sporozoite-infection of mosquitoes was confirmed by PCR using the protocol established by Padley et al. [40]. DNA extraction of the head and thorax of mosquitoes was performed with DNAzol according to the procedure described above.

PCR was performed with the T3000 thermocycler (Analytik Jena AG). Forward (5'-AAC AGA CGG GTA GTC ATG ATT GAG-3') and reverse primers (5'-GTA TCT GAT CGT CTT CAC TCC C-3') amplified a 276 bp sequence of Plasmodium falciparum. PCR was done on 2 Ī¼l of DNA with 240 nM of each primer for a final volume of 25 Ī¼l. DNA was amplified during 40 cycles with denaturation at 95 Ā°C, annealing at 60 Ā°C and extension at 72 Ā°C, each for 45 s. Infection was confirmed from the correct band of DNA by electrophoresis on a 1.8% agarose gel stained with the RedSafeTM Nucleic Acid stain (iNtRON Biotechnology).

PCR amplifications were also done on the abdomens of mosquitoes that were exposed to infected blood in order to detect oocyst infection. We followed the same procedure for the PCR and the DNA extraction as for sporozoite-infection, except that we added only 150 Ī¼l of DNAzol for each abdomen.

Mosquito body size

The wing length was measured as the distance from the distal end of the alula to the tip of the wing (vein R3) without the fringes.

Statistical analysis

We estimated repellency in two ways: (i) as the proportion of mosquitoes that did not respond to the human cue (Additional file 2: Table S5), i.e. that stayed inside the releases chamber (Fig. 2); and (ii) as the proportion of responding mosquitoes that did not pass the net to move from the central to the stimulus chamber (Additional file 2: Table S4). The former is an expression of spatial repellency, the latter could be due to a combination of spatial repellency and contact irritancy. We analysed the likelihood to respond to human cues with a GLMM where the type of net and the infection status [unexposed (fed on blood containing heat-killed parasites); exposed but not infected; or infectious] were nominal factors, wing length was a covariate and the isolate (the child from whom the parasite had been obtained) and the replicate (the group of 20 mosquitoes tested simultaneously) were random factors. We analysed the likelihood that responding mosquitoes passed through the net with a GLMM that included the type of net and the infection status as nominal factors, wing length as a covariate, and the isolate and the replicate as random factors. Mosquitoes that were infected, but harboured only oocysts, were removed from the analyses. We further analysed to variations of this basic model. First, as the true infection status of ā€˜exposed but uninfectedā€™ mosquitoes was not clear, we left them out in a second analysis; this did not change the qualitative results. Secondly, we compared the repellency of the treated and untreated nets separately.

Mosquito mortality 24 h after the test of repellency was also analysed with a GLMM (Additional file 2: Table S4). The model included whether mosquitoes died or not as the response variable, and the same factors and covariate as the analysis for repellency. We also compared the mortality of each net individually.

All analyses were performed with the software R [33] and the RstudioĀ© interface (version 1.0.143) [34]. For the mixed effect models (which were done with the package ā€œlme4ā€ [35]), the statistical significance of each factor and covariate was determined with a type 3 ANOVA from the package ā€œcarā€ [32].

Results

Plasmodium berghei and Anopheles gambiae

Of the 80 mosquitoes exposed to gametocytic mice, 35 did not harbour sporozoites and were therefore excluded from the analysis. The percentage of mosquitoes responding to human cues was 55.6% (95% confidence interval, CI: 45.3ā€“65.4%) if they were uninfected and 66.7% (95% CI: 52.1ā€“78.6%) if they harboured sporozoites. This difference was not statistically significant (Ļ‡2 = 1.52, df = 1, P = 0.22).

Including only the mosquitoes that responded let us analyse the choice of 80 mosquitoes (30 infected and 50 uninfected) when facing a permethrin-treated or a control net. While 68% (95% CI: 54.2ā€“79.2%) of the uninfected mosquitoes avoided the cage with the insecticide, only 43.3% (95% CI: 27.4ā€“60.8%) of the infected ones did (Ļ‡2 = 4.11, df = 1, P = 0.043) (Fig. 3). Neither the mosquitoā€™s wing length (Ļ‡2 = 0.46, df = 1, P = 0.5) nor its age (Ļ‡2 = 0.06, df = 1, P = 0.81) affected the choice of mosquitoes. There was also no effect of the interaction between mosquitoā€™s age and their infection status (Ļ‡2 = 0.04, df = 1, P = 0.85).

Fig. 3
figure 3

Proportion of repelled mosquitoes as a function of their infection status. Mosquitoes were considered as repelled if they were found in the cage without the permethrin-treated net. The bars represent the 95% confidence interval

Plasmodium falciparum and Anopheles arabiensis

Of the mosquitoes that survived up to 14 days after blood-feeding, we selected haphazardly 80 mosquitoes from each of 3 isolates, 160 from a fourth one and 360 from a fifth (these numbers reflect the greater survival for the latter two isolates). We tested 50 mosquitoes that had fed on blood with heat-killed parasites; none of them were infected. Of the mosquitoes that had fed on infectious blood, 29.0% harboured sporozoites (half of these also harboured oocysts; in our analyses we did not distinguish between the mosquitoes with or without oocysts, both groups were considered as infectious). A further 6.3% harboured oocysts, but had no sporozoites; we removed these from the analysis.

The percentage of mosquitoes that responded to host cues (that is, that moved out of the release chamber) was 83.7% (95% CI: 79.6ā€“87.1%) for the unexposed ones (those fed on blood containing heat-killed parasites), 79.3% (95% CI: 73.8ā€“83.9%) for exposed-but-uninfected mosquitoes and 89.1% (95% CI: 81.9ā€“93.7%) for the infectious ones, (Ļ‡2 = 7.4, df = 2, P = 0.02). The difference of the percentage of responders between the untreated and treated nets was less than 2% for each infection status (main effect of net: Ļ‡2 = 0.01, df = 1, P = 0.92; interaction between infections status and net: Ļ‡2 = 0.13, df = 2, P = 0.94). Wing length had no effect on the proportion of responders (Ļ‡2 = 1.98, df = 2, P = 0.16).

The analysis comparing the proportion of mosquitoes passing the net into the stimulus chamber included 611 responding mosquitoes, with 286 mosquitoes (151 not exposed, 89 exposed-but-uninfected and 46 infectious) that were tested on an untreated net and 325 (167 not exposed, 106 exposed-but-uninfected and 52 infectious) mosquitoes that were tested on the permethrin-treated Olyset net. About of the quarter of the responding mosquitoes remained in the central chamber, if they were tested with an untreated net. This proportion was not affected by infection status: unexposed: 25.2% (95% CI: 18.9ā€“32.9%.); exposed-but-uninfected: 25.8% (95% CI: 17.9ā€“35.8%); infectious: 26.1% (95% CI: 15.6ā€“40.3%) (Ļ‡2 = 0.05, df = 2, P = 0.97). When tested on an Olyset net, in average 40.6% of mosquitoes were repelled, with the response ranging from 46.7% (95% CI: 39.3ā€“54.3%) of the unexposed mosquitoes down to 23.1% (95% CI: 13.7ā€“36.1%) of the infectious mosquitoes (Ļ‡2 = 7.67, df = 2, P = 0.02). Pairwise comparisons showed that unexposed mosquitoes were repelled more than sporozoite-infected ones (z = 2.71, P = 0.007), and that exposed-but-uninfected mosquitoes were not significantly different from either unexposed (z = 1.23, P = 0.22) or infectious ones (z = 1.74, P = 0.08).

The full analysis, pooling the data with the two nets, confirmed that the mosquitoes were less likely to pass through the Olyset than the untreated net (Ļ‡2 = 16.06, df = 1, P < 0.001), and that the effect of the net depended on the infection status. While the interaction between net type and infection was not quite significant when we compare all three infection groups (Ļ‡2 = 4.66, df = 2, P = 0.098), it was significant when we omitted the exposed-but-uninfected mosquitoes (Ļ‡2 = 4.44, df = 1, P = 0.035). Wing length did not affect the behavioural response in either of these analyses (three infection groups: Ļ‡2 = 0.02, df = 1, P = 0.90; unexposed vs infectious Ļ‡2 = 0.45, df = 1, P = 0.51).

Overall, uninfected mosquitoes were more repelled by the Olyset than the untreated net, but infectious mosquitoes were just as likely to pass through the Olyset as the untreated net (Fig. 4a).

Fig. 4
figure 4

a Proportion of mosquitoes found in the central cage (and therefore responded to the odour but were repelled by the net) for the two types of net [Olyset permethrin-treated net (light grey) or control untreated net (dark grey)] and for the three infection statuses. b Proportion of mosquitoes that died within 24 h after the behavioural test, as a function of the type of net and infection status. These proportions represent the overall mortality, combining the mortality of the mosquitoes considered as repelled and the one considered as unrepelled. For both graphs the bars represent the 95% confidence interval

Twenty-four hours after the test with the untreated net 23.1% of mosquitoes had died, independently of their infection status: unexposed 21.9% (95% CI: 16ā€“29.1%); exposed-but-uninfected: 24.7% (95% CI: 16.9ā€“34.6%); infectious: 23.9% (95% CI: 13.9ā€“37.9%) (Ļ‡2 = 0.55, df = 2, P = 0.76). After the test with a permethrin-treated net, 38.8% of mosquitoes had died, again independently of their infection status unexposed: 36.5% (95% CI: 29.6ā€“44.1%); exposed-but-uninfected: 38.7% (95% CI: 30ā€“48.2%); infectious: 46.2% (95% CI: 33.3ā€“59.5%) (Ļ‡2 = 1.45, df = 2, P = 0.48).

In the analysis pooling the two types of net, we found that the Olyset net increased the risk of dying by a factor of 1.7 (Ļ‡2 = 15.05, df = 1, P < 0.001, Fig. 4b) and we confirmed that there was no significant effect of the infection status (Ļ‡2 = 1.48, df = 2, P = 0.48) and that neither the wing length (Ļ‡2 = 0.0002, df = 1, P = 0.99) nor the interaction between the type of net and the infection status (Ļ‡2 = 0.46, df = 2, P = 0.79) affected the mosquitoā€™s mortality.

Discussion

Our results with two experimental systems show that infection by malaria affects the behavioural response of mosquitoes to insecticides. In particular, in both systems, infection by sporozoites negated the repellency of permethrin, therefore decreasing the personal protection offered by an ITN.

The decrease of repellency was not only observed with two combinations of malaria and mosquito species, but also with two ways of measuring the repellent effect of an insecticide. In the first experiment, when given the choice between flying into (and biting in) a chamber containing either an untreated or a treated net, uninfected mosquitoes preferred to fly into the chamber that was not protected by permethrin (showing that they were repelled by the insecticide), while infectious mosquitoes showed no preference. In the second, mosquitoes were not given the choice of an insecticide-free chamber, but we measured the likelihood that mosquitoes would pass an untreated or a treated net to get to their blood meal. Again, uninfected mosquitoes were repelled, and thus less likely to fly through a treated bednet than a treated one to reach their source of blood, while infectious mosquitoes responded similarly to treated and untreated nets.

Both of these aspects of repellency could be a consequence of spatial repellency (which impede mosquitoes from approaching treated nets) or contact irritancy (which makes mosquitoes fly away once they contact a treated net). However, our second experiment gives some evidence that spatial repellency is not an important factor. If it were important, we would have expected that fewer mosquitoes move towards a treated net (and thus stay in the release chamber) than towards an untreated net. The proportions, however, did not differ. This corroborated the patterns observed in several other studies that found no spatial repellency, but strong irritancy for permethrin (reviewed in [7]).

The proportion of mosquitoes remaining in the release chamber did, however, depend on their infection: sporozoite-infected mosquitoes were more likely to move out of the chamber and towards the human cues than uninfected or exposed-but-uninfected mosquitoes. This is in line with previously observed changes of mosquitoesā€™ biting behaviour induced by malaria, in particular that sporozoite-infection increases the host attraction [41] and the motivation and persistence of mosquitoes to obtain a blood meal [19,20,21,22,23] by enhancing the responsiveness of mosquitoā€™s olfactory receptors [41].

The possible mechanism of the decreased repellency of infectious mosquitoes by permethrin may be based on this behavioural manipulation: infectious mosquitoes could be more motivated to pass the permethrin-treated net to obtain their blood meal than uninfected ones. Alternatively, the chronic infection by malaria may simply damage in some way the receptors (or their integrative pathways) of the mosquitoes, analogously to the effect of a previous exposure to an insecticide in the odour receptors of Aedes aegypti [42]. Thus the mosquitoes would be less aware of the insecticideā€™s odour or irritancy, and therefore less repelled by it. Indeed, malaria infection reduces the levels of several proteins in the head of mosquitoes, including a synaptic and a neural signaling protein, though the functions of these proteins are not known [43].

Our observations on mortality are a bit less straightforward than those on behavior. As expected, the permethrin-treated net killed more mosquitoes than the untreated net. However, although the sporozoites increased the mortality from about 37% to about 46%, this difference was far from statistically significant. This contrasts our expectation that the increased contact of infectious mosquitoes with the insecticide would increase the likelihood that they are killed. A likely explanation for the lack of statistical significance is simply a lack of power. If (as the numbers of our experiment indicate) about 20% fewer mosquitoes are repelled and these are 10% more likely to be killed, we would need a large sample size to detect the combined effect as a statistically significant result. Alternatively, infectious mosquitoes may pass through the net so quickly that their exposure is too short to kill them. It would be important to know which of these (or any other) explanation is correct, for community protection and the rate of the evolution of insecticide resistance both depend mainly on the insecticidal property of ITNs [3]. Furthermore, to fully determine the level of community protection offer by ITNs it would be important to take in account the weaker attraction or motivation of oocyst-infected mosquitoes [23, 41]. Indeed, it might increase the overall number of infectious mosquitoes, for oocyst infection might decrease the mosquitoā€™s risk of dying due to ITNs or host defence behaviours [23, 44].

A caveat of our study is that both experiments were done in laboratory conditions, which can differ considerably from natural situations. In particular, laboratory infections might differ from natural infection and thus might influence the density of the sporozoites in the mosquitoesā€™ salivary glands [45,46,47], which in turn can influence the malaria-induced change of behaviour [24]. Nevertheless, we suggest that malaria infection in mosquitoes can decrease the repellency of insecticides and may reduce the personal protection at the stage when the protection is most relevant.

Conclusions

Insecticide-based tools, like ITNs and IRS, have proven their efficacy against malaria [2, 48] and remain the most cost-effective tools for malaria control [1]. However, to fully understand and predict the impact of ITNs on the epidemiology of malaria, we need to know how malaria affects the personal and community protection offered by the nets. Our main finding is that the malaria infection, at least in the laboratory, reduces the personal protection offered by insecticides. We suggest that further studies should investigate whether our laboratory results are consistent in the field and go on to investigate whether the less effective personal protection observed here can be balanced by increased community protection.

Abbreviations

CI:

Confidence interval

GLMM:

Generalized Linear Mixed Model

IRS:

Indoor residual spraying

ITN:

Insecticide-treated net

References

  1. Hanson K, Goodman C, Lines JD, Meek S, Bradley D, Mills A, et al. The economics of malaria control interventions. Geneva: World Health Organization; 2004. http://www.who.int/iris/handle/10665/43004. Accessed 20 September 2017.

  2. Steketee RW, Campbell CC. Impact of national malaria control scale-up programmes in Africa: magnitude and attribution of effects. Malar J. 2010;9:299.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  3. Birget PLG, Koella JC. An epidemiological model of the effects of insecticide-treated bed nets on malaria transmission. PLoS One. 2015;10:e0144173.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  4. Roberts D, Tren R. The excellent powder: DDTā€™s political and scientific history. Indianapolis: Dog Ear Publishing; 2010.

    Google ScholarĀ 

  5. Gatton ML, Chitnis N, Churcher T, Donnelly MJ, Ghani AC, Godfray HCJ, et al. The importance of mosquito behavioural adaptations to malaria control in Africa. Evolution. 2013;67:1218ā€“30.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  6. Grieco JP, Achee NL, Chareonviriyaphap T, Suwonkerd W, Chauhan K, Sardelis MR, et al. A new classification system for the actions of IRS chemicals traditionally used for malaria Control. PLoS One. 2007;2:e716.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  7. Chareonviriyaphap T Behavioral responses of mosquitoes to insecticides, Insecticides - Pest Engineering, Dr. Farzana Perveen (Ed.), InTech, https://doi.org/10.5772/28773. Available from: https://mts.intechopen.com/books/insecticides-pest-engineering/behavioral-responses-of-mosquitoes-to-insecticides. Accessed 20 Sept 2017.

  8. Atieli HE, Zhou G, Afrane Y, Lee M-C, Mwanzo I, Githeko AK, et al. Insecticide-treated net (ITN) ownership, usage, and malaria transmission in the highlands of western Kenya. Parasit Vectors. 2011;4:113.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  9. Kweka EJ, Nkya WM, Mahande AM, Assenga C, Mosha FW, Lyatuu EE, et al. Mosquito abundance, bed net coverage and other factors associated with variations in sporozoite infectivity rates in four villages of rural Tanzania. Malar J. 2008;7:59.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  10. Gimnig JE, Vulule JM, Lo TQ, Kamau L, Kolczak MS, Phillips-Howard PA, et al. Impact of permethrin-treated bed nets on entomologic indices in an area of intense year-round malaria transmission. Am J Trop Med Hyg. 2003;68:16ā€“22.

    PubMedĀ  Google ScholarĀ 

  11. Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, et al. Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010;9:62.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  12. Russell TL, Lwetoijera DW, Maliti D, Chipwaza B, Kihonda J, Charlwood JD, et al. Impact of promoting longer-lasting insecticide treatment of bed nets upon malaria transmission in a rural Tanzanian setting with pre-existing high coverage of untreated nets. Malar J. 2010;9:187.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  13. Deressa W, Yihdego YY, Kebede Z, Batisso E, Tekalegne A, Dagne GA. Effect of combining mosquito repellent and insecticide-treated net on malaria prevalence in southern Ethiopia: a cluster-randomised trial. Parasit Vectors. 2014;7:132.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  14. Chen-Hussey V, Carneiro I, Keomanila H, Gray R, Bannavong S, Phanalasy S, et al. Can topical insect repellents reduce malaria? A cluster-randomised controlled trial of the insect repellent N,N-diethyl-m-toluamide (DEET) in Lao PDR. PLoS One. 2013;8:e70664.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  15. Gu W, Novak RJ. Predicting the impact of insecticide-treated bed nets on malaria transmission: the devil is in the detail. Malar J. 2009;8:256.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  16. Smallegange RC, van Gemert G-J, van de Vegte-Bolmer M, Gezan S, Takken W, Sauerwein RW, et al. Malaria-infected mosquitoes express enhanced attraction to human odor. PLoS One. 2013;8:e63602.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  17. Cator LJ, Lynch PA, Read AF, Thomas MB. Do malaria parasites manipulate mosquitoes? Trends Parasitol. 2012;28:466ā€“70.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  18. Vantaux A, DabirĆ© KR, Cohuet A, LefĆØvre T. A heavy legacy: offspring of malaria-infected mosquitoes show reduced disease resistance. Malar J. 2014;13:442.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  19. Koella JC, SĆørensen FL, Anderson RA. The malaria parasite, Plasmodium falciparum, increases the frequency of multiple feeding of its mosquito vector, Anopheles gambiae. Proc R Soc Lond B Biol Sci. 1998;265:763ā€“8.

    ArticleĀ  CASĀ  Google ScholarĀ 

  20. Wekesa JW, Copeland RS, Mwangi RW. Effect of Plasmodium falciparum on blood feeding behavior of naturally infected Anopheles mosquitoes in western Kenya. Am J Trop Med Hyg. 1992;47:484ā€“8.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  21. Anderson RA, Koella JC, Hurd H. The effect of Plasmodium yoelii nigeriensis infection on the feeding persistence of Anopheles stephensi Liston throughout the sporogonic cycle. Proc R Soc Lond B Biol Sci. 1999;266:1729ā€“33.

    ArticleĀ  CASĀ  Google ScholarĀ 

  22. Rossignol PA, Ribeiro JM, Spielman A. Increased biting rate and reduced fertility in sporozoite-infected mosquitoes. Am J Trop Med Hyg. 1986;35:277ā€“9.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  23. Koella JC, Rieu L, Paul REL. Stage-specific manipulation of a mosquitoā€™s host-seeking behavior by the malaria parasite Plasmodium gallinaceum. Behav Ecol. 2002;13:816ā€“20.

    ArticleĀ  Google ScholarĀ 

  24. Koella JC. An evolutionary view of the interactions between anopheline mosquitoes and malaria parasites. Microbes Infect. 1999;1:303ā€“8.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  25. World Health Organization. Global insecticide use for vector-borne disease control. 4th edn. Geneva: World Health Organization; 2009.

  26. Kulma K, Saddler A, Koella JC. Effects of age and larval nutrition on phenotypic expression of insecticide-resistance in Anopheles mosquitoes. PLoS One. 2013;8:e58322.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  27. Vanderberg JP, Yoeli M. Effects of temperature on sporogonic development of Plasmodium berghei. J Parasitol. 1966;52:559ā€“64.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  28. Cornet S, Nicot A, Rivero A, Gandon S. Both infected and uninfected mosquitoes are attracted toward malaria infected birds. Malar J. 2013;12:179.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  29. Rider MA, Byrd BD, Keating J, Wesson DM, Caillouet KA. PCR detection of malaria parasites in desiccated Anopheles mosquitoes is uninhibited by storage time and temperature. Malar J. 2012;11:193.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  30. Schwartz A, Koella JC. Melanization of Plasmodium falciparum and C-25 sephadex beads by field-caught Anopheles gambiae (Diptera: Culicidae) from southern Tanzania. J Med Entomol. 2002;39:84ā€“8.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  31. Koella JC, SĆørense FL. Effect of adult nutrition on the melanization immune response of the malaria vector Anopheles stephensi. Med Vet Entomol. 2002;16:316ā€“20.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  32. Fox J, Weisberg S. An R companion to applied regression, second edition. Thousand Oaks: Sage; 2011.

    Google ScholarĀ 

  33. R Core Team. R: A language and environment for statistical computing. 2013. https://www.R-project.org.

    Google ScholarĀ 

  34. RStudio Team. RStudio: Integrated development environment for R. 2016. http://www.rstudio.com.

    Google ScholarĀ 

  35. Bates D, Maechler M, Bolker B, Walker S, Fitting linear mixed-effects models using lme4. J Stat Software. 2015:67:jss.v067.i01

  36. Bousema T, Churcher TS, Morlais I, Dinglasan RR. Can field-based mosquito feeding assays be used for evaluating transmission-blocking interventions? Trends Parasitol. 2013;29:53ā€“9.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  37. Boudin C, Diop A, Gaye A, Gadiaga L, Gouagna C, Safeukui I, et al. Plasmodium falciparum transmission blocking immunity in three areas with perennial or seasonal endemicity and different levels of transmission. Am J Trop Med Hyg. 2005;73:1090ā€“5.

    PubMedĀ  Google ScholarĀ 

  38. Sangare I, Michalakis Y, Yameogo B, Dabire R, Morlais I, Cohuet A. Studying fitness cost of Plasmodium falciparum infection in malaria vectors: validation of an appropriate negative control. Malar J. 2013;12:2.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  39. Mendes AM, Schlegelmilch T, Cohuet A, Awono-Ambene P, Iorio MD, Fontenille D, et al. Conserved mosquito/parasite interactions affect development of Plasmodium falciparum in Africa. PLoS Pathog. 2008;4:e1000069.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  40. Padley D, Moody AH, Chiodini PL, Saldanha J. Use of a rapid, single-round, multiplex PCR to detect malarial parasites and identify the species present. Ann Trop Med Parasitol. 2003;97:131ā€“7.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  41. Cator LJ, George J, Blanford S, Murdock CC, Baker TC, Read AF, et al. ā€˜Manipulationā€™ without the parasite: altered feeding behaviour of mosquitoes is not dependent on infection with malaria parasites. Proc R Soc Lond B Biol Sci. 2013;280:20130711.

    ArticleĀ  Google ScholarĀ 

  42. Stanczyk NM, Brookfield JFY, Field LM, Logan JG. Aedes aegypti mosquitoes exhibit decreased repellency by DEET following previous exposure. PLoS One. 2013;8:e54438.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  43. Lefevre T, Thomas F, Schwartz A, Levashina E, Blandin S, Brizard J-P, et al. Malaria Plasmodium agent induces alteration in the head proteome of their Anopheles mosquito host. Proteomics. 2007;7:1908ā€“15.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  44. Ohm JR, Teeple J, Nelson WA, Thomas MB, Read AF, Cator LJ. Fitness consequences of altered feeding behavior in immune-challenged mosquitoes. Parasit Vectors. 2016;9:113.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  45. Awono-Ambene HP, Diawara L, Robert V. Comparison of direct and membrane feeding methods to infect Anopheles arabiensis with Plasmodium falciparum. Am J Trop Med Hyg. 2001;64:32ā€“4.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  46. Sattabongkot J, Maneechai N, Phunkitchar V, Eikarat N, Khuntirat B, Sirichaisinthop J, et al. Comparison of artificial membrane feeding with direct skin feeding to estimate the infectiousness of Plasmodium vivax gametocyte carriers to mosquitoes. Am J Trop Med Hyg. 2003;69:529ā€“35.

    PubMedĀ  Google ScholarĀ 

  47. Bousema T, Dinglasan RR, Morlais I, Gouagna LC, van WT, Awono-Ambene PH, et al. Mosquito feeding assays to determine the infectiousness of naturally infected Plasmodium falciparum gametocyte carriers. PLoS One. 2012;7:e42821.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  48. Lengeler C. Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004;2:CD000363.

    Google ScholarĀ 

  49. Wolfensohn S, Lloyd M. Handbook of laboratory animal management and welfare. 4th ed. New Jersey: Wiley-Blackwell; 2013.

    Google ScholarĀ 

Download references

Acknowledgements

We thank Professor Volker Heusslerā€™s team (University of Bern) for providing the mice and the malaria parasites, Olivier Rais for dealing with the mice in our laboratory experiments, and Susana Dos Santos for her help in the lab. We also thank the staff from the Ifakara Health Institute for their help during the work in the field and in their lab.

Funding

Financial support for the field experiment in Tanzania was provided by the Fonds des Donations from the University of NeuchĆ¢tel, the Fondation Hans Wilsdorf, the Schmidheiny foundation and the Rotary Club GenĆØve Groupe SERVIR (Marc Plojoux).

Availability of data and materials

All data generated or analysed during this study are included in this published article and its additional files.

Author information

Authors and Affiliations

Authors

Contributions

KT, SM and JCK conceived and designed the study. SM and JCK coordinated the study. KT, LH, ER and MMT performed the laboratory and field work and collected the data. KT and JCK carried out the statistical analysis. KT and JCK wrote the original draft and all authors contributed to the elaboration of the final version of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kevin ThiƩvent.

Ethics declarations

Ethics approval and consent to participate

All procedures that involved mice were conducted in accordance with the Swiss Tierschutzgesetz guidelines (TSchG; Animal Rights Laws). They were approved by the ethical committee of the University of Bern (Permit Number: BE109/13). Mice were anesthetized with an intra-peritoneal injection of 8.5 ml/kg of a mix of Xylazine XylasolĀ® (solution: 20mg/ml), Ketamine Ketasol Ā® (solution: 100 mg/ml) and PBS [49]. Obtaining and using the blood of the Tanzanian schoolchildren was approved by the Ethical Committee of the University of NeuchĆ¢tel, the Ifakara Health Research and Development Centre, the Medical Research Coordination Committee (MRCC) of the National Institute of Medical Research (NIMR) and the Tanzanian Commission for Science and Technology (COSTECH). Parents were invited to attend a sensitization meeting at the school premises. The children as well as their parents and teachers were informed about the project, and those who were interested to participate were given an informed consent form. The parents had to read and sign the form in order to allow us to enroll their children in the study.

Competing interests

The authors declare that they have no competing interests.

Publisherā€™s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. P. berghei - An. gambiae dataset. Data from the laboratory experiment with An. gambiae and P.berghei used for statistical analysis of repellency. Table S2. P. berghei - An. gambiae dataset. Data from the laboratory experiment with An. gambiae and P.berghei used for statistical analysis of the proportion of responding mosquitoes. Table S3. P. berghei - An. gambiae dataset. Data from the laboratory experiment with An. gambiae and P.Ā berghei used for statistical analysis of the side preference. (XLSX 16 kb)

Additional file 2:

Table S4. P. falciparum - An. arabiensis dataset. Data from the experiment with An. arabiensis and P. falciparum used for statistical analysis of repellency. Table S5. P. falciparum - An. arabiensis dataset. Data from the experiment with An. arabiensis and P. falciparum used for statistical analysis of the proportion of responding mosquitoes. (XLSX 50 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

ThiƩvent, K., Hofer, L., Rapp, E. et al. Malaria infection in mosquitoes decreases the personal protection offered by permethrin-treated bednets. Parasites Vectors 11, 284 (2018). https://doi.org/10.1186/s13071-018-2846-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-018-2846-0

Keywords