Subjects were recruited between 2nd and 26th June 2020 from two sites—a 5400-bed purpose-built dormitory where migrant workers were housed in large rooms holding 7–20 workers, and a community care facility (CCF) where migrant workers diagnosed with COVID-19 but not requiring acute hospital care were sent for isolation and monitoring. All subjects at the CCF are prior confirmed cases (via RT-PCR), while subjects from the dormitory comprised two groups—(1) migrant workers presenting with symptoms of acute respiratory tract infection (ARI); and (2) asymptomatic roommates of newly diagnosed COVID-19 cases.
This study was approved by the Director of Medical Services, Ministry of Health, under Singapore’s Infectious Disease Act17. Under this Act, in the event of a major outbreak, the Director may require the obtainment of such information or samples (including human samples) as deemed appropriate or necessary that will be of significant public health benefit to the country17. Informed consent was obtained from all participants, and all methods were performed in accordance to Singapore guidelines and regulations for biomedical research.
Migrant workers from the purpose-built dormitory presenting with ARI were assessed by physicians at the medical post, who made the decision for whether diagnostic NP swabs for COVID-19 testing was necessary. Those workers requiring NP swabs were immediately approached for study participation, and consent was taken where agreeable.
For the collection of SN swabs, participants were instructed to insert the swab (about 1 cm) into their nostrils (one at a time), tilt their head back slightly, and rotate the swab in a circular motion for 3 times around the nasal wall. The swab was then inserted into the collection tube. For the collection of naso-oropharyngeal saliva samples, participants were asked to tilt their head back slightly, clear their throat and nose, and spit the saliva into the collection bottle. The steps were repeated until the required volume (2 mL) was achieved. For “naso-oropharyngeal” saliva collection, instructional videos in the major native languages of the migrant workers were shown, following which these samples were collected under the supervision of a trained researcher.
For consenting subjects from CCF and asymptomatic roommates of newly diagnosed cases at the dormitory, NP swab collection procedure was performed by a trained researcher. SN swab and saliva samples were collected in the same sitting.
Each subject was tested up to three times at 2–3-days interval where possible, in order to compare the sensitivity of different samples across time. Subjects from the purpose-built dormitory who tested negative across all three samples during the first round of testing were not retested. Subjects from the CCF were not retested if all samples from the initial two rounds were negative.
NP swabs from subjects with ARI were sent dry in cooler boxes to the Singapore General Hospital (SGH) molecular laboratory as part of routine clinical testing. NP swabs and self-administered nasal swabs from other subjects were sent in 3 mL of viral transport medium, while up to 2 mL of saliva was collected in a container with 2 mL of viral RNA stabilization fluid (SAFER-Sample Stabilization Fluid, Lucence, Singapore) before transfer to Lucence. All samples were processed within the same day. Both service laboratories are the College of American Pathologists (CAP) accredited, and Lucence is CLIA-licensed.
RT-PCR at SGH was performed using the automated cobas 6800 system (Roche, Branchburg, NJ, USA) on an automated cobas 6800 system, with results inferred according to the manufacturer’s specifications. NP and saliva samples sent to Lucence Laboratory underwent RNA extraction (200 μL of the sample) (GeneAid Biotech Ltd) and were tested with a laboratory-developed RT-PCR test (CDC-LDT) based on primers published by the Division of Viral Diseases, National Center for Immunization and Respiratory Diseases, Centers for Disease Control and Prevention, Atlanta, GA, USA, while saliva and SN swabs were additionally tested using the Fortitude 2.1 kit (MiRXES, Singapore).
The analytical limit of detection of the CDC-LDT was determined to be 25 copies per reaction based on a synthetic SARS-CoV-2 genome (Twist Bioscience). Saliva was pre-processed with the addition of dithiothreitol (DTT) at 0.4–0.85% of total sample volume, vortexing, and incubation at room temperature for 15 min. Solubilization was visibly apparent post-treatment at room temperature and RNA was extracted immediately post-treatment.
A limited number of samples was selected for the initial stage of determining performance specifications for a NGS-based SARS-CoV-2 assay. Both saliva and SN swab samples were included to demonstrate compatibility of RNA extracts from samples collected in the viral RNA stabilization fluid. Thirty samples were selected including high and low viral load samples, and those that had discordant results from the two RT-PCR assays.
Six of the 30 samples were paired sets of saliva and SN swab samples from the same time point for 3 individuals, and 6 samples (4 saliva, 2 SN swabs) were collected at different time points from 3 other individuals. The remaining 18 samples comprised 10 saliva, 7 SN swab, and 1 NP swab sample from individual patients.
SARS-CoV-2 amplicon-based NGS was done using 330 primer pairs to generate amplicons (size range 130–178 bp) covering the entire virus genome (except the first 25 bases and 30 bases upstream of the final polyA tail) to establish a direct-from-sample workflow. To rule out potential non-specific amplification of other viruses related in sequence, all amplicons were verified to have limited similarity to sarbecoviruses, outside of SARS-related coronaviruses (assumed not to be present in circulation).
The threshold coverage (%) for making positive call by NGS was established by performing NGS on 5 negative (by RT-PCR) samples from this study, 11 negative NP swab samples from community testing, and 10 no-template controls (NTC). For samples with complete viral genomes (100% coverage ≥ 1 × coverage), phylogenetic analysis was performed to identify lineages based on sequence variants.
We described our data using frequencies/percentages and median/interquartile range. We assessed the comparability between sampling methods using kappa-statistic and percent agreement. STATA 13.1 (StataCorp, Texas, USA) was used for all statistical calculations.