 |
 |
 |
 |
|
|
|
|
 |
|
|
 |
|
|
|
|
| |
|||||||||||||||||
|
|
|
|
|
|
|
|||||||||||
|
|||||||||||||||||
 |
|
EID Home | Ahead of Print | Past Issues | EID Search | Contact Us | Announcements | Suggested Citation | Submit Manuscript
|
|
|
Research Ring Vaccination and Smallpox ControlMirjam Kretzschmar,*
Recently, concerns about a bioterror attack with the smallpox virus or other infectious disease agents have risen (1,2). While new vaccines with fewer adverse consequences are being developed (3), the existing vaccines, which have potential side effects and may be lethal, are the only vaccines available (4,5). In the United States during the recent voluntary smallpox vaccination program, a limited number of healthcare workers volunteered for vaccination because of the risks associated with vaccination and the low risk for infection (6). If an outbreak occurs, vaccination strategies include ring vaccination around diagnosed cases of smallpox or a mass vaccination to begin as soon as the first cases are diagnosed. Without natural smallpox infections, practical experience with ring vaccination against smallpox cannot be gained; accounts of the vaccination programs that eradicated smallpox in the 1970s are the only source of information (7). Combined with information collected during the last decades of smallpox circulation, mathematical modeling offers a tool to explore various vaccination scenarios if an outbreak occurs (8–15). We investigated which conditions are the best for effective use of ring vaccination, a strategy in which direct contacts of diagnosed cases are identified and vaccinated. We also investigated whether monitoring contacts contributes to the success of ring vaccination. We used a stochastic model that distinguished between close and casual contacts to explore the variability in the number of infected persons during an outbreak, and the time until the outbreak is over. We derived expressions for the basic reproduction number (R0) and the effective reproduction number (Rυ). We investigated how effectiveness of ring vaccination depends on the time until diagnosis of a symptomatic case, the time to identify and vaccinate contacts in the close contact and casual contact ring, and the vaccination coverage required to contain an epidemic. MethodsThe model describes the number of infected persons after one or more index cases are introduced. It simulates a stochastic process in which every infected person generates a number of new infections according to a given probability distribution. This process implies that contacts of different infected persons are independent of each other and that no saturation of the incidence occurs at higher prevalence. The model is applicable for the first few generations of infection, if the outbreak goes unchecked, and for the complete outbreak if it is contained. We summarize the main features of the model; the formal model definition is given in the appendix. Course of Infection and Transmission
The noninfectious state (incubation period plus prodromal phase) lasts 12–15 days (7,16,17) with specified probabilities per day of moving to the infectious state. The assumption that infectivity during the prodromal phase is negligible is supported by a recently published statistical analysis of outbreak data (16). The duration of the infectious state DI is 14 days (7,16,18), with variable infectiousness during that time (8,13,19). The probability of transmission per contact pτ , where τ denotes the day of the infectious period, is high at the beginning and low at the end of the infectious period (Figure 1A). At the end of the infectious period, a person either recovers or dies. The case-fatality rate is 30% (14), which is an average value for the case-fatality rate of variola major. Transmission takes place in two rings of contacts: 1) household and other close contacts, and 2) more casual face-to-face contacts. We assumed that in the close contact ring the probability of transmission is five times higher than in the casual contact ring (g = 0.2). The number of contacts on day τ of the infectious period in the close contact ring follows a Poisson distribution with mean μτ(1), and in the casual contact ring this number follows a a negative binomial distribution with mean μτ(2). The values (μτ(1) = 2 and μτ(2) = 14.9) were chosen such that the total number of contacts per day was comparable to numbers observed in empirical studies (20) (Figure 1B). For every contact, the event of transmission is determined by the infectiousness by day of the infectious period. Most people can be infected again, and smallpox can develop 10–20 years after vaccination (21,22). While residual immunity might lower the case-fatality rate, it might also lead to a later diagnosis for infected persons because disease symptoms are milder. An infection with milder symptoms is probably less infectious, but an infectious person might have more contacts with others because he or she feels less impaired by disease symptoms. Hence, the net effects of residual immunity are difficult to assess. We assumed that all persons are equally susceptible, and that no protective immunity remains in the population from previous vaccination. The basic reproduction number R0 describes the average number of secondary cases produced from contact with an infected person during the infectious period and without intervention. The number can be computed as the sum of the reproduction numbers in the close contact and casual contact ring
For the baseline parameter values given in the Table, R0
= 5.23, which, when broken down by rings of contacts, gives Ring VaccinationRing vaccination in the model includes complete isolation of symptomatic patients with diagnosed cases of smallpox and vaccination of some or all contacts of the diagnosed case-patient. In our baseline scenario, we assumed that vaccinated contacts are not isolated after vaccination and may therefore transmit the infection to others if they become infectious. In addition, we enhance the baseline intervention by including monitoring of identified contacts. The effectiveness of the intervention therefore is determined by the probability of diagnosis per day of the infectious period, the time needed to identify contacts of the close contact and casual contact ring, the vaccination coverage in the close contact and the casual contact ring, and whether contacts are monitored. Some of those parameters (speed of diagnosis and time to identifying contacts) differ between the first index case in the population and cases occurring later in the epidemic. In Figure 2, the timing of the key events in the chain of transmission and intervention is shown schematically. The index patient can cause new cases of infection between the beginning of the infectious period until diagnosis and isolation. For a secondary case, vaccination has to take place within 4 days after infection (14) to prevent disease. We denote with δτ the probability of diagnosis on day τ of the infectious period for those persons who have not been diagnosed before. From those probabilities, one can derive the probability that an infectious person is not yet diagnosed on day τ of his or her infectious period (Figure 1C). By υτ(i), we denote the probability that a contact in ring i (i = 1 or 2), who was infected on day τ of the index patient's infectious period, will be vaccinated within 4 days of being infected. In the appendix, υτ depends on the diagnosis probabilities, the time needed for contact tracing, and the vaccination coverage. Throughout, we assume that the vaccine efficacy is 100%. We can now determine an effective reproduction number Rυ that describes the number of secondary cases caused by an index patient in a situation with intervention:
A special strategy included in this formula is an intervention where only case isolation is performed without vaccination of contacts. The vaccination coverage c(i) is set to zero. Equivalently, it describes the situation that no window period exists (24). The reproduction number can then be calculated as
If vaccination is ineffective, but contacts are monitored, the monitoring will have the same effect on R0 as an effective vaccination, because the contacts will not be able to disseminate the virus any further (assuming a fully effective monitoring). Therefore, Rv can be computed with the formula including vaccination, where the window period is now set to w = 15, of the full duration of the infectious period. This assumption means that regardless of when the index patient's condition is diagnosed, contacts can effectively be excluded from further transmission. If monitoring of contacts is not 100% effective, the parameter c(i) for the coverage can be used to express the extent of successful monitoring. The outbreak can be controlled if Rυ <1. In the Table, the model parameters and their baseline values are listed. In the Appendix, the formal model definition is given. ResultsBaseline Parameter SetAn epidemic starting with one index case in a completely susceptible population without intervention grows exponentially, if it survives early extinction. The large range of possible courses of the epidemic reflects the stochastic variability (Figure 1D). If the intervention does not succeed in reducing the effective reproduction number Rυ to below 1, the epidemic will continue to grow exponentially, albeit at a lower rate. For example, if diagnosed infectious persons are isolated, but no ring vaccination is performed, the effective reproduction number is Rυ=1.65 and the epidemic cannot be contained. If the intervention succeeds in reducing the effective reproduction number Rυ to <1, the size of successive generations of infected persons declines. For the parameter values of the baseline scenario given in the Table, we have Rυ = 0.67 and 94.4% of all transmissions take place in the casual contact ring. The epidemic can then be contained and the virus eradicated. In Figure 3A and B, the distribution of the total number of infected persons (excluding those who were vaccinated in time to prevent symptomatic infection) and of the time until recovery of the last infected patient is shown for 500 simulation runs with the baseline parameters (Table). The time until the epidemic is over is quite variable: on average it takes 82 days, in some cases it takes up to 1 year (range 22–334 days). During this time, an average of 209 contacts are vaccinated (range 6–1,038 contacts), with a mean number of 13.5 vaccinated contacts per infected case. On average, 15 persons are infected (range 1–123 persons), if we exclude those infected persons who were vaccinated on time. If we include the infected contacts who were vaccinated on time, the mean number of persons infected is 29 (range 1–200 contacts). On average, 5 persons die of smallpox (range 0–30 persons). If monitoring of contacts is added to the intervention (incorporated in the model by assuming that w = 15), the effective reproduction number can be further reduced to 0.53. The average number of infected persons drops to 6 (range 1–43 persons), excluding the identified infected contacts who are vaccinated or monitored and to 14 (range 1–96 persons), including the vaccinated and monitored infected contacts. The mean time to extinction is now 53 days (range 20–288 days). The fraction of transmissions taking place in the casual contact ring is slightly lower. at 93.7%. To contrast the baseline scenario, in Figure 3C and D we show results for the case that R0 = 10.46, i.e., twice the value of baseline scenario. To contain the epidemic, we now assumed that 80% of all contacts in the casual contact ring were vaccinated in time. The effective reproduction number Rυ was 0.91. A fraction of 91.8% of transmissions took place in the casual contact ring. The mean number of infected persons during the epidemic was 101 (range 2–663), excluding the infected contacts vaccinated in time and 340 (range 2–2,175 persons) including those contacts. The mean time to extinction was 229 days (range 32 to >900 days). The time to extinction can be very long when Rυ is near 1 because the epidemic can flare up again when a case by chance produces many secondary infections. A similar picture would result if, in addition to vaccinating casual contacts with a coverage of 55%, those contacts are monitored. The effective reproduction number is then 0.96; 93.0% of transmissions are in the casual contact ring. Sensitivity AnalysisInitial Phase of EpidemicThe initial phase of the epidemic (time before discovery of the first case) is determined by the number of index patients that start the epidemic outbreak and by the time it takes to diagnose the first case. We varied those two variables separately while assuming that after diagnosis intervention took place within the parameters defined in the Table, i.e., with an Rυ of 0.67 (Figure 4). While the number of cases increases almost linearly with the number of index cases, the dependency on the time to diagnosis shows the influence of the variable infectiousness during the infectious period. In the beginning, when infectiousness is high, the number of infected persons increases rapidly. Another rise occurs toward the end of the index patient's infectious period because the second-generation patients become infectious and produce the third generation of infected persons. Once the second generation of infected persons has the opportunity to disseminate the infection further, the range of possible outcomes increases greatly (range 6–221 infected persons). A similar picture emerges for the time needed to extinguish the outbreak. The duration of the outbreak increases when diagnosis is delayed during the first few days of infectiousness, then stays on a stable level, and finally increases again when diagnosis is delayed towards the end of the infectious period (Figure 4D). Therefore, diagnosing the first index case before the second generation of infected persons start transmitting the virus is important. If diagnosing the first case at the beginning of its infectious period is possible, the number of cases during the epidemic can be kept at a low level. Effectiveness of ResponseAmong others, the value of R0 determines whether ring vaccination as defined above can contain an epidemic or not. As the value of R0 is uncertain (2,9,19), we studied Rυ as a function of R0 (Figure 5). In Figure 5A, an intervention without monitoring of contacts is considered with various assumptions on how long tracing and vaccinating casual contacts take. In this case, ring vaccination can contain the epidemic if R0 is <7, and contacts can be traced within 3 days. In Figure 5B, monitoring of contacts is added to the intervention. In that instance, the epidemic can be contained up to an R0 of 10. In Figure 5A and B, we assumed that 50% of all casual contacts can be identified and vaccinated or monitored. In Figure 6, we show how the critical vaccination coverage needed to control the epidemic depends on R0, or, more specifically, on the average number of daily contacts (Figure 6C). In addition, we varied the baseline assumption about the time to diagnosis by shifting the probability of being diagnosed by n days towards a later time in the symptomatic period. In Figure 6A, without monitoring of vaccinated contacts, a shift by 1 day greatly increases the coverage needed to contain the epidemic. If diagnosis is delayed by >1 day, the chances of controlling the epidemic diminish greatly. If vaccinated contacts are monitored, the situation improves (Figure 6B), and a high vaccination coverage in the casual contact ring ensures that the epidemic stays under control. Time to Extinction Depending on RυFinally, we looked at how differences in intervention effectiveness influence the duration of the epidemic and the cumulative number of infected persons. The effective reproduction number Rυ was varied by decreasing the vaccination coverage in the casual contact ring stepwise to 0.1. The effective reproduction number increased up to 0.96 (starting from the baseline value of 0.67). Figures 7A and B show how the cumulative number of infected persons and the time to extinction increase with increasing Rυ. The mean time until extinction approximately doubles to almost 200 days, and the range of possible outcomes increases with maximum possible durations of >2 years. The mean number of infected persons increases by a factor of 5, and the range of possible outcomes increases such that epidemics with several hundreds of infected persons are possible. Hence, if Rυ is slightly <1, the epidemic might take a long time to control, and the number of persons who become infected and die might be high. DiscussionOur simulation results show that a smallpox epidemic starting from a small number of index cases can be contained by ring vaccination provided the intervention measures are very effective. The time to diagnosis has proven to be an essential and sensitive parameter in determining the intervention effectiveness. The speed of diagnosis is less essential if identified contacts are isolated to prevent them from transmitting further if their vaccination fails. The time window limiting the success of vaccination then loses its importance for determining the effectiveness of intervention. The time to diagnosis of cases and the fraction of contacts found by contact tracing are then the key parameters. Quick contact tracing would be even more essential if substantial transmission would take place during the prodromal period of infection as is assumed by some authors (11,12). Some limitations of our modeling approach should be kept in mind. First, we only consider epidemics that are started by a small number of index cases. The branching process approach does not allow for overlapping rings of contact, but we implicitly include such an effect by varying the effective transmission probability such that the distribution of transmissions over the infectious period agrees with empiric findings (8). In other words, the decreasing probability of contacting new susceptible persons during the infectious period is incorporated in the decreasing transmission probability per contact. For larger numbers of index cases, our approach can be viewed as a worst-case scenario. Second, we assume that the population is completely susceptible, i.e., no residual immunity from vaccination in the pre-eradication era exists. This lack of immunity means that if previously vaccinated persons cannot become infectious to others, our results are too pessimistic, whereas if they become infectious with mitigated symptoms, our results might be too optimistic. In the recent literature, other models, both stochastic and deterministic, of smallpox outbreaks have been introduced to analyze the effects of ring and mass vaccination (8,10–14). On the basis of a low estimate for R0, Meltzer et al. (8) concluded that even quarantine alone can control the epidemic, as can vaccination alone, if the transmission rate is reduced sufficiently. However, the model does not allow analysis of how intervention parameters determine the reduction of the transmission rate. On the other hand, Kaplan et al. (10) conclude that with a large number of initial cases, mass vaccination will prevent more deaths than will a vaccination strategy based on contact tracing. Those authors explicitly take into account the limited resources available for tracing and vaccinating contacts, a limitation that is an important factor in large outbreaks. Also, they assume that infectivity is high during the prodromal phase. As in the model by Kaplan et al. (10), our model takes into account that the intervals between infection and diagnosis of an index case, and between diagnosis of an index case and tracing of the contact, may exceed the time window in which vaccination has to take place. This window limits the possible effectiveness of contact vaccination—a phenomenon termed "race to trace" by Kaplan et al. While the model by Kaplan et al. is based on differential equations with exponentially distributed sojourn times in different compartments, our model is a stochastic model that is able to deal with more realistic distributions for sojourn times in different disease states. Also, our model can provide estimates for variability in outcomes (Discussion in [25]). In a study by Halloran et al. (11), a stochastic model for smallpox outbreaks in small, structured communities is described. In some respects the model is similar to ours, namely, that there is a distinction between household contacts and other contacts in the community with differing transmission probabilities. The values of most biologic parameters are choices similar to ours, with the exception of Hallorans assumption that persons are highly infectious during the prodromal phase. An important difference between the models is the natural limitation of the number of infected persons in any epidemic attributable to the rather small community size in Halloran's model. The nonlinear effects of saturation play a rather large role in determining the outbreak size, especially for less effective intervention measures. Also, Halloran's model seems to be too complex to derive an explicit formula for the basic reproduction number R0, and thus makes sensitivity analysis of the results based on that quantity much more tedious. One of the main differences, however, is the way that vaccination is incorporated into the model does not allow investigation of the effects of those parameters that largely determine the success of intervention, namely the time to diagnosis of new cases and the time needed to trace contacts. In a small closed community as the one described in the model, tracing of contacts is not that difficult, but in modern society with its increasing mobility the tracing of casual contacts can pose a big problem. The main difference between the modeling approach of Bozette et al. (12,13) and our model is that in Bozette's R0 and Rυ are set to prescribed values, while in our model those numbers can be derived from measurable quantities inherent to the transmission and intervention process. In comparison to the estimate of Rυ of 0.53 in our baseline scenario with vaccination and monitoring of contacts, Bozette et al. assume a much lower value of 0.1. So compared to our results, their results are rather optimistic, but they cannot relate the assumed value of Rυ to R0 or to parameters describing the process of contact tracing and vaccination. Finally, Eichner (14) recently published a modeling study that uses a simulation model to assess the effectiveness of case isolation and contact tracing. Modeling approach and choice of parameter values resemble our approach, but the intervention is modeled in a more phenomenologic way by defining the outcomes of intervention without explicitly including intervention-related parameters into the model. This does not allow for an explicit calculation of the effective reproduction number based on intervention parameters as is possible in our approach. With respect to preparing for a smallpox outbreak, alertness and ability to diagnose quickly are important. Physicians and nurses need to be educated and the public needs to be more aware. Also, since we know little about the timing and effectiveness of identifying infectious persons and their contacts in case of a bioterror attack, obtaining more empirical information about contact patterns and contact tracing will be helpful. Recently, some useful data about contact patterns have been collected during severe acute respiratory syndrome outbreaks, but a more systematic investigation of contact tracing is advisable. Considering the uncertainties connected to all parameter values, we conclude that any contingency plan for use of ring vaccination must also identify the criteria under which switching to large-scale mass vaccination is justified. Acknowledgments
AppendixWe used a stochastic model to describe the number of infected persons following the introduction of one or more index cases. The model simulates a discrete time-branching process. Formal Model DefinitionBy Et,τ we denote the number of persons infected at time t-τ who are still latent at time t, 1 < τ < DE. By It,τ we denote the number of persons who became symptomatic at time t-τ and are not yet in isolation at time t, 1 < τ < DI. By Qt,τ we denote the number of persons who became symptomatic at time t-τ and are in isolation at time t, 1 < τ < DI. Sτ denotes a Bernoulli distributed random variable with mean γτ, 1 < τ < DE, where γτ is the probability to move to the infectious state on day τ of the latent period. DE is the maximum duration of the latent period. Tτ denotes a Bernoulli distributed random variable
with mean pτ, 1 < τ < DI,
where pτ is the transmission probability per contact
on day τ of the infectious period. DI is the duration
of the infectious period. We assume that pτ can
be well described by the functional form Cτ(1) denotes a Poissson distributed
random variable with mean Cτ(2) denotes a random variable with
a negative binomial distribution with parameters nτ
and qτ, 1 < τ < DI
describing the number of contacts per day in the casual contact ring.
The mean number of contacts per day in the casual contact ring is The transition from being undiagnosed infectious to isolation depends on the probability of diagnosis Δτ on day τ after the start of the infectious period. Δτ is a Bernoulli distributed random variable with parameter δτ, τ=1,...,DI. Vτ(1) and Vτ(2) are
Bernoulli distributed random variables with parameters υτ(1)
and υτ(1), respectively, with τ=1,...,
DI. They describe the probability that a contact in ring
1 or 2 , respectively, will be effectively vaccinated within the time
window of 4 days after infection. The subscript τ refers to day of
the infectious period of the index case at the moment of transmission.
The vaccination probabilities υτ(i) for
i=1,2 depend on the probability that the index case is diagnosed and on
the vaccination coverage as follows. The probability
The probability for the infected contact to be vaccinated in time depends on the number of days w after infection that vaccination can still be effective, the number of days r(i) that are needed for tracing the contact, and the coverage c(i) (the fraction of contacts in ring i that are effectively immunized) for i=1,2. Then we get
To describe the probability of diagnosis per day of the infectious period we use the functional form
with b1 > 0, and an integer b2 between 0 and DI. In the simulations different values for those parameters are chosen for the first index patient, who starts the epidemic, and later cases assuming that it takes longer to diagnose the first case as there is not yet that much alertness of the public health system. Also, the time needed to find contacts and the vaccination coverage may vary between the first index case and later cases. The transitions through the different stages in time is described by the system of difference equations.
The inflow of new latent and symptomatically infected persons is given by
The initial conditions for the case that the epidemic is started by one infected index case entering the population at the beginning of his latent period are given by E0,1 = 1, E0,τ = 0 for τ = 2,…, DE , I0,τ = 0 for τ = 1,…, DI , and Q0,τ = 0 for τ = 1,…, DI . At the end of the infectious period an infected patient either recovers and becomes immune or dies. Death occurs in a fraction f of all cases, i.e., f denotes the case fatality. This implies that the number of deaths Mt+1 at time t+1 is given by
where Z is a Bernoulli distributed random variable with parameter f
, the case-fatality rate. The cumulative mortality from the start of the
epidemic up to time t is given by The model was implemented and run in Mathematica 4.2. References
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
||||||||||||
|
||||||||||||
|
|
|
|
|
EID Home | Top of Page | Ahead-of-Print | Past Issues | Suggested Citation | EID Search | Contact Us | Accessibility | Privacy Policy Notice | CDC Home | CDC Search | Health Topics A-Z |
|
|
|
|
This page posted April
1, 2004 |
|
|
|
|
Emerging
Infectious Diseases Journal |
|
|
and
, i.e., 40.2% of all transmissions take
place in the close contact ring. We use the parameter a1
in the function describing the transmission probability (
with
a1, a2 > 0.
, 1 < τ <
DI describing the number of contacts per day in the close
contact ring.
and the standard deviation
is 
.
an infectious person who has transmitted
to a contact on day τ of the infectious period is diagnosed on day
τ+j can be computed as
for j=1,..., DI-τ+1.
for i=1,2.
for τ > b2,
and 
for τ < b2
,
.
