Adequacy of SEIR models when epidemics have spatial structure: Ebola in Sierra Leone.

Dynamic SEIR (Susceptible, Exposed, Infectious, Removed) compartmental models provide a tool for predicting the size and duration of both unfettered and managed outbreaks-the latter in the context of interventions such as case detection, patient isolation, vaccination and treatment. The reliability of this tool depends on the validity of key assumptions that include homogeneity of individuals and spatio-temporal homogeneity. Although the SEIR compartmental framework can easily be extended to include demographic (e.g. age) and additional disease (e.g. healthcare workers) classes, dependence of transmission rates on time, and metapopulation structure, fitting such extended models is hampered by both a proliferation of free parameters and insufficient or inappropriate data. This raises the question of how effective a tool the basic SEIR framework may actually be. We go some way here to answering this question in the context of the 2014-2015 outbreak of Ebola in West Africa by comparing fits of an SEIR time-dependent transmission model to both country- and district-level weekly incidence data. Our novel approach in estimating the effective-size-of-the-populations-at-risk ( Neff) and initial number of exposed individuals ( E0) at both district and country levels, as well as the transmission function parameters, including a time-to-halving-the-force-of-infection ( tf/2) parameter, provides new insights into this Ebola outbreak. It reveals that the estimate R0 ≈ 1.7 from country-level data appears to seriously underestimate R0 ≈ 3.3 - 4.3 obtained from more spatially homogeneous district-level data. Country-level data also overestimate tf/2 ≈ 22 weeks, compared with 8-10 weeks from district-level data. Additionally, estimates for the duration of individual infectiousness is around two weeks from spatially inhomogeneous country-level data compared with 2.4-4.5 weeks from spatially more homogeneous district-level data, which estimates are rather high compared with most values reported in the literature. This article is part of the theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes'. This issue is linked with the subsequent theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control'

Mesoscale movement and recursion behaviors of Namibian black rhinos.

Background:Understanding rhino movement behavior, especially their recursive movements, holds significant promise for enhancing rhino conservation efforts, and protecting their habitats and the biodiversity they support. Here we investigate the daily, biweekly, and seasonal recursion behavior of rhinos, to aid conservation applications and increase our foundational knowledge about these important ecosystem engineers. Methods:Using relocation data from 59 rhinos across northern Namibia and 8 years of sampling efforts, we investigated patterns in 24-h displacement at dawn, dusk, midday, and midnight to examine movement behaviors at an intermediate scale and across daily behavioral modes of foraging and resting. To understand recursion patterns across animals' short and long-term ranges, we built T-LoCoH time use grids to estimate recursive movement by each individual. Comparing these grids to contemporaneous MODIS imagery, we investigated productivity's influence on short-term space use and recursion. Finally, we investigated patterns of recursion within a year's home range, measuring the time to return to the most intensively used patches. Results:Twenty four-hour displacements at dawn were frequently smaller than 24-h displacements at dusk or at midday and midnight resting periods. Recursion analyses demonstrated that short-term recursion was most common in areas of median rather than maximum NDVI values. Investigated across a full year, recursion analysis showed rhinos most frequently returned to areas within 8-21 days, though visits were also seen separated by months likely suggesting seasonality in range use. Conclusions:Our results indicate that rhinos may frequently stay within the same area of their home ranges for days at a time, and possibly return to the same general area days in a row especially during morning foraging bouts. Recursion across larger time scales is also evident, and likely a contributing mechanism for maintaining open landscapes and browsing lawns of the savanna

Consumer-Resource Dynamics: Quantity, Quality, and Allocation

