Modelling the spread of Wolbachia in spatially heterogeneous environments

The endosymbiont Wolbachia infects a large number of insect species and is capable of rapid spread when introduced into a novel host population. The bacteria spread by manipulating their hosts' reproduction, and their dynamics are influenced by the demographic structure of the host population and patterns of contact between individuals. Reaction–diffusion models of the spatial spread of Wolbachia provide a simple analytical description of their spatial dynamics but do not account for significant details of host population dynamics. We develop a metapopulation model describing the spatial dynamics of Wolbachia in an age-structured host insect population regulated by juvenile density-dependent competition. The model produces similar dynamics to the reaction–diffusion model in the limiting case where the host's habitat quality is spatially homogeneous and Wolbachia has a small effect on host fitness. When habitat quality varies spatially, Wolbachia spread is usually much slower, and the conditions necessary for local invasion are strongly affected by immigration of insects from surrounding regions. Spread is most difficult when variation in habitat quality is spatially correlated. The results show that spatial variation in the density-dependent competition experienced by juvenile host insects can strongly affect the spread of Wolbachia infections, which is important to the use of Wolbachia to control insect vectors of human disease and other pests.

[1]  J. Werren Biology of Wolbachia. , 2003, Annual review of entomology.

[2]  P. Hammerstein,et al.  Wolbachia‐induced unidirectional cytoplasmic incompatibility and the stability of infection polymorphism in parapatric host populations , 2007, Journal of evolutionary biology.

[3]  S. Ritchie,et al.  The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations , 2011, Nature.

[4]  V. S. Nam,et al.  Characterizing the Aedes aegypti Population in a Vietnamese Village in Preparation for a Wolbachia-Based Mosquito Control Strategy to Eliminate Dengue , 2009, PLoS neglected tropical diseases.

[5]  T. Scott,et al.  Skeeter Buster: A Stochastic, Spatially Explicit Modeling Tool for Studying Aedes aegypti Population Replacement and Population Suppression Strategies , 2009, PLoS neglected tropical diseases.

[6]  H. Laven,et al.  Cytoplasmic Inheritance in Culex , 1956, Nature.

[7]  A. Lloyd,et al.  Density-Dependent Intraspecific Competition in the Larval Stage of Aedes aegypti (Diptera: Culicidae): Revisiting the Current Paradigm , 2009, Journal of medical entomology.

[8]  C. Codeço,et al.  Temporal Distribution of Aedes aegypti in Different Districts of Rio De Janeiro, Brazil, Measured by Two Types of Traps , 2009, Journal of medical entomology.

[9]  A. R. Barr,et al.  New Hypothesis of the Cause of Cytoplasmic Incompatibility in Culex pipiens L. , 1971, Nature.

[10]  N. Barton,et al.  Spatial Waves of Advance with Bistable Dynamics: Cytoplasmic and Genetic Analogues of Allee Effects , 2011, The American Naturalist.

[11]  T. Scott,et al.  Influence of container size, location, and time of day on oviposition patterns of the dengue vector, Aedes aegypti, in Thailand. , 2008, Vector borne and zoonotic diseases.

[12]  N. Barton,et al.  The dynamics of hybrid zones , 1979, Heredity.

[13]  S. Levin,et al.  Diffusion and Ecological Problems: Modern Perspectives , 2013 .

[14]  Christopher Dye,et al.  Models for the population dynamics of the yellow fever mosquito, Aedes aegypti , 1984 .

[15]  Peter Schlattmann,et al.  How many species are infected with Wolbachia? – a statistical analysis of current data , 2008, FEMS microbiology letters.

[16]  M. Turelli EVOLUTION OF INCOMPATIBILITY‐INDUCING MICROBES AND THEIR HOSTS , 1994, Evolution; international journal of organic evolution.

[17]  Masanobu Shinozuka,et al.  Simulation of Stochastic Fields by Statistical Preconditioning , 1990 .

[18]  Bodil N. Cass,et al.  Stable Introduction of a Life-Shortening Wolbachia Infection into the Mosquito Aedes aegypti , 2009, Science.

[19]  P. M. Sheppard,et al.  The dynamics of an adult population of Aedes aegypti in relation to dengue haemorrhagic fever in Bangkok. , 1969 .

[20]  Zhiyong Xi,et al.  The Endosymbiotic Bacterium Wolbachia Induces Resistance to Dengue Virus in Aedes aegypti , 2010, PLoS pathogens.

[21]  Peter A. Ryan,et al.  A Wolbachia Symbiont in Aedes aegypti Limits Infection with Dengue, Chikungunya, and Plasmodium , 2009, Cell.

