Prey group defense and hunting cooperation among generalist-predators induce complex dynamics: a mathematical study.

Group defense in prey and hunting cooperation in predators are two important ecological phenomena and can occur concurrently. In this article, we consider cooperative hunting in generalist predators and group defense in prey under a mathematical framework to comprehend the enormous diversity the model could capture. To do so, we consider a modified Holling-Tanner model where we implement Holling type IV functional response to characterize grazing pattern of predators where prey species exhibit group defense. Additionally, we allow a modification in the attack rate of predators to quantify the hunting cooperation among them. The model admits three boundary equilibria and up to three coexistence equilibrium points. The geometry of the nontrivial prey and predator nullclines and thus the number of coexistence equilibria primarily depends on a specific threshold of the availability of alternative food for predators. We use linear stability analysis to determine the types of hyperbolic equilibrium points and characterize the non-hyperbolic equilibrium points through normal form and center manifold theory. Change in the model parameters leading to the occurrences of a series of local bifurcations from non-hyperbolic equilibrium points, namely, transcritical, saddle-node, Hopf, cusp and Bogdanov-Takens bifurcation; there are also occurrences of global bifurcations such as homoclinic bifurcation and saddle-node bifurcation of limit cycles. We observe two interesting closed 'bubble' form induced by global bifurcations due to change in the strength of hunting cooperation and the availability of alternative food for predators. A three dimensional bifurcation diagram, concerning the original system parameters, captures how the alternation in model formulation induces gradual changes in the bifurcation scenarios. Our model highlights the stabilizing effects of group or gregarious behaviour in both prey and predator, hence supporting the predator-herbivore regulation hypothesis. Additionally, our model highlights the occurrence of "saltatory equilibria" in ecological systems and capture the dynamics observed for lion-herbivore interactions.

A two-timescale model of plankton-oxygen dynamics predicts formation of oxygen minimum zones and global anoxia.

Decline of the dissolved oxygen in the ocean is a growing concern, as it may eventually lead to global anoxia, an elevated mortality of marine fauna and even a mass extinction. Deoxygenation of the ocean often results in the formation of oxygen minimum zones (OMZ): large domains where the abundance of oxygen is much lower than that in the surrounding ocean environment. Factors and processes resulting in the OMZ formation remain controversial. We consider a conceptual model of coupled plankton-oxygen dynamics that, apart from the plankton growth and the oxygen production by phytoplankton, also accounts for the difference in the timescales for phyto- and zooplankton (making it a "slow-fast system") and for the implicit effect of upper trophic levels resulting in density dependent (nonlinear) zooplankton mortality. The model is investigated using a combination of analytical techniques and numerical simulations. The slow-fast system is decomposed into its slow and fast subsystems. The critical manifold of the slow-fast system and its stability is then studied by analyzing the bifurcation structure of the fast subsystem. We obtain the canard cycles of the slow-fast system for a range of parameter values. However, the system does not allow for persistent relaxation oscillations; instead, the blowup of the canard cycle results in plankton extinction and oxygen depletion. For the spatially explicit model, the earlier works in this direction did not take into account the density dependent mortality rate of the zooplankton, and thus could exhibit Turing pattern. However, the inclusion of the density dependent mortality into the system can lead to stationary Turing patterns. The dynamics of the system is then studied near the Turing bifurcation threshold. We further consider the effect of the self-movement of the zooplankton along with the turbulent mixing. We show that an initial non-uniform perturbation can lead to the formation of an OMZ, which then grows in size and spreads over space. For a sufficiently large timescale separation, the spread of the OMZ can result in global anoxia.

Does mutual interference stabilize prey-predator model with Bazykin-Crowley-Martin trophic function?

