Immuno-epidemiological Model Research Articles

An immuno-epidemiological model with non-exponentially distributed disease stage on complex networks

Most of epidemic models assume that duration of the disease phase is distributed exponentially for the simplification of model formulation and analysis. Actually, the exponentially distributed assumption on the description of disease stages is hard to accurately approximate the interplay of drug concentration and viral load within host. In this article, we formulate an immuno-epidemiological epidemic model on complex networks, which is composed of ordinary differential equations and integral equations. The linkage of within- and between-host is connected by setting that the death caused by the disease is an increasing function in viral load within host. Mathematical analysis of the model includes the existence of the solution to the epidemiological model on complex networks, the existence and stability of equilibrium, which are completely determined by the basic reproduction number of the between-host system. Numerical analysis are shown that the non-exponentially distributions and the topology of networks have significant roles in the prediction of epidemic patterns.

Modeling and analysis of a two-strain immuno-epidemiological model with reinfection

In this paper, we formulate a two-strain model with reinfection that combines immunological and epidemiological dynamics across scales, using the COVID-19 pandemic as a case study. Firstly, we conduct a qualitative analysis of both within-host and between-host models. For the within-host model, we prove the existence and stability of equilibria, and Hopf bifurcation occurs from the infection equilibrium with immune response. This implies that, under specific immune states, the virus within an infected individual may persist, and its concentration may also oscillate periodically. For the between-host model, the disease-free equilibrium always exists and is locally asymptotically stable when the epidemiological basic reproduction number ℜ0<1. In addition, the model can have boundary equilibria of strain 1 or strain 2, which are locally asymptotically stable under specific conditions. However, the co-existence equilibrium does not exist. Secondly, to explore the infection and transmission mechanisms of two strain models and obtain reliable parameter values, we utilize statistical data to fit the immuno-epidemiological model. Simultaneously, we conduct an identifiability analysis of the immuno-epidemiological model to ensure the robustness of the fitted parameters. The results demonstrate the reliable estimation of parameter ranges for structurally unidentifiable parameters with minor measurement errors using the affine invariant ensemble Markov Chain Monte Carlo algorithm (GWMCMC). Moreover, simulations illustrate that enhancing treatment of patients infected with BA.2 strains to inhibit the number of viruses released by infected cells can significantly reduce disease spread.

An immuno-epidemiological model with waning immunity after infection or vaccination

In epidemics, waning immunity is common after infection or vaccination of individuals. Immunity levels are highly heterogeneous and dynamic. This work presents an immuno-epidemiological model that captures the fundamental dynamic features of immunity acquisition and wane after infection or vaccination and analyzes mathematically its dynamical properties. The model consists of a system of first order partial differential equations, involving nonlinear integral terms and different transfer velocities. Structurally, the equation may be interpreted as a Fokker-Planck equation for a piecewise deterministic process. However, unlike the usual models, our equation involves nonlocal effects, representing the infectivity of the whole environment. This, together with the presence of different transfer velocities, makes the proved existence of a solution novel and nontrivial. In addition, the asymptotic behavior of the model is analyzed based on the obtained qualitative properties of the solution. An optimal control problem with objective function including the total number of deaths and costs of vaccination is explored. Numerical results describe the dynamic relationship between contact rates and optimal solutions. The approach can contribute to the understanding of the dynamics of immune responses at population level and may guide public health policies.

Optimal control of a multi-scale HIV-opioid model.

In this study, we apply optimal control theory to an immuno-epidemiological model of HIV and opioid epidemics. For the multi-scale model, we used four controls: treating the opioid use, reducing HIV risk behaviour among opioid users, entry inhibiting antiviral therapy, and antiviral therapy which blocks the viral production. Two population-level controls are combined with two within-host-level controls. We prove the existence and uniqueness of an optimal control quadruple. Comparing the two population-level controls, we find that reducing the HIV risk of opioid users has a stronger impact on the population who is both HIV-infected and opioid-dependent than treating the opioid disorder. The within-host-level antiviral treatment has an effect not only on the co-affected population but also on the HIV-only infected population. Our findings suggest that the most effective strategy for managing the HIV and opioid epidemics is combining all controls at both within-host and between-host scales.

