Search

Transmission heterogeneity, whereby a disproportionate fraction of pathogen transmission events result from a small number of individuals or geographic locations, is an inherent property of many, if not most, infectious disease systems. For vector-borne diseases, transmission heterogeneity is inferred from the distribution of the number of vectors per host, which could lead to significant bias in situations where vector abundance and transmission risk at the household do not correlate, as is the case with dengue virus (DENV). We used data from a contact tracing study to quantify the distribution of DENV acute infections within human activity spaces (AS), the collection of residential locations an individual routinely visits, and quantified measures of virus transmission heterogeneity from two consecutive dengue outbreaks (DENV-4 and DENV-2) that occurred in the city of Iquitos, Peru. Negative-binomial distributions and Pareto fractions showed evidence of strong overdispersion in the number of DENV infections by AS and identified super-spreading units (SSUs): i.e. AS where most infections occurred. Approximately 8% of AS were identified as SSUs, contributing to more than 50% of DENV infections. SSU occurrence was associated more with DENV-2 infection than with DENV-4, a predominance of inapparent infections (74% of all infections), households with high Aedes aegypti mosquito abundance, and high host susceptibility to the circulating DENV serotype. Marked heterogeneity in dengue case distribution, and the role of inapparent infections in defining it, highlight major challenges faced by reactive interventions if those transmission units contributing the most to transmission are not identified, prioritized, and effectively treated.

Dengue virus (DENV) transmission from humans to mosquitoes is a poorly documented, but critical component of DENV epidemiology. Magnitude of viremia is the primary determinant of successful human-to-mosquito DENV transmission. People with the same level of viremia, however, can vary in their infectiousness to mosquitoes as a function of other factors that remain to be elucidated. Here, we report on a field-based study in the city of Iquitos, Peru, where we conducted direct mosquito feedings on people naturally infected with DENV and that experienced mild illness. We also enrolled people naturally infected with Zika virus (ZIKV) after the introduction of ZIKV in Iquitos during the study period. Of the 54 study participants involved in direct mosquito feedings, 43 were infected with DENV-2, two with DENV-3, and nine with ZIKV. Our analysis excluded participants whose viremia was detectable at enrollment but undetectable at the time of mosquito feeding, which was the case for all participants with DENV-3 and ZIKV infections. We analyzed the probability of onward transmission during 50 feeding events involving 27 participants infected with DENV-2 based on the presence of infectious virus in mosquito saliva 7-16 days post blood meal. Transmission probability was positively associated with the level of viremia and duration of extrinsic incubation in the mosquito. In addition, transmission probability was influenced by the day of illness in a non-monotonic fashion; i.e., transmission probability increased until 2 days after symptom onset and decreased thereafter. We conclude that mildly ill DENV-infected humans with similar levels of viremia during the first two days after symptom onset will be most infectious to mosquitoes on the second day of their illness. Quantifying variation within and between people in their contribution to DENV transmission is essential to better understand the biological determinants of human infectiousness, parametrize epidemiological models, and improve disease surveillance and prevention strategies.

The mosquito Aedes aegypti is the vector of a number of medically-important viruses, including dengue virus, yellow fever virus, chikungunya virus, and Zika virus, and as such vector control is a key approach to managing the diseases they cause. Understanding the impact of vector control on these diseases is aided by first understanding its impact on Ae. aegypti population dynamics. A number of detail-rich models have been developed to couple the dynamics of the immature and adult stages of Ae. aegypti. The numerous assumptions of these models enable them to realistically characterize impacts of mosquito control, but they also constrain the ability of such models to reproduce empirical patterns that do not conform to the models' behavior. In contrast, statistical models afford sufficient flexibility to extract nuanced signals from noisy data, yet they have limited ability to make predictions about impacts of mosquito control on disease caused by pathogens that the mosquitoes transmit without extensive data on mosquitoes and disease. Here, we demonstrate how the differing strengths of mechanistic realism and statistical flexibility can be fused into a single model. Our analysis utilizes data from 176,352 household-level Ae. aegypti aspirator collections conducted during 1999-2011 in Iquitos, Peru. The key step in our approach is to calibrate a single parameter of the model to spatio-temporal abundance patterns predicted by a generalized additive model (GAM). In effect, this calibrated parameter absorbs residual variation in the abundance time-series not captured by other features of the mechanistic model. We then used this calibrated parameter and the literature-derived parameters in the agent-based model to explore Ae. aegypti population dynamics and the impact of insecticide spraying to kill adult mosquitoes. The baseline abundance predicted by the agent-based model closely matched that predicted by the GAM. Following spraying, the agent-based model predicted that mosquito abundance rebounds within about two months, commensurate with recent experimental data from Iquitos. Our approach was able to accurately reproduce abundance patterns in Iquitos and produce a realistic response to adulticide spraying, while retaining sufficient flexibility to be applied across a range of settings.

