Journal of Theoretical Biology最新文献

Insect outbreaks can cause large scale defoliation of forest trees or destruction of crops, leading to ecosystem degradation and economic losses. Some outbreaks occur simultaneously across large geographic scales and some outbreaks occur periodically every few years across space. Parasitoids are a natural enemy of these defoliators and could help mitigate these pest outbreaks. A holistic understanding of the host-parasitoid interactions in a spatial context would thus enhance our ability to understand, predict and prevent these outbreaks. We use a discrete time deterministic model of the host parasitoid system with populations migrating between 2 patches to elucidate features of spatial host outbreaks. We show that whenever populations persist indefinitely, host outbreaks in both patches can occur alternatively (out of phase) at low migration between patches whereas host outbreaks always occur simultaneously (in phase) in both patches at high migration between patches. We show that our results are robust across a large range of parameters across different modelling approaches used typically to model intraspecific competition among hosts and parasitism, in the host-parasitoid literature. We give an analytical expression for the period of oscillations when the migration is low i.e., when host outbreaks in both patches are out of phase, show it is in agreement with numerical results. We end our paper by showing that we get the same results whether we include the biologically rooted formulations from May et al. (1981) or a general cellular automata model with qualitative rules.

Biomathematical models of fatigue capture the physiology of sleep/wake regulation and circadian rhythmicity to predict changes in neurobehavioral functioning over time. We used a biomathematical model of fatigue linked to the adenosinergic neuromodulator/receptor system in the brain as a framework to predict sleep inertia, that is, the transient neurobehavioral impairment experienced immediately after awakening. Based on evidence of an adenosinergic basis for sleep inertia, we expanded the biomathematical model with novel differential equations to predict the propensity for sleep inertia during sleep and its manifestation after awakening. Using datasets from large laboratory studies of sleep loss and circadian misalignment, we calibrated the model by fitting just two new parameters and then validated the model’s predictions against independent data. The expanded model was found to predict the magnitude and time course of sleep inertia with generally high accuracy. Analysis of the model’s dynamics revealed a bifurcation in the predicted manifestation of sleep inertia in sustained sleep restriction paradigms, which reflects the observed escalation of the magnitude of sleep inertia in scenarios with sleep restriction to less than ∼ 4 h per day. Another emergent property of the model involves a rapid increase in the predicted propensity for sleep inertia in the early part of sleep followed by a gradual decline in the later part of the sleep period, which matches what would be expected based on the adenosinergic regulation of non-rapid eye movement (NREM) sleep and its known influence on sleep inertia. These dynamic behaviors provide confidence in the validity of our approach and underscore the predictive potential of the model. The expanded model provides a useful tool for predicting sleep inertia and managing impairment in 24/7 settings where people may need to perform critical tasks immediately after awakening, such as on-demand operations in safety and security, emergency response, and health care.

Natural systems show heterogeneous patchy distributions of vegetation over large landscapes. Reaction–diffusion systems can demonstrate such heterogeneity of species distributions. Here, we analyse a reaction–diffusion model of plant–herbivore interactions in two-dimensional space to illustrate non-homogeneous distributions of plants and herbivores. The non-spatial system shows bottom-up control, where herbivore density is low under low and high primary productivity but increased at intermediate productivity. In addition, the non-spatial system provides bistability between a dense vegetation state devoid of herbivores and a coexisting state of plants and herbivores. In the spatiotemporal model, we give analytical conditions of occurring diffusion-driven (Turing) instability, where a novel point in our model is the relative dispersal of herbivores, which represents the movement of herbivores from a higher to a lower vegetation state in addition to the self-diffusion of both species. It is shown that heterogeneity in the population distribution does not occur if the relative dispersal of herbivores is low, but it appears in the opposite case. Due to bistability in the underlying non-spatial system, the spatiotemporal model produces initial value-dependent patterns. The two initial values make different patterns despite having the same primary productivity and relative dispersal rate. As productivity increases with a given relative herbivore dispersal, pattern transition occurs from a blend of stripes and spots of low vegetation state to a predominantly low-density vegetation state with smaller patches of densely vegetated states with one initial value. On the contrary, a discernible change in vegetation patterns from cold spots in the dense vegetation to hot stripes in the primarily low-vegetated state is noticed under the other initial population value. Furthermore, the population distributions of plants and herbivores in the entire domain after a long period are heterogeneous for both initial values, provided the relative herbivore dispersal is substantial. We estimated mean population densities to observe species fitness in the whole domain under variable productivity. When productivity is high, the mean population density of plants may go up or down, depending on the herbivore’s relative dispersal rate. In contrast to the bottom-up control dynamics of the non-spatial system, the system exhibits a top-down control under high relative dispersal, where the herbivore regulates vegetation growth under high productivity. On the other hand, herbivores are extinct under high productivity if the relative dispersal is low.

Fungal pellets are hierarchical systems that can be found in an ample variety of applications. Modeling transport phenomena in this type of systems is a challenging but necessary task to provide knowledge-based processes that improve the outcome of their biotechnological applications. In this work, an upscaled model for total mass and momentum transport in fungal pellets is implemented and analyzed, using elements of the volume averaging and adjoint homogenization methods departing from the governing equations at the microscale in the intracellular and extracellular phases. The biomass is assumed to be composed of a non-Newtonian fluid and the organelles impervious to momentum transport are modeled as a rigid solid phase. The upscaled equations contain effective-medium coefficients, which are predicted from the solution of adjoint closure problems in a three-dimensional periodic domains representative of the microstructure. The construction of these domains was performed for Laccaria trichodermophora based on observations of actual biological structures. The upscaled model was validated with direct numerical simulations in homogeneous portions of the pellets core. It is shown that no significant differences are observed when the dolipores are open or closed to fluid flow. By comparing the predictions of the average velocity in the extracellular phase resulting from the upscaled model with those from the classical Darcy equation (i.e., assuming that the biomass is a solid phase) the contribution of the intracellular fluid phase was evidenced. This work sets the foundations for further studies dedicated to transport phenomena in this type of systems.