Detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in a fourplex real-time quantitative reverse transcription-PCR assays
===============================================================================================================================================

* Mathieu Durand
* Philippe Thibault
* Simon Lévesque
* Ariane Brault
* Alex Carignan
* Louis Valiquette
* Philippe Martin
* Simon Labbé

## Abstract

The early diagnosis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is required to identify and isolate contagious patients to prevent further transmission of the coronavirus disease 2019 (COVID-19). In this study, we present a multitarget real-time TaqMan reverse transcription PCR (rRT-PCR) assay for the quantitative detection of SARS-CoV-2 and some of its circulating variants harboring mutations that give SARS-CoV-2 a selective advantage. Seven different primer-probe sets that included probes containing locked nucleic acid (LNA) nucleotides were designed to amplify specific wild-type and mutant sequences in Orf1ab, Envelope (E), Spike (S), and Nucleocapsid (N) genes. Furthermore, a newly developed primer-probe set targeted human β2-microglobulin (B2M) as a highly sensitive internal control for RT efficacy. All singleplex and fourplex assays detected ≤ 14 copies/reaction of quantified synthetic RNA transcripts, with a linear amplification range of 9 logarithmic orders. Primer-probe sets for detection of SARS-CoV-2 exhibited no false-positive amplifications with other common respiratory pathogens, including human coronaviruses NL63, 229E, OC43, and HKU-1. Given the emergence of SARS-CoV-2 variants and their rapid spread in some populations, fourplex rRT-PCR assay containing four primer-probe sets represents a reliable approach to detect multiple viral target sequences containing typical mutations of SARS-CoV-2 variants in a single reaction, allowing quicker detection of circulating relevant variants.

Keywords
*   COVID-19
*   SARS-CoV-2
*   real-time TaqMan reverse transcription PCR assays
*   molecular diagnostics
*   locked nucleic acid (LNA)

## Introduction

The human SARS-CoV-2 virus that causes the respiratory infection COVID-19 started in the Wuhan province of China in late 2019, and then spread around the world (Wu et al. 2020; Zhang and Holmes 2020; Zhu et al. 2020). This novel coronavirus belongs to the subgroup betacoronavirus and the subgenus Sarbecovirus (Chan et al. 2020; Lu et al. 2020a; Zhu et al. 2020). The SARS-CoV-2 genome is constituted of a positive single-stranded RNA with a genome size of nearly ∼29 800 to 29 900 nucleotides in length (Lu et al. 2020a). At the beginning of the outbreak, numerous patients exhibiting an atypical viral pneumonia in Wuhan City were notified to World Health Organization (Wu et al. 2020). Shortly thereafter, infected foreign residents of China and foreign travelers arriving in their countries of residence from international destinations, including visitors from Wuhan, contributed to the spread of the virus through person-to-person contact (Chan et al. 2020; Holshue et al. 2020; Oude Munnink et al. 2020). Subsequently, spreading of the virus occurred easily and rapidly from person-to-person transmission in local communities, resulting in a pandemic spread across the globe.

Although a number of vaccines against SARS-CoV-2 have been developed, they are not available to all people as their production and distribution require efficient logistic and time. Until vaccines become universally available, governments from many countries around the world strongly suggest to their citizens to maintain and sometimes enforce social distancing as they impose lockdown measures, overnight curfew and implement an obligatory quarantine in the case of international travelers flying or driving to different countries. Although rapid identification of SARS-CoV-2-infected patients is critical to trigger off their isolation, the early diagnosis of SARS-CoV-2 remains difficult for the following reasons. First, some SARS-CoV-2-infected patients exhibit no symptoms throughout the course of the infection. Second, other SARS-CoV-2-infected patients with mild symptoms could be confounded with patients who possess other kinds of atypical respiratory tract infections (Guarneri et al. 2021; Lai et al. 2020; Rivett et al. 2020; Yang et al. 2020). In response to the outbreak of COVID-19, numerous national public health agencies across the world have been proactive by setting up testing programs. Most of these agencies have recommended and used a nucleic-acid-based method for early detection of SARS-CoV-2. These nucleotide-based tests rely on the real-time reverse transcription PCR (rRT-PCR) that has been considered the “gold standard” approach for SARS-CoV-2 detection due to its hypersensitivity for correctly identifying the viral RNA genome found in nasopharyngeal swabs and sputum samples collected from infected patients (Corman et al. 2020; Esbin et al. 2020; Lu et al. 2020b).

For daily clinical specimen testing, clinical microbiology laboratories have used different primer-probe sets for detection of SARS-CoV-2 (Corman et al. 2020; Lu et al. 2020b; Nalla et al. 2020). Among them, the CDC N2 and Corman E (also called E Sarbeco) primer-probe sets have been identified to be particularly sensitive (Nalla et al. 2020). A typical singleplex reaction contains one set of primer pairs (forward and reverse) and a TaqMan probe that hybridizes to a specific targeted region of the viral genome. Short sequences within the N and E gene regions are examples of viral templates proved to be detectable with high reproducibility (Corman et al. 2020; Lu et al. 2020b). In the same master mixture, a second set of primer pairs and TaqMan probe containing a distinct fluorophore are added to target a human gene, such as the ribonuclease P gene (RNase P also denoted RP or RPP30) to monitor nucleic acid extraction and amplification (Lu et al. 2020b; Lu et al. 2014).

Since the beginning of the COVID-19 pandemic, several studies have reported genomic sequence variations of SARS-CoV-2 isolates (Leung et al. 2021; Martin et al. 2021; Rondinone et al. 2021; Tegally et al. 2021a; Tegally et al. 2021b). Among these genetic variations, a single nucleotide polymorphism (SNP) found at position 8,782 (in *orf1ab*; C instead of T) of SARS-CoV-2 genomic sequence is a hallmark of the L strain (Tang et al. 2020), whereas the unaltered nucleotide is found at the same position in the S strain. Several mutations have also been found in the coding sequence of the Spike protein (Starr et al. 2020; Yi et al. 2020). One critical mutation in which an adenine has been substituted by a guanine at position 23,403 in the genome of the Wuhan reference strain produces a mutant form of Spike containing D614G substitution (aspartic acid to glycine substitution at amino acid residue 614) (Korber et al. 2020). Another critical mutation in which an adenine has been substituted by a thymine at position 23,063 generates a mutant form of Spike harboring N501Y substitution (asparagine to tyrosine substitution at amino acid residue 501) that is found in the United Kingdom (UK) B.1.1.7 strain as well as other reported variants (Leung et al. 2021). These mutations represent few examples of the ability of SARS-CoV-2 to rapidly evolve since the beginning of the outbreak.

