- 1Biomeme, Inc., Philadelphia, PA, United States
- 2MRIGlobal, Kansas City, MO, United States
- 3Naval Health Research Center (NHRC), San Diego, CA, United States
- 4Integrated DNA Technologies, Coralville, IA, United States
- 5Defense Biological Product Assurance Office (DBPAO), Joint Program Executive Office for Chemical, Biological, Radiological and Nuclear Defense (JPEO-CBRND), Enabling Biotechnologies, Frederick, MD, United States
- 6Logistics Management Institute, Tysons, VA, United States
Many SARS-CoV-2 variants have emerged during the course of the COVID-19 pandemic. These variants have acquired mutations conferring phenotypes such as increased transmissibility or virulence, or causing diagnostic, therapeutic, or immune escape. Detection of Alpha and the majority of Omicron sublineages by PCR relied on the so-called S gene target failure due to the deletion of six nucleotides coding for amino acids 69–70 in the spike (S) protein. Detection of hallmark mutations in other variants present in samples relied on whole genome sequencing. However, whole genome sequencing as a diagnostic tool is still in its infancy due to geographic inequities in sequencing capabilities, higher cost compared to other molecular assays, longer turnaround time from sample to result, and technical challenges associated with producing complete genome sequences from samples that have low viral load and/or high background. Hence, there is a need for rapid genotyping assays. In order to rapidly generate information on the presence of a variant in a given sample, we have created a panel of four triplex RT-qPCR assays targeting 12 mutations to detect and differentiate all five variants of concern: Alpha, Beta, Gamma, Delta, and Omicron. We also developed an expanded pentaplex assay that can reliably distinguish among the major sublineages (BA.1–BA.5) of Omicron. In silico, analytical and clinical testing of the variant panel indicate that the assays exhibit high sensitivity and specificity. This panel can help fulfill the need for rapid identification of variants in samples, leading to quick decision making with respect to public health measures, as well as treatment options for individuals. Compared to sequencing, these genotyping PCR assays allow much faster turn-around time from sample to results—just a couple hours instead of days or weeks.
Background
The emergence of COVID-19 in late 2019 in Wuhan, China and subsequent spread of the virus across the globe led to the unprecedented COVID-19 pandemic and declaration as a Public Health Emergency of Global Concern on Jan 30, 2020 by the World Health Organization (WHO) (1). Identification of the etiological agent of COVID-19 as SARS-CoV-2 and release of the first whole genome sequence of SARS-CoV-2 shortly thereafter on Jan 11, 2020 was instrumental in the rapid development of diagnostics, therapeutics and vaccines (2, 3). Public health measures (such as “test, trace and isolate”) and individual health measures (such as treating symptomatic individuals as needed) rely on rapid and accurate diagnosis early in the course of infection. The first molecular diagnostic test was issued an emergency use authorization (EUA) on Feb 04, 2020. As of May 05, 2022 (~2 years later), there were 307 EUA SARS-CoV-2 molecular tests and 48 antigen tests (4). The COVID-19 pandemic response has been complicated by the fact that as many as 40% of infected individuals may be asymptomatic depending on the population cohort, while secondary transmission might still be occurring (5). Hence, PCR tests remain the gold standard for SARS-CoV-2 detection because of their exquisite sensitivity and better performance identifying asymptomatic cases when compared to rapid antigen tests (6). It has been argued that PCR tests may show positivity for a longer duration—even after the individual is no longer infectious—and hence rapid antigen tests may be more reflective of infectivity (7, 8). However, PCR tests are preferred for their sensitivity, especially in testing for asymptomatic infections and early in the infection cycle (9, 10).
Despite the availability of hundreds of different tests, testing demands and needs have not been fully met and there is always a shortage of test quantities and sites around the globe. This adversely affects rapid decision making to test, trace, and isolate infected persons to prevent secondary transmissions (11). For this reason, there is a need for rapid and accurate tests that have higher sensitivity and can be performed in less complex settings. Indeed, at-home diagnostic kits meet the turn-around time requirements, but few options are currently available for sensitive PCR-based tests compared to the less sensitive rapid antigen tests (12). In addition to testing constraints, the emergence of variants has further complicated testing capabilities.
Rapid proliferation and sustained transmission of SARS-CoV-2 in the human population resulted in waves of variants that swept through the globe (13, 14). To identify the variants, genomic surveillance has come to the forefront of the COVID-19 pandemic response, as evidenced by the tremendous explosion of whole genome sequencing data from all around the globe (15). Whole genome sequencing has become ubiquitous, and as of May 24, 2022, more than 11 million SARS-CoV-2 sequences have been generated and deposited to the GISAID database by the global sequencing community (16–18). Whole genome sequence data have been extensively used in determining the changing genetic profiles of viral variants and tracking viral mobility across communities in what has come to be known as genomic epidemiology (19–21). Rapid identification of the prevalent variants without having to wait weeks for sequence data is paramount for public health measures and decision-making, and in some instances even treatment options. For example, some antibody therapies fail against Omicron infections due to specific mutations in critical positions in the spike protein (22), and antibody evasion capabilities vary not only amongst the major lineages (23), but also amongst the Omicron sublineages (24). This prompted the U.S. Food and Drug Administration to release a statement in January of 2022 that limited the use of certain antibody treatments to only patients known to be infected/exposed to susceptible lineages (25). Hence, there is a need for RT-qPCR assays that can rapidly provide genotype information just a couple hours after sample collection, compared to the many days that are usually required before sequencing results are available (26–36).
Here we describe a panel of multiplex RT-qPCR assays to detect and differentiate all five WHO designated variants of concern: Alpha, Beta, Gamma, Delta, and Omicron. These RT-qPCR assays can be performed on a hand-held, battery-operated device, the Biomeme Franklin™ three9 Real-Time PCR Thermocycler. We chose the Biomeme platform because it is portable, easy to use, battery-powered, and hence, field deployable with a minimal footprint. Biomeme also provides multiple assay formats—including lyophilized, room-temperature stable 3-well Go-Strips and 96-well Go-Plates—that can be used in both point-of-need, low-throughput devices like the FranklinTM, as well as in traditional lab-based, high-throughput instruments.
Materials and methods
Primers, probes and PCR reagents
The primer, probe, and template sequences that were designed in this study are presented in Supplementary Tables 1, 2. Reagents and assay conditions are described under different subsections below.
Assay design for Triplexes 1–3
The whole genome sequences of SARS-CoV-2 variants were downloaded from the GISAID database and aligned to the genome of the Wuhan-Hu-1 isolate (GenBank accession no. MN908947.3) using Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/). Primers were designed at conserved regions of the spike (S) and nucleocapsid (N) genes to generate amplicons comprising the mutations of interest using Primer Quest Software (https://www.idtdna.com/PrimerQuest/Home/Index) from Integrated DNA Technologies, Inc. (IDT; Coralville, Iowa, USA) and analyzed for potential dimerization and hairpin products by the OligoAnalyzer tool (IDT). Probes containing LNA and conjugated minor groove binder (MGB) were manually designed. In particular, the assistance of IDT's in-house tools afforded the predication of melting temperatures and mismatch discrimination potentials. The S:Y144del primer and probe sequences were designed according to Vogels et al. (37). Primer candidates were then analyzed by Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) and the probes by blastn (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE_TYPE=BlastSearch) to ensure recognition of the intended targets. For suitability of multiplexing, primer combinations were predicted by MultiPLX 2.1 (38). Primers and LNA probes were purchased either from IDT or Sigma-Aldrich (St. Louis, MO, USA), with the primers purified by standard desalting and the probes by high-performance liquid chromatography (HPLC). The MGB probe was purchased from Applied Biosystems (Foster City, CA, USA) as HPLC-purified products. gBlocks were purchased from IDT. An overview of the assay design and development workflow is presented in Figure 1.
Figure 1. Biomeme assay design, testing and validation workflow. Triplex assay development workflow. Using the Wuhan-Hu-1 isolate (accession no. MN908947.3) as the reference, whole genomes of SARS-CoV-2 variants were compared. The in silico framework consisted of identification of appropriate primers and probes and their subsequent suitability in triplex combinations. The primers were designed at regions of consensus sequences, and the probes were designed manually at target mutation sites. For the locked nucleic acid (LNA) probes, melting temperature predictions served as guides for mismatch discrimination. Candidate primers and probes were then aligned using BLAST to ensure specificity, and suitable triplex combinations were subsequently predicted. Identified triplexes were then tested and optimized initially using dsDNA templates (IDT gBlocks), followed by the use of in vitro-transcribed RNA templates. The optimized conditions were formulated into freeze-dried Go-Strips and Go-Plates for use on the Biomeme FranklinTM and standard 96-well RT-qPCR instruments.
RT-qPCR on standard thermocyclers
One-step reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed to avoid potential cross-contamination. A typical 20-μL reaction containing the DNA or RNA sample was prepared in LyoRNA 2.0 Master Mix (Biomeme, Inc., Philadelphia, PA, USA) according to the manufacturer's instructions and assayed on a Bio-Rad CFX96 instrument (Bio-Rad Laboratories, Inc., Hercules, CA, USA) in a triplex mixture. The assays were run under the following conditions: reverse transcription at 55°C for 2 min, initial denaturation at 95°C for 1 min, and 45 cycles of 95°C for 3 s and annealing/extension at 62°C for 30 s. Assay development was accomplished using gBlock templates in the RT-qPCR reactions and later confirmed with in vitro-transcribed RNA templates generated by the MEGAscript T7 Transcription Kit (Invitrogen, Waltham, MA, USA) according to the manufacturer's instructions using 80 ng of gBlock templates. Triplexes were tested against both wild-type and mutant templates, and only the probes that could specifically discriminate between the two templates were selected for further testing and optimization.
RT-qPCR on the Biomeme FranklinTM Thermocycler
Incorporating the same primer and probe concentrations used for measurement on the Bio-Rad CFX96 instrument, triplex assays were assembled into freeze-dried test strips (Go-Strips and Go-Plates) using Biomeme's proprietary formulations. Each well of the Go-Strip or Go-Plate was reconstituted with 20 μL total volume of nuclease-free water and template DNA/RNA. The Go-Strips were then loaded into the portable FranklinTM Real-Time Thermocycler (39) in the appropriate orientation and monitored by the Biomeme software application on the included smartphone, while the Go-Plates were measured on the Bio-Rad CFX96 instrument. RT and cycling conditions used on all platforms were the same as those described for the LyoRNA master mix.
Lyophilized assay manufacturing
Triplex assays were lyophilized into Go-Strip and Go-Plate formats. Go-Strips contained three reactions packaged individually in a small foil pouch designed for use on the Biomeme FranklinTM. Go-Plates are 96-well plates packaged in a foil pouch meant for standard laboratory RT-qPCR systems. Go-Strips and Go-Plates are shelf-stable at room temperatures for 18 months in their original packaging.
Omicron assay design and optimization
Probes for the omicron multiplexes (Triplex 4 and the Omicron pentaplex) were manually designed to overlie the targeted mutation site, using Primer Express® Software Version 3.0.1 to estimate melting temperatures for both the standard and MGB probes. Candidate primer pairs were designed using Primer3 Plus (40), then checked against the NCBI nt database for off-target matches and on-target mismatches using the NCBI Primer-BLAST tool (41). Candidate primer pairs that passed NCBI Primer-BLAST screening and candidate probe sequences were then tested for multiplex compatibility (oligo interactions and melting temperature) using ThermoFisher's Multiple Primer Analyzer (42), and some manual adjustments to primer sequences were made to enhance compatibility. Candidate primer and probe combinations were then screened, down-selected, and optimized using 10-fold serial dilutions of template (50–50,000 copies per reaction) to evaluate sensitivity, efficiency, and signal-to-noise ratios. Candidate assays were also tested against a high-level of wild-type genomic RNA to confirm exclusivity against non-variant template. Initial testing and optimization were performed in singleplex reactions using gBlock templates, but then progressed to multiplex testing using synthetic SARS-CoV-2 gRNA commercially available from Twist Biosciences (San Francisco, CA, USA), as well as RNA extracted from patient specimens and stock cultures (once available). Because gRNA could not be obtained for BA.1.1, BA.3, BA.4, and BA.5, gBlocks encompassing the S gene and relevant portion of the N gene were designed based on reference sequences, and serial dilutions of those were used for testing.
Verification of assay signatures in silico
To evaluate inclusivity of each assay on a per-lineage basis, SARS-CoV-2 genomes designated as “complete” (N = 1005,125) were downloaded from NCBI Virus (43) on March 31, 2022, along with a corresponding metadata table containing the Pangolin (44) lineage assignments. Upon observing that representation of the BA.3 lineage was limited to a single record in the dataset drawn from NCBI Virus, we supplemented the dataset by downloading BA.3 assigned records (N = 4,370) from GISAID.org on April 15, 2022. Similarly, BA.4 sequences (N = 1,174) and BA.5 sequences (N = 663) were downloaded from GISAID.org on May 11, 2022. A complete table of author attributions for the GISAID sequences is available in Supplementary Table 3.
Inclusivity was calculated as set coverage for each assay target for each lineage. A lineage's coverage was calculated by dividing the in silico predicted PCR-positive sequences (based on cutoff criteria for primer and probe alignment mismatches) by the total number of sequences assigned that lineage in the database. Formally, the inclusivity C of each primer set is: , where B is the set of records returned by a primer pair/probe and A is the set of records for a given lineage.
In silico prediction of amplicon products for the evaluated assays was conducted using a modified version of simulate_PCR (45) re-implemented in Python. This approach relies on BLAST results of the primers along with pre-specified cutoff criteria for amplicon product sizes and mismatch criteria to aggregate candidate PCR products from a target database. This pipeline was modified to better accommodate the use of short probe sequences that may contain modifications such as minor groove binders by decomposing the problem into a BLAST-based estimate of amplicon products returned by the primers, followed by a Smith-Waterman local alignment of the probe sequences against the predicted amplicon products. Because of the computational cost of alignment with expanding degeneracies, predicted amplicons that contain “N” degenerate characters are not aligned and not considered as inclusive hits. Probe alignments that meet the designated mismatch threshold are annotated to the corresponding amplicon product. This amplicon table is then evaluated using a custom Python script for inclusivity by the criteria above, where inclusivity for a set of oligos is defined as a forward primer, reverse primer, and probe sequence are all aligned within the pre-specified mismatch cutoffs. The output from this script is a summary table that gives inclusivity estimates for each unique label, in this case, the Pangolin lineage. Cutoff criteria set inclusion for all assays was a maximum of 3 mismatches in either the forward or reverse primer, 3 mismatches for standard (non-LNA) probes, and a 0 mismatch tolerance for LNA probes.
MRIGlobal analytical testing
Assay inclusivity/exclusivity testing for all triplexes was performed using serial dilutions of RNA extracted from live SARS-CoV-2 viral stocks, whereas Omicron pentaplex testing was performed using synthetic gRNA from Twist Biosciences due to lack of available stocks for Omicron sublineages. Viral stocks were extracted using the Qiagen Viral RNA Mini Kit (Qiagen, MD, USA). A 10-fold dilution series of each extract was made and then screened for quality/quantity using the Biomeme SARS-CoV-2 Dx assay on a Bio-Rad CFX96 PCR instrument. Dilutions producing Cq values in the early 20s to late 30s (or negative results) were selected for further testing with the variant multiplexes, side-by-side with the Biomeme SARS-CoV-2 Dx assay as a control/comparator assay. This testing was performed using lyophilized versions of Triplexes 2 and 3, but “wet” versions of Triplex 1, Triplex 4, and the Omicron pentaplex were used because final lyophilized versions were not available. “Wet” versions were prepared using commercially available frozen master mix (TaqPath 1-Step RT-qPCR Master Mix, CG from ThermoFisher (St. Louis, MO, USA). When “wet” versions were used, the thermocycling protocol changed to the following: UNG incubation at 25°C for 2 min, reverse transcription at 50°C for 15 min, initial denaturation at 95°C for 2 min, and 45 cycles of 95°C for 3 s and annealing/extension at 62°C for 30 s. In addition to testing the full dilution series on the Bio-Rad CFX96, single-level testing of each extract was also performed on the handheld Biomeme FranklinTM PCR device. Viral copy numbers for gRNA dilutions were later determined by droplet digital PCR (ddPCR) on a Bio-Rad QX200 (Hercules, CA, USA) and CDC 2019-nCov N1 assay primers and probe.
MRIGlobal testing of RADx residual clinical specimens
Twenty SARS-CoV-2 residual clinical specimens were provided by the National Institutes of Health (NIH) Rapid Acceleration of Diagnostics (RADx®) program. These samples consisted of heat inactivated nasopharyngeal swabs in transport media that had previously tested positive for SARS-CoV-2 and had been sequenced to determine lineage before being provided to us for this study. Five samples each of four different lineages—Alpha, Beta, Gamma, and Delta—were provided by NIH RADx® in early December of 2021 (before Omicron samples were available). Specimens were extracted using the Qiagen Viral RNA Mini kit using 140 μL of sample as input and an elution volume of 70 μL. Extracts were tested with all four variant triplexes, as well as with the Biomeme SARS-CoV-2 Dx assay as a control/comparator assay. This testing was performed on a CFX96 instrument using 2 μL of extract per reaction, and testing duplicate PCR reactions of each sample per triplex. Sample identities were blinded to the operators until after data analysis and interpretation was performed.
Naval Health Research Center (NHRC) clinical sample testing
Clinical samples tested by NHRC consisted of nasal and nasopharyngeal (NP) swabs in either transport media or saline. For sensitivity and specificity testing, we used clinical specimens positive for SARS-CoV-2 of known lineage (determined by the providing lab via whole genome sequencing). Ten samples were selected for each variant targeted by the panel except for B.1.351, for which NHRC only had four samples. An additional 20 variants not targeted by the panel were also selected to test for specificity of the variant panels. Viral RNA extraction was performed on 100 μl samples using the Qiagen Viral RNA Mini Kit. Extracts were tested with the CDC 2019-nCoV N1 and N2 assays on an ABI 7500 real-time PCR system (Foster City, CA, USA) to confirm presence of SARS-CoV-2 RNA. Variant panel testing was also conducted on the ABI 7500 using TaqPath 1-Step RT-qPCR Master Mix and the variant panel primers and probes provided by MRIGlobal at a 10X concentration. After initial variant panel testing, high-throughput testing was performed for >1,500 samples using Triplex 4 (the Omicron triplex). During high-throughput testing, all samples were tested with Triplex 4, and in some cases, only samples that were negative for Triplex 4 assays were then tested with the CDC N1 and N2 assays to verify presence/absence of SARS-CoV-2 in the sample.
Results
Assay target selection
The Biomeme FranklinTM Real-Time Thermocycler includes three different fluorescence channels, enabling one triplex assay per well. Based on the high prevalence of unique markers among the variant sequences available at the time of assay design, we selected markers to include in each triplex. These markers were either a single nucleotide polymorphism (SNP), multiple SNPs in close proximity, or an insertion or deletion (indel) spanning the probe region. For Triplexes 1 and 2, markers were selected in a way that some are shared among variant lineages and others are unique to a specific lineage to allow for identifying three variants using six markers. For Triplexes 3 and 4, markers were selected for specific identification of a single variant per triplex (Delta and Omicron, respectively). The targeted markers and algorithm for adjudication of results and variant calling are presented in Table 1. Data were obtained from outbreak.info (46).
Prevalence and distribution of sequences carrying the markers are also shown in the last three columns of Table 1. The total number of sequences in the database for each variant of concern is displayed, as well as the total number of sequences with a given set of markers. The last column represents the percentage of sequences of the specific sublineage listed in column 2 that carries the set of expected mutations. All triplexes cover 81–96% of the specific sublineages they were designed to detect and distinguish. Since the time of assay design, many new sublineages have emerged, including Omicron sublineages BA.3, BA.4, and BA.5. With the increasing number of sublineages and diversity of markers combined with the convergent evolution of markers, the assignment to sublineages based on just two or three markers has become a challenge. In order to address this issue and to increase the confidence of calls for Omicron sublineages, we expanded the panel to a pentaplex assay described in a later section.
Assay design
The rationale for assay design was to amplify and detect the target amplicon only if the mutation was present in the template. This was accomplished by designing dual-labeled fluorescent probes that can discriminate between the variant markers and the wild-type sequence, in addition to optimizing the annealing temperature and concentration of the probes. For SNP markers in particular, specificity was accomplished using LNA probe designs. As a first step toward assay design, the primer/probe sequences were examined in silico to assess the sensitivity and specificity of the designs against whole genome sequences of SARS-CoV-2 variants obtained from GenBank and GISAID (Table 2).
The in silico results indicated very high specificity among the variants of concern, with ≤ 0.6% of sequences for a particular lineage containing any marker not typically associated with that lineage. Predicted sensitivity was >90% for most of the markers, although there were some exceptions for Delta and Omicron lineages. This reduction in sensitivity is generally a reflection of absence of the mutations (as opposed to PCR failure), as there can be differences among sequences of the same lineage and even among sublineages of the same variant. In some instances, there are shared mutations in some lineages or sublineages due to convergent evolution. For example, S:N501Y was initially designed as a shared mutation for the three variants targeted in Triplexes 1 and 2. It is not present in the Delta sequences, but the majority of Omicron sequences carry this mutation. However, this assay failed to produce detection during analytical testing with Omicron (see Table 3) due to mismatches in the forward and reverse primer regions (Supplementary Figure 1). These results were further confirmed by clinical sample testing, during which no positive results were obtained for the S:N501Y assay with Omicron samples (see Supplementary Tables 4, 5). It is noted that when primers were redesigned to incorporate mixed bases at the mismatch sites, the S:N501Y mutation in Omicron could then be reliably detected (unpublished data). Similarly, the S:69/70del mutation is present in both Alpha and Omicron BA.1; however, in the latter, there is an additional SNP in the probe binding region (Supplementary Figures 2, 3). The S:69/70del probe used in Triplex 2 can still detect Omicron albeit at reduced amplitude, while a more Omicron specific probe can be discriminatory between Alpha and Omicron templates (Supplementary Figures 4A–C).
Analytical testing of Triplexes 1–4
Initial assay testing and optimization was performed using dsDNA templates and/or in vitro transcribed RNA templates harboring the targeted regions (data not shown). Once the assays and reaction conditions were finalized, analytical and clinical testing were performed. Analytical testing of the triplex assays was performed using a serial dilution of live viral stocks after extraction as described in the Materials and Methods section. All assays exhibited the expected results for each lineage. Average Cq values (n = 2) are presented for dilutions that produced a Cq value of around 22–25 (Table 3). Comparable values were observed between the diagnostic assays (SARS-CoV-2 Dx: spike and Orf1ab) and the variant specific assays.
For each variant gRNA dilution series and the relevant assays, plots of Cq values vs. viral copies per reaction are presented in Figure 2. The PCR efficiencies are presented in Supplementary Table 6. Amplification efficiencies of all assays with the tested lineages were >80%, with two exceptions (P26S at 79%, and 156/157del at 76%).
Figure 2. Serial dilution of extracted gRNA of different SARS-CoV-2 variants were used as templates in RT-qPCR reactions. The Cq values are plotted as a function of the concentrations. PCR efficiencies were derived from the slopes of the curves.
Clinical sample testing
Retrospective testing of the four triplexes using clinical samples was performed as described in the Materials and Methods section. The results are shown in Table 4 and the detailed sample data with Cq values are provided in Supplementary Tables 4, 5. The 64 samples presented in Table 4 are from two different sources: NIH RADx® and NHRC. Overall, the panels have very high sensitivity and specificity (see Supplementary Figures 5A,B). However, in this analysis results of the S:P26S assay in Triplex 1 have been excluded due to indiscriminate signal (and hence false positives) with all SARS-CoV-2 samples analyzed on the ABI 7500 (see Supplementary Figure 5B). Therefore, positive calls for the P.1 lineage are based on only two markers instead of three.
Table 4. Positive percent agreement (PPA) and negative percent agreement (NPA) of Triplexes 1–4 for each of the targeted lineages.
In addition to the 64 samples from targeted lineages, 20 samples (1 per lineage) from lineages not targeted by the panel were also tested, with only one producing a false positive lineage call; a B.1.621 specimen was falsely identified as P.1 on the basis of having the same two mutations as P.1 when the P26S assay is ignored. This was the only false positive call obtained from the test set of 84 clinical specimens used to assess PPA and NPA of Triplexes 1–4. Difficulty distinguishing amplification from background for the S:P26S assay was seen both in the NHRC clinical sample testing on ABI 7500 systems (Supplementary Figure 5B) and in testing on the Biomeme FranklinTM (Figure 4), but was less of an issue during analytical testing and testing of the RADx clinical specimens on the Bio-Rad CFX96.
Additional clinical sample testing with Triplex 4
Due to the surge in US clinical samples and testing needs for Omicron, beginning in December of 2021, additional testing was performed on the Omicron triplex (Triplex 4). This consisted of dynamic range finding and relative limit of detection (LoD) testing with a clinical specimen, additional sensitivity and specificity testing with an expanded set of samples, and high-throughput testing of >1,500 clinical samples being screened for sequencing.
Dynamic range and relative limit of detection (LoD) of Triplex 4
In order to assess the dynamic range and limit of detection of the Triplex 4 panel of assays, a 10-fold serial dilution of a clinical specimen extract was made and PCR was performed as described in the Materials and Methods section. This specimen was known to be SARS-CoV-2 positive and of sublineage BA.1 based on prior PCR and sequencing. Triplicate reactions of Triplex 4 assays were tested for each dilution and the working LoD was determined as the lowest sample level producing 3 of 3 positive results. Working LoD was the same for all three assays in the triplex, with Cq values ranging from 34 to 37 at that level (Table 5). This working LoD was 1:1e5 dilution of the specimen extract, which produced N1 and N2 Cq values of ~18 at the undiluted level.
Sensitivity and specificity of Triplex 4
We further assessed sensitivity and specificity of the Triplex 4 using newly obtained clinical samples. Twenty of 20 Omicron (BA.1) clinical specimens tested positive for all three Omicron assays, including three specimens with “failed” sequencing results due to lower titers. The Cq values of CDC N1 and CDC N2 assays ranged from ~16–30 for these samples. Fourteen of 14 non-Omicron specimens (all of different lineages; see Supplementary Table 5 for lineages) were negative for all three Omicron assays and positive for the CDC N1 and N2 assays, with Cq values of 22–30. Based on these data for sequence-confirmed samples, the sensitivity and specificity were 100% for all three assays in the panel. Ten samples without sequence data were also tested with Triplex 4. Of these ten unknown samples, eight were positive for all three assays, and two were negative. Of the two negative samples, one was positive for the CDC N1 and N2 assays (suggesting it contained a non-Omicron SARS-CoV-2 lineage) and the other one was negative for the CDC assays. Cq data for all samples are summarized in Figure 3 and Supplementary Table 7.
Figure 3. Box plot of Ct values of omicron sample testing results. Cq values for clinical samples tested with Triplex 4 and the CDC N1 and N2 assays. The total number of samples is 1,548. The positive sample counts for each assay are: CDC N1 (1,013), CDC N2 (1,014), S:371–375 (1,255), S:142–144del (1,213), and S:ins214 (1,189). Gray circles represent outliers. Not all samples were tested with the CDC assays.
High-throughput clinical sample testing with Triplex 4
A total of 1,548 clinical samples were tested with Triplex 4. Of these, 1,255 (81%) were positive for one or more of the Omicron triplex assays. The sample and Cq values of these samples are presented in Supplementary Table 8 and the results are shown in the box plot (Figure 3).
Performance assessment of all triplexes on the Biomeme Franklin
Synthetic gRNA obtained from Twist Biosciences for six lineages (B.1.1.7, B.1.351, P.1, B.1.617.2, BA.1, and BA.2) were used to test performance of the variant triplexes on Biomeme FranklinTM devices. The Biomeme FranklinTM is a handheld, battery-powered nucleic acid testing device with 9 wells and 3-color detection that is designed for field portability in low-resource, far-forward environments (39). Negative controls were also included for each assay/run and consisted of water-only NTC reactions, as well as synthetic gRNA from the wild-type (A) lineage. We initially tested each lineage at a concentration of 500 genome equivalent copies (GEC) per reaction. At this level, all assays produced easily discernable amplification for all expected mutations in each lineage, with two exceptions: the S:P26S and S:N501Y assays in which Triplex 1 yielded only dim amplification curves that were not consistently detected (Figure 4A). When Triplex 1 was retested at a 100-fold higher concentration (50,000 GEC per reaction), these assays did produce positive results, but amplification curves remained dim (Figure 4B). These results are consistent with other results obtained on the Bio-Rad CFX96 and ABI 7500, with the S:P26S and S:N501Y assays producing dimmer curves and demonstrating less sensitive detection (higher LoDs) than the other assays in the panel, nonetheless detectable. All NTC reactions and other expected negatives were negative across all runs.
Figure 4. Amplification plots of triplexes 1–4 on the Biomeme FranklinTM device. (A) Plots of templates with 500 genomic copy equivalents of different variants. (B) Amplification plot of P1 template at 50,000 copies.
Omicron pentaplex analytical testing
Omicron pentaplex testing was performed using serial dilutions (50 to 50,000 copies per reaction) of gRNA for BA.1 and BA.2, and of gBlocks for BA.1.1, BA.3, BA.4, and BA.5 due to unavailability of gRNA for those lineages. Wild-type gRNA was used as a negative control. For comparison, these serial dilutions were also tested with Triplexes 1–4. All results were as expected based on in silico predictions and known variant mutation profiles according to outbreak.info (40). For positive samples, average Cq values (n = 2) are provided for the 50 copies per reaction test level (Table 6). Use of either the Omicron pentaplex or a combination of Triplexes 2–4 allowed for discrimination among BA.1/BA.1.1, BA.2, BA.3, and BA.4/BA.5, whereas use of only the original Omicron triplex (Triplex 4) would not have allowed for discrimination of BA.2 from BA.4 and BA.5. The BA.4 and BA.5 sublineages (as well as some sublineages of BA.2) also have an additional marker of clinical significance, S:L452R, which is part of Triplex 3 and can be used for further discriminatory power (47).
Discussion
The emergence of variants with phenotypic attributes may have impacts on the effectiveness of vaccines, therapeutics, and diagnostics. The resulting mutant proteins can lead to diminished molecular recognition for vaccine and therapeutic interventions (22–25). In regard to diagnostics, signature erosion is a term used for defining the mutational changes in molecular assay signatures such as primer/probe binding sites (48, 49). The term “signature erosion” is used here to signify potential false positive or false negative results in molecular assays due to mutations in the PCR primers/probe/amplicon target sequences (PCR signatures).
In this study, we describe a panel of multiplex RT-qPCR assays targeting a set of mutations across the major SARS-CoV-2 variants of concern; viz., Alpha, Beta, Gamma, Delta, and Omicron for use in a point-of-care PCR device as well as other standard PCR platforms. The Omicron triplex (Triplex 4) can distinguish among Omicron sublineages BA.1, BA.2, and BA.3, and the Omicron pentaplex can further delineate BA.4 and BA.5. We determined a sensitivity and specificity of >90% for these panels with sequence-confirmed SARS-CoV-2 clinical specimens. Based on these data, we posit that these panels could be valuable tools to fill variant testing needs and serve as a screening test to inform selection/prioritization of samples for whole genome sequencing. One of the drawbacks of these panels is that in the current configuration, the panels are used as non-Dx assays because of the cost and time involved in Dx development and are meant to be used as Research Use Only (RUO) assays. The general strategy is to test clinical samples using the Biomeme EUA Dx assay first, followed by screening for variants using the RUO genotyping panel. As of October 2022, there are no EUA SARS-CoV-2 genotyping RT-qPCR assays, although there are a couple with a sequencing component as part of the assay (4).
There are several distinct advantages to using SARS-CoV-2 variant genotyping PCR assays as a screening tool. (1) PCR is extremely less complex, cheaper and faster compared to sequencing, meaning a PCR panel could be a better alternative for rapid identification and tracking of variants (50). (2) A genotyping PCR panel could be very useful in resource-limited settings, since sequencing is not done at the same level in all locations/countries, and many laboratories do not have the infrastructure or expertise to perform sequencing (51). (3) Obtaining whole genome sequence data from clinical samples is very much dependent on viral load; reported PCR Cq cutoffs for obtaining full-length genome sequences range from ~25–33 (52, 53). A PCR panel could fill this gap by providing genotype information for numerous samples that cannot be sequenced due to low viral load as reflected by high Cq values (31). (4) A variant PCR panel can be used as an efficient and cost-effective screening assay to prioritize samples for sequencing, conserve resources, and minimize over/under-sampling of certain lineages for sequencing (54). (5) Mutation-specific PCR amplification can be a more precise option than a proxy based on assay failure, such as the well-known S gene target failure (SGTF). Unlike SGTF, in which negative results are used as a proxy to indicate presence of variants like Alpha and Omicron BA.1, this variant PCR panel comprises specific assays that detect only the mutant template and not the wild type, so variant detection is based on positive results rather than a lack of results. (6) During the course of the pandemic, multiple “discrete” waves of variants emerged. Hence, unique variant panels can be very useful to detect variants with no ambiguity when there are multiple variants in the community; the same is true for mixed infections (55). (7) Wastewater testing is a good predictor of impending community prevalence, with detection of new variants a few days to a week in advance of community prevalence (56), but sequencing wastewater samples to monitor the spread of genotypes is problematic/difficult due to sample quality and the presence of many different strains in the same sample, often at low-titer (57). Thus, a highly sensitive PCR panel that can detect variant markers would be ideal for wastewater monitoring (58). Similarly, the variant panel may also be a useful, cheaper alternative to sequencing and in vitro diagnostic (IVD) PCR assays for other types of environmental/non-clinical testing (59) (e.g., contamination/decontamination verification) (60), as RUO assays are typically less expensive than IVD assays. Although not the focus of this study, very limited preliminary studies have also been conducted using our RT-qPCR assays to test for SARS-CoV-2 variants in wastewater samples (unpublished data).
There are of course some limitations of using RT-qPCR panels compared to other approaches such as sequencing. The panels were designed in a way to detect specific variants using a combination of mutations indicative of the variant in the sample. However, these panels may not be useful to detect newly emerging variants with mutations that are not targeted by the panels. Similarly, these panels may not be useful for detecting recombinant variants, since these are potentially chimeras between variants. Also, convergent evolution of the same mutations in multiple lineages may also confound the accurate adjudication of PCR results and assignment of variants present in a sample. In the long run, it would be advantageous to develop assays targeting specific mutations and customize panels depending on the mutations for a given lineage or sublineage in question, and for recombinant lineages (e.g., PCR targeting the junctions or recombination breakpoints). The problem with such a flexible multiplexing approach is that there may not be enough validation data on the performance of a custom triplex with a given set of mutations, unless such data are generated during assay development. In that case, triplex combination testing will become enormous depending on the number of mutations. Hence, the panels described here were designed, tested and validated to detect these specific variants only, and may not work for other variants with a different constellation of mutations, or for recombinant variants. Also, the variants appeared in distinct waves during the course of the pandemic and got replaced by the next wave in a very short time frame. In some cases, the assay development, testing and validation took longer than the time it took for wave replacement, so when the new triplex was ready to be tested in clinical samples, the variant had already been replaced by a newer variant with a new set of mutations. Such an “arms race” between variant emergence and assay development is a challenge that needs to be addressed so that new panels do not become obsolete even before they can be fully developed and verified. Thus, the best approach may be to have assays for different mutations tested and ready to incorporate into customized triplexes, rather than targeting specific lineages that may rise and fall quickly.
Despite these limitations, these variant panels can help fulfill the need for rapid identification of the SARS-CoV-2 variant present in a sample, leading to quick decision making with respect to public health measures as well as treatment options for individuals, especially when sequencing is not an option. In some areas—for instance, many regions in Africa—access to sequencing is extremely limited, while PCR is much more widely available (61, 62). Although progress has been made toward globalizing access to sequencing technology, there are still significant disparities among regions and countries. One reason for this is the significant capability standup costs, including the high cost of next-generation sequencers, accessory lab equipment, data storage and analysis infrastructure, and non-equipment capability standup costs like staff training, protocol development, and method verification studies (63). In addition to the standup costs, another hurdle to broad/universal application of sequencing is the relatively high cost per sample of whole genome sequencing relative to PCR and other nucleic acid amplification tests (NAATs). The RT-qPCR triplexes we developed cost ~$2.83 per reaction in the “wet” (non-lyophilized) format we used for high-throughput testing of clinical specimens, whereas per-sample reagent costs for non-targeted whole genome sequencing are often over $100 depending on instrument and samples-per-run. Although it is possible to achieve per-sample sequencing reagent costs under $10 when using a targeted amplicon sequencing approach and combining a high number of samples (96 or more) on the same sequencing run (64), PCR is still more economical in most cases, especially after accounting for the cost of labor hours, equipment purchases, and other startup costs. Thus, although whole genome sequencing represents a more comprehensive and informative test, multiplexed genotyping RT-qPCR assays like those developed in this study can be invaluable tools when sequencing is not feasible/practical because of access, cost, or time constraints.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
Conceptualization: Jv, JS, MN, and SS. Experimental design: Jv, MN, EW, JS, and SS. Methodology and experimental validation: BS, MN, BP, BK, PD, RD, RP, and AF. Formal analysis: BS, MN, BP, BK, PD, and JS. Data curation and visualization: JS and BP. Resources: CM, Jv, and JS. Project administration and supervision: Jv, JS, and SS. Funding acquisition: SS and BG. Writing—original draft preparation: SS, JS, and MN. Writing—review and editing: All authors. All authors contributed to the article and approved the submitted version.
Funding
Funding for this study was provided by the Department of Defense's Joint Program Executive Office for Chemical, Biological, Radiological and Nuclear Defense (JPEO-CBRND), in collaboration with the Defense Health Agency (DHA) COVID funding initiative for this effort: Development of a Biomeme Variant Panel under the Contract (FA807518D0017) with MRIGlobal, and a Subcontract (867-111084-20) between MRIGlobal and Biomeme, Inc. Work performed by Naval Health Research Center (NHRC) was also funded via DBPAO contract (N2016).
Acknowledgments
We would like to acknowledge Eric Lai and the NIH RADx program (https://www.nih.gov/research-training/medical-research-initiatives/radx) for providing some of the clinical samples used in this study. We would also like to thank Robert Player for help generating the Figure 3 box plot.
Conflict of interest
Authors BS, MN, and Jv were employed by company Biomeme, Inc. Author EW was employed by Integrated DNA Technologies.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author disclaimer
The content of this publication are the authors' opinions and do not necessarily reflect the views or policies of the U.S. Department of Defense, the U.S. Department of the Army, or the institutions and companies affiliated with the authors, nor does the mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government. Some authors of this work are military service members or federal/contracted employees of the United States government. This work was prepared as part of their official duties. Title 17 U.S.C. 105 provides that copyright protection under this title is not available for any work of the United States Government. Title 17 U.S.C. 101 defines a U.S. Government work as work prepared by a military service member or employee of the U.S. Government as part of that person's official duties.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.1042647/full#supplementary-material
References
1. COVID-19 Public Health Emergency of International Concern (PHEIC) Global research and innovation forum. https://www.who.int/publications/m/item/covid-19-public-health-emergency-of-international-concern-(pheic)-global-research-and-innovation-forum (accessed February 12, 2020).
2. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. (2020) 382:727–33. doi: 10.1056/NEJMoa2001017
3. Wu F, Zhao S, Yu B, Chen YM, Wang W, Song ZG, et al. A new coronavirus associated with human respiratory disease in China. Nature. (2020) 579:265–9. doi: 10.1038/s41586-020-2008-3
4. In Vitro Diagnostics EUAs - Molecular Diagnostic Tests for SARS-CoV-2. Available online at: https://www.fda.gov/medical-devices/coronavirus-disease-2019-covid-19-emergency-use-authorizations-medical-devices/in-vitro-diagnostics-euas-molecular-diagnostic-tests-SARS-CoV-2 (press release May 2022).
5. Ma Q, Liu J, Liu Q, Kang L, Liu R, Jing W, et al. Global percentage of asymptomatic SARS-CoV-2 infections among the tested population and individuals with confirmed COVID-19 diagnosis: a systematic review and meta-analysis. JAMA Netw Open. (2021) 4:e2137257. doi: 10.1001/jamanetworkopen.2021.37257
6. Parikh A, Cooper L, Frogel D, Le Benger K, Cooper CK, Parvu V. Large-scale SARS-CoV-2 antigen testing with real-world specimens. Front Public Health. (2022) 10:836328. doi: 10.3389/fpubh.2022.836328
7. Larremore DB, Wilder B, Lester E, Shehata S, Burke JM, Hay JA, et al. Test sensitivity is secondary to frequency and turnaround time for COVID-19 screening. Sci Adv. (2021) 7:eabd5393. doi: 10.1126/sciadv.abd5393
8. Mina MJ, Parker R, Larremore DB. Rethinking Covid-19 test sensitivity—a strategy for containment. N Engl J Med. (2020) 383:e120. doi: 10.1056/NEJMp2025631
9. Kissler SM, Fauver JR, Mack C, Olesen SW, Tai C, Shiue KY, et al. Viral dynamics of acute SARS-CoV-2 infection and applications to diagnostic and public health strategies. PLoS Biol. (2021) 19:e3001333. doi: 10.1371/journal.pbio.3001333
10. Smith RL, Gibson LL, Martinez PP, Ke R, Mirza A, Conte M, et al. Longitudinal assessment of diagnostic test performance over the course of acute SARS-CoV-2 infection. J Infect Dis. (2021) 224:976–82. doi: 10.1093/infdis/jiab337
11. Miller MF, Shi M, Motsinger-Reif A, Weinberg CR, Miller JD, Nichols E. Community-based testing sites for SARS-CoV-2 - United States, March 2020–November 2021. MMWR Morb Mortal Wkly Rep. (2021) 70:1706–11. doi: 10.15585/mmwr.mm7049a3
12. At-Home OTC. COVID-19 Diagnostic Tests. Available online at: https://wwwfdagov/medical-devices/coronavirus-covid-19-and-medical-devices/home-otc-covid-19-diagnostic-tests (press release, 2022).
13. Lambrou AS, Shirk P, Steele MK, Paul P, Paden CR, Cadwell B, et al. Genomic surveillance for SARS-CoV-2 variants: predominance of the delta (B.1.617.2) and Omicron (B.1.1.529) variants - United States, June 2021–January 2022. MMWR Morb Mortal Wkly Rep. (2022) 71:206–11. doi: 10.15585/mmwr.mm7106a4
15. Chen Z, Azman AS, Chen X, Zou J, Tian Y, Sun R, et al. Global landscape of SARS-CoV-2 genomic surveillance and data sharing. Nat Genet. (2022) 54:499–507. doi: 10.1038/s41588-022-01033-y
16. Khare S, Gurry C, Freitas L, Schultz MB, Bach G, Diallo A, et al. GISAID's role in pandemic response. China CDC Wkly. (2021) 3:1049–51. doi: 10.46234/ccdcw2021.255
17. Elbe S, Buckland-Merrett G. Data, disease and diplomacy: GISAID's innovative contribution to global health. Glob Chall. (2017) 1:33–46. doi: 10.1002/gch2.1018
18. Shu Y, McCauley J. GISAID: global initiative on sharing all influenza data - from vision to reality. Euro Surveill. (2017) 22:30494. doi: 10.2807/1560-7917.ES.2017.22.13.30494
19. Volz E, Mishra S, Chand M, Barrett JC, Johnson R, Geidelberg L, et al. Assessing transmissibility of SARS-CoV-2 lineage B. 117 in England Nature. (2021) 593:266–9. doi: 10.1038/s41586-021-03470-x
20. Vohringer HS, Sanderson T, Sinnott M, De Maio N, Nguyen T, Goater R, et al. Genomic reconstruction of the SARS-CoV-2 epidemic in England. Nature. (2021) 600:506–11. doi: 10.1101/2021.05.22.21257633
21. Alpert T, Brito AF, Lasek-Nesselquist E, Rothman J, Valesano AL, MacKay MJ, et al. Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States. Cell. (2021) 184:2595–604 e13. doi: 10.1016/j.cell.2021.03.061
22. Harvey WT, Carabelli AM, Jackson B, Gupta RK, Thomson EC, Harrison EM, et al. SARS-CoV-2 variants, spike mutations and immune escape. Nat Rev Microbiol. (2021) 19:409–24. doi: 10.1038/s41579-021-00573-0
23. Liu L, Iketani S, Guo Y, Chan JF, Wang M, Liu L, et al. Striking antibody evasion manifested by the Omicron variant of SARS-CoV-2. Nature. (2022) 602:676–81. doi: 10.1038/s41586-021-04388-0
24. Iketani S, Liu L, Guo Y, Liu L, Chan JF, Huang Y, et al. Antibody evasion properties of SARS-CoV-2 Omicron sublineages. Nature. (2022) 604:553–6. doi: 10.1038/s41586-022-04594-4
25. Cavazzoni P. Coronavirus (COVID-19) Update: FDA Limits Use of Certain Monoclonal Antibodies to treat COVID-19 Due to the Omicron Variant. Silver Spring, MD: US Food and Drug Administration (2022).
26. Lai E, Kennedy EB, Lozach J, Hayashibara K, Davis-Turak J, Becker D, et al. A method for variant agnostic detection of SARS-CoV-2, rapid monitoring of circulating variants, and early detection of emergent variants such as omicron. J Clin Microbiol. (2022) 7:e0034222. doi: 10.1128/jcm.00342-22
27. Vogels CBF, Breban MI, Ott IM, Alpert T, Petrone ME, Watkins AE, et al. Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2. PLoS Biol. (2021) 19:e3001236. doi: 10.1371/journal.pbio.3001236
28. Neopane P, Nypaver J, Shrestha R, Beqaj SS. SARS-CoV-2 variants detection using TaqMan SARS-CoV-2 mutation panel molecular genotyping assays. Infect Drug Resist. (2021) 14:4471–9. doi: 10.2147/IDR.S335583
29. Dikdan RJ, Marras SAE, Field AP, Brownlee A, Cironi A, Hill DA, et al. Multiplex PCR assays for identifying all major severe acute respiratory syndrome coronavirus 2 variants. J Mol Diagn. (2022) 24:309–19. doi: 10.1016/j.jmoldx.2022.01.004
30. Migueres M, Lhomme S, Tremeaux P, Dimeglio C, Ranger N, Latour J, et al. Evaluation of two RT-PCR screening assays for identifying SARS-CoV-2 variants. J Clin Virol. (2021) 143:104969. doi: 10.1016/j.jcv.2021.104969
31. Wang H, Jean S, Eltringham R, Madison J, Snyder P, Tu H, et al. Mutation-specific SARS-CoV-2 PCR screen: rapid and accurate detection of variants of concern and the identification of a newly emerging variant with spike L452R mutation. J Clin Microbiol. (2021) 59:e0092621. doi: 10.1128/JCM.00926-21
32. Harper H, Burridge A, Winfield M, Finn A, Davidson A, Matthews D, et al. Detecting SARS-CoV-2 variants with SNP genotyping. PLoS ONE. (2021) 16:e0243185. doi: 10.1371/journal.pone.0243185
33. Bechtold P, Wagner P, Hosch S, Siegrist D, Ruiz-Serrano A, Gregorini M, et al. Rapid identification of SARS-CoV-2 variants of concern using a portable peakPCR platform. Anal Chem. (2021) 93:16350–9. doi: 10.1021/acs.analchem.1c02368
34. Methods for the detection and characterisation of SARS-CoV-2 variants - first update. Available online at: https://www.ecdc.europa.eu/en/publications-data/methods-detection-and-characterisation-SARS-CoV-2-variants-first-update (press release, 2021).
35. Vega-Magana N, Sanchez-Sanchez R, Hernandez-Bello J, Venancio-Landeros AA, Pena-Rodriguez M, Vega-Zepeda RA, et al. RT-qPCR assays for rapid detection of the N501Y, 69-70del, K417N, and E484K SARS-CoV-2 mutations: a screening strategy to identify variants with clinical impact. Front Cell Infect Microbiol. (2021) 11:672562. doi: 10.3389/fcimb.2021.672562
36. De Pace V, Bruzzone B, Orsi A, Ricucci V, Domnich A, Guarona G, et al. Comparative analysis of five multiplex RT-PCR assays in the screening of SARS-CoV-2 variants. Microorganisms. (2022) 10:306. doi: 10.3390/microorganisms10020306
37. Vogels CBF, Alpert T, Breban M, Fauver JR, Grubaugh ND. (2021). Available online at: https://virological.org/t/multiplexed-rt-qpcr-to-screen-for-SARS-CoV-2-b-1-1-7-variants-preliminary-results/588 (accessed March 15, 2021).
38. Kaplinski L, Andreson R, Puurand T, Remm M. MultiPLX: automatic grouping and evaluation of PCR primers. Bioinformatics. (2005) 21:1701–2. doi: 10.1093/bioinformatics/bti219
39. Zowawi HM, Alenazi TH, AlOmaim WS, Wazzan A, Alsufayan A, Hasanain RA, et al. Portable RT-PCR system: a rapid and scalable diagnostic tool for COVID-19 testing. J Clin Microbiol. (2021) 59:e03004–20. doi: 10.1128/JCM.03004-20
40. Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JA. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. (2007) 35:W71–4. doi: 10.1093/nar/gkm306
41. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. (2012) 13:134. doi: 10.1186/1471-2105-13-134
42. ThermoFisher Scientific. Multiple Primer Analyzer. Available online at: https://www.thermofisher.com/us/en/home/brands/thermo-scientific/molecular-biology/molecular-biology-learning-center/molecular-biology-resource-library/thermo-scientific-web-tools/multiple-primer-analyzer.html (accessed December 02, 2021).
43. Hatcher EL, Zhdanov SA, Bao Y, Blinkova O, Nawrocki EP, Ostapchuck Y, et al. Virus variation resource—improved response to emergent viral outbreaks. Nucleic Acids Res. (2017) 45:D482–D90. doi: 10.1093/nar/gkw1065
44. Rambaut A, Holmes EC, O'Toole A, Hill V, McCrone JT, Ruis C, et al. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat Microbiol. (2020) 5:1403–7. doi: 10.1038/s41564-020-0770-5
45. Gardner SN, Slezak T. Simulate_PCR for amplicon prediction and annotation from multiplex, degenerate primers and probes. BMC Bioinformatics. (2014) 15:237. doi: 10.1186/1471-2105-15-237
46. Latif AA, Mullen JL, Alkuzweny M, Tsueng G, Cano M, Haag E, et al. Variant reports. Available online at: https://outbreak.info/situation-reports (accessed April 29, 2022).
47. Cao Y, Yisimayi A, Jian F, Song W, Xiao T, Wang L, et al. BA.2.12.1, BA.4 and BA.5 escape antibodies elicited by Omicron infection. Nature. (2022) 602: 593–602. doi: 10.1038/s41586-022-04980-y
48. Slezak T, Kuczmarski T, Ott L, Torres C, Medeiros D, Smith J, et al. Comparative genomics tools applied to bioterrorism defence. Brief Bioinform. (2003) 4:133–49. doi: 10.1093/bib/4.2.133
49. Gardner SN, Lam MW, Smith JR, Torres CL, Slezak TR. Draft vs. finished sequence data for DNA and protein diagnostic signature development. Nucleic Acids Res. (2005) 33:5838–50. doi: 10.1093/nar/gki896
50. Methods for the detection and characterisation of SARS-CoV-2 variants - second update. Stockholm and Copenhagen: European Centre for Disease Prevention and Control/World Health Organization Regional Office for Europe (2022).
51. Ranasinghe D, Jayadas TTP, Jayathilaka D, Jeewandara C, Dissanayake O, Guruge D, et al. Comparison of different sequencing techniques for identification of SARS-CoV-2 variants of concern with multiplex real-time PCR. PLoS ONE. (2022) 17:e0265220. doi: 10.1371/journal.pone.0265220
52. SeqCOVID. Evaluating qPCR Cycle Threshold (Ct) As a Predictive Value of Sequencing Quality. (2022). Available online at: https://seqcovid.csic.es/evaluating-qpcr-cycle-threshold-ct-as-a-predictive-value-of-sequencing-quality/ (accessed February 23, 2022).
53. Doddapaneni H, Cregeen SJ, Sucgang R, Meng Q, Qin X, Avadhanula V, et al. Oligonucleotide capture sequencing of the SARS-CoV-2 genome and subgenomic fragments from COVID-19 individuals. PLoS ONE. (2021) 16:e0244468. doi: 10.1371/journal.pone.0244468
54. Borillo GA, Kagan RM, Marlowe EM. Rapid and accurate identification of SARS-CoV-2 variants using real time PCR assays. Front Cell Infect Microbiol. (2022) 12:894613. doi: 10.3389/fcimb.2022.894613
55. Thakur V, Bhola S, Thakur P, Patel SKS, Kulshrestha S, Ratho RK, et al. Waves and variants of SARS-CoV-2: understanding the causes and effect of the COVID-19 catastrophe. Infection. (2022) 50:309–25. doi: 10.1007/s15010-021-01734-2
57. Wilhelm A, Schoth J, Meinert-Berning C, Agrawal S, Bastian D, Orschler L, et al. Wastewater surveillance allows early detection of SARS-CoV-2 omicron in North Rhine-Westphalia, Germany. Sci Total Environ. (2022) 846:157375. doi: 10.1016/j.scitotenv.2022.157375
58. Gregory DA, Trujillo M, Rushford C, Flury A, Kannoly S, San KM, et al. Genetic diversity and evolutionary convergence of cryptic SARS- CoV-2 lineages detected via wastewater sequencing. PLoS Pathog. (2022) 18:e1010636. doi: 10.1371/journal.ppat.1010636
59. Maltezou HC, Tseroni M, Daflos C, Anastassopoulou C, Vasilogiannakopoulos A, Daligarou O, et al. Environmental testing for SARS-CoV-2 in three tertiary-care hospitals during the peak of the third COVID-19 wave. Am J Infect Control. (2021) 49:1435–7. doi: 10.1016/j.ajic.2021.08.022
60. Glinert I, Ben-Shmuel A, Szwartcwort-Cohen M, Beth-Din A, Laskar O, Barlev-Gross M, et al. Revisiting SARS-CoV-2 environmental contamination by patients with COVID-19: The Omicron variant does not differ from previous strains. Int J Infect Dis. (2022) 118:211–3. doi: 10.1016/j.ijid.2022.03.001
61. Adebisi YA, Oke GI, Ademola PS, Chinemelum IG, Ogunkola IO, Lucero-Prisno Iii DE. SARS-CoV-2 diagnostic testing in Africa: needs and challenges. Pan Afr Med J. (2020) 35:4. doi: 10.11604/pamj.2020.35.4.22703
62. Ondoa P, Kebede Y, Loembe MM, Bhiman JN, Tessema SK, Sow A, et al. COVID-19 testing in Africa: lessons learnt. Lancet Microbe. (2020) 1:e103–e4. doi: 10.1016/S2666-5247(20)30068-9
63. Yohe S, Thyagarajan B. Review of clinical next-generation sequencing. Arch Pathol Lab Med. (2017) 141:1544–57. doi: 10.5858/arpa.2016-0501-RA
Keywords: SARS-CoV-2, variants, PCR, Biomeme, rapid detection, screening assay, RUO
Citation: Stanhope BJ, Peterson B, Knight B, Decadiz RN, Pan R, Davis P, Fraser A, Nuth M, vanWestrienen J, Wendlandt E, Goodwin B, Myers C, Stone J and Sozhamannan S (2022) Development, testing and validation of a SARS-CoV-2 multiplex panel for detection of the five major variants of concern on a portable PCR platform. Front. Public Health 10:1042647. doi: 10.3389/fpubh.2022.1042647
Received: 12 September 2022; Accepted: 11 November 2022;
Published: 15 December 2022.
Edited by:
Sergio E. Rodriguez, Centers for Disease Control and Prevention (CDC), United StatesReviewed by:
Subhra Subhadra, Regional Medical Research Center (ICMR), IndiaConnie Lam, Western Sydney Local Health District, Australia
Copyright Copyright © 2022 Stanhope, Peterson, Knight, Decadiz, Pan, Davis, Fraser, Nuth, vanWestrienen, Wendlandt, Goodwin, Myers, Stone and Sozhamannan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The authors grant nonexclusive right to copy, distribute, adapt, and transmit the published Work for commercial or non-commercial use with proper attribution under the Creative Commons Attribution 4.0 International license (CC-BY). No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jennifer Stone, anN0b25lJiN4MDAwNDA7bXJpZ2xvYmFsLm9yZw==; Shanmuga Sozhamannan, c2hhbm11Z2Euc296aGFtYW5uYW4uY3RyJiN4MDAwNDA7YXJteS5taWw=