{"title":"Decoding cancer dynamics: The role of miRNA in E2F/Myc networks and bifurcation points","authors":"Lanjiang Luo, Jianwei Shen","doi":"10.1016/j.jtbi.2025.112219","DOIUrl":"10.1016/j.jtbi.2025.112219","url":null,"abstract":"<div><div>Cancer development is driven by critical state transitions in gene regulatory networks, particularly the miRNA-E2F/Myc pathway, yet the nonlinear dynamics underlying this process remain relatively scarce. Here, we construct a dynamic model incorporating parameters, time delays, and diffusion to elucidate how bifurcations and spatiotemporal oscillations drive carcinogenesis. The model reveals that regulatory variations can induce saddle-node, Hopf, Bogdanov-Takens, and Cusp bifurcations, with key events-such as extracellular matrix protein levels exceeding a threshold-triggering irreversible oncogenic switches. Time delays and diffusion further generate gene oscillations and spatial heterogeneity through Hopf and Turing-Hopf bifurcations; these dynamic signatures can serve as early warning indicators, and our multiscale modeling enables quantitative prediction of critical transitions, providing a theoretical basis for optimizing targeted therapy timing.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"613 ","pages":"Article 112219"},"PeriodicalIF":2.0,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144700348","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}
Diepreye Ayabina , Hasan Sevil , Adam Kleczkowski , M.Gabriela-M. Gomes
{"title":"A note on metapopulation models","authors":"Diepreye Ayabina , Hasan Sevil , Adam Kleczkowski , M.Gabriela-M. Gomes","doi":"10.1016/j.jtbi.2025.112213","DOIUrl":"10.1016/j.jtbi.2025.112213","url":null,"abstract":"<div><div>Metapopulation models are commonly used in ecology, evolution, and epidemiology. These models usually entail homogeneity assumptions within patches and study networks of migration between patches to generate insights into conservation of species, differentiation of populations, and persistence of infectious diseases. Here, focusing on infectious disease epidemiology, we take a complementary approach and study the effects of individual variation within patches while neglecting any form of disease transmission between patches. Consistently with previous work on single populations, we show how metapopulation models that neglect in-patch heterogeneity also underestimate basic reproduction numbers (<span><math><msub><mi>R</mi><mn>0</mn></msub></math></span>) and the effort required to control or eliminate infectious diseases by uniform interventions. We then go beyond this confirmatory result and introduce a scheme to infer distributions of individual susceptibility or exposure to infection based on suitable stratifications of a population into patches. We apply the resulting metapopulation models to a simple case study of the COVID-19 pandemic.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"613 ","pages":"Article 112213"},"PeriodicalIF":1.9,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144651265","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":"The role of selective extinction in promoting cooperation in a many demes model","authors":"Dhaker Kroumi","doi":"10.1016/j.jtbi.2025.112212","DOIUrl":"10.1016/j.jtbi.2025.112212","url":null,"abstract":"<div><div>This study investigates the evolution of cooperation in a deme-structured population where individuals interact via the Prisoner’s Dilemma. Reproduction is payoff-dependent, and demes contribute to a global migrant pool proportionally to their fecundity. Deme extinction occurs with a baseline probability plus an additional probability term that is proportional to the fraction of defectors in the deme, modulated by a parameter <span><math><mi>α</mi></math></span>. Extinct demes are then repopulated from the migrant pool. A two-timescale analysis, valid in the limit of a large number of demes, reveals that the system’s evolutionary dynamics are well-approximated by a continuous-time diffusion process. Using this diffusion approximation, we show that uniform extinction (<span><math><mrow><mi>α</mi><mo>=</mo><mn>0</mn></mrow></math></span>) suppresses cooperation. However, extinction targeting defector-dominated demes can significantly promote cooperation, establishing critical thresholds of <span><math><mi>α</mi></math></span> that determine when cooperation dominates. The study further explores how varying extinction rates and reproductive dynamics impact cooperation’s persistence and evolutionary success in structured populations.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"613 ","pages":"Article 112212"},"PeriodicalIF":1.9,"publicationDate":"2025-07-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144638735","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":"Modeling alternative strategies of male emergence timing","authors":"Hidaka Kubo , Shinji Nakaoka , Ryo Yamaguchi","doi":"10.1016/j.jtbi.2025.112214","DOIUrl":"10.1016/j.jtbi.2025.112214","url":null,"abstract":"<div><div>In many butterfly species, males emerge earlier than females as part of a strategy to maximize male reproductive success. Although behavioral ecological studies using mathematical models have been conducted to explain this phenomenon, certain emergence patterns remain unexplained. In the butterfly species <em>Fabriciana nerippe</em>, some males emerge at the same time as females, in addition to males that emerge earlier than the females. However, it is unclear what emergence patterns occur in populations with male dimorphism, as observed in this species. In this study, we showed the existence of male body size dimorphism in <em>Fabriciana nerippe</em> by conducting a comparative analysis of forewing lengths between males and females. In addition, we developed a comprehensive mathematical model to investigate emergence patterns in the presence of dimorphic males. By introducing a trade-off between large size and early emergence, the model considered a scenario where small early-emerging and large late-emerging males could coexist. Numerical analysis demonstrated the emergence patterns of these two male types with a switch in emergence time. Furthermore, the higher the death rate before emergence, the earlier the emergence switch. These findings suggested that the timing of the switch depends on the death rate and is influenced by environmental factors. This work contributes to ecological and theoretical studies on two types of emergence timing in life-history strategies across a broader range of species.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112214"},"PeriodicalIF":1.9,"publicationDate":"2025-07-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144621183","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":"Success, failure and collapse of an ecological system due to the introduction of a generalist super-predator","authors":"Ayanava Basak, Nandadulal Bairagi","doi":"10.1016/j.jtbi.2025.112198","DOIUrl":"10.1016/j.jtbi.2025.112198","url":null,"abstract":"<div><div>Species reintroduction is an effective strategy to save endangered species. New species are also introduced for economic purposes or to control the density of some other species. However, not all species introduction programs are successful. In this article, using mathematical models and the theory of dynamical systems, we explain a species introduction program’s success, failure and collapse due to introducing a generalist super-predator in a predator-prey system. Our study reveals that the persistence of the introduced species with the others mainly depends on two factors: the population density of the native species at the time of the third species’ introduction and the population density of the introduced species. Further, we show that the introduced species can have significant negative impacts on native species, often leading to their replacement by themselves. In extreme cases, this disruption can result in the extinction of native prey and predators along with the introduced species. This novel study may find applications in conservation biology and species control mechanisms.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112198"},"PeriodicalIF":1.9,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144546226","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":"Response of social norms to individual differences in error-proneness","authors":"Quang Anh Le , Seung Ki Baek","doi":"10.1016/j.jtbi.2025.112199","DOIUrl":"10.1016/j.jtbi.2025.112199","url":null,"abstract":"<div><div>Indirect reciprocity explains the evolution of cooperation by considering how our cooperative behavior toward someone is reciprocated by someone else who has observed us. A cohesive society has a shared norm that prescribes how to assess observed behavior as well as how to behave toward others based on the assessments, and the eight social norms that are evolutionarily stable against the invasion of mutants with different behavioral rules are referred to as the leading eight, whose member norms are called L1 to L8, respectively. Among the leading eight, L8 (also known as ‘Judging’) has been deemed mostly irrelevant due to its poor performance in maintaining cooperation when each person may have a different opinion about someone instead of forming a public consensus. In this work, we propose that L8 can nevertheless be best protected from assessment errors among the leading eight if we take into account the intrinsic heterogeneity of error proneness among individuals because the more prone to assessment errors, the more disadvantaged under the action of L8. This finding suggests that individual differences should be explicitly taken into account as quenched randomness to obtain a thorough understanding of a social norm working in a heterogeneous environment.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112199"},"PeriodicalIF":1.9,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144546225","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}
Maxime Estavoyer , Thomas Lepoutre , Marie Manceau
{"title":"Mathematical modeling of the feather follicle morphogenetic wave in birds","authors":"Maxime Estavoyer , Thomas Lepoutre , Marie Manceau","doi":"10.1016/j.jtbi.2025.112200","DOIUrl":"10.1016/j.jtbi.2025.112200","url":null,"abstract":"<div><div>During the development of the avian skin, feather follicles are produced in a medio-lateral morphogenetic wave that results in their spatial arrangement in typical patterns. This wave involves the timely acquisition of pattern-forming competence followed by a row-by-row production of feather follicles. While several mathematical models combining self-organizing systems accurately reproduced dynamics of feather follicle pattern formation, the events that control timely parameters of wave propagation remain poorly understood. Here, we built on previous modeling work to theoretically calculate the speed at which tissue competence progresses. Using a weakly non-linear analysis, we calculated the speed at which follicles emerge once competence is attained. We produced numerical simulations of our model to predict the respective influences of competence acquisition and follicle emergence on each other and on wave propagation. Our results show that the theoretical speed of follicle emergence is limited by competence acquisition, but that, in turn, competence acquisition is not constrained by follicle emergence. This modeling work provides an approximation of the timely parameters of the morphogenetic wave, and sheds light on the interplay between competence and patterning events in the developing skin.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112200"},"PeriodicalIF":1.9,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144531172","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}
Helen Kettle , Duncan J. Coston , Sacha White , Tom W. Pope , Joe M. Roberts , David Ewing
{"title":"A process-based life cycle model of the cabbage stem flea beetle on winter oilseed rape: Effects of temperature","authors":"Helen Kettle , Duncan J. Coston , Sacha White , Tom W. Pope , Joe M. Roberts , David Ewing","doi":"10.1016/j.jtbi.2025.112185","DOIUrl":"10.1016/j.jtbi.2025.112185","url":null,"abstract":"<div><div>The cabbage stem flea beetle (CSFB) is an economically important pest of oilseed rape crops that has increased in importance over recent years. Here we present the first attempt at a process-based model of this species using delay-differential equations (implemented in the R package stagePop) to simulate the abundance of each stage of the life cycle. We split the life cycle into 6 distinct stages: eggs, larvae, pupae-diapause (pupae, pre-aestivation adults and aestivating adults), pre-oviposition adults, reproductive adults and non-reproductive adults (past reproductive age), with temperature-dependent stage durations and death rates. We use adult trap data to inform on the initial conditions for the model and drive the model with a generic air temperature time series. To validate our model we compare the results with expected timings based on the scientific literature and field data on larval abundance. We use the model to investigate the effects of increasing the yearly average air temperature, and the effects of unusually warm seasons on CSFB population seasonal abundance patterns. Due to a lack of data there are key processes (e.g. temperature dependence of larval development) and parameter values that are uncertain. However, we believe the model, and its future iterations as more data are obtained, has the potential to be used as part of a decision support system for farmers and advisors, enabling more informed decisions on drilling and IPM strategies based on the within season weather conditions.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112185"},"PeriodicalIF":1.9,"publicationDate":"2025-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144523563","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":"Accurate stochastic simulation algorithm for multiscale models of infectious diseases","authors":"Yuan Yin , Jennifer A. Flegg , Mark B. Flegg","doi":"10.1016/j.jtbi.2025.112194","DOIUrl":"10.1016/j.jtbi.2025.112194","url":null,"abstract":"<div><div>In the infectious disease literature, significant effort has been devoted to studying dynamics at a single scale. For example, compartmental models describing population-level dynamics are often formulated using differential equations. In cases where small numbers or noise play a crucial role, these differential equations are replaced with memoryless Markovian models, where discrete individuals can be members of a compartment and transition stochastically. Classic stochastic simulation algorithms, such as the next reaction method, can be employed to solve these Markovian models exactly. The intricate coupling between models at different scales underscores the importance of multiscale modelling in infectious diseases. However, several computational challenges arise when the multiscale model becomes non-Markovian. In this paper, we address these challenges by developing a novel exact stochastic simulation algorithm. We apply it to a showcase multiscale system where all individuals share the same deterministic within-host model while the population-level dynamics are governed by a stochastic formulation. We demonstrate that as long as the within-host information is harvested at a reasonable resolution, the novel algorithm will always be accurate. Furthermore, our implementation is still efficient even at finer resolutions. Beyond infectious disease modelling, the algorithm is widely applicable to other multiscale systems, providing a versatile, accurate, and computationally efficient framework.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112194"},"PeriodicalIF":1.9,"publicationDate":"2025-06-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144487150","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}
Kamryn Spinelli , Chaozhen Wei , Luis Vidali , Min Wu
{"title":"Inferring exocytosis profiles from cell shapes using a dual-configuration model of walled cell tip growth","authors":"Kamryn Spinelli , Chaozhen Wei , Luis Vidali , Min Wu","doi":"10.1016/j.jtbi.2025.112196","DOIUrl":"10.1016/j.jtbi.2025.112196","url":null,"abstract":"<div><div>Tip growth in filamentous cells, such as root hairs, moss protonemata, and fungal hyphae, depends on coordinated cell wall extension driven by turgor pressure, wall mechanics, and exocytosis. We introduce a dual-configuration model that incorporates both turgid and unturgid states to describe cell wall growth as the combined effect of elastic deformation and irreversible extension. This framework infers exocytosis profiles directly from cell morphology and elastic stretches, formulated as an initial value problem based on the self-similarity condition. Applying the model to <em>Medicago truncatula</em> root hairs, moss <em>Physcomitrium patens</em> protonemata, and hyphoid-like shapes, we find that exocytosis peaks at the tip in tapered cells but shifts to an annular region away from the apex in flatter-tip cells beyond a threshold. The model generalizes previous fluid models and provides a mechanistic link between exocytosis distribution and cell shape, explaining observed variations in tip-growing cells across species.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"612 ","pages":"Article 112196"},"PeriodicalIF":1.9,"publicationDate":"2025-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144477865","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}