In this study, we report the development of novel primer-probe sets for the detection of SARS-CoV-2 variants. We use fluorogenic probes containing locked nucleic acid (LNA) for detection of single nucleotide polymorphisms (SNPs) in the genome of SARS-CoV-2 using the TaqMan-LNA rRT-PCR method. In addition to the new primer-probe sets, we have developed a fourplex rRT-PCR assay in which four sets of primer pairs and fluorogenic probes are contained in a single reaction master mixture. Furthermore, we have created and used a novel set of primers/probe that detects the ubiquitously expressed human β2-microglobulin (B2M) transcripts (instead of RNase P) as a positive control for monitoring performance of the whole procedure. These include the presence of RNAs in the collected sample, RNA extraction, reverse transcription efficacy, and real-time amplification. Thus, this multitarget rRT-PCR assay will improve our ability to narrow the differential identification of SARS-CoV-2 isolates.

## Materials and methods

### Primer and probe design

The Wuhan-Hu-1 genome reference sequence of SARS-CoV-2 (Lu et al. 2020a; Wu et al. 2020) as well as several other listed sequences of its genome (available in the repository of the National Center for Biotechnology Information) were aligned using BLASTn software to identify conserved nucleotide regions corresponding only to the SARS-CoV-2 viral genome and not found in other viral sequences. Subsequently, different primer-probe sets targeting distinct regions of the SARS-CoV-2 genome, especially within E, Orf1ab, and S genes, were designed using Primer Express 3.0.1 software. The N LSPQ primer-probe set was designed as described previously (Longtin et al. 2021). After validation of their exclusivity by comparing them to Virus Pathogen Resource (ViPR), Reference Viral DataBase (RVDB) and human GRCh38 genomic DNA/mRNA databases, the indicated TaqMan primer-probe sets (Table 1) were predicted to specifically amplify SARS-CoV-2, exhibiting the absence of non-specific homologies with other respiratory viral pathogens or human-related gene sequences. The only exception was the E Sarbeco primer-probe set that was predicted to hybridize with the E gene sequence of SARS-CoV-1 in addition to the E gene of SARS-CoV-2 (Corman et al. 2020). TaqMan probes were labeled at the 5’ end with either 6-carboxyfluorescein (6-FAM), Texas Red 615 (TEX 615), hexachlorofluorescein (HEX), or Cytiva5 (Cy5) and at the 3’ end with either Iowa Black FQ (IBFQ) or Black Hole Quencher 1 (BHQ1) (Integrated DNA Technologies, Coralville, Iowa).

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/T1)

Table 1. 
Oligonucleotide primers and fluorogenic probes used in real-time quantitative RT-PCR assays.

In the case of probes containing locked nucleic acid (LNA) residues, they were designed for detection of single nucleotide polymorphisms (SNPs) in the Orf1ab and S genes of SARS-CoV-2. Essentially, we used forward and reverse primers that hybridized with each side of the target sequence and wild-type or mutant TaqMan LNA probes that hybridized at the mutation site, which was located between the forward and reverse primers. All primers and TaqMan probes were designed to avoid self-complementarity and the formation of primer-dimer by-products and hairpins. Furthermore, all primer/probe sets were validated for their compatibility in multiplex rRT-PCR assays.

### Synthesis of RNA transcripts used as standards

gBlocks double-stranded DNA fragments encompassing the amplified region of each targeted SARS-CoV-2 region were synthesized accompanied with a T7 promoter sequence at their 5’ends (Fig. S1). A similar approach was used in the case of the DNA fragment covering the amplified region of the human β-2-microglobulin (B2M) gene. Purified gBlocks products were used for *in vitro* T7-dependent transcription, as described previously (Vijgen et al. 2005). The gBlocks DNA templates were eliminated by digestion with RNase-free DNase I for 15 min at 37°C. The RNA transcripts were purified using a MEGAclear transcription clean-up kit (Invitrogen) and were quantified spectrophotometrically at 260 nm. After RNA levels quantification, measurements of transcripts were converted to the molecule number per µl as described previously (Fronhoffs et al. 2002). To validate calculations of viral RNA copies per reaction, we have established standard curves from 10-fold dilutions of RNA isolated from known number of copies from a prototype viral RNA preparation that was commercially available (Exact Diagnostics, Fort Worth, TX; Vircell, Granada, Spain).

### Real-time quantitative reverse transcription PCR (rRT-PCR)

Reactions were performed in a 10-l reaction mixture including 7 µl extracted RNA-containing samples or RNA transcripts used as standards, 2.5 µl Reliance One-Step Multiplex rRT-PCR Supermix (4X concentrated) (BioRad), and 0.5 µl of forward and reverse primers (500 nM) and the indicated TaqMan or TaqMan-LNA probe (250 nM). One-step rRT-PCR amplification and detection were performed in a CFX96 Touch Real-Time PCR System (BioRad) under the following conditions. Reverse transcription of 10 min at 50°C was followed by PCR activation at 95°C for 10 min and 50 cycles of amplification (10 s at 95°C and 30 s at 60°C or for some indicated reactions at 63°C). In singleplex reaction mixtures, each reaction mixture contained a single primer pair and TaqMan or TaqMan-LNA probe. In the case of fourplex reactions, they contained all four sets of primers (500 nM per set) and probes (250 nM per probe) within a single reaction master mixture. Fluorescence measurements were monitored after each amplification round and the threshold cycle (Ct) value for each sample was calculated by assessing the point at which fluorescence crossed the threshold line, exhibiting an increase in fluorescence above the calculated background levels. The result was considered valid if two or more of the target-specific fluorescent signal showed the Ct value ≤ 37 cycles and all positive and negative control reactions gave a successful and no amplification, respectively.

### BioFire FilmArray assays

