Global analysis of a model of plasmid-bearing,plasmid-free competition in a chemostat

A model of competition between plasmid-bearing and plasmid-free organisms in a chemostat was proposed in a paper of Stephanopoulis and Lapidus. The model was in the form of a system of nonlinear ordinary differential equations. Such models are relevant to commercial production by genetically altered organisms in continuous culture. The analysis there was local (using index arguments). This paper provides a mathematically rigorous analysis of the global asymptotic behavior of the governing equations in the case of uninhibited specific growth rate.Research supported by the National Council of Science, Republic of ChinaResearch supported by National Science Foundation Grant, DMS-9204490Research supported by the Natural Science and Engineering Council of Canada. This author's contribution was made while on research leave visiting the Department of Ecology and Evolutionary Biology at Princeton University. She would especially like to thank Simon Levin for his guidance as well as for providing an exceptional working environment     

The periodically forced Droop model for phytoplankton growth in a chemostat

 It is proved that the periodically forced Droop model for phytoplankton growth in a chemostat has precisely two dynamic regimes depending on a threshold condition involving the dilution rate. If the dilution rate is such that the sub-threshold condition holds, the phytoplankton population is washed out of the chemostat. If the super-threshold condition holds, then there is a unique periodic solution, having the same period as the forcing, characterized by the presence of the phytoplankton population, to which all solutions approach asymptotically. Furthermore, this result holds for a general class of models with monotone growth rate and monotone uptake rate, the latter possibly depending on the cell quota. Received 10 October 1995; received in revised form 26 March 1996     

Global analysis of a model of plasmid-bearing,plasmid-free competition in a chemostat with inhibitions

A model of competition between plasmid-bearing and plasmid-free organisms in a chemostat was proposed in a paper of Stephanopoulis and Lapidus. The model was in the form of a system of nonlinear ordinary differential equations. Such models were relevant to commercial production by genetically altered organisms in continuous culture. The analysis there was local. The rigorous global analysis was done in a paper of Hsu, Waltman and Wolkowicz in the case of the uninhibited specific growth rates. This paper provides a mathematically rigorous analysis of the global asymptotic behavior of the governing equations in the cases of combinations of inhibited and uninhibited specific growth rates.Research Supported by the National Council of Science, Republic of China     

Long run coexistence in the chemostat with multiple species

In this work we analyze the transient behavior of the dynamics of multiple species competing in a chemostat for a single resource, presenting slow/fast characteristics. We prove that coexistence among a subset of species, with growth functions close to each other, can last for a substantially long time. For these cases, we also show that the proportion of non-dominant species can be increasing before decreasing, under certain conditions on the initial distribution.     

Global dynamics of a chemostat competition model with distributed delay

 We study the global dynamics of n-species competition in a chemostat with distributed delay describing the time-lag involved in the conversion of nutrient to viable biomass. The delay phenomenon is modelled by the gamma distribution. The linear chain trick and a fluctuation lemma are applied to obtain the global limiting behavior of the model. When each population can survive if it is cultured alone, we prove that at most one competitor survives. The winner is the population that has the smallest delayed break-even concentration, provided that the orders of the delay kernels are large and the mean delays modified to include the washout rate (which we call the virtual mean delays) are bounded and close to each other, or the delay kernels modified to include the washout factor (which we call the virtual delay kernels) are close in L ¹-norm. Also, when the virtual mean delays are relatively small, it is shown that the predictions of the distributed delay model are identical with the predictions of the corresponding ODEs model without delay. However, since the delayed break-even concentrations are functions of the parameters appearing in the delay kernels, if the delays are sufficiently large, the prediction of which competitor survives, given by the ODEs model, can differ from that given by the delay model. Received: 9 August 1997 / Revised version: 2 July 1998     

Coexistence in the chemostat as a result of metabolic by-products

Classical chemostat models assume that competition is purely exploitative and mediated via a common, limiting and single resource. However, in laboratory experiments with pathogens related to the genetic disease Cystic Fibrosis, species specific properties of production, inhibition and consumption of a metabolic by-product, acetate, were found. These assumptions were implemented into a mathematical chemostat model which consists of four nonlinear ordinary differential equations describing two species competing for one limiting nutrient in an open system. We derive classical chemostat results and find that our basic model supports the competitive exclusion principle, the bistability of the system as well as stable coexistence. The analytical results are illustrated by numerical simulations performed with experimentally measured parameter values. As a variant of our basic model, mimicking testing of antibiotics for therapeutic treatments in mixed cultures instead of pure ones, we consider the introduction of a lethal inhibitor, which cannot be eliminated by one of the species and is selective for the stronger competitor. We discuss our theoretical results in relation to our experimental model system and find that simulations coincide with the qualitative behavior of the experimental result in the case where the metabolic by-product serves as a second carbon source for one of the species, but not the producer.     

