Persistence in epidemic metapopulations: quantifying the rescue effects for measles, mumps, rubella and whooping cough

Metapopulation rescue effects are thought to be key to the persistence of many acute immunizing infections. Yet the enhancement of persistence through spatial coupling has not been previously quantified. Here we estimate the metapopulation rescue effects for four childhood infections using global WHO reported incidence data by comparing persistence on island countries vs all other countries, while controlling for key variables such as vaccine cover, birth rates and economic development. The relative risk of extinction on islands is significantly higher, and approximately double the risk of extinction in mainland countries. Furthermore, as may be expected, infections with longer infectious periods tend to have the strongest metapopulation rescue effects. Our results quantitate the notion that demography and local community size controls disease persistence

From regional pulse vaccination to global disease eradication: insights from a mathematical model of Poliomyelitis

Mass-vaccination campaigns are an important strategy in the global fight against poliomyelitis and measles. The large-scale logistics required for these mass immunisation campaigns magnifies the need for research into the effectiveness and optimal deployment of pulse vaccination. In order to better understand this control strategy, we propose a mathematical model accounting for the disease dynamics in connected regions, incorporating seasonality, environmental reservoirs and independent periodic pulse vaccination schedules in each region. The effective reproduction number, $R_e$, is defined and proved to be a global threshold for persistence of the disease. Analytical and numerical calculations show the importance of synchronising the pulse vaccinations in connected regions and the timing of the pulses with respect to the pathogen circulation seasonality. Our results indicate that it may be crucial for mass-vaccination programs, such as national immunisation days, to be synchronised across different regions. In addition, simulations show that a migration imbalance can increase $R_e$ and alter how pulse vaccination should be optimally distributed among the patches, similar to results found with constant-rate vaccination. Furthermore, contrary to the case of constant-rate vaccination, the fraction of environmental transmission affects the value of $R_e$ when pulse vaccination is present.Comment: Added section 6.1, made other revisions, changed titl

Predators reduce extinction risk in noisy metapopulations

Background Spatial structure across fragmented landscapes can enhance regional population persistence by promoting local “rescue effects.” In small, vulnerable populations, where chance or random events between individuals may have disproportionately large effects on species interactions, such local processes are particularly important. However, existing theory often only describes the dynamics of metapopulations at regional scales, neglecting the role of multispecies population dynamics within habitat patches. Findings By coupling analysis across spatial scales we quantified the interaction between local scale population regulation, regional dispersal and noise processes in the dynamics of experimental host-parasitoid metapopulations. We find that increasing community complexity increases negative correlation between local population dynamics. A potential mechanism underpinning this finding was explored using a simple population dynamic model. Conclusions Our results suggest a paradox: parasitism, whilst clearly damaging to hosts at the individual level, reduces extinction risk at the population level

The impact of the demographic transition on dengue in Thailand: Insights from a statistical analysis and mathematical modeling

Background: An increase in the average age of dengue hemorrhagic fever (DHF) cases has been reported in Thailand. The cause of this increase is not known. Possible explanations include a reduction in transmission due to declining mosquito populations, declining contact between human and mosquito, and changes in reporting. We propose that a demographic shift toward lower birth and death rates has reduced dengue transmission and lengthened the interval between large epidemics. Methods and Findings: Using data from each of the 72 provinces of Thailand, we looked for associations between force of infection (a measure of hazard, defined as the rate per capita at which susceptible individuals become infected) and demographic and climactic variables. We estimated the force of infection from the age distribution of cases from 1985 to 2005. We find that the force of infection has declined by 2% each year since a peak in the late 1970s and early 1980s. Contrary to recent findings suggesting that the incidence of DHF has increased in Thailand, we find a small but statistically significant decline in DHF incidence since 1985 in a majority of provinces. The strongest predictor of the change in force of infection and the mean force of infection is the median age of the population. Using mathematical simulations of dengue transmission we show that a reduced birth rate and a shift in the population's age structure can explain the shift in the age distribution of cases, reduction of the force of infection, and increase in the periodicity of multiannual oscillations of DHF incidence in the absence of other changes. Conclusions: Lower birth and death rates decrease the flow of susceptible individuals into the population and increase the longevity of immune individuals. The increase in the proportion of the population that is immune increases the likelihood that an infectious mosquito will feed on an immune individual, reducing the force of infection. Though the force of infection has decreased by half, we find that the critical vaccination fraction has not changed significantly, declining from an average of 85% to 80%. Clinical guidelines should consider the impact of continued increases in the age of dengue cases in Thailand. Countries in the region lagging behind Thailand in the demographic transition may experience the same increase as their population ages. The impact of demographic changes on the force of infection has been hypothesized for other diseases, but, to our knowledge, this is the first observation of this phenomenon

Optimizing Metapopulation Sustainability through a Checkerboard Strategy

The persistence of a spatially structured population is determined by the rate of dispersal among habitat patches. If the local dynamic at the subpopulation level is extinction-prone, the system viability is maximal at intermediate connectivity where recolonization is allowed, but full synchronization that enables correlated extinction is forbidden. Here we developed and used an algorithm for agent-based simulations in order to study the persistence of a stochastic metapopulation. The effect of noise is shown to be dramatic, and the dynamics of the spatial population differs substantially from the predictions of deterministic models. This has been validated for the stochastic versions of the logistic map, the Ricker map and the Nicholson-Bailey host-parasitoid system. To analyze the possibility of extinction, previous studies were focused on the attractiveness (Lyapunov exponent) of stable solutions and the structure of their basin of attraction (dependence on initial population size). Our results suggest that these features are of secondary importance in the presence of stochasticity. Instead, optimal sustainability is achieved when decoherence is maximal. Individual-based simulations of metapopulations of different sizes, dimensions and noise types, show that the system's lifetime peaks when it displays checkerboard spatial patterns. This conclusion is supported by the results of a recently published Drosophila experiment. The checkerboard strategy provides a technique for the manipulation of migration rates (e.g., by constructing corridors) in order to affect the persistence of a metapopulation. It may be used in order to minimize the risk of extinction of an endangered species, or to maximize the efficiency of an eradication campaign

