Subharmonic resonance and chaos in forced excitable systems

 Forced excitable systems arise in a number of biological and physiological applications and have been studied analytically and computationally by numerous authors. Existence and stability of harmonic and subharmonic solutions of a forced piecewise-linear Fitzhugh-Nagumo-like system were studied in Othmer ad Watanabe (1994) and in Xie et al. (1996). The results of those papers were for small and moderate amplitude forcing. In this paper we study the existence of subharmonic solutions of this system under large-amplitude forcing. As in the case of intermediate-amplitude forcing, bistability between 1 : 1 and 2 : 1 solutions is possible for some parameters. In the case of large-amplitude forcing, bistability between 2 : 2 and 2 : 1 solutions, which does not occur in the case of intermediate-amplitude forcing, is also possible for some parameters. We identify several new canonical return maps for a singular system, and we show that chaotic dynamics can occur in some regions of parameter space. We also prove that there is a direct transition from 2 : 2 phase-locking to chaos after the first period-doubling bifurcation, rather than via the infinite sequence of period doublings seen in a smooth quadratic interval map. Coexistence of chaotic dynamics and stable phase-locking can also occur. Received: 6 July 1998 / Revised version: 2 October 1998     

Overcompensatory recruitment and generation delay in discrete age-structured population models

 The effect of overcompensatory recruitment and the combined effect of overcompensatory recruitment and generation delay in discrete nonlinear age-structured population models is studied. Considering overcompensatory recruitment alone, we present formal proofs of the supercritical nature of bifurcations (both flip and Hopf) as well as an extensive analysis of dynamics in unstable parameter regions. One important finding here is that in case of small and moderate year to year survival probabilities there are large regions in parameter space where the qualitative behaviour found in a general n+1 dimensional model is retained already in a one-dimensional model. Another result is that the dynamics at or near the boundary of parameter space may be very complicated. Generally, generation delay is found to act as a destabilizing effect but its effect on dynamics is by no means unique. The most profound effect occurs in the n-generation delay cases. In these cases there is no stable equilibrium X ^* at all, but whenever X ^* small, a stable cycle of period n+1 where the periodic points in the cycle are on a very special form. In other cases generation delay does not alter the dynamics in any substantial way. Received 25 April 1995; received in revised form 21 November 1995     

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     

Communication networks: social environments for receiving and signalling behaviour

Communication and social behaviour are inextricably linked, with communication mediating important social behaviours such as resource defence and mate attraction. However, the social environment in which communication occurs is often ignored in discussions of communication behaviour. We argue that networks of several individuals are the common social environment for communication behaviour. The consequences for receivers and signallers of communicating in a network environment are the main subjects of this review. Eavesdropping is a receiving behaviour that is only possible in the environment of a network and therefore we concentrate on this behaviour. The main effect of communication networks on signallers is to create competition with other signallers for receiver attention. We discuss the consequences of such competition. To conclude, we explore the role of signals and signalling interactions as sources of information that animals exploit to direct their behaviour. Received: 1 December 1999 / Received in revised form: 18 January 2000 / Accepted: 18 January 2000     

Optimization models for the force and energy in competitive running

 In [2] the author has developed an optimization model for the force and energy in competitive running. In this paper the energy processes in the muscle were described by a three-compartment hydraulic model. Here this is reviewed briefly and applied to the current world records in order to determine the key parameters, maximal force, energy reserves and oxygen uptake. These values agree well with those given in the literature and those obtained by other means. The velocity profiles for 100 m sprints are described equally well. The model is then applied to older world records to deduce a relation between the force and energy by linear regression. Finally the fully parameterized model is used to compute the effects of adverse wind and altitude. Inasmuch as there are data available, there is a good agreement. Received 19 July 1995; received in revised form 27 February 1996     

Convergence in almost periodic Fisher and Kolmogorov models

 We study convergence of positive solutions for almost periodic reaction diffusion equations of Fisher or Kolmogorov type. It is proved that under suitable conditions every positive solution is asymptotically almost periodic. Moreover, all positive almost periodic solutions are harmonic and uniformly stable, and if one of them is spatially homogeneous, then so are others. The existence of an almost periodic global attractor is also discussed. Received: 11 November 1996 / Revised version: 8 January 1998     

Quasi-steady-state approximation for chemical reaction networks

 The parameter embedding leading to the quasi-steady-state approximation of Heinrich [9] is investigated within the theory of invariant manifolds of Fenichel [4] in order to clarify the essential assumptions needed for this reduction to a low dimensional system. In particular, the concept of pool-variables can be avoided in this generalized approach. Moreover, the dominating influence of the slow subnetwork over the complementary fast subnetwork is interpreted geometrically and in chemical terms and this can be seen as an “enslaving” of the fast subsystem by the slow subsystem. Finally, the results are applied to a system of slime mould communication [6, 7, 13] and to a maltose transport system [2, 3]. Received: 16 June 1997     

