Joint effects of mitosis and intracellular delay on viral dynamics: two-parameter bifurcation analysis

To understand joint effects of logistic growth in target cells and intracellular delay on viral dynamics in vivo, we carry out two-parameter bifurcation analysis of an in-host model that describes infections of many viruses including HIV-I, HBV and HTLV-I. The bifurcation parameters are the mitosis rate r of the target cells and an intracellular delay τ in the incidence of viral infection. We describe the stability region of the chronic-infection equilibrium E* in the two-dimensional (r, τ) parameter space, as well as the global Hopf bifurcation curves as each of τ and r varies. Our analysis shows that, while both τ and r can destabilize E* and cause Hopf bifurcations, they do behave differently. The intracellular delay τ can cause Hopf bifurcations only when r is positive and sufficiently large, while r can cause Hopf bifurcations even when τ = 0. Intracellular delay τ can cause stability switches in E* while r does not.     

Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations I. Global organization of bistable periodic solutions

 The Hodgkin-Huxley equations (HH) are parameterized by a number of parameters and shows a variety of qualitatively different behaviors depending on the parameter values. We explored the dynamics of the HH for a wide range of parameter values in the multiple-parameter space, that is, we examined the global structure of bifurcations of the HH. Results are summarized in various two-parameter bifurcation diagrams with I _ext (externally applied DC current) as the abscissa and one of the other parameters as the ordinate. In each diagram, the parameter plane was divided into several regions according to the qualitative behavior of the equations. In particular, we focused on periodic solutions emerging via Hopf bifurcations and identified parameter regions in which either two stable periodic solutions with different amplitudes and periods and a stable equilibrium point or two stable periodic solutions coexist. Global analysis of the bifurcation structure suggested that generation of these regions is associated with degenerate Hopf bifurcations. Received: 23 April 1999 / Accepted in revised form: 24 September 1999     

Periodicity in an epidemic model with a generalized non-linear incidence

We develop and analyze a simple SIV epidemic model including susceptible, infected and perfectly vaccinated classes, with a generalized non-linear incidence rate subject only to a few general conditions. These conditions are satisfied by many models appearing in the literature. The detailed dynamics analysis of the model, using the Poincaré index theory, shows that non-linearity of the incidence rate leads to vital dynamics, such as bistability and periodicity, without seasonal forcing or being cyclic. Furthermore, it is shown that the basic reproductive number is independent of the functional form of the non-linear incidence rate. Under certain, well-defined conditions, the model undergoes a Hopf bifurcation. Using the normal form of the model, the first Lyapunov coefficient is computed to determine the various types of Hopf bifurcation the model undergoes. These general results are applied to two examples: unbounded and saturated contact rates; in both cases, forward or backward Hopf bifurcations occur for two distinct values of the contact parameter. It is also shown that the model may undergo a subcritical Hopf bifurcation leading to the appearance of two concentric limit cycles. The results are illustrated by numerical simulations with realistic model parameters estimated for some infectious diseases of childhood.     

A disease transmission model in a nonconstant population

A general SIRS disease transmission model is formulated under assumptions that the size of the population varies, the incidence rate is nonlinear, and the recovered (removed) class may also be directly reinfected. For a class of incidence functions it is shown that the model has no periodic solutions. By contrast, for a particular incidence function, a combination of analytical and numerical techniques are used to show that (for some parameters) periodic solutions can arise through homoclinic loops or saddle connections and disappear through Hopf bifurcations.Supported in part by NSERC grant A-8965, the University of Victoria Committee on Faculty Research & Travel, and the Institute for Mathematics and its Applications, Minneapolis, MN, with funds provided by NSF     

一类具有时滞的传染病模型的稳定性分析

研究了一类具有时滞的传染病生物模型．首先研究了该模型的线性稳定性,并给出了一列Hopf分支值,然后利用中心流形定理和正规型方法,给出了确定分支周期解的分支方向与稳定性的计算公式．     

