Stochastic dynamical aspects of neuronal activity

A stochastic model of neuronal activity is proposed. Some stochastic differential equations based on jump processes are used to investigate the behavior of the membrane potential at a time scale small with respect to the neuronal states time evolution. A model for learning, implying short memory effects, is described.     

Long-term evolution of an ecosystem with spontaneous periodicity of mass extinctions

Twenty years ago, after analysing palaeontological data, Raup and Sepkoski suggested that mass extinctions on Earth appear cyclically in time with a period of approximately 26 million years (My). To explain the 26 My period, a number of proposals were made involving, e.g., astronomical effects, increased volcanic acitivity, or the Earth's magnetic field reversal, none of which, however, has been confirmed. Here we study a spatially extended discrete model of an ecosystem and show that the periodicity of mass extinctions might be a natural feature of the ecosystem's dynamics and not the result of a periodic external perturbation. In our model, periodic changes of the diversity of an ecosystem and some of its other characteristics are induced by the coevolution of species. In agreement with some palaeontological data, our results show that the longevity of a species depends on the evolutionary stage at which the species is created. Possible further tests of our model are also discussed.     

Stochastic payoff evaluation increases the temperature of selection

We study stochastic evolutionary game dynamics in populations of finite size. Moreover, each individual has a randomly distributed number of interactions with other individuals. Therefore, the payoff of two individuals using the same strategy can be different. The resulting "payoff stochasticity" reduces the intensity of selection and therefore increases the temperature of selection. A simple mean-field approximation is derived that captures the average effect of the payoff stochasticity. Correction terms to the mean-field theory are computed and discussed.     

Evolution of Defence Portfolios in Exploiter–Victim Systems

Some organisms maintain a battery of defensive strategies against their exploiters (predators, parasites or parasitoids), while others fail to employ a defence that seems obvious. In this paper, we shall investigate the circumstances under which defence strategies might be expected to evolve. Brood parasites and their hosts provide our main motivation, and we shall discuss why the reed warbler Acrocephalus scirpaceus has evolved an egg-rejection but not a chick-rejection strategy as a defence against the common (Eurasian) cuckoo Cuculus canorus, while the superb fairy-wren Malurus cyaneus has evolved a chick-rejection but not an egg-rejection strategy as a defence against Horsfield's bronze-cuckoo Chrysococcyx basalis. We suggest that the answers lie in strategy-blocking, where one strategy (the blocking strategy) prevents the appearance of another (the blocked strategy) that would be adaptive in its absence. This may be common in exploiter–victim systems.     

Stochastic P systems and the simulation of biochemical processes with dynamic compartments

We introduce a sequential rewriting strategy for P systems based on Gillespie's stochastic simulation algorithm, and show that the resulting formalism of stochastic P systems makes it possible to simulate biochemical processes in dynamically changing, nested compartments. Stochastic P systems have been implemented using the spatially explicit programming language MGS. Implementation examples include models of the Lotka-Volterra auto-catalytic system, and the life cycle of the Semliki Forest virus.     

Prisoner's dilemma posed by fitness-associated recombination strategies

Genetic recombination is a central and repeated topic of study in the evolution of life. However, along with the influence of recombination on evolution, we understand surprisingly little of how selection shapes the nature of recombination. One explanation for recombination is that it allows organisms to escape from perilous situations where they experience very low fitness. As a corollary, it has been suggested that selection should favor recombination at low fitness and not at high fitness (fitness-associated recombination, FAR), and theory suggests that such strategies can indeed be selected. Here we develop models to further investigate the evolution of FAR. Consistent with previous works, we find that FAR can invade and dominate over a strategy of uniform recombination that is independent of fitness. However, our simulation results suggest that extreme FAR strategies, known as group-elitism, are not necessarily superior to other FAR strategies. Moreover, we argue that FAR domination will often occur with a net loss of mean population fitness. Interestingly, this suggests that the strategy of not recombining at high fitness will sometimes be analogous to a defector strategy from the famous "prisoner's dilemma" game: a selfish strategy that is selected but leads to a loss of mean fitness for all players.     

Strategy selection in structured populations