Evolutionary Game Theory and Social Learning Can Determine How Vaccine Scares Unfold

Immunization programs have often been impeded by vaccine scares, as evidenced by the measles-mumps-rubella (MMR) autism vaccine scare in Britain. A “free rider” effect may be partly responsible: vaccine-generated herd immunity can reduce disease incidence to such low levels that real or imagined vaccine risks appear large in comparison, causing individuals to cease vaccinating. This implies a feedback loop between disease prevalence and strategic individual vaccinating behavior. Here, we analyze a model based on evolutionary game theory that captures this feedback in the context of vaccine scares, and that also includes social learning. Vaccine risk perception evolves over time according to an exogenously imposed curve. We test the model against vaccine coverage data and disease incidence data from two vaccine scares in England & Wales: the whole cell pertussis vaccine scare and the MMR vaccine scare. The model fits vaccine coverage data from both vaccine scares relatively well. Moreover, the model can explain the vaccine coverage data more parsimoniously than most competing models without social learning and/or feedback (hence, adding social learning and feedback to a vaccine scare model improves model fit with little or no parsimony penalty). Under some circumstances, the model can predict future vaccine coverage and disease incidence—up to 10 years in advance in the case of pertussis—including specific qualitative features of the dynamics, such as future incidence peaks and undulations in vaccine coverage due to the population's response to changing disease incidence. Vaccine scares could become more common as eradication goals are approached for more vaccine-preventable diseases. Such models could help us predict how vaccine scares might unfold and assist mitigation efforts

Global Patterns in Seasonal Activity of Influenza A/H3N2, A/H1N1, and B from 1997 to 2005: Viral Coexistence and Latitudinal Gradients

Despite a mass of research on the epidemiology of seasonal influenza, overall patterns of infection have not been fully described on broad geographic scales and for specific types and subtypes of the influenza virus. Here we provide a descriptive analysis of laboratory-confirmed influenza surveillance data by type and subtype (A/H3N2, A/H1N1, and B) for 19 temperate countries in the Northern and Southern hemispheres from 1997 to 2005, compiled from a public database maintained by WHO (FluNet). Key findings include patterns of large scale co-occurrence of influenza type A and B, interhemispheric synchrony for subtype A/H3N2, and latitudinal gradients in epidemic timing for type A. These findings highlight the need for more countries to conduct year-round viral surveillance and report reliable incidence data at the type and subtype level, especially in the Tropics

Bursts of vertex activation and epidemics in evolving networks

The dynamic nature of contact patterns creates diverse temporal structures. In particular, empirical studies have shown that contact patterns follow heterogeneous inter-event time intervals, meaning that periods of high activity are followed by long periods of inactivity. To investigate the impact of these heterogeneities in the spread of infection from a theoretical perspective, we propose a stochastic model to generate temporal networks where vertices make instantaneous contacts following heterogeneous inter-event intervals, and may leave and enter the system. We study how these properties affect the prevalence of an infection and estimate , the number of secondary infections of an infectious individual in a completely susceptible population, by modeling simulated infections (SI and SIR) that co-evolve with the network structure. We find that heterogeneous contact patterns cause earlier and larger epidemics in the SIR model in comparison to homogeneous scenarios for a vast range of parameter values, while smaller epidemics may happen in some combinations of parameters. In the case of SI and heterogeneous patterns, the epidemics develop faster in the earlier stages followed by a slowdown in the asymptotic limit. For increasing vertex turnover rates, heterogeneous patterns generally cause higher prevalence in comparison to homogeneous scenarios with the same average inter-event interval. We find that is generally higher for heterogeneous patterns, except for sufficiently large infection duration and transmission probability

Genetic Diversity in the SIR Model of Pathogen Evolution

We introduce a model for assessing the levels and patterns of genetic diversity in pathogen populations, whose epidemiology follows a susceptible-infected-recovered model (SIR). We model the population of pathogens as a metapopulation composed of subpopulations (infected hosts), where pathogens replicate and mutate. Hosts transmit pathogens to uninfected hosts. We show that the level of pathogen variation is well predicted by analytical expressions, such that pathogen neutral molecular variation is bounded by the level of infection and increases with the duration of infection. We then introduce selection in the model and study the invasion probability of a new pathogenic strain whose fitness (R0(1+s)) is higher than the fitness of the resident strain (R0). We show that this invasion probability is given by the relative increment in R0 of the new pathogen (s). By analyzing the patterns of genetic diversity in this framework, we identify the molecular signatures during the replacement and compare these with those observed in sequences of influenza A

Model for in vivo progression of tumors based on co-evolving cell population and vasculature

With countless biological details emerging from cancer experiments, there is a growing need for minimal mathematical models which simultaneously advance our understanding of single tumors and metastasis, provide patient-personalized predictions, whilst avoiding excessive hard-to-measure input parameters which complicate simulation, analysis and interpretation. Here we present a model built around a co-evolving resource network and cell population, yielding good agreement with primary tumors in a murine mammary cell line EMT6-HER2 model in BALB/c mice and with clinical metastasis data. Seeding data about the tumor and its vasculature from in vivo images, our model predicts corridors of future tumor growth behavior and intervention response. A scaling relation enables the estimation of a tumor's most likely evolution and pinpoints specific target sites to control growth. Our findings suggest that the clinically separate phenomena of individual tumor growth and metastasis can be viewed as mathematical copies of each other differentiated only by network structure