Simulating dynamic insecticide selection pressures for resistance management in mosquitoes assuming polygenic resistance

Insecticide resistance management (IRM) is critical to maintain the operational effectiveness of insecticides used in public health vector control. Evaluating IRM strategies rests primarily on computational models. Most models assume monogenic resistance, but polygenic resistance may be a more appropriate assumption. Conventionally, polygenic models assume selection differentials are constant over successive generations. We present a dynamic method for calculating the selection differentials accounting for the level of resistance and insecticide efficacy. This allows the inclusion of key parameters namely insecticide dosing, insecticide decay and cross resistance, increasing biological and operational realism. Two methods for calculating the insecticide selection differential were compared: truncation (only the most resistant individuals in the population survive) and probabilistic (individual survival depends on their level of resistance). The probabilistic calculation is extendable to multiple gonotrophic cycles, whereby mosquitoes may encounter different insecticides over their life span. A range of IRM strategies of direct policy relevance can be simulated, including the implication of reduced dose mixtures. We describe in detail the calculation and calibration of these models. We demonstrate the ability of the models to simulate a variety of IRM strategies and implications of including these features of the models. In simple IRM strategy evaluations, the truncation and probabilistic models give comparable results to each other and against published polygenic and monogenic models. Analysis of model simulations indicates there is often little difference between sequences or rotations of insecticides. Full-dose mixtures remain the best evaluated IRM strategy. Consistency between models increases confidence in their predictions especially when demonstrating model assumptions do not significantly impact key operational decisions. Using the multiple-gonotrophic cycle model we calculate the age distributions of mosquitoes which provides a framework to link resistance management with disease transmission. Future applications will investigate more scenario-specific evaluations of IRM strategies to inform public health policy.