The collection of respiratory viral and bacterial pathogens was obtained from ZeptoMetrix (cat #NATRVP-IDI). The collection is called NATtrol Respiratory Verification Panel (NATRVP) and contains 20 vials x 0.6 ml, each containing viral and bacterial targets listed in Table 4. A master mixture that was constituted of an aliquot of each pathogen of NATRVP was injected into the BioFire Respiratory Panel 2.1 (RP2.1) Pouch in accordance with the manufacturer’s instructions for analyze using the BioFire FilmArray 2.0 System (BioFire Diagnostics) as described previously (Poritz et al. 2011). The loaded RP2.1 Pouch inserted into the FilmArray instrument contained all necessary reagents for automated nucleic acid extraction, reverse transcription, and two consecutive multiplex PCR amplification runs. Furthermore, the FilmArray instrument has the potential to undergo gene target melt curve analysis with each target in a valid run reported as detected or not detected.

## Results

### Limits of detection (LoD) with SARS-CoV-2 RNA transcripts in singleplex reactions

gBlocks DNA templates containing a 5’ T7 RNA polymerase promoter sequence were used to produce viral transcripts that corresponded to specific SARS-CoV-2 genomic regions (Fig. S1). These RNA molecules were used as standards for the generation of standard curves to determine the limit of detection (LoD) of each primer-probe set. Oligonucleotide primer pairs and fluorogenic probes are listed in Table 1. In the case of the E Sarbeco probe and primers, their nucleotide sequences were identical to those described previously (Corman et al. 2020). The other primer pairs and fluorogenic probes were designed from the SARS-CoV-2 complete genome (Lu et al. 2020a; Wu et al. 2020). In the cases of SARS-CoV-2 L and S strains, specific TaqMan locked nucleic acid (LNA) probes were designed to discriminate between the single nucleotide polymorphism (SNP) found at position 8,782 (in *orf1ab*; C instead of T) of SARS-CoV-2 genomic sequence as described previously (Tang et al. 2020). Similarly, a specific TaqMan LNA probe was synthesized to detect a guanine instead of an adenine at position 23,403 in the Wuhan reference strain (Lu et al. 2020a; Wu et al. 2020), allowing the identification of the Spike variant harboring a glycine residue instead of an aspartic acid residue at position 614 (Zhang et al. 2020). An additional mutation found in the UK strain B.1.1.7 was of interest to probe with a specific TaqMan LNA oligonucleotide to discriminate between the SNP found at position 23,063 (in *Spike*; thymine instead of adenine) where a tyrosine residue instead of an asparagine residue at position 501 was found (Kemp et al. 2021). The sensitivity of each set of primer pair and TaqMan probe was first evaluated in singleplex reactions. Tenfold serial dilutions of viral transcript ranging from 1.3 to 1.4 × 109 copies per reaction mixture were tested in triplicate by real-time quantitative RT-PCR (rRT-PCR) assays. Results were analyzed in terms of the Ct value that was defined as the threshold cycle in which a target viral sequence was first detected. Samples with Ct values ≤ 37 were considered to be positive in comparison with background cross-reactivity of the primers and probes in non-template control (NTC) reactions. The highest dilution of transcript that gave a significative Ct value was defined as the limit of detection (LoD) for a targeted RNA transcript. Results showed that LoD values ranged from 1.3 to 13 or 14 RNA transcript copies / reaction (Table 2). Linear regression curves were achieved over a 9-log dynamic range, from 1.3 to 1.3 × 109 copies or 1.4 to 1.4 × 109 copies per reaction for all nine probes, with calculated efficient values of 83.3% to 95.1% (Fig. 1).

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/T2)

Table 2. 
Limits of detection of SARS-CoV-2-specific primers and probes using synthetic RNA transcripts in singleplex rRT-PCR assays.

![Fig. 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/03/26/2021.03.23.21254196/F1.medium.gif)

[Fig. 1.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/F1)

Fig. 1. Standard curves to determine synthetic RNA copy number.
*A – I*. Tenfold serial dilutions ranging from 10 to 109 copies of Orf1ab L. Orf1ab S. Spike 614G. Spike D614. Spike N501. Spike 501Y. E Sarbeco. N LSPQ. and B2M synthetic RNA transcripts were analyzed by rRT-PCR assays. Each representative graph was generated by plotting the Ct values (y-axis) and the log of initial number of copies of synthetic transcripts (x-axis). Calculated linear correlation coefficients (R2) and percentage of amplification efficiencies are indicated for each primer-probe set. The graphs represent quantification of the results of three independent experiments.

### Limits of detection (LoD) in the case of fourplex rRT-PCR assays

Following analytical sensitivity of the primers and probes in singleplex reactions, we sought to develop a fourplex rRT-PCR assay that contained four sets of primers and probes within one reaction master mixture. Limits of detection in the indicated fourplex assays were determined to verify whether their sensitivities were comparable to those of singleplex reactions (Tables 2 and 3). Fourplex reactions that were prioritized included TaqMan probes for detection of single nucleotide polymorphisms in the genome of the SARS-CoV-2 strain that had been found in patient samples (Korber et al. 2020; Leung et al. 2021). With the goal of identifying distinct SARS-CoV-2 virus variants, different combinations of primer-probe sets were mixed within one reaction master mixture to specifically detect viral RNA genes that have undergone genetic variations. As an example, a fourplex master mixture contained all four sets of primers and probes to determine whether a SARS-CoV-2 isolate possessed viral RNA that encoded a Spike variant harboring either an Asp614 or a Gly614 (also denoted D614G mutation) (Korber et al. 2020). Similarly, additional fourplex master mixtures contained sets of primers and probes to determine whether a SARS-CoV-2 isolate belongs to the L or S strain (Tang et al. 2020). Other examples of fourplex reactions included primer-probe sets to identify the SNP found in *Spike gene* at position 23,063 (thymine instead of adenine) that results in an amino acid substitution N501Y that triggers a stronger interaction of the Spike receptor-binding domain with the human cell-surface receptor angiotensin-converting enzyme 2 (ACE2) (Leung et al. 2021; Martin et al. 2021). This multiplex approach could facilitate tracking of SARS-CoV-2 variants that carry mutations. Serial 10-fold dilutions of SARS-CoV-2 transcripts were produced as described above and were tested using different combinations of sets of primer pairs and fluorogenic probes in one reaction mixture as described in Table 3. Linear amplification was performed over a 9-log dynamic range, from 1.4 to 1.4 × 109 copies per tetraplex reaction for all combined probes (Table 3). Results showed linear regression curves that exhibited calculated efficient values of 85.1% to 110.4% (Fig. 2).

View this table:
[Table 3.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/T3)

Table 3. 
Sensitivities of SARS-CoV-2-specific primers and probes in fourplex rRT-PCR assays.

View this table:
[Table 4.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/T4)

Table 4. 
Specificities of fluorogenic SARS-CoV-2 probes with a number of other viruses and bacteria.

![Fig. 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/03/26/2021.03.23.21254196/F2.medium.gif)

[Fig. 2.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/F2)

