Infection-age structured epidemic models with behavior change or treatment

We formulate infection-age structured susceptible-infective-removed (SIR) models with behavior change or treatment of infections. Individuals change their behavior or have treatment after they are infected. Using infection age as a continuous variable, and dividing infectives into discrete groups with different infection stages, respectively, we formulate a partial differential equation model and an ordinary differential equation model with behavior change or treatment. We derive explicit formulas for the reproductive number by linear stability analysis of the infection-free equilibrium, and explicit formulas for the unique endemic equilibrium, when it exists, for both models. These formulas provide mathematical theoretical frameworks for analysis of impact of behavior change or treatment of infection to the transmission dynamics of infectious diseases. We study several special cases and provide sensitivity analysis for the reproductive numbers with respect to model parameters based on those formulas.     

Differential susceptibility epidemic models

We formulate compartmental differential susceptibility (DS) susceptible-infective-removed (SIR) models by dividing the susceptible population into multiple subgroups according to the susceptibility of individuals in each group. We analyze the impact of disease-induced mortality in the situations where the number of contacts per individual is either constant or proportional to the total population. We derive an explicit formula for the reproductive number of infection for each model by investigating the local stability of the infection-free equilibrium. We further prove that the infection-free equilibrium of each model is globally asymptotically stable by qualitative analysis of the dynamics of the model system and by utilizing an appropriately chosen Liapunov function. We show that if the reproductive number is greater than one, then there exists a unique endemic equilibrium for all of the DS models studied in this paper. We prove that the endemic equilibrium is locally asymptotically stable for the models with no disease-induced mortality and the models with contact numbers proportional to the total population. We also provide sufficient conditions for the stability of the endemic equilibrium for other situations. We briefly discuss applications of the DS models to optimal vaccine strategies and the connections between the DS models and predator-prey models with multiple prey populations or host-parasitic interaction models with multiple hosts are also given.This research was partially supported by the Department of Energy under contracts W-7405-ENG-36 and the Applied Mathematical Sciences Program KC-07-01-01.     

Mathematical analysis of delay differential equation models of HIV-1 infection

Models of HIV-1 infection that include intracellular delays are more accurate representations of the biology and change the estimated values of kinetic parameters when compared to models without delays. We develop and analyze a set of models that include intracellular delays, combination antiretroviral therapy, and the dynamics of both infected and uninfected T cells. We show that when the drug efficacy is less than perfect the estimated value of the loss rate of productively infected T cells, delta, is increased when data is fit with delay models compared to the values estimated with a non-delay model. We provide a mathematical justification for this increased value of delta. We also provide some general results on the stability of non-linear delay differential equation infection models.     

On local bifurcations in neural field models with transmission delays

Neural field models with transmission delays may be cast as abstract delay differential equations (DDE). The theory of dual semigroups (also called sun-star calculus) provides a natural framework for the analysis of a broad class of delay equations, among which DDE. In particular, it may be used advantageously for the investigation of stability and bifurcation of steady states. After introducing the neural field model in its basic functional analytic setting and discussing its spectral properties, we elaborate extensively an example and derive a characteristic equation. Under certain conditions the associated equilibrium may destabilise in a Hopf bifurcation. Furthermore, two Hopf curves may intersect in a double Hopf point in a two-dimensional parameter space. We provide general formulas for the corresponding critical normal form coefficients, evaluate these numerically and interpret the results.     

An alternative formulation for a delayed logistic equation

We derive an alternative expression for a delayed logistic equation, assuming that the rate of change of the population depends on three components: growth, death, and intraspecific competition, with the delay in the growth component. In our formulation, we incorporate the delay in the growth term in a manner consistent with the rate of instantaneous decline in the population given by the model. We provide a complete global analysis, showing that, unlike the dynamics of the classical logistic delay differential equation (DDE) model, no sustained oscillations are possible. Just as for the classical logistic ordinary differential equation (ODE) growth model, all solutions approach a globally asymptotically stable equilibrium. However, unlike both the logistic ODE and DDE growth models, the value of this equilibrium depends on all of the parameters, including the delay, and there is a threshold that determines whether the population survives or dies out. In particular, if the delay is too long, the population dies out. When the population survives, i.e., the attracting equilibrium has a positive value, we explore how this value depends on the parameters. When this value is positive, solutions of our DDE model seem to be well approximated by solutions of the logistic ODE growth model with this carrying capacity and an appropriate choice for the intrinsic growth rate that is independent of the initial conditions.     