Current knowledge of dengue virus (DENV) transmission provides only a partial understanding of a complex and dynamic system yielding a public health track record that has more failures than successes. An important part of the problem is that the foundation for contemporary interventions includes a series of longstanding, but untested, assumptions based on a relatively small portion of the human population; i.e., people who are convenient to study because they manifest clinically apparent disease. Approaching dengue from the perspective of people with overt illness has produced an extensive body of useful literature. It has not, however, fully embraced heterogeneities in virus transmission dynamics that are increasingly recognized as key information still missing in the struggle to control the most important insect-transmitted viral infection of humans. Only in the last 20 years have there been significant efforts to carry out comprehensive longitudinal dengue studies. This manuscript provides the rationale and comprehensive, integrated description of the methodology for a five-year longitudinal cohort study based in the tropical city of Iquitos, in the heart of the Peruvian Amazon. Primary data collection for this study was completed in 2019. Although some manuscripts have been published to date, our principal objective here is to support subsequent publications by describing in detail the structure, methodology, and significance of a specific research program. Our project was designed to study people across the entire continuum of disease, with the ultimate goal of quantifying heterogeneities in human variables that affect DENV transmission dynamics and prevention. Because our study design is applicable to other Aedes transmitted viruses, we used it to gain insights into Zika virus (ZIKV) transmission when during the project period ZIKV was introduced and circulated in Iquitos. Our prospective contact cluster investigation design was initiated by detecttion of a person with a symptomatic DENV infection and then followed that person's immediate contacts. This allowed us to monitor individuals at high risk of DENV infection, including people with clinically inapparent and mild infections that are otherwise difficult to detect. We aimed to fill knowledge gaps by defining the contribution to DENV transmission dynamics of (1) the understudied majority of DENV-infected people with inapparent and mild infections and (2) epidemiological, entomological, and socio-behavioral sources of heterogeneity. By accounting for factors underlying variation in each person's contribution to transmission we sought to better determine the type and extent of effort needed to better prevent virus transmission and disease.

The importance of mosquitoes in human pathogen transmission has motivated major research efforts into mosquito biology in pursuit of more effective vector control measures. Aedes aegypti is a particular concern in tropical urban areas, where it is the primary vector of numerous flaviviruses, including the yellow fever, Zika, and dengue viruses. With an anthropophilic habit, Ae. aegypti prefers houses, human blood meals, and ovipositioning in water-filled containers. We hypothesized that this relatively simple ecological niche should allow us to predict the impacts of insecticidal control measures on mosquito populations. To do this, we use Skeeter Buster 2 (SB2), a stochastic, spatially explicit, mechanistic model of Ae. aegypti population biology. SB2 builds on Skeeter Buster, which reproduced equilibrium dynamics of Ae. aegypti in Iquitos, Peru. Our goal was to validate SB2 by predicting the response of mosquito populations to perturbations by indoor insecticidal spraying and widespread destructive insect surveys. To evaluate SB2, we conducted two field experiments in Iquitos, Peru: a smaller pilot study in 2013 (S-2013) followed by a larger experiment in 2014 (L-2014). Here, we compare model predictions with (previously reported) empirical results from these experiments. In both simulated and empirical populations, repeated spraying yielded substantial yet temporary reductions in adult densities. The proportional effects of spraying were broadly comparable between simulated and empirical results, but we found noteworthy differences. In particular, SB2 consistently over-estimated the proportion of nulliparous females and the proportion of containers holding immature mosquitoes. We also observed less temporal variation in simulated surveys of adult abundance relative to corresponding empirical observations. Our results indicate the presence of ecological heterogeneities or sampling processes not effectively represented by SB2. Although additional empirical research could further improve the accuracy and precision of SB2, our results underscore the importance of non-linear dynamics in the response of Ae. aegypti populations to perturbations, and suggest general limits to the fine-grained predictability of its population dynamics over space and time.

