{"title":"The power of Allee effects: inducing multistability and oscillations in a stoichiometric producer-herbivore system.","authors":"Zhiwei Zhu, Tao Feng","doi":"10.1007/s00285-025-02197-x","DOIUrl":"https://doi.org/10.1007/s00285-025-02197-x","url":null,"abstract":"<p><p>Understanding producer-herbivore dynamics is fundamental for maintaining ecosystem stability and biodiversity. This study proposes a novel stoichiometric producer-herbivore model that incorporates positive density dependence induced by demographic factors. We conduct a rigorous mathematical analysis of the proposed model, covering well-posedness, nullcline analysis, and system stability. This analysis is expanded through numerical bifurcation analysis to explore the effects of critical biological parameters, including light intensity, on producer-herbivore interactions. Our findings reveal that variations in the severity of the Allee effect significantly influence these interactions, driving multistability and periodic oscillations. Severe Allee effects lead to complex dynamics, including four forms of bistability and three forms of tristability. Severe Allee effects can also lead to the extinction of both producer and herbivore populations due to positive density dependence. Intermediate levels of parameters such as light intensity, producer growth rate, herbivore loss rate, saturation levels of the Allee effect, total phosphorus, and sufficiently high production efficiency can lead to system instability and oscillations. Conversely, in scenarios with low-severity Allee effects, the system shows relatively simpler dynamics, with three types of bistability. Low producer growth rate and herbivore loss rate, moderate saturation levels of the Allee effect, light intensity, and sufficiently high herbivore production efficiency and total phosphorus levels can induce periodic oscillations. These findings emphasize the importance of managing Allee effect severity in conservation efforts to sustain biodiversity and prevent undesirable state transitions.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 4","pages":"35"},"PeriodicalIF":2.2,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143517173","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Dynamics of an epidemic model arising in a spatial segregation control strategy.","authors":"Zhiguo Wang, Hua Nie, Sanyi Tang","doi":"10.1007/s00285-025-02195-z","DOIUrl":"https://doi.org/10.1007/s00285-025-02195-z","url":null,"abstract":"<p><p>In this paper, we propose a free boundary problem to model the spread of an epidemic by introducing a spatial segregation control strategy. The model consists of two coupled reaction-diffusion equations along with an ordinary differential equation, while the free boundary is described by an integro-differential equation. The results reveal a trichotomy in which the epidemic can shrink, reach equilibrium, or expand spatially. Moreover, we establish the final size of the cumulative number of infected populations and characterize the threshold phenomenon of epidemic outbreak using the principal eigenvalue of an elliptic operator. Additionally, we apply this model to simulate the spatial spread of the COVID-19 epidemic in Xi'an, China, during 2021-2022. This study provides valuable model references for dynamically designing spatial isolation control strategies for newly emerging major infectious diseases.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 4","pages":"34"},"PeriodicalIF":2.2,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143505777","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convective stability of the critical waves of an FKPP-type model for self-organized growth.","authors":"Florian Kreten","doi":"10.1007/s00285-025-02189-x","DOIUrl":"10.1007/s00285-025-02189-x","url":null,"abstract":"<p><p>We construct the traveling wave solutions of an FKPP growth process of two densities of particles, and prove that the critical traveling waves are locally stable in a space where the perturbations can grow exponentially at the back of the wave. The considered reaction-diffusion system was introduced by Hannezo et al. (Cell 171(1):242-255, 2017) in the context of branching morphogenesis: active, branching particles accumulate inactive particles, which do not react. Thus, the system features a continuum of steady state solutions, complicating the analysis. We adopt a result by Faye and Holzer (J Differ Equ 269(9):6559-6601, 2020) for proving the stability of the critical traveling waves, and modify the semi-group estimates to spaces with unbounded weights. We use a Feynman-Kac formula to get an exponential a priori estimate for the tail of the PDE, a novel and simple approach.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"33"},"PeriodicalIF":2.2,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11832597/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143442624","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Basic concepts for the Kermack and McKendrick model with static heterogeneity.","authors":"Hisashi Inaba","doi":"10.1007/s00285-025-02187-z","DOIUrl":"10.1007/s00285-025-02187-z","url":null,"abstract":"<p><p>In this paper, we consider the infection-age-dependent Kermack-McKendrick model, where host individuals are distributed in a continuous state space. To provide a mathematical foundation for the heterogeneous model, we develop a <math><msup><mi>L</mi> <mn>1</mn></msup> </math> -framework to formulate basic epidemiological concepts. First, we show the mathematical well-posedness of the basic model under appropriate conditions allowing for unbounded structural variables in an unbounded domain. Next, we define the basic reproduction number and prove pandemic threshold results. We then present a systematic procedure to compute the effective reproduction number and the herd immunity threshold. Finally, we give some illustrative examples and concrete results by using the separable mixing assumption.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"32"},"PeriodicalIF":2.2,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143442621","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Description of chemical systems by means of response functions.","authors":"E Franco, B Kepka, J J L Velázquez","doi":"10.1007/s00285-025-02191-3","DOIUrl":"10.1007/s00285-025-02191-3","url":null,"abstract":"<p><p>In this paper we introduce a formalism that allows to describe the response of a part of a biochemical system in terms of renewal equations. In particular, we examine under which conditions the interactions between the different parts of a chemical system, described by means of linear ODEs, can be represented in terms of renewal equations. We show also how to apply the formalism developed in this paper to some particular types of linear and non-linear ODEs, modelling some biochemical systems of interest in biology (for instance, some time-dependent versions of the classical Hopfield model of kinetic proofreading). We also analyse some of the properties of the renewal equations that we are interested in, as the long-time behaviour of their solution. Furthermore, we prove that the kernels characterising the renewal equations derived by biochemical system with reactions that satisfy the detail balance condition belong to the class of completely monotone functions.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"31"},"PeriodicalIF":2.2,"publicationDate":"2025-02-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11830649/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143434299","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The effect of pathogens from environmental breeding and accumulative release by the infected individuals on spread dynamics of a SIRP epidemic model.","authors":"Ning Wang, Long Zhang, Zhidong Teng","doi":"10.1007/s00285-025-02194-0","DOIUrl":"10.1007/s00285-025-02194-0","url":null,"abstract":"<p><p>In this paper, a SIRP epidemic model is proposed, wherein the pathogens derive from two ways, i.e., environmental breeding, and accumulative excretion by the infected individuals. The former is characterized by Logistic growth, while the latter is in the form of infinite integral. First, the positivity and ultimate boundedness of solutions are obtained. Second, the basic reproduction number <math><msub><mi>R</mi> <mn>0</mn></msub> </math> is obtained, by which the model is analyzed if either the intrinsic growth rate of environmental pathogens is lower or higher than its clearance rate. For the first case, the disease-free equilibrium is globally asymptotically stable when <math> <mrow><msub><mi>R</mi> <mn>0</mn></msub> <mo><</mo> <mn>1</mn></mrow> </math> , while the endemic equilibrium is globally asymptotically stable when <math> <mrow><msub><mi>R</mi> <mn>0</mn></msub> <mo>></mo> <mn>1</mn></mrow> </math> . Conversely, if the growth rate exceeds the removal rate, the disease-free equilibrium is always unstable, meanwhile, the uniform persistence of the model indicates that there could exist one or multi-endemic equilibria, and it is globally asymptotically stable if the endemic equilibrium is unique. Finally, the theoretical results are illustrated by numerical simulations. We find that the accumulative release of pathogens by the infected individuals in the form of infinite integral is more realistic and consistent with the disease spread than that of linear form by real data.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"30"},"PeriodicalIF":2.2,"publicationDate":"2025-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143411308","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jianing Xu, Huimin Hu, Gregory Ellison, Lili Yu, Christopher C Whalen, Liang Liu
{"title":"Bayesian estimation of transmission networks for infectious diseases.","authors":"Jianing Xu, Huimin Hu, Gregory Ellison, Lili Yu, Christopher C Whalen, Liang Liu","doi":"10.1007/s00285-025-02193-1","DOIUrl":"10.1007/s00285-025-02193-1","url":null,"abstract":"<p><p>Reconstructing transmission networks is essential for identifying key factors like superspreaders and high-risk locations, which are critical for developing effective pandemic prevention strategies. This study presents a Bayesian transmission model that combines genomic and temporal data to reconstruct transmission networks for infectious diseases. The Bayesian transmission model incorporates the latent period and distinguishes between symptom onset and actual infection time, improving the accuracy of transmission dynamics and epidemiological models. It also assumes a homogeneous effective population size among hosts, ensuring that the coalescent process for within-host evolution remains unchanged, even with missing intermediate hosts. This allows the model to effectively handle incomplete samples. Simulation results demonstrate the model's ability to accurately estimate model parameters and transmission networks. Additionally, our proposed hypothesis test can reliably identify direct transmission events. The Bayesian transmission model was applied to a real dataset of Mycobacterium tuberculosis genomes from 69 tuberculosis cases. The estimated transmission network revealed two major groups, each with a superspreader who transmitted M. tuberculosis, either directly or indirectly, to 28 and 21 individuals, respectively. The hypothesis test identified 16 direct transmissions within the estimated network, demonstrating the Bayesian model's advantage over a fixed threshold by providing a more flexible criterion for identifying direct transmissions. This Bayesian approach highlights the critical role of genetic data in reconstructing transmission networks and enhancing our understanding of the origins and transmission dynamics of infectious diseases.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"29"},"PeriodicalIF":2.2,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143400543","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A theoretical and computational study of heteroclinic cycles in Lotka-Volterra systems.","authors":"M C Bortolan, P Kalita, J A Langa, R O Moura","doi":"10.1007/s00285-025-02190-4","DOIUrl":"10.1007/s00285-025-02190-4","url":null,"abstract":"<p><p>In general, global attractors are composed of isolated invariant sets and the connections between them. This structure can possibly be highly complex, encompassing attraction basins, repeller sets and invariant sets that, collectively, form a dynamical landscape. Lotka-Volterra systems have long been pivotal as preliminary models for dynamics in complex networks exhibiting pairwise interactions. In scenarios involving Volterra-Lyapunov (VL) stable matrices, the dynamics is simplified in such a way that the positive solutions converge to a single, globally asymptotically stable stationary point as time tends to infinity, thereby excluding the existence of periodic solutions. In this work, we conduct a systematic study on the emergence of heteroclinic cycles within Lotka-Volterra systems characterized by Volterra-Lyapunov stable matrices. Although VL stability of the matrix implies that <math><mi>ω</mi></math> -limit sets of solutions are always stationary points, our analysis of <math><mi>α</mi></math> -limit sets reveals finite sets of stationary points interconnected by global trajectories, forming structures referred to as heteroclinic cycles. Our findings indicate that even within the framework of VL stable matrices, such structures are more prevalent than previously thought in literature, driven by the interplay between the symmetric and antisymmetric components of the model matrix. This understanding also reinforces our comprehension of the classical three-dimensional May-Leonard model, which is known to be the unique case exhibiting heteroclinic cycle within the VL framework in dimension three, while also pointing to a surprising richness in the dynamics of these structures in higher dimensions.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"28"},"PeriodicalIF":2.2,"publicationDate":"2025-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11805806/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143371431","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Feedback-delay dependence of the stability of cluster periodic orbits in populations of degrade-and-fire oscillators with common activator.","authors":"Bastien Fernandez, Matteo Tanzi","doi":"10.1007/s00285-024-02169-7","DOIUrl":"10.1007/s00285-024-02169-7","url":null,"abstract":"<p><p>Feedback delay has been identified as a key ingredient in the quorum sensing synchronization of synthetic gene oscillators. While this influence has been evidenced at the theoretical level in a simplified system of degrade-and-fire oscillators coupled via a common activator protein, full mathematical certifications remained to be provided. Here, we prove from a rigorous mathematical viewpoint that, for the very same model, the synchronized degrade-and-fire oscillations are 1/ unstable with respect to out-of-sync perturbations in absence of delay, and 2/ are otherwise asymptotically stable in presence of delay, no matter how small is its amplitude. To that goal, we proceed to an extensive study of the population dynamics in this system, which in particular identifies the mechanisms of, and related criteria for, the delay-dependent stability of periodic orbits with respect to out-of-sync perturbations. As an additional outcome, the analysis also reveals that, depending on the parameters, multiple stable partially synchronized periodic orbits can coexist with the fully synchronized one.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"27"},"PeriodicalIF":2.2,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11802712/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143257085","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Evolutionary rescue of bacterial populations by heterozygosity on multicopy plasmids.","authors":"Ian Dewan, Hildegard Uecker","doi":"10.1007/s00285-025-02182-4","DOIUrl":"10.1007/s00285-025-02182-4","url":null,"abstract":"<p><p>Bacterial plasmids and other extrachromosomal DNA elements frequently carry genes with important fitness effects for their hosts. Multicopy plasmids can additionally carry distinct alleles of host-fitness-relevant genes on different plasmid copies, allowing for heterozygosity not possible for loci on haploid chromosomes. Plasmid-mediated heterozygosity may increase the fitness of bacterial cells in circumstances where there is an advantage to having multiple distinct alleles (heterozyogote advantage); however, plasmid-mediated heterozygosity is also subject to constant loss due to random segregation of plasmid copies on cell division. We analyze a multitype branching process model to study the evolution and maintenance of plasmid-mediated heterozygosity under a heterozygote advantage. We focus on an evolutionary rescue scenario in which a novel mutant allele on a plasmid must be maintained together with the wild-type allele to allow population persistance (although our results apply more generally to the maintenance of heterozygosity due to heterozygote advantage). We determine the probability of rescue and derive an analytical expression for the threshold on the fitness of heterozygotes required to overcome segregation and make rescue possible; this threshold decreases with increasing plasmids copy number. We further show that the formation of cointegrates from the fusion of plasmid copies increases the probability of rescue. Overall, our results provide a rigorous quantitative assessment of the conditions under which bacterial populations can adapt to multiple stressors through plasmid-mediated heterozygosity. Many of the results are furthermore applicable to the related problem of the maintenance of incompatible plasmids in the same cell under selection for both.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"90 3","pages":"26"},"PeriodicalIF":2.2,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11799102/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143256898","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}