Alexandre Benatti , Henrique Ferraz De Arruda , Luciano Da Fontoura Costa
{"title":"Interrelating neuronal morphology by coincidence similarity networks","authors":"Alexandre Benatti , Henrique Ferraz De Arruda , Luciano Da Fontoura Costa","doi":"10.1016/j.jtbi.2025.112104","DOIUrl":"10.1016/j.jtbi.2025.112104","url":null,"abstract":"<div><div>The study of neuronal morphology presents potential not only for identifying possible relationship with neuronal dynamics, but also as a means to characterize and classify types of neuronal cells and compare them among species, organs, and conditions. In the present work, we approach this problem by using the concept of coincidence similarity index, considering a methodology for mapping datasets into similarity networks. The adopted similarity presents some specific interesting properties, including more strict comparisons. A set of 20 morphological features has been considered, and coincidence similarity networks estimated respectively to 735 considered neuronal cells from 8 groups of Drosophila melanogaster.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"606 ","pages":"Article 112104"},"PeriodicalIF":1.9,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143733228","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}
Adrita Ghosh , Parthasakha Das , Tanujit Chakraborty , Pritha Das , Dibakar Ghosh
{"title":"Developing cholera outbreak forecasting through qualitative dynamics: Insights into Malawi case study","authors":"Adrita Ghosh , Parthasakha Das , Tanujit Chakraborty , Pritha Das , Dibakar Ghosh","doi":"10.1016/j.jtbi.2025.112097","DOIUrl":"10.1016/j.jtbi.2025.112097","url":null,"abstract":"<div><div>Cholera, an acute diarrheal disease, is a serious concern in developing and underdeveloped areas. A qualitative understanding of cholera epidemics aims to foresee transmission patterns based on reported data and mechanistic models. The mechanistic model is a crucial tool for capturing the dynamics of disease transmission and population spread. However, using real-time cholera cases is essential for forecasting the transmission trend. This prospective study seeks to furnish insights into transmission trends through qualitative dynamics followed by machine learning-based forecasting. The Monte Carlo Markov Chain approach is employed to calibrate the proposed mechanistic model. We identify critical parameters that illustrate the disease’s dynamics using partial rank correlation coefficient-based sensitivity analysis. The basic reproduction number as a crucial threshold measures asymptotic dynamics. Furthermore, forward bifurcation directs the stability of the infection state, and Hopf bifurcation suggests that trends in transmission may become unpredictable as societal disinfection rates rise. Further, we develop epidemic-informed machine learning models by incorporating mechanistic cholera dynamics into autoregressive integrated moving averages and autoregressive neural networks. We forecast short-term future cholera cases in Malawi by implementing the proposed epidemic-informed machine learning models to support this. We assert that integrating temporal dynamics into the machine learning models can enhance the capabilities of cholera forecasting models. The execution of this mechanism can significantly influence future trends in cholera transmission. This evolving approach can also be beneficial for policymakers to interpret and respond to potential disease systems. Moreover, our methodology is replicable and adaptable, encouraging future research on disease dynamics.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112097"},"PeriodicalIF":1.9,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143694664","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":"An optimal network that promotes the spread of an advantageous variant in an SIR epidemic","authors":"Samuel Lopez , Natalia L. Komarova","doi":"10.1016/j.jtbi.2025.112095","DOIUrl":"10.1016/j.jtbi.2025.112095","url":null,"abstract":"<div><div>In the course of epidemics, the pathogen may mutate to acquire a higher fitness. At the same time, such a mutant is automatically in an unfavorable position because the resident virus has a head start in accessing the pool of susceptible individuals. We considered a class of tunable small-world networks, where a parameter, <span><math><mi>p</mi></math></span> (the rewiring probability), characterizes the prevalence of non-local connections, and we asked, whether the underlying network can influence the fate of a mutant virus. Under an SIR model, we considered two measures of mutant success: the expected height of the peak of mutant infected individuals, and the total number of recovered from mutant individuals at the end of the epidemic. Using these measures, we have found the existence of an optimal (for an advantageous mutant virus) rewiring probability that promotes a larger infected maximum and a larger total recovered population corresponding to the advantageous pathogen strain. This optimal rewiring probability decreases as mean degree and the infectivity of the wild type are increased, and it increases with the mutant advantage. The non-monotonic behavior of the advantageous mutant as a function of rewiring probability may shed light into some of the complex patterns in the size of mutant peaks experienced by different countries during the COVID-19 pandemic.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112095"},"PeriodicalIF":1.9,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143665383","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 model of tubulin removal and exchange caused by kinesin motor walking on microtubule lattices","authors":"Ping Xie","doi":"10.1016/j.jtbi.2025.112088","DOIUrl":"10.1016/j.jtbi.2025.112088","url":null,"abstract":"<div><div>The kinesin motor walking on microtubule lattices can cause disassembly of GDP-tubulins, generating defects, and repair the defects by incorporating GTP-tubulins. To explore the underlying mechanism, a model is proposed here. On the basis of the model, the dynamics of the defect generation, defect repair and tubulin exchange induced by the kinesin motor is studied theoretically. The theoretical results explain well the available experimental data. Moreover, predicted results are provided.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112088"},"PeriodicalIF":1.9,"publicationDate":"2025-03-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644146","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 spatial aggregation of phytophagous insects driven by the evolution of preference for plant chemicals","authors":"Haruna Ohsaki , Akira Yamawo , Yuuya Tachiki","doi":"10.1016/j.jtbi.2025.112094","DOIUrl":"10.1016/j.jtbi.2025.112094","url":null,"abstract":"<div><div>Ecologists have shown considerable interest in the spatial patterns of organism distribution and the processes responsible for their formation and maintenance. The phytophagous insects typically use chemicals in plants as host-finding cues. Because nonvolatile chemicals remain near the source, the spatial structure of plant community determines the local distribution of insects. In addition, the plant chemical accumulation due to plant–plant interaction also influences the distribution of insects. In <em>Rumex obtusifolius</em>, for example, the production of phenolics is mediated by conspecific interaction. <em>Rumex</em> plants with high phenolic concentrations are preferred by the leaf beetle <em>Gastrophysa atrocyanea</em>, resulting in its spatial aggregation. Although this preference of beetles for nonvolatile chemicals should be beneficial in finding host plants, there is also a cost in terms of intraspecific competition among the beetles due to aggregation on certain chemical-rich hosts. To investigate the evolutionary significance of preference for nonvolatile chemicals and the ecological consequence of spatial distribution in leaf beetles, we constructed a mathematical model for the joint evolution of two preferences for plant size and chemical condition. In the model, beetles choose a resource based on the size and chemical concentrations of plants and are exposed to resource competition. Host plants accumulate the chemicals when they interact with neighboring conspecifics, and hence the level of chemical accumulation varies depending on the species composition and spatial distribution of the plant community. As a result, beetles became more sensitive to chemicals when the host species was rare and sparsely distributed in the community. The evolution of high chemical preference caused the aggregation of beetles and hence population size declined. We proposed a potential mechanism that underlies aggregated distribution in phytophagous insects, driven by the evolution of chemical preferences in response to plant community structure.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112094"},"PeriodicalIF":1.9,"publicationDate":"2025-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143639853","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":"Absorbing Markov chain model of PrEP drug adherence to estimate adherence decay rate and probability distribution in clinical trials","authors":"Renee Dale , Hongyu He , Yingqing Chen","doi":"10.1016/j.jtbi.2025.112086","DOIUrl":"10.1016/j.jtbi.2025.112086","url":null,"abstract":"<div><div>Pre-exposure prophylaxis (PrEP) is increasingly used to prevent the transmission of H.I.V. in at-risk populations. However, PrEP users may discontinue use of the medicine due to side effects, lower perceived risk, or other reasons. The usage metrics of 594 individuals was tracked over 350 days using the Wisepill electronic monitoring system. We model the PrEP drug adherence level using an absorbing Markov chain with a unique absorbing state. The transition matrix <span><math><mi>T</mi></math></span> obtained from the Wisepill data will have a trivial eigenvector (eigendistribution) associated with the first (i.e., largest) eigenvalue 1. The 2nd eigenvalue(s) then become important in determining the asymptotic behavior of the Markov chain, dictating how fast the Markov chain <em>decays</em> to the absorbing state. Under a fairly general assumption, we prove that the second positive eigenvalue is unique and the corresponding eigenvector will have nonnegative entries with exceptions at absorbing states. In addition, we define the asymptotic half life of the absorbing Markov chain directly from the 2nd eigenvalue. We then determine the 2nd eigenvalue of <span><math><mi>T</mi></math></span> and the asymptotic half life of the Markov chain, which turns out to be very close to the real half life of the Markov chain. Finally, we interpret the 2nd eigenvector as the relative probability distribution of <span><math><msub><mi>X</mi><mi>∞</mi></msub></math></span> with respect to the decay rate of the 2nd eigenvalue. By applying these methods to the Wisepill data, we estimate the half-life of population adherence to be 46 weeks. The bi-weekly decay rate observed in these data from 90 to 100 % adherence is 3 %. This work produces an estimate at which adherence falls over time, given no external intervention is applied. These results suggest an eigenvector-based approach to estimate adherence trends, as well as the timing of interventions to improve adherence.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"604 ","pages":"Article 112086"},"PeriodicalIF":1.9,"publicationDate":"2025-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143611685","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":"Branching architecture affects genetic diversity within an individual tree","authors":"Sou Tomimoto , Yoh Iwasa , Akiko Satake","doi":"10.1016/j.jtbi.2025.112093","DOIUrl":"10.1016/j.jtbi.2025.112093","url":null,"abstract":"<div><div>While a tree grows over many years, somatic mutations accumulate and form genetic variation among branches within an individual. Trees can transmit such mutations to subsequent generations, potentially enhancing the genetic diversity of the population. We study a mathematical model to understand the relationship between within-individual genetic variation and branching architecture. We generate branching architecture by repeatedly adding two new branches (main and lateral daughter branches) to each terminal branch (mother branch). The architecture is characterized by two key parameters: main-lateral ratio (ML) and daughter-mother ratio (DM). During branch elongation, somatic mutations accumulate in the stem cells of a shoot apical meristem (SAM) at the tip of each branch. In branching, all the stem cells are passed on from the mother to the main daughter branch, but only one stem cell is chosen for the lateral daughter branch. We evaluate genetic variation by <span><math><mover><mrow><mi>Z</mi></mrow><mrow><mo>¯</mo></mrow></mover></math></span>, the mean genetic differences between all pairs of branches of a tree, and examine how <span><math><mover><mrow><mi>Z</mi></mrow><mrow><mo>¯</mo></mrow></mover></math></span> varies with DM and ML while keeping the total branch length constant. As a result, (1) <span><math><mover><mrow><mi>Z</mi></mrow><mrow><mo>¯</mo></mrow></mover></math></span> increases monotonically with ML; (2) <span><math><mover><mrow><mi>Z</mi></mrow><mrow><mo>¯</mo></mrow></mover></math></span> attains the maximum for an intermediate DM, when stem cells in a SAM are genetically homogeneous; (3) <span><math><mover><mrow><mi>Z</mi></mrow><mrow><mo>¯</mo></mrow></mover></math></span> decreases monotonically with DM when stem cells are heterogeneous. The effect of branching architecture varies significantly depending on the genetic heterogeneity within a SAM, which results from the behavior of stem cells during growth. Our study sheds light on the overlooked role of branching architecture in storing genetic diversity.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112093"},"PeriodicalIF":1.9,"publicationDate":"2025-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143634975","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":"Conditions for the establishment of creole languages from an evolutionary game theoretic perspective","authors":"Raiki Nakano, Hisashi Ohtsuki","doi":"10.1016/j.jtbi.2025.112090","DOIUrl":"10.1016/j.jtbi.2025.112090","url":null,"abstract":"<div><div>Language contact refers to a situation in which speakers of different languages meet and interact. There are three consequences of language contact. The first one is language shift, where one language is replaced by another. The second one is language maintenance, where languages coexist without any language change. The question of how language shift takes place and under what conditions coexistence of languages is possible has drawn wide attention of linguists as well as of theoretical researchers in other disciplines. Most of the previous mathematical studies of language dynamics have shown that there is majority advantage in the competition and that the population eventually becomes dominated by a single language that was spoken by a majority when the competition started. However, it is known that the contact of different languages could lead to the third consequence: the emergence of a new language, such as creole languages. Our study aims to model the spread of creole languages by using evolutionary game theory and determine the factors that drive the spread of creoles in the population. We consider an evolutionary game with individuals who speak either one of two ancestral languages or the creole that resulted from language contact between those two ancestral languages. Due to the formation process of creole, the creole speakers can communicate with speakers of other languages to some, but not full, extent. We assume that their frequencies change according to their payoffs that reflect how well they can communicate with others. Under this basic assumption, we varied (i) the population size and (ii) the difficulty of language change, to see their effects on the spread of the creole. Results show that the creole is more likely to spread if the population size is small enough for random drift to work well and/or if two ancestral languages are linguistically distant. We also find that finiteness of the population can either favor or disfavor the spread of creole depending on detailed conditions. Our results suggest that the contact between two distant languages and random drift drive the spread of creoles. It is known that many creole languages appeared in plantations where European owners brought labor forces from Africa as slaves. Our results match the homestead phase of these plantations, where European languages are distant from African ones, and where the population size was small.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112090"},"PeriodicalIF":1.9,"publicationDate":"2025-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143631065","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 short comment on the individual response to ionizing radiation","authors":"Krzysztof Wojciech Fornalski","doi":"10.1016/j.jtbi.2025.112092","DOIUrl":"10.1016/j.jtbi.2025.112092","url":null,"abstract":"","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"604 ","pages":"Article 112092"},"PeriodicalIF":1.9,"publicationDate":"2025-03-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143594244","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":"Emergence of food webs with a multi-trophic hierarchical structure driven by nonlinear trait-matching","authors":"Christophe Laplanche, Benjamin Pey, Robin Aguilée","doi":"10.1016/j.jtbi.2025.112091","DOIUrl":"10.1016/j.jtbi.2025.112091","url":null,"abstract":"<div><div>Food webs are a central subject in community ecology, because consumption supports the flow of matter through the system, which is at the base of many of its functions. Identifying the mechanisms that are at the origin of food web structure is useful, e.g., for restoration purposes. We investigated the extent to which trait-matching, which contributes to defining the strength of trophic interactions, can cause the emergence of food webs with a non-trivial, multi-trophic, hierarchical structure. We compared for that purpose the structural properties of food webs simulated by four food web model variants, depending whether trait-matching was linear or nonlinear and whether population dynamics and evolution were accounted for (dynamical model) or not (static model). Nonlinear trait-matching can restrict interactions in phenotypic space so as to obtain localized interactions (i.e., each species interact with a small subset of species), which is a key element for food web formation. In the static case, nonlinear trait-matching allowed for the emergence of food webs, at a relatively low connectance as with random graphs. In the dynamical case, nonlinear trait-matching combined with population dynamics and evolution allowed for the formation of groups of phenotypically close species, resulting in food webs with a multi-trophic, hierarchical structure.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"605 ","pages":"Article 112091"},"PeriodicalIF":1.9,"publicationDate":"2025-03-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143588243","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}