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.     

Equilibrium structure and stability in a frequency-dependent,two-population diploid model

We investigate the equilibrium structure for an evolutionary genetic model in discrete time involving two monoecious populations subject to intraspecific and interspecific random pairwise interactions. A characterization for local stability of an equilibrium is found, related to the proximity of this equilibrium with evolutionarily stable strategies (ESS). This extends to a multi-population framework a principle initially proposed for single populations, which states that the mean population strategy at a locally stable equilibrium is as close as possible to an ESS.     

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).     

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.     

A new explicit stability criterion for human periodic breathing<--RID="*"--><--ID="*" This research was supported by grants from DRED (96-920).-->

The aim of this paper is to carry out a stability analysis for periodic breathing in humans that incorporates the dynamic characteristics of ventilation control. A simple CO₂ model that takes into account the main elements of the respiratory system, i.e. the lungs and the ventilatory controller with its dynamic properties, is presented. This model results in a three-dimensional non-linear delay differential system for which there exists a unique equilibrium point. Our stability analysis of this equilibrium point leads to the definition of a new explicit stability criterion and to the demonstration of the existence of a Hopf bifurcation. Numerical simulations illustrate the influence of physiological parameters on the stability of ventilation, and particularly the major role of the dynamic characteristics of the respiratory controller. Received: 2 February 1999 / Revised version: 18 June 1999 / Published online: 23 October 2000     

Models of central pattern generators for quadruped locomotion I. Primary gaits

In this paper we continue the analysis of a network of symmetrically coupled cells modeling central pattern generators for quadruped locomotion proposed by Golubitsky, Stewart, Buono, and Collins. By a cell we mean a system of ordinary differential equations and by a coupled cell system we mean a network of identical cells with coupling terms. We have three main results in this paper. First, we show that the proposed network is the simplest one modeling the common quadruped gaits of walk, trot, and pace. In doing so we prove a general theorem classifying spatio-temporal symmetries of periodic solutions to equivariant systems of differential equations. We also specialize this theorem to coupled cell systems. Second, this paper focuses on primary gaits; that is, gaits that are modeled by output signals from the central pattern generator where each cell emits the same waveform along with exact phase shifts between cells. Our previous work showed that the network is capable of producing six primary gaits. Here, we show that under mild assumptions on the cells and the coupling of the network, primary gaits can be produced from Hopf bifurcation by varying only coupling strengths of the network. Third, we discuss the stability of primary gaits and exhibit these solutions by performing numerical simulations using the dimensionless Morris-Lecar equations for the cell dynamics.     

Traveling waves and pulses in a one-dimensional network of excitable integrate-and-fire neurons

 We study the existence and stability of traveling waves and pulses in a one-dimensional network of integrate-and-fire neurons with synaptic coupling. This provides a simple model of excitable neural tissue. We first derive a self-consistency condition for the existence of traveling waves, which generates a dispersion relation between velocity and wavelength. We use this to investigate how wave-propagation depends on various parameters that characterize neuronal interactions such as synaptic and axonal delays, and the passive membrane properties of dendritic cables. We also establish that excitable networks support the propagation of solitary pulses in the long-wavelength limit. We then derive a general condition for the (local) asymptotic stability of traveling waves in terms of the characteristic equation of the linearized firing time map, which takes the form of an integro-difference equation of infinite order. We use this to analyze the stability of solitary pulses in the long-wavelength limit. Solitary wave solutions are shown to come in pairs with the faster (slower) solution stable (unstable) in the case of zero axonal delays; for non-zero delays and fast synapses the stable wave can itself destabilize via a Hopf bifurcation. Received: 27 October 1998     

A simple SIS epidemic model with a backward bifurcation

It is shown that an SIS epidemic model with a non-constant contact rate may have multiple stable equilibria, a backward bifurcation and hysteresis. The consequences for disease control are discussed. The model is based on a Volterra integral equation and allows for a distributed infective period. The analysis includes both local and global stability of equilibria.     

Regions of Stable Equilibria for Models of Differential Selection in the Two Sexes under Random Mating