Models of competition in the chemostat with instantaneous and delayed nutrient recycling

Two models for competition of two populations in a chemostat environment with nutrient recycling are considered. In the first model, the recycling is instantaneous, whereas in the second, the recycling is delayed. For each model an equilibrium analysis is carried out, and persistence criteria are obtained. This paper extends the work done by Beretta et al. (1990) for a single species.Research partially supported by the Natural Sciences and Engineering Research Council of Canada, Grant NSERC A4823Research carried out at the University of Alberta while on a Canada-China Scholarly Exchange Program     

Analysis of a delayed diffusive model with Beddington-DeAngelis functional response

In this paper, a delayed diffusive phytoplankton-zooplankton model with Beddington DeAngelis functional response and toxins is investigated. Existence of equilibria of the system are solved. The global asymptotic stability of the zooplankton-free equilibrium is obtained. The local stability of the coexistent equilibrium and existence of Hopf bifurcation are discussed. In addition, the properties of the Hopf bifurcation are studied based on the center manifold and normal form theory for partial differential equations. Finally, some numerical simulations are also carried out to confirm our theoretical analysis.     

Practical coexistence of two species in the chemostat - A slow-fast characterization

We show that the chemostat model with two species having different but close break-even concentrations exhibits a slow-fast dynamics. Considering small perturbations about the dilution rate for which break-even concentrations are identical, we use the Fenichel theory to show the coexistence of species for large times. Then we determine the reduced dynamics, which is non-trivial and characterized by the slopes of the growth functions about their break-even concentrations.     

Time-delayed model of RNA interference

RNA interference (RNAi) is a fundamental cellular process that inhibits gene expression through cleavage and destruction of target mRNA. It is responsible for a number of important intracellular functions, from being the first line of immune defence against pathogens to regulating development and morphogenesis. In this paper we consider a mathematical model of RNAi with particular emphasis on time delays associated with two aspects of primed amplification: binding of siRNA to aberrant RNA, and binding of siRNA to mRNA, both of which result in the expanded production of dsRNA responsible for RNA silencing. Analytical and numerical stability analyses are performed to identify regions of stability of different steady states and to determine conditions on parameters that lead to instability. Our results suggest that while the original model without time delays exhibits a bi-stability due to the presence of a hysteresis loop, under the influence of time delays, one of the two steady states with the high (default) or small (silenced) concentration of mRNA can actually lose its stability via a Hopf bifurcation. This leads to the co-existence of a stable steady state and a stable periodic orbit, which has a profound effect on the dynamics of the system.     

The analysis and regulation for the dynamics of a temperate bacteriophage model

The purpose of this paper is to study the asymptotical behavior of a temperate bacteriophage model in chemostat, which was first proposed by Levin et al. [B.R. Levin, F.M. Stewart, L. Chao, Resource-limited growth, competition and predation: A model and experiment studies with bacteria and bacteriophage, Am. Nat. 125 (1977) 3]. Firstly, a classification for the equilibria of the model and their stability are obtained; secondly, sufficient conditions for uniform persistence are obtained; thirdly, sufficient conditions for the global asymptotic behavior are given, and simulations for the model are presented. The theoretical results show that there are more than eight cases for the classification of the model, and that the decrease (increase) of the nutrient concentration or average lytic time (flow rate) is beneficial to the survival of the sensitive cells. Both the simulated and theoretical results show that there is a possibility of switch phenomena or a periodical outburst of the phages and the lysogens, which is caused by the internal factors rather than by some external environment. Finally, the simulation and regulation of the dynamics of the model with experimental data are presented.     

Limit cycles in a chemostat model for a single species with age structure

We model an age-structured population feeding on an abiotic resource by combining the Gurtin-MacCamy [Math. Biosci. 43 (1979) 199] approach with a standard chemostat model. Limit cycles arise by Hopf bifurcations at low values of the chemostat dilution rate, even for simple maternity functions for which the original Gurtin-MacCamy model has no oscillatory solutions. We find the exact location in parameter space of the Hopf bifurcations for special cases of our model. The onset of cycling is largely independent of both the form of the resource uptake function and the shape of the maternity function.     