This study describes the evolution of knockdown resistance (kdr) haplotypes in Aedes aegypti in response to pyrethroid insecticide use over the course of 18 years in Iquitos, Peru. Based on the duration and intensiveness of sampling (~10,000 samples), this is the most thorough study of kdr population genetics in Ae. aegypti to date within a city. We provide evidence for the direct connection between programmatic citywide pyrethroid spraying and the increase in frequency of specific kdr haplotypes by identifying two evolutionary events in the population. The relatively high selection coefficients, even under infrequent insecticide pressure, emphasize how quickly Ae. aegypti populations can evolve. In our examination of the literature on mosquitoes and other insect pests, we could find no cases where a pest evolved so quickly to so few exposures to low or nonresidual insecticide applications. The observed rapid increase in frequency of resistance alleles might have been aided by the incomplete dominance of resistance-conferring alleles over corresponding susceptibility alleles. In addition to dramatic temporal shifts, spatial suppression experiments reveal that genetic heterogeneity existed not only at the citywide scale, but also on a very fine scale within the city.

Heterogeneous exposure to mosquitoes determines an individual's contribution to vector-borne pathogen transmission. Particularly for dengue virus (DENV), there is a major difficulty in quantifying human-vector contacts due to the unknown coupled effect of key heterogeneities. To test the hypothesis that the reduction of human out-of-home mobility due to dengue illness will significantly influence population-level dynamics and the structure of DENV transmission chains, we extended an existing modeling framework to include social structure, disease-driven mobility reductions, and heterogeneous transmissibility from different infectious groups. Compared to a baseline model, naïve to human pre-symptomatic infectiousness and disease-driven mobility changes, a model including both parameters predicted an increase of 37% in the probability of a DENV outbreak occurring; a model including mobility change alone predicted a 15.5% increase compared to the baseline model. At the individual level, models including mobility change led to a reduction of the importance of out-of-home onward transmission (R, the fraction of secondary cases predicted to be generated by an individual) by symptomatic individuals (up to -62%) at the expense of an increase in the relevance of their home (up to +40%). An individual's positive contribution to R could be predicted by a GAM including a non-linear interaction between an individual's biting suitability and the number of mosquitoes in their home (>10 mosquitoes and 0.6 individual attractiveness significantly increased R). We conclude that the complex fabric of social relationships and differential behavioral response to dengue illness cause the fraction of symptomatic DENV infections to concentrate transmission in specific locations, whereas asymptomatic carriers (including individuals in their pre-symptomatic period) move the virus throughout the landscape. Our findings point to the difficulty of focusing vector control interventions reactively on the home of symptomatic individuals, as this approach will fail to contain virus propagation by visitors to their house and asymptomatic carriers.

