PLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, USA06-PONE-RA-00078R110.1371/journal.pone.0000012Research ArticleComputational BiologyInfectious Diseases/Epidemiology and Control of Infectious DiseasesPublic Health and Epidemiology/Infectious DiseasesThe Effectiveness of Contact Tracing in Emerging EpidemicsThe Effectiveness of TracingKlinkenbergDon^{1}^{*}FraserChristophe^{2}HeesterbeekHans^{1}Theoretical Epidemiology, Department of Farm Animal Health, Faculty of Veterinary Medicine, Utrecht University, Utrecht, NetherlandsDepartment of Infectious Disease Epidemiology, Imperial College London, London, United KingdomGetzWayneAcademic EditorUniversity of California, Berkeley, United States of America* To whom correspondence should be addressed. E-mail: d.klinkenberg@vet.uu.nl
Conceived and designed the experiments: CF HH DK. Performed the experiments: DK. Analyzed the data: DK. Wrote the paper: CF HH DK
The authors have declared that no competing interests exist.
201220062012200611e12258200615920062006Klinkenberg et alThis is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Background
Contact tracing plays an important role in the control of emerging infectious diseases, but little is known yet about its effectiveness. Here we deduce from a generic mathematical model how effectiveness of tracing relates to various aspects of time, such as the course of individual infectivity, the (variability in) time between infection and symptom-based detection, and delays in the tracing process. In addition, the possibility of iteratively tracing of yet asymptomatic infecteds is considered. With these insights we explain why contact tracing was and will be effective for control of smallpox and SARS, only partially effective for foot-and-mouth disease, and likely not effective for influenza.
Methods and Findings
We investigate contact tracing in a model of an emerging epidemic that is flexible enough to use for most infections. We consider isolation of symptomatic infecteds as the basic scenario, and express effectiveness as the proportion of contacts that need to be traced for a reproduction ratio smaller than 1. We obtain general results for special cases, which are interpreted with respect to the likely success of tracing for influenza, smallpox, SARS, and foot-and-mouth disease epidemics.
Conclusions
We conclude that (1) there is no general predictive formula for the proportion to be traced as there is for the proportion to be vaccinated; (2) variability in time to detection is favourable for effective tracing; (3) tracing effectiveness need not be sensitive to the duration of the latent period and tracing delays; (4) iterative tracing primarily improves effectiveness when single-step tracing is on the brink of being effective.
The research of D.K. is supported by the EU Sixth Framework Programme for research for policy support (contract SP22-CT-2004-511066). The research of H.H. is supported by the Netherlands Organisation for Scientific Research (NWO grant 918.56.620).Introduction
Control of epidemics of (emerging) infectious diseases, such as SARS, pandemic influenza, or foot-and-mouth disease, always faces the difficulty that some infectives are not yet observed. By concentrating control measures only on observed cases (treatment, isolation, culling), resources are used efficiently but control is often not effective enough. On the other hand, by directing control to the whole population (mass vaccination, prophylactic treatment, preventive culling), epidemics are most likely contained, but at high cost.
Contact tracing of symptomatic infecteds works on an intermediate level: treatment or quarantine of contactees (by contact we mean the possible transmission event, by contactee the individual that is contacted) may be effective because unidentified infecteds are most likely to be found among contactees, and efficient because the resources can be directed towards individuals at risk only. Tracing and quarantine has been successfully applied for smallpox control, where the term ‘Leicester method’ refers to exactly this strategy, with the establishment of specific smallpox hospitals [1]. It was also successful during the more recent SARS epidemic [2], but not during the British foot-and-mouth epidemic [3]. Here we wish to investigate how these differences can be accounted for by the tracing process, thereby distinguishing tracing and quarantine as having separate effects which can be explained as follows.
Symptoms and detection divide the reproduction ratio R (the average number of secondary infecteds per primary infected in a susceptible population) into a part before detection and isolation and a part after [4]:R = R_{0}^{pre}+c (R_{0}−R_{0}^{pre}),where c is the reduction factor due to isolation, R_{0} is the basic reproduction ratio, when no efforts are made to isolate, and R_{0}^{pre} is that part of R_{0} occurring prior to detection. Contact tracing will lead to earlier prevention of transmission due to quarantine of traced infecteds, thereby reducing the uncontrolled transmission R_{0}^{pre} to R^{pre}, and the reproduction ratio R to R_{q}:R_{q} = R^{pre}+c(R_{0}−R^{pre}).The effect of decreasing R_{0}^{pre} to R^{pre} is distinct from the effect of isolation and quarantine which reduce c. For tracing to be effective, R_{q} should be smaller than 1, so R^{pre} should be smaller than (1−c)/(1−cR_{0}), a threshold value determined by c and R_{0}[4]–[6]. For comparison of contact tracing in different situations, we will choose c = 0 (as in [6]–[10]), because it is easy and as arbitrary as any other value. The value c = 0 is valid for foot-and-mouth disease where traced herds are culled.
Several tracing studies have been published, some focussing on specific infections [10], [11] and some with the objective to obtain more general insight, e.g. in a general Markov-type SIR model [7], in a model incorporating symptom development [4], or in simulation models with specific network contact structures [8], [9]. Different assumptions were made with regard to the possibility of only tracing contacts of symptomatic infecteds (single-step tracing [4], [10], [11]), or of iteratively tracing the contacts of traced infecteds (iterative tracing [6]–[9]). However, no general insight has yet been obtained in reducing R^{pre} in relation to time-related characteristics of the infection and the tracing process, such as the latent and infectious periods, the time of symptom-based detection, and delays in the tracing process.
In this paper we consider detection and isolation as an autonomous process most likely governed by detection via clinical symptoms. This allows us to concentrate on contact tracing and quarantine, which are initiated by this autonomous process. Also we only regard transmission prior to control: we will consider tracing effective if R^{pre}<1.
The Model
We study an epidemic in its initial phase and use a branching process for its description, which means that the epidemic can be regarded as a growing tree. Each branch (contact) connects two nodes (contactees) with a hierarchical relation, the infector having infected the infectee. If one of the contactees is notified as being infected, by becoming symptomatic or by being traced via a secondary contact, the contact is traced with a probability p_{c} (all parameters and functions of the model are given in Table 1). Each contact may be traced only once, either from infector to infectee (forwards tracing) or from infectee to infector (backwards tracing), so its traceability can be determined at the time of transmission. Thus, infection trees with traceable and untraceable contacts arise (Figure 1).
10.1371/journal.pone.0000012.g001Single-step and iterative tracing on an epidemic tree, developing from left to right.
Nodes are infecteds, lines are contacts, contactees that were not infected are not represented on the tree. Grey infecteds are asymptomatic, white infecteds are symptomatic, infecteds with a thick border are isolated or quarantined.
Solid lines are traceable contacts, dotted lines are untraceable contacts.
A. Single-step tracing.
In A1a-c, a symptomatic infected is isolated and his traceable contactees are quarantined.
In A2a-b (some time later), one of the quarantined infecteds got symptomatic and his traceable contactees are quarantined.
B. Iterative tracing.
In B1a-c, a symptomatic infected is isolated and all infecteds directly or indirectly linked to this infected by traceable contacts are quarantined.
All quarantined infecteds form a traceable cluster.
10.1371/journal.pone.0000012.t001List of most important parameters.
parameter
description
p_{c}
probability of a contact being traced
δ
duration of the time delay in each tracing step along the epidemic tree
R, R_{q}, R_{0}
reproduction ratio: expected number of secondary infecteds produced by a single infected in a susceptible population; R_{0} is used in the absence of control measures; R_{q} is with tracing and quarantine in effect
R_{0}^{pre}
pre-detection part of R_{0}: expected number of secondary infecteds produced prior to symptom-based detection
τ
time since infection of an infected
τ_{lat}
latent period: time between infection and the start of the infectious period
τ_{inf}
time between infection and the end of the infectious period
α
shape parameter of the detection time distribution (Gamma distribution). If α = 1, the time to detection is exponentially distributed; if α = ∞, it is fixed
p_{c}*
critical tracing probability for which R_{q} = 1
On these infection trees, two types of contact tracing are distinguished. The first type is single-step contact tracing, in which all traceable contactees of a symptomatic case are quarantined, but tracing only continues when an undiscovered infected, quarantined in this way, is detected itself (Figure 1A). The second type is iterative tracing, in which tracing of traceable contactees is continued, and a whole cluster of infecteds that is linked through traceable contacts is quarantined (Figure 1B). Such continuation would be possible if a test were available to determine the infection status of traced contactees.
The underlying model for infection dynamics is based on the framework of [4]. In our model, τ measures time since infection of an individual, which starts with a latent period until τ = τ_{lat} without transmission of the pathogen. During the infectious period, lasting from τ_{lat} to τ_{inf}, infecteds give rise to β new infecteds per unit of time, as long as they are not detected. As we are interested in the effectiveness of tracing only, we do not model possible transmission that might occur while being isolated, so transmission ‘ceases’ after detection of the infected. This leaves us with the basic reproduction ratio prior to detection, defined as the expected number of secondary infections per infected in a susceptible population before detection, R_{0}^{pre}. Here we have adjusted the interpretation of the model in [4], where R_{0} is divided in an asymptomatic and a symptomatic part, by adding a detection delay after becoming symptomatic. The part R_{0}^{pre} is determined by τ_{lat}, τ_{inf}, β, and the distribution of the detection time, a Gamma-distribution with mean 1 (so time is measured relative to the mean detection time) and shape parameter α. Throughout our analyses, β will be scaled accordingly to achieve a desired value of R_{0}^{pre}. This model construction allows a flexible way of exploring different assumptions about the time to detection, the infectious period, and their overlap, and it enables us to evaluate tracing effectiveness for most infections.
For full understanding of contact tracing in our model, we analyzed four special cases regarding the infectious period and the detection time distribution. The infectious period was assumed to be either very short (all transmission occurs instantaneously, so τ_{inf} = τ_{lat}) or very long (of infinite duration, so τ_{inf} = ∞). (In our model, infinite duration can be assumed because each infected will stop spreading the infection after detection. If some infecteds would never be detected, some large τ_{inf}<∞ should be taken.) The detection time was assumed to be either fixed (α = ∞) or highly variable (α = 1, i.e. exponentially distributed). As a control we analyzed intermediate cases (results not shown) and four examples of real infections (influenza, SARS, smallpox, and foot-and-mouth disease), of which the parameter values are listed in Table 2. These parameter values were obtained from literature [2], [3], [12]–[14], assuming that the time to detection consists of the incubation period (time to symptom onset) plus a symptom-to-detection delay, which we assumed to be distributed as observed in the SARS epidemic (average 3.67 days, [2]).
10.1371/journal.pone.0000012.t002Parameter values for the real infections
Infection
R_{0}^{pre}
Detection time1
Infectious period2
References
mean
variance
α
τ_{lat}
τ_{inf}
Influenza
1.5
1.48+3.67
0.22+11.3
2.43
0.097
0.971
12,14
SARS
1.5
3.81+3.67
8.34+11.3
2.94
0.668
3.61
[2]
Smallpox
3
15.5+3.67
4.08+11.3
24.9
0.781
1.35
[13]
FMD (82%)3
1.6
7.83
2.033
303
0.33
∞
[3]
FMD (18%)3
14.93
6.123
43
Sum of incubation period and time from symptom onset to detection as observed during the SARS epidemic [2]
τ_{lat} and τ_{inf} are measured relative to the mean detection time
The detection time distribution of FMD was modeled as a mixed Gamma distribution as in [3]
In our analyses tracing effectiveness will be expressed as the critical tracing probability p_{c}*, defined as the proportion of contacts that need to be traced to achieve R^{pre} = 1. If at least that many contacts are traced, epidemics will certainly die out if transmission during isolation or quarantine is prevented or limited to a small number of health-care workers that do not re-introduce the infection into the general community (see also [5]). If such re-introductions cannot be excluded, a lower R^{pre} may be aimed for. Because of this threshold of 1, the R_{0}^{pre} values (without tracing) were assumed to have some value larger than 1 (otherwise tracing would not be needed at all) and less than published R_{0} values for the specific cases (Table 2). First we study p_{c}* as a function of R_{0}^{pre} and τ_{lat} (with tracing delay δ = 0), and second as a function of τ_{lat} and δ (with R_{0}^{pre} = 1.5).
ResultsSingle-step tracing without delays
In an epidemic where single-step tracing is applied, asymptomatic infecteds can spread the infection until they are detected and isolated. Detected infecteds are asked to name their contactees, and a proportion p_{c} of all contactees will be reported and quarantined. Only when quarantined infecteds are detected themselves, tracing is continued (Figure 1A).
We determined the critical tracing probability p_{c}* as a function of the latent period for three values of R_{0}^{pre} (1.5, 2, and 3). The results are shown in Figure 2. If the detection time is fixed (α = ∞), a too large latent period (larger than the detection time) results in a situation where every infected is detected before transmitting the infection, so tracing need not prevent any transmission; hence the maximum τ_{lat} of 1 in panels 2B and 2D. If we locate the approximate position of real infections (parameter values in Table 2) in Figure 2, we observe that the long infectious period will be the best approximation for most infections, because we only regard pre-detection transmission and infecteds will often still be infectious at the time of detection. This does not entirely hold for influenza, which therefore fits best between the short and long infectious period. Because of the detection period distributions, smallpox and FMD are best described by panel 2D, whereas SARS and influenza fit best into panel 2C (influenza also into 2A).
10.1371/journal.pone.0000012.g002The effectiveness of single-step contact tracing without tracing delays.
Effectiveness is expressed as the minimum proportion of contacts that need to be traced for effective control (critical tracing probability p_{c}*).
The plots show p_{c}* as a function of the latent period relative to the mean time to detection (τ_{lat}).
There are four special cases: A. Short infectious period and variable time to detection; B. Short infectious period and fixed detection time; C. Long infectious period and variable time to detection; and D. Long infectious period and fixed detection time.
The three curves denote p_{c}* for different values of the pre-isolation reproduction ratio R_{0}^{pre}.
Indicated by dashed lines are the average τ_{lat} for four infections, in the panels with closest correspondence to the actual parameter values (Table 2).
Influenza appears in two panels with long and short infectious period, because it corresponds to both parameter sets equally.
In three special cases (panels 2A,B,D), with fixed incubation period and/or short infectious period, the proportion of contacts to be traced is at least 1−1/R_{0}^{pre}. This lower limit 1−1/R_{0}^{pre} is due to forwards tracing, when all infecteds that are traceable via their infector, are quarantined before the end of their latent period. Then, the effective reproduction ratio will be equal to the untraceable proportion of R_{0}^{pre}, i.e. R_{q} = (1−p_{c})R_{0}^{pre}, resulting in p_{c}* = 1−1/R_{0}^{pre}. This is likely to be case with smallpox (panel 2D). However, if the latent period is short, as seen for influenza and possibly FMD, quarantine will occur too late to prevent all infections and more contacts need to be traced.
In the fourth panel (2C), with variable detection time and long infectious period, effective contact tracing requires a proportion of contacts smaller than 1−1/R_{0}^{pre} to be traced, if the latent period is large enough (like SARS). This can be explained as follows. If the detection time is variable and the latent and infectious periods are large, many infecteds will be detected before becoming infectious. If many infecteds are not infectious before being detected, the few that are should be very infectious (“superspreaders”) to attain a given R_{0}^{pre} (the average number of secondary infections before detection). Because many of the infectees of these superspreaders will be detected early (variable detection time), the superspreaders will be quarantined after backwards tracing, which adds to the effect of forwards tracing preventing infectees to reach their infectious period.
Effectiveness may be very sensitive to the latent period, especially if there is little variation in the detection time. This is most apparent in the sharp transition in panel 2B, where tracing is only effective if all infectors are detected (at τ = 1−τ_{lat}) before the infectious period (at τ = τ_{lat}), so τ_{lat}>0.5. The high sensitivity can be a problem for assessing the tracing effectiveness for a specific infection: the conclusion may largely depend on how correct the available estimates for the latent and incubation periods are, the incubation period determining the time to detection. Figure 2 indicates that this might be a problem for FMD.
Single-step tracing with delays
In the previous section tracing was an instantaneous process: detection and isolation of a given case were followed by quarantine of traceable contactees of that case without any delay. The effect of a delay of duration δ is that contactees of detected infecteds can continue transmitting the infection for an extra δ time units, unless they are detected themselves during this interval.
The effect of a delay in each tracing step was studied by determining the critical tracing probability p_{c}* for R_{0}^{pre} = 1.5, as a function of the latent period and tracing delay. Figure 3 shows the resulting contour plots for the four special cases with single-step tracing, and it indicates the approximate positions of the four real infections. Along the τ_{lat} axis (δ = 0) lie the R_{0}^{pre} = 1.5 curves of Figure 2.
10.1371/journal.pone.0000012.g003The effectiveness of single-step contact tracing with tracing delays, with the pre-detection reproduction ratio <italic>R</italic>
<sub>0</sub>
<italic>
<sup>pre</sup>
</italic> = 1.5.
Effectiveness is expressed as the minimum proportion of contacts that need to be traced for effective control (critical tracing probability p_{c}*).
The contour plots show p_{c}* as a function of the tracing delay δ and the latent period τ_{lat}, measured relative to the mean detection time, for four special cases: A. Short infectious period and variable incubation period; B. Short infectious period and fixed incubation period; C. Long infectious period and variable incubation period; and D. Long infectious period and fixed incubation period.
Dark grey shadows indicate areas where tracing is ineffective, light grey shadows indicate areas where p_{c}* = 0.33.
Indicated by dashed lines are the average τ_{lat} for four infections, in the panels with closest correspondence to the actual parameter values (Table 2).
Influenza appears in two panels with long and short infectious period, because it corresponds to both parameter sets equally.
It appears that the iso-p_{c}* contours in Figure 3 are linear in three of the four special cases (panel 3A,B,D), and approximately linear in the fourth case, with long infectious period and variable incubation time (panel 3C: τ_{inf} = ∞; α = 1). The slopes of all (approximate) lines are always between 0.5 and 1.5. This means that the effect of tracing delays is comparable to the sensitivity to the latent period as observed in Figure 2, so plots of the critical tracing probability as a function of the delay will resemble the plots in Figure 2, only mirrored (as in Figure 4).
10.1371/journal.pone.0000012.g004The effectiveness of single-step and iterative contact tracing for control of influenza, smallpox, SARS, and foot-and-mouth disease.
Effectiveness is expressed as the minimum proportion of contacts that need to be traced for effective control (critical tracing probability p_{c}*); p_{c}* is plotted as a function of the relative delay (δ, proportion of the incubation period) or the absolute delay (days).
In the contour plots in Figure 3, tracing is ineffective in the dark grey areas (small latent period or large delay), so smallpox and SARS control are predicted to be able to cope with some delays, whereas it might be more of a problem with FMD and influenza. If forwards tracing is maximally effective and not affected by the delay, p_{c}* = 1−1/R_{0}^{pre} = 0.33 as indicated by the light grey area (smallpox if δ<0.5). If the detection time is variable, tracing may be effective even if a proportion less than 1−1/R_{0}^{pre} is traced (SARS if δ<0.4, panel 3C). Because backwards tracing is needed to attain this result, it is only possible when the tracing delay is shorter than the infectious period (not shown in the Figures).
The contour plots show that for some combinations of latent period and delay, the sensitivity to small changes in the delay is large (contours lie close to one another). Then, small changes in the delay can have a major impact when the time of quarantine is shifted from just before to just after the end of the latent period, similar to the sensitivity to τ_{lat} as observed in Figure 2. This effect is most dramatic if there is little variation in the detection time (smallpox: δ≈0.7, and FMD: δ≈0.2).
Figure 4 shows the effectiveness of single-step tracing for SARS, smallpox, influenza, and FMD using the parameter values listed in Table 2. Panel 4A shows the relation between p_{c}* and the relative tracing delay δ, from which it appears that the general cases shown in Figures 2 and [3] are good predictors for assessing tracing effectiveness. Influenza can be placed somewhere between long and short infectious periods (panels 2A,C, and 3A,C): it is hardly effective without delay, and ineffective already if δ = 0.1. SARS control requires a tracing probability p_{c}*<1−1/R_{0}^{pre} which is relatively insensitive to delays. Smallpox requires p_{c}* = 1−1/R_{0}^{pre}, and is insensitive to delays up to some point (δ≈0.6) where tracing becomes quickly ineffective. Finally, FMD can be traced effectively only if δ is small, but it is sensitive to delays already if δ is small and requires a tracing probability p_{c}*>1−1/R_{0}^{pre}. Measured in real time (panel 4B), effectiveness of tracing appears to be highly dependent on the actual generation time of the infection. Influenza control is hardly possible, FMD control will be difficult, whereas tracing is likely to be successful for smallpox and SARS.
Iterative tracing
Suppose a test is available to determine the infection status of traced contactees. We can then continue tracing iteratively before infected contactees are detected due to symptoms, until no further infecteds are found (Figure 1B).
We evaluated iterative tracing with and without delays for the same cases as single-step tracing, resulting in Figures similar to Figures 2 and [3] (see Supporting Information). As expected, the universal effect of iterative tracing is to lower p_{c}*, although the lower limits 1−1/R_{0}^{pre} as observed in three of the four special cases remain intact. The largest difference between single-step and iterative tracing is observed when single-step tracing is on the brink of being effective, so that is the only situation where iterative tracing will make a difference (see also Figure 4). Without delays this difference is rather intuitive: single-step tracing is not effective if the latent period is small, but iterative tracing will always be effective, because if all contacts are traceable, the first detection is immediately followed by quarantine of every infected. This is the case with influenza, although the improvement will probably not be sufficient to making tracing effective for influenza control. On the other hand, when the latent period is large and single-step tracing is effective, iterative tracing only improves on this if there are significant delays, providing one or two more days to reach the required p_{c}* (smallpox, SARS, and FMD).
Discussion
We studied the effectiveness of contact tracing in a model for the start of an epidemic, that is flexible enough to use for most infections. Effectiveness of contact tracing was expressed as the minimum proportion of contacts that need to be traced to obtain a reproduction ratio prior to control of 1 (p_{c}*). Other threshold values may be chosen if more is known about transmission of the infection to the general community while infectious individuals are being isolated or quarantined. For instance, if isolated infecteds still cause an average of 0.3 new cases in the general community after isolation or quarantine, p_{c}* should be redefined such that R^{pre} = 0.7.
The first conclusion from our model is that, for a given R_{0}^{pre}, the critical tracing probability can take any value depending on all infectious disease characteristics in the model: the latent period, the infectious period, and the detection time distribution. In contrast to some earlier publications on contact tracing [7], [9], there exists no general expression for p_{c}* as there is for the proportion to vaccinate for effective control in a well-mixed population. For smallpox the relation p_{c}* = 1−1/R_{0}^{pre} holds reasonably well, but for SARS it is smaller, and for influenza and FMD it is larger.
The second conclusion is that a variable detection time improves tracing effectiveness, possibly even resulting in p_{c}*<1−1/R_{0}^{pre}. This does not mean, of course, that one should aim at late detection of some infecteds (which would increase R_{0}^{pre}), but that apparent variability is an argument to use tracing. The reason is that the few infecteds that are detected late (or not at all, which is essentially the same) will be discovered by backwards tracing which is an additional effect to forwards tracing that in itself may already result in p_{c}* = 1−1/R_{0}^{pre}. It was earlier found that p_{c}* can be smaller than 1−1/R_{0}^{pre}[9], but by a different mechanism, namely the presence of shared contacts in a network.
The third conclusion is that the sensitivity of tracing effectiveness to the latent period and tracing delay may be large, especially in the case of single-step tracing. If this is the case already with a small delay (influenza and to a lesser extent FMD), reliability of parameter estimates will be crucial to establish whether tracing is advantageous. If it is only the case with larger delays (smallpox, SARS), tracing may be effective as long as it is done quick enough.
The fourth conclusion is that in most situations single-step and iterative tracing are almost equally effective. A considerable difference can only be expected if single-step tracing is not or hardly effective, which is also when the sensitivity to the latent period and tracing delay is largest (see above). Thus, influenza control will benefit from iterative tracing already without delays, whereas with the other infections one or two days may be gained (Figure 4C,D), thus having more time to achieve the required p_{c}*. In the real world, the choice between single-step and iterative tracing will be based on what is possible. Capacity problems may reduce the effectiveness of iterative tracing if effort is directed towards secondary contactees prior to primary contactees; on the other hand, if not quarantine but vaccination is applied, it might be worthwhile to traced contactees of contactees even without diagnostic tests [6].
For determining the cost-effectiveness of contact tracing, not only the critical tracing probability p_{c}* is required, but also how easily that p_{c}* can be achieved. A key aspect will be the possibility to distinguish relevant contacts: for control of sexually transmitted diseases relevant contacts are easily identified (that is, if people are willing to co-operate), but in case of respiratory pathogens like SARS or influenza the notion of relevant contacts is rather diffuse. Good insight will be obtained by not only regarding the sensitivity of tracing (p_{c}), but also the positive predictive value: what proportion of traced individuals is actually infected, and if we put more effort into increasing p_{c}, how will it decrease the positive predictive value?
A relatively easy way to increase the proportion of traced contacts for respiratory pathogens like influenza might be quarantine of households. This can be effective if households are the prime location for spread of the infection; specific models taking into account the contact structures within and between households will be better-suited to study this strategy.
In our model we only regarded transmission of the infection before tracing or isolation, allowing us to focus on the characteristics of the contact tracing itself. Regarding the control of animal disease epidemics, e.g. foot-and-mouth disease or avian influenza, it is not unreasonable to assume that there will be no other transmission: all infected and suspected farms will be culled and they really will stop being infectious. However, with human infections it is very likely that isolation and quarantine will be less effective. This results in a complicated situation, because transmission will not only be reduced after detection, the contact structure will also change, with only a limited number of people (health-care workers) having contact with multiple isolated patients. For the case of SARS, a more detailed model on exactly this aspect was studied by Lloyd-Smith et al [5], who indeed concluded that minimizing within-hospital transmission might be crucial, especially if tracing and isolation occur rather slowly. Further studies are required for this complex problem. The present paper provides a good understanding of the principles of contact tracing, and how infectious disease characteristics determine the effectiveness.
Methods
The effect of single-step tracing can be measured by the effective reproduction ratio R_{q}, the average number of new infections per infected. Whereas the rate of being traced backwards at time τ since infection (due to detection of any infectee) is equal for all infecteds, the rate of being traced forwards (due to detection of the infector) depends on the generation-time distribution, which in turn depends on the number of traceable generations backwards in the transmission tree (traceable ancestors). Hence, infectives need to be typed according to the number of traceable ancestors j, and R_{q} is the largest eigenvalue of the next-generation matrix K with entries k_{ij} being the expected number of type-i infecteds (with i traceable ancestors) produced per type-j infected [15]. Because type-j infecteds only produce type-0 and type-j+1 infecteds, by untraceable and traceable contacts respectively, all entries other than k_{0j
} and k_{j}_{+1,j
} are equal to 0.
Even though the matrix is infinitely large, we conjecture from all observed numerical calculations that the positive elements k_{j}_{+1,j
} converge in the sense that |k_{j}_{+1,j
}−k_{j}_{+2,j+1}|→0 as j→∞. Therefore, for numerical evaluation we ‘closed’ the matrix to an (n+1)^{2} matrix with k_{nn} = k_{n}_{+1,n
}. For the four special cases, R_{q} could be calculated numerically in Mathematica®, but for other cases considered, the entries of the matrix had to be determined by numerical simulation. Analysis of single-step tracing with delay could be done by calculating R from the next-generation matrix, as for the model without delay (see Supporting Information for all details).
With iterative tracing, isolation of a single infected results in quarantine of a cluster of infecteds, all mutually linked by traceable contacts (Figure 1B). By recognizing these clusters of infectives, an epidemic of infected individuals can be regarded as an epidemic of traceable clusters. Each untraceable contact in the transmission tree results in infection of a new cluster index case, so the average number of untraceable contact infections caused by a single cluster, the cluster reproduction ratio R_{c}, determines the effectiveness of iterative contact tracing. If Y(p_{c}) denotes the expected cumulative infectiousness of a cluster at the time of cluster quarantine, then
As with single-step tracing, effectiveness of iterative tracing is determined by considering the critical tracing probability (for achieving R_{c} = 1):If p_{c} = 0, then each infective will be a separate cluster, and R_{c} = Y(0) = R_{0}^{pre} by definition. For p_{c}>0 the cluster size will be larger than one, but cluster infectiousness need not be larger, as backwards tracing can reduce the infectious period of superspreaders with large detection time as it does in single-step tracing.
For two of the four special cases (α = ∞), and for the case considered by Müller et al [7] (α = 1, τ_{inf} = ∞, τ_{lat} = 0), Y(p_{c}) can be calculated numerically in Mathematica®. For all other cases, including the real infections, stochastic simulations were needed.
Incorporating delays into iterative tracing does not generally change the concept of R_{c}, but it makes the calculation of Y(p_{c}) very complicated, because there is no longer a single time of cluster quarantine. Two complications arise: first, delays may cause contactees in the chain emanating from a symptomatic infected to become symptomatic themselves before the tracing process reaches them. They then initiate a new tracing process of their own within the same cluster. Second, if the infection process is faster than the tracing process, the cluster size grows infinitely large and iterative tracing becomes ineffective. We determined p_{c}* by stochastic simulation of clusters until quarantine of the final infected (see Supporting Information for all details).
Supporting Information
Supporting information for the paper “The Effectiveness of Contact Tracing in Emerging Epidemics”
(0.60 MB PDF)
We thank one anonymous reviewer for his or her significant contribution to the present paper.
ReferencesFennerFHendersonDAAritaIJežekZLadnyiID1988DonnellyCAGhaniACLeungGMHedleyAJFraserC2003Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong.FergusonNMDonnellyCAAndersonRM2001Transmission intensity and impact of control policies on the foot and mouth epidemic in Great Britain.FraserCRileySAndersonRMFergusonNM2004Factors that make an infectious disease outbreak controllable.Lloyd-SmithJOGalvaniAPGetzWM2003Curtailing transmission of severe acute respiratory syndrome within a community and its hospital.PorcoTHolbrookKFernyakSPortnoyDReiterR2004Logistics of community smallpox control through contact tracing and ring vaccination: a stochastic network model.MüllerJKretzschmarMDietzK2000Contact tracing in stochastic and deterministic epidemic models.KissIZGreenDMKaoRR2005Disease contact tracing in random and clustered networks.EamesKTDKeelingMJ2003Contact tracing and disease control.EichnerM2003Case isolation and contact tracing can prevent the spread of smallpox.BeckerNGGlassKLiZAldisGK2005Controlling emerging infectious diseases like SARS.Group WHOW2006Nonpharmaceutical interventions for pandemic influenza, international measures.AldisGKRobertsMG2005An integral equation model for the control of a smallpox outbreak.FergusonNMCummingsDATCauchemezSFraserCRileyS2005Strategies for containing an emerging influenza pandemic in southeast Asia.DiekmannOHeesterbeekJAP2000Mathematical epidemiology of infectious diseases.