Stability of traveling wave solutions for a nonlocal Lotka-Volterra model.

In this paper, we studied the stability of traveling wave solutions of a two-species Lotka-Volterra competition model in the form of a coupled system of reaction diffusion equations with nonlocal intraspecific and interspecific competitions in space at times. First, the uniform upper bounds for the solutions of the model was proved. By using the anti-weighted method and the energy estimates, the asymptotic stability of traveling waves with large wave speeds of the system was established.

Stereotaxic Viral Injection and Gradient-Index Lens Implantation for Deep Brain In Vivo Calcium Imaging.

A miniature fluorescence microscope (miniscope) is a potent tool for in vivo calcium imaging from freely behaving animals. It offers several advantages over conventional multi-photon calcium imaging systems: (1) compact; (2) light-weighted; (3) affordable; and (4) allows recording from freely behaving animals. This protocol describes brain surgeries for deep brain in vivo calcium imaging using a custom-developed miniscope recording system. The preparation procedure consists of three steps, including (1) stereotaxically injecting the virus at the desired brain region of a mouse brain to label a specific subgroup of neurons with genetically encoded calcium sensor; (2) implantation of gradient-index (GRIN) lens that can relay calcium image from deep brain region to the miniscope system; and (3) affixing the miniscope holder over the mouse skull where miniscope can be attached later. To perform in vivo calcium imaging, the miniscope is fastened onto the holder, and neuronal calcium images are collected along with simultaneous behavior recordings. The present surgery protocol is compatible with any commercial or custom-built single-photon and two-photon imaging systems for deep brain in vivo calcium imaging.

Mathematical modeling of population structure in bioreactors seeded with light-controllable microbial stem cells.

Industrial bioreactors use microbial organisms as living factories to produce a wide range of commercial products. For most applications, yields eventually become limited by the proliferation of "escape mutants" that acquire a growth advantage by losing the ability to make product. The goal of this work is to use mathematical models to determine whether this problem could be addressed in continuous flow bioreactors that include a "stem cell" population that multiplies rapidly and could be used to compete against the emergence of cheater mutants. In this system, external stimuli can be used to induce stem cell multiplication through symmetric cell division, or to limit stem cell multiplication and induce higher production through an asymmetric cell division that produces one stem cell and one new product-producing "factory cell". Our results show product yields from bioreactors with microbial stem cells can be increased by 18% to 127% over conventional methods, and sensitivity analysis shows that yields could be improved over a broad range of parameter space.

Multiple Dose Pharmacokinetic Models Predict Bioavailability of Toxins in Vertebrate Herbivores.

In this paper, compartmental pharmacokinetic models are built to predict the concentration of toxic phytochemical in the gastrointestinal tract and blood following oral intake by an individual vertebrate herbivore. The existing single and multiple dose pharmacokinetic models are extended by inclusion of impulsive differential equations which account for an excretion factor whereby unchanged toxins are excreted in the feces due to gastrointestinal mobility. An index α is defined to measure the fraction of bioavailability attributed to the excretion factor of gastrointestinal motility. Sensitivity analysis was conducted and suggests, for any toxin, the bioavailability index α depends mostly on absorption rate of toxin from gastrointestinal tract into the blood, frequency of elimination due to gastrointestinal motility, and the frequency of toxin intake, under the model assumptions.

Negative frequency-dependent foraging behaviour in a generalist herbivore (Alces alces) and its stabilizing influence on food web dynamics.

Resource selection is widely appreciated to be context-dependent and shaped by both biological and abiotic factors. However, few studies have empirically assessed the extent to which selective foraging behaviour is dynamic and varies in response to environmental conditions for free-ranging animal populations. Here, we assessed the extent that forage selection fluctuated in response to different environmental conditions for a free-ranging herbivore, moose (Alces alces), in Isle Royale National Park, over a 10-year period. More precisely, we assessed how moose selection for coniferous versus deciduous forage in winter varied between geographic regions and in relation to (a) the relative frequency of forage types in the environment (e.g. frequency-dependent foraging behaviour), (b) moose abundance, (c) predation rate (by grey wolves) and (d) snow depth. These factors are potentially important for their influence on the energetics of foraging. We also built a series of food-chain models to assess the influence of dynamic foraging strategies on the stability of food webs. Our analysis indicates that moose exhibited negative frequency dependence, by selectively exploiting rare resources. Frequency-dependent foraging was further mediated by density-dependent processes, which are likely to be predation, moose abundance or some combination of both. In particular, frequency dependence was weaker in years when predation risk was high (i.e. when the ratio of moose to wolves was relatively low). Selection for conifers was also slightly weaker during deep snow years. The food-chain analysis indicates that the type of frequency-dependent foraging strategy exhibited by herbivores had important consequences for the stability of ecological communities. In particular, the dynamic foraging strategy that we observed in the empirical analysis (i.e. negative frequency dependence being mediated by density-dependent processes) was associated with more stable food web dynamics compared to fixed foraging strategies. The results of this study indicated that forage selection is a complex ecological process, varying in response to both biological (predation and moose density) and abiotic factors (snow depth) and over relatively small spatial scales (between regions). This study also provides a useful framework for assessing the influence of other aspects of foraging behaviour on the stability of food web dynamics.

