Periodic coexistence in the chemostat with three species competing for three essential resources

A chemostat model of three species of microorganisms competing for three essential, growth-limiting nutrients is considered. J. Husiman and F.J. Weissing [Nature 402 (1999) 407] show numerically that this model can generate periodic oscillations. The present contribution is concerned with rigorous analysis regarding the existence of periodic oscillations in this model. Our analysis is based on the following observation made by Huisman and Weissing: there is a cyclic replacement of species, if each species becomes limited by the resource for which it is the intermediate competitor. Using a permanence theory, an index theory, and a Poincaré-Bendixson theory for three-dimensional competitive systems, we analytically succeed to give sufficient conditions for the existence of periodic orbits in the limit sets in this model. The results in this paper suggest that with a wide range of parameter values, sustained periodic oscillations of species abundances for the model are possible, without involving external disturbances. Our results also suggest that competition is not necessarily destructive, i.e., in the case of existence of sustained periodic oscillations, if one of three competitors is absent, one of the other two rivals cannot survive.     

Periodic difference equations, population biology and the Cushing-Henson conjectures

We show that for a k-periodic difference equation, if a periodic orbit of period r is globally asymptotically stable (GAS), then r must be a divisor of k. Moreover, if r divides k we construct a non-autonomous dynamical system having minimum period k and which has a GAS periodic orbit with minimum period r. Our method uses the technique of skew-product dynamical systems. Our methods are then applied to prove two conjectures of Cushing and Henson concerning a non-autonomous Beverton-Holt equation which arises in the study of the response of a population to a periodically fluctuating environmental force such as seasonal fluctuations in carrying capacity or demographic parameters like birth or death rates. We give an equality linking the average population with the growth rates and carrying capacities (in the 2-periodic case) which shows that out-of-phase oscillations in these quantities always have a deleterious effect on the average population. We give an example where in-phase oscillations cause the opposite to occur.     

A two variable delay model for the circadian rhythm of Neurospora crassa

A two variable model with delay in both the variables, is proposed for the circadian oscillations of protein concentrations in the fungal species Neurospora crassa. The dynamical variables chosen are the concentrations of FRQ and WC-1 proteins. Our model is a two variable simplification of the detailed model of Smolen et al. (J. Neurosci. 21 (2001) 6644) modeling circadian oscillations with interlocking positive and negative feedback loops, containing 23 variables. In our model, as in the case of Smolen's model, a sustained limit cycle oscillation takes place in both FRQ and WC-1 protein in continuous darkness, and WC-1 is anti-phase to FRQ protein, as observed in experiments. The model accounts for various characteristic features of circadian rhythms such as entrainment to light dark cycles, phase response curves and robustness to parameter variation and molecular fluctuations. Simulations are carried out to study the effect of periodic forcing of circadian oscillations by light-dark cycles. The periodic forcing resulted in a rich bifurcation diagram that includes quasiperiodicity and chaotic oscillations, depending on the magnitude of the periodic changes in the light controlled parameter. When positive feedback is eliminated, our model reduces to the generic one dimensional delay model of Lema et al. (J. Theor. Biol. 204 (2000) 565), delay model of the circadian pace maker with FRQ protein as the dynamical variable which represses its own production. This one-dimensional model also exhibits all characteristic features of circadian oscillations and gives rise to circadian oscillations which are reasonably robust to parameter variations and molecular noise.     

High frequency spikes in long period blood cell oscillations

Several hematological diseases are characterised by oscillations of various blood cell populations. Two of these are a variant of chronic myelogenous leukemia (CML) and cyclical neutropenia (CN). These oscillations typically have long periods ranging from 20 to 60 days, despite the fact that the stem cell cycling time is thought to be of the order of 2–3 days. Clinical data from humans and laboratory data from the grey collie animal model of CN is suggestive of the idea that these long period oscillations may also contain higher frequency spiky oscillations. We show how such oscillations can be understood in the context of slow periodic stem cell oscillations, by analysing a two component differential-delay equation model of stem cell and neutrophil populations.For Karl Hadeler, on his 70th birthday, leader, teacher, colleague and friend.     