Final and peak epidemic sizes of immuno-epidemiological SIR models

Final and peak epidemic sizes of immuno-epidemiological SIR models

Rich dynamics of a bidirectionally linked immuno-epidemiological model for cholera.

Cholera is an environmentally driven disease where the human hosts both ingest the pathogen from polluted environment and shed the pathogen to the environment, generating a two-way feedback cycle. In this paper, we propose a bidirectionally linked immuno-epidemiological model to study the interaction of within- and between-host cholera dynamics. We conduct a rigorous analysis for this multiscale model, with a focus on the stability and bifurcation properties of each feasible equilibrium. We find that the parameter that represents the bidirectional connection is a key factor in shaping the rich dynamics of the system, including the occurrence of the backward bifurcation and Hopf bifurcation. Numerical results illustrate a practical application of our model and add new insight into the prevention and intervention of cholera epidemics.

Stability Switches, Hopf Bifurcation and Chaotic Dynamics in Simple Epidemic Model with State-Dependent Delay

Dynamic behavior investigations of infectious disease models are central to improve our understanding of emerging characteristics of model states interaction. Here, we consider a Susceptible-Infected (SI) model with a general state-dependent delay, which covers an immuno-epidemiological model of pathogen transmission, developed in our early study, using a threshold delay to examine the effects of multiple exposures to a pathogen. The analysis in the previous work showed the appearance of forward as well as backward bifurcations of endemic equilibria when the basic reproductive ratio [Formula: see text] is less than unity. The analysis, in the present work, of the endemically infected equilibrium behavior, through the study of a second order exponential polynomial characteristic equation, concludes the existence of a Hopf bifurcation on the upper branch of the backward bifurcation diagram and gives the criteria for stability switches. Furthermore, the inclusion of state-dependent delays is shown to entirely change the dynamics of the SI model and give rise to rich behaviors including periodic, torus and chaotic dynamics.

An immuno-epidemiological model for transient immune protection: A case study for viral respiratory infections.

The dynamics of infectious disease in a population critically involves both within-host pathogen replication and between host pathogen transmission. While modeling efforts have recently explored how within-host dynamics contribute to shaping population transmission, fewer have explored how ongoing circulation of an epidemic infectious disease can impact within-host immunological dynamics. We present a simple, influenza-inspired model that explores the potential for re-exposure during a single, ongoing outbreak to shape individual immune response and epidemiological potential in non-trivial ways. We show how even a simplified system can exhibit complex ongoing dynamics and sensitive thresholds in behavior. We also find epidemiological stochasticity likely plays a critical role in reinfection or in the maintenance of individual immunological protection over time.

An age-dependent immuno-epidemiological model with distributed recovery and death rates.

The work is devoted to a new immuno-epidemiological model with distributed recovery and death rates considered as functions of time after the infection onset. Disease transmission rate depends on the intra-subject viral load determined from the immunological submodel. The age-dependent model includes the viral load, recovery and death rates as functions of age considered as a continuous variable. Equations for susceptible, infected, recovered and dead compartments are expressed in terms of the number of newly infected cases. The analysis of the model includes the proof of the existence and uniqueness of solution. Furthermore, it is shown how the model can be reduced to age-dependent SIR or delay model under certain assumptions on recovery and death distributions. Basic reproduction number and final size of epidemic are determined for the reduced models. The model is validated with a COVID-19 case data. Modelling results show that proportion of young age groups can influence the epidemic progression since disease transmission rate for them is higher than for other age groups.

An immuno-epidemiological model linking between-host and within-host dynamics of cholera.