We investigated a system of ordinary differential equations that describes the dynamics of prey and predator populations, taking into account the Allee effect affecting the reproduction of the predator population, and mutual interference amongst predators, which is modeled with the Bazykin-Crowley-Martin (BCM) trophic function. Bifurcation analysis revealed a rich spectrum of bifurcations occurring in the system. In particular, analytical conditions for the saddle-node, Hopf, cusp, and Bogdanov-Takens bifurcations were derived for the model parameters, quantifying the strength of the predator interference, the Allee effect, and the predation efficiency. Numerical simulations verify and illustrate the analytical findings. The main purpose of the study was to test whether the mutual interference in the model with BCM trophic function provides a stabilizing or destabilizing effect on the system dynamics. The obtained results suggest that the model demonstrates qualitatively the same pattern concerning varying the interference strength as other predator-dependent models: both low and very high interference levels increase the risk of predator extinction, while moderate interference has a favorable effect on the stability and resilience of the prey-predator system.

Combinatorial Cooperativity in miR200-Zeb Feedback Network can Control Epithelial-Mesenchymal Transition.

Carcinomas often utilize epithelial-mesenchymal transition (EMT) programs for cancer progression and metastasis. Numerous studies report SNAIL-induced miR200/Zeb feedback circuit as crucial in regulating EMT by placing cancer cells in at least three phenotypic states, viz. epithelial (E), hybrid (h-E/M), mesenchymal (M), along the E-M phenotypic spectrum. However, a coherent molecular-level understanding of how such a tiny circuit controls carcinoma cell entrance into and residence in various states is lacking. Here, we use molecular binding data and mathematical modeling to report that the miR200/Zeb circuit can essentially utilize combinatorial cooperativity to control E-M phenotypic plasticity. We identify minimal combinatorial cooperativities that give rise to E, h-E/M, and M phenotypes. We show that disrupting a specific number of miR200 binding sites on Zeb as well as Zeb binding sites on miR200 can have phenotypic consequences-the circuit can dynamically switch between two (E, M) and three (E, h-E/M, M) phenotypes. Further, we report that in both SNAIL-induced and SNAIL knock-out miR200/Zeb circuits, cooperative transcriptional feedback on Zeb as well as Zeb translation inhibition due to miR200 are essential for the occurrence of intermediate h-E/M phenotype. Finally, we demonstrate that SNAIL can be dispensable for EMT, and in the absence of SNAIL, the transcriptional feedback can control cell state transition from E to h-E/M, to M state. Our results thus highlight molecular-level regulation of EMT in miR200/Zeb circuit and we expect these findings to be crucial to future efforts aiming to prevent EMT-facilitated dissemination of carcinomas.

Hunting cooperation in a prey-predator model with maturation delay.

We investigate the dynamics of a prey-predator model with cooperative hunting among specialist predators and maturation delay in predator growth. First, we consider a model without delay and explore the effect of hunting time on the coexistence of predator and their prey. When the hunting time is long enough and the cooperation rate among predators is weak, prey and predator species tend to coexist. Furthermore, we observe the occurrences of a series of bifurcations that depend on the cooperation rate and the hunting time. Second, we introduce a maturation delay for predator growth and analyse its impact on the system's dynamics. We find that as the delay becomes larger, predator species become more likely to go extinct, as the long maturation delay hinders the growth of the predator population. Our numerical exploration reveals that the delay causes shifts in both the bifurcation curves and bifurcation thresholds of the non-delayed system.

Delay epidemic models determined by latency, infection, and immunity duration.

We propose new single and two-strain epidemic models represented by systems of delay differential equations and based on the number of newly exposed individuals. Transitions between exposed, infectious, recovered, and back to susceptible compartments are determined by the corresponding time delays. Existence and positiveness of solutions are proved. Reduction of delay differential equations to integral equations allows the analysis of stationary solutions and their stability. In the case of two strains, they compete with each other, and the strain with a larger individual basic reproduction number dominates the other one. However, if the basic reproduction number exceeds some critical values, stationary solution loses its stability resulting in periodic time oscillations. In this case, both strains are present and their dynamics is not completely determined by the basic reproduction numbers but also by other parameters. The results of the work are illustrated by comparison with data on seasonal influenza.

