Theoretical perspectives on the infectiousness of Ebola virus disease
Theoretical Biology and Medical Modellingvolume 12, Article number: 1 (2015)
Ebola virus disease (EVD) has generated a large epidemic in West Africa since December 2013. This mini-review is aimed to clarify and illustrate different theoretical concepts of infectiousness in order to compare the infectiousness across different communicable diseases including EVD.
We employed a transmission model that rests on the renewal process in order to clarify theoretical concepts on infectiousness, namely the basic reproduction number, R0, which measures the infectiousness per generation of cases, the force of infection (i.e. the hazard rate of infection), the intrinsic growth rate (i.e. infectiousness per unit time) and the per-contact probability of infection (i.e. infectiousness per effective contact).
Whereas R0 of EVD is similar to that of influenza, the growth rate (i.e. the measure of infectiousness per unit time) for EVD was shown to be comparatively lower than that for influenza. Moreover, EVD and influenza differ in mode of transmission whereby the probability of transmission per contact is lower for EVD compared to that of influenza.
The slow spread of EVD associated with the need for physical contact with body fluids supports social distancing measures including contact tracing and case isolation. Descriptions and interpretations of different variables quantifying infectiousness need to be used clearly and objectively in the scientific community and for risk communication.
An epidemic of Ebola virus disease (EVD) centred in three West African countries has been ongoing since December 2013, with limited international spread to other countries in Africa, Europe and the USA . It is likely that the duration of this EVD epidemic, associated with a high case fatality risk (CFR) estimated at ~70% [2, 3], will extend well into 2015. To investigate the ongoing EVD transmission dynamics and consider a range of possible countermeasures, it is vital to understand the natural history and epidemiological dynamics of this disease.
Owing to the rapid progression of the EVD epidemic in West Africa, attempts have been made to clarify the fundamental epidemiological characteristics of EVD [1, 2, 4]. For instance, several studies have reported statistical estimates of the reproduction number, i.e., the average number of secondary cases generated by a single primary case, as a measure of the transmission potential of EVD [2, 5–12]. Despite substantial progress, it remains unclear how measures of infectiousness (or the transmissibility) of EVD should be communicated to the public and interpreted in light of the set of control interventions that could be considered in practical settings. Hence, the purpose of this mini-review is to comprehensively classify different theoretical aspects of infectiousness using a basic transmission model formulated in terms of a renewal process. This approach allows us to compare different measures of infectiousness across different communicable diseases and design possible countermeasures.
Here we briefly review the definition of the basic reproduction number, R0 using the renewal process model . Let i(t) represent the incidence (i.e. the transient number of new cases) at calendar time t. Assuming that the contribution of initial cases to the dynamics is negligible, the renewal process is written as
where A(s) is the rate of secondary transmission per single primary case at its infection-age (i.e., the time since infection) s. Using A(s), one can model the dependency of the transmission dynamics on infection-age . By far the most commonly used measure of infectiousness is the basic reproduction number, R0, which is computed as
and it can be interpreted as the number of secondary cases produced by a single primary case throughout its entire course of infection in a completely susceptible population. Although the concept of R0 is well-known, it is important to note from (2) that R0 results from the integration over all infection-ages. It is well known that the mathematical definition of R0 in a heterogeneously mixing population is described by using the multivariate version of (1) and the next-generation matrix that maps secondary transmissions between and within sub-populations. R0 is defined as the largest eigenvalue of this matrix [15, 16]. Similarly, the definition of R0 can be adapted to the situation of periodic infectious diseases by handling the seasonal dynamics using a vector and employing Floquet theory (see e.g., ).
Although R0 is clearly a dimensionless quantity, the conceptual interpretation from the renewal process (1) permits us to regard R0 as the average number of infected cases produced “per generation”. For this reason, R0 could also be referred to as the basic reproductive ratio, as it could be calculated as the ratio of secondary to primary cases.
where s(t) is the fraction susceptible at time t, β(s) the rate of transmission per single infected individual at infection-age s, and Γ(s) the survivorship of infectiousness at infection-age s. Here we define the force of infection, λ(t) as
which yields a measure of the risk of infection in a susceptible population. The force of infection can be interpreted as the hazard of infection in statistical sense -- the rate at which susceptible individuals are infected . In the classical Kermack and McKendrick epidemic model, λ(t) is modelled as proportional to the disease prevalence . The force of infection is useful for the analysis of incidence data.
Comparison of three communicable diseases
Table 1 shows empirical estimates of R0 and the mean generation time for three different infectious diseases that are characterized by significantly different transmissibility and natural history parameters, i.e., measles, influenza H1N1-2009 and EVD [1, 19, 20]. The mean generation time, Tg can be mathematically derived from the transmission kernel in the renewal process (1), i.e.,
The mode of transmission greatly differs for three diseases considered (Table 1). Measles is transmitted efficiently through the air while the transmission of influenza mostly occurs via droplet although airborne transmission is also possible in a confined setting . In contrast, transmission of EVD is greatly constrained to physical contacts via body fluids . Despite the differences in the mode of transmission for these diseases, it is important to note that the estimates of R0 for H1N1-2009 and EVD are not too different (Table 1). Does that indicate that influenza (H1N1-2009) and Ebola are similarly infectious?
While the average R0 for influenza and Ebola are similar, here we underscore that their underlying transmission dynamics show fundamental differences. This can be understood by analysing the intrinsic growth rate r for both diseases. Assuming that the early growth of each disease follows an exponential form, i.e., i(t) = i0exp(rt) (where i0 is a constant), the renewal equation (1) is rewritten as the so-called Euler-Lotka equation. Replacing i(t) in both sides of (1) by i0exp(rt) and cancelling exp(rt) from both sides, we obtain
yielding the relationship between R0 and the generation time,
where g(s) is the probability density function of the generation time. Equation (7) frequently appears in discussions of mathematical demography  and theoretical epidemiology , which is useful to describe how the relationship is determined between R0 and the intrinsic growth rate r as a function of the generation time distribution. For instance, if the generation time distribution follows the exponential distribution or Kronecker delta function, we obtain the well-known estimators of R0, i.e., R0 = 1 + rTg and R0 = exp(rTg), respectively . Assuming that g(s) follows a gamma distribution with the coefficient of variation k, we have
It should be noted that it is possible that the right-tail of g(s) for EVD might have been underestimated if there were substantial number of secondary transmissions from deceased persons during funerals. Adopting the values of R0 and Tg given in Table 1, and assuming that the coefficient of variation of the generation time at 50%, the intrinsic growth rate of influenza H1N1-2009 is calculated as 0.125 per day, while that of EVD is calculated as 0.038 per day. Figure 1A compares the growth rates (r) of three representative communicable diseases for different values of the coefficient of variation of the generation time. An epidemic of measles appears to grow the fastest followed by one of influenza while an outbreak of EVD is expected to grow the slowest. Whereas the R0 for EVD is similar to that of influenza, the growth rate of EVD is far smaller than that of influenza. This is because each disease generation in the context of EVD transmission takes approximately two weeks, while each generation of new influenza cases occurs on a much shorter time scale - every 3 days on average. Moreover, EVD spreads comparatively slowly mainly by physical contact. This feature indicates that social distancing measures including contact tracing and case isolation could be powerful options for controlling EVD assuming that public health infrastructure exists for these interventions to be feasible .
Thus, based on the infectiousness as measured by the growth rate of cases per unit time, it is very encouraging that EVD is far less dispersible than influenza. Although static countermeasures (e.g. mass vaccination at a certain age) can be planned using R0, the feasibility to deploy dynamic countermeasures, such as contact tracing and case isolation rests on the competition between the growth of cases and public health control, and in this context, the key parameter of infectiousness to assess the feasibility of control interventions is the intrinsic growth rate of cases.
Per contact risk of infection
We further decompose the rate of secondary transmission per single primary case in the renewal equation (3) into the product of the contact rate c(s) and the per-contact probability of infection p(s), i.e.,
Assuming that the per contact probability of infection, p is independent of infection-age, we have
The interpretation of p is straightforward, i.e., it can be regarded as the risk of successful secondary transmission given an infectious contact to a susceptible individual. Assuming that everyone is susceptible at time zero, R0 in (2) is rewritten as
As mentioned above, R0 for EVD is similar to that of influenza. Nevertheless, the infectious period, modelled by Γ(s) for EVD is longer than that of influenza. Assuming an identical contact rate, c, between EVD and influenza, equation (11) indicates that the per-contact probability of infection for EVD is smaller than that for influenza.
The mode of transmission differs across communicable diseases. Figure 2 illustrates the physical range of “contact” that can potentially lead to infection for three representative infectious diseases. Measles causes airborne transmission, and thus, it can lead to secondary infections across different rooms (or sometimes even across buildings). The extent of contact for EVD is very limited as it is highly constrained to physical contacts with body fluids. Hence, effective contact for EVD is limited to close contacts that might be unavoidable among healthcare workers and household members of cases.
We have comparatively discussed concepts of infectiousness for EVD in relation to other communicable diseases from a mathematical modelling point of view. The measure of infectiousness per generation of cases is R0. R0 offers a threshold principle and we have discussed that this measure is important for planning some static countermeasures such as mass vaccination. Based on R0, the overall infectiousness of EVD may be perceived to be similar to that of influenza. Nevertheless, the infectiousness per unit time for EVD was shown to be comparatively lower than influenza. The slow spread of EVD supports social distancing measures including contact tracing and case isolation. Moreover, the per-contact probability of infection for EVD is lower than that for influenza, and the mode of transmission also differs. These findings should also be regarded as encouraging news for healthcare workers who would have to have unavoidable and protected contact with EVD cases. In summary, there is a need for the use of clear and objective descriptions and interpretations of different variables quantifying infectiousness among the scientific community and for risk communication.
- R 0 :
The basic reproduction number
Ebola virus disease
Case fatality risk.
Chowell G, Nishiura H: Transmission dynamics and control of Ebola virus disease (EVD): a review.BMC Med 2014, 12:196.
WHO Ebola Response Team: Ebola virus disease in West Africa–the first 9 months of the epidemic and forward projections.N Engl J Med 2014, 371:1481–1495.
Kucharski AJ, Edmunds WJ: Case fatality rate for Ebola virus disease in west Africa.Lancet 2014, 384:1260.
Incident Management System Ebola Epidemiology Team, CDC; Ministries of Health of Guinea, Sierra Leone, Liberia, Nigeria, and Senegal; Viral Special Pathogens Branch, National Center for Emerging and Zoonotic Infectious Diseases, CDC: Ebola virus disease outbreak - West Africa, September 2014.Morb Mortal Wkly Rep 2014, 63:865–866.
Nishiura H, Chowell G: Early transmission dynamics of Ebola virus disease (EVD), West Africa, March to August 2014.Euro Surveill 2014., 19: Available online: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=20894
Althaus CL: Estimating the reproduction number of Zaire ebolavirus (EBOV) during the 2014 outbreak in West Africa.PLOS Curr Outbreaks 2014. Sep 2. Edition 1. doi:10.1371/currents.outbreaks.91afb5e0f279e7f29e7056095255b288
Fisman D, Khoo E, Tuite A: Early epidemic dynamics of the West African 2014 Ebola outbreak: estimates derived with a simple two-parameter model.PLOS Curr Outbreaks 2014. Sep 8. Edition 1. doi:10.1371/currents.outbreaks.89c0d3783f36958d96ebbae97348d571
Gomes MF, Piontti AP, Rossi L, Chao D, Longini I, Halloran ME, et al.: Assessing the international spreading risk associated with the 2014 West African Ebola outbreak.PLOS Curr Outbreaks 2014. Sep 2. Edition 1. doi:10.1371/currents.outbreaks.cd818f63d40e24aef769dda7df9e0da5
Towers S, Patterson-Lomba O, Castillo-Chavez C: Temporal variations in the effective reproduction number of the 2014 West Africa Ebola outbreak.PLOS Curr Outbreaks 2014. Sep 18. Edition 1. doi:10.1371/currents.outbreaks.9e4c4294ec8ce1adad283172b16bc908
Yamin D, Gertler S, Ndeffo-Mbah ML, Skrip LA, Fallah M, Nyenswah TG, et al.: Effect of Ebola progression on transmission and control in Liberia.Ann Intern Med 2014. doi:10.7326/M14–2255. in press
Fasina F, Shittu A, Lazarus D, Tomori O, Simonsen L, Viboud C, et al.: Transmission dynamics and control of Ebola virus disease outbreak in Nigeria, July to September 2014.Eur Surveill 2014., 19: Available online: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=20920
Althause CL, Gsteiger S, Low N: Ebola virus disease outbreak in Nigeria: lessons to learn.Peer J PrePrints 2014, 2:e569v1.
Diekmann O, Heesterbeek JAP: Mathematical epidemiology of infectious diseases: model building, analysis and interpretation. Chichester: Wiley; 2000.
Nishiura H: Time variations in the generation time of an infectious disease: implications for sampling to appropriately quantify transmission potential.Math Biosci Eng 2010, 7:851–869.
Diekmann O, Heesterbeek JA, Metz JA: On the definition and the computation of the basic reproduction ratioR 0 in models for infectious diseases in heterogeneous populations.J Math Biol 1990, 28:365–382.
Nishiura H, Chowell G, Safan M, Castillo-Chavez C: Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009.Theor Biol Med Model 2010, 7:1. 10.1186/1742-4682-7-1
Bacaër N, Ait Dads el H: On the biological interpretation of a definition for the parameter R 0 in periodic population models.J Math Biol 2012, 65:601–621. 10.1007/s00285-011-0479-4
Farrington CP: Modelling forces of infection for measles, mumps and rubella.Stat Med 1990, 9:953–967. 10.1002/sim.4780090811
Fine PE: Herd immunity: history, theory, practice.Epidemiol Rev 1993, 15:265–302.
Boëlle PY, Ansart S, Cori A, Valleron AJ: Transmission parameters of the A/H1N1 (2009) influenza virus pandemic: a review.Influenza Other Respir Viruses 2011, 5:306–316. 10.1111/j.1750-2659.2011.00234.x
Inaba H, Nishiura H: The state-reproduction number for a multistate class age structured epidemic system and its application to the asymptomatic transmission model.Math Biosci 2008, 216:77–89. 10.1016/j.mbs.2008.08.005
Eichner M, Dowell SF, Firese N: Incubation period of Ebola hemorrhagic virus subtype Zaire.Osong Public Health Res Persptect 2011, 2:3–7. 10.1016/j.phrp.2011.04.001
Cowling BJ, Ip DK, Fang VJ, Suntarattiwong P, Olsen SJ, Levy J, et al.: Aerosol transmission is an important mode of influenza A virus spread.Nat Commun 2013, 4:1935.
Keyfitz BL, Keyfitz N: The McKendrick partial differential equation and its uses in epidemiology and population study.Math Comp Model 1997, 26:1–9.
Wallinga J, Lipsitch M: How generation intervals shape the relationship between growth rates and reproductive numbers.Proc R Soc Lond Ser B 2007, 274:599–604. 10.1098/rspb.2006.3754
Roberts MG, Heesterbeek JA: Model-consistent estimation of the basic reproduction number from the incidence of an emerging infection.J Math Biol 2007, 55:803–816. 10.1007/s00285-007-0112-8
Eichner M, Dietz K: Transmission potential of smallpox: estimates based on detailed data from an outbreak.Am J Epidemiol 2003, 158:110–117. 10.1093/aje/kwg103
HN received funding support from the Japan Science and Technology Agency (JST) CREST program, RISTEX program for Science of Science, Technology and Innovation Policy, and St Luke’s Life Science Institute Research Grant for Clinical Epidemiology Research 2014. GC acknowledges financial support from the NSF grant 1414374 as part of the joint NSF-NIH-USDA Ecology and Evolution of Infectious Diseases program, UK Biotechnology and Biological Sciences Research Council grant BB/M008894/1, and the Division of International Epidemiology and Population Studies, The Fogarty International Center, US National Institutes of Health.
The authors declare that they have no competing interests.
HN conceived of the study. HN conducted mathematical analyses and drafted the manuscript. HN and GC drafted figures and table together and revised the manuscript. All authors read and approved the final manuscript.