Fig. 2. Linear regression curves to assess synthetic RNA copy number in fourplex rRT-PCR assays.
*A – G*. For each representative graph. four primer-probe sets within one reaction mixture were used to detect serial 10-fold dilutions of four synthetic RNA transcripts ranging from 10 to 109 copies per reaction. The y and x axes of each graph are the Ct values and the log of the input of synthetic RNA transcripts. respectively. For each fourplex reaction. R2 indicates calculated linear coefficients and Eff. shows percentage of amplification efficiencies for each primer-probe set that is present with three other sets of primers and probes within one reaction master mixture. Fluorogenic probe color codes are as follows: blue. E Sarbeco; red. Orf1ab L; violet. N LSPQ; black. B2M; green. Spike 614G; orange. Orf1ab S; brown. Spike D614; pink. Spike 501Y; and gold. Spike N501. The graphs represent quantification of the results of three independent experiments.

When we compared the sensitivity measured as the lowest Ct value of the TaqMan probes (Orf1ab L, Orf1ab S, Spike 614G, Spike D614, Spike N501, Spike 501Y, E Sarbeco, N LSPQ, and B2M) in singleplex reactions (Table2) and fourplex reactions (Table 3), results showed that sensitivities of these probes between singleplex and fourplex reactions were comparable. In the cases of E Sarbeco, Orf1ab L, and B2M, their limit of detection was generally 1.3 or 1.4 copies of transcript per reaction in both singleplex and fourplex rRT-PCR assays. In the cases of Orf1ab S, Spike D614, Spike 614G, Spike N501, Spike 501Y, and N LSPQ probes, a limit of detection of 13 or 14 copies per reaction was observed with much more consistency in singleplex and fourplex assays. Taken together, all designed primer-probe sets performed comparably in singleplex and fourplex assays that reproducibly detected as few as 1.4 to 14 copies of target sequences per reaction.

### Specificity of primers and probes in fourplex reactions

In the course of designing primer and probe sequences used for detection of SARS-CoV-2, we performed BLAST analyses to verify the absence of significant sequence homologies with other respiratory viruses and human genome sequences to avoid the possibility of false positive results. To test the specificity of the primers and probes that we have used in fourplex rRT-PCR reactions, we performed several assays using aliquots of the entire collection of NATtrol Respiratory Verification Panel (NATRVP) from ZeptoMetrix, which contains purified intact respiratory virus and bacteria particles that have been inactivated to make them non-infectious. These microbial particles were supplied in liquid samples containing a specialized matrix that included human cells to mimic the composition of a true clinical specimen. All SARS-CoV-2 primer and probe sets used in fourplex reactions (see Table 3 for all combinations of four sets of primer and probes per one reaction master mixture that have been tested) showed the absence of nonspecific amplification against the ZeptoMetrix NATRVP preparation that includes 19 respiratory pathogens such as influenza A and B, different types of parainfluenza, respiratory syncytial virus A, and other coronaviruses such as NL63, 229E, OC43, and HKU-1 (Table 4). In contrast, the B2M primer-probe set yielded positive results for the same biological samples due to the presence of human cells that are associated with the ZeptoMetrix NATtrol Respiratory pathogen samples. As positive assays for nucleic acid detection of the 19 respiratory pathogens of the ZeptoMetrix NATRVP, we performed BioFire tests using the BioFire FilmArray 2.0 system (BioFire Diagnostics, Salt Lake City, UT). A master mixture containing all viral and bacterial particles of ZeptoMetrix was injected into the BioFire Respiratory Panel 2.1 Pouch that contained all the necessary reagents for automated nuclei acid extraction, reverse transcription, two stages PCR amplification, and detection of multiple respiratory pathogen targets in a single assay. Results showed that the FilmArray instrument in combination with BioFire Pouch reagents including multiple independent sets of primers allowed multiplex-tandem PCR detection of the microbial panel members listed in Table 4. Therefore, FilmArray runs validated the presence of viral and bacterial nucleic acids in the ZeptoMetrix NATRVP preparation. Taken together, the results showed that fourplex assays are specific for the SARS-CoV-2 virus target and that there is an absence of false-positive signals with other respiratory viral or bacterial pathogens.

### LNA probes allow differential identification of specific SARS-CoV-2 transcripts in fourplex reactions

We tested different groups of differentially labeled four-color sets of probes that allow discrimination of single-nucleotide mutations found in genetic variants of SARS-CoV-2 to validate specific LNA probe-target interaction for SNPs detection in the genome of SARS-CoV-2 (Table 5). Results showed that the Orf1ab L probe containing a “C LNA” at the polymorphic site (position 8782 in SARS-CoV-2 sequence; Table 1) generated a positive fluorogenic signal that was specific for transcripts produced from the new circulating L variant strain, whereas no fluorogenic signal was detected in the presence of RNA obtained from the older S variant strain (Table 5; *e.g*. probes in fourplexes 1 and 2). In the case of the Orf1ab S probe containing a “T LNA” at the polymorphic site (Table 1), this probe exhibited a positive amplification signal in the presence of transcripts produced from the S variant strain and failed to detect transcripts of the L variant strain (Table 5; *e.g*. probes in fourplexes 3 and 4).

View this table:
[Table 5.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/T5)

Table 5. 
Real-time genotyping with oligonucleotide probes containing locked nucleic acids (LNA) for detection of SARS-CoV-2 variants.

Two oligonucleotide probes containing LNA residues were created for specific hybridization of target sequences containing nucleotide substitutions that confer amino acid substitution in Spike at Asp614 (also denoted D614) and Gly614 (also called 614G) (Table1). The LNA probe for the detection of the Wuhan reference D614 strain carried an “A LNA” at the SNP site (position 23,403 in the viral genome), whereas the fluorogenic probe specific for the detection of the circulating Spike variant (614G) contained a “G LNA” at this position (Table 1). In the case of the Spike D614 LNA probe, results showed that it was specific for the detection of transcripts from the Spike D614 gene and did not generate amplification signal in the presence of transcripts of the Spike 614G variant gene (Table 5; *e.g*. probes in fourplexes 3 and 5). In contrast, results showed that the Spike 614G LNA probe exhibited high specificity for detection of transcripts from the Spike 614G gene but failed to amplify transcripts from the Spike D614 coding sequence (Table 5; *e.g*. probes in fourplex 5). Similarly, experiments by fourplex rRT-PCR assays verified that the LNA probe for the detection of Spike 501Y transcripts (nucleotide change A23,063T) (Tables 1 and 2) failed to detect the Spike N501 transcripts (Table 5, *e.g*. probes in fourplex 6). Consistently, results showed that the LNA probe for detection of the Spike N501 transcripts (nucleotide A at position 23,063) (Tables 1 and 2) did not detect the Spike 501Y transcripts (Table 5, *e.g*. probes in fourplex 7). In contrast, results showed that the Spike N501 and Spike 501Y LNA probes exhibited high specificity for detection of transcripts from the Spike N501 and Spike 501Y coding sequences, respectively (Table 5, *e.g*. probes in fourplex 8). Taken together, the results showed that the oligonucleotide probes containing LNA residues can be used for differential identification of SARS-CoV-2 viral templates in which mutations have occurred.

