Coevolutionary interactions between a haploid species and a diploid species

We investigate a general model describing coevolutionary interaction between a haploid population and a diploid population, each with two alleles at a single locus. Both species are allowed to evolve, with the fitness of the genotypes of each species assumed to depend linearly on the frequencies of the genotypes of the other species. We explore the resulting outcomes of these interactions, in particular determining the location of equilibria under various conditions. The coevolution here is much more complex than that between two haploid populations and allows for the possibility of two polymorphic equilibria. To allow for further analysis, we construct a semi-symmetric model. The variety of outcomes possible even in this second model provides support for the geographic mosaic theory of coevolution by suggesting the possibility of small local populations coevolving to very different outcomes, leading to a shifting geographic mosaic as neighboring populations interact with each other through migration.     

Two SIS epidemiologic models with delays

 The SIS epidemiologic models have a delay corresponding to the infectious period, and disease-related deaths, so that the population size is variable. The population dynamics structures are either logistic or recruitment with natural deaths. Here the thresholds and equilibria are determined, and stabilities are examined. In a similar SIS model with exponential population dynamics, the delay destabilized the endemic equilibrium and led to periodic solutions. In the model with logistic dynamics, periodic solutions in the infectious fraction can occur as the population approaches extinction for a small set of parameter values. Received: 10 January 1997 / 18 November 1997     

The influence of drug treatment on the maintenance of schistosome genetic diversity

Drug treatment of patients with schistosomiasis may select for drug-resistant parasites. In this article, we formulate a deterministic model with multiple strains of schistosomes (helminth parasites with a two-host life cycles) in order to explore the role of drug treatment in the maintenance of a polymorphism of parasite strains that differ in their resistance levels. The basic reproductive numbers for all strains are computed, and are shown to determine the stabilities of equilibria of the model and consequently the distribution of parasite phenotypes with different levels of drug tolerance. Analysis of our model shows that the likelihood that resistant strains will increase in frequency depends on the interplay between their relative fitness, the cost of resistance, and the degree of selection pressure exerted by the drug treatments.     

Structured population on two patches: modeling dispersal and delay

We derive from the age-structured model a system of delay differential equations to describe the interaction of spatial dispersal (over two patches) and time delay (arising from the maturation period). Our model analysis shows that varying the immature death rate can alter the behavior of the homogeneous equilibria, leading to transient oscillations around an intermediate equilibrium and complicated dynamics (in the form of the coexistence of possibly stable synchronized periodic oscillations and unstable phase-locked oscillations) near the largest equilibrium.     

Brucellosis, botflies, and brainworms: the impact of edge habitats on pathogen transmission and species extinction

Ecological interactions between species that prefer different habitat types but come into contact in edge regions at the interfaces between habitat types are modeled via reaction-diffusion systems. The primary sort of interaction described by the models is competition mediated by pathogen transmission. The models are somewhat novel because the spatial domains for the variables describing the population densities of the interacting species overlap but do not coincide. Conditions implying coexistence of the two species or the extinction of one species are derived. The conditions involve the principal eigenvalues of elliptic operators arising from linearizations of the model system around equilibria with only one species present. The conditions for persistence or extinction are made explicit in terms of the parameters of the system and the geometry of the underlying spatial domains via estimates of the principal eigenvalues. The implications of the models with respect to conservation and refuge design are discussed. Received: 10 June 1999 / Revised version: 7 July 2000 / Published online: 20 December 2000     

Global analysis of the Michaelis-Menten-type ratio-dependent predator-prey system

The recent broad interest on ratio-dependent based predator functional response calls for detailed qualitative study on ratio-dependent predator-prey differential systems. A first such attempt is documented in the recent work of Kuang and Beretta(1998), where Michaelis-Menten-type ratio-dependent model is studied systematically. Their paper, while contains many new and significant results, is far from complete in answering the many subtle mathematical questions on the global qualitative behavior of solutions of the model. Indeed, many of such important open questions are mentioned in the discussion section of their paper. Through a simple change of variable, we transform the Michaelis-Menten-type ratio-dependent model to a better studied Gause-type predator-prey system. As a result, we can obtain a complete classification of the asymptotic behavior of the solutions of the Michaelis-Menten-type ratio-dependent model. In some cases we can determine how the outcomes depend on the initial conditions. In particular, open questions on the global stability of all equilibria in various cases and the uniqueness of limit cycles are resolved. Biological implications of our results are also presented.     