Population equilibria and stability in plant-animal pollination systems

Differential equation models of the change in population size with change in time have provided valuable insight into nature and have served as tools for management of both predator-prey and competition systems. Presented here is an extension of the differential equation models for application to pollination ecology plant-animal interactions. Analyses define conditions for equilibria and for equilibrium stability for a one-plant one-animal pollinator system. Predictions about plant speciation and pollinator foraging behavior are discussed.     

Dynamics of malaria transmission model with sterile mosquitoes

In this paper, a malaria transmission model with sterile mosquitoes is considered. We first formulate a simple SEIR malaria transmission model as our baseline model. Then sterile mosquitoes are introduced into the baseline model. We consider the case that the release rate of sterile mosquitoes is proportional to the wild mosquito population size. To investigate the impact of releasing sterile mosquitoes on the malaria transmission, the dynamics of the baseline model and the models with the sterile mosquitoes are discussed. We derive formulas of the reproductive numbers and explore the existence of endemic equilibrium as the reproductive number is more than unity for these models. It is shown that both the baseline model and the models with the sterile mosquitoes undergo backward bifurcations. Based on theoretical analysis and numerical simulation, we investigate the impact of releasing sterile mosquitoes on malaria transmission.     

A PDMP model of the epithelial cell turn-over in the intestinal crypt including microbiota-derived regulations

Mechanistic models are a powerful tool to gain insights into biological processes. The parameters of such models, e.g. kinetic rate constants, usually cannot be measured directly but need to be inferred from experimental data. In this article, we study dynamical models of the translation kinetics after mRNA transfection and analyze their parameter identifiability. That is, whether parameters can be uniquely determined from perfect or realistic data in theory and practice. Previous studies have considered ordinary differential equation (ODE) models of the process, and here we formulate a stochastic differential equation (SDE) model. For both model types, we consider structural identifiability based on the model equations and practical identifiability based on simulated as well as experimental data and find that the SDE model provides better parameter identifiability than the ODE model. Moreover, our analysis shows that even for those parameters of the ODE model that are considered to be identifiable, the obtained estimates are sometimes unreliable. Overall, our study clearly demonstrates the relevance of considering different modeling approaches and that stochastic models can provide more reliable and informative results.

     

Structured and Unstructured Continuous Models for <Emphasis Type="Italic">Wolbachia</Emphasis> Infections

We introduce and investigate a series of models for an infection of a diplodiploid host species by the bacterial endosymbiont Wolbachia. The continuous models are characterized by partial vertical transmission, cytoplasmic incompatibility and fitness costs associated with the infection. A particular aspect of interest is competitions between mutually incompatible strains. We further introduce an age-structured model that takes into account different fertility and mortality rates at different stages of the life cycle of the individuals. With only a few parameters, the ordinary differential equation models exhibit already interesting dynamics and can be used to predict criteria under which a strain of bacteria is able to invade a population. Interestingly, but not surprisingly, the age-structured model shows significant differences concerning the existence and stability of equilibrium solutions compared to the unstructured model.     

The role of seasonality in the dynamics of deer tick populations

In this paper, we formulate a nonlinear system of difference equations that models the three-stage life cycle of the deer tick over four seasons. We study the effect of seasonality on the stability and oscillatory behavior of the tick population by comparing analytically the seasonal model with a non-seasonal one. The analysis of the models reveals the existence of two equilibrium points. We discuss the necessary and sufficient conditions for local asymptotic stability of the equilibria and analyze the boundedness and oscillatory behavior of the solutions. A main result of the mathematical analysis is that seasonality in the life cycle of the deer tick can have a positive effect, in the sense that it increases the stability of the system. It is also shown that for some combination of parameters within the stability region, perturbations will result in a return to the equilibrium through transient oscillations. The models are used to explore the biological consequences of parameter variations reflecting expected environmental changes.     

Phylogenetic distances for neighbour dependent substitution processes

