Skip to main content

Extracting key information from historical data to quantify the transmission dynamics of smallpox

Abstract

Background

Quantification of the transmission dynamics of smallpox is crucial for optimizing intervention strategies in the event of a bioterrorist attack. This article reviews basic methods and findings in mathematical and statistical studies of smallpox which estimate key transmission parameters from historical data.

Main findings

First, critically important aspects in extracting key information from historical data are briefly summarized. We mention different sources of heterogeneity and potential pitfalls in utilizing historical records. Second, we discuss how smallpox spreads in the absence of interventions and how the optimal timing of quarantine and isolation measures can be determined. Case studies demonstrate the following. (1) The upper confidence limit of the 99th percentile of the incubation period is 22.2 days, suggesting that quarantine should last 23 days. (2) The highest frequency (61.8%) of secondary transmissions occurs 3–5 days after onset of fever so that infected individuals should be isolated before the appearance of rash. (3) The U-shaped age-specific case fatality implies a vulnerability of infants and elderly among non-immune individuals. Estimates of the transmission potential are subsequently reviewed, followed by an assessment of vaccination effects and of the expected effectiveness of interventions.

Conclusion

Current debates on bio-terrorism preparedness indicate that public health decision making must account for the complex interplay and balance between vaccination strategies and other public health measures (e.g. case isolation and contact tracing) taking into account the frequency of adverse events to vaccination. In this review, we summarize what has already been clarified and point out needs to analyze previous smallpox outbreaks systematically.

Background

Smallpox epidemiology has the longest and richest history in investigating the mechanisms of spread and in evaluating the effectiveness of vaccination [1, 2]. Modern epidemiological methods have developed in parallel with smallpox control practice, and consequently, the disease had already been eradicated before statistical and epidemiological techniques for analyzing infectious disease outbreaks had sufficiently matured.

Although the world is free from smallpox, researchers continue revisiting smallpox epidemiology and virology with recent techniques. In the aftermath of the 9-11-2001 attack, the awareness of the threat of bioterrorism has grown significantly [3]. Mathematical models and computer simulations have been developed to design and optimize public health measures against re-introduced variola virus, the causative agent of smallpox [417]. These models are based on different epidemiological assumptions of smallpox. For example, assumptions about the number of secondary transmissions before onset of illness had not been not carefully validated in earlier mathematical modelling studies [16, 17]. Accordingly, the policy implications of these models differed widely, and thus the necessity arose to capture the basic mechanisms of smallpox transmission precisely [6, 18]. To date, it has been demonstrated that transmission dynamics and intervention strategies cannot be modelled without sufficiently quantifying the detailed intrinsic mechanisms by using observed data [6, 19, 20]. Because of the global eradication, we have had to maximize the use of historical data to estimate nearly all biological and epidemiological parameters that are needed to optimize interventions [21].

This review article has two purposes. The first purpose is to summarize the issues that have been clarified in recent mathematical and statistical studies and to discuss the relevant policy implications. The second purpose is to specify what important aspects of smallpox epidemiology remain unknown and to suggest how these could be addressed by analyzing historical records. In the following section, we first give a technical overview of the use of historical data and then present some examples of quantification. Subsequently, we summarize the basic concept and interpretation of the transmission potential and the resulting implications for vaccination strategies. The paper concludes with a summary of the findings, emphasizing the importance of systematically analyzing historical datasets.

Review

Extracting key information from historical data

Although historical data have frequently been revisited using modern statistical techniques to identify epidemiological determinants of smallpox, many key issues remain unknown in spite of great efforts. To clarify important aspects of smallpox epidemiology, it remains necessary to maximize the use of historical data. To understand their usefulness and to avoid common pitfalls, we briefly discuss technical issues in utilizing historical publications.

Issues to consider when looking at historical smallpox data

In the following, we list key points to be remembered whenever we statistically extract information from the historical literature. As we may not be able to find all the answers to the following questions in a single historical data set, we may have to combine different data sets or to merge in information from laboratory experiments.

(A) Were all cases caused by variola virus?

As cases could not be confirmed virologically before the middle of the 20th century, it is crucial to know on what observations historical diagnoses were based. It was not uncommon to misdiagnose chickenpox cases as smallpox [22, 23]. In the older literature, it sometimes even remains unclear which kind of "plague" was being described [24, 25]. Ascertainment of diagnostic methods is one of the biggest challenges in utilizing historical outbreak data.

(B) Clinical documentations and time-varying medical trends

Similarly, clinical classifications of smallpox have been revised over time [1, 2628]. The definition of severe smallpox has varied greatly even in the 20th century. Vaccines have continuously been improved [29], and we still do not even know from where the vaccinia virus emerged and when it started to be used as a smallpox vaccine [30]. It is necessary to identify and to select the most useful sources of literature, in order to understand which classification in a given publication was adopted and which type of vaccine was most likely to have been used in the population described.

(C) Pathogenicity and virulence of the variola virus

Classically, smallpox was classified into two different types according to the observed case fatality. The traditional form of smallpox, referred to as variola major, was believed to have a case fatality of 20% or more. A milder form of the disease, variola minor, with a case fatality of 1% or less was first reported in the late 19th century in South Africa, then it was observed in European countries and finally in Brazil [1, 3133]. Variola minor accounted for the majority of cases in the early 20th century in the United States, where it remained the only form of smallpox from the 1930s until its eradication [34]. The epidemiology of variola minor and its interplay with variola major have only partly been clarified [35, 36]. There are clear genetic differences between variola major and minor, supporting the taxonomic distinction; recently, the virulence of variola virus has also been attributed to detailed genomic information [37, 38]. However, if case fatality was a major criterion in determining the virulence of variola virus, the outcome of these laboratory studies may have been distorted by the vaccination history of cases and maybe also by other factors. Epidemiological clarification of this point still remains an open question.

(D) Definition of the reported events

When extracting information on the incubation and infectious periods (or similar parameters describing the epidemiological characteristics), it is crucial to know how the time of infection (which cannot be observed directly) and the onset of disease were defined. There were two traditional ways to define the onset of smallpox: onset of fever or appearance of rash. If the period from onset to recovery is documented, it is important furthermore to identify what "recovery" stands for (e.g. recovery from pyrexia or solidification/disappearance of rash).

Extracting data from historical publications

The foregoing list does not cover all common pitfalls. Tackling historical data requires not only statistical techniques, but also understanding of the social history and the background of the cases. Moreover, as noted above, we often have to draw conclusions with implications for public health decision-making using combined data from different sources. Identifying the most useful and important data and addressing key questions are major parts of an essential process to shed light on the mechanisms of transmission and spread of smallpox. In the next two sections, we review studies on parameter estimation that use historical records and predominantly originate from our previous studies. The following case studies were conducted, carefully accounting for common technical problems as listed above when looking at the historical data.

Intrinsic transmission process of smallpox

To understand the spread of smallpox, it is essential to know the intrinsic transmission process, i.e. after what time symptoms appear, when secondary transmission occurs, and how severe the disease will be. Although basic, such knowledge of the intrinsic transmission process already allows us to assess whether public health interventions in the event of a bioterrorist attack can contain smallpox by means of mass vaccination or by a combination of contact tracing, quarantine and isolation [6, 19]. As practical examples, here we briefly discuss basic methodologies and recent findings concerning the incubation period, the infectious period and the case fatality.

Incubation period of smallpox

The incubation period is defined as the time from infection to onset of disease [39]. Usually, symptoms of smallpox appear 10–14 days after infection [40]. The knowledge of the incubation period distribution enables us to determine the appropriate length of quarantine [41]. 'Quarantine' here refers to physical separation of healthy individuals who were exposed to cases. In the practice of outbreak investigations, the time of exposure is sometimes determined by contact tracing. Historically, the suggested length of smallpox quarantine tended to be 14–16 days, based on professional experience and an accumulation of epidemiological data, but not on an explicit statistical analysis of the incubation period distribution.