Mathematical modelling of the intravenous glucose tolerance test

 Several attempts at building a satisfactory model of the glucose-insulin system are recorded in the literature. The minimal model, which is the model currently mostly used in physiological research on the metabolism of glucose, was proposed in the early eighties for the interpretation of the glucose and insulin plasma concentrations following the intravenous glucose tolerance test. It is composed of two parts: the first consists of two differential equations and describes the glucose plasma concentration time-course treating insulin plasma concentration as a known forcing function; the second consists of a single equation and describes the time course of plasma insulin concentration treating glucose plasma concentration as a known forcing function. The two parts are to be separately estimated on the available data. In order to study glucose-insulin homeostasis as a single dynamical system, a unified model would be desirable. To this end, the simple coupling of the original two parts of the minimal model is not appropriate, since it can be shown that, for commonly observed combinations of parameter values, the coupled model would not admit an equilibrium and the concentration of active insulin in the “distant” compartment would be predicted to increase without bounds. For comparison, a simple delay-differential model is introduced, is demonstrated to be globally asymptotically stable around a unique equilibrium point corresponding to the pre-bolus conditions, and is shown to have positive and bounded solutions for all times. The results of fitting the delay-differential model to experimental data from ten healthy volunteers are also shown. It is concluded that a global unified model is both theoretically desirable and practically usable, and that any such model ought to undergo formal analysis to establish its appropriateness and to exclude conflicts with accepted physiological notions. Received: 22 June 1998 / Revised version: 24 February 1999     

A generalized transport model for biased cell migration in an anisotropic environment

 A generalized transport model is derived for cell migration in an anisotropic environment and is applied to the specific cases of biased cell migration in a gradient of a stimulus (taxis; e.g., chemotaxis or haptotaxis) or along an axis of anisotropy (e.g., contact guidance). The model accounts for spatial or directional dependence of cell speed and cell turning behavior to predict a constitutive cell flux equation with drift velocity and diffusivity tensor (termed random motility tensor) that are explicit functions of the parameters of the underlying random walk model. This model provides the connection between cell locomotion and the resulting persistent random walk behavior to the observed cell migration on longer time scales, thus it provides a framework for interpreting cell migration data in terms of underlying motility mechanisms. Received: 8 April 1999     

Parametric analysis of the ratio-dependent predator–prey model

We present a complete parametric analysis of stability properties and dynamic regimes of an ODE model in which the functional response is a function of the ratio of prey and predator abundances. We show the existence of eight qualitatively different types of system behaviors realized for various parameter values. In particular, there exist areas of coexistence (which may be steady or oscillating), areas in which both populations become extinct, and areas of "conditional coexistence" depending on the initial values. One of the main mathematical features of ratio-dependent models, distinguishing this class from other predator-prey models, is that the Origin is a complicated equilibrium point, whose characteristics crucially determine the main properties of the model. This is the first demonstration of this phenomenon in an ecological model. The model is investigated with methods of the qualitative theory of ODEs and the theory of bifurcations. The biological relevance of the mathematical results is discussed both regarding conservation issues (for which coexistence is desired) and biological control (for which extinction is desired).     

Stability analysis of the partial selfing selection model

We undertake a detailed study of the one-locus two-allele partial selfing selection model. We show that a polymorphic equilibrium can exist only in the cases of overdominance and underdominance and only for a certain range of selfing rates. Furthermore, when it exists, we show that the polymorphic equilibrium is unique. The local stability of the polymorphic equilibrium is investigated and exact analytical conditions are presented. We also carry out an analysis of local stability of the fixation states and then conclude that only overdominance can maintain polymorphism in the population. When the linear local analysis is inconclusive, a quadratic analysis is performed. For some sets of selective values, we demonstrate global convergence. Finally, we compare and discuss results under the partial selfing model and the random mating model.     

The interaction of thin-film flow, bacterial swarming and cell differentiation in colonies of Serratia liquefaciens

 The rate of expansion of bacterial colonies of S. liquefaciens is investigated in terms of a mathematical model that combines biological as well as hydrodynamic processes. The relative importance of cell differentiation and production of an extracellular wetting agent to bacterial swarming is explored using a continuum representation. The model incorporates aspects of thin film flow with variable suspension viscosity, wetting, and cell differentiation. Experimental evidence suggests that the bacterial colony is highly sensitive to its environment and that a variety of mechanisms are exploited in order to proliferate on a variety of surfaces. It is found that a combination of effects are required to reproduce the variation of bacterial colony motility over a large range of nutrient availability and medium hardness. Received: 29 April 1999     