Recent years have seen rising incidence of dengue and large outbreaks of Zika and chikungunya, which are all caused by viruses transmitted by Aedes aegypti mosquitoes. In most settings, the primary intervention against Aedes-transmitted viruses is vector control, such as indoor, ultra-low volume (ULV) spraying. Targeted indoor residual spraying (TIRS) has the potential to more effectively impact Aedes-borne diseases, but its implementation requires careful planning and evaluation. The optimal time to deploy these interventions and their relative epidemiological effects are, however, not well understood. We used an agent-based model of dengue virus transmission calibrated to data from Iquitos, Peru to assess the epidemiological effects of these interventions under differing strategies for deploying them. Specifically, we compared strategies where spray application was initiated when incidence rose above a threshold based on incidence in recent years to strategies where spraying occurred at the same time(s) each year. In the absence of spraying, the model predicted 361,000 infections [inter-quartile range (IQR): 347,000-383,000] in the period 2000-2010. The ULV strategy with the fewest median infections was spraying twice yearly, in March and October, which led to a median of 172,000 infections [IQR: 158,000-183,000], a 52% reduction from baseline. Compared to spraying once yearly in September, the best threshold-based strategy utilizing ULV had fewer median infections (254,000 vs. 261,000), but required more spraying (351 vs. 274 days). For TIRS, the best strategy was threshold-based, which led to the fewest infections of all strategies tested (9,900; [IQR: 8,720-11,400], a 94% reduction), and required fewer days spraying than the equivalent ULV strategy (280). Although spraying twice each year is likely to avert the most infections, our results indicate that a threshold-based strategy can become an alternative to better balance the translation of spraying effort into impact, particularly if used with a residual insecticide.

Despite estimates that, each year, as many as 300 million dengue virus (DENV) infections result in either no perceptible symptoms (asymptomatic) or symptoms that are sufficiently mild to go undetected by surveillance systems (inapparent), it has been assumed that these infections contribute little to onward transmission. However, recent blood-feeding experiments with Aedes aegypti mosquitoes showed that people with asymptomatic and pre-symptomatic DENV infections are capable of infecting mosquitoes. To place those findings into context, we used models of within-host viral dynamics and human demographic projections to (1) quantify the net infectiousness of individuals across the spectrum of DENV infection severity and (2) estimate the fraction of transmission attributable to people with different severities of disease. Our results indicate that net infectiousness of people with asymptomatic infections is 80% (median) that of people with apparent or inapparent symptomatic infections (95% credible interval (CI): 0-146%). Due to their numerical prominence in the infectious reservoir, clinically inapparent infections in total could account for 84% (CI: 82-86%) of DENV transmission. Of infections that ultimately result in any level of symptoms, we estimate that 24% (95% CI: 0-79%) of onward transmission results from mosquitoes biting individuals during the pre-symptomatic phase of their infection. Only 1% (95% CI: 0.8-1.1%) of DENV transmission is attributable to people with clinically detected infections after they have developed symptoms. These findings emphasize the need to (1) reorient current practices for outbreak response to adoption of pre-emptive strategies that account for contributions of undetected infections and (2) apply methodologies that account for undetected infections in surveillance programs, when assessing intervention impact, and when modeling mosquito-borne virus transmission.

BackgroundAedes aegypti is a primary vector of dengue, chikungunya, Zika, and urban yellow fever viruses. Indoor, ultra low volume (ULV) space spraying with pyrethroid insecticides is the main approach used for Ae. aegypti emergency control in many countries. Given the widespread use of this method, the lack of large-scale experiments or detailed evaluations of municipal spray programs is problematic.Methodology/principal findingsTwo experimental evaluations of non-residual, indoor ULV pyrethroid spraying were conducted in Iquitos, Peru. In each, a central sprayed sector was surrounded by an unsprayed buffer sector. In 2013, spray and buffer sectors included 398 and 765 houses, respectively. Spraying reduced the mean number of adults captured per house by ~83 percent relative to the pre-spray baseline survey. In the 2014 experiment, sprayed and buffer sectors included 1,117 and 1,049 houses, respectively. Here, the sprayed sector's number of adults per house was reduced ~64 percent relative to baseline. Parity surveys in the sprayed sector during the 2014 spray period indicated an increase in the proportion of very young females. We also evaluated impacts of a 2014 citywide spray program by the local Ministry of Health, which reduced adult populations by ~60 percent. In all cases, adult densities returned to near-baseline levels within one month.Conclusions/significanceOur results demonstrate that densities of adult Ae. aegypti can be reduced by experimental and municipal spraying programs. The finding that adult densities return to approximately pre-spray densities in less than a month is similar to results from previous, smaller scale experiments. Our results demonstrate that ULV spraying is best viewed as having a short-term entomological effect. The epidemiological impact of ULV spraying will need evaluation in future trials that measure capacity of insecticide spraying to reduce human infection or disease.