具有时滞和密度制约基于比率的捕食-被捕食系统的稳定性

考虑了一类具有时滞和密度制约基于比率的捕食-被捕食系统的稳定性.证明了系统在一定条件下的一致持久性,并且得到了系统正平衡位置局部稳定和全局稳定的充分条件     

Competing predators for a prey in a chemostat model with periodic nutrient input

A system of ordinary differential equations is used to model the interactions of n competing predators on a single prey population in a chemostat environment with a periodic nutrient input. In the case of one or no predators, criteria for the existence of periodic solutions are given. In the general case, conditions for all populations to persist are derived.Research is in part based on a Ph.D. thesis submitted to the Faculty of Graduate Studies, University of AlbertaResearch is partly supported by the Natural Sciences and Engineering Research Council of Canada, Grant No. NSERC A4823     

A functional response model of a predator population foraging in a patchy habitat

1. Functional response models (e.g. Holling's disc equation) that do not take the spatial distributions of prey and predators into account are likely to produce biased estimates of predation rates. 2. To investigate the consequences of ignoring prey distribution and predator aggregation, a general analytical model of a predator population occupying a patchy environment with a single species of prey is developed. 3. The model includes the density and the spatial distribution of the prey population, the aggregative response of the predators and their mutual interference. 4. The model provides explicit solutions to a number of scenarios that can be independently combined: the prey has an even, random or clumped distribution, and the predators show a convex, sigmoid, linear or no aggregative response. 5. The model is parameterized with data from an acarine predator-prey system consisting of Phytoseiulus persimis and Tetranychus urticae inhabiting greenhouse cucumbers. 6. The model fits empirical data quite well and much better than if prey and predators were assumed to be evenly distributed among patches, or if the predators were distributed independently of the prey. 7. The analyses show that if the predators do not show an aggregative response it will always be an advantage to the prey to adopt a patchy distribution. On the other hand, if the predators are capable of responding to the distribution of prey, then it will be an advantage to the prey to be evenly distributed when its density is low and switch to a more patchy distribution when its density increases. The effect of mutual interference is negligible unless predator density is very high. 8. The model shows that prey patchiness and predator aggregation in combination can change the functional response at the population level from type II to type III, indicating that these factors may contribute to stabilization of predator-prey dynamics.     

Structural and functional characterization of a stable, 4-chlorosalicylic-acid-degrading,bacterial community in a chemostat

A mixed, stable microbial community, obtained by continuous enrichment of a sediment core using 4-chlorosalicylic acid as sole source of carbon and energy, contained 10 different bacterial species, including Klebsiella pneumonia, Pseudomonas fluorescens, P. mendocina and P. cichorii. The members of the community were grown separately on various chlorinated compounds which were readily degraded.The author is with the Department of Environmental Sciences, College of Basic Sciences and Humanities, Gobind Ballav Pant University of Agriculture and Technology, Pantnagar, Naini-tal-263145, India     

一类含有无限时滞和离散时滞的Logistic系统的全局渐近稳定性

考虑了一类含有无限时滞和离散时滞的非自治Logistic系统。通过运用时滞泛函微分方程的基本原理及新的计算技巧表明系统是持久的,且系统在持久性条件下也是全局渐近稳定的。结果表明时滞是“有害的”,获得了一些新的结果并改进了一些已有的结果。     

Dynamical behavior of a hepatitis B virus transmission model with vaccination

Hepatitis B virus (HBV) infection is a globally health problem. In 2005, the WHO Western Pacific Regional Office set a goal of reducing chronic HBV infection rate to less than 2% among children five years of age by 2012, as an interim milestone towards the final goal of less than 1%. Many countries made some plans (such as free HBV vaccination program for all neonates in China now) to control the transmission HBV. We develop a model to explore the impact of vaccination and other controlling measures of HBV infection. The model has simple dynamical behavior which has a globally asymptotically stable disease-free equilibrium when the basic reproduction number R₀≤1, and a globally asymptotically stable endemic equilibrium when R₀>1. Numerical simulation results show that the vaccination is a very effective measure to control the infection and they also give some useful comments on controlling the transmission of HBV.