This article has Open Peer Review reports available.
A large-scale immuno-epidemiological simulation of influenza A epidemics
- Sarah Lukens^{1, 9}Email author,
- Jay DePasse^{2},
- Roni Rosenfeld^{3},
- Elodie Ghedin^{4},
- Ericka Mochan^{7},
- Shawn T Brown^{2},
- John Grefenstette^{6},
- Donald S Burke^{8},
- David Swigon^{1} and
- Gilles Clermont^{5}
https://doi.org/10.1186/1471-2458-14-1019
© Lukens et al.; licensee BioMed Central Ltd. 2014
Received: 28 July 2014
Accepted: 18 September 2014
Published: 29 September 2014
Abstract
Background
Agent based models (ABM) are useful to explore population-level scenarios of disease spread and containment, but typically characterize infected individuals using simplified models of infection and symptoms dynamics. Adding more realistic models of individual infections and symptoms may help to create more realistic population level epidemic dynamics.
Methods
Using an equation-based, host-level mathematical model of influenza A virus infection, we develop a function that expresses the dependence of infectivity and symptoms of an infected individual on initial viral load, age, and viral strain phenotype. We incorporate this response function in a population-scale agent-based model of influenza A epidemic to create a hybrid multiscale modeling framework that reflects both population dynamics and individualized host response to infection.
Results
At the host level, we estimate parameter ranges using experimental data of H1N1 viral titers and symptoms measured in humans. By linearization of symptoms responses of the host-level model we obtain a map of the parameters of the model that characterizes clinical phenotypes of influenza infection and immune response variability over the population. At the population-level model, we analyze the effect of individualizing viral response in agent-based model by simulating epidemics across Allegheny County, Pennsylvania under both age-specific and age-independent severity assumptions.
Conclusions
We present a framework for multi-scale simulations of influenza epidemics that enables the study of population-level effects of individual differences in infections and symptoms, with minimal additional computational cost compared to the existing population-level simulations.
Background
Models of infectious disease epidemics, such as due to influenza A virus, have proven useful for understanding dynamics of disease transmission, vaccination strategies, and are frequently used as a decision support tool by public health officials. Single scale models of influenza A virus have been developed both on the host-level and on the population level in order to study epidemiological, immune response and viral characteristics of influenza A infection [1].
On the population scale, continuous models based on SIR (susceptible/infected/removed) approaches have been used for large-scale epidemiological predictions [2]. Continuous SIR models assume homogeneous populations and any model assumption about heterogeneity (such as human characteristics, network structure, and environmental factors) increases the number of equations and parameters, necessitating changing the model structure for any additional assumption. An alternative discrete approach to the study of disease epidemics is through agent-based models (ABM), where individuals are represented as autonomous agents whose infectious status is followed in time. Population density, age-structure, travel patterns and inter-individual contact patterns are derived or inferred from available census and demographic data. The infection propagates in the population according to disease transmission and duration rules, where inter-individual variability is represented in the form of statistical distributions. Because ABMs evolve individual agents, simulations require considerable computational cost. Many large-scale collaborative networks have developed and use large platform ABMs to make epidemiological predictions [3–6] and the impact of a variety of mitigation strategies. ABMs have been used in studies of school closure strategies [7, 8], in determining the role of subway travel in an epidemic [9], in the role of presenteeism on disease transmission in the workplace [10], and in vaccine allocation [11, 12].
Immuno-epidemiological models explore how host-level immune response affects population-level epidemiological patterns [13] and allow much finer tuning of host-behavior, compared to population-level models. The nested modeling approach embeds a mechanistic model of host-pathogen disease dynamics into an epidemiological model of infection by linking epidemiological parameters such as transmission rate or duration of infectiousness [14–16]. Immuno-epidemiological models consider assumptions about acquired immunity to study spread of infection [17–19] or host-parasite co-evolution [14, 20]. To account for heterogeneity in immune response, SIR models may be refined to incorporate individual host response through different linking mechanisms. Understanding the biological mechanisms underlying these variations in host response is important for estimating disease parameters and for designing optimal prevention strategies [21, 22].
Symptoms of influenza infection include fever, sore throat, chills, and cough, and have been attributed to the immune response to viral infection [23–26]. Clinical markers related to disease and disease-severity, largely due to host-immunity factors and prior exposures, are referred to as clinical phenotypes. Mathematical models of within-host influenza A infection have been used in estimation of kinematic parameters such as infection and clearance rates of the virus, as well as in prediction of dynamics and drug therapy strategies for individual infection [27, 28], yet have not been used to provide a mechanistic basis for varying disease phenotypes. Immuno-epidemiological modeling could be used to account for variability in clinical phenotype and evaluate the impact of symptom-based, biologically-based, and potentially behavior-based mitigation strategies at the population level.
A key challenge in deriving a realistic description of within-host responses is lack of data at the host-level. Human data is available from volunteer challenge studies where volunteers are inoculated intranasally, tend to be young (ages 18-35), healthy, and pre-screened for existing immunity to the experimental IAV strain. For H3N2 it has been shown that illness following intranasal inoculation is milder than illness following aerosol inoculation, manifested by shorter duration of cough and fever [29]. Mechanisms of transmission of naturally acquired infection include contact transmission, droplet spray transmission, and aerosol transmission. Therefore, using a host-level model to predict individual response to influenza infection requires altering model parameters based on healthy volunteer data in a biologically meaningful way.
Morbidity and mortality resulting from influenza infection is highest among the elderly (>65) and children. In a study of pandemic H1N1 in Ontario Canada, recovery was faster among patients under 18 years of age [30]. To gain qualitative insight of how host-level dynamics are varied as a function of age, animal studies prove to be useful. A study between aged and adult mice infected with sublethal doses of influenza virus (A/Puerto Rico/8/1934) showed slower recovery and delay in immune system activation and virus clearance in aged mice [31]. Differences in mouse response as a function of age provide insight in how to vary baseline responses in a within-host model across different age groups, while synthetic population data sets provide reasonable age distributions across a population [32, 33].
Naive integration of a within-host model in a populationlevel ABM requires the evaluation of a nonlinear system of differential equations for each individual infected, which is computationally prohibitive for realistically large simulations. In this manuscript, we present a prototypical hybrid immuno-epidemiological ABM model linking an equation-based within-host model and providing mechanistically-based host variability to an agent-based population level model. For outputs of the within-host model we choose the infectivity and symptoms scores during the time course of infection, which reflect likelihood of transmission and stay-at-home behavior of an infected individual. In this fashion, we can replicate existing results from simulated epidemics, but we introduce additional flexibility in studying variation of host-response at no additional computational cost. To eliminate computational cost, we replace the exact response with an approximation obtained by analyzing the most significant response modes of the within-host model. Further, the within-host model is simple to program and implement, and can be used to represent within-host dynamics in a variety of population-level platforms. We apply the hybrid ABM model to study questions arising from individual heterogeneity that cannot be addressed using a single scale model. In particular, we examine the role of varying severity of disease as a function of age on the dynamics of simulated epidemics. We find that epidemiological estimates such as attack rate and incidence are reduced when differences in host response across age groups are considered. Use of the immuno-epidemiological ABM model opens up many possibilities regarding the realistic impact of influenza on populations with different contact patterns or on mitigation strategies which depend on host phenotype or behavior.
Methods
Human volunteer data
The data set consists of viral load and symptom data from five previously published human volunteer studies conducted between 1995 and 1999 for 84 total individuals infected with A/Texas/36/91 (H1N1) [23, 24, 34–36]. Of these studies, only one [35] (n =17) does not contain symptom data. The viral titers across the five studies are all reported as averages across individuals in the respective studies in units of l o g _{10} TCID _{50}/mL. To adjust the data across studies, we used mixed-effects modeling to compute a vertical shift for each study (random-effect), assuming that the distance between each curve is minimized. We computed weighted means and standard deviations (fixed-effect) on the shifted data.
In all studies, symptom scores were measured twice a day, and averaged. Across studies, 65 of the 67 individuals with symptoms data experienced at least one symptom. While similar symptoms were scored, their measurements were on different scales across studies. As done in [37], we scaled each study curve to its maximum clinical score. The summary curves for both titer and symptoms were calculated as a weighted average with weights as the number of individuals considered in each study. Data is summarized in the Additional file 1.
Overview of the immuno-epidemiological ABM model
Estimation of parameters for virus/symptoms model
A schematic of the model is shown in the left panel of Figure 1(a). An individual is infected by an initial viral load V(0)=V _{0}, which infects target cells (H) of the respiratory tract with infection rate β. Virus replicates within latent phase cells (J), which after average rate of time k become infectious (I) and shed virus at a rate p. Infected cells die at rate δ, primarily through cell necrosis, which triggers productions of cytokines such as type I interferons (IFN), interleukin(IL)-6, tumor necrosis factor (TNF)- α and IL-1 [39]. The novel contribution of this model is the inclusion of a variable Ψ quantifying the intensity of systemic symptoms of influenza, such as malaise and fever. Of all the variables of the model, it is most reasonable to have systemic symptoms dependent on infected cells, since both flu-like symptoms and severity of systemic symptoms have been shown to be correlated to cytokine levels [40, 41], particularly IL-6 and IFN- α/β[23, 24, 42], and inflammatory cytokines are, in turn, produced by infected cells. It has also been noticed experimentally that the number of days of viral shedding is correlated with cytokine levels, severity of symptoms, and mucus weight [43, 44]. While resolution of symptoms is a complex process, we assume a constant decay rate a. Model trajectories include viral load and symptoms scores, which are used as determinants of infectivity and symptomaticity in the population level simulation and represent state variables corresponding to data. We use Bayesian inference to construct an ensemble of models consistent with the data in order to explore the parameter space associated with our model and data. We sample the Bayesian posterior density P(α|D) using a Markov Chain Monte Carlo (MCMC) method with parallel tempering [45–48]. Ensemble modeling results in a collection of trajectories consistent with the data, in which each trajectory corresponds to a parameter set from the posterior distribution.
Abstraction of the within-host model for population-level evaluation
Population level model
The within-host representation is linked to the Framework for Replication of Epidemiological Dynamics (FRED) ABM software, an open source population-level large scale modeling system developed by the University of Pittsburgh Models of infectious Disease Agents Study (MIDAS) Center of Excellence. We use data from the synthetic population database, which is freely available and based on 2005-2009 U.S. census data [32] from Allegheny County, Pennsylvania. The data from the 1,164,879 agents assigns characteristics and behaviors such as occupation and workplace or school [32, 33] to each individual. Each household is assigned a latitude/longitude coordinate, income, size, and the sex, race, and age of household occupants, representing the distribution of U.S. households. Behavior and daily routine of agents are based on a set of assumptions that describe contact patterns and movements of individuals. Individuals interact daily with schools or workplaces, household, and neighborhood networks. In each scenario, individuals contact infection based on a fixed mean number of people per day. For detailed descriptions on how agents interact and perform daily routines, the reader is referred to [4, 7, 9–12, 51].
At any time in FRED, each agent is in one of four states: susceptible (S), exposed (E), infectious (I) and recovered (R). We assume that individuals transition from the exposed state to the infectious state when the computed infectivity crosses a defined threshold. All agents are initially susceptible to infection, and at t=0 days, 100 random agents are initialized in the exposed state. Contact with an infectious person has an assigned probability of transmission, which is further modulated by the individual’s infectivity (See Additional file 1).
Results
Model ensemble
Parameter values
V _{0} | β | p | c | δ | k | θ | a | |
---|---|---|---|---|---|---|---|---|
Description | Inoculum | Infection | Viral | Viral | Infected | Latent | Rate of | Symptoms |
size | rate | production | clearance | cell decay | cell | symptoms | clearance | |
rate | rate | rate | rate | onset | rate | |||
Units | TCID _{50} mL | (TCID${}_{50}^{-1}$ mL day ^{-1}) | TCID${}_{50}^{-1}$ mL day ^{-1} | (day ^{-1}) | (day ^{-1}) | (day ^{-1}) | day ^{-1}[S]^{-1} | (day ^{-1}) |
Baseline (α _{0}) | 0.02 | 0.083 | 4000 | 3 | 1.3 | 5 | 2.5 | 0.55 |
Lower Bounds | 0.75e-5 | β _{0}/35 | p _{0}/35 | 8 | 0.5 | 2 | θ _{0}/35 | a _{0}/35 |
Upper Bounds | 1 | 35 β _{0} | 88 p _{0} | 40 | 2 | 6 | 35 θ _{0} | 35 a _{0} |
Best fit | 7.50e-6 | 0.0674 | 40356 | 8.00 | 1.364 | 3.684 | 2.75 | 0.498 |
Representing biological phenotypes
The second direction d _{2}=d _{2}(α) (Figure 4a) corresponds to onset dynamics which occur in the first two days of infection; virus that quickly reaches a high viral load corresponds to earlier symptoms than with those with a viral load that is slower to peak. Correlation between time of viral and symptoms peak was observed between the different experimental studies [23, 24, 34, 36], as well as between individual volunteers in an H3N2 study [52].
Within-host epidemiological parameters derived in a population-level model
As a starting point for the full population-scale model, we assign two uniform random numbers (x _{ j }∈U(0,1),j=1,2) to every individual in the population that becomes infected. These numbers correspond to coordinates in the response plane and are used to compute trajectories describing infectivity and symptomaticity of that agent during infection. We assume infection is symptomatic only if the respective trajectory is above a scalar threshold value, which also provide a metric for estimating the duration of latent and incubation periods in the model. The threshold values can be varied according to the particular purpose of the model, as we illustrate in the age-severity study below. Sample thresholds are plotted in Figure 5 as solid black lines.
Symptoms are assigned values 0 (asymptomatic) if below the threshold value or 1 (symptomatic) if above. The symptoms threshold S ^{∗}=0.717 is selected so that 33% of individuals fall below it in order to maintain a percentage of asymptomatic cases consistent with literature estimates [37]. The viral load threshold is selected at V ^{∗}=1.35, which gives a mean duration of infectivity of 5.2 days. The mean onset to infectivity calculated from evaluating the response surfaces is 0.83 days. This threshold was selected so the model would be comparable with estimates in the literature, and a comparison is found in the Discussion [6, 30, 43, 53–59].
Estimates of epidemiological parameters across response surfaces
Total | Symptomatic | Asymptomatic | |
---|---|---|---|
(st.d.) | (st.d.) | (st.d.) | |
Latent period (days) | 0.83 (0.40) | 0.69 (0.35) | 1.11 (0.34) |
Incubation period (days) | - | 2.14 (0.53) | None |
Duration of infectivity (days) | 5.17 (1.54) | 5.87 (1.38) | 3.71 (0.47) |
Duration of symptoms (days) | 4.18 (2.32) | 4.18 (2.32) | None |
Peak infectivity time (days) | 2.10 (0.51) | 2.00 (0.51) | 2.29 (0.45) |
Peak infectivity value | 3.15 (0.31) | 3.30 (0.22) | 2.85 (0.23) |
Peak symptoms time (days) | 3.04 (0.58) | 2.97 (0.60) | 3.19 (0.51) |
Peak symptoms score | 0.86 (0.25) | 0.99 (0.20) | 0.60 (0.075) |
Comparison of the baseline model with a reference ABM within-host model
We first show that the baseline within-host model produces similar results to a standard within-host model used in ABM studies, which we call the “reference” ABM model [4, 10]. In the reference within-host model in FRED, the mean duration of infectivity (4.1 days) and mean latent period (1.2 days) can be evaluated directly using the discrete probability distributions explicitly used in the model, and give an average total duration of infection as 5.3 days from inoculation to resolution [4]. Each infected agent is either symptomatic or not, based on a Bernoulli random variable.
In both cases in this comparison, all individuals who become infectious have an equal probability of contracting the same response, meaning that when an individual becomes infected they are assigned random numbers to represent illness independent of information such as a demographic or pre-existing conditions. Therefore the observed epidemiological differences are due only to the differences in the models of within-host response. The main difference between the two within-host models is that symptomatic individuals have a higher infectivity lasting over a shorter period and asymptomatic individuals have a lower infectivity lasting over a longer period in the reference ABM within-host model, compared to the new baseline within-host model.
Comparison of epidemiological parameters between reference ABM model and baseline
Reference ABM (st.d.) | Baseline (st.d.) | |
---|---|---|
Attack rate (AR) | 49.6588 (0.105) | 50.024 (0.109) |
Symptomatic AR | 33.271 (0.0879) | 33.6062 (0.083) |
% Symptomatic | 67.00% | 67.18% |
R _{0} | 2.14 (0.38) | 1.98 (0.37) |
Peak incidence (Days) | 31.12 | 25.10 |
Peak incidence (# People) | 36,811 | 41,883 |
Average run time | 80 seconds | 79 seconds |
Age-severity study
Aging alters both the innate and adaptive branches of the immune system, influencing response to influenza virus infection [60]. We assume that age of an infected individual maps to a scalar input value along d _{1}, reflecting longer and more severe infection and symptoms in older individuals. To illustrate the flexibility of our immune-epidemiological approach, we conduct a study where we compare the simulated epidemics between three different age-severity maps: a linear increasing function, and two functions based on mortality data from past epidemics in 1918 and 1957 [61]. Mortality patterns associated with epidemics traditionally form a “U” curve, reflecting higher mortality and more severe disease in the very young and the elderly. A distinctive feature of the 1918-19 influenza pandemic was a “W” shaped pattern in the mortality curve, showing a secondary peak in mortality rate for the 25-35 years age group [62, 63]. The age-severity model thus assumes mortality rates are indicators of severity of disease experienced by different age groups. Piece-wise linear functions map age to [0,1] for all three maps, and mortality data is scaled so the area under all three curves (for a g e∈[0,100]) is preserved. In this way the model considers two factors linking age with disease severity: deterioration of the immune system with age and immunological memory. Immunological memory refers to the adaptive immune response to pathogen the body has previously encountered, and contributes to explaining the “W” shape in the 1918 pandemic mortality curve [64].
Summary of demographics
Allegheny County | Franklin County | |
---|---|---|
Population size | 1,164,879 | 1,069,386 |
Land area (miles ^{2}) | 730.07 | 532.19 |
Median age | 41.3 | 33.4 |
Average household size | 2.23 | 2.38 |
Attack rates, reproductive number R _{ 0 } , and percent symptomatic for Allegheny and Franklin Counties
Baseline | Linear | 'U’ curve | 'W’ curve | |
---|---|---|---|---|
AR: Allegheny County (PA) | 50.02 (0.11) | 43.87 (0.15) | 42.66 (0.10) | 45.11 (0.12) |
AR: Franklin County (OH) | 54.84 (0.11) | 46.28 (0.12) | 45.98 (0.13) | 40.43 (0.11) |
R _{0}: Allegheny County (PA) | 2.14 (0.38) | 1.92 (0.32) | 1.93 (0.28) | 2.00 (0.33) |
R _{0}: Franklin County (OH) | 2.28 (0.34) | 2.05 (0.29) | 2.05 (0.27) | 2.16 (0.30) |
%S: Allegheny County (PA) | 67% | 54.23% | 49.25% | 62.16% |
%S: Franklin County (OH) | 67% | 48.21% | 44.69% | 63.16% |
We plot the fold difference in the attack rates in the age-severity models to their baselines in order to make a direct comparison between the models (Figure 7(d)). We observe a larger fold change for the W-model in Allegheny County, showing that for these counties, the decreased illness severity in the 45-65 ranges has a larger impact than the increased illness severity in the 25-35 age range. The larger fold change in the linear and U-shaped models of Franklin County implies a more dramatic change in epidemic outcome when age-severity assumptions are made in the younger age distribution scenario. The model assumes symptomatic cases by determining threshold maximum symptoms score. Less severe illness is more likely to correspond to asymptomatic cases, so the age-severity model also has tendencies to have an age dependence on symptomatic cases. We report the percentage of symptomatic patients in each scenario in Table 5.
Discussion
The strategy for incorporating the host-level model described here could be employed with other deterministic population level models as well as other well-established ABM simulation platforms such as EpiSimS [3]. Our aim here is to develop a methodology to bridge the gap between within-host models of virus/symptoms dynamics and population-level models that may also capture human movement and demographics. FRED provides a ready to use framework for population-level modeling using the synthetic population data set, allowing easy extension of the within-host model to simulations of populations with different demographic descriptions (here by distribution of age). We note that integration of our within-host model into a deterministic population level model may have similar results with less complexity, and we provide the coefficients for the response surfaces in the supplemental for straightforward use.
Our model implicitly uses age-structured contact patterns in the synthetic data set and ABM, and the results from the age-severity study reflect this relationship between distribution of age and epidemic severity. Previous studies have derived location specific contact matrices, either through use of synthetic populations with ABM network structures [67, 68], or from surveys of European communities [69]. In these cases, two mixing patterns emerge: a strong contact structure between people of similar age and a weaker structure between children and middle-aged adults, most likely due to parent-child contact. Use of contact matrices based on social networks better capture age-specific infection patterns of infectious diseases which helps in estimation of the basic reproductive number R _{0}. In our age-severity study, we provide a method to further fine-tune these estimates depending on host-response to a circulating strain. The computed values of R _{0} in this study are derived from assumptions about one particular experimental strain, and results may vary depending on characteristics unique to different circulating strains.
When integrating our derived within-host model into the population-level framework of FRED, further assumptions were required on scaling and thresholds. The symptoms threshold was set so approximately 1/3 of individuals in the baseline scenario are asymptomatic, consistent with accepted public health estimates [37]. However, for the infectivity threshold I ^{∗}, there was flexibility. Different values of I ^{∗}, correspond to different latent periods, impacting the duration of time an individual is in the exposed category. The threshold (1.35) was selected such that the average baseline latent periods are 0.83 days, and mean duration of infectivity is 5.17 days. Estimates for latent period are experimentally unknown due to difficulty in measuring, but a mathematical model predicts 0.4-1.5 days [53]. Duration of infectivity has been estimated at 4.8 days [37], and various studies have reported estimates at 3.38 days [30], 6.6 days [54], 5.0 days [55], and 3.1 days [43]. Our mean value of 4.98 days thus is consistent with existing estimates. This assumption reflects that in the baseline model, the incubation period exceeds the latent period, so there is a period in which an individual is not symptomatic and shedding virus. In one modeling study, a lag of 1.9 days between latent and incubation periods was assumed [6].
In our study, the incubation period is estimated at 2.14 days, duration of symptoms about 4 days for the symptomatic cases, in which both estimates rely on our choice of symptoms threshold. Incubation period has been estimated anywhere from 1.4 days (1.3 - 1.5 95% CI) [56], 1.5-2 days [57], 2-3 days [58], to 4 days [30]. Estimates for duration of symptoms vary from 4.5-5 days [37], 6 days [59], 5.6 days [55]. We remark that lowering the symptoms threshold would coincide with a shorter incubation period. One estimate that is well agreed upon is that viral load peaks at around 2 days post inoculation, and symptoms peak at around 3 days p.i. [37], which is consistent with our simulation measurements as well.
We selected the mathematical model of Baccam et al.[38] to describe viral and symptoms dynamics for three reasons. First, structural identifiability analysis showed parameter values could be estimated given the data. Second, it is the simplest model that generated the required inputs for a population level simulation (infectivity and symptoms score), and last, the model fit the data. Our method also provides a method for varying host-level responses from experimental baseline measurements in order to predict various illness scenarios under different host assumptions. As a starting point, we used a nominal response calibrated to one non-pandemic seasonal human influenza strain (A/Texas/36/91). As more data allows more subtle description of disease phenotypes characterized by more complex within-host models, our method could be easily scaled to such models as, although computation of the response-surfaces is more complex, it is done a priori and thus requires no additional resources for the population-level simulation.
Community based studies also give insight corresponding to our assumption that infection in the elderly may be characterized by longer viral load, and longer and higher symptoms scores. One study suggests a median duration of symptoms of 7 days for the entire community, compared to 8 days for older people [30]. Another study estimates that duration of infectiousness/viral shedding is 6.3 days for children, and 6.7 days for adults [54]. Using a linear mapping from age to infection severity is a simplistic assumption without biological motivation, which is why three different mappings were studied.
Our model estimates of duration and intensity of infectivity and symptoms show asymptomatic infections correspond to reduced viral load, which is consistent with previous studies of asymptomatic disease and transmission, but the degree is unknown [70]. In a community-based study of 2009 H1N1 in Hong Kong, peak viral load averaged over the asymptomatic subjects was 3.2e3 copies/mL compared to 3.6e7 copies/mL averaged over the symptomatic subjects [57]. In a recent human volunteer study of H3N2/Wisconsin, 50% of asymptomatic subjects had evident viral shedding and a reported a significantly reduced viral load in the asymptomatic cases [52]. However, the exact mechanism contributing to asymptomatic infection is not incorporated in this model, and a study of these mechanisms and how they influence an epidemic are important future considerations.
This study makes the assumption that infectivity scales linearly with contact rate and is proportional to the log10 of the viral load as measured by nasopharyngeal swabs. Although this is plausible, the exact relationship between viral load and infectivity on disease transmission is unknown. One model, which studies the relationship between virus/host dynamics and a population-level model, multiplied viral load by the amount of nasal discharge as the estimate of infectivity [71]. Another model of measles used the area under the infected cell curve to model transmission strength by 'amount’ of infection which was also used as a determinant for symptomatic or asymptomatic infection [72]. In a household study of influenza transmission in Hong Kong, models of molecular viral shedding and log10 molecular viral shedding were studied to compare viral load with transmission [57]. The case in which viral load is used as an indicator of transmission leads to nearly all transmission occurring in the first 1-2 days after inoculation, implying that intervention strategies must be rapid. Using the log10 model of viral shedding to transmission corresponds to significant transmission lasting strongly up to 3-4 days after viral inoculation. Therefore, knowing the exact relationship between infectivity and disease transmission would improve modeling and have significant implication for mitigation strategies.
In the mathematical model describing the host immune response, symptoms scores are generated from the biological mechanisms described by the ODE model, which greatly simplifies complex immune system dynamics. The immune system is composed of innate and adaptive components, which respond to viral infection through inflammation and creation of virus-specific antibodies and effector cells to help resolve the infection [73]. Previous models have been studied within the context of improved understanding of population-level dynamics have incorporated interactions of the innate immune system [53, 74]. Incorporating features such as inflammation and virus-specific antibody production would yield more realistic within-host dynamics, and could explicitly account for factors such as preexisting immunity. Importantly, more complex models could be attempted using response surface methodology at little additional computational cost for large-scale simulations. Further, the ensemble methodology allows the computation of the thickness of the response surface, which would add an additional stochastic element better representing our incomplete knowledge of individual responses.
The added flexibility of a mechanistic within-host to population-level models might be particularly relevant when evaluating the impact of host-level containment measures and additional factors known to demonstrate significant within-host effect depending on immune status. The ability to include this flexibility at minimal computational cost also offers exiting possibilities when simulating new variants, co-infection with bacterial pathogens, for providing a biological basis for morbidity, undoubtedly a major driver of agent behavior and cost of disease, and for mortality related to disease, which is currently lacking in existing population-level simulators.
Conclusions
We have developed a multi-scale immuno-epidemiological ABM model linking an equation-based within-host model to an agent-based population level model. The key innovation is the use of response surface analysis to simplify the computation of output trajectories of the equation-based within-host model, which makes the immuno-epidemiological model just as efficient and scalable as the underlying population-level ABM. The hybrid ABM model replicates results of simulated epidemics with simpler within-host assumptions such as the reference FRED within-host model, and provides additional flexibility to explore questions relevant to host heterogeneity at no extra computational cost.
Declarations
Acknowledgments
Research reported in this publication was supported by the National Institute Of General Medical Sciences of the National Institutes of Health under Award Number U54 GM088491, R01-GM83602, and NSF RTG DMS0739261. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
The authors wish to acknowledge David Galloway and Molly Eggleston at the Public Health Dynamics Laboratory in the University of Pittsburgh Graduate School of Public Health for their help during the preparation of this manuscript.
Authors’ Affiliations
References
- Murillo LN, Murillo MS, Perelson AS: Towards multiscale modeling of influenza infection. J Theor Biol. 2013, 332: 267-290.View ArticlePubMedPubMed CentralGoogle Scholar
- Keeling MJ, Rohani P: Modeling infectious diseases in humans and animals. 2008, New Jersey: Princeton University PressGoogle Scholar
- Eubank S, Guclu H, Kumar VA, Marathe MV, Srinivasan A, Toroczkai Z, Wang N: Modelling disease outbreaks in realistic urban social networks. Nature. 2004, 429 (6988): 180-184. 10.1038/nature02541.View ArticlePubMedGoogle Scholar
- Grefenstette JJ, Brown ST, Rosenfeld R, DePasse J, Stone NT, Cooley PC, Wheaton WD, Fyshe A, Galloway DD, Sriram A, Guclu H, Abraham T, Burke DS: FRED (A Framework for Reconstructing Epidemic Dynamics): an open-source software system for modeling infectious diseases and control strategies using census-based populations. BMC Public Health. 2013, 13: 940-10.1186/1471-2458-13-940.View ArticlePubMedPubMed CentralGoogle Scholar
- Halloran ME, Ferguson NM, Eubank S, Longini IM, Cummings DA, Lewis B, Xu S, Fraser C, Vullikanti A, Germann TC, Wagener D, Beckman R, Kadau K, Barrett C, Macken CA, Burke DS, Cooley P: Modeling targeted layered containment of an influenza pandemic in the United States. Proc Natl Acad Sci. 2008, 105 (12): 4639-4644. 10.1073/pnas.0706849105.View ArticlePubMedPubMed CentralGoogle Scholar
- Longini I, Nizam A, Xu S, Ungchusak K, Hanshaoworakul W, Cummings D, Halloran M: Containing pandemic influenza at the source. Science. 2005, 309 (5737): 1083-1087. 10.1126/science.1115717.View ArticlePubMedGoogle Scholar
- Lee B, Brown S, Cooley P, Potter M, Wheaton W, Voorhees R, Stebbins S, Grefenstette J, Zimmer S, Zimmerman R, Assi T-M, Bailey RR, Wagener DK, Burke DS: Simulating school closure strategies to mitigate an influenza epidemic. J Public Health Manag Pract: JPHMP. 2010, 16 (3): 252-10.1097/PHH.0b013e3181ce594e.View ArticlePubMedPubMed CentralGoogle Scholar
- Brown ST, Tai JH, Bailey RR, Cooley PC, Wheaton WD, Potter MA, Voorhees RE, LeJeune M, Grefenstette JJ, Burke DS, McGlone SM, Lee BY: Would school closure for the 2009 H1N1 influenza epidemic have been worth the cost?: A computational simulation of Pennsylvania. BMC Public Health. 2011, 11: 353-10.1186/1471-2458-11-353.View ArticlePubMedPubMed CentralGoogle Scholar
- Cooley P, Brown S, Cajka J, Chasteen B, Ganapathi L, Grefenstette J, Hollingsworth C, Lee B, Levine B, Wheaton WD, Wagener DK: The role of subway travel in an influenza epidemic: A New York City simulation. J Urban Health. 2011, 88 (5): 982-995. 10.1007/s11524-011-9603-4.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumar S, Grefenstette J, Galloway D, Albert S, Burke D: Policies to reduce influenza in the workplace: impact assessments using an agent-based model. Am J Public Health. 2013, 103 (8): 1406-1411. 10.2105/AJPH.2013.301269.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee B, Brown S, Korch G, Cooley P, Zimmerman R, Wheaton W, Zimmer S, Grefenstette J, Bailey R, Assi T-M, Burke DS: A computer simulation of vaccine prioritization, allocation, and rationing during the 2009 H1N1 influenza pandemic. Vaccine. 2010, 28 (31): 4875-4879. 10.1016/j.vaccine.2010.05.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee B, Brown S, Bailey R, Zimmerman R, Potter M, McGlone S, Cooley P, Grefenstette J, Zimmer S, Wheaton W, Quinn SC, Voorhees RE, Burke DS: The benefits to all of ensuring equal and timely access to influenza vaccines in poor communities. Health Aff. 2011, 30 (6): 1141-1150. 10.1377/hlthaff.2010.0778.View ArticleGoogle Scholar
- Hellriegel B: Immunoepidemiology–bridging the gap between immunology and epidemiology. Trends Parasitol. 2001, 17 (2): 102-106. 10.1016/S1471-4922(00)01767-0.View ArticlePubMedGoogle Scholar
- Gilchrist MA, Sasaki A: Modeling host–parasite coevolution: a nested approach based on mechanistic models. J Theor Biol. 2002, 218 (3): 289-308. 10.1006/jtbi.2002.3076.View ArticlePubMedGoogle Scholar
- Mideo N, Alizon S, Day T: Linking within-and between-host dynamics in the evolutionary epidemiology of infectious diseases. Trends Ecol Evol. 2008, 23 (9): 511-517. 10.1016/j.tree.2008.05.009.View ArticlePubMedGoogle Scholar
- Gandolfi A, Pugliese A, Sinisgalli C: Epidemic dynamics and host immune response: a nested approach. J Math Biol. 2014, 1: 1-37.Google Scholar
- Kostova T: Persistence of viral infections on the population level explained by an immunoepidemiological model. Math Biosci. 2007, 206 (2): 309-319. 10.1016/j.mbs.2005.08.003.View ArticlePubMedGoogle Scholar
- Heffernan J, Keeling MJ: Implications of vaccination and waning immunity. Proc R Soc B: Biol Sci. 2009, 276 (1664): 2071-2080. 10.1098/rspb.2009.0057.View ArticleGoogle Scholar
- Martcheva M, Pilyugin SS: An epidemic model structured by host immunity. J Biol Syst. 2006, 14 (02): 185-203. 10.1142/S0218339006001787.View ArticleGoogle Scholar
- Gilchrist MA, Coombs D: Evolution of virulence: interdependence, constraints, and selection using nested models. Theor Popul Biol. 2006, 69 (2): 145-153. 10.1016/j.tpb.2005.07.002.View ArticlePubMedGoogle Scholar
- Snider D, Bridges C, Weissman D: Meeting summary of the workshop 'approaches to better understand human influenza transmission. Approaches to better understand human influenza transmission. Centers for Disease Control and Prevention (United States of America). 4–5 November 2010; Atlanta, Georgia. 2010, [http://www.cdc.gov/influenzatransmissionworkshop2010/],Google Scholar
- Weinstein RA, Bridges CB, Kuehnert MJ, Hall CB: Transmission of influenza: implications for control in health care settings. Clin Infect Dis. 2003, 37 (8): 1094-1101. 10.1086/378292.View ArticleGoogle Scholar
- Hayden F, Fritz R, Lobo M, Alvord W, Strober W, Straus S: Local and systemic cytokine responses during experimental human influenze A virus infection. J Clin Invest. 1998, 101: 643-649. 10.1172/JCI1355.View ArticlePubMedPubMed CentralGoogle Scholar
- Fritz R, Hayden F, Calfee D: Nasal cytokine and chemokine responses in experimental influenza A virus infection: results of a placebo-controlled trial of intravenous zanamivir treatment. J Infect Dis. 1999, 180: 586-593. 10.1086/314938.View ArticlePubMedGoogle Scholar
- Julkunen I, Melen K, Nyqvist M, Pirhonen J, Sareneva T, Matikainen S: Inflammatory responses in influenza A virus infection. Vaccine. 2001, 19: S32-S37.View ArticleGoogle Scholar
- Eccles R: Understanding the symptoms of the common cold and influenza. Lancet Infect Dis. 2005, 5: 718-725. 10.1016/S1473-3099(05)70270-X.View ArticlePubMedGoogle Scholar
- Smith A, Perelson A: Influenza A virus infection kinetics: quantitative data and models. Wiley Interdisciplinary Rev: Syst Biol Med. 2011, 3: 429-445. 10.1002/wsbm.129. [http://dx.doi.org/10.1002/WSBM.129],Google Scholar
- Beauchemin C, Handel A: A review of mathematical models of influenza A infections within a host or cell culture: lessons learned and challenges ahead. BMC Public Health. 2011, 11 (Suppl 1): S7-10.1186/1471-2458-11-S1-S7.View ArticlePubMedPubMed CentralGoogle Scholar
- Little JW, Gordon RD, Hall WJ, Roth FK: Attenuated influenza produced by experimental intranasal inoculation. J Med Virol. 1979, 3 (3): 177-188. 10.1002/jmv.1890030303.View ArticlePubMedGoogle Scholar
- Tuite A, Greer A, Whelan M, Winter A, Lee B, Yan P, Wu J, Moghadas S, Buckeridge D, Pourbohloul B, Fisman DN: Estimated epidemiologic parameters and morbidity associated with pandemic H1N1 influenza. Can Med Assoc J. 2010, 182 (2): 131-136. 10.1503/cmaj.091807.View ArticleGoogle Scholar
- Toapanta F, Ross TM: Impaired immune responses in the lungs of aged mice following influenza infection. Respir Res. 2009, 10: 112-131. 10.1186/1465-9921-10-112.View ArticlePubMedPubMed CentralGoogle Scholar
- Wheaton W: 2005-2009 U.S. Synthetic Population Ver. 2. RTI International. 2012, [https://www.epimodels.org/midas/Rpubsyntdata1.do],Google Scholar
- Cajka J, Cooley P, Wheaton W: Attribute assignment to a synthetic population in support of agent-based disease modeling. Methods Rep (RTI Press). 2010, 19 (1009): 1-Google Scholar
- Hayden F, Treanor J, Betts RF, Lobo M, Esinhart JD, Hussey EK: Safety and efficacy of the neuraminidase inhibitor GG167 in experimental human influenza. J Am Med Ass. 1996, 275: 295-299. 10.1001/jama.1996.03530280047035.View ArticleGoogle Scholar
- Barroso L, Treanor J, Gubareva L, Hayden F: Efficacy and tolerability of the oral neuraminidase inhibitor peramivir in experimental human influenza: randomized, controlled trials for prophylaxis and treatment. Antivir Ther. 2005, 10: 901-910.PubMedGoogle Scholar
- Murphy A, Platts-Mills T, Lobo M, Hayden F: Respiratory nitric oxide levels in experimental human influenza. Chest. 1998, 114: 452-456. 10.1378/chest.114.2.452.View ArticlePubMedGoogle Scholar
- Carrat F, Vergu E, Ferguson N, Lemaitre M, Cauchemex S, Leach S, Valleron A: Time lines of infection and disease in human influenza: a review of volunteer challenge studies. Am J Epidemiol. 2008, 167: 775-785. 10.1093/aje/kwm375.View ArticlePubMedGoogle Scholar
- Baccam P, Beauchemin C, Macken C, Hayden F, Perelson A: Kinetics of influenza A virus infection in humans. J Virol. 2006, 80: 7590-7599. 10.1128/JVI.01623-05.View ArticlePubMedPubMed CentralGoogle Scholar
- La Gruta NL, Kedzierska K, Stambas J, Doherty PC: A question of self-preservation: immunopathology in influenza virus infection. Immunol Cell Biol. 2007, 85 (2): 85-92. 10.1038/sj.icb.7100026.View ArticlePubMedGoogle Scholar
- Conti B, Taberean I, Andrei C, Bartfai T: Cytokines and fever. Front Biosci. 2004, 9: 1433-1449. 10.2741/1341.View ArticlePubMedGoogle Scholar
- Descotes J, Vial T: Flu-like syndrome and cytokines. Cytokines in Human Health. 2007, New Jersey: Humana PressGoogle Scholar
- Cohen S, Doyle WJ, Skoner DP: Psychological stress, cytokine production, and severity of upper respiratory illness. Psychosom Med. 1999, 61 (2): 175-180. 10.1097/00006842-199903000-00009.View ArticlePubMedGoogle Scholar
- Skoner DP, Gentile DA, Patel A, Doyle WJ: Evidence for cytokine mediation of disease expression in adults experimentally infected with influenza A virus. J Infect Dis. 1999, 180: 10-14. 10.1086/314823.View ArticlePubMedGoogle Scholar
- Gentile D, Doyle W, Whiteside T, Fireman P, Hayden FG, Skoner D: Increased interleukin-6 levels in nasal lavage samples following experimental influenza A virus infection. Clin Diagn Lab Immunol. 1998, 5 (5): 604-608.PubMedPubMed CentralGoogle Scholar
- Metropolis N, Rosenbluth A, Rosenbluth M, Teller A, Teller E: Equations of state calculations by fast computing machines. J Chem Phys. 1953, 21: 1087-1092. 10.1063/1.1699114.View ArticleGoogle Scholar
- Hastings W: Monte Carlo sampling methods using Markov chains and their applications. Biometrika. 1970, 57: 97-109. 10.1093/biomet/57.1.97.View ArticleGoogle Scholar
- Gammerman D, Lopas H: Markov chain Monte Carlo stochastic simulation for Bayesian inference. 2006, Boca Raton, FL: Chapman & Hall/CRC, Taylor and Francis GroupGoogle Scholar
- Earl D, Deem M: Parallel tempering: theory, applictions and new perspectives. Phys Chem Chem Phys. 2005, 7: 3910-3916. 10.1039/b509983h.View ArticlePubMedGoogle Scholar
- Brown K, Sethna J: Statistical mechanical approaches to models with many poorly known parameters. Phys Rev E. 2003, 68: 021904-View ArticleGoogle Scholar
- Gutenkunst R, Waterfall J, Casey F, Brown K, Myers C, Sethna J: Universally sloppy parameter sensitivities in systems biology models. PLoS Comput Biol. 2007, 3: e189-10.1371/journal.pcbi.0030189.View ArticlePubMed CentralGoogle Scholar
- Lee H, Topham D, Park S, Hollenbaugh J, Treanor J, Mosmann T, Jin X, Ward B, Miao H, Holden-Wiltse J, Perelson A, Zand M, Wu H: Simulation and prediction of the adaptive immune response to Influenza A virus infection. J Virol. 2009, 83: 7151-7165. 10.1128/JVI.00098-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang Y, Zaas A, Rao A, Dobigeon N, Woolf P, Velman T, ien NO, McClain M, Varkey J, Nicholson B, Carin L, Kingsmore S, Woods C, Ginsburg G, III AH: Temporal dynamics of host molecular responses to differentiate symptomatic and asymptomatic influenza A infection. PLoS Genet. 2011, 7: e100234-Google Scholar
- Canini L, Carrat F: Population modeling of influenza A/H1N1 virus kinetics and symptom dynamics. J Virol. 2011, 85: 2764-2770. 10.1128/JVI.01318-10.View ArticlePubMedGoogle Scholar
- Seuss T, Buchholz U, Dupke S, Grunow R, an der Heiden M, Heider A, Biere B, Schweiger B, Haas W, Krause G: Shedding and transmission of novel influenza virus A/H1N1 infection in households - Germany, 2009. Am J Epidemiol. 2009, 11: 1157-64.Google Scholar
- Witkop C, Duffy M, Macias E, Gibbons T, Escobar J, Burwell K, Knight K: Novel influenza A (H1N1) outbreak at the US Air Force Academy: epidemiology and viral shedding duration. Am J Prev Med. 2010, 38 (2): 121-126. 10.1016/j.amepre.2009.10.005.View ArticlePubMedGoogle Scholar
- Lessler J, Reich N, Brookmeyer R, Perl T, Nelson K, Cummings D: Incubation periods of acute respiratory viral infections: a systematic review. Lancet Infect Dis. 2009, 9 (5): 291-300. 10.1016/S1473-3099(09)70069-6.View ArticlePubMedPubMed CentralGoogle Scholar
- Lau L, Cowling B, Chan KH, Lau E, Lipsitch M, Cheng C, Houck P, Uyeki T, Peiris J, Leung G: Viral shedding and clinical illness in naturally acquired Influenza A virus infections. J Infect Dis. 2010, 201: 1509-1516. 10.1086/652241.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamagishi T, Matsui T, Nakamura N, Oyama T, Taniguchi K, Aoki T, Hirakawa K, Okabe N: Onset and duration of symptoms and timing of disease transmission of 2009 influenza A (H1N1) in an outbreak in Fukuoka, Japan, June 2009. Jpn J Infect Dis. 2010, 63 (5): 327-31.PubMedGoogle Scholar
- Lessler J, Reich N, Cummings D: Outbreak of 2009 pandemic influenza A (H1N1) at a New York City school. N Engl J Med. 2009, 361 (27): 2628-2636. 10.1056/NEJMoa0906089.View ArticlePubMedGoogle Scholar
- Lambert ND, Ovsyannikova IG, Pankratz VS, Jacobson RM, Poland GA: Understanding the immune response to seasonal influenza vaccination in older adults: a systems biology approach. Expert Rev Vaccines. 2012, 11: 985-994. 10.1586/erv.12.61.View ArticlePubMedPubMed CentralGoogle Scholar
- Luk J, Gross P, Thompson WW: Observations on mortality during the 1918 influenza pandemic. Clin Infect Dis. 2001, 33 (8): 1375-1378. 10.1086/322662.View ArticlePubMedGoogle Scholar
- Dauer C, Serfling R: Mortality from influenza, 1957–1958 and 1959–1960. Am Rev Respir Dis. 1961, 83 (Suppl 2): 15-26.Google Scholar
- Palese P: Influenza: old and new threats. Nat Med. 2004, 10: S82-S87. 10.1038/nm1141.View ArticlePubMedGoogle Scholar
- Ahmed R, Oldstone MB, Palese P: Protective immunity and susceptibility to infectious diseases: lessons from the 1918 influenza pandemic. Nat Immunol. 2007, 8 (11): 1188-1193. 10.1038/ni1530.View ArticlePubMedGoogle Scholar
- Howden LM, Meyer JA: Age and sex composition: 2010:2010 census briefs. 2011, Washington, DC: U.S. Census Bureau, Available at: [http://www.census.gov/prod/cen2010/briefs/c2010br-03.pdf]Google Scholar
- Bureau, USCensus: American fact finder. Website. 2013, [http://www.census.gov/prod/cen2010/briefs/c2010br-03.pdf],Google Scholar
- Del Valle SY, Hyman J, Hethcote HW, Eubank SG: Mixing patterns between age groups in social networks. Social Netw. 2007, 29 (4): 539-554. 10.1016/j.socnet.2007.04.005.View ArticleGoogle Scholar
- Fumanelli L, Ajelli M, Manfredi P, Vespignani A, Merler S: Inferring the structure of social contacts from demographic data in the analysis of infectious diseases spread. PLoS Comput Biol. 2012, 8 (9): e1002673-10.1371/journal.pcbi.1002673.View ArticlePubMedPubMed CentralGoogle Scholar
- Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, Heijne J, Sadkowska-Todys M, Rosinka M, Edmunds WJ: Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5 (3): e74-10.1371/journal.pmed.0050074.View ArticlePubMedPubMed CentralGoogle Scholar
- Patrozou E, Mermel LA: Does influenza transmission occur from asymptomatic infection or prior to symptom onset?. Public Health Rep. 2009, 124 (2): 193-PubMedPubMed CentralGoogle Scholar
- Handel A, Longini IM Jr, Antia R: Neuraminidase inhibitor resistance in influenza: assessing the danger of its generation and spread. PLoS Comput Biol. 2007, 3: e240-10.1371/journal.pcbi.0030240.View ArticlePubMedPubMed CentralGoogle Scholar
- Heffernan J, Keeling MJ: An in-host model of acute infection: Measles as a case study. Theor Popul Biol. 2008, 73: 134-147. 10.1016/j.tpb.2007.10.003.View ArticlePubMedGoogle Scholar
- Tamura S, Kurata T: Defense mechanisms against Influenza virus infection in the respiratory tract mucosa. Jpn J Infect Dis. 2004, 57: 236-247.PubMedGoogle Scholar
- Handel A, Longini IM Jr, Antia R: Towards a quantitative understanding of the within-host dynamics of influenza A infections. J R Soc Interface. 2010, 7: 35-47. 10.1098/rsif.2009.0067.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/14/1019/prepub
Pre-publication history
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.