Evolutionary game theory studies frequency dependent selection. The fitness of a strategy is not constant, but depends on the relative frequencies of strategies in the population. This type of evolutionary dynamics occurs in many settings of ecology, infectious disease dynamics, animal behavior and social interactions of humans. Traditionally evolutionary game dynamics are studied in well-mixed populations, where the interaction between any two individuals is equally likely. There have also been several approaches to study evolutionary games in structured populations. In this paper we present a simple result that holds for a large variety of population structures. We consider the game between two strategies, A and B, described by the payoff matrix . We study a mutation and selection process. For weak selection strategy A is favored over B if and only if σa+b>c+σd. This means the effect of population structure on strategy selection can be described by a single parameter, σ. We present the values of σ for various examples including the well-mixed population, games on graphs, games in phenotype space and games on sets. We give a proof for the existence of such a σ, which holds for all population structures and update rules that have certain (natural) properties. We assume weak selection, but allow any mutation rate. We discuss the relationship between σ and the critical benefit to cost ratio for the evolution of cooperation. The single parameter, σ, allows us to quantify the ability of a population structure to promote the evolution of cooperation or to choose efficient equilibria in coordination games.     

The pace of evolution across fitness valleys

How fast does a population evolve from one fitness peak to another? We study the dynamics of evolving, asexually reproducing populations in which a certain number of mutations jointly confer a fitness advantage. We consider the time until a population has evolved from one fitness peak to another one with a higher fitness. The order of mutations can either be fixed or random. If the order of mutations is fixed, then the population follows a metaphorical ridge, a single path. If the order of mutations is arbitrary, then there are many ways to evolve to the higher fitness state. We address the time required for fixation in such scenarios and study how it is affected by the order of mutations, the population size, the fitness values and the mutation rate.     

Adaptive dynamics for physiologically structured population models

We develop a systematic toolbox for analyzing the adaptive dynamics of multidimensional traits in physiologically structured population models with point equilibria (sensu Dieckmann et al. in Theor. Popul. Biol. 63:309–338, 2003). Firstly, we show how the canonical equation of adaptive dynamics (Dieckmann and Law in J. Math. Biol. 34:579–612, 1996), an approximation for the rate of evolutionary change in characters under directional selection, can be extended so as to apply to general physiologically structured population models with multiple birth states. Secondly, we show that the invasion fitness function (up to and including second order terms, in the distances of the trait vectors to the singularity) for a community of N coexisting types near an evolutionarily singular point has a rational form, which is model-independent in the following sense: the form depends on the strategies of the residents and the invader, and on the second order partial derivatives of the one-resident fitness function at the singular point. This normal form holds for Lotka–Volterra models as well as for physiologically structured population models with multiple birth states, in discrete as well as continuous time and can thus be considered universal for the evolutionary dynamics in the neighbourhood of singular points. Only in the case of one-dimensional trait spaces or when N = 1 can the normal form be reduced to a Taylor polynomial. Lastly we show, in the form of a stylized recipe, how these results can be combined into a systematic approach for the analysis of the (large) class of evolutionary models that satisfy the above restrictions.      

Coping with an unpredictable and stressful environment: the life history and metabolic response to variable food and host availability in a polyphagous tephritid fly

The way energy resources are used under variable environmental conditions lies at the heart of our understanding of resource management and opportunism in many organisms. Here we sought to determine how a time-limited, synovigenic and polyphagous insect with a high reproductive-potential (Anastrephaludens), copes behaviourally and metabolically with environmental unpredictability represented by constant and variable regimes of host availability and variation in food quality. We hypothesized that an adaptive response to a windfall of nutritious food would be the rapid accumulation of energy metabolites (whole body lipids, glycogen and proteins) in the female. We also studied patterns of oogenesis as an indicator of egg-reabsorption under stressful environmental conditions. As predicted, patterns of energy metabolites were mainly driven by the quality and temporal pattern of food availability. In contrast, patterns of host availability had a lower impact upon metabolites. When given constant access to high quality nutrients, after an initial increase early in life, whole body lipids and glycogen were regulated downward to a steady-state level and somatic protein levels did not vary. In contrast, when food uncertainty was introduced, whole body lipid, glycogen and protein oscillated sharply with peaks associated with pulses of high-quality food. Production of eggs was highest when offered continuous access to hosts and high quality food. Importantly, females fully recovered their reproductive capacity when fruit became available following a period of host deprivation. With no evidence of egg resorption and high levels of egg dumping, it appears that egg dumping may favour the continuous production of eggs such that the female’s reproductive tissues are ready to respond to rapid changes in the availability of hosts. Our results exemplify the capacity of insects to maximize reproduction under variable and stressful environmental conditions.     