The equilibrium structure of models of differential selection in the sexes is investigated. It is shown that opposing additive selection leads to stable polymorphic equilibria for only a restricted set of selection intensities, and that for weak selection the selection intensities must be of approximately the same magnitude in the sexes. General models of opposing directional selection, with arbitrary dominance, are investigated by considering simultaneously the stability properties of the trivial equilibria and the curve along which multiple roots appear. Numerical calculations lead us to infer that the average degree of dominance determines the equilibrium characteristics of models of opposing selection. It appears that if the favored alleles are, on the average, recessive, there may be multiple polymorphic equilibria, whereas only a single polymorphic equilibrium can occur when the favored alleles are, on the average, dominant. The principle that the average degree of dominance controls equilibrium behavior is then extended to models allowing directional selection in one sex with overdominance in the other sex, by showing that polymorphism is maintained if and only if the average fitness in heterozygotes exceeds one.     

Linear stability criteria for population models with periodically perturbed delays

We examine some simple population models that incorporate a time delay which is not a constant but is instead a known periodic function of time. We examine what effect this periodic variation has on the linear stability of the equilibrium states of scalar population models and of a simple predator prey system. The case when the delay differs from a constant by a small amplitude periodic perturbation can be treated analytically by using two-timing methods. Of particular interest is the case when the system is initially marginally stable. The introduction of variation in the delay can then have either a stabilising effect or a destabilizing one, depending on the frequency of the periodic perturbation. The case when the periodic perturbation has large amplitude is studied numerically. If the fluctuation is large enough the effect can be stabilising.     

An inbreeding model of associative overdominance during a population bottleneck

Associative overdominance, the fitness difference between heterozygotes and homozygotes at a neutral locus, is classically described using two categories of models: linkage disequilibrium in small populations or identity disequilibrium in infinite, partially selfing populations. In both cases, only equilibrium situations have been considered. In the present study, associative overdominance is related to the distribution of individual inbreeding levels (i.e., genomic autozygosity). Our model integrates the effects of physical linkage and variation in inbreeding history among individual pedigrees. Hence, linkage and identity disequilibrium, traditionally presented as alternatives, are summarized within a single framework. This allows studying nonequilibrium situations in which both occur simultaneously. The model is applied to the case of an infinite population undergoing a sustained population bottleneck. The effects of bottleneck size, mating system, marker gene diversity, deleterious genomic mutation parameters, and physical linkage are evaluated. Bottlenecks transiently generate much larger associative overdominance than observed in equilibrium finite populations and represent a plausible explanation of empirical results obtained, for instance, in marine species. Moreover, the main origin of associative overdominance is random variation in individual inbreeding whereas physical linkage has little effect.     

Polymorphic equilibria with assortative mating and selection in subdivided populations

Two modes of assortative mating, partial selfing and assorting by phenotypic classes, are studied in a subdivided population. Differential viability is allowed and the selection intensities and assorting tendencies are permitted to vary among the habitats. There exists a symmetric polymorphism; we delimit its level of heterozygosity and stability nature (dependent on the selection intensities and assorting propensities). This complements studies of the fixation states and thereby provides further insight into the global equilibrium structure in subdivided populations. Circumstances are given where the fixation states and symmetric polymorphism comprise the global equilibrium structure. Examples are also given where migration engenders stable polymorphic equilibria and stable polymorphic equilibrium cycles which are absent in single demes without migration.     

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     

Numerical bifurcation analysis of delay differential equations arising from physiological modeling

This paper has a dual purpose. First, we describe numerical methods for continuation and bifurcation analysis of steady state solutions and periodic solutions of systems of delay differential equations with an arbitrary number of fixed, discrete delays. Second, we demonstrate how these methods can be used to obtain insight into complex biological regulatory systems in which interactions occur with time delays: for this, we consider a system of two equations for the plasma glucose and insulin concentrations in a diabetic patient subject to a system of external assistance. The model has two delays: the technological delay of the external system, and the physiological delay of the patient's liver. We compute stability of the steady state solution as a function of two parameters, compare with analytical results and compute several branches of periodic solutions and their stability. These numerical results allow to infer two categories of diabetic patients for which the external system has different efficiency.     

The evolution of a Müllerian mimic in a spatially distributed community

