Durable Resistance to Crop Pathogens: An Epidemiological Framework to Predict Risk under Uncertainty

Increasing the durability of crop resistance to plant pathogens is one of the key goals of virulence management. Despite the recognition of the importance of demographic and environmental stochasticity on the dynamics of an epidemic, their effects on the evolution of the pathogen and durability of resistance has not received attention. We formulated a stochastic epidemiological model, based on the Kramer-Moyal expansion of the Master Equation, to investigate how random fluctuations affect the dynamics of an epidemic and how these effects feed through to the evolution of the pathogen and durability of resistance. We focused on two hypotheses: firstly, a previous deterministic model has suggested that the effect of cropping ratio (the proportion of land area occupied by the resistant crop) on the durability of crop resistance is negligible. Increasing the cropping ratio increases the area of uninfected host, but the resistance is more rapidly broken; these two effects counteract each other. We tested the hypothesis that similar counteracting effects would occur when we take account of demographic stochasticity, but found that the durability does depend on the cropping ratio. Secondly, we tested whether a superimposed external source of stochasticity (for example due to environmental variation or to intermittent fungicide application) interacts with the intrinsic demographic fluctuations and how such interaction affects the durability of resistance. We show that in the pathosystem considered here, in general large stochastic fluctuations in epidemics enhance extinction of the pathogen. This is more likely to occur at large cropping ratios and for particular frequencies of the periodic external perturbation (stochastic resonance). The results suggest possible disease control practises by exploiting the natural sources of stochasticity.     

Chaos and the (un)predictability of evolution in a changing environment

Among the factors that may reduce the predictability of evolution, chaos, characterized by a strong dependence on initial conditions, has received much less attention than randomness due to genetic drift or environmental stochasticity. It was recently shown that chaos in phenotypic evolution arises commonly under frequency‐dependent selection caused by competitive interactions mediated by many traits. This result has been used to argue that chaos should often make evolutionary dynamics unpredictable. However, populations also evolve largely in response to external changing environments, and such environmental forcing is likely to influence the outcome of evolution in systems prone to chaos. We investigate how a changing environment causing oscillations of an optimal phenotype interacts with the internal dynamics of an eco‐evolutionary system that would be chaotic in a constant environment. We show that strong environmental forcing can improve the predictability of evolution by reducing the probability of chaos arising, and by dampening the magnitude of chaotic oscillations. In contrast, weak forcing can increase the probability of chaos, but it also causes evolutionary trajectories to track the environment more closely. Overall, our results indicate that, although chaos may occur in evolution, it does not necessarily undermine its predictability.     

种群统计随机性和环境随机性对种群绝灭的影响

影响种群绝灭的随机干扰可分为种群统计随机性、环境随机性和随机灾害三大类。在相对稳定的环境条件下和相对较短的时间内,以前两类随机干扰对种群绝灭的影响为生态学家关注的焦点。但是,由于自然种群动态及其影响因子的复杂特征,进一步深入研究随机干扰对种群绝灭的作用在理论上和实践上都必须发展新的技术手段。本文回顾了种群统计随机性与环境随机性的概念起源与发展,系统阐述了其分析方法。归纳了两类随机性在种群绝灭研究中的应用范围、作用方式和特点的异同和区别方法。各类随机作用与种群动态之间关系的理论研究与对种群绝灭机理的实践研究紧密相关。根据理论模型模拟和自然种群实际分析两方面的研究现状,作者提出了进一步深入研究随机作用与种群非线性动态方法的策略。指出了随机干扰影响种群绝灭过程的研究的方向：更多的研究将从单纯的定性分析随机干扰对种群动力学简单性质的作用,转向结合特定的种群非线性动态特征和各类随机力作用特点具体分析绝灭极端动态的成因,以期做出精确的预测。     

Establishment versus population growth in spatio-temporally varying environments

