Journal of Mathematical Biology

To maintain renewing epithelial tissues in a healthy, homeostatic state, cell divisions and differentiation need to be tightly regulated. Mechanisms of homeostatic regulation often rely on crowding feedback control: cells are able to sense the cell density in their environment, via various molecular and mechanosensing pathways, and respond by adjusting division, differentiation, and cell state transitions appropriately. Here, we determine, via a mathematically rigorous framework, which general conditions for the crowding feedback regulation (i) must be minimally met, and (ii) are sufficient, to allow the maintenance of homeostasis in renewing tissues...

Emerging and re-emerging pathogens are latent threats in our society with the risk of killing millions of people worldwide, without forgetting the severe economic and educational backlogs. From COVID-19, we learned that self isolation and quarantine restrictions (confinement) were the main way of protection till availability of vaccines. However, abrupt lifting of social confinement would result in new waves of new infection cases and high death tolls. Here, inspired by how an extracellular solution can make water move into or out of a cell through osmosis, we define confinement tonicity...

COVID-19 is a respiratory disease triggered by an RNA virus inclined to mutations. Since December 2020, variants of COVID-19 (especially Delta and Omicron) continuously appeared with different characteristics that influenced death and transmissibility emerged around the world. To address the novel dynamics of the disease, we propose and analyze a dynamical model of two strains, namely native and mutant, transmission dynamics with mutation and imperfect vaccination. It is also assumed that the recuperated individuals from the native strain can be infected with mutant strain through the direct contact with individual or contaminated surfaces or aerosols...

We consider stochastic dynamics of a population which starts from a small colony on a habitat with large but limited carrying capacity. A common heuristics suggests that such population grows initially as a Galton-Watson branching process and then its size follows an almost deterministic path until reaching its maximum, sustainable by the habitat. In this paper we put forward an alternative and, in fact, more accurate approximation which suggests that the population size behaves as a special nonlinear transformation of the Galton-Watson process from the very beginning...

We study a viral infection model incorporating both cell-to-cell infection and immune chemokines. Based on experimental results in the literature, we make a standing assumption that the cytotoxic T lymphocytes (CTL) will move toward the location with more infected cells, while the diffusion rate of CTL is a decreasing function of the density of infected cells. We first establish the global existence and ultimate boundedness of the solution via a priori energy estimates. We then define the basic reproduction number of viral infection <mml:math xmlns:mml="https://www...

In the Antarctic, the whale population had been reduced dramatically due to the unregulated whaling. It was expected that Antarctic krill, the main prey of whales, would grow significantly as a consequence and exploratory krill fishing was practiced in some areas. However, it was found that there has been a substantial decline in abundance of krill since the end of whaling, which is the phenomenon of krill paradox. In this paper, to study the krill-whale interaction we revisit a harvested predator-prey model with Holling I functional response...

Clinical and pre-clinical data suggest that treating some tumors at a mild, patient-specific dose might delay resistance to treatment and increase survival time. A recent mathematical model with sensitive and resistant tumor cells identified conditions under which a treatment aiming at tumor containment rather than eradication is indeed optimal. This model however neglected mutations from sensitive to resistant cells, and assumed that the growth-rate of sensitive cells is non-increasing in the size of the resistant population...

Phylogenetic diversity is a popular measure for quantifying the biodiversity of a collection Y of species, while phylogenetic diversity indices provide a way to apportion phylogenetic diversity to individual species. Typically, for some specific diversity index, the phylogenetic diversity of Y is not equal to the sum of the diversity indices of the species in Y. In this paper, we investigate the extent of this difference for two commonly-used indices: Fair Proportion and Equal Splits. In particular, we determine the maximum value of this difference under various instances including when the associated rooted phylogenetic tree is allowed to vary across all rooted phylogenetic trees with the same leaf set and whose edge lengths are constrained by either their total sum or their maximum value...

The presence or absence of synaptic plasticity can dramatically influence the collective behavior of populations of coupled neurons. In this work, we consider spike-timing dependent plasticity (STDP) and its resulting influence on phase cohesion in computational models of heterogeneous populations of conductance-based neurons. STDP allows for the influence of individual synapses to change over time, strengthening or weakening depending on the relative timing of the relevant action potentials. Using phase reduction techniques, we derive an upper bound on the critical coupling strength required to retain phase cohesion for a network of synaptically coupled, heterogeneous neurons with STDP...

In this paper, we study a time-delayed nonlocal reaction-diffusion model of within-host viral infections. We introduce the basic reproduction number <mml:math xmlns:mml="https://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>R</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> and show that the infection-free steady state is globally asymptotically stable when <mml:math xmlns:mml="https://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>R</mml:mi> <mml:mn>0</mml:mn></mml:msub> <mml:mo>≤</mml:mo> <mml:mn>1</mml:mn></mml:mrow> </mml:math> , while the disease is uniformly persistent when <mml:math xmlns:mml="https://www...