We consider models of nucleotidic substitution processes where the rate of substitution at a given site depends on the state of the neighbours of the site. We first estimate the time elapsed between an ancestral sequence at stationarity and a present sequence. Second, assuming that two sequences are issued from a common ancestral sequence at stationarity, we estimate the time since divergence. In the simplest non-trivial case of a Jukes-Cantor model with CpG influence, we provide and justify mathematically consistent estimators in these two settings. We also provide asymptotic confidence intervals, valid for nucleotidic sequences of finite length, and we compute explicit formulas for the estimators and for their confidence intervals. In the general case of an RN model with YpR influence, we extend these results under a proviso, namely that the equation defining the estimator has a unique solution.     

Backward bifurcation and oscillations in a nested immuno-eco-epidemiological model

This paper introduces a novel partial differential equation immuno-eco-epidemiological model of competition in which one species is affected by a disease while another can compete with it directly and by lowering the first species' immune response to the infection, a mode of competition termed stress-induced competition. When the disease is chronic, and the within-host dynamics are rapid, we reduce the partial differential equation model (PDE) to a three-dimensional ordinary differential equation (ODE) model. The ODE model exhibits backward bifurcation and sustained oscillations caused by the stress-induced competition. Furthermore, the ODE model, although not a special case of the PDE model, is useful for detecting backward bifurcation and oscillations in the PDE model. Backward bifurcation related to stress-induced competition allows the second species to persist for values of its invasion number below one. Furthermore, stress-induced competition leads to destabilization of the coexistence equilibrium and sustained oscillations in the PDE model. We suggest that complex systems such as this one may be studied by appropriately designed simple ODE models.     

Amplitude Metrics for Cellular Circadian Bioluminescence Reporters

Bioluminescence rhythms from cellular reporters have become the most common method used to quantify oscillations in circadian gene expression. These experimental systems can reveal phase and amplitude change resulting from circadian disturbances, and can be used in conjunction with mathematical models to lend further insight into the mechanistic basis of clock amplitude regulation. However, bioluminescence experiments track the mean output from thousands of noisy, uncoupled oscillators, obscuring the direct effect of a given stimulus on the genetic regulatory network. In many cases, it is unclear whether changes in amplitude are due to individual changes in gene expression level or to a change in coherence of the population. Although such systems can be modeled using explicit stochastic simulations, these models are computationally cumbersome and limit analytical insight into the mechanisms of amplitude change. We therefore develop theoretical and computational tools to approximate the mean expression level in large populations of noninteracting oscillators, and further define computationally efficient amplitude response calculations to describe phase-dependent amplitude change. At the single-cell level, a mechanistic nonlinear ordinary differential equation model is used to calculate the transient response of each cell to a perturbation, whereas population-level dynamics are captured by coupling this detailed model to a phase density function. Our analysis reveals that amplitude changes mediated at either the individual-cell or the population level can be distinguished in tissue-level bioluminescence data without the need for single-cell measurements. We demonstrate the effectiveness of the method by modeling experimental bioluminescence profiles of light-sensitive fibroblasts, reconciling the conclusions of two seemingly contradictory studies. This modeling framework allows a direct comparison between in vitro bioluminescence experiments and in silico ordinary differential equation models, and will lead to a better quantitative understanding of the factors that affect clock amplitude.     

Effect of media-induced social distancing on disease transmission in a two patch setting

We formulate an SIS epidemic model on two patches. In each patch, media coverage about the cases present in the local population leads individuals to limit the number of contacts they have with others, inducing a reduction in the rate of transmission of the infection. A global qualitative analysis is carried out, showing that the typical threshold behavior holds, with solutions either tending to an equilibrium without disease, or the system being persistent and solutions converging to an endemic equilibrium. Numerical analysis is employed to gain insight in both the analytically tractable and intractable cases; these simulations indicate that media coverage can reduce the burden of the epidemic and shorten the duration of the disease outbreak.     

The differential infectivity and staged progression models for the transmission of HIV