Formulas for intrinsic noise evaluation in oscillatory genetic networks

The linear noise approximation is a useful method for stochastic noise evaluations in genetic regulatory networks, where the covariance equation described as a Lyapunov equation plays a central role. We discuss the linear noise approximation method for evaluations of an intrinsic noise in autonomously oscillatory genetic networks; in such oscillatory networks, the covariance equation becomes a periodic differential equation that provides generally an unbounded covariance matrix, so that the standard method of noise evaluation based on the covariance matrix cannot be adopted directly. In this paper, we develop a new method of noise evaluation in oscillatory genetic networks; first, we investigate structural properties, e.g., orbital stability and periodicity, of the solutions to the covariance equation given as a periodic Lyapunov differential equation by using the Floquet-Lyapunov theory, and propose a global measure for evaluating stochastic amplitude fluctuations on the periodic trajectory; we also derive an evaluation formula for the period fluctuation. Finally, we apply our method to a model of circadian oscillations based on negative auto-regulation of gene expression, and show validity of our method by comparing the evaluation results with stochastic simulations.     

Transient oscillations induced by delayed growth response in the chemostat

In this paper, in order to try to account for the transient oscillations observed in chemostat experiments, we consider a model of single species growth in a chemostat that involves delayed growth response. The time delay models the lag involved in the nutrient conversion process. Both monotone response functions and nonmonotone response functions are considered. The nonmonotone response function models the inhibitory effects of growth response of certain nutrients when concentrations are too high. By applying local and global Hopf bifurcation theorems, we prove that the model has unstable periodic solutions that bifurcate from unstable nonnegative equilibria as the parameter measuring the delay passes through certain critical values and that these local periodic solutions can persist, even if the delay parameter moves far from the critical (local) bifurcation values.When there are two positive equilibria, then positive periodic solutions can exist. When there is a unique positive equilibrium, the model does not have positive periodic oscillations and the unique positive equilibrium is globally asymptotically stable. However, the model can have periodic solutions that change sign. Although these solutions are not biologically meaningful, provided the initial data starts close enough to the unstable manifold of one of these periodic solutions they may still help to account for the transient oscillations that have been frequently observed in chemostat experiments. Numerical simulations are provided to illustrate that the model has varying degrees of transient oscillatory behaviour that can be controlled by the choice of the initial data.Mathematics Subject Classification: 34D20, 34K20, 92D25Research was partially supported by NSERC of Canada.This work was partly done while this author was a postdoc at McMaster.     

Entrainment to Periodic Initiation and Transition Rates in a Computational Model for Gene Translation

Periodic oscillations play an important role in many biomedical systems. Proper functioning of biological systems that respond to periodic signals requires the ability to synchronize with the periodic excitation. For example, the sleep/wake cycle is a manifestation of an internal timing system that synchronizes to the solar day. In the terminology of systems theory, the biological system must entrain or phase-lock to the periodic excitation. Entrainment is also important in synthetic biology. For example, connecting several artificial biological systems that entrain to a common clock may lead to a well-functioning modular system. The cell-cycle is a periodic program that regulates DNA synthesis and cell division. Recent biological studies suggest that cell-cycle related genes entrain to this periodic program at the gene translation level, leading to periodically-varying protein levels of these genes. The ribosome flow model (RFM) is a deterministic model obtained via a mean-field approximation of a stochastic model from statistical physics that has been used to model numerous processes including ribosome flow along the mRNA. Here we analyze the RFM under the assumption that the initiation and/or transition rates vary periodically with a common period . We show that the ribosome distribution profile in the RFM entrains to this periodic excitation. In particular, the protein synthesis pattern converges to a unique periodic solution with period . To the best of our knowledge, this is the first proof of entrainment in a mathematical model for translation that encapsulates aspects such as initiation and termination rates, ribosomal movement and interactions, and non-homogeneous elongation speeds along the mRNA. Our results support the conjecture that periodic oscillations in tRNA levels and other factors related to the translation process can induce periodic oscillations in protein levels, and may suggest a new approach for re-engineering genetic systems to obtain a desired, periodic, protein synthesis rate.     

