The evolution of parasitic and mutualistic plant-virus symbioses through transmission-virulence trade-offs.

Virus-plant interactions range from parasitism to mutualism. Viruses have been shown to increase fecundity of infected plants in comparison with uninfected plants under certain environmental conditions. Increased fecundity of infected plants may benefit both the plant and the virus as seed transmission is one of the main virus transmission pathways, in addition to vector transmission. Trade-offs between vertical (seed) and horizontal (vector) transmission pathways may involve virulence, defined here as decreased fecundity in infected plants. To better understand plant-virus symbiosis evolution, we explore the ecological and evolutionary interplay of virus transmission modes when infection can lead to an increase in plant fecundity. We consider two possible trade-offs: vertical seed transmission vs infected plant fecundity, and horizontal vector transmission vs infected plant fecundity (virulence). Through mathematical models and numerical simulations, we show (1) that a trade-off between virulence and vertical transmission can lead to virus extinction during the course of evolution, (2) that evolutionary branching can occur with subsequent coexistence of mutualistic and parasitic virus strains, and (3) that mutualism can out-compete parasitism in the long-run. In passing, we show that ecological bi-stability is possible in a very simple discrete-time epidemic model. Possible extensions of this study include the evolution of conditional (environment-dependent) mutualism in plant viruses

A model of population dynamics of TB in a prison system and application to South Africa

BACKGROUND: Tuberculosis (TB) continues to spread in South African prisons in particular, as prisons are over-capacitated and have poor ventilation. The awaiting trial detainees are not screened on admission and are at high risk of getting infected with TB. RESULTS: We propose a compartmental model to describe the population dynamics of TB disease in prisons. Our model considers the inflow of susceptible, exposed and TB infectives into the prison population. Removal of individuals out of the prison population can be either by death or by being released from prison, as compared to a general population in which removal is only by death. We describe conditions, including non-inflow of infectives into the prison, which will ensure that TB can be eradicated from the prison population. The model is calibrated for the South African prison system, by using data in existing literature. The model can be used to make quantitative projections of TB prevalence and to measure the effect of interventions. Illustrative simulations in this regard are presented. The model can be used for other prison populations too, if data is available to calculate the model parameters. CONCLUSIONS: Various simulations generated with our model serve to illustrate how it can be utilized in making future projections of the levels of prevalence of TB, and to quantify the effect of interventions such as screening, treatment or reduction of transmission parameter values through improved living conditions for inmates. This makes it particularly useful as there are various targets set by the World Health Organization and by governments, for reduction of TB prevalence and ultimately its eradication. Towards eradication of TB from a prison system, the theorem on global stability of the disease-free state is a useful indicator

Effect of spatial configuration of an extended nonlinear Kierstead-Slobodkin reaction-transport model with adaptive numerical scheme

In this paper, we consider the numerical simulations of an extended nonlinear form of Kierstead-Slobodkin reaction-transport system in one and two dimensions. We employ the popular fourth-order exponential time differencing Runge-Kutta (ETDRK4) schemes proposed by Cox and Matthew (J Comput Phys 176:430-455, 2002), that was modified by Kassam and Trefethen (SIAM J Sci Comput 26:1214-1233, 2005), for the time integration of spatially discretized partial differential equations. We demonstrate the supremacy of ETDRK4 over the existing exponential time differencing integrators that are of standard approaches and provide timings and error comparison. Numerical results obtained in this paper have granted further insight to the question "What is the minimal size of the spatial domain so that the population persists?" posed by Kierstead and Slobodkin (J Mar Res 12:141-147, 1953 ), with a conclusive remark that the popula- tion size increases with the size of the domain. In attempt to examine the biological wave phenomena of the solutions, we present the numerical results in both one- and two-dimensional space, which have interesting ecological implications. Initial data and parameter values were chosen to mimic some existing patternsScopus 201

Spatiotemporal Infectious Disease Modeling: A BME-SIR Approach

This paper is concerned with the modeling of infectious disease spread in a composite space-time domain under conditions of uncertainty. We focus on stochastic modeling that accounts for basic mechanisms of disease distribution and multi-sourced in situ uncertainties. Starting from the general formulation of population migration dynamics and the specification of transmission and recovery rates, the model studies the functional formulation of the evolution of the fractions of susceptible-infected-recovered individuals. The suggested approach is capable of: a) modeling population dynamics within and across localities, b) integrating the disease representation (i.e. susceptible-infected-recovered individuals) with observation time series at different geographical locations and other sources of information (e.g. hard and soft data, empirical relationships, secondary information), and c) generating predictions of disease spread and associated parameters in real time, while considering model and observation uncertainties. Key aspects of the proposed approach are illustrated by means of simulations (i.e. synthetic studies), and a real-world application using hand-foot-mouth disease (HFMD) data from China.J.M. Angulo and A.E. Madrid have been partially supported by grants MTM2009-13250 and MTM2012-32666 of SGPI, and P08-FQM-3834 of the Andalusian CICE, Spain. H-L Yu has been partially supported by a grant from National Science Council of Taiwan (NSC101-2628-E-002-017-MY3 and NSC102-2221-E-002-140-MY3). A. Kolovos was supported by SpaceTimeWorks, LLC. G. Christakos was supported by a Yongqian Chair Professorship (Zhejiang University, China)

A hierarchical kinetic theory of birth, death, and fission in age-structured interacting populations

We study mathematical models describing the evolution of stochastic age-structured populations. After reviewing existing approaches, we develop a complete kinetic framework for age-structured interacting populations undergoing birth, death and fission processes in spatially dependent environments. We define the full probability density for the population-size age chart and find results under specific conditions. Connections with more classical models are also explicitly derived. In particular, we show that factorial moments for non-interacting processes are described by a natural generalization of the McKendrick-von Foerster equation, which describes mean-field deterministic behavior. Our approach utilizes mixed-type, multidimensional probability distributions similar to those employed in the study of gas kinetics and with terms that satisfy BBGKY-like equation hierarchies