Cholera, a severe gastrointestinal infection caused by the bacterium Vibrio cholerae, remains a major threat to public health, with a yearly estimated global burden of 2.9 million cases. Although most existing models for the disease focus on its population dynamics, the disease evolves from within-host processes to the population, making it imperative to link the multiple scales of the disease to gain better perspectives on its spread and control. In this study, we propose an immuno-epidemiological model that links the between-host and within-host dynamics of cholera. The immunological (within-host) model depicts the interaction of the cholera pathogen with the adaptive immune response. We distinguish pathogen dynamics from immune response dynamics by assigning different time scales. Through a time-scale analysis, we characterise a single infected person by their immune response. Contrary to other within-host models, this modelling approach allows for recovery through pathogen clearance after a finite time. Then, we scale up the dynamics of the infected person to construct an epidemic model, where the infected population is structured by individual immunological dynamics. We derive the basic reproduction number ($ \mathcal{R}_0 $) and analyse the stability of the equilibrium points. At the disease-free equilibrium, the disease will either be eradicated if $ \mathcal{R}_0 < 1 $ or otherwise persists. A unique endemic equilibrium exists when $ \mathcal{R}_0 > 1 $ and is locally asymptotically stable without a loss of immunity.

The Enemy of my Enemy is my Friend: Immune-Mediated Facilitation Contributes to Fitness of Co-Infecting Helminths.

The conceptual understanding of immune-mediated interactions between parasites is rooted in the theory of community ecology. One of the limitations of this approach is that most of the theory and empirical evidence has focused on resource or immune-mediated competition between parasites and yet there is ample evidence of positive interactions that could be generated by immune-mediated facilitation. We developed an immuno-epidemiological model and applied it to longitudinal data of two gastrointestinal helminths in two rabbit populations to investigate, through model testing, how immune-mediated mechanisms of parasite regulation could explain the higher intensities of both helminths in rabbits with dual than single infections. The model framework was selected and calibrated on rabbit population A and then validated on the nearby rabbit population B to confirm the consistency of the findings and the generality of the mechanisms. Simulations suggested that the higher intensities in rabbits with dual infections could be explained by a weakened or low species-specific IgA response and an asymmetrical IgA cross-reaction. Simulations also indicated that rabbits with dual infections shed more free-living stages that survived for longer in the environment, implying greater transmission than stages from hosts with single infections. Temperature and humidity selectively affected the free-living stages of the two helminths. These patterns were comparable in the two rabbit populations and support the hypothesis that immune-mediated facilitation can contribute to greater parasite fitness and local persistence.

An Epidemic Model with Time Delay Determined by the Disease Duration

Immuno-epidemiological models with distributed recovery and death rates can describe the epidemic progression more precisely than conventional compartmental models. However, the required immunological data to estimate the distributed recovery and death rates are not easily available. An epidemic model with time delay is derived from the previously developed model with distributed recovery and death rates, which does not require precise immunological data. The resulting generic model describes epidemic progression using two parameters, disease transmission rate and disease duration. The disease duration is incorporated as a delay parameter. Various epidemic characteristics of the delay model, namely the basic reproduction number, the maximal number of infected, and the final size of the epidemic are derived. The estimation of disease duration is studied with the help of real data for COVID-19. The delay model gives a good approximation of the COVID-19 data and of the more detailed model with distributed parameters.

Sensitivity Analysis in an Immuno-Epidemiological Vector-Host Model.

Sensitivity Analysis (SA) is a useful tool to measure the impact of changes in model parameters on the infection dynamics, particularly to quantify the expected efficacy of disease control strategies. SA has only been applied to epidemic models at the population level, ignoring the effect of within-host virus-with-immune-system interactions on the disease spread. Connecting the scales from individual to population can help inform drug and vaccine development. Thus the value of understanding the impact of immunological parameters on epidemiological quantities. Here we consider an age-since-infection structured vector-host model, in which epidemiological parameters are formulated as functions of within-host virus and antibody densities, governed by an ODE system. We then use SA for these immuno-epidemiological models to investigate the impact of immunological parameters on population-level disease dynamics such as basic reproduction number, final size of the epidemic or the infectiousness at different phases of an outbreak. As a case study, we consider Rift Valley Fever Disease utilizing parameter estimations from prior studies. SA indicates that 1% increase in within-host pathogen growth rate can lead up to 8% increase in mathcal R_0, up to 1 % increase in steady-state infected host abundance, and up to 4% increase in infectiousness of hosts when the reproduction number mathcal R_0 is larger than one. These significant increases in population-scale disease quantities suggest that control strategies that reduce the within-host pathogen growth can be important in reducing disease prevalence.

