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.     

Impulsive control strategies in biological control of pesticide

By presenting and analyzing the pest-predator model under insecticides used impulsively, two impulsive strategies in biological control are put forward. The first strategy: the pulse period is fixed, but the proportional constant E(1) changes, which represents the fraction of pests killed by applying insecticide. For this scheme, two thresholds, E(1)(**) and E(1)(*) for E(1) are obtained. If E(1)>or=E(1)(*), both the pest and predator (natural enemies) populations go to extinction. If E(1)(**)相似文献   

Linked selected and neutral loci in heterogeneous environments

We analyze a system of ordinary differential equations modeling haplotype frequencies at a physically linked pair of loci, one selected and one neutral, in a population consisting of two demes with divergent selection regimes. The system is singularly perturbed, with the migration rate m between the demes serving as a small parameter. We use geometric singular perturbation theory to show that when m is sufficiently small, each solution not initially fixed for the same selected allele in both demes approaches one of a 1-dimensional continuum of equilibria. We then obtain asymptotic expansions of the solutions and show their validity on arbitrarily long finite time intervals. From these expansions we obtain formulas for the transient dynamics of F _ST (a measure of population structure) at both loci, as well as for the rate of genotyping error if the allelic state at the selected locus is inferred from that at the neutral (marker) locus. We examine two cases in detail, one modeling two populations in secondary contact after a period of evolution in allopatry, and the other modeling the origination and spread of a resistance allele.Electronic supplementary material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00285-006-0038-6 and is accessible for authorized users.     

Multiple limit cycles in the chemostat with variable yield

The global asymptotic behavior of solutions of the variable yield model is determined. The model generalizes the classical Monod model and it assumes that the yield is an increasing function of the nutrient concentration. In contrast to the Monod model, it is demonstrated that the variable yield model exhibits sustained oscillations. Moreover, it is shown that the variable yield model may undergo a subcritical Hopf bifurcation and feature at least two distinct limit cycles. Implications for the coexistence of competing populations are discussed.     

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.     

Competition between plasmid-bearing and plasmid-free organisms in a chemostat with nutrient recycling and an inhibitor

The asymptotic behavior of solutions of a model for competition between plasmid-bearing and plasmid-free organisms in the chemostat with two distributed delays and an external inhibitor is considered. The model presents a refinement of the one considered by Lu and Hadeler [Z. Lu, K.P. Hadeler, Model of plasmid-bearing plasmid-free competition in the chemostat with nutrient recycling and an inhibitor, Math. Biosci. 167 (2000) p. 177]. The delays model the fact that the nutrient is partially recycled after the death of the biomass by bacterial decomposition. Furthermore, it is assumed that there is inter-specific competition between the plasmid-bearing and plasmid-free organisms as well as intra-specific competition within each population. Conditions for boundedness of solutions and existence of non-negative equilibrium are given. Analysis of the extinction of the organisms, including plasmid-bearing and plasmid-free organisms, and the uniform persistence of the system are also carried out. By constructing appropriate Liapunov-like functionals, some sufficient conditions of global attractivity to the extinction equilibria are obtained and the combined effects of the delays and the inhibitor are studied.     

Dynamic multipopulation and density dependent evolutionary games related to replicator dynamics. A metasimplex concept

This paper contains the basic extensions of classical evolutionary games (multipopulation and density dependent models). It is shown that classical bimatrix approach is inconsistent with other approaches because it does not depend on proportion between populations. The main conclusion is that interspecific proportion parameter is important and must be considered in multipopulation models. The paper provides a synthesis of both extensions (a metasimplex concept) which solves the problem intrinsic in the bimatrix model. It allows us to model interactions among any number of subpopulations including density dependence effects. We prove that all modern approaches to evolutionary games are closely related. All evolutionary models (except classical bimatrix approaches) can be reduced to a single population general model by a simple change of variables. Differences between classic bimatrix evolutionary games and a new model which is dependent on interspecific proportion are shown by examples.     

Measurement and modeling of peripheral detection and discrimination thresholds

This report describes experimental measurements of threshold contrasts as a function of the angle to the visual axis (peripheral threshold contrasts). The visual tasks consist in detection (perception of presence) and discrimination (perception of a form feature) of simple visual signs during a fixation period realistic observing conditions being chosen. Proceeding from the experimental findings a model for forecasting off-axis threshold contrast functions on different visual conditions is developed based upon spatial frequency filters. Further with the aid of a known model visibility fields are calculated.     

Plasmid-bearing,plasmid-free organisms competing for two complementary nutrients in a chemostat