CITATION: Getz, W. M. & Owen-Smith, N. 2011. Consumer-resource dynamics : quantity, quality, and allocation. PLoS ONE, 6(1): e14539, doi:10.1371/journal.pone.0014539.The original publication is available at http://journals.plos.org/plosoneBackground: The dominant paradigm for modeling the complexities of interacting populations and food webs is a system of coupled ordinary differential equations in which the state of each species, population, or functional trophic group is represented by an aggregated numbers-density or biomass-density variable. Here, using the metaphysiological approach to model consumer-resource interactions, we formulate a two-state paradigm that represents each population or group in a food web in terms of both its quantity and quality. Methodology and Principal Findings: The formulation includes an allocation function controlling the relative proportion of extracted resources to increasing quantity versus elevating quality. Since lower quality individuals senesce more rapidly than higher quality individuals, an optimal allocation proportion exists and we derive an expression for how this proportion depends on population parameters that determine the senescence rate, the per-capita mortality rate, and the effects of these rates on the dynamics of the quality variable. We demonstrate that oscillations do not arise in our model from quantity-quality interactions alone, but require consumer-resource interactions across trophic levels that can be stabilized through judicious resource allocation strategies. Analysis and simulations provide compelling arguments for the necessity of populations to evolve quality-related dynamics in the form of maternal effects, storage or other appropriate structures. They also indicate that resource allocation switching between investments in abundance versus quality provide a powerful mechanism for promoting the stability of consumer-resource interactions in seasonally forcing environments. Conclusions/Significance: Our simulations show that physiological inefficiencies associated with this switching can be favored by selection due to the diminished exposure of inefficient consumers to strong oscillations associated with the wellknown paradox of enrichment. Also our results demonstrate how allocation switching can explain observed growth patterns in experimental microbial cultures and discuss how our formulation can address questions that cannot be answered using the quantity-only paradigms that currently predominate. © 2011 Getz, Owen-Smith.http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0014539Publisher's versio

Effect of the SOS Response on the Mean Fitness of Unicellular Populations: A Quasispecies Approach

The goal of this paper is to develop a mathematical model that analyzes the selective advantage of the SOS response in unicellular organisms. To this end, this paper develops a quasispecies model that incorporates the SOS response. We consider a unicellular, asexually replicating population of organisms, whose genomes consist of a single, double-stranded DNA molecule, i.e. one chromosome. We assume that repair of post-replication mismatched base-pairs occurs with probability , and that the SOS response is triggered when the total number of mismatched base-pairs is at least . We further assume that the per-mismatch SOS elimination rate is characterized by a first-order rate constant . For a single fitness peak landscape where the master genome can sustain up to mismatches and remain viable, this model is analytically solvable in the limit of infinite sequence length. The results, which are confirmed by stochastic simulations, indicate that the SOS response does indeed confer a fitness advantage to a population, provided that it is only activated when DNA damage is so extensive that a cell will die if it does not attempt to repair its DNA

Selection against Accumulating Mutations in Niche-Preference Genes Can Drive Speciation

Our current understanding of sympatric speciation is that it occurs primarily through disruptive selection on ecological genes driven by competition, followed by reproductive isolation through reinforcement-like selection against inferior intermediates/heterozygotes. Our evolutionary model of selection on resource recognition and preference traits suggests a new mechanism for sympatric speciation. We find speciation can occur in three phases. First a polymorphism of functionally different phenotypes is established through evolution of specialization. On the gene level, regulatory functions have evolved in which some alleles are conditionally switched off (i.e. are silent). These alleles accumulate harmful mutations that potentially may be expressed in offspring through recombination. Second mating associated with resource preference invades because harmful mutations in parents are not expressed in the offspring when mating assortatively, thereby dividing the population into two pre-zygotically isolated resource-specialist lineages. Third, silent alleles that evolved in phase one now accumulate deleterious mutations over the following generations in a Bateson-Dobzhansky-Muller fashion, establishing a post-zygotic barrier to hybridization

Coevolution of exploiter specialization and victim mimicry can be cyclic and saltational

Darwin’s Principle of Divergence explains sympatric speciation as gradual and directional. Contradicting evidence suggests that species’ traits evolve saltationally. Here, we model coevolution in exploiter-victim systems. Victims (resource population) have heritable, mutable cue phenotypes with different levels of defense. Exploiters have heritable, mutable perceptual phenotypes. Our simulations reveal coevolution of victim mimicry and exploiter specialization in a saltational and reversible cycle. Evolution is gradual and directional only in the specialization phase of the cycle thereby implying that specialization itself is saltational in such systems. Once linked to assortative mating, exploiter specialization provides conditions for speciation