## Discussion

The ongoing COVID-19 pandemic characterized by successive waves of infection is a driving force for optimization and constant development of SARS-CoV-2 testing assays (Esbin et al. 2020; Graham et al. 2021). Here, we report the development of a fourplex TaqMan-LNA rRT-PCR assay that included per single reaction: four sets of primer pairs and four distinct TaqMan probes labelled with different fluorophores. Results showed that a tetraplex reaction in a single master mixture reached a detection limit of 14 viral transcript copies per reaction (Table 3). The amplification efficiencies of tetraplex reactions (values of 85.1% to 110.4%) were equivalent to those in the singleplex reactions (values of 83.3% and 95.1%) in which a single viral primer-probe set was used within one reaction master mixture (Table 2). The primer-probe set targeting a single nucleotide substitution in the Orf1ab L transcript was more sensitive than primer-probe sets targeting single nucleotide substitutions in the Orf1ab S, Spike 614G, Spike D614, Spike N501, and Spike 501Y transcripts. This set could detect transcript at levels as low as 1.4 copies per reaction in singleplex and fourplex rRT-PCR assays (Tables 2 and 3). In the cases of primer-probe sets targeting single nucleotide substitutions in the Orf1ab S, Spike 614G, Spike D614, Spike N501, and Spike 501Y transcripts, a detection limit of 14 transcript copies was reached when they were used in fourplex assays on three independent experiments (Table 3). In the case of the primer-probe set targeting the E gene that is denoted E Sarbeco (Corman et al. 2020), it was more sensitive than the N LSPQ probe targeting the N gene (Tables 2 and 3). The E Sarbeco primer-probe set could detect transcripts levels as low as 1.3 to 1.4 copies per reaction, whereas 14 transcript copies per reaction was reached with the N LSPQ primer-probe set under our experimental conditions.

During the course of studies, we have found that the human β2-microglobulin (B2M) gene encoding an ubiquitously expressed protein could be a positive control of choice for monitoring efficiency of RNA extraction from swab samples and subsequent reverse transcription reaction reliability. To create the B2M primer-probe set and ensure that its resulting product originates from RNA/complementary DNA (cDNA) amplification and not from genomic DNA amplification, we designed the reverse primer on exon2-exon3 junction of the B2M gene. Because this exon2-exon3 boundary sequence is only present in spliced B2M mRNA and not in genomic DNA, this approach eliminated the possibility that the B2M rRT-PCR signal resulted from an amplification of genomic DNA. Results showed that the B2M primer-probe set was highly sensitive when tested in singleplex assays, exhibiting an efficiency of 95.1% (r2 = 0.996) and a detection limit of 1.3 transcript copies per reaction (Fig. 1 and Table 2). When the B2M primer-probe set was used in fourplex assays, it remained highly sensitive with high levels of rRT-PCR efficiency (95.5% to 106.5%) and a similar detection limit (1.4 transcript copies per reaction) as compared to that observed in the case of singleplex assays (Figs 1 and 2) (Tables 2 and 3). The vast majority of COVID-19 tests are currently targeting the human RNase P gene (denoted RPP30 or h-RP) as a positive control to monitor RNA extraction and its quality to be reverse transcribed for its subsequent PCR amplification (Lu et al. 2020b; Lu et al. 2014). However, problems with the use of one CDC-recommended h-RP primer-probe set have been reported previously (Dekker et al. 2020). One issue is that the reverse primer hybridizes in the same exon as the forward primer and the TaqMan probe, therefore allowing amplification from genomic DNA instead of exclusively from spliced mRNA (Dekker et al. 2020). In that case, a positive result with this specific primer and probe set fails to validate the quality of the clinical sample specimen with respect to its RNA content.

The use of multiplex rRT-PCR assays allowed simultaneous targeting of different regions of the SARS-CoV-2 genome in a single master mixture. Among the targeted regions, we often included the E Sarbeco probe due to its reported successful use and high sensitivity (Corman et al. 2020; Nalla et al. 2020). In this manner, Ct values of the E Sarbeco probe in fourplex assays served as a point of reference to assess the significance of Ct values of the other tested probes. Interestingly, we have observed that utilization of multiplex SARS-CoV-2 primer/probe sets with distinct fluorophores shed light on which regions of the virus genome are more sensitive for a rRT-PCR-based SARS-CoV-2-specific identification assay. Another advantage of using a multitarget rRT-PCR method is the fact that it reduces the amount of non-enzymatic reagents, enzyme mixtures and plasticware required to analyze all the samples on a daily basis.

RNA viruses exhibit high mutation rates that may result in self beneficial properties by increasing their ability to infect host cells (Cuevas et al. 2015; Geller et al. 2016; Lauring et al. 2013). In the case of SARS-CoV-2, several mutations have been identified since the start of the pandemic (Fontanet et al. 2021; Mascola et al. 2021; Tegally et al. 2021b). For instance, two linked SNPs in the genome of SARS-CoV-2 have defined two major lineage of the virus that are designated S and L strains (Tang et al. 2020). Another example consists of nucleotide mutations found in the viral sequence of the gene encoding Spike that result in increased virus infectivity, especially nucleotide changes resulting in amino acid changes at Asp614 (D614G) and Asn501 (N501Y) (Korber et al. 2020; Leung et al. 2021). Here, we have used forward and reverse primers that hybridize to each side of hot spot mutational sites found in Orf1ab (nucleotide change T8782C, generating lineage S versus L), and at different locations in Spike such as nucleotide change A23403G (resulting in Spike variation D614G) and nucleotide change A23063T (giving rise Spike variation N501Y). These pairs of primers were used in combination with wild-type and mutant-LNA probes that hybridize at the mutation site. The fourplex real-time RT-PCR assay using LNA-based TaqMan probes exploits the 5’-3’ nuclease activity of Taq DNA polymerase, allowing direct detection of each of the four fluorogenic qPCR products by the uncoupling of a reporter dye from its quencher dye during qPCR. Furthermore, LNA modification of TaqMan probes increases their base pairing stability and creates a highly favorable context of duplex formation between them and their target sequences under more stringent conditions (Braasch and Corey 2001; McTigue et al. 2004; Petersen and Wengel 2003; Ugozzoli et al. 2004). Therefore, these properties allow LNA-based probes to discriminate between distinct viral genomic sequences that differ by a single nucleotide at a precise hot spot mutational site. This was performed, in three different places on the genome of SARS-CoV-2 at the same time into one reaction. In this way, it is possible to detect a range of representative genotypes that could reveal the presence of different variants of SARS-CoV-2 geographically. This method for differential diagnosis of SARS-CoV-2 virus isolates could be of help for effective surveillance. Furthermore, this multitarget rRT-PCR assay enables rapid follow-up of the most clinically relevant variants across the population in real-time and effectively guide the selection of viral strains that disserve further analysis by their whole-genome sequencing using next-generation sequencing methods.

