Screening plans for SARS-CoV-2 based on sampling and rotation: an example in the school setting =============================================================================================== * Michela Baccini * Giulia Cereda ## Abstract Screening plans for prevention and containment of SARS-CoV-2 infection should take into account the epidemic context, the fact that undetected infected individuals may transmit the disease, and that the infection spreads through outbreaks, creating clusters in the population. In this paper, we compare the performance of six screening plans based on poorly sensitive individual tests, in detecting infection outbreaks at the level of single classes in a school context. The performance evaluation is done by simulating different epidemic dynamics within the class during the five weeks following the day of the first infection. The plans have different costs in terms of number of individual tests required for the screening and are based on recurrent evaluations on all students or subgroups of students in rotation. Especially in scenarios where the rate of contagion is high, at an equal cost, testing half of the class in rotation every week appears to be better in terms of sensitivity than testing all students every two weeks. Similarly, testing one-fourth of the students every week is comparable with testing all students every two weeks, despite the first one is a much cheaper strategy. In the presence of natural clusters in the population, testing subgroups of individuals belonging to the same cluster in rotation may have a better performance than testing all the individuals less frequently. The proposed simulations approach can be extended to evaluate more complex screening plans than those presented in the paper. ## INTRODUCTION Since the beginning of the COVID-19 emergency in the early 2020, the importance of implementing extensive screening procedures to prevent or slow down the spread of the SARS-CoV-2 infection has been emphasized [1, 2] and, in light of the threat of new variants of the virus that could be more widespread and the critical issues related to the rapid implementation of vaccination plans [3, 4], it still seems early to consider extensive surveillance strategies on the population (or specific subgroups of it) no longer necessary. Pharmaceutical industries have produced tests of various nature and cost, which have been proposed and used in screening plans aimed at early detection of asymptomatic or paucisymptomatic individuals in specific populations. The ability of these tests to correctly classify the single patient as infected or not has been widely discussed and debated, sometimes overshadowing the necessity that screening plans account for strengths and limitations of the used tests and are tailored to the specific context in which they are applied [5, 6]. When dealing with a screening plan for an infectious disease, two points should not be overlooked: * undetected infected subjects can transmit the disease; * the infection usually spreads in small outbreaks, creating clusters in the population (families, classes, work colleagues). For these two reasons, screening plans similar to those implemented in the case of non-communicable diseases could lead to suboptimal results in terms of cost-benefit ratio. Furthermore, it is crucial that the screening procedures are assessed accounting for the actual epidemic context, including the strength of contagion [7]. In this paper, we compare through simulations alternative strategies designed for screening in schools, based on repeated tests to be performed at regular time intervals on all students, or on tests to be performed in rotation on subgroups. The comparison, which refers to a single class of *N* = 24 students - the class is here defined as a group of students who attend the same course each day at university or school - takes into account the epidemic context and is performed under alternative scenarios of epidemic spread. Although the simulation analysis refers to school settings, the methods adopted and the results obtained are valid for any context in which there are natural clusters of individuals, within which contagion could spread starting from a single initial infection. ## METHODS Our objective is to early identify infections in school settings, in order to quarantine the classes where there is at least one infected student. Let us suppose that individual tests are performed on all students in the class or on a subset of them. The class is considered positive if at least one of the individual tests is positive; it is considered negative if no individual tests are positive. Therefore, at the class level, the sensitivity of the test is defined as the probability that the class is positive given that there is at least one infected student present. Specificity is the probability that the class is negative given that there are no infected students in it. In our analysis, we assume that the individual tests used in the screening procedure have maximum specificity and sensitivity equal to *p*. The assumption of maximum specificity of the individual test implies maximum specificity at the class level and rules out those situations where false positives may lead to quarantining classes when not necessary. This assumption allows us to simplify subsequent calculations without compromising evaluations regarding the ability of the proposed plans to detect outbreaks. ### Screening plans We consider six screening plans (Figure 1), which differ from each other for the time interval between consecutive evaluations on the class and number of students involved in each of them: ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/12/2021.02.10.21251502/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/02/12/2021.02.10.21251502/F1) Figure 1. Example of an epidemic curve and its intersection with the screening plans (A1, A2, B1, B2, C, D). Dashed and solid red lines: evaluations carried out on the class according to the screening plans; Solid blue line: beginning of the epidemic; Grey area: time window on which the assessment is made (4 weeks from the beginning of the epidemic); Dashed blue lines: 1, 2, 3 and 4 weeks from the beginning of the epidemic; Solid red lines: evaluations carried out on the class within four weeks from the beginning of the epidemic. * Plan A1. Individual tests on all students of the class every week; * Plan A2. Individual tests on all students of the class every 2 weeks; * Plan B1. Individual tests every week on 1/2 of the students of the class, in rotation; * Plan B2. Individual tests every 2 weeks on half of the students of the class, in rotation; * Plan C. Individual tests every 10 days on 1/3 of the students of the class, in rotation; * Plan D. Individual tests every week on 1/4 of the students of the class, in rotation. Plan A1 guarantees the best performance in terms of surveillance but requires more resources compared to the others. Plans B1, B2, C, and D have an additional element of risk compared to plan A2 because they test each time sub-samples of the students of the class. Plans B1, C and D allow more frequent monitoring of the class compared to plan A2 and B2. We assume *p*=0.7, which is close to the average sensitivity reported by OECD for antigenic tests [6]. Then, in order to perform a sensitivity analysis, we assume *p*=0.9 as well. ### Epidemic scenarios We assume that an infected student in the class may generate new infections among his/her classmates. Let us suppose that one of the *N* students is infected on day 1. The epidemic dynamic within the class can be simulated by using a compartmental model [8], where the contagion strength depends on the average time of infectivity *T* and on the basic reproduction number *R*, which is the number of secondary infections generated from the first infected student in the class. In particular, we assume that, on average, each infected student may spread the contagion in the class for *T* days, still remaining detectable as infected for 4 weeks [9]. We consider different epidemic contexts, characterized by different combinations of *R* and *T* (*R*=1.1, 1.5, 3, 5 and *T*=7, 14, 21, for a total of 12 scenarios) [10], and we apply on each of them the six screening plans. The ratio *β*=*R**/T* is proportional to the contagion rate; combinations of *R* and *T* which result in the same *β* generate the same epidemic dynamic, net of stochastic variability. This is the case for *R*=2, *T*=14 and *R*=3, *T*=21. ### Simulations Separate simulations, for a minimum of 7000, are performed for each combination of *R* and *T* and each screening plan. We assume that the number of new infections and the number of infected that become not infectious at time *t* in the class, *I**new*(*t*) and *R**new*(*t*), follow Binomial distributions: ![Formula][1] where *S*(*t* − 1) and *I*(*t* − 1) are the number of susceptible students and the number of infectious students at time *t* − 1, respectively. We further assume that the groups, when required by the screening plan, are randomly generated, that the probability of becoming infected for a susceptible subject does not depend on the group to which he/she belongs, and that the new infections at each time are randomly distributed among the groups. Let us suppose that *g*=2 groups of students *G*1 and *G*2 have been created and the susceptible individuals in the two groups at time *t* are *S*1(*t*), *S*2(*t*). The number of new infections in *G*1 is sampled from a Hypergeometric distribution: ![Formula][2] where *M* is the population size, *K* is the number of successes in the population, and *n* is the number of draws. Then the number of new infections in *G*2 is obtained as difference: ![Graphic][3]. In general, if *g*>2 the number of new infections in the groups is obtained by sampling from a Multivariate Hypergeometric distribution. The results of the individual tests on different individuals are assumed to be independent and follow a Bernoulli distribution with parameter *p*. According to the assumption of maximum specificity, the result of the test on the not infected students is always negative. In the simulations, we allow the epidemic to originate at any time between two consecutive assessments of the class and we focus on a time window of 4 weeks from the first infection (Figure 1) to evaluate the performance of the screening plan in terms of: * probability of detecting the outbreak at 7, 14, 21 and 28 days since its onset; * total number of infection-days which are left undetected by the screening plan in the 4 weeks time window. An infection-day is here defined as a day spent by a subject in the infectious status, thus a day in which he/she can spread the contagion. ## RESULTS Figure 2 reports, for each scenario, the probabilities of a positive result on the class within 7, 14, 21 and 28 days from the beginning of the epidemic. The curve describes the overall performance of the screening plans in detecting the presence of infections in the class. As expected, plan A1 guarantees the best performance in terms of epidemic detection. The cumulative probabilities for plans A2 and B1 are very similar in scenarios where the infection spreads slowly within the class (upper left quadrant), while plan B1 seems to have better relative performance compared to A2 in high-epidemic spread contexts (bottom right quadrant). Screening plans B2, C and D have the worst performance if the infection spread is low, but they reach good results in high-risk scenarios. Plan D seems to detect infections within 2 weeks from the beginning of the epidemic with a probability between 70% and 80% in scenarios where *β* ≥ 0.36. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/12/2021.02.10.21251502/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/02/12/2021.02.10.21251502/F2) Figure 2. Probability of detecting the infection within t days from the beginning of the epidemic in the class. Comparison of the performance of the six screening plans (A1, A2, B1, B2, C, D), under different epidemic scenarios, assuming that individual tests have sensitivity 0.7 and maximum specificity. The probabilities that the screening plans do not detect the infection within 4 weeks from the beginning of the epidemic are reported in Table 1. The risk of not detecting the outbreak at all decreases as the rate of infection within the class increases. For plan A1, the probability of a false negative is always negligible (lower than 0.4%). For plans A2 and B1 probabilities are very similar (from 1.0% when *β* = 0.71 to 5.5% when *β* = 0.05), as well as for plans B2, C and D (from 3.4% when *β* = 0.71 to 24.4% when *β* = 0.05). View this table: [Table 1.](http://medrxiv.org/content/early/2021/02/12/2021.02.10.21251502/T1) Table 1. Probability of not detecting the outbreak in the class within 35 days from the beginning of the epidemic, by screening plan (A1, A2, B1, B2, C, D), under different epidemic scenarios (*R*0 and *T*), assuming that individual tests have sensitivity 0.7 and maximum specificity. Figure 3 shows, for each scenario and screening plan, the distribution of the number of lost infection-days over a time window of 4 weeks, arisen from the simulations. The colors indicate the cost of the plans in terms of individual tests. If *x* is the cost of plans B2, C and D (green), plan A2 and B1 cost 2*x* (yellow) and plan A1 costs 4*x* (red). The number of lost infection-days increases as *β* increases. After A1, B1 is the plan which assures the lowest number of lost infection-days. If we exclude scenarios of low epidemic spread, when plan B2, C and D are comparable in terms of lost infection-days, of the three cheaper plans, plan D seems the one that guarantees the lowest number of lost infection-days. Interestingly, plan D is equivalent or better than A2 in scenarios where *β* ≥ 0.36. Plan B1 leaves less undetected infection-days than plan A2. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/12/2021.02.10.21251502/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/02/12/2021.02.10.21251502/F3) Figure 3. Boxplots of the lost infection-days by screening plan (A1, A2, B1, B2, C, D), under different epidemic scenarios, assuming that individual tests have sensitivity 0.7 and maximum specificity. In green the less expensive plans, in yellow the medium and in red the most expensive. In Figure 4, focusing on plan A1 and on plans B1 and D, which are the best ones within their cost range, we compare the cumulative probabilities of a positive result when assuming *p*=0.7 and 0.9, under the scenarios characterized by the larger and the lower rates of contagion. Increasing the sensitivity of the individual test to 0.9, the performance of the screening plans increases, but their relative accuracy remains similar. This result arises also from the comparison of Tables S1 and S2 (Supplemental Material), which report averages and 90th percentiles of the number of lost infection-days for the six screening plans and the 12 epidemic scenarios, for *p*=0.7 and 0.9, respectively. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/12/2021.02.10.21251502/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2021/02/12/2021.02.10.21251502/F4) Figure 4. Probability of detecting the infection within *t* days from the beginning of the epidemic in the class. Comparison of the performance of the screening plans A1, B1, D, under the scenarios characterized by the lower (left panel) and the largest (right panel) rates of contagion, assuming that individual tests have sensitivity 0.7 (solid line) and 0.9 (dashed line) and maximum specificity. ## DISCUSSION In this paper, we perform simulations to compare the performance of six screening plans based on individual tests with low sensitivity and maximum specificity in detecting infection outbreaks at the class level in schools or, more in general, at the cluster level in a population. We account for uncertainties around the epidemic dynamic in the class through a stochastic compartmental model and around the screening plan implementation through random generation of groups, when required by the plan, and random allocation of the new infections across them. Our work is not far from others which used mathematical models or simulations to investigated relevant issues during the COVID-19 emergency, such as the definition of optimal quarantine strategies [11] or optimal pool size in pooled testing [12, 13]. The compared plans have different costs and are based on recurrent evaluations on all students or on subgroups of students on rotation. Among all possible plans that perform assessments on the students at time intervals greater than one week, the best option obviously consists in testing all students at weekly intervals. This option assures the maximum sensitivity and the lowest number of infection-days left undetected. However, it is very expensive, hence the need of exploring cheaper strategies. Without claiming to be exhaustive, we consider and compare through simulations five alternatives to the best option. At an equal cost, testing half of the class on rotation every week proves to be better than testing all students every two weeks, because it allows to earlier detect the presence of infections, especially in scenarios of a high rate of contagion in the class. If resource constraints are even tighter, less expensive plans can be considered: testing half of the class once every two weeks, testing one-third of the class every 10 days, testing one-fourth of the class every week. These plans, which have same cost, perform similarly in case of a low rate of contagion. However, in case of a high rate, testing one-fourth of the class every week seems to be the best option. This again suggest that reducing times between successive evaluations may largely balance the risk related to a lower class coverage at each assessment. Moreover, considering that strength and speed of virus transmission within each single class are not known in advance, it seems reasonable to focus on high-transmission scenarios, potentially more dangerous in terms of infection spread within and outside the class [14]. Interestingly, the plan which tests one-fourth of the students every week results to be comparable to the plan that tests all students every two weeks in terms of sensitivity and turns out to be better in terms of lost infection-days. This suggests that reducing costs by simply increasing the interval between successive assessments should not be considered a priori a good option, because much less expensive strategies based on testing subgroups on rotation could end up performing similarly. The absolute performance of the screening plans increases with the sensitivity of the individual tests used, but their relative performance remains unchanged. It should be noticed that the sensitivity value of 0.9 roughly corresponds to the one that would be obtained by repeating a test with a sensibility of 0.7 twice on the same individual, in a procedure of double testing [6, 15]. Our simulations can be extended also to screening plans which combine rotation and sample pooling strategies, that could be less costly than those based on individual tests [16]. Each group of six students in plans C or D could be tested in pool with only one molecular test on mixed material from individual swabs. Similarly, the screening on the entire class could be done by performing only one pooled test on the whole class or by defining *k*≥1 subgroups of students, then performing *k* pooled tests. However, the sensitivity of pooled testing for various pool sizes and its relationship with viral loads distribution in the population should be considered essential inputs to assess and compare the performance of plans which involve sample pooling [17]. A final remark concerns our assumption of maximum specificity. This assumption affects the calculation of the sensitivity at the class level in a conservative sense: if a student is falsely declared positive in a class where there is at least one infected, the probability of a positive result at the class level increases. On the other hand, a suboptimal, even if high, specificity could induce a not negligible number of inappropriate quarantines, with the resulting social costs. We do not address this issue in the paper, but it is worth noting that the risk of inappropriate quarantine is lower if tests are performed on smaller groups of students. For example, if the individual test has a specificity of 0.99, the specificity on the class at each assessment is 0.996=0.94 under plan C and 0.9924=0.79 under plan A1. In conclusion, our simulations emphasize the importance of taking into account the epidemic context when comparing alternative screening strategies aimed at contagion prevention. Combining calculation of test accuracy with contagion dynamic models is an easy way to do this, even when the plans to be evaluated are more complex than those considered in this paper. We show that, in the presence of clusters in the population, less costly strategies that exploit the correlation between subjects may have good performance, in particular in detecting outbreaks at high-rate of contagion, and that the time interval between successive evaluations on the same cluster is a very relevant input. The closer the assessments, the lower the number of infection-days left undetected, with a reduction of the risk that infectious subjects spread the contagion within the cluster and outside it. In real applications, economic and social costs as well as positive and negative predictive values of the proposed plans, should be evaluated making assumptions about the prevalence of infections in the population. ## Supporting information Supplemental Material [[supplements/251502_file02.pdf]](pending:yes) ## Data Availability The paper is based on simulations. ## Competing Interests The authors have no competing interests. * Received February 10, 2021. * Revision received February 10, 2021. * Accepted February 12, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].CDC. CDC Guidance for Expanded Screening Testing to Reduce Silent Spread of SARS-CoV-2; 21 January 2021. Available from: [https://www.cdc.gov/coronavirus/2019-ncov/php/open-america/expanded-screening-testing.html](https://www.cdc.gov/coronavirus/2019-ncov/php/open-america/expanded-screening-testing.html). 2. [2].OECD. OECD Policy Responses to Coronavirus (COVID-19). Testing for COVID-19: A way to lift confinement restrictions; May 2020. Available from: [http://www.oecd.org/coronavirus/policy-responses/testing-for-covid-19-a-way-to-lift-confinement-restrictions-89756248](http://www.oecd.org/coronavirus/policy-responses/testing-for-covid-19-a-way-to-lift-confinement-restrictions-89756248). 3. [3].COVID-19 vaccines: no time for complacency. Editorial. The Lancet. 2021;396:1607. 4. [4].Burki T. Understanding variants of SARS-CoV-2. The Lancet. 2021;397:462. 5. [5].Cheng MP, Papenburg J, Desjardin M, Kanjilal S, Quach C, Libman M, et al. Diagnostic Testing for Severe Acute Respiratory Syndrome–Related Coronavirus 2. Annals of Internal Medicine. 2020;172:726–734. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F12%2F2021.02.10.21251502.atom) 6. [6].OECD. OECD Policy Responses to Coronavirus (COVID-19). Testing for COVID-19: How to best use the various tests? Available from: [http://www.oecd.org/coronavirus/policy-responses/testing-for-covid-19-how-to-best-use-the-various-tests-c76df201](http://www.oecd.org/coronavirus/policy-responses/testing-for-covid-19-how-to-best-use-the-various-tests-c76df201). 7. [7].Baccini M, Mattei A, Rocco E, Vannucci G, Mealli F. Evaluating a SARS-CoV-2 screening strategy based on serological tests. Epidemiologia & Prevenzione. 2020;44:184– 192. 8. [8].Brauer F, van den Driessche P, Wu J. Compartmental Models in Epidemiology. In: Mathematical Epidemiology. Springer Berlin Heidelberg. 2008:19–79. 9. [9].Walsh K, Jordan K, Clyne B, Rohde D, Drummond L, Byrne P, et al. SARS-CoV-2 detection, viral load and infectivity over the course of an infection. J Infect. 2020;81:357– 371. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2020.06.067&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F12%2F2021.02.10.21251502.atom) 10. [10].Inglesby T. Public Health Measures and the Reproduction Number of SARS-CoV-2. JAMA. 2020;323:2186–2187. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.7878&link_type=DOI) 11. [11].Wells CR, Townsend JP, Pandey A. et al. Optimal COVID-19 quarantine and testing strategies. Nat Commun. 2021;12:356. [https://doi.org/10.1038/s41467-020-20742-8](https://doi.org/10.1038/s41467-020-20742-8) 12. [12].Rewley, Jeffrey. Specimen Pooling to Conserve Additional Testing Resources When Persons’ Infection Status Is Correlated. Epidemiology. 2020;31:832–835. 13. [13].Regen F, Eren N, Heuser I, Hellmann-Regen J. A simple approach to optimum pool size for pooled SARS-CoV-2 testing. International Journal of Infectious Diseases. 2020;100:324–326. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.08.063&link_type=DOI) 14. [14].Majra D, Benson J, Pitts J, Stebbing J. SARS-CoV-2 (COVID-19) super-spreader events. J Infect. 2021;82:36–40. 15. [15].Griffin DWJ, Bramwell E, Cheng AC, Spelman DW. Diagnostic testing for COVID-19: demonstrated utility of duplicate testing for inpatients in a low incidence setting. Internal Medicine Journal. 50:1594–1595. 16. [16].Hogan C, Sahoo M, Pinsky B. Sample Pooling as a Strategy to Detect Community Transmission of SARS-CoV-2. JAMA. 2020;323(19):1967–1969. 17. [17].Nguyen N, Aprahamian H, Bish E, Bish D. A methodology for deriving the sensitivity of pooled testing, based on viral load progression and pooling dilution. J Transl Med. 2019;17(1):252. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12967-019-1992-2&link_type=DOI) [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/inline-graphic-1.gif