Gray Wolves as Climate Change Buffers in Yellowstone

Understanding the mechanisms by which climate and predation patterns by top predators co-vary to affect community structure accrues added importance as humans exert growing influence over both climate and regional predator assemblages. In Yellowstone National Park, winter conditions and reintroduced gray wolves (Canis lupus) together determine the availability of winter carrion on which numerous scavenger species depend for survival and reproduction. As climate changes in Yellowstone, therefore, scavenger species may experience a dramatic reshuffling of food resources. As such, we analyzed 55 y of weather data from Yellowstone in order to determine trends in winter conditions. We found that winters are getting shorter, as measured by the number of days with snow on the ground, due to decreased snowfall and increased number of days with temperatures above freezing. To investigate synergistic effects of human and climatic alterations of species interactions, we used an empirically derived model to show that in the absence of wolves, early snow thaw leads to a substantial reduction in late-winter carrion, causing potential food bottlenecks for scavengers. In addition, by narrowing the window of time over which carrion is available and thereby creating a resource pulse, climate change likely favors scavengers that can quickly track food sources over great distances. Wolves, however, largely mitigate late-winter reduction in carrion due to earlier snow thaws. By buffering the effects of climate change on carrion availability, wolves allow scavengers to adapt to a changing environment over a longer time scale more commensurate with natural processes. This study illustrates the importance of restoring and maintaining intact food chains in the face of large-scale environmental perturbations such as climate change

Cooperation and Stability through Periodic Impulses

Basic games, where each individual chooses between two strategies, illustrate several issues that immediately emerge from the standard approach that applies strategic reasoning, based on rational decisions, to predict population behavior where no rationality is assumed. These include how mutual cooperation (which corresponds to the best outcome from the population perspective) can evolve when the only individually rational choice is to defect, illustrated by the Prisoner's Dilemma (PD) game, and how individuals can randomize between two strategies when neither is individually rational, illustrated by the Battle of the Sexes (BS) game that models male-female conflict over parental investment in offspring. We examine these questions from an evolutionary perspective where the evolutionary dynamics includes an impulsive effect that models sudden changes in collective population behavior. For the PD game, we show analytically that cooperation can either coexist with defection or completely take over the population, depending on the strength of the impulse. By extending these results for the PD game, we also show that males and females each evolve to a single strategy in the BS game when the impulsive effect is strong and that weak impulses stabilize the randomized strategies of this game

LoCoH: nonparameteric kernel methods for constructing home ranges and utilization distributions.

Parametric kernel methods currently dominate the literature regarding the construction of animal home ranges (HRs) and utilization distributions (UDs). These methods frequently fail to capture the kinds of hard boundaries common to many natural systems. Recently a local convex hull (LoCoH) nonparametric kernel method, which generalizes the minimum convex polygon (MCP) method, was shown to be more appropriate than parametric kernel methods for constructing HRs and UDs, because of its ability to identify hard boundaries (e.g., rivers, cliff edges) and convergence to the true distribution as sample size increases. Here we extend the LoCoH in two ways: "fixed sphere-of-influence," or r-LoCoH (kernels constructed from all points within a fixed radius r of each reference point), and an "adaptive sphere-of-influence," or a-LoCoH (kernels constructed from all points within a radius a such that the distances of all points within the radius to the reference point sum to a value less than or equal to a), and compare them to the original "fixed-number-of-points," or k-LoCoH (all kernels constructed from k-1 nearest neighbors of root points). We also compare these nonparametric LoCoH to parametric kernel methods using manufactured data and data collected from GPS collars on African buffalo in the Kruger National Park, South Africa. Our results demonstrate that LoCoH methods are superior to parametric kernel methods in estimating areas used by animals, excluding unused areas (holes) and, generally, in constructing UDs and HRs arising from the movement of animals influenced by hard boundaries and irregular structures (e.g., rocky outcrops). We also demonstrate that a-LoCoH is generally superior to k- and r-LoCoH (with software for all three methods available at http://locoh.cnr.berkeley.edu)