We consider situations where repeated invasion attempts occur from a source population into a receptor population over extended periods of time. The receptor population contains two locations that provide different expected offspring numbers to invaders. There is demographic stochasticity in offspring numbers. In addition, temporal variation causes local invader fitnesses to vary. We show that effects of environmental autocorrelation on establishment success depend on spatial covariance of the receptor subpopulations. In situations with a low spatial covariance this effect is positive, whereas high spatial covariance and/or high migration probabilities between the subpopulations causes the effect to be negative. This result reconciles seemingly contradictory results from the literature concerning effects of temporal variation on population dynamics with demographic stochasticity. We study an example in the context of genetic introgression, where invasions of cultivar plant genes occur through pollen flow from a source population into wild-type receptor populations, but our results have implications in a wider range of contexts, such as the spread of exotic species, metapopulation dynamics and epidemics.     

Using mechanistic models to understand synchrony in forest insect populations: the North American gypsy moth as a case study

In many forest insects, subpopulations fluctuate concurrently across large geographical areas, a phenomenon known as population synchrony. Because of the large spatial scales involved, empirical tests to identify the causes of synchrony are often impractical. Simple models are, therefore, a useful aid to understanding, but data often seem to contradict model predictions. For instance, chaotic population dynamics and limited dispersal are not uncommon among synchronous forest defoliators, yet both make it difficult to achieve synchrony in simple models. To test whether this discrepancy can be explained by more realistic models, we introduced dispersal and spatially correlated stochasticity into a mechanistic population model for the North American gypsy moth Lymantria dispar. The resulting model shows both chaotic dynamics and spatial synchrony, suggesting that chaos and synchrony can be reconciled by the incorporation of realistic dynamics and spatial structure. By relating alterations in model structure to changes in synchrony levels, we show that the synchrony is due to a combination of spatial covariance in environmental stochasticity and the origins of chaos in our multispecies model.     

Spatially heterogeneous stochasticity and the adaptive diversification of dormancy

Diversified bet‐hedging, a strategy that leads several individuals with the same genotype to express distinct phenotypes in a given generation, is now well established as a common evolutionary response to environmental stochasticity. Life‐history traits defined as diversified bet‐hedging (e.g. germination or diapause strategies) display marked differences between populations in spatial proximity. In order to find out whether such differences can be explained by local adaptations to spatially heterogeneous environmental stochasticity, we explored the evolution of bet‐hedging dormancy strategies in a metapopulation using a two‐patch model with patch differences in stochastic juvenile survival. We found that spatial differences in the level of environmental stochasticity, restricted dispersal, increased fragmentation and intermediate survival during dormancy all favour the adaptive diversification of bet‐hedging dormancy strategies. Density dependency also plays a major role in the diversification of dormancy strategies because: (i) it may interact locally with environmental stochasticity and amplify its effects; however, (ii) it can also generate chaotic population dynamics that may impede diversification. Our work proposes new hypotheses to explain the spatial patterns of bet‐hedging strategies that we hope will encourage new empirical studies of this topic.     

Chemotactic Sensing towards Ambient and Secreted Attractant Drives Collective Behaviour of E. coli

We simulate the dynamics of a suspension of bacterial swimmers, which chemotactically sense gradients in either ambient or self-secreted attractants (e.g. nutrient or aspartate respectively), or in both. Unlike previous mean field models based on a set of continuum partial differential equations, our model resolves single swimmers and therefore incorporates stochasticity and effects due to fluctuations in the bacterial density field. The algorithm we use is simple enough that we can follow the evolution of colonies of up to over a million bacteria for timescales relevant to pattern formation for E. coli growing in semisolid medium such as agar, or in confined geometries. Our results confirm previous mean field results that the patterns observed experimentally can be reproduced with a model incorporating chemoattractant secretion, chemotaxis (towards gradients in the chemoattractant field), and bacterial reproduction. They also suggest that further experiments with bacterial strains chemotactically moving up both nutrient and secreted attractant field may yield yet more dynamical patterns.     

Intraindividual variability in behavior shapes fitness landscapes

