Journal of Mathematical Biology

In this paper, we use an integrodifference equation model and pairwise invasion analysis to find what dispersal strategies are evolutionarily stable strategies (also known as evolutionarily steady or ESS) when there is spatial heterogeneity and possibly seasonal variation in habitat suitability. In that case there are both advantages and disadvantages of dispersing. We begin with the case where all spatial locations can support a viable population, and then consider the case where there are non-viable regions in the habitat. If the viable regions vary seasonally, and the viable regions in summer and winter do not overlap, dispersal may really be necessary for sustaining a population. Our findings generally align with previous findings in the literature that were based on other modeling frameworks, namely that dispersal strategies associated with ideal free distributions are evolutionarily stable. In the case where only part of the habitat can sustain a population, we show that a partial occupation ideal free distribution that occupies only the viable region is associated with a dispersal strategy that is evolutionarily stable. As in some previous works, the proofs of these results make use of properties of line sum symmetric functions, which are analogous to those of line sum symmetric matrices but applied to integral operators.

In this paper we present a model for erythropoiesis under the basic assumption that sufficient iron availability is guaranteed. An extension of the model including a sub-model for the iron dynamics in the body is topic of present research efforts. The model gives excellent results for a number of important situations: recovery of the red blood cell mass after blood donation, adaptation of the number of red blood cells to changes in the altitude of residence and, most important, the reaction of the body to different administration regimens of erythropoiesis stimulating agents, as for instance in the case of pre-surgical administration of Epoetin-α. The simulation results concerning the last item show that choosing an appropriate administration regimen can reduce the total amount of the administered drug considerably. The core of the model consists of structured population equations for the different cell populations which are considered. A key feature of the model is the incorporation of neocytolysis.

Cholera is an environmentally driven disease where the human hosts both ingest the pathogen from polluted environment and shed the pathogen to the environment, generating a two-way feedback cycle. In this paper, we propose a bidirectionally linked immuno-epidemiological model to study the interaction of within- and between-host cholera dynamics. We conduct a rigorous analysis for this multiscale model, with a focus on the stability and bifurcation properties of each feasible equilibrium. We find that the parameter that represents the bidirectional connection is a key factor in shaping the rich dynamics of the system, including the occurrence of the backward bifurcation and Hopf bifurcation. Numerical results illustrate a practical application of our model and add new insight into the prevention and intervention of cholera epidemics.

A selection model which comprises models of additive fecundities as well as models of viability, fecundity, or differential mating selection acting only in one sex, is investigated for an autosomal gene locus in a population reproducing in nonoverlapping generations. The recurrence equations and basic properties of the genotypic population trajectories and equilibrium points are formulated for the multiallelic case. For the diallelic case, the trajectory development is discussed in more detail, and it is proven that every population trajectory converges to a Hardy-Weinberg equilibrium point.

A model for the sterile insect release method of pest control in which the target species is under predatory or parasitic regulation is analyzed. The equations are nondimensionalized and the rescaled parameters are interpreted. There are four types of equilibria, whose existence and stability depend on which of ten regions of parameter space contain the rescaled parameters, and in turn give minimal release rates to achieve eradication of the pest. In at least one region, Hopf bifurcation theory shows the existence of limit cycles, but they are found to be unstable. In addition, the optimal release rate to minimize a total cost functional for pest control by the sterile release method is studied. Both approaches show that when predation accounts for a large fraction of the natural deaths, the necessary release rate and total cost are higher than for weak predation. If the predators are removed without being replaced by any other source of mortality, the cost rises in all cases but rises much more dramatically for cases with strong predation. A definite danger of the sterile release method when some predatory control exists is that the predators are frequently driven extinct before the prey, so that the target species could explode to much higher levels and be more difficult to eradicate again after the sterile release is terminated.

We present a model for the coevolution of seed size and germination time within a season when both affect the ability of the seedlings to compete for space. We show that even in the absence of a morphological or physiological constraint between the two traits, a correlation between seed size and germination time is nevertheless likely to evolve. This raises the more general question to what extent a correlation between any two traits should be considered as an a priori constraint or as an evolved means (or “instrument”) to actually implement a beneficial combination of traits. We derive sufficient conditions for the existence of a positive or a negative correlation. We develop a toy model for seed and seedling survival and seedling growth and use this to illustrate in practice how to determine correlations between seed size and germination time.

 In a class of model vascular trees having distensible blood vessels, we prove that flow partitioning throughout the tree remains constant, independent of the nonzero driving flow (or nonzero inlet to terminal outlet pressure difference). Underlying assumptions are: (1) every vessel in the tree exhibits the same distensibility relationship given by $D/D_0 = f(P)$ where $D$ is the diameter which results from distending pressure $P$ and $D_0$ is the diameter of the individual vessel at zero pressure (each vessel may have its own individual $D_0$). The choice of $f(P)$ includes distensibilities often used in vessel biomechanics modeling, e.g., $f(P) = 1 + \alpha P$ or $f(P) = b + (1-b) \exp(-c P)$, as well as $f(P)$ which exhibit autoregulatory behavior. (2) Every terminal vessel in the tree is subjected to the same terminal outlet pressure. (3) Bernoulli effects are ignored. (4) Flow is nonpulsatile. (5) Blood viscosity within any individual vessel is constant. The results imply that for a vascular tree consistent with assumptions 2–5, the flow distribution calculations based on a rigid geometry, e.g., $D=D_0$, also gives the flow distribution when assuming the common distensibility relationships.

This paper investigates the partial differential equation for the evolving distribution of prostate-specific antigen (PSA) levels following radiotherapy. We also present results on the behavior of moments for the evolving distribution of PSA levels and estimate the probability of long-term treatment success and failure related to values of treatment and disease parameters. Results apply to a much wider range of parameter values than was considered in earlier studies, including parameter combinations that are patient specific.

Consider a positively regular, slightly supercritical branching process withK types. An approximation to the probability of survival of a line descended from a single individual of typei has recently been derived by Hoppe. IfK is large, however, this approximation may not be easy to compute. A further approximation that is easily computable is given. The result is used to estimate probabilities of survival of an allele A that is originally present in one male or one female in a large, random mating, age-structured population. Both autosomal and sex-linked loci are considered. Another application of the approximation is also discussed.