Bifurcation analysis of piecewise smooth ecological models

The aim of this paper is the study of the long-term behavior of population communities described by piecewise smooth models (known as Filippov systems). Models of this kind are often used to describe populations with selective switching between alternative habitats or diets or to mimic the evolution of an exploited resource where harvesting is forbidden when the resource is below a prescribed threshold. The analysis is carried out by performing the bifurcation analysis of the model with respect to two parameters. A relatively simple method, called the puzzle method, is proposed to construct the complete bifurcation diagram step-by-step. The method is illustrated through four examples concerning the exploitation and protection of interacting populations.     

具有比率依赖型功能反应函数的食饵-捕食者征税模型分析

本文研究了一个具有功能反应函数的食饵一捕食者征税模型,得到了该系统正平衡点的存在性、局部渐近稳定性和全局渐近稳定性的条件,并利用Pontrjagin最大值原理得到了最优税收策略．该文为资源管理者制定合理的管理政策提供理论依据．     

Influence of stochastic perturbation on prey-predator systems

We analyse the influence of various stochastic perturbations on prey-predator systems. The prey-predator model is described by stochastic versions of a deterministic Lotka-Volterra system. We study long-time behaviour of both trajectories and distributions of the solutions. We indicate the differences between the deterministic and stochastic models.     

The stabilizing effect of a random environment

It is shown that the lottery competition model permits coexistence in a stochastic environment, but not in a constant environment. Conditions for coexistence and competitive exclusion are determined. Analysis of these conditions shows that the essential requirements for coexistence are overlapping generations and fluctuating birth rates which ensure that each species has periods when it is increasing. It is found that a species may persist provided only that it is favored sufficiently by the environment during favorable periods independently of the extent to which the other species is favored during its favorable periods.Coexistence is defined in terms of the stochastic boundedness criterion for species persistence. Using the lottery model as an example this criterion is justified and compared with other persistence criteria. Properties of the stationary distribution of population density are determined for an interesting limiting case of the lottery model and these are related to stochastic boundedness. An attempt is then made to relate stochastic boundedness for infinite population models to the behavior of finite population models.     

Evolution probabilities and phylogenetic distance of dinucleotides

We develop here an analytical evolution model based on a dinucleotide mutation matrix 16 x 16 with six substitution parameters associated with the three types of substitutions in the two dinucleotide sites. It generalizes the previous models based on the nucleotide mutation matrices 4 x 4. It determines at some time t the exact occurrence probabilities of dinucleotides mutating randomly according to these six substitution parameters. Furthermore, several properties and two applications of this model allow to derive 16 evolutionary analytical solutions of dinucleotides and also a dinucleotide phylogenetic distance. Finally, based on this mathematical model, the SED (Stochastic Evolution of Dinucleotides) web server has been developed for deriving evolutionary analytical solutions of dinucleotides.     

Mutation-selection equilibrium in games with multiple strategies

In evolutionary games the fitness of individuals is not constant but depends on the relative abundance of the various strategies in the population. Here we study general games among n strategies in populations of large but finite size. We explore stochastic evolutionary dynamics under weak selection, but for any mutation rate. We analyze the frequency dependent Moran process in well-mixed populations, but almost identical results are found for the Wright-Fisher and Pairwise Comparison processes. Surprisingly simple conditions specify whether a strategy is more abundant on average than 1/n, or than another strategy, in the mutation-selection equilibrium. We find one condition that holds for low mutation rate and another condition that holds for high mutation rate. A linear combination of these two conditions holds for any mutation rate. Our results allow a complete characterization of n×n games in the limit of weak selection.