Immuno-epidemiological co-affection model of HIV infection and opioid addiction.

In this paper, we present a multi-scale co-affection model of HIV infection and opioid addiction. The population scale epidemiological model is linked to the within-host model which describes the HIV and opioid dynamics in a co-affected individual. CD4 cells and viral load data obtained from morphine addicted SIV-infected monkeys are used to validate the within-host model. AIDS diagnoses, HIV death and opioid mortality data are used to fit the between-host model. When the rates of viral clearance and morphine uptake are fixed, the within-host model is structurally identifiable. If in addition the morphine saturation and clearance rates are also fixed the model becomes practical identifiable. Analytical results of the multi-scale model suggest that in addition to the disease-addiction-free equilibrium, there is a unique HIV-only and opioid-only equilibrium. Each of the boundary equilibria is stable if the invasion number of the other epidemic is below one. Elasticity analysis suggests that the most sensitive number is the invasion number of opioid epidemic with respect to the parameter of enhancement of HIV infection of opioid-affected individual. We conclude that the most effective control strategy is to prevent opioid addicted individuals from getting HIV, and to treat the opioid addiction directly and independently from HIV.

A network immuno-epidemiological model of HIV and opioid epidemics.

In this paper, we introduce a novel multi-scale network model of two epidemics: HIV infection and opioid addiction. The HIV infection dynamics is modeled on a complex network. We determine the basic reproduction number of HIV infection, $ \mathcal{R}_{v} $, and the basic reproduction number of opioid addiction, $ \mathcal{R}_{u} $. We show that the model has a unique disease-free equilibrium which is locally asymptotically stable when both $ \mathcal{R}_{u} $ and $ \mathcal{R}_{v} $ are less than one. If $ \mathcal{R}_{u} > 1 $ or $ \mathcal{R}_{v} > 1 $, then the disease-free equilibrium is unstable and there exists a unique semi-trivial equilibrium corresponding to each disease. The unique opioid only equilibrium exist when the basic reproduction number of opioid addiction is greater than one and it is locally asymptotically stable when the invasion number of HIV infection, $ \mathcal{R}^{1}_{v_i} $ is less than one. Similarly, the unique HIV only equilibrium exist when the basic reproduction number of HIV is greater than one and it is locally asymptotically stable when the invasion number of opioid addiction, $ \mathcal{R}^{2}_{u_i} $ is less than one. Existence and stability of co-existence equilibria remains an open problem. We performed numerical simulations to better understand the impact of three epidemiologically important parameters that are at the intersection of two epidemics: $ q_v $ the likelihood of an opioid user being infected with HIV, $ q_u $ the likelihood of an HIV-infected individual becoming addicted to opioids, and $ \delta $ recovery from opioid addiction. Simulations suggest that as the recovery from opioid use increases, the prevalence of co-affected individuals, those who are addicted to opioids and are infected with HIV, increase significantly. We demonstrate that the dependence of the co-affected population on $ q_u $ and $ q_v $ are not monotone.

Immuno-Epidemiological Model-Based Prediction of Further Covid-19 Epidemic Outbreaks Due to Immunity Waning

We develop a new data-driven immuno-epidemiological model with distributed infectivity, recovery and death rates determined from the epidemiological, clinical and experimental data. Immunity in the population is taken into account through the time-dependent number of vaccinated people with different numbers of doses and through the acquired immunity for recovered individuals. The model is validated with the available data. We show that for the first time from the beginning of pandemic COVID-19 some countries reached collective immunity. However, the epidemic continues because of the emergence of new variant BA.2 with a larger immunity escape or disease transmission rate than the previous BA.l variant. Large epidemic outbreaks can be expected several months later due to immunity waning. These outbreaks can be restrained by an intensive booster vaccination.

Quantification of the Tradeoff between Test Sensitivity and Test Frequency in a COVID-19 Epidemic-A Multi-Scale Modeling Approach.