Strong positive density-dependence should lead to a loss of diversity, but warning-colour and Müllerian mimicry systems show extraordinary levels of diversity. Here, we propose an analytical model to explore the dynamics of two forms of a Müllerian mimic in a heterogeneous environment with two alternative model species. Two connected populations of a dimorphic, chemically defended mimic are allowed to evolve and disperse. The proportions of the respective model species vary spatially. We use a nonlinear approximation of Müller's number-dependent equations to model a situation where the mortality for either form of the mimic decreases hyberbolically when its local density increases. A first non-spatial analysis confirms that the positive density-dependence makes coexistence of mimetic forms unstable in a single isolated patch, but shows that mimicry of the rarer model can be stable once established. The two-patch analysis shows that when models have different abundance in different places, local mimetic diversity in the mimic is maintained only if spatial heterogeneity is strong, or, more interestingly, if the mimic is not too strongly distasteful. Therefore, mildly toxic species can become polymorphic in a wider range of ecological settings. Spatial dynamics thus reveal a region of Müllerian polymorphism separating classical Batesian polymorphism and Müllerian monomorphism along the mimic's palatability spectrum. Such polymorphism-palatability relationship in a spatial environment provides a parsimonious hypothesis accounting for the observed Müllerian polymorphism that does not require quasi-Batesian dynamics. While the stability of coexistence depends on all factors, only the migration rate and strength of selection appear to affect the level of diversity at the polymorphic equilibrium. Local adaptation is predicted in most polymorphic cases. These results are in very good accordance with recent empirical findings on the polymorphic butterflies Heliconius numata and H. cydno.     

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     

Random binding of dimers to chains

We develop a probabilistic model for the binding of a small linear polymer to a larger chain. We assume that we can approximate the energy of interaction of the two chains by summing the pairwise interactions between subunits. Because the energy of interaction between a pair of subunits can depend on neighboring subunits, which we assume vary along the chain, we assign the pairwise energies of interactions according to a specified probability distribution. Thus we develop a statistical model for the binding of two molecules. While such models may not be appropriate for studying the interaction of a particular pair of molecules, they can provide insight into questions that deal with populations of molecules, such as why do MHC molecules bind peptides of a certain size? Here we analyze in detail the special case of a heterodimer binding to a polymer.     

Partial male sterility and the evolution of nuclear gynodioecy in plants

Gynodioecy, a genetic dimorphism of females and hermaphrodites, is pertinent to an understanding of the evolution of plant gender, mating and genetic variability. Classical models of nuclear gynodioecy attribute the maintenance of the dimorphism to frequency-dependent selection in which the female phenotype has a fitness advantage at low frequency owing to a doubled ovule fertility. Here, I analyse explicit genetic models of nuclear gynodioecy that expand on previous work by allowing partial male sterility in combination with either fixed or dynamically evolving mutational inbreeding depression. These models demonstrate that partial male sterility causes fitness underdominance at the mating locus, which can prevent the spread of females. However, if partial male sterility is compensated by a change in selfing rate, overdominance at the mating locus can cause the spread of females. Overdominance at introduction of the male sterility allele can be caused by high inbreeding depression and a lower selfing rate in the heterozygote, by purging of mutations by a higher selfing rate in the heterozygote, and by low inbreeding depression and a higher selfing rate in the heterozygote. These processes might be of general importance in the maintenance of mating polymorphisms in plants.     

The diffusion model for migration and selection in a plant population

 The diffusion approximation is derived for migration and selection at a multiallelic locus in a partially selfing plant population subdivided into a lattice of colonies. Generations are discrete and nonoverlapping; both pollen and seeds disperse. In the diffusion limit, the genotypic frequencies at each point are those determined at equilibrium by the local rate of selfing and allelic frequencies. If the drift and diffusion coefficients are taken as the appropriate linear combination of the corresponding coefficients for pollen and seeds, then the migration terms in the partial differential equation for the allelic frequencies have the standard form for a monoecious animal population. The selection term describes selection on the local genotypic frequencies. The boundary conditions and the unidimensional transition conditions for a geographical barrier and for coincident discontinuities in the carrying capacity and migration rate have the standard form. In the diallelic case, reparametrization renders the entire theory of clines and of the wave of advance of favorable alleles directly applicable to plant populations. Received 30 August 1995; received in revised form 23 February 1996     

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