Despite estimates that, each year, as many as 300 million dengue virus (DENV) infections result in either no perceptible symptoms (asymptomatic) or symptoms that are sufficiently mild to go undetected by surveillance systems (inapparent), it has been assumed that these infections contribute little to onward transmission. However, recent blood-feeding experiments with Aedes aegypti mosquitoes showed that people with asymptomatic and pre-symptomatic DENV infections are capable of infecting mosquitoes. To place those findings into context, we used models of within-host viral dynamics and human demographic projections to (1) quantify the net infectiousness of individuals across the spectrum of DENV infection severity and (2) estimate the fraction of transmission attributable to people with different severities of disease. Our results indicate that net infectiousness of people with asymptomatic infections is 80% (median) that of people with apparent or inapparent symptomatic infections (95% credible interval (CI): 0-146%). Due to their numerical prominence in the infectious reservoir, clinically inapparent infections in total could account for 84% (CI: 82-86%) of DENV transmission. Of infections that ultimately result in any level of symptoms, we estimate that 24% (95% CI: 0-79%) of onward transmission results from mosquitoes biting individuals during the pre-symptomatic phase of their infection. Only 1% (95% CI: 0.8-1.1%) of DENV transmission is attributable to people with clinically detected infections after they have developed symptoms. These findings emphasize the need to (1) reorient current practices for outbreak response to adoption of pre-emptive strategies that account for contributions of undetected infections and (2) apply methodologies that account for undetected infections in surveillance programs, when assessing intervention impact, and when modeling mosquito-borne virus transmission.

Most host-parasite systems exhibit remarkable heterogeneity in the contribution to transmission of certain individuals, locations, host infectious states, or parasite strains. While significant advancements have been made in the understanding of the impact of transmission heterogeneity in epidemic dynamics and parasite persistence and evolution, the knowledge base of the factors contributing to transmission heterogeneity is limited. We argue that research efforts should move beyond considering the impact of single sources of heterogeneity and account for complex couplings between conditions with potential synergistic impacts on parasite transmission. Using theoretical approaches and empirical evidence from various host-parasite systems, we investigate the ecological and epidemiological significance of couplings between heterogeneities and discuss their potential role in transmission dynamics and the impact of control.

BackgroundPathogen transmission by mosquitos is known to be highly sensitive to mosquito bionomic parameters. Mosquito mark-release-recapture (MMRR) experiments are a standard method for estimating such parameters including dispersal, population size and density, survival, blood feeding frequency and blood meal host preferences.MethodsWe assembled a comprehensive database describing adult female MMRR experiments. Bibliographic searches were used to build a digital library of MMRR studies and selected data describing the reported outcomes were extracted.ResultsThe resulting database contained 774 unique adult female MMRR experiments involving 58 vector mosquito species from the three main genera of importance to human health: Aedes, Anopheles and Culex. Crude examination of these data revealed patterns associated with geography as well as mosquito genus, consistent with bionomics varying by species-specific life history and ecological context. Recapture success varied considerably and was significantly different amongst genera, with 8, 4 and 1% of adult females recaptured for Aedes, Anopheles and Culex species, respectively. A large proportion of experiments (59%) investigated dispersal and survival and many allowed disaggregation of the release and recapture data. Geographic coverage was limited to just 143 localities around the world.ConclusionsThis MMRR database is a substantial contribution to the compilation of global data that can be used to better inform basic research and public health interventions, to identify and fill knowledge gaps and to enrich theory and evidence-based ecological and epidemiological studies of mosquito vectors, pathogen transmission and disease prevention. The database revealed limited geographic coverage and a relative scarcity of information for vector species of substantial public health relevance. It represents, however, a wealth of entomological information not previously compiled and of particular interest for mosquito-borne pathogen transmission models.

