Abstract
The test-trace-isolate-quarantine (TTIQ) strategy is used to break chains of transmission during a disease outbreak. Confirmed-positive pathogen carriers are isolated from the community to prevent onward transmission and their recent close contacts are identified and pre-emptively quarantined. TTIQ, along with mask wearing and social distancing, make up the non-pharmaceutical interventions that are utilised to suppress the ongoing SARS-CoV-2 pandemic. The efficacy of the TTIQ strategy depends on the probability of isolating a case, the fraction of contacts quarantined, and the delays in these processes. Here we use empirical distributions of the timing of SARS-CoV-2 transmission to quantify how these parameters individually contribute to the reduction of onwards infection. We show that finding and isolating index cases, and doing so with minimal delay after symptom onset, have the largest effects on case reduction, and that contact tracing can make up for deficiencies in testing coverage and delays. These results can be used to assess how TTIQ can be improved and optimised. We provide an online application to assess the efficacy as a function of these parameters.
Code is publicly available at http://github.com/ashcroftp.
1 Introduction
Individuals who are confirmed as infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogen are isolated from the population to prevent further transmission. The individuals who have been in recent close contact with an infected individual have an increased risk of being infected themselves. By identifying the potentially-infected contacts through contact tracing, and eventually quarantining them, transmission chains can be broken. Thus contact tracing is an essential public health tool for controlling epidemics (WHO, 2020). The strategy of testing to identify infected cases, isolating them to prevent further transmission, and tracing & quarantining their recent close contacts is known as test-trace-isolate-quarantine (TTIQ) (Salathé et al., 2020). This strategy is a fundamental non-pharmaceutical intervention which is used globally to control the ongoing SARS-CoV-2 pandemic (Kucharski et al., 2020).
Testing typically occurs once an individual develops symptoms of coronavirus disease 2019 (COVID-19). As presymptomatic transmission makes up approximately 40% of total onward transmission (He et al., 2020; Ashcroft et al., 2020a; Ferretti et al., 2020a), it would be possible for the number of secondary cases to be more than halved if infected individuals are isolated from the community at the time of symptom onset. However, as testing follows from symptoms, the testing & isolating strategy without subsequent contact tracing & quarantine is unlikely to capture asymptomatic cases which make up 20% of cases (Buitrago-Garcia et al., 2020), and thus isolating 100% of cases would not be possible.
Contact tracing & quarantine have the potential to be effective interventions against the spread of COVID-19 because of the high frequency of pre-symptomatic or asymptomatic transmission from recently-infected individuals (Moghadas et al., 2020). Potentially-infected contacts can be identified and quarantined before they would be isolated as a result of developing symptoms and/or receiving a positive test result, such that their onward transmission is reduced. This is exemplified in the light of the high dispersion of the offspring distribution and frequency of super-spreader events (Riou & Althaus, 2020; Endo et al., 2020; Adam et al., 2020), where large numbers of potentially-infected contacts can be quarantined to prevent widespread community transmission. Tracing & quarantine does not depend on symptom development. Hence, this strategy is capable of reducing onward transmission even from asymptomatically infected individuals.
TTIQ strategies are not perfect: each stage in the process is subject to delays and uncertainties and it would be impossible to prevent all onward transmission through TTIQ alone (Ferretti et al., 2020b; Kucharski et al., 2020; Kretzschmar et al., 2020; Quilty et al., 2020; Ashcroft et al., 2020b). Furthermore, in the presence of widespread community transmission the contact tracers may be overwhelmed by the volume of cases. In this scenario it is important to optimise the resources (i.e. the person hours of the contact tracers) to minimise onward transmission.
In a previous study of TTIQ efficacy, Ferretti et al. (2020b) used an approach based on the empirically-observed timing of transmission events – but with sub-stantial approximations around the TTIQ process – to get to an analytically tractable prediction of the impact of TTIQ on SARS-CoV-2 transmission. They concluded that widespread digital contact tracing (with minimal delay between index case identification and quarantine of secondary cases) would be necessary to reduce the effective reproduction number, Re, below one to bring an outbreak under control. Kucharski et al. (2020) used an agent-based model with detailed contact structures to simulate intervention strategies. While the TTIQ process is more accurately described than in Ferretti et al. (2020b), they did not use empirical data about the timing of transmission, which is crucial for quantifying the impact of isolation and quarantine. Kretzschmar et al. (2020) opted for a discrete-time branching process model of transmission and TTIQ. While they explicitly accounted for the timing of infection events and accurately described the TTIQ process, they predominantly focussed on assessing the role of digital contact tracing based on mobile applications. In this paper we develop an analytical approach which builds on our previous work in which we have quantified the impact of quarantine duration and high-lighted the optimal use of test-and-release strategies (Ashcroft et al., 2020b). Briefly, we use the empirically-observed distributions of transmission timing [Fig. 2; Ferretti et al. (2020a)] to determine when infections occur (Fig. 1). We then introduce five parameters to describe the TTIQ process: i) f, the probability that an index case is isolated from the population and is interviewed by contact tracers; ii) Δ1, the time delay between symptom onset and isolation of the index case; iii) τ, the duration prior to symptom onset in which contacts are identifiable; iv) g, the fraction of identifiable contacts that are quarantined; and v) Δ2, the delay between isolation of the index case and the start of quarantine for the contacts. We compute the expected number of tertiary cases per index case under the TTIQ interventions, with the aim being to reduce this number below one to suppress the growth of the epidemic (see Methods for details). We systematically explore this parameter space, first for the “testing & isolation” intervention in the absence of contact tracing (Fig. 1A), and then with additional “tracing & quarantine” (Fig. 1B).
2 Methods
Our primary goal is to quantify the reduction of transmission by isolating individuals who test positive for SARS-CoV-2 and by quarantining their recent close contacts with an increased risk of infection. We refer to the initial confirmed case as the index case, and the infected contacts as secondary cases. We know that the index case developed symptoms at time , but the time at which they were infected, t1, is generally unknown. Secondary cases will be infected by the index case at some time t2 (t2 > t1), and develop symptoms at time (Fig. 2A).
2.1 Generation times, infectivity profiles, and incubation periods
The relationships between the times are determined by: the generation time distribution, q(t2 − t1|θq), describing the time interval between the infection of an index case and secondary case (Fig. 2B); the infectivity profile, , describing the time interval between the onset of symptoms in the index case and infection of the secondary case (Fig. 2C); and the incubation period distribution, , describing the time between the infection of an individual and the onset of their symptoms (Fig. 2D). For these distributions, we use empirical estimates from Ferretti et al. (2020a) which are based on a large set of transmission pairs and minimal assumptions about the relationship between infectiousness and symptoms, which would otherwise bias the resulting generation time distribution (Lehtinen et al., 2020).
2.2 Quantifying the number of secondary cases
Consider an index case who develops symptoms of COVID-19 at time . The time of infection, , is generally unknown. Without any TTIQ intervention this individual would contact and infect k1 individuals during the course of the infection, where this number of contacts is distributed as across individuals in the population. Note that this number of contacts depends on the current level of other non-pharmaceutical interventions, such as mask wearing and social distancing. The number of secondary infections up to a time T after developing symptoms would then be where p(t|θp) is the infectivity profile and is the cumulative infectivity profile.
Index cases who develop symptoms and/or test positive for SARS-CoV-2 should be isolated from the population. This occurs in a fraction f of index cases who are isolated at a time , where Δ1 > 0 is the delay between symptom onset and isolation. The parameter Δ1 can be interpreted as the delay of taking a test after symptom onset, waiting for the result, and entering isolation, or alternatively as the delay between symptom onset and self-isolation. The remaining 1 − f index cases are not isolated (T → ∞). We can compute the number of secondary infections, n2, as a function of testing coverage f and delay Δ1, as shown in Fig. S1. For a given symptom onset time and degree k1 of the index case, we have Averaging over k1, which is distributed as , and keeping fixed as the reference time point, we arrive at where Re = ⟨k1⟩ is the mean of , i.e. the average number of secondary infections in the absence of testing & isolation (f = 0).
2.3 Quantifying the number of tertiary cases
Each secondary case has some potential to cause further infections, which will be the tertiary cases of the index case. The number of tertiary infections caused by a secondary case who is infected at t2 and isolated at time T, will be where k2 is the number of contacts of the secondary case, t3 is the infection time of the tertiary cases, q(t|θq) is the generation time distribution, and is the cumulative generation time distribution. Note that we use the generation time distribution here, as our reference point is the time of infection (t2), whereas in Eq. (3) the reference point was the time of symptom onset .
Under TTIQ interventions, the index and secondary cases can be isolated following a positive test result and/or self-isolation after symptom onset. If an index case is confirmed positive, then contact tracing can be used to identify and quarantine individuals who have recently been exposed to the confirmed case. Quarantining these individuals prevents the onward infection of tertiary cases (Fig. 1B). We introduce three further parameters to quantify contact tracing: i) τ > 0, the duration of lookback prior to symptom onset of the index case; ii) 0 ≤ g ≤ 1, the probability to identify and quarantine a secondary contact that was infected within the contact tracing window; and iii) Δ2 > 0, the delay between isolating the index case and quarantining the identified secondary contacts.
There are many permutations of events that contribute to the number of tertiary cases under TTIQ, as shown in Fig. S2. The index case may not be detected (1 − f), and hence contact tracing is not possible. If the index case is detected (f), then a fraction g of the secondary cases that were infected within the contact tracing window are quarantined at time (Fig. 1B). The remaining fraction 1 − g, as well as the secondary cases that were infected outside of the contact tracing window , are not quarantined. However, the non-traced contacts may themselves be tested and become index cases that are isolated at time , where is the symptom onset time of the secondary case. By considering these different scenarios, we arrive at an expression for the number of tertiary cases per index case under TTIQ, We now have to average Eq. (5) over , and k2 to obtain the expected number of tertiary cases per index case under TTIQ. We first note that for incubation period γ ≥ 0. Hence we can write where g(γ) is the incubation period distribution. Note that we have assumed the independence between symptom onset and infectivity, which may lead to an overestimation of the fraction of tertiary cases prevented. Keeping fixed as the reference time point, averaging Eq. (5) over , and k2 gives the expected number of tertiary cases per index case under TTIQ: where we have substituted such that Eq. (7) can be further simplified to Finally, in the absence of contact tracing (g = 0), the number of tertiary cases under testing & isolation only is given by
2.4 Confidence intervals
The primary sources of uncertainty in the outcomes of this model come from the generation time distribution and infectivity profile, which are inferred from empirical serial interval distributions (Ferretti et al., 2020a). Following Ferretti et al. (2020a), we use a likelihood ratio test to extract sample parameter sets for each distribution that lie within the 95% confidence interval.
Concretely, we first identify the maximum likelihood parameter sets and for the infectivity profile and generation time distribution, respectively. We then randomly sample the parameter space of each distribution, and keep 1,000 parameter sets whose likelihood satisfies , where λn is the 95% quantile of a χ2 distribution with n degrees of freedom. The infectivity profile is described a shifted Student’s t-distribution, which has n = 3 parameters, while the generation time is described by a Weibull distribution with n = 2 parameters.
We then use these sampled parameter sets to generate the number of secondary and tertiary cases, and the extrema of cases across all of these parameter sets determines the 95% confidence interval for the number of cases. For the estimate of the number of secondary cases under testing & isolation [Eq. (3)], we only have to consider the uncertainty of the parameters of the infectivity profile θp. Under the full TTIQ strategy, we need to use estimates of both θp and θq. We assume parameter independence, and keep all (θp, θq) combinations whose joint likelihood satisfies .
2.5 Interactive app
To complement the results in this manuscript, and to allow readers to investigate different TTIQ parameter settings, we have developed an online interactive application. This can be found on the CH Covid-19 Dashboard at https://ibz-shiny.ethz.ch/covidDashboard/.
3 Results
3.1 Reducing cases by testing & isolating
The efficacy of testing & isolating is determined by two parameters: the probability f to find and isolate an infected individual; and the time delay Δ1 between symptom onset and isolation of the index case. The expected number of secondary or tertiary cases [Eqs. (3) & (10)] is also dependent on the current intensity of the epidemic, Re, which is the expected number of secondary cases per infected in the absence of testing & isolating (f = 0). This effective reproduction number depends on the current suppression measures against SARS-CoV-2 transmission (social distancing, mask wearing, home office, etc.), as well as seasonality and levels of immunity/vaccination.
Epidemics can be controlled by testing & isolating if this intervention reduces the expected number of secondary or tertiary cases per index case to below one. We here focus on the number of tertiary cases, but results for the number of secondary cases are qualitatively equivalent (Fig. S3).
The region of (f, Δ1) parameter space in which the number of tertiary cases is less than one, i.e. the region in which the epidemic is controlled by testing & isolating, is shrinking for higher Re epidemics (Fig. 3A). Higher testing & isolation coverage (f) or shortened delays between symptom onset and isolation (Δ1) are required to control SARS-CoV-2 outbreaks as Re increases. Increasing the fraction of infecteds that are isolated buys more time to isolate them, but with diminishing returns.
A SARS-CoV-2 outbreak with Re = 1.1 can be controlled by isolating as few as 18% [95% confidence interval (CI): 15%,22%] of infected cases at the time of symptom onset (Δ1 = 0 days) (Fig. 3B). If the infected index or secondary cases wait Δ1 = 2 days after symptom onset before isolating (i.e. they wait for a test result), then 39% [CI: 28%,60%] of infecteds would have to be isolated for the epidemic to be controlled. Isolating after Δ1 = 4 days would be insufficient to control the epidemic even if all cases were isolated [CI: 65%,n.a.]. For faster-spreading SARSCoV-2 outbreaks (Re = 1.5), we would require 65% [CI: 55%,81%] of infecteds to be isolated immediately after they develop symptoms (Δ1 = 0 days) to control the epidemic. With a delay Δ1 ≥ 2 days, testing & isolating would be insufficient to control the epidemic even if 100% of infecteds are isolated. We note that the frequency of asymptomatic cases (20%; Buitrago-Garcia et al. (2020)) means that we would not be able to isolate 100% of infecteds if we wait for symptoms to develop.
3.2 Reducing cases by additional contact tracing & quarantine
The efficacy of tracing & quarantine is determined by three further parameters: the duration of the contact tracing window prior to symptom onset in the index case τ; the probability to identify and quarantine a secondary case that was infected by an index case within the contact tracing window g; and the delay between isolating the index case and quarantining the secondary cases Δ2. The expected number of tertiary cases [Eq. (7)] is also dependent on the intensity of the epidemic in the absence of TTIQ, Re, as well as the probability (f) and delay (Δ1) of finding and isolating an index case.
The impact that contact tracing has on epidemic control can be seen by varying the parameter g. For g = 0, no contacts are traced & quarantined, and hence we return to the testing & isolation strategy (Fig. 3). By increasing g, we expand the parameter space in which n3 < 1 (Fig. 4), i.e. contract tracing allows an epidemic to be controlled for lower fractions of index cases found (f) and/or longer delays to isolating the index case after they develop symptoms (Δ1).
To visualise the impact of each parameter on the number of tertiary cases, we consider focal parameter sets for the five TTIQ parameters, (f, g, Δ1, Δ2, τ). We then calculate the expected number of tertiary cases when we perturb each single parameter, keeping the remaining four parameters fixed (Fig. 5).
Modifying the fraction of index cases that are identified and isolated (f) has the largest effect of all parameter changes. By identifying more index cases (increasing f), we not only prevent the onward transmission to new secondary cases through isolation, but we also allow infected contacts to be traced and quarantined.
Increasing the fraction of secondary cases that are quarantined (g) has a smaller return than increasing f. If only 30% of index cases are identified, then increasing g results in a small reduction of the number of tertiary cases and for Re = 1.5 the epidemic cannot be controlled even if all secondary cases (g = 1) of known index cases are quarantined (Figs. 5A & B). However, if a large fraction of index cases are identified (f = 0.7), then increasing g can control an epidemic that would be out of control in the absence of contact tracing (Figs. 5C & D).
After increasing f, the next most effective control strategy is to reduce the delay between symptom onset and isolation of the index case (Δ1). Reducing the time taken to quarantine secondary cases has a lesser effect on the total number of tertiary cases. Finally, looking back further while contact tracing (increasing τ) allows more secondary cases to be traced & quarantined. However, this does not translate into a substantial reduction in the number of tertiary cases as the extra cases which are traced have already been infectious for a long time, and will thus have less remaining infectivity potential. Hence increasing τ comes with diminishing returns.
To check the robustness of these effects across all parameter combinations (not just varying a single parameter), we performed uniform parameter sampling and used linear discriminant analysis (LDA) to capture the impact that each parameter has on the number of tertiary cases (Fig. 6). We find that f is the dominant parameter to determine the number of tertiary cases, followed by Δ1, g, Δ2, and finally τ has the smallest impact (Fig. 6B).
4 Discussion
By combining empirically well-supported estimates of the infection timing of SARSCoV-2 with a simple model of transmission, we have calculated the impact of test-trace-isolate-quarantine (TTIQ) interventions against the spread of COVID-19. Overall, we find that TTIQ has the potential to control epidemics with an Re of up to 1.5. This would be practically infeasible under testing & isolation alone, which would require 65% of positive cases to isolate immediately after the time of symptom onset to be effective. By increasing the fraction of contacts that are identified and quarantined, we can successfully suppress an epidemic even if fewer index cases are isolated or if isolation is delayed by up to 2 days. Based on this analysis, we find that the greatest impact comes from increased identification of index cases and reduction of delay between symptom onset and isolation. These parameters have a compound effect on overall transmission as they contribute to the direct reduction of onward infection from an index case, and they allow more contacts to be traced earlier.
Increasing the duration of the contact tracing window by looking back further in time has limited return under our model of forward contact tracing (identifying who is infected by the index case). However, if we were interested in identifying the source of infection (backwards contact tracing), then increasing the duration of the contact tracing window could lead to the identification of transmission clusters.
When comparing to the findings of Ferretti et al. (2020b), we find that contact tracing has less impact on epidemic suppression, and that the speed of contact tracing is of secondary importance to the speed of isolating index cases. This difference can be attributed to Ferretti et al. (2020b)’s approach to model contact tracing and isolation as independent events (i.e. tracing an index cases’ contacts says nothing about whether the index case has been isolated), which leads to an overestimation of contact tracing’s impact (Fraser et al., 2004).
In Kretzschmar et al. (2020) – this time with contact tracing dependent on testing & isolation – they concluded that reducing the delay to isolation after symptom onset has the greatest impact on TTIQ effectiveness. This conclusion was made without systematic analysis of all parameters, and we now find that changing testing & isolation coverage has a greater effect on the number of tertiary cases.
Our approach and results are crucially dependent on the distribution of infection times (generation time and infectivity profile) and although we have used well-supported estimates, there’s inherent limitations to deriving these distributions based on transmission pairs. These transmission pairs are representative of symptomatic cases, but the infectiousness profiles for fully asymptomatic cases are unknown (Ferretti et al., 2020a). We have assumed that asymptomatic cases have the same infectiousness profiles as symptomatic cases, but if asymptomatic cases are infectious for a shorter duration, or have a lower probability of transmission during a contact (Buitrago-Garcia et al., 2020), then we would overestimate the transmission prevented by quarantining these cases. We do account for uncertainty in the infection time distributions, and this uncertainty is carried through into our analysis and is captured by the confidence intervals shown in the figures and reported in the text.
In terms of modelling the TTIQ process, we have assumed that identified index cases are isolated and have their contracts traced. If the index case fails to adhere to the isolation protocol, then we will overestimate the amount of transmission prevented by isolation. However, uncertainty in whether contacts adhere to quarantine protocols, or whether contact tracers actually identify contacts, is contained in the parameter g. Lower adherence to quarantine or missed cases due to overwhelmed contact tracers is captured by lowering g.
Here we have shown through systematic analysis that TTIQ processes can be optimised to bring the effective reproductive number below one. Crucially, contact tracing & quarantine adds security to testing & isolating strategies, where high coverage and short delays are necessary to control an epidemic. By improving the testing & isolation coverage and reducing the delay to index case isolation, we can greatly increase the efficacy of the overall TTIQ strategy.
Data Availability
Code is publicly available at http://github.com/ashcroftp.