Analysis of an excitable dendritic spine with an activity-dependent stem conductance

 Dendritic spines are the major target for excitatory synaptic inputs in the vertebrate brain. They are tiny evaginations of the dendritic surface consisting of a bulbous head and a tenuous stem. Spines are considered to be an important locus for plastic changes underlying memory and learning processes. The findings that synaptic morphology may be activity-dependent and that spine head membrane may be endowed with voltage-dependent (excitable) channels is the motivation for this study. We first explore the dynamics, when an excitable, yet morphologically fixed spine receives a constant current input. Two parameter Andronov–Hopf bifurcation diagrams are constructed showing stability boundaries between oscillations and steady-states. We show how these boundaries can change as a function of both the spine stem conductance and the conductance load of the attached dendrite. Building on this reference case an idealized model for an activity-dependent spine is formulated and analyzed. Specifically we examine the possibility that the spine stem resistance, the tunable “synaptic weight” parameter identified by Rall and Rinzel, is activity-dependent. In the model the spine stem conductance depends (slowly) on the local electrical interactions between the spine head and the dendritic cable; parameter regimes are found for bursting, steady states, continuous spiking, and more complex oscillatory behavior. We find that conductance load of the dendrite strongly influences the burst pattern as well as other dynamics. When the spine head membrane potential exhibits relaxation oscillations a simple model is formulated that captures the dynamical features of the full model. Received: 10 January 1997/Revised version: 25 March 1997     

Global stability in consumer-resource cascades

 Models of population growth in consumer-resource cascades (serially arranged containers with a dynamic consumer population, v, receiving a flow of resource, u, from the previous container) with a functional response of the form h(u/v ^b ) are investigated. For b∈[0, 1], it is shown that these models have a globally stable equilibrium. As a result, two conclusions can be drawn: (1) Consumer density dependence in the functional or in the per-capita numerical response can result in persistence of the consumer population in all containers. (2) In the absence of consumer density dependence, the consumer goes extinct in all containers except possibly the first. Several variations of this model are discussed including replacing discrete containers by a spatial continuum and introducing a dynamic resource. Received 25 February 1995 / received in revised form 27 July 1995     

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     

Discrete chiasma formation models and their associated high order interference

. We introduce some special chiasma formation processes. First a family of discrete chiasma formation processes is introduced and we determine the nature of higher order interference associated with those processes. Secondly we consider a two-stage chiasma formation process, where the associated recombination frequency between two markers depends not only on their map distance but also on their location along the chromosomes. We characterise under this process, in some cases, the nature of interference between two segments. Received: 22 January 1996 / Revised version: 17 September 1997     

The effects of the functional response on the bifurcation behavior of a mite predator–prey interaction model

 A predator–prey interaction model based on a system of differential equations with temperature-dependent parameters chosen appropriately for a mite interaction on apple trees is analyzed to determine how the type of functional response influences bifurcation and stability behavior. Instances of type I, II, III, and IV functional responses are considered, the last of which incorporates prey interference with predation. It is shown that the model systems with the type I, II, and III functional responses exhibit qualitatively similar bifurcation and stability behavior over the interval of definition of the temperature parameter. Similar behavior is found in the system with the type IV functional response at low levels of prey interference. Higher levels of interference are destabilizing, as illustrated by the prevalence of bistability and by the presence of three attractors for some values of the model parameters. All four systems are capable of modeling population oscillations and outbreaks. Received 12 March 1996; received in revised form 25 October 1996     

Sampling distributions, biases, variances, and confidence intervals for genetic correlations

