Krisztian Magori

Skeeter Buster: a stochastic, spatially explicit modeling tool for studying Aedes aegypti population replacement and population suppression strategies.

Background Dengue is the most important mosquito-borne viral disease affecting humans. The only prevention measure currently available is the control of its vectors, primarily Aedes aegypti. Recent advances in genetic engineering have opened the possibility for a new range of control strategies based on genetically modified mosquitoes. Assessing the potential efficacy of genetic (and conventional) strategies requires the availability of modeling tools that accurately describe the dynamics and genetics of Ae. aegypti populations. Methodology/Principal findings We describe in this paper a new modeling tool of Ae. aegypti population dynamics and genetics named Skeeter Buster. This model operates at the scale of individual water-filled containers for immature stages and individual properties (houses) for adults. The biology of cohorts of mosquitoes is modeled based on the algorithms used in the non-spatial Container Inhabiting Mosquitoes Simulation Model (CIMSiM). Additional features incorporated into Skeeter Buster include stochasticity, spatial structure and detailed population genetics. We observe that the stochastic modeling of individual containers in Skeeter Buster is associated with a strongly reduced temporal variation in stage-specific population densities. We show that heterogeneity in container composition of individual properties has a major impact on spatial heterogeneity in population density between properties. We detail how adult dispersal reduces this spatial heterogeneity. Finally, we present the predicted genetic structure of the population by calculating FST values and isolation by distance patterns, and examine the effects of adult dispersal and container movement between properties. Conclusions/Significance We demonstrate that the incorporated stochasticity and level of spatial detail have major impacts on the simulated population dynamics, which could potentially impact predictions in terms of control measures. The capacity to describe population genetics confers the ability to model the outcome of genetic control methods. Skeeter Buster is therefore an important tool to model Ae. aegypti populations and the outcome of vector control measures.

Regional Differences in the Association Between Land Cover and West Nile Virus Disease Incidence in Humans in the United States

West Nile virus (WNV) is generally considered to be an urban pathogen in the United States, but studies associating land cover and disease incidence, seroprevalence, or infection rate in humans, birds, domesticated and wild mammals, and mosquitoes report varying and sometimes contradictory results at an array of spatial extents. Human infection can provide insight about basic transmission activity; therefore, we analyzed data on the incidence of WNV disease in humans to obtain a comprehensive picture of how human disease and land cover type are associated across the United States. Human WNV disease incidence in Northeastern regions was positively associated with urban land covers, whereas incidence in the Western United States was positively associated with agricultural land covers. We suggest that these regional associations are explained by the geographic distributions of prominent WNV vectors: Culex pipiens complex (including Cx. pipiens and Cx. quinquefasciatus) in the Northeast and Cx. tarsalis in the Western United States.

Spread of white-nose syndrome on a network regulated by geography and climate

Wildlife and plant diseases can reduce biodiversity, disrupt ecosystem services and threaten human health. Emerging pathogens have displayed a variety of spatial spread patterns due to differences in host ecology, including diffusive spread from an epicentre (West Nile virus), jump dispersal on a network (foot-and-mouth disease), or a combination of these (Sudden oak death). White-nose syndrome is a highly pathogenic infectious disease of bats currently spreading across North America. Understanding how bat ecology influences this spread is crucial to management of infected and vulnerable populations. Here we show that white-nose syndrome spread is not diffusive but rather mediated by patchily distributed habitat and large-scale gradients in winter climate. Simulations predict rapid expansion and infection of most counties with caves in the contiguous United States by winter 2105-2106. Our findings show the unique pattern of white-nose syndrome spread corresponds to ecological traits of the host and suggest hypotheses for transmission mechanisms acting at the local scale.

Genetically Engineered Underdominance for Manipulation of Pest Populations: A Deterministic Model

We theoretically investigate the potential for introgressing a desired engineered gene into a pest population by linking the desired gene to DNA constructs that exhibit underdominance properties. Our deterministic model includes two independently segregating engineered constructs that both carry a lethal gene, but suppress each other. Only genotypes containing both or neither construct are viable. Both constructs also carry the desired gene with an independent regulatory mechanism. We examine the minimal number of individuals of an engineered strain that must be released into a natural population to successfully introgress the desired gene. We compare results for strains carrying single and multiple insertions of the constructs. When there are no fitness costs associated with the inserted constructs (when the lethal sequences are not expressed), the number of individuals that must be released decreases as the number of insertions in the genome of the released strain increases. As fitness costs increase, the number of individuals that must be released increases at a greater rate for release strains with more insertions. Under specific conditions this results in the strain with only a single insertion of each construct being the most efficient for introgressing the desired gene. We discuss practical implications of our findings.