Understanding the implications of under-reporting, vaccine efficiency and social behavior on the post-pandemic spread using physics informed neural networks: A case study of China.

In late 2019, the emergence of COVID-19 in Wuhan, China, led to the implementation of stringent measures forming the zero-COVID policy aimed at eliminating transmission. Zero-COVID policy basically aimed at completely eliminating the transmission of COVID-19. However, the relaxation of this policy in late 2022 reportedly resulted in a rapid surge of COVID-19 cases. The aim of this work is to investigate the factors contributing to this outbreak using a new SEIR-type epidemic model with time-dependent level of immunity. Our model incorporates a time-dependent level of immunity considering vaccine doses administered and time-post-vaccination dependent vaccine efficacy. We find that vaccine efficacy plays a significant role in determining the outbreak size and maximum number of daily infected. Additionally, our model considers under-reporting in daily cases and deaths, revealing their combined effects on the outbreak magnitude. We also introduce a novel Physics Informed Neural Networks (PINNs) approach which is extremely useful in estimating critical parameters and helps in evaluating the predictive capability of our model.

Analytical detection of stationary and dynamic patterns in a prey-predator model with reproductive Allee effect in prey growth.

Allee effect in population dynamics has a major impact in suppressing the paradox of enrichment through global bifurcation, and it can generate highly complex dynamics. The influence of the reproductive Allee effect, incorporated in the prey's growth rate of a prey-predator model with Beddington-DeAngelis functional response, is investigated here. Preliminary local and global bifurcations are identified of the temporal model. Existence and non-existence of heterogeneous steady-state solutions of the spatio-temporal system are established for suitable ranges of parameter values. The spatio-temporal model satisfies Turing instability conditions, but numerical investigation reveals that the heterogeneous patterns corresponding to unstable Turing eigenmodes act as a transitory pattern. Inclusion of the reproductive Allee effect in the prey population has a destabilising effect on the coexistence equilibrium. For a range of parameter values, various branches of stationary solutions including mode-dependent Turing solutions and localized pattern solutions are identified using numerical bifurcation technique. The model is also capable to produce some complex dynamic patterns such as travelling wave, moving pulse solution, and spatio-temporal chaos for certain range of parameters and diffusivity along with appropriate choice of initial conditions. Judicious choices of parametrization for the Beddington-DeAngelis functional response help us to infer about the resulting patterns for similar prey-predator models with Holling type-II functional response and ratio-dependent functional response.

An epidemic model with time delays determined by the infectivity and disease durations.

We propose an epidemiological model with distributed recovery and death rates. It represents an integrodifferential system of equations for susceptible, exposed, infectious, recovered and dead compartments. This model can be reduced to the conventional ODE model under the assumption that recovery and death rates are uniformly distributed in time during disease duration. Another limiting case, where recovery and death rates are given by the delta-function, leads to a new point-wise delay model with two time delays corresponding to the infectivity period and disease duration. Existence and positiveness of solutions for the distributed delay model and point-wise delay model are proved. The basic reproduction number and the final size of the epidemic are determined. Both, the ODE model and the delay models are used to describe COVID-19 epidemic progression. The delay model gives a better approximation of the Omicron data than the conventional ODE model from the point of view of parameter estimation.

Saddle-node bifurcation of periodic orbit route to hidden attractors.

Hidden attractors are present in many nonlinear dynamical systems and are not associated with equilibria, making them difficult to locate. Recent studies have demonstrated methods of locating hidden attractors, but the route to these attractors is still not fully understood. In this Research Letter, we present the route to hidden attractors in systems with stable equilibrium points and in systems without any equilibrium points. We show that hidden attractors emerge as a result of the saddle-node bifurcation of stable and unstable periodic orbits. Real-time hardware experiments were performed to demonstrate the existence of hidden attractors in these systems. Despite the difficulties in identifying suitable initial conditions from the appropriate basin of attraction, we performed experiments to detect hidden attractors in nonlinear electronic circuits. Our results provide insights into the generation of hidden attractors in nonlinear dynamical systems.