A model of competition for two complementary nutrients between plasmid-bearing and plasmid-free organisms in a chemostat is proposed. A rigorous mathematical analysis of the global asymptotic behavior of the model is presented. The work extends the model of competition for a single-limited nutrient studied by Stephanopoulos and Lapidus [Chem. Engng. Sci. 443 (1988) 49] and Hsu, Waltman and Wolkowicz [J. Math. Biol. 32 (1994) 731].     

Estimation of effective vaccination rate: pertussis in New Zealand as a case study

In some cases vaccination is unreliable. For example vaccination against pertussis has comparatively high level of primary and secondary failures. To evaluate efficiency of vaccination we introduce the idea of effective vaccination rate and suggest an approach to estimate it. We consider pertussis in New Zealand as a case study. The results indicate that the level of immunity failure for pertussis is considerably higher than was anticipated.     

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.     

An impulsively controlled pest management model with n predator species and a common prey

This paper investigates the dynamics of a competitive single-prey n-predators model of integrated pest management, which is subject to periodic and impulsive controls, from the viewpoint of finding sufficient conditions for the extinction of prey and for prey and predator permanence. The per capita death rates of prey due to predation are given in abstract, unspecified forms, which encompass large classes of death rates arising from usual predator functional responses, both prey-dependent and predator-dependent. The stability and permanence conditions are then expressed as balance conditions between the cumulative death rate of prey in a period, due to predation from all predator species and to the use of control, and to the cumulative birth rate of prey in the same amount of time. These results are then specialized for the case of prey-dependent functional responses, their biological significance being also discussed.     

Mathematical modeling and qualitative analysis of insulin therapies

Several insulin therapies are widely in clinical use with the basic strategy that mimics insulin secretion in a normal glucose-insulin endocrine metabolic regulatory system. In this paper, we model the insulin therapies using a delay differential equation model. We study the dynamics of the model both qualitatively and quantitatively. The analytical results show the existence and uniqueness of a stable periodic solution that corresponds to ultradian insulin secretion oscillations. Numerically we simulate the insulin administration based on our model. The numerical simulation results are in agreement with findings of clinical studies.     

Cluster formation for multi-strain infections with cross-immunity

Many infectious diseases exist in several pathogenic variants, or strains, which interact via cross-immunity. It is observed that strains tend to self-organise into groups, or clusters. The aim of this paper is to investigate cluster formation. Computations demonstrate that clustering is independent of the model used, and is an intrinsic feature of the strain system itself. We observe that an ordered strain system, if it is sufficiently complex, admits several cluster structures of different types. Appearance of a particular cluster structure depends on levels of cross-immunity and, in some cases, on initial conditions. Clusters, once formed, are stable, and behave remarkably regularly (in contrast to the generally chaotic behaviour of the strains themselves). In general, clustering is a type of self-organisation having many features in common with pattern formation.     

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.     

Rabbits killing birds revisited

We formulate and study a three-species population model consisting of an endemic prey (bird), an alien prey (rabbit) and an alien predator (cat). Our model overcomes several model construction problems in existing models. Moreover, our model generates richer, more reasonable and realistic dynamics. We explore the possible control strategies to save or restore the bird by controlling or eliminating the rabbit or the cat when the bird is endangered. We confirm the existence of the hyperpredation phenomenon, which is a big potential threat to most endemic prey. Specifically, we show that, in an endemic prey-alien prey-alien predator system, eradication of introduced predators such as the cat alone is not always the best solution to protect endemic insular prey since predator control may fail to protect the indigenous prey when the control of the introduced prey is not carried out simultaneously.     

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.     

The asymptotic behavior of solutions of the buffered bistable system

In this paper, we study a model for calcium buffering with bistable nonlinearity. We present some results on the stability of equilibrium states and show that there exists a threshold phenomenon in our model. In comparing with the model without buffers, we see that stationary buffers cannot destroy the asymptotic stability of the associated equilibrium states and the threshold phenomenon. Moreover, we also investigate the propagation property of solutions with initial data being a disturbance of one of the stable states which is confined to a half-line. We show that the more stable state will eventually dominate the whole dynamics and that the speed of this propagation (or invading process) is positive.     

Introduction to special issue on RNA

In this introduction to the special issue on RNA, we provide a brief overview of some of the novel and exciting biological discoveries concerning diverse roles played by RNA, and subsequently we give a rapid summary of some algorithmic aspects of RNA structure and alignment. Each of the contributions to this special issue is briefly described. This work is funded in part by the National Science Foundation under Grant No. DBI-0543506.