ABSTRACT
On November 26, 2021, the B.1.1.529 COVID-19 variant was named as the Omicron variant of concern. Reports of higher transmissibility and potential immune evasion triggered flight bans and heightened health control measures across the world to stem its distribution. Wastewater-based surveillance has demonstrated to be a useful complement for community-based tracking of SARS-CoV-2 variants. Using design principles of our previous assays that detect SARS-CoV-2 variants (Alpha and Delta), here we report an allele-specific RT-qPCR assay that simultaneously targets mutations Q493R, G496S and Q498R for quantitative detection of the Omicron variant in wastewater. This method is open-sourced and can be implemented using commercially available RT-qPCR protocols, and would be an important tool for tracking the spread and introduction of the Omicron variant in communities for informed public health responses.
INTRODUCTION
The COVID-19 disease was first detected late December 2019 and rapidly spread globally, leading to the WHO declaring it a global pandemic by March 2020 (WHO, 2020). The etiological cause of this disease is the single-stranded RNA virus SARS-CoV-2 of the genus Betacoronavirus which also contains other human respiratory pathogens such as SARS-CoV and MERS-CoV (Pal et al., 2020). SARS-CoV-2 infection in the respiratory and gastrointestinal tract is mediated by the binding between viral spike protein (S) and human angiotensin-converting enzyme 2 (ACE-2) receptor (Zhou et al., 2020). Over the course of the pandemic, multiple SARS-CoV-2 variants emerged due to genetic mutation of the viral genome. While many variants contain inconsequential mutations, some obtained mutations that confer higher fitness through increased transmissibility and the ability to elude medical countermeasures such as vaccines (Harvey et al., 2021). These variants are deemed by the WHO as variants of concerns (VOCs) due to their increased risk (CDC, 2021). To date, WHO has classified five variants as VOCs. These include Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1) and the Delta (B.1.617.2) variant, which emerged from India in October 2020 and quickly became the predominant strain by mid-2021 (WHO, 2021a). The fifth VOC, Omicron (B.1.1.529) and descendent lineages were reported in multiple countries in November 2021. The Omicron variant draws attention due to its high number of mutations, with 26-32 of them in the spike protein (WHO, 2021b). At the time of writing, the Omicron variant only makes up 10% of all COVID-19 cases but their rapid expansion in the South African region and global distribution amidst early reports of higher transmissibility and immune escape compared to the Delta variant has raised alarm for researchers and public health officials alike (Pulliam et al., 2021).
In light of the emergence of VOCs, surveillance efforts tracking their introduction and spread in both naive and vaccinated populations become important to fight the pandemic. Current widely-used method for variant surveillance involves sequencing of clinical samples. However, the accessibility of genomic sequencing is limited and unsustainable due to the high capital and operational costs and the specialised infrastructure requirements associated (Gwinn et al., 2019). Furthermore, a significant fraction of confirmed cases needs to be sequenced to generate a meaningful dataset. A complementary method of surveillance that has been drawing a lot of attention during this pandemic is wastewater-based surveillance (WBS). WBS promises a low-cost, real-time, and unbiased snapshot of entire populations within its catchment range (Polo et al., 2020; Thompson et al., 2020). It has been shown to be effective at determining SARS-CoV-2 circulation trends during this COVID-19 pandemic across different countries (Medema et al., 2020; Randazzo et al., 2020; Wu et al., 2021). Further, WBS has been used to track the introduction and spread of SARS-CoV-2 strains. Most of the efforts for variant tracking in wastewater rely on enriching and sequencing the environmental SARS-CoV-2 genome (Crits-Christoph et al., 2021; Fontenele et al., 2021; Napit et al., 2021). However, poor sensitivity towards low-frequency variants and the lack of quantitative modelling of the data generated limits widespread application (Van Poelvoorde et al., 2021). RT-qPCR-based methods targeting variant-specific mutants of the SARS-CoV-2 genome have been demonstrated for variant identification in clinical samples (Wang et al., 2021) and we and others have recently adapted and validated use of such assays for quantification of multiple variants in wastewater samples (Graber et al., 2021; Lee et al., 2021a, 2021b, 2021c; Yaniv et al., 2021). RT-qPCR methods enable differentiation of specific variant-linked mutations and are more sensitive than sequence-based approaches, allowing for their quantitation, and providing readily interpretable results within hours. However, designing assays for VOCs is a moving target since emerging VOCs may possess unique mutations that necessitate the development of new assays.
Our group has previously designed RT-qPCR-based methods for the detection of mutations associated with the Alpha variant and the Delta variant in wastewater (Lee et al., 2021a, 2021b, 2021c). The assays are designed based on principles of allele-specific (AS) qPCR (Petruska et al., 1988; Wu et al., 1989). Here, we develop a new set of AS RT-qPCR primers for the detection of the Omicron variant, simultaneously targeting loci Q493R, G496S and Q498R, which are currently specific for Omicron variant. This newly designed assay can be run as a singleplex or as a duplex with an earlier assay designed to target HV69-70del (Lee et al., 2021b, 2021c) for improved convenience.
RESULTS AND DISCUSSION
In this work we explore and develop allele-specific (AS) RT-qPCR primers that could be used to detect and quantify the Omicron variant in wastewater. We have developed two approaches for Omicron detection, including (1) a new assay that is currently specific for mutations present in the spike gene in Omicron (Figure 1), and (2) targeting the locus 69-70 (HV69-70del) in Omicron, also in the spike gene, that is shared with the Alpha variant (GISAID, 2021a) (Lee et al., 2021b, 2021c). In the new Omicron specific assay, we simultaneously target Omicron mutations Q493R, G496S and Q498R. This stretch of mutations are present in most of the sequences (97.5%) classified as Omicron, and can be used to differentiate Omicron from other variants (Table 1). Initial validation with DNA and RNA templates demonstrate good sensitivity and specificity of this assay.
SARS-CoV-2 sequences obtained within the last two months (19th Oct 2021 - 19th Dec 2021) were obtained from the GISAID database. The percentage of Omicron and wild-type (WT) (which includes all non-Omicron) genomes containing each target mutation were used as a gauge to determine reliability of assay at identifying the Omicron variant.
(a) Each assay requires a common reverse primer (black) and probe (pink), paired with an Omicron-specific (dark blue) or WT (light blue) primer (dark blue) in the forward direction. Those allele-specific primers are designed to target the mutations of interest (red), with the mutations near the 3′ end to enhance assay specificity.
Primer design
Specificity and cross-reactivity of AS RT-qPCR primers against DNA containing WT and Omicron sequences
Assays for locus 69-70 were validated in an earlier work (Lee et al., 2021b, 2021c). Here we describe the validation of the approach that simultaneously targets loci Q493R, G496S and Q498R. The AS RT-qPCR assay developed was validated using synthetic DNA containing WT and Omicron sequences at these loci. We demonstrate the specificity for their respective WT and mutant targets at locus 493-498 by screening them against synthetic DNA containing sequences representing Omicron and non-Omicron (WT) genotypes (Figure 2). Cross-reactivity was observed at only 104 copies of DNA of the opposite genotype suggesting good specificity and selectivity towards the respective genotypes. (Figure 2). Amplification efficiencies were between 90% and 94% for the correct DNA (i.e., Omicron assay for DNA sequences representative of Omicron and WT assay for WT sequences).
Purple diamonds and squares represent tests against the matching genotype (WT-specific primers to WT DNA and Omicron-specific primers to Omicron DNA), and grey diamonds and squares denote tests against DNA of the opposite genotype. The presented data reflect two sets of independent measurements taken on different days.
Besides the use of DNA templates, we tested our assay using synthetic WT RNA. As commercial Omicron variant RNA was not available in the laboratory at the time of this work, the primer-probe set developed was validated only against synthetic WT RNA (represented by Delta variant) (Figure 3). We use the performance (specificity and sensitivity) of the WT assay against WT RNA as a proxy to gauge how the Omicron assay would perform against Omicron RNA, given that the WT assay differs from the Omicron assay by only three nucleotides. The amplification efficiency for the WT RNA is 85.2% and 100% in two independent tests, respectively. The Omicron-specific assay does not cross-react with WT RNA at 103 copies and below and confers sufficient specificity for our work given that the number of copies of SARS-CoV-2 RNA in each reaction containing RNA template from wastewater is typically below 1000 (Duvallet et al., 2021; Wu et al., 2022, 2021, 2020). Further, using cycle thresholds (Ct) as a proxy for the sensitivity of the assay (Figure 4), we compared Ct values to those of the U.S. CDC N1 assay and found comparable sensitivities for the same amount of input RNA across a 10,000-fold difference in target RNA concentrations.
Purple diamonds represent tests against the matching genotype (WT-specific primers to WT RNA), and grey squares denote tests against RNA of the opposite genotype. The presented data reflect two sets of independent measurements taken on different days.
Comparison of Ct values suggests similar sensitivity of AS RT-qPCR primers (purple diamonds) in comparison to US CDC N1 primers (blue diamonds). Primer-probe sets were tested with tenfold dilutions of their respective full-length synthetic WT SARS-CoV-2 RNA (represented by Delta RNA) in nuclease-free water. The data shown reflect two sets of independent measurements taken on different days.
DISCUSSION
The study presented here proposes an AS RT-qPCR assay that has been developed for specific detection and quantitation of the SARS-CoV-2 Omicron variant of concern in wastewater samples. This assay simultaneously targets loci Q493R, G496S and Q498R. Both the WT (non-Omicron) and Omicron assays are highly specific and do not cross-react with DNA of the opposite genotype up to 103 and 104 copies, respectively. When tested against SARS-CoV-2 WT RNA, the WT assay has shown good amplification efficiency of above 85% and appears to have comparable sensitivity as the widely used US CDC N1 assay. The Omicron assay does not show non-specific cross reactions with WT RNA in concentrations of up to 103 copies. In all, the assay targeting locus 493-498 should possess sufficient sensitivity and specificity for accurate detection of the Omicron variant in wastewater.
We also mention that an allele-specific qPCR assay which has been previously developed and described for the Alpha variant (Lee et al., 2021b,c) could be used to indicate presence of the Omicron variant. This deletion at locus 69-70 is found in the Omicron variant and in areas with low prevalence of the Alpha variant, can be used to track dynamics of the Omicron variant, corroborating results from the Omicron-specific assay. While there may be ambiguity if a positive read out with this assay would signify the Alpha or the Omicron variant, as of late Dec 2021, less than 0.2% of newly reported sequences is attributed to the Alpha variant (GISAID, 2021b). This assay targeting 69-70 is highly specific and does not cross react with up to 105 copies of RNA of the opposite (WT) genotype (Lee et al., 2021b,c), while possessing comparable sensitivities to widely accepted US CDC N1/N2 assays. While not demonstrated in this work, assay 69-70 can be combined with the 493-498 assay described here, into a duplex for greater convenience of handling. Potential duplex combinations are shown in Figure 5.
The assay for the locus 69-70 has been developed in an earlier study on Alpha variant quantification (Lee et al., 2021a, 2021b). Two individual RT-qPCRs are needed to assay for variant frequencies at each locus, though the two assays can be combined in the same reaction as a duplex (dark blue box) and the two WT-specific assays in another duplex (light blue box).
While variant detection only requires performing the reaction that targets the mutated locus in the SARS-CoV-2 variant of concern, we present assays that target both the WT and the mutant at the same loci. This way, quantification of both WT and mutant loci in wastewater allows us to determine the proportion of WT to the variant sequence at the target loci. Further, since at the time of this work, commercial Omicron variant RNA was not yet available in our laboratory, the performance of the Omicron assay is approximated using the WT assay as a proxy. WT and Omicron assays differ by only three nucleotides and hybridization kinetics are likely to be similar. This, however, constitutes a caveat for this work. Assays for the Omicron variant were validated on synthetic DNA carrying WT and Omicron sequences, and against synthetic WT SARS-CoV-2 RNA but not Omicron variant RNA. Further, these assays have not yet been validated against Omicron sequences in real wastewater samples. However, given the rapid increase in circulation of the Omicron VOC, we are sharing these assays and their initial validations for the benefit of the wider wastewater surveillance community. Detailed analytical validation against the Omicron variant RNA will be performed and reported as an update to this preprint.
MATERIALS AND METHODS
Determination of frequencies of targeted mutations in Omicron and (WT) non-Omicron strains
A global masked alignment of SARS-CoV-2 sequences from GISAID aligned to the GISAID reference sequence WIV04 (Zhou et al., 2020) were downloaded on GISAID (19th Dec 2021). We analysed sequences from the recent two months (19th Oct 2021 - 19th Dec 2021). Among 1,082,005 sequences, 0.81% sequences were assigned as omicron strains (B.1.1.529 + BA.*) and 99.2% were assigned as non-omicron strains.
Assay design
We designed AS RT-qPCR reactions to detect a stretch of mutations in the SARS-CoV-2 spike gene. Primers and probes were designed following our previous work (Lee et al., 2021) and using the Integrated DNA Technologies (IDT)’s PrimerQuest Tool. Target mutations were placed near the 3’ end of the forward primer. All primers were designed to have a melting temperature in the range of 59–65°C and the probes in the range of 64–72°C. Probes were designed to anneal to the same strand as the allele-specific primer, with the probe as close to the 3′-end of the Allele Specific (AS) primers as possible. Guanines are avoided at the 5′-end of the probe. WT primers are designed to bind to all non-Omicron sequences. The probe for 493-498 is designed with mixed bases at two positions to enable binding to both WT (non-Omicron) and Omicron sequences. All primers and probes were purchased from IDT (Table 2).
Allele-specific nucleotides are marked bold. Gen denotes general (common primers), and are used in both WT and omicron-specific reactions. Om denotes Omicron, and is used in Omicron-specific reactions. WT denotes wild type, and is used for WT (non-Omicron)-specific reactions. If the assays were to be duplexed, one of the probes should be synthesised with a different fluorophore i.e. HEX.
Analysis of assays against DNA standards by qPCR
qPCR reactions were performed using TaqMan Fast Advanced Master Mix (ThermoFisher #4444556) at a final volume of 10 μL in duplicates, according to the manufacturer’s recommendations with a final primer concentration of 0.5 μM and probe of 0.2 μM with 1 μL of DNA template, at 20 s at 95 °C, 40 cycles of 3 s at 95 °C and 30 s at 60 °C. A single reverse primer and probe was used with each allele-specific forward primer (Table 3).
Yellow highlight demarcates binding of AS forward primer. Pink highlight demarcates binding of probe. Green highlight demarcates binding of reverse primer. Bold underlined denotes positions where nucleotides differ between the two sequences.
RNA Standards and their quantification by RT-ddPCR
Twist synthetic SARS-CoV-2 RNA control 23 (Delta, B.1.617.2) was used as the RNA standard, as it contains the wild-type (WT) sequence at the three targeted mutant loci. RNA standards were prepared as single-use aliquots. Control 23 was quantified by digital droplet PCR (ddPCR) to be 3.96 × 105 copies/μL. Quantification was performed using One-Step RT-ddPCR Advanced Kit for Probes #1864022 (Bio-Rad) following manufacturer’s recommendations.
Analysis of assays against RNA standards by RT-qPCR
AS RT-qPCR was performed using the Taqman Virus 1-Step master mix (Thermofisher #4444434) with technical duplicates, at a final volume of 10 µL, according to the manufacturer’s recommendations. A single reverse primer and probe was used with each allele-specific forward primer (Table 3). The final concentration of the AS RT-qPCR primers were 500 nM, probe at 200 nM, with 1 µL of template. No template controls were included for each assay and none of them amplified. The reactions are setup using electronic pipettes (Eppendorf) and performed on a Bio-Rad CFX384 real-time PCR instrument under the following conditions, 5 min at 50 °C and 20 s at 95 °C, followed by 45 cycles of 3 s at 95 °C and 30 s at 60 °C.
Data analysis
Frequency of targeted mutations was analysed using Python. qPCR data was analysed using Microsoft Excel and Graphpad Prism. Graphs were presented using Graphpad Prism.
Data availability
Source data will be made available upon request.
Declaration of competing interests
EJA is an advisor to Biobot Analytics and holds shares in the company.
Funding Statement
This research is supported by the National Research Foundation, Prime Minister’s Office, Singapore, under its Campus for Research Excellence and Technological Enterprise (CREATE) program funding to the Singapore-MIT Alliance for Research and Technology (SMART) Antimicrobial Resistance Interdisciplinary Research Group (AMR IRG) and the Intra-CREATE Thematic Grant (Cities) grant NRF2019-THE001-0003a to JT and EJA and funding from the Singapore Ministry of Education and National Research Foundation through an RCE award to Singapore Centre for Environmental Life Sciences Engineering (SCELSE).
Contributions
EJA and JT conceptualized the project. WLL designed the experiments. WLL, XG, FA, FC, HC, HC, FW, AX, ML, FJDC and GWCK analyzed the data. WLL, JYRT and CYJL performed experiments. All authors contributed to writing the manuscript. WLL, JT and EJA supervised the project. All authors read and approved the manuscript.
Acknowledgements
We thank members of the Biobot Analytics, Inc team for helpful discussions.
REFERENCES
Subject Area
Reviews and Context
0
Comment
0
TRIP Peer Reviews
0
Community Reviews
0
Automated Services
1
Blogs/Media
Author Videos