INTRODUCING DESIRABLE TRANSGENES INTO INSECT POPULATIONS USING Y-LINKED MEIOTIC DRIVE—A THEORETICAL ASSESSMENT

Abstract The use of genetic drive mechanisms to replace native mosquito genotypes with individuals bearing antipathogen transgenes is a potential strategy for repressing insect transmission of human diseases such as malaria and dengue. Antipathogen transgenes have been developed and tested, but efficient gene drive mechanisms are lacking. Here we theoretically assess the feasibility of introducing antipathogen genes into wild Aedes aegypti populations by using a naturally occurring meiotic drive system. We consider the release of males having both a Y-linked meiotic drive gene and an X-linked drive-insensitive response allele to which an antipathogen gene is linked. We use mathematical models and computer simulations to determine how the post-introduction dynamics of the antipathogen gene are affected by specific genetic characteristics of the system. The results show that when the natural population is uniformly sensitive to the meiotic drive gene, the antipathogen gene may be driven close to fixation if the fitness costs of the drive gene, the insensitive response allele, and the antipathogen gene are low. However, when the natural population has a small proportion of an X-linked insensitive response allele or an autosomal gene that strongly reduces the effect of the drive gene, the antipathogen gene does not spread if it has an associated fitness cost. Our modeling results provide a theoretical foundation for further experimental tests.

Impact of Herbivore-induced Plant Volatiles on Parasitoid Foraging Success: A Spatial Simulation of the Cotesia rubecula, Pieris rapae, and Brassica oleracea System

Many parasitoids are known to use herbivore-induced plant volatiles as cues to locate hosts. However, data are lacking on how much of an advantage a parasitoid can gain from following these plant cues and which factors can limit the value of these cues to the parasitoid. In this study, we simulate the Cotesia rubecula–Pieris rapae–Brassica oleracea system, and ask how many more hosts can a parasitoid attack in a single day of foraging by following plant signals versus randomly foraging. We vary herbivore density, plant response time, parasitoid flight distance, and available host stages to see under which conditions parasitoids benefit from herbivore-induced plant cues. In most of the parameter combinations studied, parasitoids that responded to cues attacked more hosts than those that foraged randomly. Parasitoids following plant cues attacked up to ten times more hosts when they were able to successfully attack herbivores older than first instar; however, if parasitoids were limited to first instar hosts, those following plant cues were at a disadvantage when plants took longer than a day to respond to herbivory. At low herbivore densities, only parasitoids with a larger foraging radius could take advantage of plant cues. Although preference for herbivore-induced volatiles was not always beneficial for a parasitoid, under the most likely natural conditions, the model predicts that C. rubecula gains fitness from following plant cues.

Evaluation of Location-Specific Predictions by a Detailed Simulation Model of Aedes aegypti Populations

Background Skeeter Buster is a stochastic, spatially explicit simulation model of Aedes aegypti populations, designed to predict the outcome of vector population control methods. In this study, we apply the model to two specific locations, the cities of Iquitos, Peru, and Buenos Aires, Argentina. These two sites differ in the amount of field data that is available for location-specific customization. By comparing output from Skeeter Buster to field observations in these two cases we evaluate population dynamics predictions by Skeeter Buster with varying degrees of customization. Methodology/Principal Findings Skeeter Buster was customized to the Iquitos location by simulating the layout of houses and the associated distribution of water-holding containers, based on extensive surveys of Ae. aegypti populations and larval habitats that have been conducted in Iquitos for over 10 years. The model is calibrated by adjusting the food input into various types of containers to match their observed pupal productivity in the field. We contrast the output of this customized model to the data collected from the natural population, comparing pupal numbers and spatial distribution of pupae in the population. Our results show that Skeeter Buster replicates specific population dynamics and spatial structure of Ae. aegypti in Iquitos. We then show how Skeeter Buster can be customized for Buenos Aires, where we only had Ae. aegypti abundance data that was averaged across all locations. In the Argentina case Skeeter Buster provides a satisfactory simulation of temporal population dynamics across seasons. Conclusions This model can provide a faithful description of Ae. aegypti populations, through a process of location-specific customization that is contingent on the amount of data available from field collections. We discuss limitations presented by some specific components of the model such as the description of food dynamics and challenges that these limitations bring to model evaluation.