Although intraindividual variability (IIV) in behavior is fundamental to ecological dynamics, the factors that contribute to the expression of IIV are poorly understood. Using an individual‐based model, this study examined the effects of stochasticity on the evolution of IIV represented by the residual variability of behavior. The model describes a population of prey with nonoverlapping generations, in which prey take refuge upon encountering a predator. The strategy of a prey is characterized by the mean and IIV (i.e., standard deviation) of hiding duration. Prey with no IIV will spend the same duration hiding in a refuge at each predator encounter, while prey with IIV will have variable hiding durations among encounters. For the sources of stochasticity, within‐generation stochasticity (represented by random predator encounters) and between‐generation stochasticity (represented by random resource availability) were considered. Analysis of the model indicates that individuals with high levels of IIV are maintained in a population in the presence of between‐generation stochasticity even though the optimal strategy in each generation is a strategy with no IIV, regardless of the presence or absence of within‐generation stochasticity. This contradictory pattern emerges because the mean behavioral trait and IIV do not independently influence fitness (e.g., the sign of the selection gradient with respect to IIV depends on the mean trait). Consequently, even when evolution eventually leads toward a strategy with no IIV (i.e., the optimal strategy), greater IIV may be transiently selected. Between‐generation stochasticity consistently imposes such transient selection and maintain high levels of IIV in a population.     

The interplay between determinism and stochasticity in childhood diseases

An important issue in the history of ecology has been the study of the relative importance of deterministic forces and processes noise in shaping the dynamics of ecological populations. We address this question by exploring the temporal dynamics of two childhood infections, measles and whooping cough, in England and Wales. We demonstrate that epidemics of whooping cough are strongly influenced by stochasticity; fully deterministic approaches cannot achieve even a qualitative fit to the observed data. In contrast, measles dynamics are extremely well explained by a deterministic model. These differences are shown to be caused by their contrasting responses to dynamical noise due to different infectious periods.     

Interspecific differences in stochastic population dynamics explains variation in Taylor's temporal power law

Taylor’s power law, i.e. that the slope for the increase in variance with mean population size is between 1 and 2 at a logarithmic scale, provides one of the few quantitative relationships in population ecology, yet the underlying ecological mechanisms are only poorly understood. Stochastic theory of population dynamics predicts that demographic and environmental stochasticity will affect the slope differently. In a stable environment under the influence of demographic stochasticity alone the slope will be equal to 1. In large populations in which demographic variance will have a negligible effect on the dynamics the slope will approach 2. In addition, the slope will also be influenced by how the strength of density dependence is related to mean population size. To disentangle the relative contribution of these processes we estimate the mean‐variance relationship for a large number of populations of British birds. The variance in population size of most species decreased with the mean due to decreased influence of demographic stochasticity at larger population sizes. Interspecific differences in demographic stochasticity was the main factor influencing variation in slopes of Taylor’s power law among species through a significant negative relationship between the slope and demographic variance. In addition, slopes were influenced by interspecific variation in life history parameters such as adult survival and clutch size. These analyses show that Taylor’s power law is generated from an interplay between stochastic and density dependent factors, modulated by life history.     

The effect of intrinsic stochasticity on transmitted HIV drug resistance patterns

Estimates of transmitted HIV drug-resistance prevalence vary widely among and within epidemiological surveys. Interpretation of trends from available survey data is therefore difficult. Because the emergence of drug-resistance involves small populations of infected drug-resistant individuals, the role of stochasticity (chance events) is likely to be important. The question addressed here is: how much variability in transmitted HIV drug-resistance prevalence patterns arises due to intrinsic stochasticity alone, i.e., if all starting conditions in the different epidemics surveyed were identical? This ‘thought experiment’ gives insight into the minimum expected variabilities within and among epidemics. A simple stochastic mathematical model was implemented. Our results show that stochasticity alone can generate a significant degree of variability and that this depends on the size and variation of the pool of new infections when drug treatment is first introduced. The variability in transmitted drug-resistance prevalence within an epidemic (i.e., the temporal variability) is large when the annual pool of all new infections is small (fewer than 200, typical of the HIV epidemics in Central European and Scandinavian countries) but diminishes rapidly as that pool grows. Epidemiological surveys involving hundreds of new infections annually are therefore needed to allow meaningful interpretation of temporal trends in transmitted drug-resistance prevalence within individual epidemics. The stochastic variability among epidemics shows a similar dependence on the pool of new infections if treatment is introduced after endemic equilibrium is established, but can persist even when there are more than 10,000 new infections annually if drug therapy is introduced earlier. Stochastic models may therefore have an important role to play in interpreting differences in transmitted drug-resistance prevalence trends among epidemiological surveys.     