A dynamic model for the p53 stress response networks under ion radiation

Summary. P53 controls the cell cycle arrest and cell apoptosis through interaction with the downstream genes and their signal pathways. To stimulate the investigation into the complicated responses of p53 under the circumstance of ion radiation (IR) in the cellular level, a dynamic model for the p53 stress response networks is proposed. The model can be successfully used to simulate the dynamic processes of generating the double-strand breaks (DSBs) and their repairing, ataxia telangiectasia mutated (ATM) activation, as well as the oscillations occurring in the p53-MDM2 feedback loop.     

An alternative stochastic model of generation of oligodendrocytes in cell culture

According to our previous model, oligodendrocyte – type 2 (O-2A) astrocyte progenitor cells become competent for differentiation in vitro after they complete a certain number of critical mitotic cycles. After attaining the competency to differentiate, progenitor cells divide with fixed probability p in subsequent cycles. The number of critical cycles is random; analysis of data suggests that it varies from zero to two. The present paper presents an alternative model in which there are no critical cycles, and the probability that a progenitor cell will divide again decreases gradually to a plateau value as the number of completed mitotic cycles increases. In particular all progenitor cells have the ability to differentiate from the time of plating. The Kiefer-Wolfowitz procedure is used to fit the new model to experimental data on the clonal growth of purified O-2A progenitor cells obtained from the optic nerves of 7 day old rats. The new model is shown to fit the experimental data well, indicating that it is not possible to determine whether critical cycles exist on the basis of these experimental data. In contrast to the fit of the previous model, which suggested that the addition of thyroid hormone increased the limiting probability of differentiation as the number of mitotic cycles increases, the fit of the new model suggests that the addition of thyroid hormone has almost no effect on the limiting probability of differentiation. Received: 6 March 2000 / Revised version: 18 September 2000 / Published online: 30 April 2001     

Mathematical modeling of the onset of capillary formation initiating angiogenesis

It is well accepted that neo-vascular formation can be divided into three main stages (which may be overlapping): (1) changes within the existing vessel, (2) formation of a new channel, (3) maturation of the new vessel. In this paper we present a new approach to angiogenesis, based on the theory of reinforced random walks, coupled with a Michaelis-Menten type mechanism which views the endothelial cell receptors as the catalyst for transforming angiogenic factor into proteolytic enzyme in order to model the first stage. In this model, a single layer of endothelial cells is separated by a vascular wall from an extracellular tissue matrix. A coupled system of ordinary and partial differential equations is derived which, in the presence of an angiogenic agent, predicts the aggregation of the endothelial cells and the collapse of the vascular lamina, opening a passage into the extracellular matrix. We refer to this as the onset of vascular sprouting. Some biological evidence for the correctness of our model is indicated by the formation of teats in utero. Further evidence for the correctness of the model is given by its prediction that endothelial cells will line the nascent capillary at the onset of capillary angiogenesis. Received: 27 May 1999 / Revised version: 28 December 1999 / Published online: 16 February 2001     

Programmed cell death: similarities and differences in animals and plants. A flower paradigm

Summary. After an overview of the criteria for the definition of cell death in the animal cell and of its different types of death, a comparative analysis of PCD in the plant cell is reported. The cytological characteristics of the plant cell undergoing PCD are described. The role of plant hormones and growth factors in the regulation of this event is discussed with particular emphasis on PCD activation or prevention by polyamine treatment (doses, timing and developmental stage of the organism) in a Developmental cell death plant model: the Nicotiana tabacum (tobacco) flower corolla. Some of the effects of polyamines might be mediated by transglutaminase catalysis. The activity of this enzyme was examined in different parts of the corolla during its life span showing an acropetal trend parallel to the cell death wave. The location of transglutaminase in some sub-cellular compartments suggests that it exerts different functions in the corolla DCD.     

Influence of nitric oxide synthase activity on amino acid concentration in the quinolinate lesioned rat striatum: A microdialysis and histochemical study