Infectious disease models play a key role in public health planning. These models rely on accurate estimates of key transmission parameters such as the force of infection (FoI), which is the per-capita risk of a susceptible person being infected. The FoI captures the fundamental dynamics of transmission and is crucial for gauging control efforts, such as identifying vaccination targets. Dengue virus (DENV) is a mosquito-borne, multiserotype pathogen that currently infects ∼390 million people a year. Existing estimates of the DENV FoI are inaccurate because they rely on the unrealistic assumption that risk is constant over time. Dengue models are thus unreliable for designing vaccine deployment strategies. Here, we present to our knowledge the first time-varying (daily), serotype-specific estimates of DENV FoIs using a spline-based fitting procedure designed to examine a 12-y, longitudinal DENV serological dataset from Iquitos, Peru (11,703 individuals, 38,416 samples, and 22,301 serotype-specific DENV infections from 1999 to 2010). The yearly DENV FoI varied markedly across time and serotypes (0-0.33), as did daily basic reproductive numbers (0.49-4.72). During specific time periods, the FoI fluctuations correlated across serotypes, indicating that different DENV serotypes shared common transmission drivers. The marked variation in transmission intensity that we detected indicates that intervention targets based on one-time estimates of the FoI could underestimate the level of effort needed to prevent disease. Our description of dengue virus transmission dynamics is unprecedented in detail, providing a basis for understanding the persistence of this rapidly emerging pathogen and improving disease prevention programs.

Mosquito-borne diseases pose some of the greatest challenges in public health, especially in tropical and sub-tropical regions of the world. Efforts to control these diseases have been underpinned by a theoretical framework developed for malaria by Ross and Macdonald, including models, metrics for measuring transmission, and theory of control that identifies key vulnerabilities in the transmission cycle. That framework, especially Macdonald's formula for R0 and its entomological derivative, vectorial capacity, are now used to study dynamics and design interventions for many mosquito-borne diseases. A systematic review of 388 models published between 1970 and 2010 found that the vast majority adopted the Ross-Macdonald assumption of homogeneous transmission in a well-mixed population. Studies comparing models and data question these assumptions and point to the capacity to model heterogeneous, focal transmission as the most important but relatively unexplored component in current theory. Fine-scale heterogeneity causes transmission dynamics to be nonlinear, and poses problems for modeling, epidemiology and measurement. Novel mathematical approaches show how heterogeneity arises from the biology and the landscape on which the processes of mosquito biting and pathogen transmission unfold. Emerging theory focuses attention on the ecological and social context for mosquito blood feeding, the movement of both hosts and mosquitoes, and the relevant spatial scales for measuring transmission and for modeling dynamics and control.

Objective: Surveillance is critical for the rapid implementation of control measures for diseases caused by aerially dispersed plant pathogens, but such programs can be resource-intensive, especially for epidemics caused by long-distance dispersed pathogens. The current cucurbit downy mildew platform for monitoring, predicting and communicating the risk of disease spread in the United States is expensive to maintain. In this study, we focused on identifying sites critical for surveillance and treatment in an attempt to reduce disease monitoring costs and determine where control may be applied to mitigate the risk of disease spread. Methods: Static networks were constructed based on the distance between fields, while dynamic networks were constructed based on the distance between fields and wind speed and direction, using disease data collected from 2008 to 2016. Three strategies were used to identify highly connected field sites. First, the probability of pathogen transmission between nodes and the probability of node infection were modeled over a discrete weekly time step within an epidemic year. Second, nodes identified as important were selectively removed from networks and the probability of node infection was recalculated in each epidemic year. Third, the recurring patterns of node infection were analyzed across epidemic years. Results: Static networks exhibited scale-free properties where the node degree followed a power-law distribution. Betweenness centrality was the most useful metric for identifying important nodes within the networks that were associated with disease transmission and prediction. Based on betweenness centrality, field sites in Maryland, North Carolina, Ohio, South Carolina and Virginia were the most central in the disease network across epidemic years. Removing field sites identified as important limited the predicted risk of disease spread based on the dynamic network model. Conclusions: Combining the dynamic network model and centrality metrics facilitated the identification of highly connected fields in the southeastern United States and the mid-Atlantic region. These highly connected sites may be used to inform surveillance and strategies for controlling cucurbit downy mildew in the eastern United States. [ABSTRACT FROM AUTHOR]