Spatiotemporal dynamics in marginal populations

Population dynamics across a mortality gradient at an ecological margin are investigated using a novel modeling approach that allows direct comparison of stochastic spatially explicit simulation results with deterministic mean field models. The results show that demographic stochasticity has a large effect at population margins such that density profiles fall off more sharply than predicted by mean field models. Substantial spatial structure emerges at the margin, and spatial correlations (measured parallel to the margin) exhibit a sharp maximum in the tail of the density profile, indicating that spatial substructuring is greatest at an intermediate point across the ecological gradient. Such substructuring may have a substantial impact on Allee effects and evolutionary processes in marginal populations.     

Reconstructing Local Population Dynamics in Noisy Metapopulations—The Role of Random Catastrophes and Allee Effects

Reconstructing the dynamics of populations is complicated by the different types of stochasticity experienced by populations, in particular if some forms of stochasticity introduce bias in parameter estimation in addition to error. Identification of systematic biases is critical when determining whether the intrinsic dynamics of populations are stable or unstable and whether or not populations exhibit an Allee effect, i.e., a minimum size below which deterministic extinction should follow. Using a simulation model that allows for Allee effects and a range of intrinsic dynamics, we investigated how three types of stochasticity—demographic, environmental, and random catastrophes— affect our ability to reconstruct the intrinsic dynamics of populations. Demographic stochasticity aside, which is only problematic in small populations, we find that environmental stochasticity—positive and negative environmental fluctuations—caused increased error in parameter estimation, but bias was rarely problematic, except at the highest levels of noise. Random catastrophes, events causing large-scale mortality and likely to be more common than usually recognized, caused immediate bias in parameter estimates, in particular when Allee effects were large. In the latter case, population stability was predicted when endogenous dynamics were actually unstable and the minimum viable population size was overestimated in populations with small or non-existent Allee effects. Catastrophes also generally increased extinction risk, in particular when endogenous Allee effects were large. We propose a method for identifying data points likely resulting from catastrophic events when such events have not been recorded. Using social spider colonies (Anelosimus spp.) as models for populations, we show that after known or suspected catastrophes are accounted for, reconstructed growth parameters are consistent with intrinsic dynamical instability and substantial Allee effects. Our results are applicable to metapopulation or time series data and are relevant for predicting extinction in conservation applications or the management of invasive species.     

Length of intervals between epidemics: evaluating the influence of maternal transfer of immunity

The length of intervals between epidemic outbreaks of infectious diseases is critical in epidemiology. In several species of marine mammals and birds, it is pivotal to also consider the life history of the species of concern, as the contact rate between individuals can have a seasonal flux, for example, due to aggregations during the breeding season. Recently, particular interest has been given to the role of the dynamics of immunity in determining the intervals between epidemics in wild animal populations. One potentially powerful, but often neglected, process in this context is the maternal transfer of immunity. Here, we explore theoretically how the transfer of maternal antibodies can delay the recurrence of epidemics using Phocine Distemper in harbor seals as an example of a system in which epidemic outbreaks are followed by pathogen extinction. We show that the presence of temporarily protected newborns can significantly increase the predicted interval between epidemics, and this effect is strongly dependent on the degree of synchrony in the breeding season. Furthermore, we found that stochasticity in the onset of epidemics in combination with maternally acquired immunity increases the predicted intervals between epidemics even more. These effects arise because newborns with maternal antibodies temporarily boost population level immunity above the threshold of herd immunity, particularly when breeding is synchronous. Overall, our results show that maternal antibodies can have a profound influence on the dynamics of wildlife epidemics, notably in gregarious species such as many marine mammals and seabirds.     