Genetic correlations are frequently estimated from natural and experimental populations, yet many of the statistical properties of estimators of are not known, and accurate methods have not been described for estimating the precision of estimates of Our objective was to assess the statistical properties of multivariate analysis of variance (MANOVA), restricted maximum likelihood (REML), and maximum likelihood (ML) estimators of by simulating bivariate normal samples for the one-way balanced linear model. We estimated probabilities of non-positive definite MANOVA estimates of genetic variance-covariance matrices and biases and variances of MANOVA, REML, and ML estimators of and assessed the accuracy of parametric, jackknife, and bootstrap variance and confidence interval estimators for MANOVA estimates of were normally distributed. REML and ML estimates were normally distributed for but skewed for and 0.9. All of the estimators were biased. The MANOVA estimator was less biased than REML and ML estimators when heritability (H), the number of genotypes (n), and the number of replications (r) were low. The biases were otherwise nearly equal for different estimators and could not be reduced by jackknifing or bootstrapping. The variance of the MANOVA estimator was greater than the variance of the REML or ML estimator for most H, n, and r. Bootstrapping produced estimates of the variance of close to the known variance, especially for REML and ML. The observed coverages of the REML and ML bootstrap interval estimators were consistently close to stated coverages, whereas the observed coverage of the MANOVA bootstrap interval estimator was unsatisfactory for some H, n, and r. The other interval estimators produced unsatisfactory coverages. REML and ML bootstrap interval estimates were narrower than MANOVA bootstrap interval estimates for most H, and r. Received: 6 July 1995 / Accepted: 8 March 1996     

Population dynamics and competition in chemostat models with adaptive nutrient uptake

 The standard Monod model for microbial population dynamics in the chemostat is modified to take into consideration that cells can adapt to the change of nutrient concentration in the chemostat by switching between fast and slow nutrient uptake and growing modes with asymmetric thresholds for transition from one mode to another. This is a generalization of a modified Monod model which considers adaptation by transition between active growing and quiescent cells. Global analysis of the model equations is obtained using the theory of asymptotically autonomous systems. Transient oscillatory population density and hysteresis growth pattern observed experimentally, which do not occur for the standard Monod model, can be explained by such adaptive mechanism of the cells. Competition between two species that can switch between fast and slow nutrient uptake and growing modes is also considered. It is shown that generically there is no coexistence steady state, and only one steady state, corresponding to the survival of at most one species in the chemostat, is a local attractor. Numerical simulations reproduce the qualitative feature of some experimental data which show that the population density of the winning species approaches a positive steady state via transient oscillations while that of the losing species approaches the zero steady state monotonically. Received 4 August 1995; received in revised form 15 December 1995     

Taurine and neural cell damage

Summary. The inhibitory amino acid taurine is an osmoregulator and neuromodulator, also exerting neuroprotective actions in neural tissue. We review now the involvement of taurine in neuron-damaging conditions, including hypoxia, hypoglycemia, ischemia, oxidative stress, and the presence of free radicals, metabolic poisons and an excess of ammonia. The brain concentration of taurine is increased in several models of ischemic injury in vivo. Cell-damaging conditions which perturb the oxidative metabolism needed for active transport across cell membranes generally reduce taurine uptake in vitro, immature brain tissue being more tolerant to the lack of oxygen. In ischemia nonsaturable diffusion increases considerably. Both basal and K⁺-stimulated release of taurine in the hippocampus in vitro is markedly enhanced under cell-damaging conditions, ischemia, free radicals and metabolic poisons being the most potent. Hypoxia, hypoglycemia, ischemia, free radicals and oxidative stress also increase the initial basal release of taurine in cerebellar granule neurons, while the release is only moderately enhanced in hypoxia and ischemia in cerebral cortical astrocytes. The taurine release induced by ischemia is for the most part Ca²⁺-independent, a Ca²⁺-dependent mechanism being discernible only in hippocampal slices from developing mice. Moreover, a considerable portion of hippocampal taurine release in ischemia is mediated by the reversal of Na⁺-dependent transporters. The enhanced release in adults may comprise a swelling-induced component through Cl⁻ channels, which is not discernible in developing mice. Excitotoxic concentrations of glutamate also potentiate taurine release in mouse hippocampal slices. The ability of ionotropic glutamate receptor agonists to evoke taurine release varies under different cell-damaging conditions, the N-methyl-D-aspartate-evoked release being clearly receptor-mediated in ischemia. Neurotoxic ammonia has been shown to provoke taurine release from different brain preparations, indicating that the ammonia-induced release may modify neuronal excitability in hyperammonic conditions. Taurine released simultaneously with an excess of excitatory amino acids in the hippocampus under ischemic and other neuron-damaging conditions may constitute an important protective mechanism against excitotoxicity, counteracting the harmful effects which lead to neuronal death. The release of taurine may prevent excitation from reaching neurotoxic levels. Received January 25, 2000/Accepted January 31, 2000     