Mathematical models of mosquito-borne pathogen transmission originated in the early twentieth century to provide insights into how to most effectively combat malaria. The foundations of the Ross-Macdonald theory were established by 1970. Since then, there has been a growing interest in reducing the public health burden of mosquito-borne pathogens and an expanding use of models to guide their control. To assess how theory has changed to confront evolving public health challenges, we compiled a bibliography of 325 publications from 1970 through 2010 that included at least one mathematical model of mosquito-borne pathogen transmission and then used a 79-part questionnaire to classify each of 388 associated models according to its biological assumptions. As a composite measure to interpret the multidimensional results of our survey, we assigned a numerical value to each model that measured its similarity to 15 core assumptions of the Ross-Macdonald model. Although the analysis illustrated a growing acknowledgement of geographical, ecological and epidemiological complexities in modelling transmission, most models during the past 40 years closely resemble the Ross-Macdonald model. Modern theory would benefit from an expansion around the concepts of heterogeneous mosquito biting, poorly mixed mosquito-host encounters, spatial heterogeneity and temporal variation in the transmission process.

Recently there have been significant advances in research on genetic strategies to control populations of disease-vectoring insects. Some of these strategies use the gene drive properties of selfish genetic elements to spread physically linked anti-pathogen genes into local vector populations. Because of the potential of these selfish elements to spread through populations, control approaches based on these strategies must be carefully evaluated to ensure a balance between the desirable spread of the refractoriness-conferring genetic cargo and the avoidance of potentially unwanted outcomes such as spread to non-target populations. There is also a need to develop better estimates of the economics of such releases. We present here an evaluation of two such strategies using a biologically realistic mathematical model that simulates the resident Aedes aegypti mosquito population of Iquitos, Peru. One strategy uses the selfish element Medea, a non-limited element that could permanently spread over a large geographic area; the other strategy relies on Killer-Rescue genetic constructs, and has been predicted to have limited spatial and temporal spread. We simulate various operational approaches for deploying these genetic strategies, and quantify the optimal number of released transgenic mosquitoes needed to achieve definitive spread of Medea-linked genes and/or high frequencies of Killer-Rescue-associated elements. We show that for both strategies the most efficient approach for achieving spread of anti-pathogen genes within three years is generally to release adults of both sexes in multiple releases over time. Even though females in these releases should not transmit disease, there could be public concern over such releases, making the less efficient male-only release more practical. This study provides guidelines for operational approaches to population replacement genetic strategies, as well as illustrates the use of detailed spatial models to assist in safe and efficient implementation of such novel genetic strategies.

Suppression of dengue and malaria through releases of genetically engineered mosquitoes might soon become feasible. Aedes aegypti mosquitoes carrying a conditionally lethal transgene have recently been used to suppress local vector populations in small-scale field releases. Prior to releases of transgenic insects on a wider scale, however, most regulatory authorities will require additional evidence that suppression will be effective in natural heterogeneous habitats. We use a spatially explicit stochastic model of an Ae. aegypti population in Iquitos, Peru, along with an uncertainty analysis of its predictions, to quantitatively assess the outcome of varied operational approaches for releases of transgenic strains with conditional death of females. We show that population elimination might be an unrealistic objective in heterogeneous populations. We demonstrate that substantial suppression can nonetheless be achieved if releases are deployed in a uniform spatial pattern using strains combining multiple lethal elements, illustrating the importance of detailed spatial models for guiding genetic mosquito control strategies.

BackgroundSkeeter 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 findingsSkeeter 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.ConclusionsThis 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.

BackgroundDengue 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 findingsWe 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 F(ST) values and isolation by distance patterns, and examine the effects of adult dispersal and container movement between properties.Conclusions/significanceWe 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.