Control strategies that employ real time polymerase chain reaction (RT-PCR) tests for the diagnosis and surveillance of COVID-19 epidemic are inefficient in fighting the epidemic due to high cost, delays in obtaining results, and the need of specialized personnel and equipment for laboratory processing. Cheaper and faster alternatives, such as antigen and paper-strip tests, have been proposed. They return results rapidly, but have lower sensitivity thresholds for detecting virus. To quantify the effects of the tradeoffs between sensitivity, cost, testing frequency, and delay in test return on the overall course of an outbreak, we built a multi-scale immuno-epidemiological model that connects the virus profile of infected individuals with transmission and testing at the population level. We investigated various randomized testing strategies and found that, for fixed testing capacity, lower sensitivity tests with shorter return delays slightly flatten the daily incidence curve and delay the time to the peak daily incidence. However, compared with RT-PCR testing, they do not always reduce the cumulative case count at half a year into the outbreak. When testing frequency is increased to account for the lower cost of less sensitive tests, we observe a large reduction in cumulative case counts, from to as low as half a year into the outbreak. The improvement is preserved even when the testing budget is reduced by one half or one third. Our results predict that surveillance testing that employs low-sensitivity tests at high frequency is an effective tool for epidemic control.

Epidemiological and evolutionary considerations of SARS-CoV-2 vaccine dosing regimes.

Given vaccine dose shortages and logistical challenges, various deployment strategies are being proposed to increase population immunity levels to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Two critical issues arise: How timing of delivery of the second dose will affect infection dynamics and how it will affect prospects for the evolution of viral immune escape via a buildup of partially immune individuals. Both hinge on the robustness of the immune response elicited by a single dose as compared with natural and two-dose immunity. Building on an existing immuno-epidemiological model, we find that in the short term, focusing on one dose generally decreases infections, but that longer-term outcomes depend on this relative immune robustness. We then explore three scenarios of selection and find that a one-dose policy may increase the potential for antigenic evolution under certain conditions of partial population immunity. We highlight the critical need to test viral loads and quantify immune responses after one vaccine dose and to ramp up vaccination efforts globally.

A Network Immuno-Epidemiological HIV Model.

In this paper we formulate a multi-scale nested immuno-epidemiological model of HIV on complex networks. The system is described by ordinary differential equations coupled with a partial differential equation. First, we prove the existence and uniqueness of the immunological model and then establish the well-posedness of the multi-scale model. We derive an explicit expression of the basic reproduction number [Formula: see text] of the immuno-epidemiological model. The system has a disease-free equilibrium and an endemic equilibrium. The disease-free equilibrium is globally stable when [Formula: see text] and unstable when [Formula: see text]. Numerical simulations suggest that [Formula: see text] increases as the number of nodes in the network increases. Further, we find that for a scale-free network the number of infected individuals at equilibrium is a hump-like function of the within-host reproduction number; however, the dependence becomes monotone if the network has predominantly low connectivity nodes or high connectivity nodes.

Infection severity across scales in multi-strain immuno-epidemiological Dengue model structured by host antibody level

Infection by distinct Dengue virus serotypes and host immunity are intricately linked. In particular, certain levels of cross-reactive antibodies in the host may actually enhance infection severity leading to Dengue hemorrhagic fever (DHF). The coupled immunological and epidemiological dynamics of Dengue calls for a multi-scale modeling approach. In this work, we formulate a within-host model which mechanistically recapitulates characteristics of antibody dependent enhancement in Dengue infection. The within-host scale is then linked to epidemiological spread by a vector-host partial differential equation model structured by host antibody level. The coupling allows for dynamic population-wide antibody levels to be tracked through primary and secondary infections by distinct Dengue strains, along with waning of cross-protective immunity after primary infection. Analysis of both the within-host and between-host systems are conducted. Stability results in the epidemic model are formulated via basic and invasion reproduction numbers as a function of immunological variables. Additionally, we develop numerical methods in order to simulate the multi-scale model and assess the influence of parameters on disease spread and DHF prevalence in the population.