[22]  K. Ruegg Genetic, Morphological, and Ecological Characterization of a Hybrid Zone That Spans a Migratory Divide , 2008, Evolution; international journal of organic evolution.

[23]  V. Jansen,et al.  Stochastic spread of Wolbachia , 2008, Proceedings of the Royal Society B: Biological Sciences.

[24]  M. Turelli,et al.  Rapid spread of an inherited incompatibility factor in California Drosophila , 1991, Nature.

[25]  N. Moran,et al.  Variation in resistance to parasitism in aphids is due to symbionts not host genotype. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[26]  A. Ellis,et al.  Linking movement and oviposition behaviour to spatial population distribution in the tree hole mosquito Ochlerotatus triseriatus. , 2008, The Journal of animal ecology.

[27]  S. Sinkins,et al.  Immune Activation by Life-Shortening Wolbachia and Reduced Filarial Competence in Mosquitoes , 2009, Science.

[28]  T. Scott,et al.  Impact of population age structure on Wolbachia transgene driver efficacy: ecologically complex factors and release of genetically modified mosquitoes. , 2004, Insect biochemistry and molecular biology.

[29]  N. Barton,et al.  Adaptation, speciation and hybrid zones , 1989, Nature.

[30]  Serena Matucci,et al.  AGGREGATIVE MOVEMENT AND FRONT PROPAGATION FOR BI-STABLE POPULATION MODELS , 2007 .

[31]  S. Ritchie,et al.  Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission , 2011, Nature.

[32]  Arthur Getis,et al.  Spatial and temporal patterns in pupal and adult production of the dengue vector Aedes aegypti in Kamphaeng Phet, Thailand. , 2008, The American journal of tropical medicine and hygiene.

[33]  M. Kot,et al.  Integrodifference equations, Allee effects, and invasions , 2002, Journal of mathematical biology.

[34]  Karyn N. Johnson,et al.  Solving the Wolbachia Paradox: Modeling the Tripartite Interaction between Host, Wolbachia, and a Natural Enemy , 2011, The American Naturalist.

[35]  L. Lehmann,et al.  The spread of incompatibility-inducing parasites in sub-divided host populations , 2008, BMC Evolutionary Biology.

[36]  Craig R Williams,et al.  The development of predictive tools for pre‐emptive dengue vector control: a study of Aedes aegypti abundance and meteorological variables in North Queensland, Australia , 2010, Tropical medicine & international health : TM & IH.

[37]  Penelope A. Hancock,et al.  Population Dynamic Models of the Spread of Wolbachia , 2011, The American Naturalist.

[38]  T. Scott,et al.  Longitudinal Studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: Population Dynamics , 2000, Journal of medical entomology.

[39]  T. Scott,et al.  Oviposition Site Selection by the Dengue Vector Aedes aegypti and Its Implications for Dengue Control , 2011, PLoS neglected tropical diseases.

[40]  Karyn N. Johnson,et al.  Wolbachia and Virus Protection in Insects , 2008, Science.

[41]  Françoise Lescourret,et al.  Factors and mechanisms explaining spatial heterogeneity: a review of methods for insect populations , 2011 .

[42]  M. Kimura,et al.  An introduction to population genetics theory , 1971 .

[43]  N. Shigesada,et al.  Biological Invasions: Theory and Practice , 1997 .

[44]  A. Himler,et al.  Rapid Spread of a Bacterial Symbiont in an Invasive Whitefly Is Driven by Fitness Benefits and Female Bias , 2011, Science.

[45]  S. Ritchie,et al.  Mark–release–recapture study to measure dispersal of the mosquito Aedes aegypti in Cairns, Queensland, Australia , 2005, Medical and veterinary entomology.

[46]  P. Schofield Spatially explicit models of Turelli-Hoffmann Wolbachia invasive wave fronts. , 2002, Journal of theoretical biology.

[47]  R. J. Gutiérrez,et al.  Genetic structure, introgression, and a narrow hybrid zone between northern and California spotted owls (Strix occidentalis) , 2005, Molecular ecology.

[48]  Marcus S. C. Blagrove,et al.  Wolbachia Stimulates Immune Gene Expression and Inhibits Plasmodium Development in Anopheles gambiae , 2010, PLoS pathogens.

[49]  T. Nagylaki,et al.  Conditions for the existence of clines. , 1975, Genetics.

[50]  Penelope A. Hancock,et al.  Strategies for Introducing Wolbachia to Reduce Transmission of Mosquito-Borne Diseases , 2011, PLoS neglected tropical diseases.