Splitting the tail of the displacement kernel shows the unimportance of kurtosis

Delays in Epidemic Outbreak Control Cost Disproportionately Large Treatment Footprints to Offset

Epidemic outbreak control often involves a spatially explicit treatment area (quarantine, inoculation, ring cull) that covers the outbreak area and adjacent regions where hosts are thought to be latently infected. Emphasis on space however neglects the influence of treatment timing on outbreak control. We conducted field and in silico experiments with wheat stripe rust (WSR), a long-distance dispersed plant disease, to understand interactions between treatment timing and area interact to suppress an outbreak. Full-factorial field experiments with three different ring culls (outbreak area only to a 25-fold increase in treatment area) at three different disease control timings (1.125, 1.25, and 1.5 latent periods after initial disease expression) indicated that earlier treatment timing had a conspicuously greater suppressive effect than the area treated. Disease spread computer simulations over a broad range of influential epidemic parameter values (R0, outbreak disease prevalence, epidemic duration) suggested that potentially unrealistically large increases in treatment area would be required to compensate for even small delays in treatment timing. Although disease surveillance programs are costly, our results suggest that treatments early in an epidemic disease outbreak require smaller areas to be effective, which may ultimately compensate for the upfront costs of proactive disease surveillance programs.

Spatiotemporal Variation in Distance Dependent Animal Movement Contacts: One Size Doesn't Fit All

The structure of contacts that mediate transmission has a pronounced effect on the outbreak dynamics of infectious disease and simulation models are powerful tools to inform policy decisions. Most simulation models of livestock disease spread rely to some degree on predictions of animal movement between holdings. Typically, movements are more common between nearby farms than between those located far away from each other. Here, we assessed spatiotemporal variation in such distance dependence of animal movement contacts from an epidemiological perspective. We evaluated and compared nine statistical models, applied to Swedish movement data from 2008. The models differed in at what level (if at all), they accounted for regional and/or seasonal heterogeneities in the distance dependence of the contacts. Using a kernel approach to describe how probability of contacts between farms changes with distance, we developed a hierarchical Bayesian framework and estimated parameters by using Markov Chain Monte Carlo techniques. We evaluated models by three different approaches of model selection. First, we used Deviance Information Criterion to evaluate their performance relative to each other. Secondly, we estimated the log predictive posterior distribution, this was also used to evaluate their relative performance. Thirdly, we performed posterior predictive checks by simulating movements with each of the parameterized models and evaluated their ability to recapture relevant summary statistics. Independent of selection criteria, we found that accounting for regional heterogeneity improved model accuracy. We also found that accounting for seasonal heterogeneity was beneficial, in terms of model accuracy, according to two of three methods used for model selection. Our results have important implications for livestock disease spread models where movement is an important risk factor for between farm transmission. We argue that modelers should refrain from using methods to simulate animal movements that assume the same pattern across all regions and seasons without explicitly testing for spatiotemporal variation.

The ecological basis of morphogenesis: branching patterns in swarming colonies of bacteria

Understanding how large-scale shapes in tissues, organs and bacterial colonies emerge from local interactions among cells and how these shapes remain stable over time are two fundamental problems in biology. Here we investigate branching morphogenesis in an experimental model system, swarming colonies of the bacterium Pseudomonas aeruginosa. We combine experiments and computer simulation to show that a simple ecological model of population dispersal can describe the emergence of branching patterns. In our system, morphogenesis depends on two counteracting processes that act on different length-scales: (1) colony expansion, which increases the likelihood of colonizing a patch at a close distance and (2) colony repulsion, which decreases the colonization likelihood over a longer distance. The two processes are included in a kernel based mathematical model using an integro-differential approach borrowed from ecological theory. Computer simulations show that the model can indeed reproduce branching, but only for a narrow range of parameter values, suggesting that P. aeruginosa has a fine-tuned physiology for branching. Simulations further show that hyperswarming, a process where highly dispersive mutants reproducibly arise within the colony and disrupt branching patterns, can be interpreted as a change in the spatial kernel.

A bayesian approach for modeling cattle movements in the United States: scaling up a partially observed network

Networks are rarely completely observed and prediction of unobserved edges is an important problem, especially in disease spread modeling where networks are used to represent the pattern of contacts. We focus on a partially observed cattle movement network in the U.S. and present a method for scaling up to a full network based on Bayesian inference, with the aim of informing epidemic disease spread models in the United States. The observed network is a 10% state stratified sample of Interstate Certificates of Veterinary Inspection that are required for interstate movement; describing approximately 20,000 movements from 47 of the contiguous states, with origins and destinations aggregated at the county level. We address how to scale up the 10% sample and predict unobserved intrastate movements based on observed movement distances. Edge prediction based on a distance kernel is not straightforward because the probability of movement does not always decline monotonically with distance due to underlying industry infrastructure. Hence, we propose a spatially explicit model where the probability of movement depends on distance, number of premises per county and historical imports of animals. Our model performs well in recapturing overall metrics of the observed network at the node level (U.S. counties), including degree centrality and betweenness; and performs better compared to randomized networks. Kernel generated movement networks also recapture observed global network metrics, including network size, transitivity, reciprocity, and assortativity better than randomized networks. In addition, predicted movements are similar to observed when aggregated at the state level (a broader geographic level relevant for policy) and are concentrated around states where key infrastructures, such as feedlots, are common. We conclude that the method generally performs well in predicting both coarse geographical patterns and network structure and is a promising method to generate full networks that incorporate the uncertainty of sampled and unobserved contacts.

