Over the course of the pandemic, the Ministry of Health, Singapore (MOH) had over time fine-tuned and adjusted the local COVID-19 suspect case definitions along with the growth and evolution of global evidence base. Under the Infectious Diseases Act in Singapore, suspect COVID-19 cases are required to undergo medical investigation14 [Supplementary Table 1] while all pneumonia cases admitted in public hospitals will also be tested for SARS-CoV-2. Furthermore, doctors can conduct discretionary tests on patients based on clinical suspicion or epidemiological risk factors. A confirmed case of acute SARS-CoV-2 infection is defined as a person with respiratory sample positive for SARS-CoV-2 using a laboratory-based reverse transcription polymerase chain reaction (RT-PCR) test.15
On March 30, 2020, MOH identified a COVID-19 outbreak in a foreign worker dormitory in Singapore and all epidemiological investigations and outbreak containment measures were implemented under the Infectious Diseases Act, which grants the use of outbreak investigation data for analysis and evaluation.
Epidemiological investigation and public health measures
The outbreak response in the dormitory was divided into two phases based on the extent of transmission in the dormitory. In the pre-lockdown phase in the dormitory (before April 6, 2020), all laboratory-confirmed cases were interviewed to collect data on demographic characteristics, clinical symptoms, and activity patterns for the fourteen days preceding symptom onset or notification date until isolation in hospital.
Contact tracing was initiated to identify close contacts within the dormitory (persons who stay in the same room as a confirmed case) and in the work and social contexts (persons who spend at least 30 minutes within 2 metres of a confirmed case). These contacts were placed under quarantine for fourteen days from last exposure to the case at designated government quarantine facilities located outside the dormitory. The health status of all persons under quarantine was monitored daily and those who developed symptoms would undergo medical investigation as part of active case finding.
To facilitate outbreak investigations and management, details of the dormitory layout and all dormitory residents were requested from the dormitory operator. Employers of a confirmed case and the dormitory operators were also advised to monitor the health condition of all other workers and dormitory residents respectively daily and to advise any person who was unwell to seek medical attention immediately.
The dormitory was locked down on April 6, 2020. During the lockdown phase, healthcare and other workers were deployed to set up medical posts and perform active case finding via swabbing exercises, case isolation and quarantine of close contacts in situ, implementation of safe distancing measures while ensuring the welfare of the dormitory residents. All personnel involved in the dormitory operations were required to observe proper infection control measures at all times and to monitor their health. Residents who were unwell were advised to seek medical attention immediately and were investigated for SARS-CoV-2 infection.
Sequencing and phylogenetic analysis
All primary samples or residual extracted nucleic acid tested positive for SARS-CoV-2 by RT-PCR at diagnostic laboratories were forwarded to Singapore’s National Public Health Laboratory. Available samples from the dormitory cases and other samples linked to the dormitory outbreak (e.g. from the workplaces or social interactions) between April 1 and June 6, 2020 were randomly selected for next generation sequencing.
Selected residual diagnostic RNA were subjected to tiled amplicon PCR using ARTIC nCoV-2019 version 3 panel,16 where One-Step RT-PCR was performed using the SuperScript™ III One-Step RT-PCR System with Platinum™ Taq DNA Polymerase (ThermoFisher Scientific). Sequencing libraries were prepared using the Nextera XT and sequenced on MiSeq (Illumina) to generate 300 bp paired-end reads. The reads were subjected to a hard-trim of 50 bp on each side to remove primer artifacts using BBMap17 prior to consensus sequence generation by Burrows-Wheeler Aligner-MEM v0.7.17, with default settings. Only sequences with ≥ 98% genome coverage and supported by an average depth of 100x were included for phylogenetic analysis. The generated consensus sequences were shared via GISAID.18 To place these sequences into global context, we searched for closely related strains using BLASTN19 against all genomes in the GISAID database and retained representative hits with 99.99% identity and matching the time window of our sequences [Acknowledgements in Supplementary Information]. The sequences were merged with hCoV-19/Wuhan/WIV04/2019 (accession: EPI_ISL_402124) as reference and root for the tree and aligned using MAFFT (v7.427)20. The alignment was manually inspected and trimmed at the 5’ and 3’ ends using Jalview.21 A maximum likelihood phylogenetic tree was created with IQ-TREE v1.6.122 using ModelFinder23 for estimating the best fit model (TN + F in this case) and 1,000 steps of ultrafast bootstrapping24 with zero length branches collapsed in the final tree, visualized with Figtree.25
Seroprevalence survey to determine extent of undetected infections
To determine the extent of undetected infections within the dormitory, we undertook a prospective cross-sectional seroprevalence survey in a convenience sample of 7,367 dormitory residents who had no travel history 14 days prior to the onset of the first dormitory case and were not previously identified as a confirmed case of acute SARS-CoV-2 infection. Blood samples were collected between May 13 to June 1, 2020 (52 to 79 days since the onset of the first case; 37 to 70 days since the lockdown of the dormitory) were tested for SARS-CoV-2 Immunoglobin G (IgG) using either Abbott Architect SARS-CoV-2 IgG assay or Roche Anti-SARS-CoV-2 assay. For the overall seroprevalence, we computed the 95% confidence intervals (CI) for binomial proportions using Wilson’s method.26
Transmission model
Despite the strengthening of public health measures during the lockdown phase, cases continued to rise. This could be attributed to the dense contact networks and living conditions resulting in rapid transmission, the presence of pre- or asymptomatic transmissions, and variable health seeking behaviour of symptomatic persons that could have resulted in a delay or failure to isolate cases who went on to transmit the virus to others. Hence, to estimate the outbreak trajectory within the dormitory and to evaluate the effectiveness of the outbreak control measures, we used an individual-based model of COVID-19 transmission in a simulated population of 12,091 individuals (scaled based on number of individuals in a room) residing in a dormitory with a similar number of blocks, levels and rooms as the dormitory under study [Supplementary Table 2]. We assumed that the entire dormitory population was naïve to SARS-CoV-2 infection and disease transmission parameters such as the infectiousness over time and incubation period were modelled based on assumed distributions as elaborated upon in the Supplementary Information.
Model fitting
We hypothesized that a diverse range of parameters could drive similar outbreak trajectory in the dormitory. We generated 50,000 random parameter combinations containing parameters related to disease transmission, contact network of the residents within the dormitory, health seeking behaviour or the effectiveness of public health measures (Table 1).
Table 1
Range of values for each parameter in a parameter set
Parameter category | Parameter | Minimum | Maximum |
Disease transmission | Initial number of cases | 3 | 20 |
Proportion of asymptomatic cases (%) | 30 | 90 |
Relative infectiousness of an asymptomatic case (%) | 0 | 50 |
Probability of infection inside a room | 0.5 | 1 |
Probability of infection outside a room | 0 | 1 |
Contact network | Mean number of random contacts form on the same level | 0 | 20 |
Mean number of random contacts form on different levels but same block | 0 | 10 |
Mean number of random contacts form in different blocks | 0 | 5 |
Health seeking behaviour | Proportion of symptomatic cases seeking medical attention (%) | 0 | 100 |
Effectiveness of public health measures | Probability that contacts with persons on the same level remains after social distancing | 0 | 1 |
Probability that contacts with persons on different levels of the same block remains after social distancing | 0 | 1 |
Probability that contacts with persons in different block remains after social distancing | 0 | 1 |
Reduction in probability of infection outside a household | 0 | 0.5 |
Days since deployment of ground officers when reduction in probability of infection outside a household occurred | 1 | 7 |
The outbreak trajectory in each iteration of a parameter combination was fitted against (i) the number of cases and number of affected locations since the earliest observed onset date (day 1) till day 12 of the outbreak (prior to the lockdown of the dormitory) and (ii) the overall serology outcomes in persons with no recent travel history and no laboratory confirmation of SARS-CoV-2 infection tested by day 79 of the outbreak. For each parameter combination, the full outbreak was simulated ten times to generate 500,000 outputs. The fit of the model against the observed data in the early and late phases of the outbreak was determined by computing the likelihood [Supplementary Information]. Observed case counts in the lockdown phase were not used for model fitting as symptom onset dates were not routinely collected.
Outbreak interventions scenarios
Parameter combinations were assigned a weight based on the corresponding likelihood and weighted sampling of the parameter combinations with replacement was performed 10,000 times. We simulated the current outbreak scenario with interventions including the deployment of ground teams to expedite case isolation, quarantine of roommates, enforcement of social distancing and reduction in probability of infection outside the room in the lockdown phase of the outbreak [Supplementary Information]. Using the same disease transmission and contact network parameters, we also simulated alternative scenarios of the dormitory outbreak: (i) baseline scenario: only case isolation and quarantine of roommates were implemented; (ii) enhanced response and physical distancing scenario: ground teams were deployed to expedite case isolation and enforce physical distancing measures but no strengthening of measures was applied to reduce probability of infections occurring outside a room in the lockdown phase of the outbreak; (iii) modified dormitory setting scenario: dormitory layout with reduced number of residents per room and, en suite bathroom, shower and cooking facilities to ensure that persons under quarantine did not leave their rooms and movement restrictions across different levels for the remaining residents when the dormitory was under lockdown.6 Details of the respective interventions in each scenario are described in the Supplementary Information.
For each parameter combination, disease progression was tracked over time and generations (\({g}_{x}\) where the subscript \(x\) indicates the respective generation). The reproduction number of the \(x\)th generation—i.e. the ratio of cases in consecutive generations \(\left(\text{i.e. }\frac{{g}_{x+1}}{{g}_{x}}\right)\) which provides an indication of the growth of an outbreak—was determined. All analyses were done using R version 3.5.1.27 A full description of the model is available in the supplementary information.