The effect of predator limitation on the dynamics of simple food chains

We investigate the influence of competition between predators on the dynamics of bitrophic predator–prey systems and of tritrophic food chains. Competition between predators is implemented either as interference competition, or as a density-dependent mortality rate. With interference competition, the paradox of enrichment is reduced or completely suppressed, but otherwise, the dynamical behavior of the systems is not fundamentally different from that of the Rosenzweig–MacArthur model, which contains no predator competition and shows only continuous transitions between fixed points or periodic oscillations. In contrast, with density-dependent predator mortality, the system shows a surprisingly rich dynamical behavior. In particular, decreasing the density regulation of the predator can induce catastrophic shifts from a stable fixed point to a large oscillation where the predator chases the prey through a cycle that brings both species close to the threshold of extinction. Other catastrophic bifurcations, such as subcritical Hopf bifurcations and saddle-node bifurcations of limit cycles, do also occur. In tritrophic food chains, we find again that fixed points in the model with predator interference become unstable only through Hopf bifurcations, which can also be subcritical, in contrast to the bitrophic situation. The model with a density limitation shows again catastrophic destabilization of fixed points and various nonlocal bifurcations. In addition, chaos occurs for both models in appropriate parameter ranges.     

Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations II. Singularity theoretic approach and highly degenerate bifurcations

In the Hodgkin-Huxley equations (HH), we have identified the parameter regions in which either two stable periodic solutions with different amplitudes and periods and an equilibrium point or two stable periodic solutions coexist. The global structure of bifurcations in the multiple-parameter space in the HH suggested that the bistabilities of the periodic solutions are associated with the degenerate Hopf bifurcation points by which several qualitatively different behaviors are organized. In this paper, we clarify this by analyzing the details of the degenerate Hopf bifurcations using the singularity theory approach which deals with local bifurcations near a highly degenerate fixed point. Received: 23 April 1999 / Accepted in revised form: 24 September 1999     

Bifurcations in a white-blood-cell production model

We study the dynamics of a model of white-blood-cell (WBC) production. The model consists of two compartmental differential equations with two discrete delays. We show that from normal to pathological parameter values, the system undergoes supercritical Hopf bifurcations and saddle-node bifurcations of limit cycles. We characterize the steady states of the system and perform a bifurcation analysis. Our results indicate that an increase in apoptosis rate of either hematopoietic stem cells or WBC precursors induces a Hopf bifurcation and an oscillatory regime takes place. These oscillations are seen in some hematological diseases.     

Hopf bifurcation and bursting synchronization in an excitable systems with chemical delayed coupling

In this paper we consider the Hopf bifurcation and synchronization in the two coupled Hindmarsh–Rose excitable systems with chemical coupling and time-delay. We surveyed the conditions for Hopf bifurcations by means of dynamical bifurcation analysis and numerical simulation. The results show that the coupled excitable systems with no delay have supercritical Hopf bifurcation, while the delayed system undergoes Hopf bifurcations at critical time delays when coupling strength lies in a particular region. We also investigated the effect of the delay on the transition of bursting synchronization in the coupled system. The results are helpful for us to better understand the dynamical properties of excitable systems and the biological mechanism of information encoding and cognitive activity.     

具有三时滞Lotka-Volterra互惠系统的Hopf分支

建立了具有三个时滞的Lotka-Volterra互惠系统;获得了正平衡点和Hopf分支存在的条件等;并对所获得的结果进行了数值模拟．     

Bifurcation analysis of the Nowak-Bangham model in CTL dynamics

