Spreading speeds as slowest wave speeds for cooperative systems

It is well known that in many scalar models for the spread of a fitter phenotype or species into the territory of a less fit one, the asymptotic spreading speed can be characterized as the lowest speed of a suitable family of traveling waves of the model. Despite a general belief that multi-species (vector) models have the same property, we are unaware of any proof to support this belief. The present work establishes this result for a class of multi-species model of a kind studied by Lui [Biological growth and spread modeled by systems of recursions. I: Mathematical theory, Math. Biosci. 93 (1989) 269] and generalized by the authors [Weinberger et al., Analysis of the linear conjecture for spread in cooperative models, J. Math. Biol. 45 (2002) 183; Lewis et al., Spreading speeds and the linear conjecture for two-species competition models, J. Math. Biol. 45 (2002) 219]. Lui showed the existence of a single spreading speed c(*) for all species. For the systems in the two aforementioned studies by the authors, which include related continuous-time models such as reaction-diffusion systems, as well as some standard competition models, it sometimes happens that different species spread at different rates, so that there are a slowest speed c(*) and a fastest speed c(f)(*). It is shown here that, for a large class of such multi-species systems, the slowest spreading speed c(*) is always characterized as the slowest speed of a class of traveling wave solutions.     

Slow and fast invasion waves in a model of acid-mediated tumour growth

This work is concerned with a reaction-diffusion system that has been proposed as a model to describe acid-mediated cancer invasion. More precisely, we consider the properties of travelling waves that can be supported by such a system, and show that a rich variety of wave propagation dynamics, both fast and slow, is compatible with the model. In particular, asymptotic formulae for admissible wave profiles and bounds on their wave speeds are provided.     

Anomalous spreading speeds of cooperative recursion systems

This work presents an example of a cooperative system of truncated linear recursions in which the interaction between species causes one of the species to have an anomalous spreading speed. By this we mean that this species spreads at a speed which is strictly greater than its spreading speed in isolation from the other species and the speeds at which all the other species actually spread. An ecological implication of this example is discussed in Sect. 5. Our example shows that the formula for the fastest spreading speed given in Lemma 2.3 of our paper (Weinberger et al. in J Math Biol 45:183–218, 2002) is incorrect. However, we find an extra hypothesis under which the formula for the faster spreading speed given in (Weinberger et al. in J Math Biol 45:183–218, 2002) is valid. We also show that the hypotheses of all but one of the theorems of (Weinberger et al. in J Math Biol 45:183–218, 2002) whose proofs rely on Lemma 2.3 imply this extra hypothesis, so that all but one of the theorems of (Weinberger et al. in J Math Biol 45:183–218, 2002) and all the examples given there are valid as they stand.     

A stochastic model for prostate-specific antigen levels

We introduce a continuous stochastic model for the prostate-specific antigen (PSA) levels following radiotherapy and derive solutions for the associated partial differential (Kolmogorov-Chapman) equation. The solutions describe the evolution of the time-dependent density for PSA levels which take into account an absorbing condition along the boundary and various initial conditions. We include implications for single-dose and multi-dose radiation treatment regimens and discuss parameter estimation and sensitivity issues.     

An age-and-cyclin-structured cell population model for healthy and tumoral tissues

We present a nonlinear model of the dynamics of a cell population divided into proliferative and quiescent compartments. The proliferative phase represents the complete cell cycle (G ₁−S−G ₂−M) of a population committed to divide at its end. The model is structured by the time spent by a cell in the proliferative phase, and by the amount of Cyclin D/(CDK4 or 6) complexes. Cells can transit from one compartment to the other, following transition rules which differ according to the tissue state: healthy or tumoral. The asymptotic behaviour of solutions of the nonlinear model is analysed in two cases, exhibiting tissue homeostasis or tumour exponential growth. The model is simulated and its analytic predictions are confirmed numerically.     