A mathematical modeling technique to understand the role of decoy receptors in ligand-receptor interaction.

The ligand-receptor interaction is fundamental to many cellular processes in eukaryotic cells such as cell migration, proliferation, adhesion, signaling and so on. Cell migration is a central process in the development of organisms. Receptor induced chemo-tactic sensitivity plays an important role in cell migration. However, recently some receptors identified as decoy receptors, obstruct some mechanisms of certain regular receptors. DcR3 is one such important decoy receptor, generally found in glioma cell, RCC cell and many various malignant cells which obstruct some mechanism including apoptosis cell-signaling, cell inflammation, cell migration. The goal of our work is to mathematically formulate ligand-receptor interaction induced cell migration in the presence of decoy receptors. We develop here a novel mathematical model, consisting of four coupled partial differential equations which predict the movement of glioma cells due to the reaction-kinetic mechanism between regular receptors CD95, its ligand CD95L and decoy receptors DcR3 as obtained in experimental results. The aim is to measure the number of cells in the chamber's filter for different concentrations of ligand in presence of decoy receptors and compute the distance travelled by the cells inside filter due to cell migration. Using experimental results, we validate our model which suggests that the concentration of ligands plays an important role in cell migration.

An age-dependent immuno-epidemiological model with distributed recovery and death rates.

The work is devoted to a new immuno-epidemiological model with distributed recovery and death rates considered as functions of time after the infection onset. Disease transmission rate depends on the intra-subject viral load determined from the immunological submodel. The age-dependent model includes the viral load, recovery and death rates as functions of age considered as a continuous variable. Equations for susceptible, infected, recovered and dead compartments are expressed in terms of the number of newly infected cases. The analysis of the model includes the proof of the existence and uniqueness of solution. Furthermore, it is shown how the model can be reduced to age-dependent SIR or delay model under certain assumptions on recovery and death distributions. Basic reproduction number and final size of epidemic are determined for the reduced models. The model is validated with a COVID-19 case data. Modelling results show that proportion of young age groups can influence the epidemic progression since disease transmission rate for them is higher than for other age groups.

An Epidemic Model with Time-Distributed Recovery and Death Rates.

A compartmental epidemiological model with distributed recovery and death rates is proposed. In some particular cases, the model can be reduced to the conventional SIR model. However, in general, the dynamics of epidemic progression in this model is different. Distributed recovery and death rates are evaluated from COVID-19 data. The model is validated by the epidemiological data for different countries, and it shows better agreement with the data than the SIR model. The time-dependent disease transmission rate is estimated.

Bifurcation analysis of the predator-prey model with the Allee effect in the predator.

The use of predator-prey models in theoretical ecology has a long history, and the model equations have largely evolved since the original Lotka-Volterra system towards more realistic descriptions of the processes of predation, reproduction and mortality. One important aspect is the recognition of the fact that the growth of a population can be subject to an Allee effect, where the per capita growth rate increases with the population density. Including an Allee effect has been shown to fundamentally change predator-prey dynamics and strongly impact species persistence, but previous studies mostly focused on scenarios of an Allee effect in the prey population. Here we explore a predator-prey model with an ecologically important case of the Allee effect in the predator population where it occurs in the numerical response of predator without affecting its functional response. Biologically, this can result from various scenarios such as a lack of mating partners, sperm limitation and cooperative breeding mechanisms, among others. Unlike previous studies, we consider here a generic mathematical formulation of the Allee effect without specifying a concrete parameterisation of the functional form, and analyse the possible local bifurcations in the system. Further, we explore the global bifurcation structure of the model and its possible dynamical regimes for three different concrete parameterisations of the Allee effect. The model possesses a complex bifurcation structure: there can be multiple coexistence states including two stable limit cycles. Inclusion of the Allee effect in the predator generally has a destabilising effect on the coexistence equilibrium. We also show that regardless of the parametrisation of the Allee effect, enrichment of the environment will eventually result in extinction of the predator population.

