Regular testing of asymptomatic healthcare workers identifies cost-efficient SARS-CoV-2 preventive measures =========================================================================================================== * Daniel Sanchez-Taltavull * Violeta Castelo-Szekely * Shaira Murugan * Tim Rollenske * Stephanie C. Ganal-Vonarburg * Isabel Büchi * Adrian Keogh * Hai Li * Lilian Salm * Daniel Spari * Bahtiyar Yilmaz * Jakob Zimmermann * UVCM-COVID researchers * Michael Gerfin * Edgar Roldan * Guido Beldi ## Abstract Protecting healthcare professionals is crucial in maintaining a functioning health-care system. The risk of infection and optimal preventive strategies for health-care workers during the COVID-19 pandemic remain poorly understood. Here we report the results of a weekly testing regime that has been performed since the beginning of the COVID-19 pandemic to identify pre- and asymptomatic healthcare workers. Based on these observations we have developed a mathematical model of SARS-CoV-2 transmission that integrates the sources of infection from inside and outside the hospital. The data were used to study how regular testing and a desynchronisation protocol are effective in preventing transmission of COVID-19 infection at work, and compared both strategies in terms of workforce availability and cost-effectiveness. We showed that case incidence among healthcare workers is higher than would be explained solely by community infection. Furthermore, while testing and desynchronisation protocols are both effective in preventing nosocomial transmission, regular testing maintains work productivity with implementation costs. ## 1 Introduction During the outbreak and spread of an epidemic such as that of SARS-CoV-2, protection of individuals relies on non-pharmaceutical interventions (NPIs) at least until a large fraction of the population is vaccinated. Common NPIs include wearing protective equipment, social distancing and working from home, as well as isolation of symptomatic individuals. However, applying some of these measures among healthcare workers is challenging, as they need to work in close contact with infected patients and with their coworkers, and as the possibilities for home-office options are limited. Therefore, the healthcare workforce could be exposed to an increased risk of infection transmission than the general population and devising protective strategies for them remains a crucial task. In addition to the above, other NPIs have been considered through modelling and case studies, notably, testing of healthcare personnel and implementation of desynchronisation to reduce worker contact. Studies of testing healthcare workers have shown differing results with regard to whether they have a reduced or increased case incidence with respect to the general population. The incidence among health workers seems to be increased [12] although it mirrored the curve in the general population [20]. Other studies showed no major occupational risk for healthcare professionals [7, 22]. Together, these results indicate that the exposure risk among healthcare workers likely depends on the specificities of the study as well as on the local magnitude of the epidemic, the availability and implementation of basic protective measures and the effective adherence to them in the studied scenario. As a caveat, most of these studies carried out testing only on symptomatic individuals and did not assess whether widespread regular testing of asymptomatic and pre-symptomatic coworkers reveals a higher incidence and can prevent transmission at the workplace in an epidemic in which a significant proportion of affected individuals can present without symptoms [5, 6]. The implementation of desynchronising the presence of workers at the work-place has been modelled in the context of the SARS-CoV-2 epidemic. A model of 4-day work, 10-day lockdown cycle showed that these cycles can allow economic activity to partially continue by exploiting the latent period of the virus and allowing those infected at work to reach maximum infectiousness while at home [14]. Similarly, organisation of laboratory or clinical staff into smaller groups working shorter shifts can lead to reduced interaction and transmission among coworkers [21, 19]. Yet such models have not been tested in with experimental data during the current pandemic. In this study we report the results of a testing pipeline of one clinical department that has been introduced from the beginning of the epidemic in Switzerland. Starting from March 23rd 2020, we implemented a RT-PCR-based testing pipeline that was applied weekly on a voluntary basis to a cohort of asymptomatic healthcare workers. Based on the results of this testing pipeline we have developed a mathematical model that includes the transmission dynamics of SARS-CoV-2 in the Swiss Canton of Bern, and that integrates the infection forces among a group of healthcare workers - including the general population, hospitalised patients and coworker - as independent sources of infection. Individuals testing positive were referred for confirmation to a second testing pipeline carried out by the official authorities, and sent into self-isolation. We used this data not only to detect and isolate pre- and asymptomatic cases, but also to inform the model on the contribution of worker-worker and worker-patient contact to case incidence among the studied healthcare group. As a comparison, we have modelled a work desynchronisation strategy in which healthcare workers would shift weekly between in-hospital and at-home work, as a way to contain infection transmission by reducing contact between colleagues. We fitted the model to the recorded number of hospitalised individuals outside the hospital to infer the time-varying infection rate of the general population, and we estimated the additional infection forces affecting the health workforce by fitting the individuals testing positive and sent to quarantine. The model allowed us to study how both regular testing of all individuals and desynchronisation strategies are efficient in preventing workers from infecting each other under the different scenarios, and we present how the model can help decision making on whether to implement such strategies. Furthermore, we studied the impact that the implementation of these interventions would have on productivity and the availability of the workforce, and compare them from an economic point view. ## 2 Methods ### 2.1 Epidemiological data We collected publicly available data on patient hospitalizations, deaths and confirmed cases in the Canton of Bern from March 1st 2020 to March 1st 2021, from the BAG (Bundesamt für Gesundheit) [1]. ### 2.2 Internal testing pipeline The Study was approved by the ethical committee of the Canton of Bern (2020-00563). A testing protocol was established and implemented weekly on a voluntary basis, starting March 23, 2020. Consenting participants could get tested either on Monday or Thursday every week, and positive cases were sent for a confirmatory test carried out by official institutions and isolated. #### Sample Collection Healthcare workers collecting samples from the participants were instructed to use personal protective equipment including FFP-2 face masks, eye protection and face shield, two layers of gloves and two layers of hairnets. Nasopharyngeal specimens were obtained by slowly and gently inserting a MiniTip swab through the nostril palate and closely to the septum until contact with the nasopharynx was reached. The swab was removed after gently rubbing and rolling at the nasopharynx. Rarely, oral swabs were taken instead by taking samples from the posterior oropharynx. Specimens taken from each participant were immediately placed into 2-ml tubes containing TRIzol (Invitrogen), an acid-guanidinium-phenol based reagent that inactivates the virus and stabilises the RNA. All collected samples were sprayed with 70% EtOH before leaving the collection site. Specimens were further processed on the same day or kept at −20°C until processed the following day. #### RNA extraction and RT-qPCR Laboratory workers handling the collected samples were instructed to use personal protective equipment including gowns, facemask and two layers of gloves. The virus was inactivated by vortexing the TRIzol tube containing the swabtips and 3.3 ng carrier RNA (Qiagen) was added. Then 200 *µ*l of chloroform was added to each sample and the tubes were inverted multiple times, followed by centrifugation for 10 min at 13,200 rpm. RNA purification and concentration was carried out using RNA easy-mini elute kit (Qiagen) according to the manufacturer’s protocol. Briefly, 270 *µ*l of the clear aqueous phase of each sample was re-suspended with 945 *µ*l of RLT buffer and 670 *µ*l of 100% EtOH. The sample was added to the column in three steps with intermittent centrifugation at 8000 rpm for 15 s allowing all sample to pass through the column. 500 *µ*l of RPE buffer was added to each column as a first washing step and centrifuged at 8000 rpm for 15 s. Subsequently, 500 *µ*l of 80% EtOH was added to each column and centrifuged at 13,200 rpm for 1 min. After drying the RNA containing membrane by centrifugation at 13,200 rpm for 3 min, 20 *µ*l RNAse-free water was added to elute the RNA. The sample was incubated for 1 min at room temperature before it was centrifuged for 1 min at 13200 rpm. One-step combined reverse transcription and qPCR using SuperScript III One-Step RT-PCR System with Platinum Taq DNA Polymerase (Invitrogen) was performed using the manufacturers recommended protocol. Briefly, 2 *µ*l of each RNA sample was amplified using primers for SARS-Cov2 E-gene, which codes for protein forming pores on the host membrane, and human RNase-P, in order to verify successful sampling of the healthcare workers. Technical triplicates were performed. The PCR was run under following conditions: 50°C for 10 min, 95°C for 3 min, followed by 45 cycles of 95°C for 15 s and 58°C for 30 s on the BioRad CFX 384 real time PCR system. ### 2.3 Mathematical modelling #### 2.3.1 Transmission dynamics in the city model We developed a population-based model that integrates the transmission dynamics in the Swiss Canton of Bern as a deterministic model, and the infection dynamics among a group of healthcare workers, studied as a stochastic system (Figure 2A). The Canton dynamics are given by the following system of differential ordinary equations (ODE): ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] Susceptible individuals in the Canton, *S**C*, get exposed at a rate *α*; exposed individuals, *E**C*, then become infectious, *I**C*, after a latency period *s*−1. A proportion *ϵ*1 of infected individuals are hospitalised, *H**C*, at a rate *h*. Infectious individuals recover, *R**C* at a rate *r*, whereas hospitalised patients recover at a rate *w*. We only considered deaths from hospitalised people: a proportion *ϵ*4 of *H**C* die at a rate *d. N* refers to the total population in the Canton of Bern. As in Althaus et al. [4], the transmission rate *α* captures the reduction in the infections following the implementation of NPIs by the government, described as: ![Formula][7] where ![Graphic][8] is the infection at the start of the epidemic, *κ* is the relative transmission reached after implementation of NPIs, *ν* is the slope of the sigmoid function, and *τ* is the midpoint of transmission reduction. #### 2.3.2 Transmission dynamics among the healthcare workers model The transmission dynamics among the HCW are formulated in terms of the Master Equation [9] ![Formula][9] where *X* = (*S, E, I, A, R, Q*) is the state vector containing the number of healthy workers susceptible to infection, *S*, the number of exposed workers who are infected but still not contagious that are working, *E*, the number of symptomatic infected workers, *I*, the number of asymptomatic infected workers, *A*, that are working, the number of recovered workers *R*, and the number of exposed and asymptomatic workers that have been identified as infected with a PCR test and put in quarantine, *Q. P* (*X, t*) is the probability of the population of having value *X* at time *t*. The transition rate per unit time of the process *j* is given by *W**j*(*X, t*), and *r**j* is the state change in *X* upon the occurrence of the process *j*. All reactions and rates are described in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2021/06/02/2021.06.01.21258124/T1) Table 1: Reaction rates of the stochastic dynamics described by Equation 8. The mean-field dynamics of the stochastic system are given by the following system of ODEs: ![Formula][10] ![Formula][11] ![Formula][12] ![Formula][13] ![Formula][14] ![Formula][15] Susceptible workers, *S*, become exposed at a rate *α* due to contact with the general population (*I**C*), at a rate *β* from being in contact with an asymptomatic coworker (*A*, see later), and at a rate *γ* from contact with infected hospitalised individuals (*H**C*). Exposed individuals (*E*) then develop the disease after a latency period *s*−1, and become either infectious and symptomatic, *I*, or infectious and asymptomatic, *A*, with the probability of being *I* or *A* given by *ϵ*2 and *1-ϵ*2 respectively. Exposed and asymptomatic can be identified by testing, where *ϵ*3 is the probability of obtaining a true positive in the test, and the details of the testing are described below. For simplicity, we write *p*(*t*) in the equations to show it is time dependent. Finally, *r* is the recovery rate. Both asymptomatic and symptomatic individuals recover at a rate *r*. We did not considered deaths among healthcare workers. *N**H* refers to the number of healthcare workers participating in the study. Due to their small number, their transmission dynamics are studied as a stochastic process using Gillespie simulations [11]. ##### Modelling a desynchronisation strategy We modelled a desynchronisation strategy in which the workforce was split into two non-mixing groups to establish a weekly rotation shift: Group 1 works in-hospital during odd weeks and at-home during even weeks, and Group 2 works at-home during odd weeks and in-hospital during even weeks. To mathematically formulate the desynchronisation, each variable is dichotomised in two non-interacting groups (*S*1, *E*1, *I*1, *A*1, *R*1, *Q*1) and (*S*2, *E*2, *I*2, *A*2, *R*2, *Q*2), each of them described by Eq. 8. When workers are doing home office, they cannot be infected by other co-workers nor by patients. That is, the parameters *β* and *γ* become time dependent, being *β*1(*t*) = 1.5*β, γ*1(*t*) = *γ, β*2(*t*) = 0, *γ*2(*t*) = 0 during the odd weeks, and exchanging values during the even weeks, *β*1(*t*) = 0, *γ*1(*t*) = 0, *β*2(*t*) = 1.5*β, γ*2(*t*) = *γ*. We assume an increase of 50% in *β* because some tasks imply cooperation between workers. ##### Modelling testing We modelled the testing as an instantaneous reaction in which *E* → *Q* or *A* → *Q*. Every *P* days, for each infected worker we generate a random number, *x* in (0, 1), if *x* < *ϵ*3*f**p*, where *f**p* is the fraction of workers tested and *ϵ*3 is the probability of detecting a true positive, the worker is sent to quarantine. #### 2.3.3 Model Fitting ##### City In Switzerland the first lockdown was started on March 16th 2020 that included reduction of elective healthcare-associated interventions such as surgery until April 27th 2020. With the onset of the second wave national measures were again installed by October 19th 2020. After a second peak in 30th October 2020 the second wave faded by end of January 2021. *τ, π, ν* and *κ* were chosen to obtain the best fit of the model to the data, by minimizing ![Formula][16] with *t**i* ∈ {0, 7, 14, …}, *O*(*t**i*) the sum of the number of cases of the last 7 days at time *t**i* starting March 1st 2020 (*t* = 0) and ![Formula][17] The fitting was done in 4 parts: 1. Wave 1: With the initial condition (*S**C*, *E**C*, *I**C*, *H**C*, *R**C*, *D**C*) = (*N*−1, 0, 1, 0, 0, 0) starting Day −8 to Day 91, with *α*(*t*) described by 7 resulting into *τ* = 17, *π* = 1.3, *κ* = 0.016, *ν* = 0.98. The optimization was performed in the space region *π* ∈ [0, 2] with step size Δ*π* = 0.05, *κ* ∈ [0, 0.15] with step size Δ*κ* = 0.002, *ν* ∈ [0, 1.5] with step size Δ*ν* = 0.02, *τ* ∈ [15, 24] with step size Δ*τ* = 1. 2. Summer: Day 91 to Day 189 with *α*(*t*) described by ![Formula][18] resulting into *π* = 0.2859. The optimization was performed in the space region *π* ∈ [0, 1] with step size Δ*π* = 0.0001 3. Wave 2: Day 189 to Day 273, with *α*(*t*) described by 7 resulting into *τ* = 239, *π* = 0.32, *κ* = 0.44, *ν* = 0.54. The optimization was performed in the space region *π* ∈ [0, 0.7] with step size Δ*π* = 0.02, *κ* ∈ [0, 1] with step size Δ*κ* = 0.02, *ν* ∈ [0, 1.5] with step size Δ*ν* = 0.02, *τ* ∈ [230, 246] with step size Δ*τ* = 1. 4. Winter: Day 273 to Day 365, with *α*(*t*) described by 7 resulting into *τ* = 299, *π* = 0.2, *κ* = 0.58, *ν* = 0.32. The optimization was performed in the space region *π* ∈ [0, 0.5] with step size Δ*π* = 0.02, *κ* ∈ [0, 1] with step size Δ*κ* = 0.02, *ν* ∈ [0, 1] with step size Δ*ν* = 0.02, *τ* ∈ [295, 305] with step size Δ*τ* = 1. ##### Hospital In order to protect workers, the department started a desynchronisation strategy on March 16th 2020 and ended on May 4th 2020. Moreover, the department implemented a weekly testing regime starting March 23rd 2020. During the desynchronisation period, only the in-hospital group were tested, on Monday and Thursday. For simplicity we assume all workers, in-hospital or doing home office, were tested once a week and the same day. In order to mimic the hospital dynamics, we simulated our system as follows. 1. Set the initial condition for (*S**C*, *E**C*, *I**C*, *H**C*, *R**C*, *D**C*) = (*N*, 0, 1, 0, 0, 0) at time *t* =−8 and integrate the ODE system 1-6 until day 15. Note that March 1st 2020 correspond to *t* = 0. 2. set (*S*1, *E*1, *A*1, *I*1, *Q*1, *R*1) = (157, 0, 0, 0, 0, 0), (*S*2, *E*2, *A*2, *I*2, *Q*2, *R*2) = (157, 0, 0, 0, 0, 0), and start the stochastic simulation. 1. Select reaction *r**i* and Δ*t* according to the Gillespie algorithm. 2. update (*S*1, *E*1, *A*1, *I*1, *Q*1, *R*1), (*S*2, *E*2, *A*2, *I*2, *Q*2, *R*2) and *t* accordingly. 3. integrate system of ODE 1-6 for a time step of Δ*t*. 4. If *t* ∈ (15, 64) 1. if floor(t/7)%2 = 0, then *β*2 = *γ*2 = 0, *β*1 = 1.5*β, γ*1 = *γ*. 2. if floor(t/7)%2 = 1, then *β*1 = *γ*1 = 0, *β*2 = 1.5*β, γ*2 = *γ*. 3. if floor(t) % 7 = 1 AND Bol=1 AND *t* > 22, we test * Bol=0 * Set Tests=E1 * for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3 ×Fraction of workers tested, then Q1=Q1+1 E1=E1-1 * Set Tests=A1 * for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3× Fraction of workers tested, then Q1=Q1+1 A1=A1-1 * Set Tests=E2 * for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3 ×Fraction of workers tested, then Q2=Q2+1 E2=E2-1 * Set Tests=A2 * for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3× Fraction of workers tested, then Q2=Q2+1 A2=A2-1 4. if floor(t) % 7 /= 1: Bol=1 5. If *t* > 64 * S1=S1+S2, S2=0 * E1=E1+E2, E2=0 * A1=A1+A2, A2=0 * I1=I1+I2, I2=0 * Q1=Q1+Q2, Q2=0 * R1=R1+R2, R2=0 * *β*1 = *β, γ*1 = *γ* * if floor(t) % 7 = 1 AND Bol=1, we test * – Tests=E1 * – for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3 ×Fraction of workers tested, then Q1=Q1+1 E1=E1-1 * – Tests=A1 * – for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3 ×Fraction of workers tested, then Q1=Q1+1 A1=A1-1 * if floor(t) % 7 /= 1: Bol=1 3. Go to (a) or END. For the second wave, the algorithm reads 1. Set the initial condition for (*S**C*, *E**C*, *I**C*, *H**C*, *R**C*) their value at time *t* = 91 and the time *t* at 91. 2. set (*S, E, A, I, Q, R*) = (307, 0, 0, 0, 0, 7), and start the stochastic simulation. 1. Select reaction *r**i* and Δ*t* according to Gillespie algorithm 2. update (*S, E, A, I, Q, R*) and *t* accordingly. 3. integrate system of ODE 1-6 for a time step of Δ*t*. 4. if floor(t) % 7 = 1 AND Bol=1, we test * Bol=0 * Tests=E * for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3× Fraction of workers tested Q=Q+1 E=E-1 * Tests=A * for (i = 0, i < Tests, i=i+1) if random(0,1) < *ϵ*3 ×Fraction of workers tested Q=Q+1 A=A-1 5. if floor(t) % 7 /= 1, then Bol=1 3. Go to (a) or END. To avoid large step-sizes selected by the Gillespie algorithm, that would lead to errors in the solution of the ODE system, we added an artificial reaction Ø → Ø at a rate 100 days−1. The fit for the first wave was done by minizing the error as in Eq. 15, but comparing the cumulative quarantined workers with the observation from the testing at times *t* = 61 and *t* = 91. Similarly, the fit for the second wave was done by minizing the function at times *t* = 191, *t* = 291 and *t* = 346. ##### Work productivity In order to study how protective strategies affect work output we defined a variable, *W*, accounting for work done by healthy individuals, and given by: ![Formula][19] where *ϵ*5 indicates the proportion of individuals in quarantine that feel healthy and can work (i.e. asymptomatic) and *µ* is the productivity of a health worker when in quarantine/working from home compared to that when working at the hospital. In particular, *µ* accounts for the fact that the tasks that can be carried out from home by most healthcare workers are limited and will thus reduce their overall normal output [3]. In the case of individuals in a desynchronisation strategy, the work done is given by: ![Formula][20] In our simulations, the worked is accounted between *t* = 15 and *t* = 364. ##### Parameters of the model All the parameters of the model are given in Table 2. View this table: [Table 2:](http://medrxiv.org/content/early/2021/06/02/2021.06.01.21258124/T2) Table 2: Parameters of the models. ### 2.4 Cost analyses In order to compare the intervention strategies from an economic point of view, we estimated the cost of both testing and desynchronisation. We estimated the cost of each test to be ~USD 42 which includes the cost of swabbing material, protective masks, and extraction and qPCR reagents, and the salary of the workers doing the tests. The cost of testing is described by ![Formula][21] The cost of home-office during a desynchronisation strategy was considered to be the proportion of the monthly salary that would correspond to the un-productive period. The average monthly wage of a healthcare professional in Switzerland was obtained from the International Labour Organisation for year 2016 (USD 7065.2) [2]. ### 2.5 Data availability Anonymized data will be released upon publication. ## 3 Results ### Implementation of an internal RT-PCR-based testing protocol on healthcare worker Starting March 23, 2020, we implemented an RT-PCR-based testing pipeline as a clinical service of the department for Visceral Surgery and Medicine at the Bern University Hospital (Fig 1A). The RT-PCR was developed within our department. Consenting participants were tested in a voluntary basis in one of the two rounds of testing performed weekly, regardless of whether they presented symptoms. Our population consists of 314 persons working in contact with patients and 58 persons who were not working in contact with patients who all attended at least five sessions or were identified positive. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/02/2021.06.01.21258124/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/06/02/2021.06.01.21258124/F1) Figure 1: Results of the testing pipeline in asymptomatic healthcare workers (A) Number of tests (left y-axis) and number of positive tests (right y-axis) for two different thresholds of the Ct value of the PCR test, 32 and 40, for the healthcare workers in contact with patients. (B) Number of official positive cases, official hospitalized cases, and positive cases for the healthcare workers. (C) Kaplan-Meier curves of infected workers healthcare workers distinguishing healthcare workers and workers without contact with patients. (D) Boxplot showing the Ct value observed in the infected healthcare workers that reported having symptoms and the healthcare workers without symptoms. Throughout the study period, the number of weekly tests performed among the 314 healthcare workers in contact with patients is shown in Figure 1A, corresponding to a weekly participation of ~65% in the first wave and a participation of ~45% during the second wave and until end of February. Included individuals were aged 18-64 and working in different healthcare roles, including medical doctors, nurses, medico-technical and laboratory staff, in potential direct contact with COVID-19 patients or samples. At the end of the reported study and considering an RT-PCR Ct cut-off of 32 (see Methods), 15 individuals in potential direct contact with COVID-19 patients or samples tested positive (Figure 1A, red bars), and with a cut-off of 40, positive tests were observed in 30 individuals. Positive cases were informed, were retested by the official testing protocol, and sent to self-isolation even if asymptomatic. The peaks of infected healthcare workers occurred with a two week delay to the peaks in the population in both waves of the COVID-19 epidemic in the Canton Bern (Figure 1B). Incidence of infections was higher in healthcare workers with patient contact compared to personnel without patient contact (Figure 1C). Low level positive results (Ct value above 32) matched with the incidence of high level positive results (Ct value below 32). High Ct values were associated with absence of symptoms (Figure 1D). ### An epidemiological model that integrates city and hospital transmission dynamics reveals an underestimation of case incidence We have developed a transmission model that describes the infection dynamics as observed in the Swiss Canton of Bern (City Model, figure 2A), and integrated the additional infection sources to which a group of healthcare workers might be exposed (Hospital Model, figure 2A). In the City Model, individuals are exposed to infection at a rate *α* and become infected after a latent period. A proportion of them are hospitalised, and both hospitalised patients and infected individuals can recover and gain immunity. In the Hospital Model we additionally considered two sources of infection: infected asymptomatic coworkers (given by *β*) and hospitalised patients (given by *γ*). After infection, we considered both asymptomatic and symptomatic as infected individuals, we considered deaths only from hospitalised individuals and we did not consider reinfection. Moreover we included a variable *Q* accounting for individuals who tested positive in the testing pipeline implemented within the department and went into quarantine and represent the fraction of detected cases before presenting symptoms. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/02/2021.06.01.21258124/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/06/02/2021.06.01.21258124/F2) Figure 2: Model calibration and integration of test results. (A) Sketch of the multi-scale epidemiological model accounting for the dynamics in the city and in the hospital. (B) Number of official cases in the canton of Bern compared with the predicted new cases shown by our model. (C) left panel: Error function (see methods) of the fit for wave 1, right panel: Observed (black curves) and simulated (50 Gillespie simulations, purple curves), and average of 1000 simulations of the 1st wave (blue curves). *γ* = 14.08, *β* = 0, *a* = 1 (D) Error function of the fit for second wave displayed as in (C). *γ* = 0, *β* = 0, *a* = 1. (E) Error function of the fit for second wave for different values of *a*. (F) Observed and simulated (50 Gillespie simulations), and average of 1000 simulations for the second wave for different values of *a*. (left)*γ* = 1.12, *β* = 0.28, (middle)*γ* = 0.24, *β* = 0.24, (right)*γ* = 0.08, *β* = 0.12. The model was fitted to the number of hospitalised patients recorded in the Swiss Canton of Bern [18] (Figure 2B), which is in line with those reported in other Swiss regions [15], as well as the infection rate in the city, *α*, which is time-dependent as a results to the implemented policies. We inferred the case incidence in the region (Figure 2B), and showed that the reported data on incidence is an underestimate, likely due to official cases accounting only for test-confirmed cases. Our model predicts a true incidence ~4.0 times higher than confirmed and reported cases during the first wave from March to May 2020 and up to ~ 1.7 times higher from June 2020 to March 2021 (Figure 2B). ### Household and community transmission alone does not explain case incidence among healthcare workers The model allowed to independently study the sources of infection that healthcare personnel might be exposed to, namely infection from the community, parameters *α* and *a*, from co-workers, *β*, and from hospitalized patients, *γ*. In order to estimate the values of *β* and *γ* that can explain the actual cases detected, we fitted the detected cases to *Q* (Quarantine variable) of our model. Of note, the fraction of detected cases corresponds to the number of pre- and asymptomatic cases detected among the workforce tested through our pipeline –i.e. the *E*, and *A* variable in our model– and does not consider individuals who tested positive at other testing centers when already presenting symptoms (i.e. the *I* variable in the model). The fit was performed for the two waves of the COVID-19 pandemic separately (Figure 2C and 2D). The results reveal that in the first wave relatively high values of *β* or *γ* are necessary to explain our testing data best (Figure 2C left panel). Figure 2C right panel shows the fit using the values of *γ* and *β* that result into the best fit obtained from 2C left panel. In the second wave, the low incidence detected in our pipeline could be explained with low values of *β* and *γ* (Figure 2D). The best fit is shown in 2D right panel. Transmission by healthcare workers therefore seemed to be lower in the second compared to the first wave. We tested if the protection from infection in healthcare workers outside the hospital is different compared to the general population and therefore would impact the results. Thus, we decided to study reduced infections from the community by reducing *a*. In Figure 2E, we observe that for a smaller values of *a*, an increase in *β* and *γ* is necessary to explain the infection rate of healthcare workers and that with decreasing *a*, the stochastic effects become more important. The best fits are shown in Figure 2F. Therefore, modelling based on the results of the continuous testing of healthcare workers allowed to determine the added risks of infection at work for the two waves, and we inferred the transmission rates at the hospital that explain the observed infections. ### Regular testing and desynchronisation prevent infection transmission between coworkers We modelled and studied two intervention measures to reduce infection transmission between coworkers: regularly testing a fraction of individuals, and desynchronisation, where the workers are split in two teams [19]. Figure 3 shows the number of infected workers as a function of *a* (Figure 3A), *β* (figure 3B) and *γ* (Figure 3C) for different strategies, namely none, testing, desynchronisation and testing+desynchronisation. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/02/2021.06.01.21258124/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/06/02/2021.06.01.21258124/F3) Figure 3: Effect of the protective measures. (A) Number of infected healthcare workers at the end of a simulation for different values of *a* for the different strategies (none, test, desynchronisation and test+desynchronisation); *β* = 0.15, *γ* = 2. (B) Number of infected healthcare workers at the end of a simulation for different values of *β* for the different strategies; *γ* = 2, *a* = 1. (C) Number of infected healthcare workers at the end of a simulation for different values of *γ* for the different strategies; *β* = 0.15, *a* = 1. (D) Fraction of detected cases in as a function of the fraction of workers being tested for different intervals of testing; *β* = 0.15, *γ* = 2, *a* = 1. Parameters *β* and *γ* are in units of days−1. test: testing strategy; des: desynchronization strategy, des+test: combining testing and desynchronization strategy The results from stochastic simulation on the effect of infection from the city, *a*, for fixed values of *γ* and *β*, are shown in Figure 3A. As expected, the number of total infected workers increase with *a*. However, we observe no remarkable differences in the behavior of the strategies. Similarly, the number of total infected workers also increases with *γ* (Figure 3C). In Figure 3B we studied the effect of infection from co-workers *β*. In the case of no infection from co-workers (*β* = 0), testing has no impact, while desynchronisation reduces the number of cases, as a result from decreased infections from hospitalized persons. As *β* increases, the number of cases and the protective effects of the measures increases. Because stochastic effects become important due to the small number of healthcare workers, the variability increases in parallel. With very large values of *β*, e.g. 0.6 this increase is reduced because herd immunity may be achieved. In Figure 3D we show the effect of the frequency of testing and of the fraction of healthcare workers (fraction of testing) that is being tested. Regimes that are included range from testing every 3 days to testing every 14 days. If a specific fraction of persons is not getting tested (e.g. because of logistic reasons) can be compensated by increasing the frequency of testing. Together, our results showed that the implementation of regular testing of healthcare workers and workforce desynchronisation are effective in preventing worker-worker transmission, and that active adherence to frequent testing is crucial for its efficiency in detecting pre-symptomatic cases. ### Regular testing outperforms desynchronisation in work output and cost-effectiveness in Switzerland Both regular testing and desynchronisation are efficient in reducing transmission between healthcare workers (Figure 3). Nevertheless, both strategies come at a cost: an impact on the overall productivity and the cost of the tests. Most healthcare workers - notably doctors, nurses, medico-technical and laboratory staff - can carry out limited work at-home. Therefore, a notable decrease in the work output is the consequence of home-office that will be added to that of infected workers. In order to account for work output of the healthcare cohort, we added a work variable, *W*, to the model that accounts for the work output of *S, E, A, R* and *Q* based on their productivity while working at-home (see Methods). In the regular testing strategy, total work output is reduced due to workers infected and in quarantine, as well as by the time devoted to the internal testing pipeline, which we estimated as two individuals full work time per day to carry out the testing protocol. In the desynchronisation strategy, work output is reduced due to workers infected and in quarantine, as well as by the decrease in productivity during home-office. Regular testing has a smaller impact on the overall work than desynchronisation for low home office productivity (Figure 4A). Indicating that testing would be highly efficient and recommended to reduce transmission with a minimal impact of workforce availability at different home-office productivity rates. On the other hand, weekly desynchronisation results in a stronger decrease in productivity, the magnitude of which depends on the proportion of work that is considered to be feasible at home. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/02/2021.06.01.21258124/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/06/02/2021.06.01.21258124/F4) Figure 4: Economical impact. (A) Work output at the end of the simulation as a function of the home office efficiency for the different strategies. (B) Total cost of the measures as a function of the home office efficiency for the different strategies. (C) Total cost of the measures as a function of the cost per test. horizontal lines represent the cost of desynchronization for different levels of home-office productivity. Average of 1000 Gillespie simulations for (A, B, C). test: testing strategy; des: desynchronization strategy; des+test: combining testing and desynchronization strategy Next we aim to include the cost of the test in the economic modelling. First we added the cost of the tests in US dollars (USD) and we estimated the cost of the productivity lost during desynchronisation to be proportional to the salary of the worker. For the costs of our in-house test we calculated ~42 USD per test. With a monthly average wage of 7065.20 USD, the cost of the testing is smaller than the cost of the desynchronisation, unless home office productivity approaches 100% (Figure 4B). However, if a commercially available test has to be applied then the requirements for home office productivity in the desynchronisation strategy would be lower to be cost efficient. In Figure 4C we compared the cost of testing with the cost of desynchronisation as a function of the cost per test. In summary, our results show that for a pandemic as observed in Bern, Switzerland, a regular internal testing regime would be the optimal strategy to prevent transmission among coworkers while keeping work output high and economic cost low. ## 4 Discussion We report the results of a testing regime of asymptomatic healthcare workers during the COVID-19 pandemic. Based on the results we have developed a multi-scale epidemiological model that integrates transmission dynamics inside and outside a hospital to describe risks of SARS-CoV-2 infection for healthcare workers. The results from the testing revealed that the peak incidence from healthcare workers lags 2 weeks behind the peak in the population. Patient contact is associated with increased incidence of infection in healthcare workers. This is supported by other studies showing high viral titers in asymptomatic healthcare workers [16]. We showed that despite the observation that test results with high Ct-values were associated with the absence of symptoms, such results can be used to predict the prevalence of infection in healthcare workers. Because of underreporting of official cases, especially during the first wave [4], we used hospitalization data to estimate the number of true infections in the city using publicly available data [1] that are in line with estimates in other Swiss studies [15, 4]. The results of the testing were used to build a model to determine sources of infection from inside and outside the hospital. The models show that household infection alone is insufficient to explain the rate of infections in healthcare workers. We detected that ~10% of healthcare workers in our cohort had been exposed (positive test with Ct value below 40) until the end of the study. This case incidence among healthcare workers is in line with that reported in other studies [7, 16]. However, estimation of the true case incidence is difficult given the seroprevalence of 1 percent healthcare workers in regions without symptomatic cases [17] up to ~2.8% [13] The parameters that are required to explain our observation indicate that for the first wave high infection rates from either co-workers or patients must occur. Conversely for the second wave, the model indicates a lower rate of infection from co-workers or patients. However, it remains to be determined if positive tests with low Ct values in asymptomatic healthcare workers are the consequence of infection from patients or from the community. Next, we used our model to study how both regular testing and desynchronisation protocols can prevent infection transmission between coworkers. Both strategies are effective in reducing transmission between coworkers, in line with other studies on the benefits of desynchronisation [19, 14, 21]. Our modelling results indicate that a reduced fraction of tested persons can be compensated for by increasing the frequency of testing. Next, we compared the consequences of the testing and desynchronisation strategies on work productivity and cost-effectiveness. Regular testing results in a minimal impact on workforce availability and work production, whereas a desynchronisation strategy would imply an important decrease on work output unless home-office productivity was virtually the same as in-hospital or in the case in which the cost of the tests is high. Given the limited tasks that can be carried out at home by health workers - estimated to be less than 50% [3] - and the access to testing resources by medical teams, we concluded that regular testing of healthcare personnel would be the intervention of choice in Switzerland for the current pandemic. However, efforts should be made for future pandemics that prepare for solutions that allow specific tasks to be performed remotely e.g. telepharmacy. In countries with lower income, desynchronisation would have a lower cost of implementation than testing. However, it would be necessary to asses the impact that halving the available personnel would have on the health system and on the epidemic management as a whole. A possible limitation is that we assumed an hospitalization rate of 2.5% throughout the epidemic, from which several model parameters are then estimated. Furthermore, we did not consider reinfections among the cohort of healthcare workers over the studied period, but they could be worth including in future models over longer periods where reinfection would be more likely. Although variations in the hospitalization rate and reinfections would lead to different parameters and in turn, case incidence outside and possibly inside the hospital, these do not affect the conclusions of our study regarding the effectiveness of preventive interventions. In summary, our study showed that frequent and widespread testing of pre- and asymptomatic healthcare workers is effective in detecting infections and preventing transmission between coworkers while optimising work output and cost-effectiveness. ## Data Availability Anonymized data will be released upon publication. ## UVCM-COVID researchers Andrew Macpherson, Daniel Candinas, Deborah Stroka, Monika Wegmüller, Sandra Wenger, Lia Bally, Andreas Melmer, Mirjam Kneubühl, Elke Beutler, Michelle Broger, Isabel Huber, Jeannine Kölliker, Kimberly König, Joseba Möri, Chiara Ziegler, Joana Freitas, Sophie Lagger, Elisabeth Leu, Sandra Tenisch, Nicole Truffer, Felix Alexander Baier, Patrick Brönnimann, Jacopo Gavini, Magdalena Eilenberg, Nicolas Melin, Daniel Rodjakovic, Annina Schmid, Riccardo Tombolini, Katharina Bacher, Marianne Berger Rentsch, Sophie Burkhalter, Marco Felber, Izzem Gemici, Dana Leuenberger, Sina Maletti, Sarah Maring, Jelena Murar, Philip Rubin, Daniela Sommer-Ezzis, Ziad Al Nabhani, Ian Young ## Acknowledgements The study was supported by the Swiss National Science Foundation with the grants 196059 (G. B., D. S.-T.) and 196641 (S. G.-V.) and by the Bern Center for Precision Medicine (G. B., D. S.-T.). * Received June 1, 2021. * Revision received June 1, 2021. * Accepted June 2, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].Federal office of public health foph. [https://www.covid19.admin.ch/en/overview](https://www.covid19.admin.ch/en/overview). Accessed: 2021-04-09. 2. [2].Statistics on Wages, [https://ilostat.ilo.org/topics/wages/](https://ilostat.ilo.org/topics/wages/). ILOSTAT, 2021. 3. [3]. Abi Adams-Prassl, Teodora Boneva, Marta Golin, and Christopher Rauh. Work That Can Be Done from Home: Evidence on Variation within and across Occupations and Industries. page 62, 2020. 4. [4]. Christian L. Althaus, Daniel Probst, Anthony Hauser, and Julien L. Riou. Time is of the essence: containment of the SARS-CoV-2 epidemic in Switzerland from February to May 2020. medRxiv, page 2020.07.21.20158014, July 2020. Publisher: Cold Spring Harbor Laboratory Press. 5. [5]. Qifang Bi, Justin Lessler, Isabella Eckerle, Stephen A. Lauer, Laurent Kaiser, Nicolas Vuilleumier, Derek AT Cummings, Antoine Flahault, Dusan Petrovic, Idris Guessous, Silvia Stringhini, Andrew S. Azman, and for the SEROCoV-POP Study Group. Household Transmission of SARS-CoV-2: Insights from a Population-based Serological Survey. medRxiv, page 2020.11.04.20225573, January 2021. Publisher: Cold Spring Harbor Laboratory Press. 6. [6]. Oyungerel Byambasuren, Magnolia Cardona, Katy Bell, Justin Clark, Mary-Louise McLaws, and Paul Glasziou. Estimating the extent of asymptomatic COVID-19 and its potential for community transmission: Systematic review and meta-analysis. Official Journal of the Association of Medical Microbiology and Infectious Disease Canada, 5(4):223–234, December 2020. Publisher: University of Toronto Press. 7. [7]. Maria Dolores Folgueira, Carmen Munoz-Ruiperez, Miguel Angel Alonso- Lopez, Rafael Delgado, and on behalf of the Hospital 12 de Octubre COVID-19 Study Groups. SARS-CoV-2 infection in Health Care Workers in a large public hospital in Madrid, Spain, during March 2020. medRxiv, page 2020.04.07.20055723, April 2020. Publisher: Cold Spring Harbor Laboratory Press. 8. [8]. Tapiwa Ganyani, Cécile Kremer, Dongxuan Chen, Andrea Torneri, Christel Faes, Jacco Wallinga, and Niel Hens. Estimating the generation interval for coronavirus disease (COVID-19) based on symptom onset data, March 2020. Eurosurveillance, 25(17):2000257, April 2020. Publisher: European Centre for Disease Prevention and Control. 9. [9]. Crispin Gardiner. Stochastic methods, volume 4. Springer Berlin, 2009. 10. [10]. K Gaythorpe, N Imai, G Cuomo-Dannenburg, M Baguelin, S Bhatia, A Boonyasiri, A Cori, Z Cucunuba Perez, A Dighe, I Dorigatti, R Fitzjohn, H Fu, W Green, J Griffin, A Hamlet, W Hinsley, N Hong, M Kwun, D Laydon, G Nedjati Gilani, L Okell, S Riley, H Thompson, S Van Elsland, R Verity, E Volz, P Walker, H Wang, Y Wang, C Walters, C Whittaker, P Winskill, X Xi, C Donnelly, A Ghani, and N Ferguson. Report 8: Symptom progression of COVID-19. Technical report, Imperial College London, March 2020. 11. [11]. Daniel T Gillespie. Exact stochastic simulation of coupled chemical reactions. The journal of physical chemistry, 81(25):2340–2361, 1977. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/j100540a008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=A1977EE49800&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F02%2F2021.06.01.21258124.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1977EE49800008&link_type=ISI) 12. [12]. Sergio Alejandro Gomez-Ochoa, Oscar H. Franco, Lyda Z. Rojas, Peter Francis Raguindin, Zayne Milena Roa-Diaz, Beatrice Minder Wyssmann, Sandra Lucrecia Romero Guevara, Luis Eduardo Echeverria, Marija Glisic, and Taulant Muka. COVID-19 in Healthcare Workers: A Living Systematic Review and Meta-analysis of Prevalence, Risk Factors, Clinical Characteristics, and Outcomes. Am J Epidemiol. 13. [13]. Aidan T Hanrath, Ina Schim Van Der Loeff, Dennis W Lendrem, Kenneth F Baker, David A Price, Peter McDowall, Kiera McDowall, Susan Cook, Peter Towns, Ulrich Schwab, et al. Sars-cov-2 testing of 11,884 healthcare workers at an acute nhs hospital trust in england: a retrospective analysis. Frontiers in medicine, 8:239, 2021. 14. [14]. Omer Karin, Yinon M. Bar-On, Tomer Milo, Itay Katzir, Avi Mayo, Yael Korem, Boaz Dudovich, Eran Yashiv, Amos J. Zehavi, Nadav Davi- dovich, Ron Milo, and Uri Alon. Adaptive cyclic exit strategies from lock-down to suppress COVID-19 and allow economic activity. medRxiv, page 2020.04.04.20053579, April 2020. Publisher: Cold Spring Harbor Labora-tory Press. 15. [15]. Joseph C. Lemaitre, Javier Perez-Saez, Andrew S. Azman, Andrea Rinaldo, and Jacques Fellay. Assessing the impact of non-pharmaceutical interventions on SARS-CoV-2 transmission in Switzerland. Swiss Medical Weekly, 150(2122), May 2020. Publisher: EMH Media. 16. [16]. M Catherine McEllistrem, Cornelius J Clancy, Deanna J Buehrle, Nina Singh, Aaron Lucas, Valerie Sirianni, and Brooke K Decker. Sars-cov-2 is associated with high viral loads in asymptomatic and recently symptomatic healthcare workers. Plos one, 16(3):e0248347. 2021. 17. [17]. Tanawin Nopsopon, Krit Pongpirul, Korn Chotirosniramit, Wutichai Jakaew, Chuenkhwan Kaewwijit, Sawan Kanchana, and Narin Hiran-suthikul. Seroprevalence of hospital staff in a province with zero covid-19 cases. PloS one, 16(4):e0238088. 2021. 18. [18]. Daniel Probst. daenuprobst/covid19-cases-switzerland, December 2020. original-date: 2020-03-12T12:09:53Z. 19. [19]. Daniel Sanchez-Taltavull, Daniel Candinas, Edgar Roldan, and Guido Beldi. Modelling strategies to organize healthcare workforce during pandemics: application to COVID-19. medRxiv, page 2020.03.23.20041863, June 2020. Publisher: Cold Spring Harbor Laboratory Press. 20. [20]. Kevin L. Schwartz, Camille Achonu, Sarah A. Buchan, Kevin A. Brown, Brenda Lee, Michael Whelan, Julie HC Wu, and Gary Garber. COVID-19 infections among Healthcare Workers and Transmission within Households. medRxiv, page 2020.06.12.20129619, June 2020. Publisher: Cold Spring Harbor Laboratory Press. 21. [21]. Chun Yee Lim, Mary Kathryn Bohn, Giuseppe Lippi, Maurizio Ferrari, Tze Ping Loh, Kwok-Yung Yuen, Khosrow Adeli, and Andrea Rita Horvath. Staff rostering, split team arrangement, social distancing (physical distancing) and use of personal protective equipment to minimize risk of workplace transmission during the COVID-19 pandemic: a simulation study. Clinical Biochemistry, September 2020. 22. [22]. Charlotte Zheng, Nema Hafezi-Bakhtiari, Victoria Cooper, Harriet David-son, Maximillian Habibi, Peter Riley, and Aodhan Breathnach. Characteristics and transmission dynamics of COVID-19 in healthcare workers at a London teaching hospital. medRxiv, page 2020.07.10.20149237, July 2020. Publisher: Cold Spring Harbor Laboratory Press. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-6.gif [7]: /embed/graphic-7.gif [8]: /embed/inline-graphic-1.gif [9]: /embed/graphic-8.gif [10]: /embed/graphic-10.gif [11]: /embed/graphic-11.gif [12]: /embed/graphic-12.gif [13]: /embed/graphic-13.gif [14]: /embed/graphic-14.gif [15]: /embed/graphic-15.gif [16]: /embed/graphic-16.gif [17]: /embed/graphic-17.gif [18]: /embed/graphic-18.gif [19]: /embed/graphic-19.gif [20]: /embed/graphic-20.gif [21]: /embed/graphic-22.gif