This paper investigates the local bifurcations of a CTL response model published by Nowak and Bangham [M.A. Nowak, C.R.M. Bangham, Population dynamics of immune responses to persistent viruses, Science 272 (1996) 74]. The Nowak-Bangham model can have three equilibria depending on the basic reproduction number, and generates a Hopf bifurcation through two bifurcations of equilibria. The main result shows a sufficient condition for the interior equilibrium to have a unique bifurcation point at which a simple Hopf bifurcation occurs. For this proof, some new techniques are developed in order to apply the method established by Liu [W.M. Liu, Criterion of Hopf bifurcations without using eigenvalues, J. Math. Anal. Appl. 182 (1) (1994) 250]. In addition, to demonstrate the result obtained theoretically, some bifurcation diagrams are presented with numerical examples.     

The computer-aided bifurcation analysis of predator-prey models

Mathematical analysis of dynamical systems can often benefit from accompanying numerical computations. This is particularly true if one has software (e.g. AUTO [6, 7]) capable of providing an automatic bifurcation analysis of such systems. Computer programs of this type now exist. We describe the application of such software to a predator-prey model. Phenomena that arise in this analysis include stationary bifurcations, limit points, Hopf bifurcations and secondary periodic bifurcations. A two-parameter numerical analysis leads quite naturally to the detection of higher order singularities.Supported in part by NSERC Canada (#4274) and FCAC Québec (#EQ1438)     

A comparison of two predator-prey models with Holling's type I functional response

In this paper, we analyze a laissez-faire predator-prey model and a Leslie-type predator-prey model with type I functional responses. We study the stability of the equilibrium where the predator and prey coexist by both performing a linearized stability analysis and by constructing a Lyapunov function. For the Leslie-type model, we use a generalized Jacobian to determine how eigenvalues jump at the corner of the functional response. We show, numerically, that our two models can both possess two limit cycles that surround a stable equilibrium and that these cycles arise through global cyclic-fold bifurcations. The Leslie-type model may also exhibit super-critical and discontinuous Hopf bifurcations. We then present and analyze a new functional response, built around the arctangent, that smoothes the sharp corner in a type I functional response. For this new functional response, both models undergo Hopf, cyclic-fold, and Bautin bifurcations. We use our analyses to characterize predator-prey systems that may exhibit bistability.     

Biomass and biodiversity in species-rich tritrophic communities

We consider a tritrophic system with one basal and one top species and a large number of primary consumers, and derive upper and lower bounds for the total biomass of the middle trophic level. These estimates do not depend on dynamical regime, holding for fixed point, periodic, or chaotic dynamics. We have two kinds of estimates, depending on whether the predator abundance is zero. All these results are uniform in a self-limitation parameter, which regulates prey diversity in the system. For strong self-limitation, diversity is large; for weak self-limitation, it is small. Diversity depends on the variance of species’ parameter values. The larger this variance, the lower the diversity, and vice versa. Moreover, variation in the parameters of the Holling type II functional response changes the bifurcation character, with the equilibrium state with nonzero predator abundance losing stability. If that variation is small then the bifurcation can lead to oscillations (the Hopf bifurcation). Under certain conditions, there exists a supercritical Hopf bifurcation. We then find a connection between diversity and Hopf bifurcations. We also show that the system exhibits top-down regulation and a hump-shaped diversity-productivity curve.We then extend the model by allowing species to experience self-regulation. For this extended model, explicit estimates of prey diversity are obtained. We study the dynamics of this system and find the following. First, diversity and system dynamics crucially depend on variation in species parameters. We show that under certain conditions, the system undergoes a supercritical Hopf bifurcation. We also establish a connection between diversity and Hopf bifurcations. For strong self-limitation, diversity is large and complex dynamics are absent. For weak self-limitation, diversity is small and the equilibrium with non-zero predator abundance is unstable.     

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.     

Asymmetry in the Presence of Migration Stabilizes Multistrain Disease Outbreaks

We study the effect of migration between coupled populations, or patches, on the stability properties of multistrain disease dynamics. The epidemic model used in this work displays a Hopf bifurcation to oscillations in a single, well-mixed population. It is shown numerically that migration between two non-identical patches stabilizes the endemic steady state, delaying the onset of large amplitude outbreaks and reducing the total number of infections. This result is motivated by analyzing generic Hopf bifurcations with different frequencies and with diffusive coupling between them. Stabilization of the steady state is again seen, indicating that our observation in the full multistrain model is based on qualitative characteristics of the dynamics rather than on details of the disease model.     

Self-sustained Oscillations in Blood Flow Through a Honeycomb Capillary Network

Numerical simulations of unsteady blood flow through a honeycomb network originating at multiple inlets and terminating at multiple outlets are presented and discussed under the assumption that blood behaves as a continuum with variable constitution. Unlike a tree network, the honeycomb network exhibits both diverging and converging bifurcations between branching capillary segments. Numerical results based on a finite difference method demonstrate that as in the case of tree networks considered in previous studies, the cell partitioning law at diverging bifurcations is an important parameter in both steady and unsteady flow. Specifically, a steady flow may spontaneously develop self-sustained oscillations at critical conditions by way of a Hopf bifurcation. Contrary to tree-like networks comprised entirely of diverging bifurcations, the critical parameters for instability in honeycomb networks depend weakly on the system size. The blockage of one or more network segments due to the presence of large cells or the occurrence of capillary constriction may cause flow reversal or trigger a transition to unsteady flow.     

A mathematical model of cell-to-cell spread of HIV-1 that includes a time delay

 We consider a two-dimensional model of cell-to-cell spread of HIV-1 in tissue cultures, assuming that infection is spread directly from infected cells to healthy cells and neglecting the effects of free virus. The intracellular incubation period is modeled by a gamma distribution and the model is a system of two differential equations with distributed delay, which includes the differential equations model with a discrete delay and the ordinary differential equations model as special cases. We study the stability in all three types of models. It is shown that the ODE model is globally stable while both delay models exhibit Hopf bifurcations by using the (average) delay as a bifurcation parameter. The results indicate that, differing from the cell-to-free virus spread models, the cell-to-cell spread models can produce infective oscillations in typical tissue culture parameter regimes and the latently infected cells are instrumental in sustaining the infection. Our delayed cell-to-cell models may be applicable to study other types of viral infections such as human T-cell leukaemia virus type 1 (HTLV-1). Received: 18 November 2000 / Published online: 28 February 2003 RID="*" ID="*" Research was partially supported by the NSERC and MITACS of Canada and a start-up fund from the College of Arts and Sciences at the University of Miami. On leave from Dalhousie University, Halifax, Nova Scotia, Canada. Current address: Department of Mathematics, Clarke College, Dubuque, Iowa 52001, USA Key words or phrases: HIV-1 – Cell-to-cell spread – Time delay – Stability – Hopf bifurcation – Periodicity     

Critical slowing down as an indicator of transitions in two-species models

Transitions in ecological systems often occur without apparent warning, and may represent shifts between alternative persistent states. Decreasing ecological resilience (the size of the basin of attraction around a stable state) can signal an impending transition, but this effect is difficult to measure in practice. Recent research has suggested that a decreasing rate of recovery from small perturbations (critical slowing down) is a good indicator of ecological resilience. Here we use analytical techniques to draw general conclusions about the conditions under which critical slowing down provides an early indicator of transitions in two-species predator-prey and competition models. The models exhibit three types of transition: the predator-prey model has a Hopf bifurcation and a transcritical bifurcation, and the competition model has two saddle-node bifurcations (in which case the system exhibits hysteresis) or two transcritical bifurcations, depending on the parameterisation. We find that critical slowing down is an earlier indicator of the Hopf bifurcation in predator-prey models in which prey are regulated by predation rather than by intrinsic density-dependent effects and an earlier indicator of transitions in competition models in which the dynamics of the rare species operate on slower timescales than the dynamics of the common species. These results lead directly to predictions for more complex multi-species systems, which can be tested using simulation models or real ecosystems.