Travelling wave phenomena in non-linear diffusion degenerate Nagumo equations

 In this paper we study the existence of one-dimensional travelling wave solutions u(x, t)=φ(x−ct) for the non-linear degenerate (at u=0) reaction-diffusion equation u _t =[D(u)u _x ] _x +g(u) where g is a generalisation of the Nagumo equation arising in nerve conduction theory, as well as describing the Allee effect. We use a dynamical systems approach to prove: 1. the global bifurcation of a heteroclinic cycle (two monotone stationary front solutions), for c=0, 2. The existence of a unique value c ^*>0 of c for which φ(x−c ^* t) is a travelling wave solution of sharp type and 3. A continuum of monotone and oscillatory fronts for c≠c ^*. We present some numerical simulations of the phase portrait in travelling wave coordinates and on the full partial differential equation. Received 15 December 1995; received in revised form 14 May 1996     

Bifurcations in haploid and diploid sequence space models

 Deterministic models of mutation and selection in the space of (binary) nucleotide-type sequences have been investigated for haploid populations during the past 25 years, and, recently, for diploid populations as well. These models, in particular their ‘error thresholds’, have mainly been analyzed by numerical methods and perturbation techniques. We consider them here by means of bifurcation theory, which improves our understanding of both equilibrium and dynamical properties. In a caricature obtained from the original model by neglecting back mutation to the favourable allele, the familiar error threshold of the haploid two-class model turns out to be a simple transcritical bifurcation, whereas its diploid counterpart exhibits an additional saddle node. This corresponds to a second error threshold. Three-class models with neutral spaces of unequal size introduce further features. Such are a global bifurcation in haploid populations, and simple examples of Hopf bifurcations (as predicted by Akin’s theorem) in the diploid case. Received 13 June 1995; received in revised form 26 July 1996     

Genealogy and subpopulation differentiation under various models of population structure

 The structured coalescent is used to calculate some quantities relating to the genealogy of a pair of homologous genes and to the degree of subpopulation differentiation, under a range of models of subdivided populations and assuming the infinite alleles model of neutral mutation. The classical island and stepping-stone models of population structure are considered, as well as two less symmetric models. For each model, we calculate the Laplace transform of the distribution of the coalescence time of a pair of genes from specified locations and the corresponding mean and variance. These results are then used to calculate the values of Wright’s coefficient F _ST , its limit as the mutation rate tends to zero and the limit of its derivative with respect to the mutation rate as the mutation rate tends to zero. From this derivative it is seen that F _ST can depend strongly on the mutation rate, for example in the case of an essentially one-dimensional habitat with many subpopulations where gene flow is restricted to neighbouring subpopulations. Received: 1 October 1997 / Revised version: 15 March 1998     

Interaction of maturation delay and nonlinear birth in population and epidemic models

 A population with birth rate function B(N) N and linear death rate for the adult stage is assumed to have a maturation delay T>0. Thus the growth equation N′(t)=B(N(t−T)) N(t−T) e⁻ d ₁ T−dN(t) governs the adult population, with the death rate in previous life stages d ₁≧0. Standard assumptions are made on B(N) so that a unique equilibrium N _e exists. When B(N) N is not monotone, the delay T can qualitatively change the dynamics. For some fixed values of the parameters with d ₁>0, as T increases the equilibrium N _e can switch from being stable to unstable (with numerically observed periodic solutions) and then back to stable. When disease that does not cause death is introduced into the population, a threshold parameter R ₀ is identified. When R ₀<1, the disease dies out; when R ₀>1, the disease remains endemic, either tending to an equilibrium value or oscillating about this value. Numerical simulations indicate that oscillations can also be induced by disease related death in a model with maturation delay. Received: 2 November 1998 / Revised version: 26 February 1999     

Gliadin polymorphism in wild and cultivated einkorn wheats

To study the relationships between different species of the Einkorn group, 408 accessions of Triticum monococcum, T. boeoticum, T. boeoticum ssp. thauodar and T. urartu were analyzed electrophoretically for their protein composition at the Gli-1 and Gli-2 loci. In all the species the range of allelic variation at the loci examined is remarkable. The gliadin patterns of T. monococcum and T. boeoticum were very similar to one another but differed substantially from those of T. urartu. Several accessions of T. boeoticum and T. monococcum were shown to share the same alleles at the Gli-1 and Gli-2 loci, confirming the recent nomenclature that considers these wheats as different subspecies of the same species, T. monococcum. The gliadin composition of T. urartu resembled that of the A genome of polyploid wheats more than did T. boeoticum or T. monococcum, supporting the hypothesis that T. urartu, rather than T. boeoticum, is the donor of the A genome in cultivated wheats. Because of their high degree of polymorphism the gliadin markers may help in selecting breeding parents from diploid wheat germ plasm collections and can be used both to search for valuable genes linked to the gliadin-coding loci and to monitor the transfer of alien genes into cultivated polyploid wheats. Received: 8 July 1996 / Accepted: 12 July 1996