Incorporating genetic selection into individual-based models of malaria and other infectious diseases

Control strategies for human infections are often investigated using individual-based models (IBMs) to quantify their impact in terms of mortality, morbidity and impact on transmission. Genetic selection can be incorporated into the IBMs to track the spread of mutations whose origin and spread are driven by the intervention and which subsequently undermine the control strategy; typical examples are mutations which encode drug resistance or diagnosis- or vaccine-escape phenotypes.

We simulated the spread of malaria drug resistance using the IBM OpenMalaria to investigate how the finite sizes of IBMs require strategies to optimally incorporate genetic selection. We make four recommendations. Firstly, calculate and report the selection coefficients, s, of the advantageous allele as the key genetic parameter. Secondly, use these values of "s" to calculate the wait time until a mutation successfully establishes itself in the pathogen population. Thirdly, identify the inherent limits of the IBM to robustly estimate small selection coefficients. Fourthly, optimize computational efficacy: when "s" is small, fewer replicates of larger IBMs may be more efficient than a larger number of replicates of smaller size.

The OpenMalaria IBM of malaria was an exemplar and the same principles apply to IBMs of other diseases.