This study discusses non-steady effects encountered in peristaltic flows in oesophagus. The purpose of this communication is to evolve a mechanism to diagnose tumor in an oesophagus mathematically. The tumor is modelled by generic bump function of certain height and width. The method of solution follows long wavelength and low-Reynolds number approximations for unsteady flow, while integrations have been performed numerically in order to plot graphs, which reveal various characteristics of the flow. The goal is to assess how pressure varies across the tumor's width...

Biochemical covalent modification networks exhibit a remarkable suite of steady state and dynamical properties such as multistationarity, oscillations, ultrasensitivity and absolute concentration robustness. This paper focuses on conditions required for a network of this type to have a species with absolute concentration robustness. We find that the robustness in a substrate is endowed by its interaction with a bifunctional enzyme, which is an enzyme that has different roles when isolated versus when bound as a substrate-enzyme complex...

We study an integro-difference equation model that describes the spatial dynamics of a species with a strong Allee effect in a shifting habitat. We examine the case of a shifting semi-infinite bad habitat connected to a semi-infinite good habitat. In this case we rigorously establish species persistence (non-persistence) if the habitat shift speed is less (greater) than the asymptotic spreading speed of the species in the good habitat. We also examine the case of a finite shifting patch of hospitable habitat, and find that the habitat shift speed must be less than the asymptotic spreading speed associated with the habitat and there is a critical patch size for species persistence...

We consider a deterministic model of HIV infection that involves macrophages as a long-term active reservoir to describe all three stages of the disease process: the acute stage, chronic infection, and the transition to AIDS. The proposed model is shown to retain crucial properties, such as the positivity of solutions, regardless of variations in model parameters. A dynamical analysis is performed to identify the local stability properties of the viral clearance steady state. This analysis illustrates how chronically infected macrophages can explain the progression to AIDS and provoke viral explosion, while previous models do not...

In the process of malaria transmission, natural recovery individuals are slightly infectious compared with infected individuals. Our concern is whether the infectivity of natural recovery category can be ignored in areas with limited medical resources, so as to reveal the epidemic pattern of malaria with simpler analysis. To achieve this, we incorporate saturated treatment into two-compartment and three-compartment models, and the infectivity of natural recovery category is only reflected in the latter. The non-spatial two-compartment model can admit backward bifurcation...

Collective cell migration is a multicellular phenomenon that arises in various biological contexts, including cancer and embryo development. 'Collectiveness' can be promoted by cell-cell interactions such as co-attraction and contact inhibition of locomotion. These mechanisms act on cell polarity, pivotal for directed cell motility, through influencing the intracellular dynamics of small GTPases such as Rac1. To model these dynamics we introduce a biased random walk model, where the bias depends on the internal state of Rac1, and the Rac1 state is influenced by cell-cell interactions and chemoattractive cues...

Fick's law and the Fokker-Planck law of diffusion are applied to manifest the cognitive dispersal of individuals in two reaction-diffusion SEIR epidemic models, where the disease transmission is illustrated by nonlocal infection mechanisms in heterogeneous environments. Building upon the well-posedness of solutions, threshold dynamics are discussed in terms of the basic reproduction numbers for the two cognitive epidemic models. The numerical investigation reveals that the Fokker-Planck law can better describe the diffusion of individuals by taking different dispersal strategies of exposed individuals in our cognitive epidemic models, and provides some insights on spatial segregation and nonpharmaceutical interventions: (i) spatial segregation occurs in the random diffusion model when the nonlocal infection radius is small, while it appears in the symmetric diffusion model when the radius is large; (ii) nonpharmaceutical interventions on restricting the dispersal of exposed and infected individuals do not contribute to reducing the infection proportion, but rather eliminate the disease in a region, which expands as the nonlocal infection radius increases...

Ontogenic resistance has been described for many plant-pathogen systems. Conversely, coffee leaf rust, a major fungal disease that drastically reduces coffee production, exhibits a form of ontogenic susceptibility, with a higher infection risk for mature leaves. To take into account stage-dependent crop response to phytopathogenic fungi, we developed an SEIR-U epidemiological model, where U stands for spores, which differentiates between young and mature leaves. Based on this model, we also explored the impact of ontogenic resistance on the sporulation rate...

Reticulations in a phylogenetic network represent processes such as gene flow, admixture, recombination and hybrid speciation. Extending definitions from the tree setting, an anomalous network is one in which some unrooted tree topology displayed in the network appears in gene trees with a lower frequency than a tree not displayed in the network. We investigate anomalous networks under the Network Multispecies Coalescent Model with possible correlated inheritance at reticulations. Focusing on subsets of 4 taxa, we describe a new algorithm to calculate quartet concordance factors on networks of any level, faster than previous algorithms because of its focus on 4 taxa...

Agent-based models (ABMs) are readily used to capture the stochasticity in tumour evolution; however, these models are often challenging to validate with experimental measurements due to model complexity. The Voronoi cell-based model (VCBM) is an off-lattice agent-based model that captures individual cell shapes using a Voronoi tessellation and mimics the evolution of cancer cell proliferation and movement. Evidence suggests tumours can exhibit biphasic growth in vivo. To account for this phenomena, we extend the VCBM to capture the existence of two distinct growth phases...