Sustained and transient oscillations and chaos induced by delayed antiviral immune response in an immunosuppressive infection model

Sustained and transient oscillations are frequently observed in clinical data for immune responses in viral infections such as human immunodeficiency virus, hepatitis B virus, and hepatitis C virus. To account for these oscillations, we incorporate the time lag needed for the expansion of immune cells into an immunosuppressive infection model. It is shown that the delayed antiviral immune response can induce sustained periodic oscillations, transient oscillations and even sustained aperiodic oscillations (chaos). Both local and global Hopf bifurcation theorems are applied to show the existence of periodic solutions, which are illustrated by bifurcation diagrams and numerical simulations. Two types of bistability are shown to be possible: (i) a stable equilibrium can coexist with another stable equilibrium, and (ii) a stable equilibrium can coexist with a stable periodic solution.     

Uncertainty and predictability in population dynamics of a bitrophic ecological model: Mixed-mode oscillations,bistability and sensitivity to parameters

We consider a two-trophic ecological model comprising of two predators competing for their common prey. We cast the model into the framework of a singular perturbed system of equations in one fast variable (prey population density) and two slow variables (predator population densities), mimicking the common observation that the per-capita productivity rate decreases from bottom to top along the trophic levels in Nature. We assume that both predators exhibit Holling II functional response with one of the predators (territorial) having a density dependent mortality rate. Depending on the system parameters, the model exhibits small, intermediate and/or large fluctuations in the population densities. The large fluctuations correspond to periodic population outbreaks followed by collapses (commonly known as cycles of “boom and bust”). The small fluctuations arise due to a singular Hopf bifurcation in the system, and are ecologically more desirable. However, more interestingly, the system exhibits mixed-mode oscillations (which are concatenations of the large amplitude oscillations and the small amplitude oscillations) that indicate the adaptability of the species to prolong the time gap between successive cycles of boom and bust. Numerical simulations are carried out to demonstrate the extreme sensitivity of the system to initial conditions (chaos and bistability of limit cycles are observed) as well as to the system parameters (here we only show the sensitivity to the density dependent mortality rate of the territorial predator). This model throws light at the uncertainties in long term behaviors that are associated with a real ecological system. We show that even very small changes in the system parameters due to natural or human-induced causes can lead to a complete different ecological phenomenon, thus affecting the predictability of the density of the prey population. In this paper, we explain the mechanisms behind the irregular fluctuations in the population sizes in an attempt to understand the dynamics occurring in a natural population and also comment on the inherent uncertainties associated with the system.     

Evidence for bile acid-evoked oscillations of Ca2(+)-dependent K+ permeability unrelated to a D-myo-inositol 1,4,5-trisphosphate effect in isolated guinea pig liver cells

In single liver cells, the D-myo-inositol 1,4,5-triphosphate (InsP3)-dependent agonists such as noradrenaline and angiotensin II evoke oscillations in intracellular calcium [Ca2+]i resulting mostly from the periodic release and reuptake of calcium from intracellular stores. In the present work, we have reexamined the effects of these agonists and investigated whether the natural bile acid taurolithocholic acid 3-sulfate (TLC-S), which permeabilizes the endoplasmic reticulum, could initiate oscillations of [Ca2+]i. Oscillations of [Ca2+]i were monitored with the Ca2(+)-dependent K+ permeability in whole-cell voltage-clamped guinea pig liver cells. Our results confirm the presence of two types of oscillations induced by hormones. They could be distinguished by their frequency periods. The fast (type I) had periods ranging from 5 to 12 s and the slow (type II) from 60 to 240 s. They have been respectively attributed to second messenger- and receptor-controlled oscillations, respectively. Our results also show that TLC-S, as noradrenaline and angiotensin II, induced the activation of this Ca(+)-dependent K+ current and was able to reproduce both types of oscillations. The bile acid effect was not blocked by intracellular perfusion of heparin known to inhibit both InsP3 binding and InsP3-evoked Ca2+ release in several tissues. In these conditions, TLC-S only evoked type I oscillations, suggesting that these fluctuations could originate from a mechanism that is independent of InsP3 and is an intrinsic property of internal Ca2+ stores.     

