Effects of pyrethroid resistance on the cost effectiveness of a mass distribution of long-lasting insecticidal nets: a modelling study
© Briët et al.; licensee BioMed Central Ltd. 2013
Received: 4 January 2013
Accepted: 22 February 2013
Published: 25 February 2013
The effectiveness of insecticide-treated nets in preventing malaria is threatened by developing resistance against pyrethroids. Little is known about how strongly this affects the effectiveness of vector control programmes.
Data from experimental hut studies on the effects of long-lasting, insecticidal nets (LLINs) on nine anopheline mosquito populations, with varying levels of mortality in World Health Organization susceptibility tests, were used to parameterize malaria models. Both simple static models predicting population-level insecticidal effectiveness and protection against blood feeding, and complex dynamic epidemiological models, where LLINs decayed over time, were used. The epidemiological models, implemented in OpenMalaria, were employed to study the impact of a single mass distribution of LLINs on malaria, both in terms of episodes prevented during the effective lifetime of the batch of LLINs, and in terms of net health benefits (NHB) expressed in disability-adjusted life years (DALYs) averted during that period, depending on net type (standard pyrethroid-only LLIN or pyrethroid-piperonyl butoxide combination LLIN), resistance status, coverage and pre-intervention transmission level.
There were strong positive correlations between insecticide susceptibility status and predicted population level insecticidal effectiveness of and protection against blood feeding by LLIN intervention programmes. With the most resistant mosquito population, the LLIN mass distribution averted up to about 40% fewer episodes and DALYs during the effective lifetime of the batch than with fully susceptible populations. However, cost effectiveness of LLINs was more sensitive to the pre-intervention transmission level and coverage than to susceptibility status. For four out of the six Anopheles gambiae sensu lato populations where direct comparisons between standard LLINs and combination LLINs were possible, combination nets were more cost effective, despite being more expensive. With one resistant population, both net types were equally effective, and with one of the two susceptible populations, standard LLINs were more cost effective.
Despite being less effective when compared to areas with susceptible mosquito populations, standard and combination LLINs are likely to (still) be cost effective against malaria even in areas with strong pyrethroid resistance. Combination nets are likely to be more cost effective than standard nets in areas with resistant mosquito populations.
KeywordsPyrethroid Insecticide Resistance LLIN ITN Modelling Piperonyl butoxide
Malaria control is based on both preventing transmission and promptly and effectively treating infection. Many countries have made significant progress in preventing malaria by focusing largely on vector control through long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) of insecticides . In areas with insecticide-susceptible mosquito populations, the insecticide on LLINs mitigates the loss of personal protection if the net becomes holed. Furthermore, the nets have a community effect by reducing the longevity of malaria vector mosquitoes . Over the last decade, many countries have significantly increased LLIN coverage with great impact: it is estimated that between 2000 and 2010, LLINs saved over 908,000 lives, with three quarters of those deaths having been prevented since 2006 . However, both IRS and LLINs face the development of physiological resistance (against insecticide) and ‘behavioural resistance’ in mosquitoes, which can reduce effectiveness of these interventions and possibly reverse the gains made in reducing malaria morbidity . Although such resistance may be inevitable with successful control programmes, new strategies need to be developed to mitigate development and spread of insecticide resistance and to preserve the efficacy of currently available insecticides and the effectiveness of malaria control interventions.
The problem of physiological resistance against insecticides is more acute for LLINs than for IRS, as LLINs rely solely on pyrethroids, whereas IRS can be done with several classes of insecticides. There is evidence of increasing levels of pyrethroid resistance  and corresponding decreases in the effectiveness of malaria control programmes that rely on pyrethroid-based interventions [6, 7]. Although LLINs containing non-pyrethroid classes of insecticides are under development, no such net currently exists. However, to mitigate the effects of pyrethroid resistance, nets using a combination of a pyrethroid and a synergist (piperonyl butoxide) have been developed.
The aim of this study is to assess the sensitivity of the effectiveness of a mass distributed batch of LLINs to insecticide resistance, both in terms of episodes averted during the effective life time of the batch, and in terms of net health benefits (NHB) expressed in disability adjusted life years (DALYs) averted.
New Bussa AG
New Bussa AA
Vand Duc A
For all these mosquito populations, except ‘Akron’ (population 3), results were available from World Health Organization (WHO) susceptibility test  assays with 0.05% deltamethrin (Table 1) conducted at the time of the experimental hut studies. However, for ‘Akron’, as well as for ‘New Bussa AA’, ‘New Bussa AG’ and ‘Yaokoffikro’, data on susceptibility to 0.75% permethrin was available, and this permitted estimating of the mortality of ‘Akron’ if exposed to 0.05% deltamethrin in WHO tests, assuming a linearb relationship between logit transformed mortality proportions of permethrin and deltamethrin. The nine mosquito populations (populations 1–9 in Table 1) varied in their susceptibility to deltamethrin over a wide range from susceptible to resistant.
Summarised experimental hut data on the numbers of unfed alive (UA), unfed dead (UD), fed alive (FA) and fed dead (FD) female mosquitoes for selected study arms are shown in Additional file 1. Some of these data were not previously published in as much detail. For populations 6 and 7, the data are for the first six weeks of data collection, whereas the reference available for this study  summarises the total of 12 weeks of data collection, and gives detailed results for ‘New Bussa AG’ (population 6 in Table 1) and a Culex quinquefasciatus population not used here, but not for ‘New Bussa AA’ (population 7 in Table 1). For population 8, the study  deviates from the standard protocol, due to ethical concerns, in that intact rather than purposely holed nets were used. For populations 1–7 and 9, studies [8–10, 12, 13] do not have study arms with intact nets. To assess the cost effectiveness of a mass distributed batch of LLINs, in the simulations, nets were presumed to be delivered intact and develop holes over time due to wear and tear. Because the studies test either artificially holed or intact nets but not both, additional data from experimental hut studies [15, 16] elsewhere (populations 10 and 11 in Additional files 1 and 2) were used to estimate effects, depending on the physical state in terms of holed area in the fabric surface.
For populations 1–9, and 10, studies do not have controls with unprotected individuals (i.e., a person without a net). For population 10, only summarised data were available and numbers for the four mosquito categories (FA, FD, UD, UA) were estimated assuming that the number of FD was the average of the number of fed mosquitoes multiplied by the proportion of FD out of fed mosquitoes and the number of dead mosquitoes multiplied by the proportion of FD out of dead mosquitoes in the corresponding arm of the study on population 11.
For the studies on populations 1–4 and 6–7, insecticide content after completion of the experimental hut studies in side panels in g/kg was converted into mg/m2 units. The conversion factor was 55/1.8 = 30.6 for 75 denier nets was and 55/1.4 = 39.3c for 100 denier nets and were derived from WHO specifications . For the study on populations 5 and 9, this was already in mg/m2. For the study on population 8, no insecticide content data were available. Therefore, these were estimated as follows: for the unwashed nets, the insecticide content was estimated as the mean target dose (P2: 1.8 g/kg; P3: 2.8 g/kg), multiplied by the average percentage retention after completion of the experimental hut essays, found in the studies on populations 1–5 (P2: 91.7%; P3: 95.7%); for the 20-times washed nets, the target dose was multiplied by the average percentage retention found in the studies on populations 1–5 after washing and completion of the experimental hut studies (P2: 27.4%; P3: 24.0%).
Parameter values for the effects of LLINs on mosquito populations needed in the mathematical models were estimated according to the results from these experimental hut studies (see Additional files 1 and 2). These parameter values were used both for analysing the relationship between insecticide resistance status and the population level insecticidal effectiveness and protection against feeding  of new and decayed LLINs, and for simulation studies on the malaria epidemiological effectiveness of LLINs.
From each scenario output, based on the number of uncomplicated episodes, severe episodes, sequelae and deaths depending on age, DALYs and health system costs (HSC) were calculated for each whole year during the simulation run (see Additional file 2). For intervention scenarios, the ‘effective lifetime’ (the period of the epidemiological effect from the year of distribution until the year where the effect had waned to just over half of the maximum impact achieved in terms of uncomplicated and severe episodes averted) as compared to non-intervention scenarios, was determined. This period is illustrated by the red arrow in Figure 1b. For each scenario and corresponding non-intervention scenario, the outcome variables during this period were averaged. This allowed calculation of the effectiveness of a mass distributed batch of LLINs in terms of episodes averted during the effective lifetime, and cost effectiveness in terms of NHB. The NHBs were calculated using a value of 1 DALY = 235.28 (2012) USD, as described in Additional file 3.
Results and discussion
Insecticide resistance status and LLIN effects in static models
The strong correlation between the population level direct protection against feeding and the mortality in WHO susceptibility assays is striking: whereas WHO susceptibility assays cannot assess behavioural resistance, results derived from experimental hut assays are affected by physiological or biochemical resistance, deterrence and dissuasion from attacking. These last two effects, resulting in avoidance of contact with the insecticide, could be seen as a form of behavioural resistance, even though it is possibly effective in reducing host-vector contact. Note that these relationships change as LLINs decay (first column of panels versus second column of panels), but that this is unimportant relative to the differences observed among populations.
Although the effectiveness of an intervention programme against malaria transmission depends both on protection against bites and the insecticidal effectiveness, more complex models, such as those implemented in OpenMalaria, are required to estimate impacts on transmission and morbidity. Nevertheless, these relatively simple measures provide some predictions against which the predictions of the more complex models can be checked. Since Figure 2 shows strong correlation with susceptibility status and both protection against feeding and insecticidal effectiveness, such a result can also be expected for transmission predictions by the more complex models. For new intact nets, where comparisons between P2 and P3 are possible, it appears that P3 is more effective than P2 in most situations, except for the fully susceptible population ‘Zeneti’: both protection against feeding and insecticidal effectiveness were different for the two LLIN types (Wilcoxon signed rank test, alpha = 0.05). For decayed nets, protection against feeding by P3 was higher than P2 in four out of seven populations, and insecticidal effectiveness of P3 was higher than P2 in five out of seven populations, but these differences were not significant.
LLIN effectiveness depending on transmission level
Insecticide resistance status and LLIN effectiveness depending on transmission level
LLIN effectiveness depending on LLIN type
The difference between the two LLIN types, up to about 15% in the most resistant population, was generally small, compared to the effect of insecticide susceptibility on LLIN effectiveness (Figure 4). The strongest differences were found at transmission levels between eight and 32 IBPAPA.
With three out of four resistant populations (‘Akron’, ‘Yaokoffikro’ and ‘Pitoa’), P3 was more effective, whereas with population ‘Kou’ there was no apparent difference. With one fully susceptible population (‘Malanville’), P3 was more effective, whereas with the fully susceptible population (‘Zeneti’), P2 was more effective.
Simple static models that predict population level insecticidal effectiveness and protection against blood feeding, and complex dynamic epidemiological models agreed that insecticide susceptibility, as measured by mortality in WHO susceptibility tests, is strongly correlated with LLIN effectiveness against malaria. An earlier modelling study by Killeen and colleagues  on the effectiveness of P3 and P2 LLINs in reducing malaria transmission provides a ‘snapshot’ of how transmission might be affected with nets in a certain physical and chemical state, when the mosquito population structure is in an equilibrium state. In the present simulation experiment, the LLINs were followed from a new and intact state along the path of attrition and physical and chemical decay, while considering the interacting dynamics of malaria in humans and mosquitoes. Nevertheless, both studies agree that, taking community level effects into account, both P3 and P2 nets are (still) effective in suppressing malaria in areas with resistant mosquito populations. Despite being more expensive, P3 was more cost effective than P2 against three out of four resistant populations. This could be due to the PBO present in the roofs of P3 nets, but effects of the higher deltamethrin content of P3 nets cannot be excluded . However, for the areas where mosquitoes are (still) fully susceptible, results were inconclusive, with P3 being more cost effective in one area and P2 being more cost effective in another. Despite resistance being an important factor in reducing the effectiveness of LLINs, both when measured in terms of averted episodes and in terms of NHB, coverage level and especially the pre-intervention transmission level appear still more important. In settings with moderate pre-intervention transmission, a mass distributed batch of LLINs was cost effective against malaria even in the presence of strong physiological or biochemical resistance. However, at pre-intervention transmission levels above 128 IBPAPA, a minority of variants of the model ensemble showed (slightly) negative net health benefits, even with fully susceptible mosquitoes. Moreover, it should be noted that in most malaria endemic countries, malaria control will not be limited to a single mass distribution of LLINs, and cost-effectiveness may alter with prolonged LLIN use.
Nevertheless, standard and combination LLINs are likely to be cost effective against malaria even in areas with strong pyrethroid resistance, and LLINs remain an effective transmission control paradigm in the fight against malaria.
aThe earlier comparative modelling study by Killeen and colleagues  used data for four populations (Pitoa, Kou, Akron, and Van Duc A, there named ‘Cameroon’, ‘Burkina Faso’, ‘Benin’, and ‘Vietnam’, respectively) of the nine used here, to study the effect of holed P3 and P2 nets on human exposure to infectious bites.
bDeltamethrin and permethrin cross resistance is likely mechanism dependent and the assumption of a linear relationship could be biased, since in population ‘Akron’, esterase, kdr and oxidase mechanisms were reported, whereas in populations ‘New Bussa AG’ and ‘Yakoffikro’ only kdr and oxidase were reported. It seems reasonable to assume that the mortality of ‘Akron’ to deltamethrin would be less than that of population ‘Yaokoffikro’, based on lower mortality of An. gambiae s.s. in Asecna and Ladji  (neighbouring Akron) than in Yaokoffikro [11, 28].
cPolyester, 75 denier netting weighs approximately 30 g/m2 +/− 10%, and 100 denier weighs approximately 40 g/m2 +/− 10%.
dNHBs are only slightly lower than DALYs averted during the effective lifetime of a mass distributed batch of LLINs, because the health systems cost of a net distribution are relatively low compared to DALYs averted if expressed in terms of DALYs (see Additional file 3).
We acknowledge our colleagues in the Swiss TPH Health Systems Research and Dynamic Modeling Unit, particularly Michael Tarantino, for help with the models and simulations. Simulations were run with the help of volunteers providing computational power. Gerry Killeen provided helpful comments on an early version of the manuscript. This study was financially supported by Vestergaard Frandsen.
- Global Malaria Programme WHO: World Malaria Report 2011. 2011, Geneva: World Health Organization
- Hawley WA, Phillips-Howard P, ter Kuile F, Terlouw DJ, Vulule JM, Ombok M, Nahlen B, Gimnig JE, Kariuki SK, Kolczak MS, Hightower AW: Community-wide effects of permethrin-treated bed nets on child mortality and malaria morbidiy in western Kenya. Am J Trop Med Hyg. 2003, 68 (Suppl 4): 121-127.PubMed
- Johansson EW, Newby H, Steketee RW: World Malaria Day 2010: Africa update. 2010, Geneva: UNICEF, [Roll Back Malaria Progress & Impact Series, vol 2.]
- Moiroux N, Gomez MB, Pennetier C, Elanga E, Djenontin A, Chandre F, Djegbe I, Guis H, Corbel V: Changes in Anopheles funestus biting behavior following universal coverage of long-lasting insecticidal nets in Benin. J Infect Dis. 2012, 206: 1622-1629. 10.1093/infdis/jis565.View ArticlePubMed
- Ranson H, N'Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V: Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control?. Trends Parasitol. 2011, 27: 91-98. 10.1016/j.pt.2010.08.004.View ArticlePubMed
- Trape JF, Tall A, Diagne N, Ndiath O, Ly AB, Faye J, Dieye-Ba F, Roucher C, Bouganali C, Badiane A, Sarr FD, Mazenot C, Touré-Baldé A, Raoult D, Druilhe P, Mercereau-Puijalon O, Rogier C, Sokhna C: Malaria morbidity and pyrethroid resistance after the introduction of insecticide-treated bednets and artemisinin-based combination therapies: a longitudinal study. Lancet Infect Dis. 2011, 11: 925-932. 10.1016/S1473-3099(11)70194-3.View ArticlePubMed
- Corbel V, Akogbeto M, Damien GB, Djenontin A, Chandre F, Rogier C, Moiroux N, Chabi J, Banganna B, Padonou GG, Henry MC: Combination of malaria vector control interventions in pyrethroid resistance area in Benin: a cluster randomised controlled trial. Lancet Infect Dis. 2012, 12: 617-626. 10.1016/S1473-3099(12)70081-6.View ArticlePubMed
- Corbel V, Chabi J, Dabire RK, Etang J, Nwane P, Pigeon O, Akogbeto M, Hougard JM: Field efficacy of a new mosaic long-lasting mosquito net (PermaNet® 3.0) against pyrethroid-resistant malaria vectors: a multi centre study in Western and Central Africa. Malar J. 2010, 9: 113. 10.1186/1475-2875-9-113.PubMed CentralView ArticlePubMed
- N'Guessan R, Asidi A, Boko P, Odjo A, Akogbeto M, Pigeon O, Rowland M: An experimental hut evaluation of PermaNet® 3.0, a deltamethrin-piperonyl butoxide combination net, against pyrethroid-resistant Anopheles gambiae and Culex quinquefasciatus mosquitoes in southern Benin. Trans R Soc Trop Med Hyg. 2010, 104: 758-765. 10.1016/j.trstmh.2010.08.008.View ArticlePubMed
- Tungu P, Magesa S, Maxwell C, Malima R, Masue D, Sudi W, Myamba J, Pigeon O, Rowland M: Evaluation of PermaNet 3.0 a deltamethrin-PBO combination net against Anopheles gambiae and pyrethroid resistant Culex quinquefasciatus mosquitoes: an experimental hut trial in Tanzania. Malar J. 2010, 9: 21. 10.1186/1475-2875-9-21.PubMed CentralView ArticlePubMed
- Koudou BG, Koffi AA, Malone D, Hemingway J: Efficacy of PermaNet® 2.0 and PermaNet® 3.0 against insecticide-resistant Anopheles gambiae in experimental huts in Côte d'Ivoire. Malar J. 2011, 10: 172. 10.1186/1475-2875-10-172.PubMed CentralView ArticlePubMed
- Adeogun AO, Olojede JB, Oduola AO, Awolola TS: Efficacy of a combination long lasting insecticidal net (PermaNet® 3.0) against pyrethroid resistant Anopheles gambiae s.s. and Culex quinquefasciatus: an experimental hut trial in Nigeria. Nig J Clin Biomed Res. 2012, 6: 37-50.
- Van Bortel W, Chinh VD, Berkvens D, Speybroeck N, Trung HD, Coosemans M: Impact of insecticide-treated nets on wild pyrethroid resistant Anopheles epiroticus population from southern Vietnam tested in experimental huts. Malar J. 2009, 8: 248. 10.1186/1475-2875-8-248.PubMed CentralView ArticlePubMed
- World Health Organization: Guidelines for Testing Mosquito Adulticides for Indoor Residual Spraying and Treatment of Mosquito Nets. 2006, Geneva: WHO/CDS/NTD/WHOPES/GCDPP/2006.3
- Lines JD, Curtis CF, Myamba J, Njau R: Tests of Repellent or Insecticide Impregnated Curtains, Bednets and Anklets Against Malaria Vectors in Tanzania. 1985, Geneva: World Health Organization, WHO/VBC/85.920
- Irish S, N'Guessan R, Boko P, Metonnou C, Odjo A, Akogbeto M, Rowland M: Loss of protection with insecticide-treated nets against pyrethroid-resistant Culex quinquefasciatus mosquitoes once nets become holed: an experimental hut study. Parasit Vectors. 2008, 1: 17. 10.1186/1756-3305-1-17.PubMed CentralView ArticlePubMed
- World Health Organization: Deltamethrin Long-lasting (Coated onto Filaments) Insecticidal Net (S)-a-cyano-3-phenoxybenzyl (1R,3R)-3-(2,2-dibromovinyl)-2, 2-dimethylcyclopropane carboxylate. 2011, Geneva: [WHO Specifications and Evaluations for Public Health Pesticides]
- Briet OJ, Smith TA, Chitnis N: Measurement of overall insecticidal effects in experimental hut trials. Parasit Vectors. 2012, 5: 256. 10.1186/1756-3305-5-256.PubMed CentralView ArticlePubMed
- OpenMalaria: A simulator of malaria epidemology and control: [http://code.google.com/p/openmalaria/]
- Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M: Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: overview. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 1-10.PubMed
- Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, Studer A, Genton B, Lengeler C, Tediosi F, de Savigny D, Tanner M: Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008, 135: 1507-1516. 10.1017/S0031182008000371.View ArticlePubMed
- Chitnis N, Hardy D, Smith T: A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012, 74: 1098-1124. 10.1007/s11538-011-9710-0.PubMed CentralView ArticlePubMed
- Briët OJT, Hardy D, Smith TA: Importance of factors determining the effective lifetime of a mass, long-lasting, insecticidal net distribution: a sensitivity analysis. Malar J. 2012, 11: 20. 10.1186/1475-2875-11-20.PubMed CentralView ArticlePubMed
- Erhart A, Thang ND, Bien TH, Tung NM, Hung NQ, Hung LX, Tuy TQ, Speybroeck N, Cong LD, Coosemans M, D'Alessandro U: Malaria epidemiology in a rural area of the Mekong Delta: a prospective community-based study. Trop Med Int Health. 2004, 9: 1081-1090. 10.1111/j.1365-3156.2004.01310.x.View ArticlePubMed
- Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, Brooks A, Penny M, Tanner M: Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012, 9: e1001157. 10.1371/journal.pmed.1001157.PubMed CentralView ArticlePubMed
- Killeen GF, Okumu FO, N'Guessan R, Coosemans M, Adeogun A, Awolola S, Etang J, Dabire RK, Corbel V: The importance of considering community-level effects when selecting insecticidal malaria vector products. Parasit Vectors. 2011, 4: 160. 10.1186/1756-3305-4-160.PubMed CentralView ArticlePubMed
- Corbel V, N'Guessan R, Brengues C, Chandre F, Djogbenou L, Martin T, Akogbeto M, Hougard JM, Rowland M: Multiple insecticide resistance mechanisms in Anopheles gambiae and Culex quinquefasciatus from Benin, West Africa. Acta Trop. 2007, 101: 207-216. 10.1016/j.actatropica.2007.01.005.View ArticlePubMed
- Koffi AA, houa Alou LP, Adja MA, Kone M, Chandre F, N'Guessan R: Update on resistance status of Anopheles gambiae s.s. to conventional insecticides at a previous WHOPES field site, "Yaokoffikro", 6 years after the political crisis in Côte d'Ivoire. Parasit Vectors. 2012, 5: 68. 10.1186/1756-3305-5-68.PubMed CentralView ArticlePubMed
- Djouaka RF, Bakare AA, Coulibaly ON, Akogbeto MC, Ranson H, Hemingway J, Strode C: Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria. BMC Genomics. 2008, 9: 538. 10.1186/1471-2164-9-538.PubMed CentralView ArticlePubMed
- Kitau J, Oxborough RM, Tungu PK, Matowo J, Malima RC, Magesa SM, Bruce J, Mosha FW, Rowland MW: Species shifts in the Anopheles gambiae complex: Do LLINs successfully control Anopheles arabiensis?. PLoS One. 2012, 7: e31481. 10.1371/journal.pone.0031481.PubMed CentralView ArticlePubMed
- Chandre F, Darrier F, Manga L, Akogbeto M, Faye O, Mouchet J, Guillet P: Status of pyrethroid resistance in Anopheles gambiae sensu lato. Bull World Health Organ. 1999, 77: 230-234.PubMed CentralPubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.