Oscillations and Pattern Formation in a Slow-Fast Prey-Predator System.

We consider the properties of a slow-fast prey-predator system in time and space. We first argue that the simplicity of the prey-predator system is apparent rather than real and there are still many of its hidden properties that have been poorly studied or overlooked altogether. We further focus on the case where, in the slow-fast system, the prey growth is affected by a weak Allee effect. We first consider this system in the non-spatial case and make its comprehensive study using a variety of mathematical techniques. In particular, we show that the interplay between the Allee effect and the existence of multiple timescales may lead to a regime shift where small-amplitude oscillations in the population abundances abruptly change to large-amplitude oscillations. We then consider the spatially explicit slow-fast prey-predator system and reveal the effect of different timescales on the pattern formation. We show that a decrease in the timescale ratio may lead to another regime shift where the spatiotemporal pattern becomes spatially correlated, leading to large-amplitude oscillations in spatially average population densities and potential species extinction.

Turing instability in an economic-demographic dynamical system may lead to pattern formation on a geographical scale.

Spatial distribution of the human population is distinctly heterogeneous, e.g. showing significant difference in the population density between urban and rural areas. In the historical perspective, i.e. on the timescale of centuries, the emergence of densely populated areas at their present locations is widely believed to be linked to more favourable environmental and climatic conditions. In this paper, we challenge this point of view. We first identify a few areas at different parts of the world where the environmental conditions (quantified by the temperature, precipitation and elevation) show a relatively small variation in space on the scale of thousands of kilometres. We then examine the population distribution across those areas to show that, in spite of the approximate homogeneity of the environment, it exhibits a significant variation revealing a nearly periodic spatial pattern. Based on this apparent disagreement, we hypothesize that there may exist an inherent mechanism that may lead to pattern formation even in a uniform environment. We consider a mathematical model of the coupled demographic-economic dynamics and show that its spatially uniform, locally stable steady state can give rise to a periodic spatial pattern due to the Turing instability, the spatial scale of the emerging pattern being consistent with observations. Using numerical simulations, we show that, interestingly, the emergence of the Turing patterns may eventually lead to the system collapse.

Pattern Formation in a Three-Species Cyclic Competition Model.

In nature, different species compete among themselves for common resources and favorable habitat. Therefore, it becomes really important to determine the key factors in maintaining the bio-diversity. Also, some competing species follow cyclic competition in real world where the competitive dominance is characterized by a cyclic ordering. In this paper, we study the formation of a wide variety of spatiotemporal patterns including stationary, periodic, quasi-periodic and chaotic population distributions for a diffusive Lotka-Volterra type three-species cyclic competition model with two different types of cyclic ordering. For both types of cyclic ordering, the temporal dynamics of the corresponding non-spatial system show the extinction of two species through global bifurcations such as homoclinic and heteroclinic bifurcations. For the spatial system, we show that the existence of Turing patterns is possible for a particular cyclic ordering, while it is not the case for the other cyclic ordering through both the analytical and numerical methods. Further, we illustrate an interesting scenario of short-range invasion as opposed to the usual invasion phenomenon over the entire habitat. Also, our study reveals that both the stationary and dynamic population distributions can coexist in different parts of a habitat. Finally, we extend the spatial system by incorporating nonlocal intra-specific competition terms for all the three competing species. Our study shows that the introduction of nonlocality in intra-specific competitions stabilizes the system dynamics by transforming a dynamic population distribution to stationary. Surprisingly, this nonlocality-induced stationary pattern formation leads to the extinction of one species and hence, gives rise to the loss of bio-diversity for intermediate ranges of nonlocality. However, the bio-diversity can be restored for sufficiently large extent of nonlocality.