Restricting the movement of exposed individuals for longer than the maximum incubation period ensures the effectiveness of quarantine measures. Unfortunately, the length of the incubation period requires knowing the exact time of infection, and thus can only be determined for cases who were exposed for a very short period of time. In addition, the maximum observed incubation period clearly depends on the sample size: the larger the sample size, the more likely we are to find cases whose incubation periods exceed the previously known maximum. The number of smallpox cases with well-known incubation period (e.g. documented patients who had been exposed for a single day only) is limited in historical records. The problem of stating a maximum incubation period can be circumvented by fitting a statistical distribution to the data. This distribution allows a time point to be determined beyond which the onset of further cases becomes extremely unlikely (e.g. the time after infection until which 99% of the patients develop symptoms). If the incubation period follows a lognormal distribution with mean, μ, and standard deviation, σ (of the variable's logarithm), the probability density of observing an incubation period of length ti is given by

f ( t i ; μ , σ 2 ) = 1 t i σ 2 π exp ( ( ln ( t i ) μ ) 2 σ 2 ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOzayMaeiikaGIaemiDaq3aaSbaaSqaaiabdMgaPbqabaGccqGG7aWocqaH8oqBcqGGSaalcqaHdpWCdaahaaWcbeqaaiabikdaYaaakiabcMcaPiabg2da9KqbaoaalaaabaGaeGymaedabaGaemiDaq3aaSbaaeaacqWGPbqAaeqaaiabeo8aZnaakaaabaGaeGOmaiJaeqiWdahabeaaaaGccyGGLbqzcqGG4baEcqGGWbaCdaqadaqaaiabgkHiTKqbaoaalaaabaGaeyOeI0YaaeWaaeaacyGGSbaBcqGGUbGBcqGGOaakcqWG0baDdaWgaaqaaiabdMgaPbqabaGaeiykaKIaeyOeI0IaeqiVd0gacaGLOaGaayzkaaWaaWbaaeqabaGaeGOmaidaaaqaaiabeo8aZnaaCaaabeqaaiabikdaYaaaaaaakiaawIcacaGLPaaaaaa@5A74@
(1)

We can estimate the parameters μ and σ from a dataset of n known incubation times ti by maximizing the likelihood function

L ( μ , σ 2 ) = i = 1 n f ( t i ; μ , σ 2 ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemitaWKaeiikaGIaeqiVd0MaeiilaWIaeq4Wdm3aaWbaaSqabeaacqaIYaGmaaGccqGGPaqkcqGH9aqpdaqeWbqaaiabdAgaMjabcIcaOiabdsha0naaBaaaleaacqWGPbqAaeqaaOGaei4oaSJaeqiVd0MaeiilaWIaeq4Wdm3aaWbaaSqabeaacqaIYaGmaaGccqGGPaqkaSqaaiabdMgaPjabg2da9iabigdaXaqaaiabd6gaUbqdcqGHpis1aaaa@48EE@
(2)

Figure 1 shows the frequency distribution of the incubation period of smallpox, which was estimated from 131 cases of smallpox who were exposed only for a single day [42]. The mean incubation period is 12.5 days (SD 2.2 days). The 99th percentile is 18.6 days with a 95% confidence interval (CI) ranging from 16.8 to 22.2 days. This indicates that a quarantine of 23 days ensures that more than 99% of infected individuals will develop symptoms before being released.

Figure 1
figure 1

Incubation period distribution of smallpox fitted to a lognormal distribution (n = 131). The vertical arrow indicates the maximum likelihood estimate of the 99th percentile of the incubation period [42]. The median and the coefficient of variation are 12.5 days and 18.0%, respectively.

Infectious period of smallpox

The infectious period has traditionally been defined as the period in which pathogens are discharged [43]. It presently refers to the period in which infected individuals are capable of generating secondary cases. Knowledge of the infectious period allows us to determine for how long known cases need to be isolated and what should be the latest time point after exposure at which newly infected individuals should be in isolation.

One approach to addressing this issue is to quantify how the pathogen load changes over time using the most sensitive microbiological techniques (e.g. polymerase chain reaction), but such observations are usually limited to the period after onset of symptoms. Several attempts have been made to measure the distribution of the virus-positive period of smallpox cases [44, 45], but sample sizes were small and only very few samples could be obtained during the early stage of illness. Moreover, linking "virus-positive" results to the probability of causing secondary transmission is difficult without further information (e.g. frequency, mode and degree of contact).

Another way of addressing this complicated issue is to determine the frequency of secondary transmission relative to disease-age, i.e. the time since onset of fever [46]. An estimate of the relative infectiousness is obtained by analyzing historical data in which it is known who acquired infection from whom. The known transmission network permits serial intervals to be extracted, i.e. the times from symptom onset in a primary case to symptom onset in a secondary case [4749]. Given the length of the serial interval s and the corresponding length of the incubation period f, the disease-age l from onset of symptoms to secondary transmission satisfies

Considering the statistical distributions for each length results in a convolution equation:

s ( t ) = 0 t l ( t τ ) f ( τ ) d τ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4CamNaeiikaGIaemiDaqNaeiykaKIaeyypa0Zaa8qmaeaacqWGSbaBcqGGOaakcqWG0baDcqGHsislcqaHepaDcqGGPaqkcqWGMbGzcqGGOaakcqaHepaDcqGGPaqkcqWGKbazcqaHepaDaSqaaiabicdaWaqaaiabdsha0bqdcqGHRiI8aaaa@455B@
(4)

The frequency l(t-τ) of secondary transmission relative to disease-age can be back-calculated by extracting the serial interval distribution s(t) from a known transmission network, and by using the incubation period distribution f(τ) given above. If we have information on the length ti of the serial interval for n cases, the likelihood function is given by

L = i = 1 n s ( t i ) = i = 1 n 0 t i l ( t i τ ) f ( τ ) d τ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemitaWKaeyypa0ZaaebCaeaacqWGZbWCcqGGOaakcqWG0baDdaWgaaWcbaGaemyAaKgabeaakiabcMcaPaWcbaGaemyAaKMaeyypa0JaeGymaedabaGaemOBa4ganiabg+GivdGccqGH9aqpdaqeWbqaamaapedabaGaemiBaWMaeiikaGIaemiDaq3aaSbaaSqaaiabdMgaPbqabaGccqGHsislcqaHepaDcqGGPaqkcqWGMbGzcqGGOaakcqaHepaDcqGGPaqkcqWGKbazcqaHepaDaSqaaiabicdaWaqaaGGabiab=bcaGiabdsha0naaBaaameaacqWGPbqAaeqaaaqdcqGHRiI8aaWcbaGaemyAaKMaeyypa0JaeGymaedabaGaemOBa4ganiabg+Givdaaaa@5AD4@
(5)

The parameters that describe the frequency of secondary transmission relative to disease-age can be estimated by maximizing this function. A similar method has been applied to estimate the number of HIV-infections from AIDS incidence [50].

Figure 2 shows the back-calculated infectiousness of smallpox relative to disease-age [46]. In the following text, day 0 represents the onset of fever. Before onset of fever (i.e. between day -5 and day -1) altogether only 2.7% of all transmissions occurred. Between day 0 and day 2 (i.e. in the prodromal period before the onset of rash) a total of 21.1% of all transmissions occurred. The daily frequency of passing on the infection was highest between day 3 and day 5, yielding a total of 61.8% of all transmissions. These estimates help determine the latest time by which cases should be in isolation. If each primary case infects on average six individuals, and if the effectiveness of isolation is 100%, the isolation of a primary case before the onset of rash reduces the expected number of victims to 6 × (0.027 + 0.211) = 1.428. Optimal isolation could, therefore, substantially reduce the number of secondary cases, and the outbreak could quickly be brought under control by additional countermeasures (e.g. contact tracing [6, 51]).

Figure 2
figure 2

Relative frequency of secondary transmissions of smallpox by disease-age. Expected daily frequency of secondary transmissions with corresponding 95% confidence intervals [46]. The percentages indicate the fraction of transmissions among all transmissions that occurred in the given intervals. The disease-age t = 0 denotes the onset of fever.

Case fatality

Case fatality is the proportion of deaths among those developing the disease. It is particularly important to understand the case fatality of smallpox in order to estimate the magnitude of the disaster in the event of a bioterrorist attack. It may also be of practical importance to predict the burden of hospital admissions and burials in such an event. The case fatality of smallpox was systematically reviewed during the Eradication Programme [52], showing extremely high crude estimates of 26% and 36% in East-Pakistan and Madras, respectively, but suggesting a wide geographical heterogeneity. Recent studies attributed part of this heterogeneity to viral genomic differences [37, 38], but many of the previously published mathematical models simply assumed an overall estimate of 30% for unvaccinated cases.

Various factors influence case fatality, most importantly previous vaccination history (which will be discussed in the Section on public health interventions) and the age at infection. Following a previous study by Dietz and Heesterbeek [53], we assume the following parametric model for the age-specific case fatality of smallpox:

where α, β, γ and δ are parameters that need to be estimated. If we have a dataset with Mi deaths and Ni survivors of age ai, the likelihood function is

L ( α , β , γ , δ ) = i c ( a i ) M i ( 1 c ( a i ) ) N i MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemitaWKaeiikaGIaeqySdeMaeiilaWIaeqOSdiMaeiilaWIaeq4SdCMaeiilaWIaeqiTdqMaeiykaKIaeyypa0ZaaebuaeaacqWGJbWycqGGOaakcqWGHbqydaWgaaWcbaGaemyAaKgabeaakiabcMcaPmaaCaaaleqabaGaemyta00aaSbaaWqaaiabdMgaPbqabaaaaOGaeiikaGIaeGymaeJaeyOeI0Iaem4yamMaeiikaGIaemyyae2aaSbaaSqaaiabdMgaPbqabaGccqGGPaqkcqGGPaqkdaahaaWcbeqaaiabd6eaonaaBaaameaacqWGPbqAaeqaaaaaaSqaaiabdMgaPbqab0Gaey4dIunaaaa@519A@
(7)

where ai is a mid-point of age group i. Figure 3 shows age-specific case fatality estimates of unvaccinated cases in Verona, Italy, from 1810–38 and Sheffield, UK, from 1887–88, respectively [54, 55]. The age-specific case fatality of smallpox can be depicted as a U-shaped curve that peaks in infancy and in old age. Smallpox case fatality also depends on other biological factors of the host such as pregnancy [56], which increases the case fatality from 12.7% (estimate for non-pregnant healthy adults; 95% CI: 11.2–14.3) to 34.3% (95% CI: 31.4–37.1) [57]. Underlying diseases (e.g. cancer, diabetes mellitus, HIV infection and medical immunosuppression for transplantation) could further increase the case fatality.

Figure 3
figure 3

Age-specific case fatality of smallpox. Observed (grey bars) and fitted (continuous line) age-specific case fatalities of unvaccinated cases in (A) Verona, Italy, 1810–1838 [53, 55] and (B) Sheffield, UK, 1887–8 [54].

Above, we have presented the three most important components of the intrinsic transmission process. Each of them plays a key role in determining the optimal intervention strategy. We showed some basic applications of utilizing likelihood functions [58], but various other statistical approaches have been taken which were motivated by similar epidemiological interests. These include the applications of non-linear models [59] and of Bayesian techniques [60].

Transmission potential

In addition to the epidemiological parameters that characterize the natural history of smallpox, we have to know the most important summary measure of transmission, the basic reproduction number, R0, in order to design and optimize public health interventions. R0 is defined as the average number of secondary cases arising from a single index case in a fully susceptible population in the absence of interventions [61, 62]. Here, we discuss the concept of R0 and its estimation, starting with its historical development. Then we use the basic reproduction number to assess the eradication threshold of smallpox by means of mass vaccination.

R0 and vaccination

Smallpox is the disease with the longest history in theoretical modelling. During the 18th century, the famous mathematician Daniel Bernoulli modelled the spread of smallpox and assessed the effectiveness of the variolation practice (variolation was the precursor of vaccination, consisting of the inoculation of variola virus) [53, 63]. Moreover, the earliest formulation and calculation of R0 may have been stimulated by smallpox [64]. The earliest concept of R0 and the relevant insights into the effectiveness of smallpox vaccination are revisited in the following.

Figure 4a shows the result of the simple mathematical model developed by Theophil Lotz in the late 19th century [64]. If a single primary case generates on average R0 = 2 secondary cases, and if we ignore for the sake of simplicity the depletion of susceptible individuals, the number of cases grows geometrically. If there are a index cases in generation 0, the expected numbers of cases in generations 1, 2, 3, ..., n will be

Figure 4
figure 4

Theoretical initial courses of smallpox outbreaks following a geometric growth. A. The infection tree (i.e. transmission network) of smallpox is shown by generation, following equation (8). For simplicity, R0 is assumed to be 2. B. Infection tree under vaccination. Vaccination is assumed to reduce the number of secondary transmission by 50%, and thus only 1 case in each generation is expected.

a R 0 , a R 0 2 , ... , a R 0 n , MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaemOuai1aaSbaaSqaaiabicdaWaqabaGccqGGSaaliiqacqWFGaaicqWGHbqycqWGsbGudaqhaaWcbaGaeGimaadabaGaeGOmaidaaOGaeiilaWIae8hiaaIae8Nla4Iae8Nla4Iae8Nla4IaeiilaWIae8hiaaIaemyyaeMaemOuai1aa0baaSqaaiabicdaWaqaaiabd6gaUbaakiabcYcaSiab=bcaGaaa@4292@
(8)

respectively. Following Lotz's example of R0 = 2, and assuming a single index case (a = 1), we expect 2, 4, 8, ..., 2n cases in the subsequent generations. Although the model ignores variations in the number of secondary transmissions (which are deemed particularly important for directly transmitted diseases [65]), the process described reasonably captures the essential dynamics of transmission during the early stages of an epidemic.

We now move on to describe various attempts to estimate R0, summarized in Table 1 together with the underlying key assumptions. Whereas an analysis of long-term temporal dynamics using a mathematical model suggested widely varying estimates of R0, ranging from 3.5 to 6.0 [66], stochastic models assuming a homogeneous pattern of spread, but ignoring the pre-existing immunity level in the afflicted population, grossly underestimated R0 as slightly above unity [67]. A revised estimate by a model that accounts for the detailed intrinsic dynamics of smallpox in an initially partially immune population suggested that R0 is in the order of 6.9 (95% CI: 4.5, 10.1) [18]. This roughly corresponds to an R0 for which 80–90% of vaccination coverage would allow sufficient herd immunity to be achieved [68, 69] (i.e., population-based protection of unvaccinated individuals due to the presence of vaccinated individuals), similarly to Bernoulli's early model, which yielded an estimate of the force of infection that can be translated to R0 = 6.7 [53].

Table 1 Reported estimates of R0 for smallpox

R0 plays a key role in determining the critical vaccination coverage in a randomly mixing population [70]. If v = 80% of individuals are protected by vaccination, the average number of secondary cases is reduced to 20%. Following the model of Lotz, the number of cases in each generation is

a ( 1 v ) R 0 , a ( 1 v ) 2 R 0 2 , ... , a ( 1 v ) n R 0 n MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaeiikaGIaeGymaeJaeyOeI0IaemODayNaeiykaKIaemOuai1aaSbaaSqaaiabicdaWaqabaGccqGGSaalcqWGHbqycqGGOaakcqaIXaqmcqGHsislcqWG2bGDcqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabdkfasnaaDaaaleaacqaIWaamaeaacqaIYaGmaaGccqGGSaalcqGGUaGlcqGGUaGlcqGGUaGlcqGGSaalcqWGHbqycqGGOaakcqaIXaqmcqGHsislcqWG2bGDcqGGPaqkdaahaaWcbeqaaiabd6gaUbaakiabdkfasnaaDaaaleaacqaIWaamaeaacqWGUbGBaaaaaa@5063@
(9)

Figure 4b shows the growth of cases when v = 50% are protected by vaccination: only a single case is expected in each generation (for R0 = 2). The number of cases decreases from one generation to the next if (1-v)R0 is less than 1 (cf. equation (9)). In line with this, we can formulate the most fundamental condition of immunization to achieve a sufficiently high herd immunity level. In order to eradicate an infectious disease by vaccination, the fraction protected by vaccination must satisfy [71]

v > 1 1 R 0 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemODayNaeyOpa4JaeGymaeJaeyOeI0scfa4aaSaaaeaacqaIXaqmaeaacqWGsbGudaWgaaqaaiabicdaWaqabaaaaaaa@3447@
(10)

If R0 is 6 for smallpox, more than 1-1/6 = 83.3% of susceptible individuals need to be successfully immunized to prevent an epidemic by vaccination alone. Although the pattern of smallpox spread is most probably non-random, equation (10) can be used as an approximation to guide policymaking. If all individuals are vaccinated, v can be interpreted as the direct effectiveness of vaccination [72, 73]. The effectiveness of smallpox vaccination remained controversial during the early 20th century, partly because of a lack of reliable estimation methods [2], but the methodologies have greatly improved since then [7478]. During the late 19th century, when vaccine quality was not always ensured, the crude overall effectiveness of vaccination seems to have been higher than 85% [1, 21].

Heterogeneity and behavioural change in relation to historical data

To predict the course and size of an epidemic appropriately, it is critically important to clarify the heterogeneity of transmission. The above-mentioned critical coverage for eradication assumes a randomly mixing population, but it has been established that smallpox spreads for example more easily within households than in the community [7982]. A theoretical approach to modelling household and community transmission separately has been described [14, 83], and a tool that allows the two different levels of transmission to be estimated has been developed [84], but it is very difficult to obtain the necessary estimates from the limited information given in historical records (e.g. detailed household transmission data are always distorted by vaccination). Age-related heterogeneity is yet another important determinant of smallpox epidemiology [85], and spatial patterns of transmission can also influence the success of interventions [86]. Unfortunately, historical records, especially those recorded during the Intensified Smallpox Eradication Programme, are considerably biased (e.g. by individual vaccination histories), and thus it is difficult to address age-related and spatial heterogeneities.

Behavioural changes during an outbreak also have to be clarified to model a bioterrorist attack realistically. It has been suggested that the frequency of contact decreases after the information on an ongoing epidemic is widely disseminated [8789]. A mathematical model that attempted to incorporate such a declining contact frequency during an epidemic suggested that even gradual and moderate behavioural changes could drastically slow the epidemic [90]. Methods incorporating such changes remain yet to be developed to help public health policy making. A generalized method could perhaps incorporate results of a psychological response survey [91].

Public health interventions

Given the basic parameters that describe the intrinsic transmission process, we are now able to examine the effectiveness of interventions. In addition to the critical level of mass vaccination that was discussed in the previous section, here we discuss further issues on vaccination strategies and other public health interventions in bioterrorism preparedness.

Duration of vaccine-induced immunity and partial protection

The degree of protection of vaccinated individuals in the present population is yet another important public health issue. Immunological studies showed that a fraction of previously vaccinated individuals still reacts to exposure with variola virus [92, 93], but it is difficult to attribute each kind of immunological response to actual protection against the disease and its severity. Thus, the degree of protection of individuals who were vaccinated 30 to 50 years ago has remained an open question.

As we have previously shown, an epidemiological model that partly addressed the effect of booster events estimated that primary vaccination protected for a median duration of 11.7 to 28.4 years against the disease [94], indicating that most vaccinated individuals in the present community may no longer be protected from contracting smallpox. However, similar estimates also indicate that vaccinated individuals are still protected against severe manifestations and death from smallpox [95]. An analysis of a statistical record of an outbreak in Liverpool from 1902–3 revealed a median duration of protection against smallpox death of 49.2 years (95% CI: 42.0–57.3) [95, 96]. This finding (of long-lasting partial protection from severe manifestations) was further supported by statistical analyses of similar historical datasets [94] and of individual case records extracted from historical outbreaks in Australia where booster events were extremely rare [97]. In the event of a bioterrorist attack in the early 21st century, residual immunity could significantly decrease the individual burden of disease. However, the persistence of partial protection does not necessarily imply a positive impact on the population level. Masked symptoms may cause difficulties in case recognition and clinical diagnosis. Although it might be virologically plausible that previously vaccinated cases are less infectious (e.g. due to low levels of virus in their nasopharynx), reduced severity may also permit movements of infectious individuals, worsening the prospects of public health control. The ripple benefit of residual immunity has yet to be clarified.

To understand the complex interplay of all partial effects of vaccination, various biological and social effects must be considered. In theory, vaccination does not only diminish the susceptibility of vaccinated individuals, but also reduces the degree and duration of infectiousness upon infection. The vaccine-induced reduction of infectiousness can be estimated using the household secondary attack rate (SAR), expressed as the ratio of the number of infected household contacts to the number of exposed household contacts [98]. Suppose that SARij denotes the household secondary attack rate where i and j, respectively, give the previous vaccination histories of the secondary and primary case (i.e. i or j = 1 represents previously vaccinated, whereas i or j = 0 represents unvaccinated individuals). Let us consider the following household transmission data, which were observed in India [79]:

The household SAR caused by unvaccinated cases among unvaccinated and vaccinated contacts were estimated to be SAR00 = 40/650 = 0.0615 and SAR10 = 11/583 = 0.0189, respectively. Those caused by vaccinated cases among unvaccinated and vaccinated household contacts were SAR01 = 10/499 = 0.0200 and SAR11 = 2/421 = 0.0048, respectively.

The crude efficacy of vaccine in reducing susceptibility VES, infectiousness VEI, and a combined effect of both VET is then estimated by

VE S = 1 SAR 10 SAR 00 = 0.693 VE I = 1 SAR 01 SAR 00 = 0.674 VE T = 1 SAR 11 SAR 00 = 0.923 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabmqaaaqaaiabbAfawjabbweafnaaBaaaleaacqWGtbWuaeqaaOGaeyypa0JaeGymaeJaeyOeI0scfa4aaSaaaeaacqqGtbWucqqGbbqqcqqGsbGudaWgaaqaaiabigdaXiabicdaWaqabaaabaGaee4uamLaeeyqaeKaeeOuai1aaSbaaeaacqaIWaamcqaIWaamaeqaaaaakiabg2da9iabicdaWiabc6caUiabiAda2iabiMda5iabiodaZaqaaiabbAfawjabbweafnaaBaaaleaacqWGjbqsaeqaaOGaeyypa0JaeGymaeJaeyOeI0scfa4aaSaaaeaacqqGtbWucqqGbbqqcqqGsbGudaWgaaqaaiabicdaWiabigdaXaqabaaabaGaee4uamLaeeyqaeKaeeOuai1aaSbaaeaacqaIWaamcqaIWaamaeqaaaaakiabg2da9iabicdaWiabc6caUiabiAda2iabiEda3iabisda0aqaaiabbAfawjabbweafnaaBaaaleaacqWGubavaeqaaOGaeyypa0JaeGymaeJaeyOeI0scfa4aaSaaaeaacqqGtbWucqqGbbqqcqqGsbGudaWgaaqaaiabigdaXiabigdaXaqabaaabaGaee4uamLaeeyqaeKaeeOuai1aaSbaaeaacqaIWaamcqaIWaamaeqaaaaakiabg2da9iabicdaWiabc6caUiabiMda5iabikdaYiabiodaZaaaaaa@7340@
(11)

If we make the simplifying assumption that the biological effect of vaccination was identical for all vaccinated individuals, vaccination reduced susceptibility by 69.3%, infectiousness by 67.4%, and the combined effect was 92.3%. Although an effect of vaccination on the duration of disease has rarely been observed and reported, historical epidemiological studies in Dalian, China, suggested that the mean symptomatic period was reduced by 13.7–48.5% if the case was previously vaccinated [21, 99].

Vaccination strategies

Given that the intrinsic dynamics as well as the effects of vaccination are sufficiently quantified, vaccination strategies against smallpox can be optimized. Three issues, of which the epidemiology has been discussed though the quantitative effect has not yet been fully clarified, are discussed in the following: revaccination, ring vaccination, and post-exposure vaccination.

After it became clear in the late 19th century that vaccine efficacy was not perfect and that vaccine-induced immunity waned over time, revaccination was put into practice. Revaccinated individuals were said to have contracted smallpox less often and had much milder manifestations, so that scheduled revaccinations became accepted in the early to mid 20th century [26], but the intervals from primary vaccination to revaccinations and the number of revaccinations varied widely within and between countries, making analytic evaluations very difficult. Accordingly, it is extremely difficult to quantify the effectiveness of revaccination in reducing the chance of smallpox even with statistical techniques in the present day. Crude estimates of the increased protection against smallpox death were obtained for several outbreaks; e.g. for Madras during the 1960s [27], where 87.1% fewer cases died in the revaccinated group than in the group who had only received the primary vaccination (770/3266 vs. 4/132 deaths, respectively), but this revaccination effect only measures what happened to people who were infected in spite of vaccination. (What makes an explicit interpretation of these findings even more difficult was the fact that vaccination in India was made using the rotary lancet, which left a scar even in the absence of "take".)

Vaccination can be combined with the practice of case finding: Ring vaccination is a surveillance containment measure that consists of vaccinating and monitoring all susceptible individuals in a prescribed area around one or several index cases [100]. This combined strategy is deemed more effective than mass vaccination [101], but combinations of vaccination and public health measures have not yet been explicitly evaluated. Ring vaccination was introduced and evaluated mainly in West and Central Africa and in Asia where it was always combined with case isolation [102]. Although it is difficult to exclude the impact of other interventions and to estimate the net effectiveness of ring vaccination explicitly (e.g. the impact of previous vaccinations can usually not be separated [103]), accumulated experience during the Intensified Eradication Programme strongly suggests that ring vaccination (accompanied by vigorous isolation) worked well [101]. The strategy is deemed logically effective in containing localized outbreaks, but it is important to ensure effective contact tracing if we are to rely on ring vaccination alone [9, 104].

Vaccination may still be protective if a person has already been exposed to the virus, a procedure referred to as post-exposure vaccination [105, 106]. Despite numerous discussions [107], the protective effect of post-exposure vaccination has remained unclear. A historical study from the early 20th century suggests that vaccination within 7 days after exposure is effective [28]. Smallpox textbooks in the 1960s and '70s claimed that 'vaccination within 72 hours almost promises protection' [26, 27], a statement roughly consistent with a more recent statistical estimate based on historical data and on several assumptions concerning the hypothetical frequencies of vaccinated and protected individuals [108], and with a laboratory study demonstrating a cell-mediated response within 4 days after exposure [109]. A similar estimate was obtained in a Delphi analysis, which concluded that post-exposure vaccination can be assumed to be 80–93% effective during the first 3 days after exposure and 2–25% thereafter [110]. However, as we have shown, a statistical exercise demonstrates that historical data, which only record cases who developed smallpox after post-exposure vaccination, hardly provide sufficient insight into the effectiveness of post-exposure vaccination [111]. Information regarding the denominator is insufficient for the majority of records (i.e. we do not know how many exposed people who were vaccinated were protected from the disease). Only the effectiveness of vaccination against severe disease upon infection can be estimated from such data: the shorter the interval between exposure and vaccination, the lower the probability of developing severe smallpox. To the best of our knowledge, only one outbreak in Leicester, UK, from 1903–04 provided insight into the protection against disease by post-exposure vaccination [112]: counting from the eruption of the index case, it was reported that none of 210 individuals (0%) who were vaccinated on the first day after exposure, 2 among 359 (0.5%) who were vaccinated on the second day, 5 among 102 (4.9%) who were vaccinated on the third day, and 10 among 116 (8.6%) who were vaccinated on the fourth day or later developed the disease. Although this seems to indicate some degree of protection, the actual efficacy of post-exposure vaccination can only be determined by comparing these findings to observations in a group of individuals who were exposed for exactly the same periods of time, but refused or were denied post-exposure vaccination.

Despite effective vaccination, pros and cons of vaccination practice always need to account for adverse events of vaccination [113]. Vaccine-strain dependent differences in the frequency of adverse events have been reported, and the risk of death due to vaccination has been analyzed in detail only recently [114, 115]. Theoretical frameworks reported to date agree with each other that we should not implement pre-attack mass vaccination in order to minimize the number of adverse events. Policy suggestions of mathematical models for post-attack vaccination strategies depend on the specific attack scenarios and need to be investigated further.

Case isolation and contact tracing

Rather than relying completely on vaccination, recent modelling studies have suggested that an outbreak could be contained by a combination of case isolation and contact tracing [6, 14], owing mainly to the characteristics of the intrinsic dynamics of smallpox (e.g. the relatively long generation time and the clear symptoms of smallpox). The importance of monitoring and controlling "contacts" has been highlighted in a historical observation [112] and was also stressed during the Eradication Programme [51, 116, 117]. A public health system's capability in conducting contact tracing may determine whether or not a smallpox outbreak can be controlled without vaccination. This should also take into account response logistics and the limited number of public health practitioners [104]. A mathematical exercise suggested that the optimal intervention also depends on the initial attack size: whereas an outbreak caused by few cases could easily be controlled by isolation and contact tracing alone, regional (targeted) mass vaccination is recommended if the initial attack size is big and R0 is large [118].

Conclusion

This article has reviewed quantifications of the transmission and spread of smallpox using historical data. Although historical data are limited and we cannot answer all questions regarding smallpox epidemiology, many publications are available from previous efforts. However, a systematic listing of surveillance data and/or outbreak reports irrespective of language (e.g. see [57]) still remains a future task. It is essential that historians, smallpox specialists and epidemiologists interact more.

Since the eradication, smallpox deaths have disappeared from the world [119], and hope has arisen that we will succeed in eradicating other infectious diseases. Owing to the conceived threat of bioterrorism, researchers nevertheless have to continue working on smallpox, and we have entered yet another round of discussing the pros and cons of smallpox vaccination. The current debates of preparedness issues are far more complex than mass vaccination, and newer vaccination strategies complicate the balance between individual and community benefits [120]. Once other infectious diseases have been eradicated, we will see similar discussions arise, but before this becomes the case, it is important to make sure that systematically collected data are aggregated and stored for posterity.

References

  1. Fenner F, Henderson DA, Arita I, Ladnyi ID: Smallpox and its eradication. 1988, Geneva: World Health Organization

    Google Scholar 

  2. Greenwood M: Epidemics and Crowd Diseases. An Introduction to the Study of Epidemiology. 1935, New York: Macmillan, (reprinted in 1977 by Arno Press, New York)

    Google Scholar 

  3. Henderson DA, Inglesby TV, Bartlett JG, Ascher MS, Eitzen E, Jahrling PB, Hauer J, Layton M, McDade J, Osterholm MT, O'Toole T, Parker G, Perl T, Russell PK, Tonat K: Smallpox as a biological weapon: medical and public health management. Working Group on Civilian Biodefense. JAMA. 1999, 281: 2127-2137. 10.1001/jama.281.22.2127.

    CAS  PubMed  Google Scholar 

  4. Aldis GK, Roberts MG: An integral equation model for the control of a smallpox outbreak. Math Biosci. 2005, 195: 1-22. 10.1016/j.mbs.2005.01.006.

    CAS  PubMed  Google Scholar 

  5. Burke DS, Epstein JM, Cummings DA, Parker JI, Cline KC, Singa RM, Chakravarty S: Individual-based computational modeling of smallpox epidemic control strategies. Acad Emerg Med. 2006, 13: 1142-1149.

    PubMed  Google Scholar 

  6. Eichner M: Case isolation and contact tracing can prevent the spread of smallpox. Am J Epidemiol. 2003, 158: 118-128. 10.1093/aje/kwg104.

    PubMed  Google Scholar 

  7. Ferguson NM, Keeling MJ, Edmunds WJ, Gani R, Grenfell BT, Anderson RM, Leach S: Planning for smallpox outbreaks. Nature. 2003, 425: 681-685. 10.1038/nature02007.

    CAS  PubMed  Google Scholar 

  8. Glasser J, Foster SO, Millar JD, Lane JM: Evaluating public health responses to reintroduced smallpox via dynamic, socially structured, and spatially distributed metapopulation models. Clin Infect Dis. 2008, 46 (suppl 3): S182-S194. 10.1086/524382.

    PubMed  Google Scholar 

  9. Kretzschmar M, Hof van den S, Wallinga J, van Wijngaarden J: Ring vaccination and smallpox control. Emerg Infect Dis. 2004, 10: 832-841.

    PubMed Central  PubMed  Google Scholar 

  10. Legrand J, Viboud C, Boelle PY, Valleron AJ, Flahault A: Modelling responses to a smallpox epidemic taking into account uncertainty. Epidemiol Infect. 2004, 132: 19-25. 10.1017/S0950268803001390.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. Longini IM, Halloran ME, Nizam A, Yang Y, Xu S, Burke DS, Cummings DA, Epstein JM: Containing a large bioterrorist smallpox attack: a computer simulation approach. Int J Infect Dis. 2007, 11: 98-108. 10.1016/j.ijid.2006.03.002.

    PubMed  Google Scholar 

  12. Meltzer MI, Damon I, LeDuc JW, Millar JD: Modeling potential responses to smallpox as a bioterrorist weapon. Emerg Infect Dis. 2001, 7: 959-969.

    PubMed Central  CAS  PubMed  Google Scholar 

  13. Nishiura H, Tang IM: Modeling for a smallpox-vaccination policy against possible bioterrorism in Japan: the impact of long-lasting vaccinal immunity. J Epidemiol. 2004, 14: 41-50. 10.2188/jea.14.41.

    PubMed  Google Scholar 

  14. Riley S, Ferguson NM: Smallpox transmission and control: spatial dynamics in Great Britain. Proc Natl Acad Sci USA. 2006, 103: 12637-12642. 10.1073/pnas.0510873103.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. Sayers BM, Angulo J: A new explanatory model of an SIR disease epidemic: a knowledge-based, probabilistic approach to epidemic analysis. Scand J Infect Dis. 2005, 37: 55-60. 10.1080/00365540410026077.

    PubMed  Google Scholar 

  16. Kaplan EH, Craft DL, Wein LM: Emergency response to a smallpox attack: the case for mass vaccination. Proc Natl Acad Sci USA. 2002, 99: 10935-10940. 10.1073/pnas.162282799.

    PubMed Central  CAS  PubMed  Google Scholar 

  17. Halloran ME, Longini IM, Nizam A, Yang Y: Containing bioterrorist smallpox. Science. 2002, 298: 1428-1432. 10.1126/science.1074674.

    CAS  PubMed  Google Scholar 

  18. 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.

    PubMed  Google Scholar 

  19. Cooper B: Poxy models and rash decisions. Proc Natl Acad Sci USA. 2006, 103: 12221-12222. 10.1073/pnas.0605502103.

    PubMed Central  CAS  PubMed  Google Scholar 

  20. Fraser C, Riley S, Anderson RM, Ferguson NM: Factors that make an infectious disease outbreak controllable. Proc Natl Acad Sci USA. 2004, 101: 6146-6151. 10.1073/pnas.0307506101.

    PubMed Central  CAS  PubMed  Google Scholar 

  21. Nishiura H, Arita I, Schwehm M, Eichner M: Epidemiological assessment of the protective effects of smallpox vaccination. Am J Infect Dis. 2006, 2: 9-17.

    Google Scholar 

  22. Lentz O, Gins HA: Handbuch der Pockenbekämpfung und Impfung. 1927, Berlin: Verlagsbuchhandlung von Richard Schoetz, (in German)

    Google Scholar 

  23. Simon J: History and practice of vaccination. Public Health Rep. 1887, 1: 167-317.

    Google Scholar 

  24. Burke A, Cunha MD: The cause of the plague of Athens:plague, typhoid, typhus, smallpox, or measles?. Infect Dis Clin North Am. 2004, 18: 29-43. 10.1016/S0891-5520(03)00100-4.

    Google Scholar 

  25. Duncan CJ, Scott S: What caused the Black Death?. Postgrad Med J. 2005, 81: 315-320. 10.1136/pgmj.2004.024075.

    PubMed Central  CAS  PubMed  Google Scholar 

  26. Dixon CW: Smallpox. 1962, London: Churchill

    Google Scholar 

  27. Rao AR: Smallpox. 1972, Bombay: Kothari Book Dept

    Google Scholar 

  28. Ricketts TF, Byles JB: The diagnosis of smallpox. 1908, London: Cassell

    Google Scholar 

  29. Regamey RH, Cohen H: International symposium on smallpox vaccine. Symposia Series in Immunobiological Standardization. 1973, Basel: Karger AS, 19:

    Google Scholar 

  30. Baxby D: The origins of vaccinia virus. J Infect Dis. 1977, 136: 453-455.

    CAS  PubMed  Google Scholar 

  31. Angulo JJ: Variola minor in Braganca Paulista County, 1956: overall description of the epidemic and of its study. Int J Epidemiol. 1976, 5: 359-366. 10.1093/ije/5.4.359.

    CAS  PubMed  Google Scholar 

  32. Chapin CV: Changes in type of infectious disease as shown by the history of smallpox in the United States, 1895–1912. J Infect Dis. 1926, 13: 171-196.

    Google Scholar 

  33. de Jong M: The alastrim epidemic in The Hague, 1953–1954. Documenta de Medicina Geographica et Tropica. 1956, 8: 207-235.

    CAS  PubMed  Google Scholar 

  34. Chapin CV, Smith J: Permanency of the mild type of smallpox. J Prev Med. 1932, 6: 273-320.

    Google Scholar 

  35. Dumbell KR, Wells DG: The pathogenicity of variola virus. A comparison of the growth of standard strains of variola major and variola minor viruses in cell cultures from human embryos. J Hyg (Lond). 1982, 89 (3): 389-397.

    CAS  Google Scholar 

  36. Dumbell KR, Huq F: The virology of variola minor. Correlation of laboratory tests with the geographical distribution and human virulence of variola isolates. Am J Epidemiol. 1986, 123: 403-415.

    CAS  PubMed  Google Scholar 

  37. Esposito JJ, Sammons SA, Frace AM, Osborne JD, Olsen-Rasmussen M, Zhang M, Govil D, Damon IK, Kline R, Laker M, Li Y, Smith GL, Meyer H, Leduc JW, Wohlhueter RM: Genome sequence diversity and clues to the evolution of variola (smallpox) virus. Science. 2006, 313: 807-812. 10.1126/science.1125134.

    PubMed  Google Scholar 

  38. Li Y, Carroll DS, Gardner SN, Walsh MC, Vitalis EA, Damon IK: On the origin of smallpox: correlating variola phylogenics with historical smallpox records. Proc Natl Acad Sci USA. 2007, 104: 15787-15792. 10.1073/pnas.0609268104.

    PubMed Central  CAS  PubMed  Google Scholar 

  39. Nishiura H: Early efforts in modeling the incubation period of infectious diseases with an acute course of illness. Emerg Themes Epidemiol. 2007, 4: 2-10.1186/1742-7622-4-2.

    PubMed Central  PubMed  Google Scholar 

  40. Mack TM: Smallpox in Europe, 1950–1971. J Infect Dis. 1972, 125: 161-169.

    CAS  PubMed  Google Scholar 

  41. Farewell VT, Herzberg AM, James KW, Ho LM, Leung GM: SARS incubation and quarantine times: when is an exposed individual known to be disease free?. Stat Med. 2005, 24: 3431-3435. 10.1002/sim.2206.

    CAS  PubMed  Google Scholar 

  42. Nishiura H: Determination of the appropriate quarantine period following smallpox exposure: an objective approach using the incubation period distribution. Int J Hyg Environ Health. 2008,

    Google Scholar 

  43. Bailey NTJ: The Mathematical Theory of Infectious Diseases and Its Applications. 1975, London: Charles Griffin, 2

    Google Scholar 

  44. Downie AW, St Vincent L, Meiklejohn G, Ratnakannan NR, Rao AR, Krishnan GN, Kempe CH: Studies on the virus content of mouth washings in the acute phase of smallpox. Bull World Health Organ. 1961, 25: 49-53.

    PubMed Central  CAS  PubMed  Google Scholar 

  45. Sakar JK, Mitra AC, Mukherjee MK, De SK, Mazumdar DG: Virus excretion in smallpox. 1. Excretion in the throat, urine, and conjunctiva of patients. Bull World Health Organ. 1973, 48: 517-522.

    Google Scholar 

  46. Nishiura H, Eichner M: Infectiousness of smallpox relative to disease age: estimates based on transmission network and incubation period. Epidemiol Infect. 2007, 135: 1145-1150.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Fine PE: The interval between successive cases of an infectious disease. Am J Epidemiol. 2003, 158: 1039-1047. 10.1093/aje/kwg251.

    PubMed  Google Scholar 

  48. Nishiura H: Epidemiology of a primary pneumonic plague in Kantoshu, Manchuria, from 1910 to 1911: statistical analysis of individual records collected by the Japanese Empire. Int J Epidemiol. 2006, 35: 1059-1065. 10.1093/ije/dyl091.

    PubMed  Google Scholar 

  49. Nishiura H, Schwehm M, Kakehashi M, Eichner M: Transmission potential of primary pneumonic plague: time inhomogeneous evaluation based on historical documents of the transmission network. J Epidemiol Community Health. 2006, 60: 640-645. 10.1136/jech.2005.042424.

    PubMed Central  PubMed  Google Scholar 

  50. Brookmeyer R, Gail MH: AIDS Epidemiology. A quantitative approach. 1994, New York: Oxford University Press

    Google Scholar 

  51. Mack T: A different view of smallpox and vaccination. N Engl J Med. 2003, 348: 460-463. 10.1056/NEJMsb022994.

    PubMed  Google Scholar 

  52. Shafa E: Case fatality ratios in smallpox. WHO Surveillance Epidemiology. 1972, 72: 35-

    Google Scholar 

  53. Dietz K, Heesterbeek JA: Daniel Bernoulli's epidemiological model revisited. Math Biosci. 2002, 180: 1-21. 10.1016/S0025-5564(02)00122-0.

    PubMed  Google Scholar 

  54. Barry FW: Report on an epidemic of small-pox at Sheffield, during 1887–88. 1889, London: Local Government Board

    Google Scholar 

  55. Rutten W: De Vreselijkste Aller Harpijen. Pokkenepidemieen en pokkenbestrijding in Nederland in de achttiende en negentiende eeuw: een social-historische en historisch-demografische studie. 1997, Wageningen, The Netherlands: Agricultural University Wageningen, (in Dutch)

    Google Scholar 

  56. Rao AR, Prahlad I, Swaminathan M, Lakshmi A: Pregnancy and smallpox. J Indian Med Assoc. 1963, 40: 353-363.

    CAS  PubMed  Google Scholar 

  57. Nishiura H: Smallpox during pregnancy and maternal outcomes. Emerg Infect Dis. 2006, 12: 1119-1121.

    PubMed Central  PubMed  Google Scholar 

  58. Pawitan Y: All Likelihood: Statistical Modelling and Inference Using Likelihood. 2001, Oxford: Oxford University Press

    Google Scholar 

  59. Lindsey JK: Nonlinear Models for Medical Statistics. 2001, Oxford: Oxford University Press, 2

    Google Scholar 

  60. Gelman A, Carlin B, Stern H, Rubin DB: Bayesian Data Analysis. 2003, Florida: Chapman and Hall/CRD, 2

    Google Scholar 

  61. Diekmann O, Heesterbeek JAP: Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. 2000, New York: Wiley

    Google Scholar 

  62. Dietz K: The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res. 1993, 2: 23-41. 10.1177/096228029300200103.

    CAS  PubMed  Google Scholar 

  63. Blower S, Bernoulli D: An attempt at a new analysis of the mortality caused by smallpox and of the advantages of inoculation to prevent it. 1766. Rev Med Virol. 2004, 14: 275-288. 10.1002/rmv.443.

    PubMed  Google Scholar 

  64. Nishiura H, Dietz K, Eichner M: The earliest notes on the reproduction number in relation to herd immunity: Theophil Lotz and smallpox vaccination. J Theor Biol. 2006, 241: 964-967.

    PubMed  Google Scholar 

  65. Lloyd-Smith JO, Schreiber SJ, Kopp PE, Getz WM: Superspreading and the effect of individual variation on disease emergence. Nature. 2005, 438: 355-359. 10.1038/nature04153.

    CAS  PubMed  Google Scholar 

  66. Gani R, Leach S: Transmission potential of smallpox in contemporary populations. Nature. 2001, 414: 748-751. 10.1038/414748a.

    CAS  PubMed  Google Scholar 

  67. Becker NG: Analysis of Infectious Disease Data. 1989, New York: Chapman and Hall

    Google Scholar 

  68. Fine PE: Herd immunity: history, theory, practice. Epidemiol Rev. 1993, 15: 265-302.

    CAS  PubMed  Google Scholar 

  69. World Health Assembly: Handbook of resolutions and decisions of the twelfth World Health Assembly and Executive Board. (12.54). 1959, Geneva: World Health Organization

    Google Scholar 

  70. Anderson RM, May RM: Directly transmitted infectious diseases: control by vaccination. Science. 1982, 215: 1053-1060. 10.1126/science.7063839.

    CAS  PubMed  Google Scholar 

  71. Smith CE: Factors in the transmission of virus infections from animal to man. Sci Basis Med Ann Rev. 1964, 125-150.

    Google Scholar 

  72. Haber M, Longini IM, Halloran ME: Measures of the effects of vaccination in a randomly mixing population. Int J Epidemiol. 1991, 20: 300-310. 10.1093/ije/20.1.300.

    CAS  PubMed  Google Scholar 

  73. Halloran ME, Haber M, Longini IM, Struchiner CJ: Direct and indirect effects in vaccine efficacy and effectiveness. Am J Epidemiol. 1991, 133: 323-331.

    CAS  PubMed  Google Scholar 

  74. Greenwood M, Yule GU: The statistics of anti-typhoid and anti-cholera inoculations, and the interpretation of such statistics in general. Proc R Soc Med. 1915, 8: 113-194.

    PubMed Central  CAS  PubMed  Google Scholar 

  75. Orenstein WA: Assessing vaccine efficacy in the field. Further observations. Epidemiol Rev. 1988, 10: 212-241.

    CAS  PubMed  Google Scholar 

  76. Fine PE, Clarkson JA, Miller E: The efficacy of pertussis vaccines under conditions of household exposure. Further analysis of the 1978–80 PHLS/ERL study in 21 area health authorities in England. Int J Epidemiol. 1988, 17: 635-642. 10.1093/ije/17.3.635.

    CAS  PubMed  Google Scholar 

  77. Farrington CP: On vaccine efficacy and reproduction numbers. Math Biosci. 2003, 185: 89-109. 10.1016/S0025-5564(03)00061-0.

    CAS  PubMed  Google Scholar 

  78. Halloran ME, Haber M, Logini IM: Interpretation and estimation of vaccine efficacy under heterogeneity. Am J Epidemiol. 1992, 136: 328-343.

    CAS  PubMed  Google Scholar 

  79. Rao AR, Jacob ES, Kamalakshi S, Appaswamy S, Bradbury : Epidemiological studies in smallpox. A study of intrafamilial transmission in a series of 254 infected families. Indian J Med Res. 1968, 56: 1826-1854.

    CAS  PubMed  Google Scholar 

  80. Angulo JJ, Walter SD: Variola minor in Braganca Paulista County, 1956: household aggregation of the disease and the influence of household size on the attack rate. J Hyg (Lond). 1979, 82 (1): 1-6.

    CAS  Google Scholar 

  81. Heiner GG, Fatima N, McCrumb FR: A study of intrafamilial transmission of smallpox. Am J Epidemiol. 1971, 94: 316-326.

    CAS  PubMed  Google Scholar 

  82. Mukherjee MK, Sarkar JK, Mitra AC: Pattern of intrafamilial transmission of smallpox in Calcutta, India. Bull World Health Organ. 1974, 51: 219-225.

    PubMed Central  CAS  PubMed  Google Scholar 

  83. Becker NG, Dietz K: The effect of household distribution on transmission and control of highly infectious diseases. Math Biosci. 1995, 127: 207-219. 10.1016/0025-5564(94)00055-5.

    CAS  PubMed  Google Scholar 

  84. Ball F, Mollison D, Scalia-Tomba G: Epidemics with two levels of mixing. Ann Appl Prob. 1997, 7: 46-89. 10.1214/aoap/1034625252.

    Google Scholar 

  85. Edmunds WJ, O'Callaghan CJ, Nokes DJ: Who mixes with whom? A method to determine the contact patterns of adults that may lead to the spread of airborne infections. Proc Roy Soc B. 1997, 264: 949-957. 10.1098/rspb.1997.0131.

    CAS  Google Scholar 

  86. Cliff AD, Haggett P, Ord JK, Versey GR: Spatial Diffusion. An Historical Geography of Epidemics in and Island Community. 1981, Cambridge: Cambridge University Press

    Google Scholar 

  87. Abbey H: An examination of the Reed-Frost theory of epidemics. Hum Biol. 1952, 24: 201-233.

    CAS  PubMed  Google Scholar 

  88. Massad E, Burattini MN, Lopez LF, Coutinho FA: Forecasting versus projection models in epidemiology: the case of the SARS epidemics. Med Hypotheses. 2005, 65: 17-22. 10.1016/j.mehy.2004.09.029.

    PubMed  Google Scholar 

  89. Nishiura H: Time variations in the transmissibility of pandemic influenza in Prussia, Germany, from 1918–19. Theor Biol Med Model. 2007, 4: 20-10.1186/1742-4682-4-20.

    PubMed Central  PubMed  Google Scholar 

  90. Del Valle S, Hethcote H, Hyman JM, Castillo-Chavez C: Effects of behavioral changes in a smallpox attack model. Math Biosci. 2005, 195: 228-251. 10.1016/j.mbs.2005.03.006.

    CAS  PubMed  Google Scholar 

  91. Hall MJ, Norwood AE, Ursano RJ, Fullerton CS: The psychological impacts of bioterrorism. Biosecurity Bioterrorism. 2003, 1: 139-144. 10.1089/153871303766275817.

    PubMed  Google Scholar 

  92. Hammarlund E, Lewis MW, Hansen SG, Strelow LI, Nelson JA, Sexton GJ, Hanifin JM, Slifka MK: Duration of antiviral immunity after smallpox vaccination. Nature Med. 2003, 9: 1131-1137. 10.1038/nm917.

    CAS  PubMed  Google Scholar 

  93. Simpson EL, Hercher M, Hammarlund EK, Lewis MW, Slifka MK, Hanifin JM: Cutaneous responses to vaccinia in individuals with previous smallpox vaccination. J Am Acad Dermatol. 2007, 57: 442-444. 10.1016/j.jaad.2007.04.022.

    PubMed  Google Scholar 

  94. Nishiura H, Schwehm M, Eichner M: Still protected against smallpox? Estimation of the duration of vaccine-induced immunity against smallpox. Epidemiology. 2006, 17: 576-581. 10.1097/01.ede.0000229196.41862.c2.

    PubMed  Google Scholar 

  95. Eichner M: Analysis of historical data suggests long-lasting protective effects of smallpox vaccination. Am J Epidemiol. 2003, 158: 717-723. 10.1093/aje/kwg225.

    PubMed  Google Scholar 

  96. Cohen J: Bioterrorism. Smallpox vaccinations: how much protection remains?. Science. 2001, 294: 985-10.1126/science.294.5544.985.

    CAS  PubMed  Google Scholar 

  97. Nishiura H, Eichner M: Estimation of the duration of vaccine-induced residual protection against severe and fatal smallpox based on secondary vaccination failure. Infection. 2006, 34: 241-246. 10.1007/s15010-006-6603-5.

    CAS  PubMed  Google Scholar 

  98. Halloran ME: Concepts of infectious disease epidemiology. Modern Epidemiology. Edited by: Rothman KJ, Greenland S. 1998, New York: Lippincott Williams and Wilkins, 529-554. 2

    Google Scholar 

  99. Toyoda T: Smallpox research. Part II. Findings based on 516 smallpox cases with special reference to hemorrhagic type and variola sine exanthemata. Tokyo Med J. 1923, 2333: 1-12. (in Japanese)

    Google Scholar 

  100. Foege WH, Millar JD, Lane JM: Selective epidemiological control in smallpox eradication. Am J Epidemiol. 1971, 94: 311-315.

    CAS  PubMed  Google Scholar 

  101. Henderson DA: Epidemiology in the global eradication of smallpox. Am J Epidemiol. 1972, 1: 25-30. 10.1093/ije/1.1.25.

    CAS  Google Scholar 

  102. Foege WH, Millar JD, Henderson DA: Smallpox eradication in West and Central Africa. Bull World Health Organ. 1975, 52: 209-222.

    PubMed Central  CAS  PubMed  Google Scholar 

  103. Kaplan EH, Wein LM: Smallpox eradication in West and Central Africa: surveillance-containment or herd immunity?. Epidemiology. 2003, 14: 90-92. 10.1097/00001648-200301000-00021.

    PubMed  Google Scholar 

  104. Porco TC, Holbrook KA, Fernyak SE, Portnoy DL, Reiter R, Aragon TJ: Logistics of community smallpox control through contact tracing and ring vaccination: a stochastic network model. BMC Public Health. 2004, 4: 34-10.1186/1471-2458-4-34.

    PubMed Central  PubMed  Google Scholar 

  105. Meltzer MI: Risks and benefits of preexposure and postexposure smallpox vaccination. Emerg Infect Dis. 2003, 9: 1363-1370.

    PubMed Central  PubMed  Google Scholar 

  106. Nishiura H: Analysis of a previous smallpox vaccination study: estimation of the time period required to acquire vaccine-induced immunity as assessed by revaccination. Southeast Asian J Trop Med Public Health. 2006, 37: 673-680.

    PubMed  Google Scholar 

  107. Mortimer PP: Can postexposure vaccination against smallpox succeed?. Clin Infect Dis. 2003, 36: 622-629. 10.1086/374054.

    PubMed  Google Scholar 

  108. Eichner M, Schwehm M: Smallpox. A vulnerable specter. Epidemiology. 2004, 15: 258-261. 10.1097/01.ede.0000121548.84183.04.

    PubMed  Google Scholar 

  109. Kennedy JS, Frey SE, Yan L, Rothman AL, Cruz J, Newman FK, Orphin L, Belshe RB, Ennis FA: Induction of human T cell-mediated immune responses after primary and secondary smallpox vaccination. J Infect Dis. 2004, 190: 1286-1294. 10.1086/423848.

    PubMed  Google Scholar 

  110. Massoudi MS, Barker L, Schwartz B: Effectiveness of postexposure vaccination for the prevention of smallpox: results of a Delphi analysis. J Infect Dis. 2003, 188: 973-976. 10.1086/378357.

    PubMed  Google Scholar 

  111. Nishiura H, Eichner M: Interpreting the epidemiology of postexposure vaccination against smallpox. Int J Hyg Env Health. 2008, 211: 219-226. 10.1016/j.ijheh.2007.01.029.

    Google Scholar 

  112. Millard CK: The vaccination question in the light of modern experience. An appeal for reconsideration. 1914, London: H.K. Lewis

    Google Scholar 

  113. Fulginiti VA, Papier A, Lane JM, Neff JM, Henderson DA: Smallpox vaccination: a review, part II. Adverse events. Clin Infect Dis. 2003, 37: 251-271. 10.1086/375825.

    PubMed  Google Scholar 

  114. Kretzschmar M, Wallinga J, Teunis P, Xing S, Mikolajczyk R: Frequency of adverse events after vaccination with different vaccinia strains. PLoS Med. 2006, 3: e272-10.1371/journal.pmed.0030272.

    PubMed Central  PubMed  Google Scholar 

  115. Mauer DM, Harrington B, Lane JM: Smallpox vaccine: contraindications, administration, and adverse reactions. Am Fam Physician. 2003, 68 (5): 889-896.

    Google Scholar 

  116. Mack TM, Thomas DB, Khan MM: Variola major in West Pakistan. J Infect Dis. 1970, 122: 479-488.

    CAS  PubMed  Google Scholar 

  117. Mack TM, Thomas DB, Muzaffar Khan M: Epidemiology of smallpox in West Pakistan. II. Determinants of intravillage spread other than acquired immunity. Am J Epidemiol. 1972, 95: 169-177.

    CAS  PubMed  Google Scholar 

  118. Hall IM, Egan JR, Barrass I, Gani R, Leach S: Comparison of smallpox outbreak control strategies using a spatial metapopulation model. Epidemiol Infect. 2007, 135: 1133-1144.

    PubMed Central  CAS  PubMed  Google Scholar 

  119. Breman JG, Arita I: The confirmation and maintenance of smallpox eradication. N Engl J Med. 1980, 303: 1263-1273.

    CAS  PubMed  Google Scholar 

  120. Bauch CT, Galvani AP, Earn DJ: Group interest versus self-interest in smallpox vaccination policy. Proc Natl Acad Sci USA. 2003, 100: 10564-10567. 10.1073/pnas.1731324100.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This review would not have been possible without technical support and input from Klaus Dietz and Isao Arita. This study was in part supported by European Union project INFTRANS (FP6 STREP; contract no. 513715). The study of HN was supported by The Netherlands Organisation for Scientific Research (NWO, ALW-IPY-NL/06-15D).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Hiroshi Nishiura or Martin Eichner.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HN reviewed the literature, analyzed the data and drafted an early version of the manuscript. ME reviewed the early version of the manuscript and assisted in editing the manuscript. SOB participated in the writing and revision of the manuscript. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Nishiura, H., Brockmann, S.O. & Eichner, M. Extracting key information from historical data to quantify the transmission dynamics of smallpox. Theor Biol Med Model 5, 20 (2008). https://doi.org/10.1186/1742-4682-5-20

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-4682-5-20

Keywords