Stochasticity in reproductive opportunity and the evolution of egg limitation in insects

Is reproduction by adult female insects limited by the finite time available to locate hosts (time limitation) or by the finite supply of eggs (egg limitation)? An influential model predicted that stochasticity in reproductive opportunity favors elevated fecundity, rendering egg limitation sufficiently rare that its importance would be greatly diminished. Here, I use models to explore how stochasticity shapes fecundity, the likelihood of egg limitation, and the ecological importance of egg limitation. The models make three predictions. First, whereas spatially stochastic environments favor increased fecundity, temporally stochastic environments favor increases, decreases, or intermediate maxima in fecundity, depending on egg costs. Second, even when spatially or temporally stochastic environments favor life histories with less‐frequent egg limitation, stochasticity still increases the proportion of all eggs laid in the population that is laid by females destined to become egg limited. This counterintuitive result is explained by noting that stochasticity concentrates reproduction in the hands of a few females that are likely to become egg limited. Third, spatially or temporally stochastic environments amplify the constraints imposed by time and eggs on total reproduction by the population. I conclude that both egg and time constraints are fundamental in shaping insect reproductive behavior and population dynamics in stochastic environments.     

Long food chains are in general chaotic

The question whether chaos exists in nature is much debated. In this paper we prove that chaotic parameter regions exist generically in food chains of length greater than three. While nonchaotic dynamics is also possible, the presence of chaotic parameter regions indicates that chaotic dynamics is likely. We show that the chaotic regions survive even at high exponents of closure. Our results have been obtained using a general food chain model that describes a large class of different food chains. The existence of chaos in models of such generality can be deduced from the presence of certain bifurcations of higher codimension.     

Invasion thresholds for fungicide resistance: deterministic and stochastic analyses

Fungicide resistance is an important practical problem, but one that is poorly understood at the population level. Here we introduce a simple nonlinear model for fungicide resistance in botanical epidemics which includes the dynamics of the chemical control agent and the host population, while also allowing for demographic stochasticity in the host-parasite dynamics. This provides a mathematical framework for analysing the risk of fungicide resistance developing by including the parameters for the amount applied, longevity and application frequency of the fungicide. The model demonstrates the existence of thresholds for the invasion of the resistant strain in the parasite population which depend on two quantities: the relative fitness of the resistant strain and the effectiveness of control. This threshold marks a change from definite elimination of the resistant strain below the threshold to a finite probability of invasion which increases above the threshold. The fungicide decay rate, the amount of fungicide applied and the period between applications affect the effectiveness of control and, consequently, they influence whether or not resistance develops and the time taken to achieve a critical frequency of resistance. All three parameters are amenable to control by the grower or by coordinating the activity of a population of growers. Providing crude estimates of the effectiveness of control and relative fitness are available, the results can be used to predict the consequences of changing these parameters for the risk of invasion and the proportion of sites at which this might be expected to occur. Although motivated for fungicide resistance, the model has broader application to herbicide, antibiotic and antiviral resistance. The modelling approach and results are discussed in the context of resistance to chemical control in general.     

Experimental support of the scaling rule for demographic stochasticity

A scaling rule of ecological theory, accepted but lacking experimental confirmation, is that the magnitude of fluctuations in population densities due to demographic stochasticity scales inversely with the square root of population numbers. This supposition is based on analyses of models exhibiting exponential growth or stable equilibria. Using two quantitative measures, we extend the scaling rule to situations in which population densities fluctuate due to nonlinear deterministic dynamics. These measures are applied to populations of the flour beetle Tribolium castaneum that display chaotic dynamics in both 20-g and 60-g habitats. Populations cultured in the larger habitat exhibit a clarification of the deterministic dynamics, which follows the inverse square root rule. Lattice effects, a deterministic phenomenon caused by the discrete nature of individuals, can cause deviations from the scaling rule when population numbers are small. The scaling rule is robust to the probability distribution used to model demographic variation among individuals.