Search

A promising strategy for reducing the transmission of dengue and other arboviral human diseases by Aedes aegypti mosquito vector populations involves field introductions of the endosymbiotic bacteria Wolbachia. Wolbachia infections inhibit viral transmission by the mosquito, and can spread between mosquito hosts to reach high frequencies in the vector population. Wolbachia spreads by maternal transmission, and spread dynamics can be variable and highly dependent on natural mosquito population dynamics, population structure and fitness components. We develop a mathematical model of an A. aegypti metapopulation that incorporates empirically validated relationships describing density‐dependent mosquito fitness components. We assume that density dependent relationships differ across subpopulations, and construct heterogeneous landscapes for which model‐predicted patterns of variation in mosquito abundance and demography approximate those observed in field populations. We then simulate Wolbachia release strategies similar to that used in field trials. We show that our model can produce rates of spatial spread of Wolbachia similar to those observed following field releases. We then investigate how different types of spatio‐temporal variation in mosquito habitat, as well as different fitness costs incurred by Wolbachia on the mosquito host, influence predicted spread rates. We find that fitness costs reduce spread rates more strongly when the habitat landscape varies temporally due to stochastic and seasonal processes. Synthesis and applications: Our empirically based modelling approach represents effects of environmental heterogeneity on the spatial spread of Wolbachia. The models can assist in interpreting observed spread patterns following field releases and in designing suitable release strategies for targeting spatially heterogeneous vector populations.

BACKGROUND: Arbovirus transmission by the mosquito Aedes aegypti can be reduced by the introduction and establishment of the endosymbiotic bacteria Wolbachia in wild populations of the vector. Wolbachia spreads by increasing the fitness of its hosts relative to uninfected mosquitoes. However, mosquito fitness is also strongly affected by population size through density-dependent competition for limited food resources. We do not understand how this natural variation in fitness affects symbiont spread, which limits our ability to design successful control strategies. RESULTS: We develop a mathematical model to predict A. aegypti-Wolbachia dynamics that incorporates larval density-dependent variation in important fitness components of infected and uninfected mosquitoes. Our model explains detailed features of the mosquito-Wolbachia dynamics observed in two independent experimental A. aegypti populations, allowing the combined effects on dynamics of multiple density-dependent fitness components to be characterized. We apply our model to investigate Wolbachia field release dynamics, and show how invasion outcomes can depend strongly on the severity of density-dependent competition at the release site. Specifically, the ratio of released relative to wild mosquitoes required to attain a target infection frequency (at the end of a release program) can vary by nearly an order of magnitude. The time taken for Wolbachia to become established following releases can differ by over 2 years. These effects depend on the relative fitness of field and insectary-reared mosquitoes. CONCLUSIONS: Models of Wolbachia invasion incorporating density-dependent demographic variation in the host population explain observed dynamics in experimental A. aegypti populations. These models predict strong effects of density-dependence on Wolbachia dynamics in field populations, and can assist in the effective use of Wolbachia to control the transmission of arboviruses such as dengue, chikungunya and