Retention time variability as a mechanism for animal mediated long-distance dispersal

Long-distance dispersal (LDD) events, although rare for most plant species, can strongly influence population and community dynamics. Animals function as a key biotic vector of seeds and thus, a mechanistic and quantitative understanding of how individual animal behaviors scale to dispersal patterns at different spatial scales is a question of critical importance from both basic and applied perspectives. Using a diffusion-theory based analytical approach for a wide range of animal movement and seed transportation patterns, we show that the scale (a measure of local dispersal) of the seed dispersal kernel increases with the organisms' rate of movement and mean seed retention time. We reveal that variations in seed retention time is a key determinant of various measures of LDD such as kurtosis (or shape) of the kernel, thinkness of tails and the absolute number of seeds falling beyond a threshold distance. Using empirical data sets of frugivores, we illustrate the importance of variability in retention times for predicting the key disperser species that influence LDD. Our study makes testable predictions linking animal movement behaviors and gut retention times to dispersal patterns and, more generally, highlights the potential importance of animal behavioral variability for the LDD of seeds.

Influence on disease spread dynamics of herd characteristics in a structured livestock industry

Studies of between-herd contacts may provide important insight to disease transmission dynamics. By comparing the result from models with different levels of detail in the description of animal movement, we studied how factors influence the final epidemic size as well as the dynamic behaviour of an outbreak. We investigated the effect of contact heterogeneity of pig herds in Sweden due to herd size, between-herd distance and production type. Our comparative study suggests that the production-type structure is the most influential factor. Hence, our results imply that production type is the most important factor to obtain valid data for and include when modelling and analysing this system. The study also revealed that all included factors reduce the final epidemic size and also have yet more diverse effects on initial rate of disease spread. This implies that a large set of factors ought to be included to assess relevant predictions when modelling disease spread between herds. Furthermore, our results show that a more detailed model changes predictions regarding the variability in the outbreak dynamics and conclude that this is an important factor to consider in risk assessment.

Should I stay or should I go? A habitat-dependent dispersal kernel improves prediction of movement

The analysis of animal movement within different landscapes may increase our understanding of how landscape features affect the perceptual range of animals. Perceptual range is linked to movement probability of an animal via a dispersal kernel, the latter being generally considered as spatially invariant but could be spatially affected. We hypothesize that spatial plasticity of an animal's dispersal kernel could greatly modify its distribution in time and space. After radio tracking the movements of walking insects (Cosmopolites sordidus) in banana plantations, we considered the movements of individuals as states of a Markov chain whose transition probabilities depended on the habitat characteristics of current and target locations. Combining a likelihood procedure and pattern-oriented modelling, we tested the hypothesis that dispersal kernel depended on habitat features. Our results were consistent with the concept that animal dispersal kernel depends on habitat features. Recognizing the plasticity of animal movement probabilities will provide insight into landscape-level ecological processes.

Bayesian analysis of animal movements related to factors at herd and between herd levels: Implications for disease spread modeling

A method to assess the influence of between herd distances, production types and herd sizes on patterns of between herd contacts is presented. It was applied on pig movement data from a central database of the Swedish Board of Agriculture. To determine the influence of these factors on the contact between holdings we used a Bayesian model and Markov chain Monte Carlo (MCMC) methods to estimate the posterior distribution of model parameters. The analysis showed that the contact pattern via animal movements is highly heterogeneous and influenced by all three factors, production type, herd size, and distance between holdings. Most production types showed a positive relationship between maximum capacity and the probability of both incoming and outgoing movements. In agreement with previous studies, holdings also differed in both the number of contacts as well as with what holding types contact occurred with. Also, the scale and shape of distance dependence in contact probability was shown to differ depending on the production types of holdings.To demonstrate how the methodology may be used for risk assessment, disease transmissions via animal movements were simulated with the model used for analysis of contacts, and parameterized by the analyzed posterior distribution. A Generalized Linear Model showed that herds with production types Sow pool center, Multiplying herd and Nucleus herd have higher risk of generating a large number of new infections. Multiplying herds are also expected to generate many long distance transmissions, while transmissions generated by Sow pool centers are confined to more local areas. We argue that the methodology presented may be a useful tool for improvement of risk assessment based on data found in central databases.

The shape of the spatial kernel and its implications for biological invasions in patchy environments

Ecological and epidemiological invasions occur in a spatial context. We investigated how these processes correlate to the distance dependence of spread or dispersal between spatial entities such as habitat patches or epidemiological units. Distance dependence is described by a spatial kernel, characterized by its shape (kurtosis) and width (variance). We also developed a novel method to analyse and generate point-pattern landscapes based on spectral representation. This involves two measures: continuity, which is related to autocorrelation and contrast, which refers to variation in patch density. We also analysed some empirical data where our results are expected to have implications, namely distributions of trees (Quercus and Ulmus) and farms in Sweden. Through a simulation study, we found that kernel shape was not important for predicting the invasion speed in randomly distributed patches. However, the shape may be essential when the distribution of patches deviates from randomness, particularly when the contrast is high. We conclude that the speed of invasions depends on the spatial context and the effect of the spatial kernel is intertwined with the spatial structure. This implies substantial demands on the empirical data, because it requires knowledge of shape and width of the spatial kernel, and spatial structure.