Multi-Annual Fluctuations in Reconstructed Historical Time-Series of a European Lobster (Homarus gammarus) Population Disappear at Increased Exploitation Levels

Through the history of ecology, fluctuations of populations have been a dominating topic, and endogenous causes of fluctuations and oscillations have been recognized and studied for more than 80 years. Here we analyzed an historical dataset, covering more than 130 years, of European lobster (Homarus gammarus) catches. The data shows periodic fluctuations, which are first dampened and then disappear over time. The disappearance of the periodicity coincided with a substantial increase in fishing effort and the oscillations have not reappeared in the time series. The shifting baseline syndrome has changed our perception of not only the status of the stock, but also the regulating pressures. We describe the transition of a naturally regulated lobster population into a heavily exploited fisheries controlled stock. This is shown by the incorporation of environmental and endogenous processes in generalized additive models, autocorrelation functions and periodicity analyses of time-series.     

Delay model of circadian gene expression with two negative feedback loops

In this theoretical paper we propose a quantitative minimal model for circadian gene expression based on two negative feedback loops. We perform numerical simulations to analyse its dynamics and parameter sensitivities in free-running conditions, and verify the entrainability by a single periodic driver. We furthermore apply two simultaneously acting external drivers, leading to aperiodic oscillations in the case of a single-loop system. These can be turned into regular periodic oscillations by introduction of a second loop. Our studies confirm the increasing evidence that multiple feedback loops increase the robustness of regulatory systems, and stress the particular situation of systems that are close to transition from free-running oscillation to steady-state behaviour. We discuss possible molecular realisations of the featured feedback loops and suggest the application of complex patterns of external stimulation as a generally useful approach to assess the functionality of models of circadian systems.     

Cannibalism in an age-structured predator-prey system

Recently, Kohlmeier and Ebenhöh showed that cannibalism can stabilize population cycles in a Lotka-Volterra type predator-prey model. Population cycles in their model are due to the interaction between logistic population growth of the prey and a hyperbolic functional response. In this paper, we consider a predator-prey system where cyclic population fluctuations are due to the age structure in the predator species. It is shown that cannibalism is also a stabilizing mechanism when population oscillations are due to this age structure. We conclude that in predator-prey systems, cannibalism by predators can stabilize both externally generated (consumer-resource) as well as internally generated (agestructure) fluctuations.     

Modelling Hematopoiesis Mediated by Growth Factors With Applications to Periodic Hematological Diseases

Hematopoiesis is a complex biological process that leads to the production and regulation of blood cells. It is based upon differentiation of stem cells under the action of growth factors. A mathematical approach of this process is proposed to understand some blood diseases characterized by very long period oscillations in circulating blood cells. A system of three differential equations with delay, corresponding to the cell cycle duration, is proposed and analyzed. The existence of a Hopf bifurcation at a positive steady-state is obtained through the study of an exponential polynomial characteristic equation with delay-dependent coefficients. Numerical simulations show that long-period oscillations can be obtained in this model, corresponding to a destabilization of the feedback regulation between blood cells and growth factors, for reasonable cell cycle durations. These oscillations can be related to observations on some periodic hematological diseases (such as chronic myelogenous leukemia, for example). ¹Research was partially supported by the INRIA Futurs, ANUBIS Team. ²Research was partially supported by the NSF and the University of Miami.     

External drive to inhibitory cells induces alternating episodes of high- and low-amplitude oscillations