Modelling Wolbachia infection in a sex-structured mosquito population carrying West Nile virus.

Wolbachia is possibly the most studied reproductive parasite of arthropod species. It appears to be a promising candidate for biocontrol of some mosquito borne diseases. We begin by developing a sex-structured model for a Wolbachia infected mosquito population. Our model incorporates the key effects of Wolbachia infection including cytoplasmic incompatibility and male killing. We also allow the possibility of reduced reproductive output, incomplete maternal transmission, and different mortality rates for uninfected/infected male/female individuals. We study the existence and local stability of equilibria, including the biologically relevant and interesting boundary equilibria. For some biologically relevant parameter regimes there may be multiple coexistence steady states including, very importantly, a coexistence steady state in which Wolbachia infected individuals dominate. We also extend the model to incorporate West Nile virus (WNv) dynamics, using an SEI modelling approach. Recent evidence suggests that a particular strain of Wolbachia infection significantly reduces WNv replication in Aedes aegypti. We model this via increased time spent in the WNv-exposed compartment for Wolbachia infected female mosquitoes. A basic reproduction number [Formula: see text] is computed for the WNv infection. Our results suggest that, if the mosquito population consists mainly of Wolbachia infected individuals, WNv eradication is likely if WNv replication in Wolbachia infected individuals is sufficiently reduced.

Delayed action insecticides and their role in mosquito and malaria control.

There is considerable interest in the management of insecticide resistance in mosquitoes. One possible approach to slowing down the evolution of resistance is to use late-life-acting (LLA) insecticides that selectively kill only the old mosquitoes that transmit malaria, thereby reducing selection pressure favoring resistance. In this paper we consider an age-structured compartmental model for malaria with two mosquito strains that differ in resistance to insecticide, using an SEI approach to model malaria in the mosquitoes and thereby incorporating the parasite developmental times for the two strains. The human population is modeled using an SEI approach. We consider both conventional insecticides that target all adult mosquitoes, and LLA insecticides that target only old mosquitoes. According to linearised theory the potency of the insecticide affects mainly the speed of evolution of resistance. Mutations that confer resistance can also affect other parameters such as mean adult life span and parasite developmental time. For both conventional and LLA insecticides the stability of the malaria-free equilibrium, with only the resistant mosquito strain present, depends mainly on these other parameters. This suggests that the main long term role of an insecticide could be to induce genetic changes that have a desirable effect on a vital parameter such as adult life span. However, when this equilibrium is unstable, numerical simulations suggest that a potent LLA insecticide can slow down the spread of malaria in humans but that the timing of its action is very important.

Spatiotemporal mutualistic model of mistletoes and birds.

A mathematical model which incorporates the spatial dispersal and interaction dynamics of mistletoes and birds is derived and studied to gain insights of the spatial heterogeneity in abundance of mistletoes. Fickian diffusion and chemotaxis are used to model the random movement of birds and the aggregation of birds due to the attraction of mistletoes, respectively. The spread of mistletoes by birds is expressed by a dispersal operator, which is typically a convolution integral with a dispersal kernel. Two different types of kernel functions are used to study the model, one is a Dirac delta function which reflects the special case that the spread behavior is local, and the other one is a general non-negative symmetric function which describes the nonlocal spread of mistletoes. When the kernel function is taken as the Dirac delta function, the threshold condition for the existence of mistletoes is given and explored in terms of parameters. For the general non-negative symmetric kernel case, we prove the existence and stability of spatially nonhomogeneous equilibria. Numerical simulations are conducted by taking specific forms of kernel functions. Our study shows that the spatial heterogeneous patterns of mistletoes are related to the specific dispersal pattern of birds which carry mistletoe seeds.

Modeling the dynamics of woody plant-herbivore interactions with age-dependent toxicity.

In this paper we study the effects that woody plant chemical defenses may have on interactions between boreal hares that in winter feed almost entirely on twigs. We focus particularly on the fact that toxin concentration often varies with the age of twig segments. The model incorporates the fact that the woody internodes of the youngest segments of the twigs of the deciduous angiosperm species that these hares prefer to eat are more defended by toxins than the woody internodes of the older segments that subtend and support the younger segments. Thus, the per capita daily intake of the biomass of the older segments of twigs by hares is much higher than their intake of the biomass of the younger segments of twigs. This age-dependent toxicity of twig segments is modeled using age-structured model equations which are reduced to a system of delay differential equations involving multiple delays in the woody plant-hare dynamics. A novel aspect of the modeling was that it had to account for mortality of non-consumed younger twig segment biomass when older twig biomass was bitten off and consumed. Basic mathematical properties of the model are established together with upper and lower bounds on the solutions. Necessary and sufficient conditions are found for the linear stability of the equilibrium in which the hare is extinct, and sufficient conditions are found for the global stability of this equilibrium. Numerical simulations confirmed the analytical results and demonstrated the existence of limit cycles over ranges of parameters reasonable for hares browsing on woody vegetation in boreal ecosystems. This showed that age dependence in plant chemical defenses has the capacity to cause hare-plant population cycles, a new result.