Recent studies of HIV RNA in infected individuals show that viral levels vary widely between individuals and within the same individual over time. Individuals with higher viral loads during the chronic phase tend to develop AIDS more rapidly. If RNA levels are correlated with infectiousness, these variations explain puzzling results from HIV transmission studies and suggest that a small subset of infected people may be responsible for a disproportionate number of infections. We use two simple models to study the impact of variations in infectiousness. In the first model, we account for different levels of virus between individuals during the chronic phase of infection, and the increase in the average time from infection to AIDS that goes along with a decreased viral load. The second model follows the more standard hypothesis that infected individuals progress through a series of infection stages, with the infectiousness of a person depending upon his current disease stage. We derive and compare threshold conditions for the two models and find explicit formulas of their endemic equilibria. We show that formulas for both models can be put into a standard form, which allows for a clear interpretation. We define the relative impact of each group as the fraction of infections being caused by that group. We use these formulas and numerical simulations to examine the relative importance of different stages of infection and different chronic levels of virus to the spreading of the disease. The acute stage and the most infectious group both appear to have a disproportionate effect, especially on the early epidemic. Contact tracing to identify super-spreaders and alertness to the symptoms of acute HIV infection may both be needed to contain this epidemic.     

The effect of delay in viral production in within-host models during early infection

ABSTRACT

Delay in viral production may have a significant impact on the early stages of infection. During the eclipse phase, the time from viral entry until active production of viral particles, no viruses are produced. This delay affects the probability that a viral infection becomes established and timing of the peak viral load. Deterministic and stochastic models are formulated with either multiple latent stages or a fixed delay for the eclipse phase. The deterministic model with multiple latent stages approaches in the limit the model with a fixed delay as the number of stages approaches infinity. The deterministic model framework is used to formulate continuous-time Markov chain and stochastic differential equation models. The probability of a minor infection with rapid viral clearance as opposed to a major full-blown infection with a high viral load is estimated from a branching process approximation of the Markov chain model and the results are confirmed through numerical simulations. In addition, parameter values for influenza A are used to numerically estimate the time to peak viral infection and peak viral load for the deterministic and stochastic models. Although the average length of the eclipse phase is the same in each of the models, as the number of latent stages increases, the numerical results show that the time to viral peak and the peak viral load increase.     

Global Stability for Delay SIR and SEIR Epidemic Models with Nonlinear Incidence Rate

In this paper, based on SIR and SEIR epidemic models with a general nonlinear incidence rate, we incorporate time delays into the ordinary differential equation models. In particular, we consider two delay differential equation models in which delays are caused (i) by the latency of the infection in a vector, and (ii) by the latent period in an infected host. By constructing suitable Lyapunov functionals and using the Lyapunov–LaSalle invariance principle, we prove the global stability of the endemic equilibrium and the disease-free equilibrium for time delays of any length in each model. Our results show that the global properties of equilibria also only depend on the basic reproductive number and that the latent period in a vector does not affect the stability, but the latent period in an infected host plays a positive role to control disease development.     

具有时滞的HBV病毒动力学模型稳定性分析

HBV（HePatitis B virus）是一种具有严重传染性的肝炎病毒,迄今为止,人们对它的免疫和慢性化的机制等方面还不甚了解。本文基于相关的病理知识,对应的建立了具有时滞的微分方程数学模型,系统地探讨了肝炎B病毒与宿主细胞之间的关系,利用Lyapunov函数方法研究了病毒动力学模型感染平衡点的局部稳定性和未感染平衡点全局稳定性,并利用数学模拟验证了理论分析。结果表明时滞的存在不会影响到感染平衡点的局部稳定性,但能影响平衡点到达的时间跨度,对于药物治疗的疗程和治疗时机的确定有参考意义。     

An intuitive formulation for the reproductive number for the spread of diseases in heterogeneous populations

The thresholds for mathematical epidemiology models specify the critical conditions for an epidemic to grow or die out. The reproductive number can provide significant insight into the transmission dynamics of a disease and can guide strategies to control its spread. We define the mean number of contacts, the mean duration of infection, and the mean transmission probability appropriately for certain epidemiological models, and construct a simplified formulation of the reproductive number as the product of these quantities. When the spread of the epidemic depends strongly upon the heterogeneity of the populations, the epidemiological models must account for this heterogeneity, and the expressions for the reproductive number become correspondingly more complex. We formulate several models with different heterogeneous structures and demonstrate how to define the mean quantities for an explicit expression for the reproductive number. In complex heterogeneous models, it seems necessary to define the reproductive number for each structured subgroup or cohort and then use the average of these reproductive numbers weighted by their heterogeneity to estimate the reproductive number for the total population.