## Data Availability

All data are included in the present manuscript. Oligonucleotides, gBlocks and plasmids used for this study are available if requested. The authors state that all results obtained for confirming the conclusions presented in the article are represented fully within the article.

## Conflict of interest

The authors declare that they have no conflict of interest with the content of this article.

## Author contributions

M.D. conceptualized, performed and analyzed rRT-PCR data. P.T. conceptualized the design of primers and TaqMan-LNA probes as well as carrying numerous bioinformatic analyses. S. Lévesque conceptualized and performed BioFire FilmArray assays. A.B. acquired and analyzed the data. A.C, L.V., and P.M. conceptualized the experimental work and analyzed data. S. Labbé conceptualized the experimental work and wrote the manuscript. The authors reviewed the results and approved the final version of the manuscript.

## Footnotes

1 Abbreviations used are: B2M, β2-microglobulin; Ct, threshold cycle; cDNA, complementary DNA; LNA, locked nuclei acid; PCR, polymerase chain reaction; NATRVP, NATtrol Respiratory verification panel; RT, reverse transcription; rRT-PCR, real-time quantitative RT-PCR; SNP, single nucleotide polymorphism.

## Supplemental figure legend

![Fig. S1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/03/26/2021.03.23.21254196/F3.medium.gif)

[Fig. S1.](http://medrxiv.org/content/early/2021/03/26/2021.03.23.21254196/F3)

Fig. S1. Sequences of the synthetic oligomers used to create the indicated gBlocks.
Each targeted SARS-CoV-2 region from the Orf1ab. Spike. E. and N genes (according to the Wuhan-Hu-1 genome reference sequence and mutant derivatives) were designed accompanied with a T7 promoter sequence (bold) at their 5’ends.

## Acknowledgments

We are grateful to Dr. Gilles Dupuis for critical reading of the manuscript and for his valuable comments. We thank Dr. Hugues Charest for his communication with respect to the N LSPQ primer-probe set and its utilization in rRT-PCR assays. This study was supported by the Ministère de l’Économie et de l’Innovation (MÉI) and, in part, by the Natural Sciences and Engineering Research Council of Canada (NSERC, grant #ALLRP 554545-20) to S.Labbé.

*   Received March 23, 2021.
*   Revision received March 23, 2021.
*   Accepted March 26, 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.  Braasch, D.A. and Corey, D.R. (2001) Locked nucleic acid (LNA): fine-tuning the recognition of DNA and RNA. Chemistry & Biology 8:1–7. [https://doi:10.1016/s1074-5521(00)00058-2](https://doi.org/10.1016/s1074-5521(00)00058-2)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1074-5521(00)00058-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11182314&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167275800002&link_type=ISI) 

2.  Chan, J.F., Kok, K.H., Zhu, Z., Chu, H., To, K.K., Yuan, S. and Yuen, K.Y. (2020) Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerging Microbes & Infections 9:221–236. [https://doi:10.1080/22221751.2020.1719902](https://doi.org/10.1080/22221751.2020.1719902)
    
    

3.  Corman, V.M., Landt, O., Kaiser, M., Molenkamp, R., Meijer, A., Chu, D.K.W., et al. (2020) Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro Surveillance 25:2000045. [https://doi:10.2807/1560-7917.es.2020.25.3.2000045](https://doi.org/10.2807/1560-7917.es.2020.25.3.2000045)
    
    

4.  Cuevas, J.M., Geller, R., Garijo, R., López-Aldeguer, J. and Sanjuán, R. (2015) Extremely high mutation rate of HIV-1 in vivo. PLoS Biology 13:e1002251. [https://doi:10.1371/journal.pbio.1002251](https://doi.org/10.1371/journal.pbio.1002251)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.1002251&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26375597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

5.  Dekker, R.J., Ensink, W.A., van Leeuwen, S., Rauwerda, H. and Breit, T.M. (2020) Overhauling a faulty control in the CDC-recommended SARS-CoV-2 RT-PCR test panel. bioRxiv:2020.06.12.147819. [https://doi:10.1101/2020.06.12.147819](https://doi.org/10.1101/2020.06.12.147819)
    
    

6.  Esbin, M.N., Whitney, O.N., Chong, S., Maurer, A., Darzacq, X. and Tjian, R. (2020) Overcoming the bottleneck to widespread testing: a rapid review of nucleic acid testing approaches for COVID-19 detection. RNA 26:771–783. [https://doi:10.1261/rna.076232.120](https://doi.org/10.1261/rna.076232.120)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoicm5hIjtzOjU6InJlc2lkIjtzOjg6IjI2LzcvNzcxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMjYvMjAyMS4wMy4yMy4yMTI1NDE5Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

7.  Fontanet, A., Autran, B., Lina, B., Kieny, M.P., Karim, S.S.A. and Sridhar, D. (2021) SARS-CoV-2 variants and ending the COVID-19 pandemic. Lancet 397:952–954. [https://doi:10.1016/s0140-6736(21)00370-6](https://doi.org/10.1016/s0140-6736(21)00370-6)
    
    

8.  Fronhoffs, S., Totzke, G., Stier, S., Wernert, N., Rothe, M., Bruning, T., et al. (2002) A method for the rapid construction of cRNA standard curves in quantitative real-time reverse transcription polymerase chain reaction. Molecular and Cellular Probes 16:99–110. [https://doi:10.1006/mcpr.2002.0405](https://doi.org/10.1006/mcpr.2002.0405)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/mcpr.2002.0405&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12030760&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176556200003&link_type=ISI) 

9.  Geller, R., Estada, Ú., Peris, J.B., Andreu, I., Bou, J.V., Garijo, R, et al. (2016) Highly heterogeneous mutation rates in the hepatitis C virus genome. Nature Microbiology 1:16045. [https://doi:10.1038/nmicrobiol.2016.45](https://doi.org/10.1038/nmicrobiol.2016.45)
    
    

10. Graham, T.G.W., Dugast-Darzacq, C., Dailey, G.M., Nguyenla, X.H., Van Dis, E., Esbin, M.N., et al. (2021) Open-source RNA extraction and RT-qPCR methods for SARS-CoV-2 detection. PloS One 16:e0246647. [https://doi:10.1371/journal.pone.0246647](https://doi.org/10.1371/journal.pone.0246647)
    
    

11. Guarneri, C., Rullo, E.V., Pavone, P., Berretta, M., Ceccarelli, M., Natale, A. and Nunnari, G. (2021) Silent COVID-19: what your skin can reveal. The Lancet Infectious Diseases 21:24–25. [https://doi:10.1016/s1473-3099(20)30402-3](https://doi.org/10.1016/s1473-3099(20)30402-3)
    
    

12. Holshue, M.L., DeBolt, C., Lindquist, S., Lofy, K.H., Wiesman, J., Bruce, H., et al. (2020) First case of 2019 novel coronavirus in the United States. New England Journal Medicine 382:929–936. [https://doi:10.1056/NEJMoa2001191](https://doi.org/10.1056/NEJMoa2001191)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001191&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

13. Kemp, S.A., Meng, B., Ferriera, I.A., Datir, R., Harvey, W.T., Papa, G., et al. (2021) Recurrent emergence and transmission of a SARS-CoV-2 spike deletion H69/V70. bioRxiv:2020.12.14.422555 doi:10.1101/2020.12.14.422555
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4xMi4xNC40MjI1NTV2NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzI2LzIwMjEuMDMuMjMuMjEyNTQxOTYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. Korber, B., Fischer, W.M., Gnanakaran, S., Yoon, H., Theiler, J., Abfalterer, W., et al. (2020) Tracking changes in SARS-CoV-2 Spike: evidence that D614G increases infectivity of the COVID-19 virus. Cell 182:812-827.e19. [https://doi:10.1016/j.cell.2020.06.043](https://doi.org/10.1016/j.cell.2020.06.043)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.06.043&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

15. Lai, C.C., Shih, T.P., Ko, W.C., Tang, H.J. and Hsueh, P.R. (2020) Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges. International Journal of Antimicrobial Agents 55:105924. [https://doi:10.1016/j.ijantimicag.2020.105924](https://doi.org/10.1016/j.ijantimicag.2020.105924)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijantimicag.2020.105924&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

16. Lauring, A.S., Frydman, J. and Andino, R. (2013) The role of mutational robustness in RNA virus evolution. Nature Reviews Microbiology 11:327–336. [https://doi:10.1038/nrmicro3003](https://doi.org/10.1038/nrmicro3003)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrmicro3003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23524517&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

17. Leung, K., Shum, M.H., Leung, G.M., Lam, T.T. and Wu, J.T. (2021) Early transmissibility assessment of the N501Y mutant strains of SARS-CoV-2 in the United Kingdom, October to November 2020. Euro Surveillance 26:2002106. [https://doi:10.2807/1560-7917.es.2020.26.1.2002106](https://doi.org/10.2807/1560-7917.es.2020.26.1.2002106)
    
    

18. Longtin, Y., Charest, H., Quach, C., Savard, P., Baz, M., Boivin, G., Fafard, J., Villeneuve, J., Roger, M., and De Serres, G. (2021) Infectivity of healthcare workers diagnosed with coronavirus disease 2019 (COVID-19) approximately 2 weeks after onset of symptoms: a cross-sectional study. Infection Control & Hospital Epidemiology. Jan 11:1-3. doi:10.1017/ice.2020.1420
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/ice.2020.1420&link_type=DOI) 

19. Lu, R., Zhao, X., Li, J., Niu, P., Yang, B., Wu, H., et al. (2020a) Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395:565–574. [https://doi:10.1016/s0140-6736(20)30251-8](https://doi.org/10.1016/s0140-6736(20)30251-8)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)302518&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32007145&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

20. Lu, X., Wang, L., Sakthivel, S.K., Whitaker, B., Murray, J., Kamili, S., et al. (2020b) US CDC real-time reverse transcription PCR panel for detection of severe acute respiratory syndrome coronavirus Emerging Infectious Diseases 26:1654–1665. [https://doi:10.3201/eid2608.201246](https://doi.org/10.3201/eid2608.201246)
    
    

21. Lu, X., Whitaker, B., Sakthivel, S.K., Kamili, S., Rose, L.E., Lowe, L., et al. (2014) Real-time reverse transcription-PCR assay panel for Middle East respiratory syndrome coronavirus. Journal of Clinical Microbiology 52:67–75. [https://doi:10.1128/jcm.02533-13](https://doi.org/10.1128/jcm.02533-13)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjc6IjUyLzEvNjciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMy8yNi8yMDIxLjAzLjIzLjIxMjU0MTk2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

22. Martin, D.P., Weaver, S., Tegally, H., San, E.J., Shank, S.D., Wilkinson, E., et al. (2021) The emergence and ongoing convergent evolution of the N501Y lineages coincides with a major global shift in the SARS-CoV-2 selective landscape. medRxiv. [https://doi:10.1101/2021.02.23.21252268](https://doi.org/10.1101/2021.02.23.21252268)
    
    

23. Mascola, J.R., Graham, B.S. and Fauci, A.S. (2021) SARS-CoV-2 viral variants-tackling a moving target. JAMA. Published On-line. [https://doi:10.1001/jama.2021.2088](https://doi.org/10.1001/jama.2021.2088)
    
    

24. McTigue, P.M., Peterson, R.J. and Kahn, J.D. (2004) Sequence-dependent thermodynamic parameters for locked nucleic acid (LNA)-DNA duplex formation. Biochemistry 43:5388–5405. https:// doi:10.1021/bi035976d
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/bi035976d&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15122905&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

25. Nalla, A.K., Casto, A.M., Huang, M.-L.W., Perchetti, G.A., Sampoleo, R., Shrestha, L., et al. (2020) Comparative performance of SARS-CoV-2 detection assays using seven different primer-probe sets and one assay kit. Journal of Clinical Microbiology 58:e00557–20. [https://doi:10.1128/jcm.00557-20](https://doi.org/10.1128/jcm.00557-20)
    
    

26. Oude Munnink, B.B., Nieuwenhuijse, D.F., Stein, M., O’Toole, A., Haverkate, M., Mollers, M., Karnga, S.K., et al. (2020) Rapid SARS-CoV-2 whole-genome sequencing and analysis for informed public health decision-making in the Netherlands. Nature Medicine 26:1405–1410. [https://doi:10.1038/s41591-020-0997-y](https://doi.org/10.1038/s41591-020-0997-y)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0997-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32678356&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

27. Petersen, M. and Wengel, J. (2003) LNA: a versatile tool for therapeutics and genomics. Trends in Biotechnology 21:74–81. [https://doi:10.1016/s0167-7799(02)00038-0](https://doi.org/10.1016/s0167-7799(02)00038-0)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0167-7799(02)00038-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12573856&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000181154900007&link_type=ISI) 

28. Poritz, M.A., Blaschke, A.J., Byington, C.L., Meyers, L., Nilsson, K., Jones, D.E., et al. (2011) FilmArray, an automated nested multiplex PCR system for multi-pathogen detection: development and application to respiratory tract infection. PloS One 6:e26047. [https://doi:10.1371/journal.pone.0026047](https://doi.org/10.1371/journal.pone.0026047)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0026047&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22039434&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

29. Rivett, L., Sridhar, S., Sparkes, D., Routledge, M., Jones, N.K., Forrest, S., et al. (2020) Screening of healthcare workers for SARS-CoV-2 highlights the role of asymptomatic carriage in COVID-19 transmission. Elife 9:e58728. [https://doi:10.7554/eLife.58728](https://doi.org/10.7554/eLife.58728)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.58728&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32392129&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

30. Rondinone, V., Pace, L., Fasanella, A., Manzulli, V., Parisi, A., Capobianchi, M.R., et al. (2021) VOC 202012/01 variant is effectively neutralized by antibodies produced by patients infected before its diffusion in Italy. Viruses 13:276. [https://doi:10.3390/v13020276](https://doi.org/10.3390/v13020276)
    
    

31. Starr, T.N., Greaney, A.J., Hilton, S.K., Ellis, D., Crawford, K.H.D., Dingens, A.S., et al. (2020) Deep mutational scanning of SARS-CoV-2 receptor binding domain reveals constraints on folding and ACE2 binding. Cell 182:1295-1310.e20. [https://doi:10.1016/j.cell.2020.08.012](https://doi.org/10.1016/j.cell.2020.08.012)
    
    

32. Tang, X., Wu, C., Li, X., Song, Y., Yao, X., Wu, X., Duan, Y., Zhang, H., Wang, Y., Qian, Z., Cui, J. and Lu, J. (2020) On the origin and continuing evolution of SARS-CoV-2. National Science Review 7:1012–1023. [https://doi:doi.org/10.1093/nsr/nwaa036](https://doi.org/doi.org/10.1093/nsr/nwaa036)
    
    

33. Tegally, H., Wilkinson, E., Giovanetti, M., Iranzadeh, A., Fonseca, V., Giandhari, J., et al. (2021a) Emergence of a SARS-CoV-2 variant of concern with mutations in spike glycoprotein. Nature. March 9, On-line. [https://doi:10.1038/s41586-021-03402-9](https://doi.org/10.1038/s41586-021-03402-9)
    
    

34. Tegally, H., Wilkinson, E., Lessells, R.J., Giandhari, J., Pillay, S., Msomi, N., et al. (2021b) Sixteen novel lineages of SARS-CoV-2 in South Africa. Nature Medicine 27:440–446. [https://doi:10.1038/s41591-021-01255-3](https://doi.org/10.1038/s41591-021-01255-3)
    
    

35. Ugozzoli, L.A., Latorra, D., Puckett R., Arar, K. and Hamby, K. (2004) Real-time genotyping with oligonucleotide probes containing locked nucleic acids. Analytical Biochemistry 324:143–152. [https://doi:10.1016/j.ab.2003.09.003](https://doi.org/10.1016/j.ab.2003.09.003)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ab.2003.09.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14654057&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187566500020&link_type=ISI) 

36. Vijgen, L., Keyaerts, E., Moës, E., Maes, P., Duson, G. and Van Ranst, M. (2005) Development of one-step, real-time, quantitative reverse transcriptase PCR assays for absolute quantitation of human coronaviruses OC43 and 229E. Journal of Clinical Microbiology 43:5452–5456. [https://doi:10.1128/jcm.43.11.5452-5456.2005](https://doi.org/10.1128/jcm.43.11.5452-5456.2005)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjEwOiI0My8xMS81NDUyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMjYvMjAyMS4wMy4yMy4yMTI1NDE5Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

37. Wu, F., Zhao, S., Yu, B., Chen, Y.M., Wang, W., Song, Z.G., et al. (2020) A new coronavirus associated with human respiratory disease in China. Nature 579:265–269. [https://doi:10.1038/s41586-020-2008-3](https://doi.org/10.1038/s41586-020-2008-3)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2008-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom) 

38. Yang, R., Gui, X. and Xiong, Y. (2020) Comparison of clinical characteristics of patients with asymptomatic vs symptomatic coronavirus disease 2019 in Wuhan, China. JAMA Network Open 3:e2010182. [https://doi:10.1001/jamanetworkopen.2020.10182](https://doi.org/10.1001/jamanetworkopen.2020.10182)
    
    

39. Yi, C., Sun, X., Ye, J., Ding, L., Liu, M., Zhang, Y., et al. (2020) Key residues of the receptor binding motif in the spike protein of SARS-CoV-2 that interact with ACE2 and neutralizing antibodies. Cellular & Molecular Immunology 17:621–630. [https://doi:10.1038/s41423-020-0458-z](https://doi.org/10.1038/s41423-020-0458-z)
    
    

40. Zhang, L., Jackson, C.B., Mou, H., Ojha, A., Peng, H.,Quinlan, B.D. et al. (2020) SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity. Nature Communications 11:6013. [https://doi:10.1038/s41467-020-19808-4](https://doi.org/10.1038/s41467-020-19808-4)
    
    

41. Zhang, Y.Z. and Holmes, E.C. (2020) A genomic perspective on the origin and emergence of SARS-CoV-2. Cell 181:223–227. [https://doi:10.1016/j.cell.2020.03.035](https://doi.org/10.1016/j.cell.2020.03.035)
    
    

42. Zhu, N., Zhang, D., Wang, W., Li, X., Yang, B., Song, J., et al. (2020) A novel coronavirus from patients with pneumonia in China, 2019. New England Journal Medicine 382:727–733. [https://doi:10.1056/NEJMoa2001017](https://doi.org/10.1056/NEJMoa2001017)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.23.21254196.atom)