Chemotactic collapse for the Keller-Segel model

 This work is concerned with the system (S) {u _t =Δu − χ∇ (u∇v) for x∈Ω, t>0Γ v _t =Δv+(u−1) for x∈Ω, t>0 where Γ, χ are positive constants and Ω is a bounded and smooth open set in ℝ². On the boundary ∂Ω, we impose no-flux conditions: (N) ∂u∂n =∂v∂n =0 for x∈∂ Ω, t>0 Problem (S), (N) is a classical model to describe chemotaxis corresponding to a species of concentration u(x, t) which tends to aggregate towards high concentrations of a chemical that the species releases. When completed with suitable initial values at t=0 for u(x, t), v(x, t), the problem under consideration is known to be well posed, locally in time. By means of matched asymptotic expansions techniques, we show here that there exist radial solutions exhibiting chemotactic collapse. By this we mean that u(r, t) →Aδ(y) as t→T for some T<∞, where A is the total concentration of the species. Received 9 March 1995; received in revised form 25 December 1995     

Mathematical analysis of a cholera model with public health interventions

Cholera, an acute gastro-intestinal infection and a waterborne disease continues to emerge in developing countries and remains an important global health challenge. We formulate a mathematical model that captures some essential dynamics of cholera transmission to study the impact of public health educational campaigns, vaccination and treatment as control strategies in curtailing the disease. The education-induced, vaccination-induced and treatment-induced reproductive numbers R(E), R(V), R(T) respectively and the combined reproductive number R(C) are compared with the basic reproduction number R(0) to assess the possible community benefits of these control measures. A Lyapunov functional approach is also used to analyse the stability of the equilibrium points. We perform sensitivity analysis on the key parameters that drive the disease dynamics in order to determine their relative importance to disease transmission and prevalence. Graphical representations are provided to qualitatively support the analytical results.     

A complex mathematical model of the human menstrual cycle

Despite the fact that more than 100 million women worldwide use birth control pills and that half of the world's population is concerned, the menstrual cycle has so far received comparatively little attention in the field of mathematical modeling. The term menstrual cycle comprises the processes of the control system in the female body that, under healthy circumstances, lead to ovulation at regular intervals, thus making reproduction possible. If this is not the case or ovulation is not desired, the question arises how this control system can be influenced, for example, by hormonal treatments. In order to be able to cover a vast range of external manipulations, the mathematical model must comprise the main components where the processes belonging to the menstrual cycle occur, as well as their interrelations. A system of differential equations serves as the mathematical model, describing the dynamics of hormones, enzymes, receptors, and follicular phases. Since the processes take place in different parts of the body and influence each other with a certain delay, passing over to delay differential equations is deemed a reasonable step. The pulsatile release of the gonadotropin-releasing hormone (GnRH) is controlled by a complex neural network. We choose to model the pulse time points of this GnRH pulse generator by a stochastic process. Focus in this paper is on the model development. This rather elaborate mathematical model is the basis for a detailed analysis and could be helpful for possible drug design.     

A theoretical study on mathematical modelling of an infectious disease with application of optimal control

In this paper, we propose and analyze an epidemic problem which can be controlled by vaccination as well as treatment. In the first part of our analysis we study the dynamical behavior of the system with fixed control for both vaccination and treatment. Basic reproduction number is obtained in all possible cases and it is observed that the simultaneous use of vaccination and treatment control is the most favorable case to prevent the disease from being epidemic. In the second part, we take the controls as time dependent and obtain the optimal control strategy to minimize both the infected populations and the associated costs. All the analytical results are verified by simulation works. Some important conclusions are given at the end of the paper.     

Singular perturbation analysis of cAMP signalling in Dictyostelium discoideum aggregates

In this paper, we use singular perturbation methods to study the structure of travelling waves for some reaction-diffusion models obtained from the Martiel-Goldbeter and Goldbeter-Segel's models of cAMP signalling in Dictyostelium discoideum. As a consequence, we derive analytic formulae for quantities like wave speed, maximum concentration and other magnitudes in terms of the different biochemical constants that appear in the model.     

Analysis of a competitive prey–predator system with a prey refuge

Gauss's competitive exclusive principle states that two competing species having analogous environment cannot usually occupy the same space at a time but in order to exploit their common environment in a different manner, they can co-exist only when they are active in different times. On the other hand, several studies on predators in various natural and laboratory situations have shown that competitive coexistence can result from predation in a way by resisting any one prey species from becoming sufficiently abundant to outcompete other species such that the predator makes the coexistence possible. It has also been shown that the use of refuges by a fraction of the prey population exerts a stabilizing effect in the interacting population dynamics. Further, the field surveys in the Sundarban mangrove ecosystem reveal that two detritivorous fishes, viz. Liza parsia and Liza tade (prey population) coexist in nature with the presence of the predator fish population, viz. Lates calcarifer by using refuges.     

Persistence in reaction diffusion models with weak allee effect

We study the positive steady state distributions and dynamical behavior of reaction-diffusion equation with weak Allee effect type growth, in which the growth rate per capita is not monotonic as in logistic type, and the habitat is assumed to be a heterogeneous bounded region. The existence of multiple steady states is shown, and the global bifurcation diagrams are obtained. Results are applied to a reaction-diffusion model with type II functional response, and also a model with density-dependent diffusion of animal aggregation. J. S. is partially supported by United States NSF grants DMS-0314736 and EF-0436318, College of William and Mary summer grants, and a grant from Science Council of Heilongjiang Province, China.     

Spreading speeds of spatially periodic integro-difference models for populations with nonmonotone recruitment functions

An idea used by Thieme (J. Math. Biol. 8, 173-187, 1979) is extended to show that a class of integro-difference models for a periodically varying habitat has a spreading speed and a formula for it, even when the recruitment function R(u, x) is not nondecreasing in u, so that overcompensation occurs. Numerical simulations illustrate the behavior of solutions of the recursion whose initial values vanish outside a bounded set.     

Stability analysis for HIV infection delay model with protease inhibitor

In this article, we considered a model of HIV-1 infection with a protease inhibitor therapy and three delays. The frequency of the bifurcating periodic solution as well as the threshold value is approximated numerically using realistic parameter. The estimated threshold value is realistic and the frequency of the oscillations is consistent with that of the observed viral blips.     

Oscillations in a patchy environment disease model

For a single patch SIRS model with a period of immunity of fixed length, recruitment-death demographics, disease related deaths and mass action incidence, the basic reproduction number R(0) is identified. It is shown that the disease-free equilibrium is globally asymptotically stable if R(0)<1. For R(0)>1, local stability of the endemic equilibrium and Hopf bifurcation analysis about this equilibrium are carried out. Moreover, a practical numerical approach to locate the bifurcation values for a characteristic equation with delay-dependent coefficients is provided. For a two patch SIRS model with travel, it is shown that there are several threshold quantities determining its dynamic behavior and that travel can reduce oscillations in both patches; travel may enhance oscillations in both patches; or travel can switch oscillations from one patch to another.     

Sustainability and optimal control of an exploited prey predator system through provision of alternative food to predator

In the present paper, we develop a simple two species prey-predator model in which the predator is partially coupled with alternative prey. The aim is to study the consequences of providing additional food to the predator as well as the effects of harvesting efforts applied to both the species. It is observed that the provision of alternative food to predator is not always beneficial to the system. A complete picture of the long run dynamics of the system is discussed based on the effort pair as control parameters. Optimal augmentations of prey and predator biomass at final time have been investigated by optimal control theory. Also the short and large time effects of the application of optimal control have been discussed. Finally, some numerical illustrations are given to verify our analytical results with the help of different sets of parameters.     

Optimal control strategies and cost-effectiveness analysis of a malaria model

The aim of this paper is to investigate the effectiveness and cost-effectiveness of three malaria preventive measures (use of treated bednets, spray of insecticides and a possible treatment of infective humans that blocks transmission to mosquitoes). For this, we consider a mathematical model for the transmission dynamics of the disease that includes these measures. We first consider the constant control parameters’ case, we calculate the basic reproduction number and investigate the existence and stability of equilibria; the model is found to exhibit backward bifurcation. We then assess the relative impact of each of the constant control parameters measures by calculating the sensitivity index of the basic reproductive number to the model's parameters. In the time-dependent constant control case, we use Pontryagin's Maximum Principle to derive necessary conditions for the optimal control of the disease. We also calculate the Infection Averted Ratio (IAR) and the Incremental Cost-Effectiveness Ratio (ICER) to investigate the cost-effectiveness of all possible combinations of the three control measures. One of our findings is that the most cost-effective strategy for malaria control, is the combination of the spray of insecticides and treatment of infective individuals. This strategy requires a 100% effort in both treatment (for 20 days) and spray of insecticides (for 57 days). In practice, this will be extremely difficult, if not impossible to achieve. The second most cost-effective strategy which consists of a 100% use of treated bednets and 87% treatment of infective individuals for 42 and 100 days, respectively, is sustainable and therefore preferable.     

Mum, why do you keep on growing?Impacts of environmental variability on optimal growth and reproduction allocation strategies of annual plants

In their 1990 paper Optimal reproductive efforts and the timing of reproduction of annual plants in randomly varying environments, Amir and Cohen considered stochastic environments consisting of i.i.d. sequences in an optimal allocation discrete-time model. We suppose here that the sequence of environmental factors is more generally described by a Markov chain. Moreover, we discuss the connection between the time interval of the discrete-time dynamic model and the ability of the plant to rebuild completely its vegetative body (from reserves). We formulate a stochastic optimization problem covering the so-called linear and logarithmic fitness (corresponding to variation within and between years), which yields optimal strategies. For ``linear maximizers', we analyse how optimal strategies depend upon the environmental variability type: constant, random stationary, random i.i.d., random monotonous. We provide general patterns in terms of targets and thresholds, including both determinate and indeterminate growth. We also provide a partial result on the comparison between ``linear maximizers' and ``log maximizers'. Numerical simulations are provided, allowing to give a hint at the effect of different mathematical assumptions.     

Optimal isolation control strategies and cost-effectiveness analysis of a two-strain avian influenza model

The most important and effective measures against disease outbreaks in the absence of valid medicines or vaccine are quarantine and isolation strategies. In this paper optimal control theory is applied to a system of ordinary differential equation describing a two-strain avian influenza transmission via the Pontryagin's Maximum Principle. To this end, a pair of control variables representing the isolation strategies for individuals with avian and mutant strains were incorporated into the transmission model. The infection averted ratio (IAR) and the incremental cost-effectiveness ratio (ICER) were calculated to investigate the cost-effectiveness of all possible combinations of the control strategies. The simulation results show that the implementation of the combination strategy during the epidemic is the most cost-effective strategy for avian influenza transmission. This is followed by the control strategy involving isolation of individuals with the mutant strain. Also observed was the fact that low mutating and more virulent virus results in an increased control effort of isolating individuals with the avian strain; and high mutating with more virulent virus results in increased efforts in isolating individuals with the mutant strain.     

Mathematical analysis of the global dynamics of a model for HTLV-I infection and ATL progression

Mathematical analysis is carried out that completely determines the global dynamics of a mathematical model for the transmission of human T-cell lymphotropic virus I (HTLV-I) infection and the development of adult T-cell leukemia (ATL). HTLV-I infection of healthy CD4(+) T cells takes place through cell-to-cell contact with infected T cells. The infected T cells can remain latent and harbor virus for several years before virus production occurs. Actively infected T cells can infect other T cells and can convert to ATL cells, whose growth is assumed to follow a classical logistic growth function. Our analysis establishes that the global dynamics of T cells are completely determined by a basic reproduction number R(0). If R(0)< or =1, infected T cells always die out. If R(0)>1, HTLV-I infection becomes chronic, and a unique endemic equilibrium is globally stable in the interior of the feasible region. We also show that the equilibrium level of ATL-cell proliferation is higher when the HTLV-I infection of T cells is chronic than when it is acute.