Multi‐scale model of epidemic fade‐out: Will local extirpation events inhibit the spread of white‐nose syndrome?

White-nose syndrome (WNS) is an emerging infectious disease that has resulted in severe declines of its hibernating bat hosts in North America. The ongoing epidemic of white-nose syndrome is a multi-scale phenomenon becau.se it causes hibernaculum-level extirpations, while simultaneously spreading over larger spatial scales. We investigate a neglected topic in ecological epidemiology: how local pathogen-driven extirpations impact large-scale pathogen spread. Previous studies have identified risk factors for propagation of WNS over hibernaculum and landscape scales but none of these have tested the hypothesis that separation of spatial scales and disease-induced mortality at the hibernaculum level might slow or halt its spread. To test this hypothesis, we developed a mechanistic multi-scale model parameterized using white-nose syndrome.county and site incidence data that connects hibernaculum-level susceptible-infectious-removed (SIR) epidemiology to the county-scale contagion process. Our key result is that hibernaculum-level extirpations will not inhibit county-scale spread of WNS. We show that over 80% of counties of the contiguous USA are likely to become infected before the current epidemic is over and that geometry of habitat connectivity is such that host refuges are exceedingly rare. The macroscale spatiotemporal infection pattern that emerges from local SIR epidemiological processes falls within a narrow spectrum of possible outcomes, suggesting that recolonization, rescue effects, and multi-host complexities at local scales are not important to forward propagation of WNS at large spatial scales. If effective control measures are not implemented, precipitous declines in bat populations are likely, particularly in cave-dense regions that constitute the main geographic corridors of the USA, a serious concern for bat conservation.

Decelerating Spread of West Nile Virus by Percolation in a Heterogeneous Urban Landscape

Vector-borne diseases are emerging and re-emerging in urban environments throughout the world, presenting an increasing challenge to human health and a major obstacle to development. Currently, more than half of the global population is concentrated in urban environments, which are highly heterogeneous in the extent, degree, and distribution of environmental modifications. Because the prevalence of vector-borne pathogens is so closely coupled to the ecologies of vector and host species, this heterogeneity has the potential to significantly alter the dynamical systems through which pathogens propagate, and also thereby affect the epidemiological patterns of disease at multiple spatial scales. One such pattern is the speed of spread. Whereas standard models hold that pathogens spread as waves with constant or increasing speed, we hypothesized that heterogeneity in urban environments would cause decelerating travelling waves in incipient epidemics. To test this hypothesis, we analysed data on the spread of West Nile virus (WNV) in New York City (NYC), the 1999 epicentre of the North American pandemic, during annual epizootics from 2000–2008. These data show evidence of deceleration in all years studied, consistent with our hypothesis. To further explain these patterns, we developed a spatial model for vector-borne disease transmission in a heterogeneous environment. An emergent property of this model is that deceleration occurs only in the vicinity of a critical point. Geostatistical analysis suggests that NYC may be on the edge of this criticality. Together, these analyses provide the first evidence for the endogenous generation of decelerating travelling waves in an emerging infectious disease. Since the reported deceleration results from the heterogeneity of the environment through which the pathogen percolates, our findings suggest that targeting control at key sites could efficiently prevent pathogen spread to remote susceptible areas or even halt epidemics.

When More Transmission Equals Less Disease: Reconciling the Disconnect between Disease Hotspots and Parasite Transmission

The assumed straightforward connection between transmission intensity and disease occurrence impacts surveillance and control efforts along with statistical methodology, including parameter inference and niche modeling. Many infectious disease systems have the potential for this connection to be more complicated–although demonstrating this in any given disease system has remained elusive. Hemorrhagic disease (HD) is one of the most important diseases of white-tailed deer and is caused by viruses in the Orbivirus genus. Like many infectious diseases, the probability or severity of disease increases with age (after loss of maternal antibodies) and the probability of disease is lower upon re-infection compared to first infection (based on cross-immunity between virus strains). These broad criteria generate a prediction that disease occurrence is maximized at intermediate levels of transmission intensity. Using published US field data, we first fit a statistical model to predict disease occurrence as a function of seroprevalence (a proxy for transmission intensity), demonstrating that states with intermediate seroprevalence have the highest level of case reporting. We subsequently introduce an independently parameterized mechanistic model supporting the theory that high case reporting should come from areas with intermediate levels of transmission. This is the first rigorous demonstration of this phenomenon and illustrates that variation in transmission rate (e.g. along an ecologically-controlled transmission gradient) can create cryptic refuges for infectious diseases.