Pub Date : 2026-01-21Epub Date: 2025-11-07DOI: 10.1016/j.jtbi.2025.112303
Zhao Song , Chen Shen , The Anh Han
Non-binding communication is common in daily life and crucial for fostering cooperation, even though it has no direct payoff consequences. However, despite robust empirical evidence, its evolutionary basis remains poorly understood. Here, we develop a game-theoretic model in which individuals can signal an intention to cooperate before playing a Donation game. Strategies differ in how they respond to these signals, ranging from unconditional to conditional types, with the latter incurring a cognitive cost for deliberation. Through evolutionary analysis, we show that non-binding communication alone cannot sustain cooperation in well-mixed, anonymous populations, consistent with empirical observations. In contrast, structured populations support the emergence of cooperation, with conditional cooperators acting as catalysts that protect unconditional cooperators through context-dependent patterns of cyclic dominance. These findings offer an evolutionary explanation for how non-binding communication promotes cooperation and provide a modelling framework for exploring its effects in diverse social settings.
{"title":"Network reciprocity turns cheap talk into a force for cooperation","authors":"Zhao Song , Chen Shen , The Anh Han","doi":"10.1016/j.jtbi.2025.112303","DOIUrl":"10.1016/j.jtbi.2025.112303","url":null,"abstract":"<div><div>Non-binding communication is common in daily life and crucial for fostering cooperation, even though it has no direct payoff consequences. However, despite robust empirical evidence, its evolutionary basis remains poorly understood. Here, we develop a game-theoretic model in which individuals can signal an intention to cooperate before playing a Donation game. Strategies differ in how they respond to these signals, ranging from unconditional to conditional types, with the latter incurring a cognitive cost for deliberation. Through evolutionary analysis, we show that non-binding communication alone cannot sustain cooperation in well-mixed, anonymous populations, consistent with empirical observations. In contrast, structured populations support the emergence of cooperation, with conditional cooperators acting as catalysts that protect unconditional cooperators through context-dependent patterns of cyclic dominance. These findings offer an evolutionary explanation for how non-binding communication promotes cooperation and provide a modelling framework for exploring its effects in diverse social settings.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112303"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145474952","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}
Pub Date : 2026-01-21Epub Date: 2025-10-17DOI: 10.1016/j.jtbi.2025.112286
Michael Napoli , Rifat Sipahi , Maurizio Porfiri
Thirty years of research into activity rhythms of ant colonies have contributed an improved understanding of this fascinating form of collective behavior. Yet, little is known about how the colony size influences the intensity and tempo of the rhythms. Here, we address this knowledge gap through a two-pronged approach, combining the re-evaluation of published experimental observations on California ant, Temnothorax rudis, with the formulation of a novel mathematical model. From the analysis of published data, we discover that the period of the activity bursts is nearly independent of the colony size, while the number of ants activated during the bursts scales hypometrically with the colony size. In search of the biological mechanisms underpinning this evidence, we put forward a compartmental model consisting of three classes: active, inactive, and refractory ants. The study of the resulting system of nonlinear delay-differential equations explains the emergence of activity rhythms as stable limit cycles. The period of these limit cycles is controlled by the refractory delay for the resting phase, independently of the colony size, but their amplitude allometrically varies with the colony size due to social deactivation. During activity rhythms, ants spontaneously deactivate when socially interacting with active individuals.
{"title":"The role of colony size on activity rhythms of ants","authors":"Michael Napoli , Rifat Sipahi , Maurizio Porfiri","doi":"10.1016/j.jtbi.2025.112286","DOIUrl":"10.1016/j.jtbi.2025.112286","url":null,"abstract":"<div><div>Thirty years of research into activity rhythms of ant colonies have contributed an improved understanding of this fascinating form of collective behavior. Yet, little is known about how the colony size influences the intensity and tempo of the rhythms. Here, we address this knowledge gap through a two-pronged approach, combining the re-evaluation of published experimental observations on California ant, <em>Temnothorax rudis</em>, with the formulation of a novel mathematical model. From the analysis of published data, we discover that the period of the activity bursts is nearly independent of the colony size, while the number of ants activated during the bursts scales hypometrically with the colony size. In search of the biological mechanisms underpinning this evidence, we put forward a compartmental model consisting of three classes: active, inactive, and refractory ants. The study of the resulting system of nonlinear delay-differential equations explains the emergence of activity rhythms as stable limit cycles. The period of these limit cycles is controlled by the refractory delay for the resting phase, independently of the colony size, but their amplitude allometrically varies with the colony size due to social deactivation. During activity rhythms, ants spontaneously deactivate when socially interacting with active individuals.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112286"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145330676","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}
Pub Date : 2026-01-21Epub Date: 2025-10-17DOI: 10.1016/j.jtbi.2025.112300
Jian-Zhong Gao , Feng Zhang , Derek W. Dunn , Hao Wang , K.Charlotte Jandér , Rui-Wu Wang
A fundamental problem in ecology is to understand how mutualisms remain stable. The density-dependent regulations within interacting species potentially impact the persistence of these interspecific relationships. Yet few studies explore such intraspecific regulations’ role in stabilizing mutualisms. In addition, partner species often gain unequal benefits in mutualisms. To what extent such an interspecific asymmetry affects the stability of mutualisms is also poorly understood. We here developed a dynamic model for the asymmetric interaction between plants and their pollinators in nursery mutualisms, considering the intraspecific competition of each mutualist. We found that (i) a mutualism can be stabilized only if both mutualists are subject to the regulation of intraspecific competition; (ii) stabilizing the system also requires that the degree of asymmetry in benefits between mutualists must be limited to a range of ‘tolerance’, which narrows as intraspecific competition increases and even fades away with strong competition within both mutualistic species; (iii) when intraspecific competition within a species increases, the tolerant range is compressed from the side beneficial for it, with thus its partner species gaining relatively more benefit allocation; (iv) if the plant-pollinator interaction initiates from a small host plant population, these host plants must offer pollinators high levels of benefits, that can be subsequently reduced to favor plants once the mutualism has been successfully established. The agreement of empirical data to theoretical predictions suggests model reliability. These results highlight the role of intraspecific competition and the degree of benefit asymmetry between host plants and symbionts in stabilizing mutualisms.
{"title":"Intraspecific competition can stabilize asymmetric nursery pollination mutualisms","authors":"Jian-Zhong Gao , Feng Zhang , Derek W. Dunn , Hao Wang , K.Charlotte Jandér , Rui-Wu Wang","doi":"10.1016/j.jtbi.2025.112300","DOIUrl":"10.1016/j.jtbi.2025.112300","url":null,"abstract":"<div><div>A fundamental problem in ecology is to understand how mutualisms remain stable. The density-dependent regulations within interacting species potentially impact the persistence of these interspecific relationships. Yet few studies explore such intraspecific regulations’ role in stabilizing mutualisms. In addition, partner species often gain unequal benefits in mutualisms. To what extent such an interspecific asymmetry affects the stability of mutualisms is also poorly understood. We here developed a dynamic model for the asymmetric interaction between plants and their pollinators in nursery mutualisms, considering the intraspecific competition of each mutualist. We found that (i) a mutualism can be stabilized only if both mutualists are subject to the regulation of intraspecific competition; (ii) stabilizing the system also requires that the degree of asymmetry in benefits between mutualists must be limited to a range of ‘tolerance’, which narrows as intraspecific competition increases and even fades away with strong competition within both mutualistic species; (iii) when intraspecific competition within a species increases, the tolerant range is compressed from the side beneficial for it, with thus its partner species gaining relatively more benefit allocation; (iv) if the plant-pollinator interaction initiates from a small host plant population, these host plants must offer pollinators high levels of benefits, that can be subsequently reduced to favor plants once the mutualism has been successfully established. The agreement of empirical data to theoretical predictions suggests model reliability. These results highlight the role of intraspecific competition and the degree of benefit asymmetry between host plants and symbionts in stabilizing mutualisms.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112300"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145330722","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}
Pub Date : 2026-01-21Epub Date: 2025-11-04DOI: 10.1016/j.jtbi.2025.112305
Quang Anh Le , Seung Ki Baek
A social norm defines what is good and what is bad in social contexts, as well as what to do based on such assessments. A stable social norm should be maintained against errors committed by its players. In addition, individuals may have different probabilities of errors in following the norm, and a social norm would be unstable if it benefited those who do not follow the norm carefully. In this work, we show that Simple Standing, which has been known to resist errors and mutants successfully, actually exhibits threshold behavior. That is, in a population of individuals playing the donation game according to Simple Standing, the residents can suppress the invasion of mutants with higher error proneness only if the residents’ own error proneness is sufficiently low. Otherwise, the population will be invaded by mutants that commit assessment errors more frequently, and a series of such invasions will eventually undermine the existing social norm. This study suggests that the stability analysis of a social norm may have a different picture if the probability of error itself is regarded as an individual attribute.
{"title":"Threshold behavior of a social norm in response to error proneness","authors":"Quang Anh Le , Seung Ki Baek","doi":"10.1016/j.jtbi.2025.112305","DOIUrl":"10.1016/j.jtbi.2025.112305","url":null,"abstract":"<div><div>A social norm defines what is good and what is bad in social contexts, as well as what to do based on such assessments. A stable social norm should be maintained against errors committed by its players. In addition, individuals may have different probabilities of errors in following the norm, and a social norm would be unstable if it benefited those who do not follow the norm carefully. In this work, we show that Simple Standing, which has been known to resist errors and mutants successfully, actually exhibits threshold behavior. That is, in a population of individuals playing the donation game according to Simple Standing, the residents can suppress the invasion of mutants with higher error proneness only if the residents’ own error proneness is sufficiently low. Otherwise, the population will be invaded by mutants that commit assessment errors more frequently, and a series of such invasions will eventually undermine the existing social norm. This study suggests that the stability analysis of a social norm may have a different picture if the probability of error itself is regarded as an individual attribute.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112305"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145460631","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}
Pub Date : 2026-01-21Epub Date: 2025-10-06DOI: 10.1016/j.jtbi.2025.112289
Rena Hayashi , Akane Hara , Yoh Iwasa
Human papillomavirus (HPV), a DNA virus, causes cervical cancer, which is the most common cancer among Japanese women in their forties. Upon infection, HPV temporarily proliferates but is usually eliminated by the immune system. However, if the virus enters the nuclei of epithelial cells, it can evade immune detection and establish a persistent infection. In this state, HPV inhibits apoptosis and allows genomic mutations to accumulate. Over many years, this can lead to dysplasia, genetic abnormalities, and eventually, invasive cancer with metastasis. While many individuals with persistent HPV infections experience spontaneous remission, a small proportion develop cervical cancer. In this study, we aim to understand the sharp contrast between cervical cancer and other solid tumors (cancers of epithelial tissues). We analyze a mathematical model for stochastic transitions between infection states, where the likelihood of persistent infection is proportional to the cumulative viral load, influenced by viral dynamics, immune effectors, and immune memory. We derive formulas for total cancer incidence, mean age at diagnosis, and age variance. Key parameters were estimated from data using the MCMC method. We conclude that major characteristics of cervical cancer arise from the strong age-dependence of viral genome incorporated into the epithelial tissue — shaped by the human sexual behavior — and from the very high rate of spontaneous remission.
{"title":"Human papillomavirus driving cervical cancer: A mathematical model with persistent infection, cancer progression, and spontaneous remission","authors":"Rena Hayashi , Akane Hara , Yoh Iwasa","doi":"10.1016/j.jtbi.2025.112289","DOIUrl":"10.1016/j.jtbi.2025.112289","url":null,"abstract":"<div><div>Human papillomavirus (HPV), a DNA virus, causes cervical cancer, which is the most common cancer among Japanese women in their forties. Upon infection, HPV temporarily proliferates but is usually eliminated by the immune system. However, if the virus enters the nuclei of epithelial cells, it can evade immune detection and establish a persistent infection. In this state, HPV inhibits apoptosis and allows genomic mutations to accumulate. Over many years, this can lead to dysplasia, genetic abnormalities, and eventually, invasive cancer with metastasis. While many individuals with persistent HPV infections experience spontaneous remission, a small proportion develop cervical cancer. In this study, we aim to understand the sharp contrast between cervical cancer and other solid tumors (cancers of epithelial tissues). We analyze a mathematical model for stochastic transitions between infection states, where the likelihood of persistent infection is proportional to the cumulative viral load, influenced by viral dynamics, immune effectors, and immune memory. We derive formulas for total cancer incidence, mean age at diagnosis, and age variance. Key parameters were estimated from data using the MCMC method. We conclude that major characteristics of cervical cancer arise from the strong age-dependence of viral genome incorporated into the epithelial tissue — shaped by the human sexual behavior — and from the very high rate of spontaneous remission.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112289"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145253781","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}
Pub Date : 2026-01-21Epub Date: 2025-10-30DOI: 10.1016/j.jtbi.2025.112306
Deepak K Pattanaik , Nachieketa K Sharma , Vasudevan Lakshminarayanan
Background
The development of an oxidative stress model has been used to explain the decrease in the ERG’s a-wave amplitude in the presence of iron ions. Reactive oxygen species (ROS), which are produced by iron ions, are thought to reduce photoreceptor hyperpolarization. In Retinitis Pigmentosa (RP), oxygen builds up inside the cells of the rods, where superoxide radicals are produced. This provides an explanation for the decrease in scotopic a-wave amplitude in RP patients.
Methods
The oxidative stress model is applied here since oxygen buildup produces ROS (superoxide radicals) that are comparable to those produced by iron ions. The loss of rod photoreceptors is a further consequence of ROS. A component denoting the fraction of rods present is incorporated into the formula for the a-wave voltage to account for the effect of rod deaths on amplitude.
Results
The ROS generated in RP eyes were able to draw calcium ions, increasing the calcium influx. The voltage vs. time graph can be determined with the aid of a factor that represents the percentage of rods present. The loss of rod photoreceptors and calcium ions both result in a decrease in the amplitude of the a-wave at any time t. The theoretical results compare well with experimental results.
Conclusion
Due to oxidative stress, which causes rod photoreceptors to die and cause tunnel vision and peripheral visual field loss, the decline in scotopic a-wave amplitude in RP patients’ eyes is explained.
{"title":"The effect of oxidative stress on the leading edge of the a-wave in retinitis pigmentosa","authors":"Deepak K Pattanaik , Nachieketa K Sharma , Vasudevan Lakshminarayanan","doi":"10.1016/j.jtbi.2025.112306","DOIUrl":"10.1016/j.jtbi.2025.112306","url":null,"abstract":"<div><h3>Background</h3><div>The development of an oxidative stress model has been used to explain the decrease in the ERG’s <em>a</em>-wave amplitude in the presence of iron ions. Reactive oxygen species (ROS), which are produced by iron ions, are thought to reduce photoreceptor hyperpolarization. In Retinitis Pigmentosa (RP), oxygen builds up inside the cells of the rods, where superoxide radicals are produced. This provides an explanation for the decrease in scotopic <em>a</em>-wave amplitude in RP patients.</div></div><div><h3>Methods</h3><div>The oxidative stress model is applied here since oxygen buildup produces ROS (superoxide radicals) that are comparable to those produced by iron ions. The loss of rod photoreceptors is a further consequence of ROS. A component denoting the fraction of rods present is incorporated into the formula for the <em>a</em>-wave voltage to account for the effect of rod deaths on amplitude.</div></div><div><h3>Results</h3><div>The ROS generated in RP eyes were able to draw calcium ions, increasing the calcium influx. The voltage vs. time graph can be determined with the aid of a factor that represents the percentage of rods present. The loss of rod photoreceptors and calcium ions both result in a decrease in the amplitude of the <em>a</em>-wave at any time <em>t</em>.<!--> <!-->The theoretical results compare well with experimental results.</div></div><div><h3>Conclusion</h3><div>Due to oxidative stress, which causes rod photoreceptors to die and cause tunnel vision and peripheral visual field loss, the decline in scotopic <em>a</em>-wave amplitude in RP patients’ eyes is explained.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112306"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145427203","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}
Pub Date : 2026-01-21Epub Date: 2025-10-14DOI: 10.1016/j.jtbi.2025.112281
Saswati Biswas, Sudeshna Sinha
While ecosystems may experience sudden transitions to a degraded state under intensified exploitation, the impact of additional food provision in exploited patchy environments remains largely unexplored. This study investigates the trade-off between connectivity and resource allocation in mitigating tipping points that could lead to metacommunity-level population collapse. We first explore predator-prey dynamics within an isolated patch, investigating the effects of predator harvesting and additional food provision on population persistence. Our results reveal that, while additional food can rescue predators in scarcity, excessive provisioning may disrupt the trophic balance. Strategic harvesting helps mitigate this risk, but multistability across harvesting intensities complicates ecological management. Extending our analysis to a two-patch system with diffusive coupling, we find that a carefully calibrated food share ratio between patches is essential for long-term steady-state coexistence, with the required ratio modulated by coupling strength. However, beyond a critical dispersal threshold, stability can be maintained without strict adherence to a specific supply ratio. While dispersal aids in local predator rescue, higher flow can trigger a tipping point, resulting in catastrophic predator collapses across the metacommunity. Our findings reveal a potential rescue mechanism in which maintaining adequate food quality – ensuring uniformity across patches – is crucial to preventing abrupt population extinction, especially under strong connectivity. Overall, our study underscores the importance of integrating dispersal dynamics and the resource allocation mechanism in shaping ecosystem resilience, providing insight into strategies to mitigate population collapses in fragmented habitats.
{"title":"Strategic synergies: Dispersal and resource allocation in mitigating tipping cascades","authors":"Saswati Biswas, Sudeshna Sinha","doi":"10.1016/j.jtbi.2025.112281","DOIUrl":"10.1016/j.jtbi.2025.112281","url":null,"abstract":"<div><div>While ecosystems may experience sudden transitions to a degraded state under intensified exploitation, the impact of additional food provision in exploited patchy environments remains largely unexplored. This study investigates the trade-off between connectivity and resource allocation in mitigating tipping points that could lead to metacommunity-level population collapse. We first explore predator-prey dynamics within an isolated patch, investigating the effects of predator harvesting and additional food provision on population persistence. Our results reveal that, while additional food can rescue predators in scarcity, excessive provisioning may disrupt the trophic balance. Strategic harvesting helps mitigate this risk, but multistability across harvesting intensities complicates ecological management. Extending our analysis to a two-patch system with diffusive coupling, we find that a carefully calibrated food share ratio between patches is essential for long-term steady-state coexistence, with the required ratio modulated by coupling strength. However, beyond a critical dispersal threshold, stability can be maintained without strict adherence to a specific supply ratio. While dispersal aids in local predator rescue, higher flow can trigger a tipping point, resulting in catastrophic predator collapses across the metacommunity. Our findings reveal a potential rescue mechanism in which maintaining adequate food quality – ensuring uniformity across patches – is crucial to preventing abrupt population extinction, especially under strong connectivity. Overall, our study underscores the importance of integrating dispersal dynamics and the resource allocation mechanism in shaping ecosystem resilience, providing insight into strategies to mitigate population collapses in fragmented habitats.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112281"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145309863","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}
Pub Date : 2026-01-21Epub Date: 2025-10-14DOI: 10.1016/j.jtbi.2025.112282
Leonid Bunimovich , Athulya Ram
The paper continues the study of the phenomenon of local immunodeficiency in viral cross-immunoreactivity networks, with a focus on the roles and interactions between central and persistent viral variants. As usual, only the state of stable (i.e. observable) local immunodeficiency is analyzed. First, we show that a single central viral variant has an upper limit for the number of persistent viral variants that it can support. Our findings reveal that in viral cross-immunoreactivity networks, central viruses act essentially autonomously from each other. Namely, connections between central viruses change neither their qualitative roles nor the quantitative values of the strengths of their connections in the cross-immunoreactivity networks. In other words, each central virus does exactly the same actions, and has the same strengths with or without specific structural features such as central viruses. This indicates that local immunodeficiency can arise purely from the network structure. However, having more central viruses allows to keep the sizes of populations of persistent viruses at higher levels. Likewise, the strength of the immune response against any central virus remains at the same constant level regardless of how many persistent viruses this central virus supports (i.e. shields from the immune response of the host’s immune system). It is also shown that viruses strongly compete with each other in order to become persistent in the state of stable local immunodeficiency. We also present an (quite unexpected) example of a cross-immunoreactivity network with stable local immunodeficiency that only consists of persistent viral variants, which shows that persistent viruses may demonstrate a kind of self-consistency.
{"title":"Antigenic cooperation in viral populations: Maximal load on viruses and self-sufficiency of persistent viruses","authors":"Leonid Bunimovich , Athulya Ram","doi":"10.1016/j.jtbi.2025.112282","DOIUrl":"10.1016/j.jtbi.2025.112282","url":null,"abstract":"<div><div>The paper continues the study of the phenomenon of local immunodeficiency in viral cross-immunoreactivity networks, with a focus on the roles and interactions between central and persistent viral variants. As usual, only the state of stable (i.e. observable) local immunodeficiency is analyzed. First, we show that a single central viral variant has an upper limit for the number of persistent viral variants that it can support. Our findings reveal that in viral cross-immunoreactivity networks, central viruses act essentially autonomously from each other. Namely, connections between central viruses change neither their qualitative roles nor the quantitative values of the strengths of their connections in the cross-immunoreactivity networks. In other words, each central virus does exactly the same actions, and has the same strengths with or without specific structural features such as central viruses. This indicates that local immunodeficiency can arise purely from the network structure. However, having more central viruses allows to keep the sizes of populations of persistent viruses at higher levels. Likewise, the strength of the immune response against any central virus remains at the same constant level regardless of how many persistent viruses this central virus supports (i.e. shields from the immune response of the host’s immune system). It is also shown that viruses strongly compete with each other in order to become persistent in the state of stable local immunodeficiency. We also present an (quite unexpected) example of a cross-immunoreactivity network with stable local immunodeficiency that only consists of persistent viral variants, which shows that persistent viruses may demonstrate a kind of self-consistency.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112282"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145310056","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}
Pub Date : 2026-01-21Epub Date: 2025-10-31DOI: 10.1016/j.jtbi.2025.112302
Santana Mondal , Ludek Berec , Subhas Khajanchi
The predator-prey interaction is likely the most studied relationship in ecology. Predator consumption rate is modulated by more efficient hunting of prey by predators or, conversely, by more efficient hiding from predators by prey. Here, we study predator and prey adaptations in these traits using a common modeling framework, consisting of juvenile prey, adult prey, and an unstructured predator. What distinguishes the cases with predator and prey adaptations is the trade-offs. While more efficient predator hunting is assumed to be negatively correlated to its mortality rate (hunting trade-off), more efficient prey hiding is assumed to reduce prey reproduction (hiding trade-off). Natural selection benefits predators with high hunting efficiency if the hunting trade-off is strongly concave. For prey hiding efficiency, concave trade-offs give rise to a continuously stable strategy, whereas evolutionary branching is observed for convex trade-offs. In the latter case, the distance separating two branches on the evolutionary tree grows with increasing prey maturation rate; one branch eventually develops toward maximal hiding efficiency and minimal reproduction rate, while the other branch evolves towards minimal hiding efficiency and maximal reproduction rate. The concave trade-off thus serves a dual function in prey and predator evolution. For prey-hiding evolution, it increases prey protection. For predator-hunting evolution, it decreases hunting efficiency and thus promotes prey protection, too. Thus, the trade-off form proves to be a crucial determinant of evolutionary outcomes.
{"title":"Evolution of hunting efficiency of predator and hiding efficiency of prey in a stage-structured prey-predator model","authors":"Santana Mondal , Ludek Berec , Subhas Khajanchi","doi":"10.1016/j.jtbi.2025.112302","DOIUrl":"10.1016/j.jtbi.2025.112302","url":null,"abstract":"<div><div>The predator-prey interaction is likely the most studied relationship in ecology. Predator consumption rate is modulated by more efficient hunting of prey by predators or, conversely, by more efficient hiding from predators by prey. Here, we study predator and prey adaptations in these traits using a common modeling framework, consisting of juvenile prey, adult prey, and an unstructured predator. What distinguishes the cases with predator and prey adaptations is the trade-offs. While more efficient predator hunting is assumed to be negatively correlated to its mortality rate (hunting trade-off), more efficient prey hiding is assumed to reduce prey reproduction (hiding trade-off). Natural selection benefits predators with high hunting efficiency if the hunting trade-off is strongly concave. For prey hiding efficiency, concave trade-offs give rise to a continuously stable strategy, whereas evolutionary branching is observed for convex trade-offs. In the latter case, the distance separating two branches on the evolutionary tree grows with increasing prey maturation rate; one branch eventually develops toward maximal hiding efficiency and minimal reproduction rate, while the other branch evolves towards minimal hiding efficiency and maximal reproduction rate. The concave trade-off thus serves a dual function in prey and predator evolution. For prey-hiding evolution, it increases prey protection. For predator-hunting evolution, it decreases hunting efficiency and thus promotes prey protection, too. Thus, the trade-off form proves to be a crucial determinant of evolutionary outcomes.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112302"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145433050","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}
Pub Date : 2026-01-21Epub Date: 2025-10-15DOI: 10.1016/j.jtbi.2025.112287
Alessandra Cambi , Diane S. Lidke , Mariya Ptashnyk , Willemijn Smit , Stefanie Sonner
G protein-coupled receptors EP2 and EP4 are both activated by the lipid messenger Prostaglandin E2 (PGE2) and induce the intracellular production of cyclic AMP (cAMP), ultimately affecting gene expression. Changes in cellular responses to PGE2 can have important consequences on immunity and disease, yet a detailed understanding of the EP2-EP4 signaling network is lacking. EP2 and EP4 are often co-expressed in cells but their specific contribution to cAMP production is poorly understood. Experimental data have shown that cAMP levels differ depending on whether PGE2 triggers EP2 or EP4, or both. To better understand the underlying mechanisms and predict cellular responses to PGE2, we developed mathematical models for EP2 and EP4 cAMP signaling, including receptor crosstalk. The mathematical models qualitatively reproduce the experimentally observed cAMP levels and provide mechanistic insight into both the differences and commonalities in EP2/EP4 signaling. We found that ligand binding dynamics play a crucial role for both single-receptor signaling and inter-receptor crosstalk. Inhibition of PGE2 signaling via receptor antagonists is gaining increasing attention in tumor immunology. These mathematical models could therefore contribute to the design of more effective anti-tumor therapies targeting EP2 and EP4.
{"title":"Mathematical models for the EP2 and EP4 signaling pathways and their crosstalk","authors":"Alessandra Cambi , Diane S. Lidke , Mariya Ptashnyk , Willemijn Smit , Stefanie Sonner","doi":"10.1016/j.jtbi.2025.112287","DOIUrl":"10.1016/j.jtbi.2025.112287","url":null,"abstract":"<div><div>G protein-coupled receptors EP2 and EP4 are both activated by the lipid messenger Prostaglandin E2 (PGE2) and induce the intracellular production of cyclic AMP (cAMP), ultimately affecting gene expression. Changes in cellular responses to PGE2 can have important consequences on immunity and disease, yet a detailed understanding of the EP2-EP4 signaling network is lacking. EP2 and EP4 are often co-expressed in cells but their specific contribution to cAMP production is poorly understood. Experimental data have shown that cAMP levels differ depending on whether PGE2 triggers EP2 or EP4, or both. To better understand the underlying mechanisms and predict cellular responses to PGE2, we developed mathematical models for EP2 and EP4 cAMP signaling, including receptor crosstalk. The mathematical models qualitatively reproduce the experimentally observed cAMP levels and provide mechanistic insight into both the differences and commonalities in EP2/EP4 signaling. We found that ligand binding dynamics play a crucial role for both single-receptor signaling and inter-receptor crosstalk. Inhibition of PGE2 signaling via receptor antagonists is gaining increasing attention in tumor immunology. These mathematical models could therefore contribute to the design of more effective anti-tumor therapies targeting EP2 and EP4.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"617 ","pages":"Article 112287"},"PeriodicalIF":2.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145314168","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}