Summary. The influence of nitric oxide synthase (NOS) activity on the KCl-evoked amino acid concentrations was investigated by in vivo microdialysis in the striatum in a rat model of excitotoxic lesion. Basal microdialysate levels of amino acids decreased during the quinolinic acid-induced neurodegeneration process, except for glutamine that increased initially and returned to control values 30 days after quinolinic acid exposure. KCl-evoked increase of extracellular amino acid concentration was reduced due to NOS activity in the striatum of both controls and lesioned animals, except for 120 days after quinolinic acid injection. These changes of amino acid concentrations in microdialysates correlated with the known biochemistry of the consecutive domineered cell types during the lesion process as revealed by histochemistry for NOS, NADPH-diaphorase, GFAP and isolectin B4. The present data provide direct evidence that NOS activity can modulate extracellular amino acid concentrations in the striatum not only under physiological conditions, but also during a pharmacologically induced lesion process and, thus, suggests that nitric oxide affects neurodegeneration via this pathway. Received October 20, 1999; Accepted February 25, 2000     

A mathematical model of the stress induced during avascular tumour growth

In this paper a mathematical model is developed to describe the effect of nonuniform growth on the mechanical stress experienced by cells within an avascular tumour. The constitutive law combines the stress-strain relation of linear elasticity with a growth term that is derived by analogy with thermal expansion. To accommodate the continuous nature of the growth process, the law relates the rate of change of the stress tensor to the rate of change of the strain (rather than relating the stress to the strain directly). By studying three model problems which differ in detail, certain characteristic features are identified. First, cells near the tumour boundary, where nutrient levels and cell proliferation rates are high, are under compression. By contrast, cells towards the centre of the tumour, where nutrient levels are low and cell death dominant, are under tension. The implications of these results and possible model developments are also discussed. Received: 15 November 1999 / Published online: 5 May 2000     

Stochastic models of a parasitic infection, exhibiting three basic reproduction ratios

Two closely related stochastic models of parasitic infection are investigated: a non-linear model, where density dependent constraints are included, and a linear model appropriate to the initial behaviour of an epidemic. Host-mortality is included in both models. These models are appropriate to transmission between homogeneously mixing hosts, where the amount of infection which is transferred from one host to another at a single contact depends on the number of parasites in the infecting host. In both models, the basic reproduction ratio R0 can be defined to be the lifetime expected number of offspring of an adult parasite under ideal conditions, but it does not necessarily contain the information needed to separate growth from extinction of infection. In fact we find three regions for a certain parameter where different combinations of parameters determine the behavior of the models. The proofs involve martingale and coupling methods.     

Selenomethionine induces polyamine biosynthesis in regenerating rat liver tissue

Summary. Our study was undertaken to elucidate the effects of selenomethionine (SeMet) on polyamine metabolism in regenerating rat liver tissue, as useful model of rapidly growing normal tissue. We have examined the levels of spermine, spermidine and putrescine in liver tissue. At the same time we have evaluated the activities of polyamine oxidase (PAO) and diamine oxidase (DAO), the catabolic enzymes of polyamine metabolism. The obtained results suggest that polyamine levels in regenerating liver tissue, at 7^th day after two-thirds partial hepatectomy, were higher in comparison with control group. The administration of selenomethionine to hepatectomized animals during seven days, in a single daily dose of 2.5 μg/100 g body weight, increases the amount of spermine and spermidine; the level of putrescine does not change under the influence of SeMet in regenerating rat liver tissue. PAO activity is lower in regenerating hepatic tissue than in control group. Supplementation of hepatectomized animals with SeMet significantly decreases the activity of this enzyme. DAO activity was significantly higher in hepatectomized and in operated animals treated with SeMet compared to the sham-operated and control ones. The differential sensitivity observed in our model of highly proliferating normal tissue to SeMet, compared with the reported anticancer activity of this molecule is discussed.     

Stochastic epidemics in dynamic populations: quasi-stationarity and extinction

Empirical evidence shows that childhood diseases persist in large communities whereas in smaller communities the epidemic goes extinct (and is later reintroduced by immigration). The present paper treats a stochastic model describing the spread of an infectious disease giving life-long immunity, in a community where individuals die and new individuals are born. The time to extinction of the disease starting in quasi-stationarity (conditional on non-extinction) is exponentially distributed. As the population size grows the epidemic process converges to a diffusion process. Properties of the limiting diffusion are used to obtain an approximate expression for τ, the mean-parameter in the exponential distribution of the time to extinction for the finite population. The expression is used to study how τ depends on the community size but also on certain properties of the disease/community: the basic reproduction number and the means and variances of the latency period, infectious period and life-length. Effects of introducing a vaccination program are also discussed as is the notion of the critical community size, defined as the size which distinguishes between the two qualitatively different behaviours. Received: 14 February 2000 / Revised version: 5 June 2000 / Published online: 24 November 2000