Electrical oscillations in neuronal network activity are ubiquitous in the brain and have been associated with cognition and behavior. Intriguingly, the amplitude of ongoing oscillations, such as measured in EEG recordings, fluctuates irregularly, with episodes of high amplitude alternating with episodes of low amplitude. Despite the widespread occurrence of amplitude fluctuations in many frequency bands and brain regions, the mechanisms by which they are generated are poorly understood. Here, we show that irregular transitions between sub-second episodes of high- and low-amplitude oscillations in the alpha/beta frequency band occur in a generic neuronal network model consisting of interconnected inhibitory and excitatory cells that are externally driven by sustained cholinergic input and trains of action potentials that activate excitatory synapses. In the model, we identify the action potential drive onto inhibitory cells, which represents input from other brain areas and is shown to desynchronize network activity, to be crucial for the emergence of amplitude fluctuations. We show that the duration distributions of high-amplitude episodes in the model match those observed in rat prefrontal cortex for oscillations induced by the cholinergic agonist carbachol. Furthermore, the mean duration of high-amplitude episodes varies in a bell-shaped manner with carbachol concentration, just as in mouse hippocampus. Our results suggest that amplitude fluctuations are a general property of oscillatory neuronal networks that can arise through background input from areas external to the network.     

Sources of Variability in a Synthetic Gene Oscillator

Synthetic gene oscillators are small, engineered genetic circuits that produce periodic variations in target protein expression. Like other gene circuits, synthetic gene oscillators are noisy and exhibit fluctuations in amplitude and period. Understanding the origins of such variability is key to building predictive models that can guide the rational design of synthetic circuits. Here, we developed a method for determining the impact of different sources of noise in genetic oscillators by measuring the variability in oscillation amplitude and correlations between sister cells. We first used a combination of microfluidic devices and time-lapse fluorescence microscopy to track oscillations in cell lineages across many generations. We found that oscillation amplitude exhibited high cell-to-cell variability, while sister cells remained strongly correlated for many minutes after cell division. To understand how such variability arises, we constructed a computational model that identified the impact of various noise sources across the lineage of an initial cell. When each source of noise was appropriately tuned the model reproduced the experimentally observed amplitude variability and correlations, and accurately predicted outcomes under novel experimental conditions. Our combination of computational modeling and time-lapse data analysis provides a general way to examine the sources of variability in dynamic gene circuits.     

The manipulation of calcium oscillations by harnessing self-organisation

This paper investigates how self-organisation might be harnessed for the manipulation and control of calcium oscillations. Calcium signalling mechanisms are responsible for a number of important functions within biological systems, such as fertilization, secretion, contraction, neuronal signalling and learning. In this paper, calcium oscillations are investigated as a biological periodic process. Within biological systems such periodic behaviour is one of the outcomes from self-organisation. The understanding of periodic processes in living systems can enable more accurate diagnosis and physiologically suitable clinical therapies to be proposed, for diseases such as cancer, epilepsy, cardiac diseases and other dynamic diseases. In this paper these ideas are investigated by means of the calcium-induced calcium release (CICR) model and a number of representative simulations of intra and inter-cellular calcium oscillations are used to illustrate the manipulation and control of these oscillations in normal and pathological situations.     

Nonlinear time-series modeling of vole population fluctuations

A central goal of population ecology is to understand and predict fluctuations in population numbers. Until recently, much of the debate focused on the issue of population regulation by density-dependent factors. In this paper, I describe an approach to nonlinear modeling of time-series data that is designed to go beyond this question by investigating the possibility of complex population dynamics, characterized by lags in regulation and periodic or chaotic oscillations. The questions motivating this approach are: what are relative contributions of endogenous vs. exogenous components of dynamics? Is the irregular component in fluctuations entirely due to exogenous noise, or do nonlinearities contribute to it, too? I describe the philosophy and the technical details of the nonlinear modeling approach, and then apply it to a collection of time-series data on vole population fluctuations in northern Europe. The results suggest that population dynamics of European voles undergo a latitudinal shift from stability to chaos. Dynamics in northern Fennoscandia are characterized by positive Lyapunov exponent estimates, and a high degree of short-term (one year ahead) predictability, suggesting a strong endogenous component. In more southerly populations estimated Lyapunov exponents are negative, and there is no one-step ahead predictability, suggesting that fluctuations are driven by exogenous factors.