Dynamics of a plant-herbivore-predator system with plant-toxicity.

A system of ordinary differential equations is considered that models the interactions of two plant species populations, an herbivore population, and a predator population. We use a toxin-determined functional response to describe the interactions between plant species and herbivores and use a Holling Type II functional response to model the interactions between herbivores and predators. In order to study how the predators impact the succession of vegetation, we derive invasion conditions under which a plant species can invade into an environment in which another plant species is co-existing with a herbivore population with or without a predator population. These conditions provide threshold quantities for several parameters that may play a key role in the dynamics of the system. Numerical simulations are conducted to reinforce the analytical results. This model can be applied to a boreal ecosystem trophic chain to examine the possible cascading effects of predator-control actions when plant species differ in their levels of toxic defense.

Spatiotemporal variation of mistletoes: a dynamic modeling approach.

Mistletoes are common aerial stem-parasites and their seeds are dispersed by fruit-eating birds. In the mutually beneficial relationships between mistletoes and bird species that disperse mistletoes' seeds, the preference of birds for infected trees influences the spread of mistletoes and the spatiotemporal pattern formation of mistletoes. We formulate a deterministic model to describe the dynamics of mistletoes in an isolated patch containing an arbitrary number of trees. We establish concrete criterions, expressed in terms of the model parameters, for mistletoes establishing in this area. We conduct numerical simulations based on a field study to reinforce and expand our results.

Plant-herbivore interactions mediated by plant toxicity.

We explore the impact of plant toxicity on the dynamics of a plant-herbivore interaction, such as that of a mammalian browser and its plant forage species, by studying a mathematical model that includes a toxin-determined functional response. In this functional response, the traditional Holling Type 2 response is modified to include the negative effect of toxin on herbivore growth, which can overwhelm the positive effect of biomass ingestion at sufficiently high plant toxicant concentrations. Two types of consumption decisions of the herbivore are considered. One of these (Case 1) incorporates the adaptation of the herbivore to control its rate of consumption of plant items when that is likely to lead to levels of toxicity that more than offset the marginal gain to the herbivore of consuming more plant biomass, while the other (Case 2) simply assumes that, although the herbivore's rate of ingestion of plant biomass is negatively affected by increasing ingestion of toxicant relative to the load it can safely deal with, the herbivore is not able to prevent detrimental or even lethal levels of toxicant intake. A primary result of this work is that these differences in behavior lead to dramatically different outcomes, summarized in bifurcation diagrams. In Case 2, a wide variety of dynamics may occur due to the interplay of Holling Type 2 dynamics and the effect of the plant toxicant. These dynamics include the occurrence of bistability, in which both a periodic solution and the herbivore-extinction equilibrium are attractors, as well the possibility of a homoclinic bifurcation. Whether the herbivore goes to extinction in the bistable case depends on initial conditions of herbivore and plant biomasses. For relatively low herbivore resource acquisition rates, the toxicant effect increases the likelihood of 'paradox of enrichment' type limit cycle oscillations, but at higher resource acquisition rates, the toxicant may decrease the likelihood of these cycles.

Eradicating vector-borne diseases via age-structured culling.

We derive appropriate mathematical models to assess the effectiveness of culling as a tool to eradicate vector-borne diseases. The model, focused on the culling strategies determined by the stages during the development of the vector, becomes either a system of autonomous delay differential equations with impulses (in the case where the adult vector is subject to culling) or a system of nonautonomous delay differential equations where the time-varying coefficients are determined by the culling times and rates (in the case where only the immature vector is subject to culling). Sufficient conditions are derived to ensure eradication of the disease, and simulations are provided to compare the effectiveness of larvicides and insecticide sprays for the control of West Nile virus. We show that eradication of vector-borne diseases is possible by culling the vector at either the immature or the mature phase, even though the size of the vector is oscillating and above a certain level.

Modeling spatial spread of west nile virus and impact of directional dispersal of birds.

A patchy model for the spatial spread of West Nile virus is formulated and analyzed. The basic reproduction number is calculated and com- pared for different long-range dispersal patterns of birds, and simulations are carried out to demonstrate discontinuous or jump spatial spread of the virus when the birds' long-range dispersal dominates the nearest neighborhood interaction and diffusion of mosquitoes and birds.