key: cord-0967718-pgvwjufy
authors: Boršová, Kristína; Paul, Evan D.; Kováčová, Viera; Radvánszka, Monika; Hajdu, Roman; Čabanová, Viktória; Sláviková, Monika; Ličková, Martina; Lukáčiková, Ľubomíra; Belák, Andrej; Roussier, Lucia; Kostičová, Michaela; Líšková, Anna; Maďarová, Lucia; Štefkovičová, Mária; Reizigová, Lenka; Nováková, Elena; Sabaka, Peter; Koščálová, Alena; Brejová, Broňa; Staroňová, Edita; Mišík, Matej; Vinař, Tomáš; Nosek, Jozef; Čekan, Pavol; Klempa, Boris
title: Surveillance of SARS-CoV-2 lineage B.1.1.7 in Slovakia using a novel, multiplexed RT-qPCR assay
date: 2021-10-14
journal: Sci Rep
DOI: 10.1038/s41598-021-99661-7
sha: 2b5214e3970701de1eb320fb10adeb19bcd10c18
doc_id: 967718
cord_uid: pgvwjufy

The emergence of a novel SARS-CoV-2 B.1.1.7 variant sparked global alarm due to increased transmissibility, mortality, and uncertainty about vaccine efficacy, thus accelerating efforts to detect and track the variant. Current approaches to detect B.1.1.7 include sequencing and RT-qPCR tests containing a target assay that fails or results in reduced sensitivity towards the B.1.1.7 variant. Since many countries lack genomic surveillance programs and failed assays detect unrelated variants containing similar mutations as B.1.1.7, we used allele-specific PCR, and judicious placement of LNA-modified nucleotides to develop an RT-qPCR test that accurately and rapidly differentiates B.1.1.7 from other SARS-CoV-2 variants. We validated the test on 106 clinical samples with lineage status confirmed by sequencing and conducted a country-wide surveillance study of B.1.1.7 prevalence in Slovakia. Our multiplexed RT-qPCR test showed 97% clinical sensitivity and retesting 6,886 SARS-CoV-2 positive samples obtained during three campaigns performed within one month, revealed pervasive spread of B.1.1.7 with an average prevalence of 82%. Labs can easily implement this test to rapidly scale B.1.1.7 surveillance efforts and it is particularly useful in countries with high prevalence of variants possessing only the ΔH69/ΔV70 deletion because current strategies using target failure assays incorrectly identify these as putative B.1.1.7 variants.

Identification of RT-qPCR targets by bioinformatic analysis. Our analysis of 1,136 spike gene sequences (spanning 1-21 December 2020) revealed 228 sequences (20%) that contained both the ΔH69/ΔV70 and ΔY144 deletions (for country of origin, see Supplementary Table S1 ). The shorter deletion (ΔY144) always co-occurred with the longer deletion (ΔH69/ΔV70), whereas the (ΔH69/ΔV70) deletion occurs independently in 17 sequences (1.5%). Pearson's correlation coefficient of the deletions is 0.953.

We analysed over 900,000 SARS-CoV-2 genomes to determine the prevalence of both ΔH69/ΔV70 and ΔY144 deletions in lineage B.1.1.7 and lineages other than B.1.1.7 and found a total of 257,953 sequences that possess both deletions. Based on the metadata file, we identified SARS-CoV-2 lineages across all called sequences with both deletions. Only 1,230 sequences (0.48%) out of 257,953 sequences are not labelled as B.1.1.7. In other words, 99.52% of sequences containing both deletions belong to lineage B.1.1.7, highlighting the notion that these two deletions are highly specific for the B.1.1.7 variant and make ideal targets for primer/probe design ( Table 1) . www.nature.com/scientificreports/ two assays targeting either the wild-type SARS-CoV-2 spike gene or the ΔH69/ΔV70 and ΔY144 deletions that are highly specific to lineage B.1.1.7 (for primer/probe locations and sequences see Fig. 1 , Table 2 , respectively). We began by designing a general S gene primer/probe assay (SARS-CoV-2 S gene) that could be used for screening purposes and would detect the most common strains of SARS-CoV-2 as well as variants containing ΔH69/ΔV70, including lineage B.1.1.7. We assessed the performance of a series of primers flanking the ΔH69/ ΔV70 deletion ( Supplementary Fig. S1a ,b; for primer/probe sequences, see Supplementary Table S2 ). After selecting the optimal primer/probe combination, we introduced an additional hydrolysis probe, identically labelled with the same reporter and quencher dyes, that would hybridize in tandem (i.e., on the same strand) with the first hydrolysis probe. Consistent with other reports 12, 13 , this dual probe approach enhanced sensitivity and specifity ( Supplementary Fig. S1c,d) .

For our assay targeting lineage B.1.1.7, we leveraged the highly specific co-occurrence of the ΔH69/ΔV70 and ΔY144 deletions in lineage B.1.1.7 (99.52%, Table 1 ). By designing a series of forward primers to target the ΔH69/ ΔV70 deletion, we differentiated wild type template from ΔH69/ΔV70 template ( Supplementary Fig. S2a ). Since other SARS-CoV-2 variants share the ΔH69/ΔV70 deletion (e.g., B.1.1.298, B.1.160, B.1.177, B.1.258, B.1.375, B.1.525), we designed a series of reverse primers to target the second, three base pair deletion (bp: 21,991-21,993; ΔY144) and utilized allele-specific PCR approaches [14] [15] [16] and judicious placement of LNA-modified nucleotides to enhance the specificity of the assay ( Supplementary Fig. S2b,c) . This approach enabled us to differentiate B.1.1.7 variants that contain both the ΔH69/ΔV70 and ΔY144 deletions from SARS-CoV-2 variants that contain only the ΔH69/ΔV70 deletion, provided that a second reaction is ran in parallel using the SARS-CoV-2 S gene set that can be used as a benchmark to assess the relative sensitivity. If the B.1.1.7 primer set amplifies the sample within five Ct cycles of the SARS-CoV-2 S gene primer set, then the sample is B.1.1.7 positive. Alternatively, if the B.1.1.7 primer set amplifies the sample in 8 or more Ct cycles relative to the SARS-CoV-2 S gene primer set, than the sample likely belongs to a variant that contains the ΔH69/ΔV70 deletion, but not the ΔY144 deletion, and hence is B.1.1.7 negative.

We compared three different versions of B.1.1.7 primer/probe sets using a selected set of 46 samples, some of which were sequenced to confirm lineage status. Given our interpretation criterion (Supplementary Table S4 ), we determined that the V3 primer/probe set performed the best since it correctly identified all B.1.1.7 and ΔH69/ΔV70 deletion samples, with the exception of one ΔH69/ΔV70 deletion sample that was interpreted as inconclusive ( Supplementary Fig. S3 ).

Analytical sensitivity and clinical evaluation of lineage B.1.1.7 S gene primer/probe set. With our final primer/probe sets for SARS-CoV-2 S gene and B.1.1.7 (Table 2) , we multiplexed each assay with the US CDC human RNase P primer/probe set (for sequences, see 17 ) as an internal control to assess RNA extraction and assay performance. We then determined the analytical sensitivity using serial dilutions of RNA extracted from a B.1.1.7 positive sample. Both assays displayed high sensitivity (Fig. 2a) with our SARS-CoV-2 S gene and www.nature.com/scientificreports/ B.1.1.7 assays reliably detecting down to only 2 copies/reaction (0.4 copies/μl) and 10 copies/reaction (2 copies/ μl), respectively, placing them among the most sensitive SARS-CoV-2 RT-qPCR assays available. We evaluated the clinical performance of our SARS-CoV-2 S gene and B.1.1.7 assays on 106 clinical samples that underwent sequencing to identify lineage status using interpretation criterion outlined in Supplementary  Table S4 . A detailed flow chart of the patient samples is depicted in Supplementary Fig. S4 . Our SARS-CoV-2 S gene assay detected all 106 clinical samples regardless of lineage (Table 3) confirming its utility as a general screening assay for the most common SARS-CoV-2 variants. Out of 67 clinical samples classified as lineage B.1.1.7 by sequencing, our B.1.1.7 assay positively identified 65 samples, while one sample (Sample 75) was not detected by the B.1.1.7 assay and another sample was deemed inconclusive. The ΔCt of this sample was slightly greater than five cycles (e.g., Sample 40, ΔCt = 5.7) relative to the Ct value for the SARS-CoV-2 S gene assay, which exceeded our cut-off for a positive identification (Fig. 2b) .

Notably, our assay was also capable of identifying samples carrying the ΔH69/ΔV70 deletion such as those belonging to the B. www.nature.com/scientificreports/ not detected by the B.1.1.7 assay possibly due to relatively high Ct values in the SARS-CoV-2 S gene assay for two samples (Sample 33, Ct = 30.1 and Sample 65, Ct = 28.9), making confirmation of the ΔH69/ΔV70 status impossible with our cut-off criterion. One sample had a ΔCt outside the criterion for ΔH69/ΔV70 deletion confirmation and was deemed inconclusive. Neither the B.1.1.7 assay nor the common SARS-CoV-2 S gene assay yielded a positive signal when testing the 94 samples that were previously identified as SARS-CoV-2 negative, while RNase P, in both multiplexed assays, was positive for all 94 samples (Table 3 ). Both assays did not amplifiy any of the respiratory pathogens that were tested in a cross reactivity experiment (Supplementary Table S6 ). Overall, the clinical evaluation confirmed the diagnostic utility of both our SARS-CoV-2 S gene and B.1.1.7 assays, which showed 100% (106/106) and 97% (65/67) diagnostic sensitivity, respectively. The assay also showed diagnostic utility for identifying variants containing only the ΔH69/ΔV70 deletion by detecting 83.3% (20/24) of lineage B.1.258 samples. When considering B.1.258 samples containing high viral loads (Ct ≤ 28), which is necessary to identify variants with only the ΔH69/ΔV70 deletion, the diagnostic sensitivity reached 91.7% (22/24) . For an overview of the clinical evaluation data, lineage of each sample, and GISAID information, see Supplementary Table S7 . We have provided a decision tree (Fig. 2c ) that users may follow to implement this test to directly detect the presence of the B.1.1.7 variant.

To determine the prevalence of lineage B.1.1.7 and its spread throughout the Slovak Republic, regional public health authorities used this test to screen nearly 7,000 samples previously identified as SARS-CoV-2 positive during three screening rounds over a one month period. The first surveillance screening for lineage B.1.1.7 on February 2nd, 2021 revealed a region-specific prevalence of B.1.1.7 ranging from 52% (Žilina) to 85% (Trnava) with an overall B.1.1.7 prevalance of 75% ( Fig. 3a -c, Supplementary Table S8 ). During the second round of screening held on February 17th, 2021, the majority of regions (5 out of 8) showed increased prevalence, although due to reduced prevalence in Banská Bystrica, Košice, and Nitra, the overall prevalence remained at 74%. Closer scrutiny of the data suggests the reduced prevalence in some regions was caused by samples originating from several "non-B.1.1.7 clusters" (e.g., in social care homes), where all positive samples were likely derived from a common infection source. Exclusion of these clusters results in a slight increase in prevalence of 78%. A final screening on March 3rd, 2021 showed increased B.1.1.7 prevalence in all regions and an overall prevalence of 85%. Taking into consideration estimates of growth rates (i.e., reproduction number) for lineage B.1.1.7 in the UK 6,8 , we used the prevalences obtained on the first screening round (February 2nd, 2021) to estimate historical and future prevalences using a range of two weeks and applying spread factors (3.3 and 5.4) estimated in other countries (Fig. 3d ). Although the expected data did not fit with our observed data, a number of factors likely influence regional spread, including the presence and virulence of competing SARS-CoV-2 lineages and selection biases from cluster outbreaks.

The recent emergence of a novel SARS-CoV-2 variant called lineage B.1.1.7 has sparked global consternation as it has now been confirmed in over 70 countries and threatens to exacerbate an already dire pandemic. To mitigate the spread of the B.1.1.7 variant, it is imperative that countries have diagnostic tools that can quickly and accurately detect and track the prevalence of the variant in order to implement the appropriate epidemiological measures. Here we report a novel RT-qPCR test to differentiate the B.1.1.7 variant from other SARS-CoV-2 Table 3 . Clinical performance of SARS-CoV-2 S gene and B.1.1.7 primer/probes sets. 18, 20, [22] [23] [24] [25] . These RT-qPCR tests contain multigene assays, with at least one assay targeting the spike gene, and during routine testing a "drop-out" in the spike gene assay may occur (often termed as S gene target failure, or SGTF), while other gene targets yield positive signals. This SGTF can indicate the presence of the B.1.1.7 variant and flag samples for confirmation by sequencing. It is important to note, however, that SGTPs are produced by other variants that contain the ΔH69/ΔV70 deletion, including the B. 1.1, B1.258, B.1.525, [27] [28] [29] [30] [31] [32] . These tests are largely based on classical SNP genotyping methods using either probe-based genotyping or melting curve analyses and typically focus on mutations that are shared between many variants (e.g., N501Y, E484K, and ΔH69/ ΔV70 deletion), making it difficult to distinguish between variants unless running multiple SNP assays and then making complex comparisons of melt curves. Similar to SGTF tests, these approaches, while providing a rapid snapshot of the presence of SARS-CoV-2 variants, often require follow up genomic sequencing to identify the particular variant.

Several groups have described publicly available RT-qPCR protocols for detection of lineage B.1.1.7 that can be divided into SNP genotyping assays using either SYBR 33 or probe-based 34 melting curve analyses, multiplexed probe-based RT-qPCR [35] [36] [37] [38] [39] , and a combination of target drop-out tests 40 . While these open source RT-qPCR protocols offer cost-effective, rapid strategies to directly detect multiple SARS-CoV-2 variants, some are limited by only assessing a small number of mutations that preclude identification of specific variants and most have yet to be tested in real world surveillance scenarios.

To differentiate B.1.1.7, we took an alternative approach by targeting both the ΔH69/ΔV70 and ΔY144 deletions using allele-specific PCR methods combined with judicious placement of LNA oligonucleotides. Together, these modifications provided us with a primer/probe set that retained specificity for B.1.1.7 variants and reduced specificity to other variants containing only the ΔH69/ΔV70 deletion. To highlight the specificity of this assay, our analysis of all GISAID sequences containing both the ΔH69/ΔV70 and ΔY144 deletions revealed that a staggering 99.52% of all these sequences belong to lineage B.1.1.7, ensuring that users can have high confidence that a positive B.1.1.7 assay result is a true positive. Our test, instead of relying on target failures to identify putative variants, provides a positive signal in the presence of B.1.1.7 and ΔH69/ΔV70 deletion variants that can easily be differentiated by comparing their relative Ct values to a common SARS-CoV-2 S gene primer/probe set that serves as a benchmark.

Since this test depends on a relative comparison of Ct values between the B.1.1.7 and common SARS-CoV-2 S gene assays (i.e., B.1.1.7 positive samples must amplify within 5 Ct cycles of the common SARS-CoV-2 S gene), it is important that samples have sufficient viral load to ensure accurate interpretation. For low viral load samples (e.g. Ct > 35), accurate detection may be compromised due to the constraints required by the relative Ct comparison as well as the assay's limit of detection. This caveat, however, is less stringent than the viral load constraints needed to successfully sequence SARS-CoV-2 samples, which in practice require samples with Ct ≤ 30.

We successfully monitored the dynamics of lineage B.1.1.7 prevalence by retesting nearly 7,000 SARS-CoV-2 positive samples in three successive testing campaigns in Slovakia. This mass surveillance effort provided invaluable information about the spread and prevalence of lineage B.1.1.7 without having to conduct expensive and time-consuming genomic sequencing. Although our observed data did not fit models of previously established replication numbers for lineage B.1.1.7 6,8 , we attribute this discrepancy to region specific presence of other competing SARS-CoV-2 variants. Indeed, a recently described B.1.258 variant that was extensively circulating throughout the Czech Republic and Slovakia 21 contains mutations in the spike protein (N439K and ΔH69/ΔV70 deletion) that result in higher viral loads and increased transmissibility 7, 21, 40 . It is plausible that the reproduction number of B.1.1.7 is altered in a time and region-dependent manner that is associated with circulation of B. Although the assays described here can accurately detect and differentiate the B.1.1.7 variant, a limitation is that it does not detect other emerging variants of concern (VOCs) such as B.1.351 (Beta), P.1 (Gamma), and B.1.617.2 (Delta). This is particularly important in geographical regions where B.1.1.7 is not prevalent or is being superseded by other VOCs. An important future objective will be to design additional assays that could be multiplexed or used in parallel with our B.1.1.7 assay to detect other VOCs (or mutations of concern). These panels could be tailored to specific geographical regions to provide rapid snapshots of VOC prevalence and their dynamics over time in the absence of sequencing. Alternatively, these panels could be used to quickly identify putative variant samples that can be marked for targeted sequencing. As an example of this approach, we recently developed an additional assay targeting the ORF1a (nsp6) ΔS3675/ΔG3676/ΔF3677 deletion that is shared by B. 1.1.7 40 .

Although the main use of this assay is to directly detect the B.1.1.7 variant, the multiplexed SARS-CoV-2 S gene and RNase P assay could also be used as an affordable and ultrasensitive diagnostic test. Indeed, with an ultrasensitive LoD of only 0.4 copies/μl, this assay is among the most sensitive open-source and commercial tests with EUA approval 41 . This notion is in accordance with several independent comparisons of open-source and commercially available SARS-CoV-2 tests that show our S gene/RNase P assay is a highly sensitive test [42] [43] [44] [45] [46] .

Given that most open-source protocols either require multiple singleplex reactions and/or lack quality control assays, our multiplexed S gene/RNase P assay can provide a streamlined, affordable, and ultrasensitive alternative to currently available open-source and commercial SARS-CoV-2 diagnostic tests.

We have provided interested users with the primer and probe sequences to implement this B.1.1.7 assay in their own laboratories with the hope this can rapidly scale the ability of countries to identify the B.1.1.7 variant and implement epidemiological measures to mitigate its spread. This test can provide labs with a powerful tool to directly confirm the presence of the B.1.1.7 variant in a sample previously determined SARS-CoV-2 positive by an approved screening test, thus avoiding the use of target gene failure assays that can be plagued with low specificity and obviating the need to conduct burdensome and costly genomic sequencing. This is particularly important for countries that are experiencing extensive circulation of variants harbouring only the ΔH69/ΔV70 deletion as current RT-qPCR assays that rely on SGTFs erroneously classify these samples as presumptive B.1.1.7 variants.

Identification of RT-qPCR targets by bioinformatic analysis. To identify suitable targets for primer/ probe design, we downloaded 1,136 sequences from the GISAID repository filtered during a collection time spanning 1-21 December 2020. We focused on the spike gene because lineage B.1.1.7 contains a number of spike gene mutations, including two deletions (ΔH69/ΔV70 and ΔY144) that are ideal for designing a specific assay. We cut the locus encoding the spike protein and used the MAFFT alignment tool (with the parameterauto) 47 to align all the sequences against the WUHAN reference (NCBI ID: NC_045512.2). Twelve sequences (1.06%) contained ambiguous signal in the loci of deletions and were not used in the downstream analysis. We separated sequences into two groups: (1) those with the ΔH69/ΔV70 and ΔY144 deletions and (2) those without the deletions (Supplementary Table S1 ). Using SeaView 48 , we called 95% consensus sequences for the ΔH69/ ΔV70 and ΔY144 group and the No deletions group that were subsequently used to design primer and probe sets specific to either B.1.1.7 or all other SARS-CoV-2 variants, respectively.

In a separate analysis to determine the prevalence of the ΔH69/ΔV70 and ΔY144 deletions in lineages other than B.1.1.7, we used the multiple sequence alignment (MSA) dataset from GISAID and the metadata file containing 922,292 sequences (of which 922,190 were determined to be unique) collected from the beginning of the pandemic through 6 April 2021. Using custom python scripts, we searched at the nucleotide level for loci with both ΔH69/ΔV70 and ΔY144 deletions and sorted hits according to Pangolin lineage ID. All commands and scripts are available here: https:// github. com/ Multi plexDX/ B117-RT-qPCR-design/ blob/ main/ B1117. ipynb.

Primer design and synthesis. We designed primers and probes using the 95% consensus sequences to target the S gene of the common SARS-CoV-2 (called SARS-CoV-2 S gene assay). To differentiate the B.1.1.7 variant from all other SARS-CoV-2 variants, we also designed primers and probes to target the S gene of SARS-CoV-2 variants containing either the ΔH69/ΔV70 deletion or the ΔY144 deletion, or both deletions (called B.1.1.7 assay). As an internal control, we synthesized a primer/probe set for human RNase P published by the US CDC 17 . We incorporated locked nucleic acid (LNA)-modified bases into some primers and probes following general guidelines in order to normalize melting temperatures, increase sensitivity, and enhance specificity 49 www.nature.com/scientificreports/ Following primer/probe design, we conducted in silico analyses using the IDT OligoAnalyzer™ tool (https:// www. idtdna. com/ pages/ tools/ oligo analy zer) to verify melting temperature (Tm), GC content, and potential to form homo-/hetero-dimers as well as the mFold server 52 (http:// www. bioin fo. rpi. edu/ appli catio ns/ mfold/) to identify problematic secondary structures or necessary hairpin formation for TaqMan probes. Probes for both SARS-CoV-2 S gene and B.1.1.7 were labelled with a 5'-FAM (6-carboxyfluorescein) reporter dye and 3'-BHQ-1 (black hole quencher 1) dye. In multiplexed reactions, the probe for human RNase P was labelled with a 5'-Cy5 (cyanine 5) reporter dye and 3'-BHQ-3 dye. Primers and probes were synthesized at MultiplexDX, s.r.o. (Bratislava, Slovakia; https:// www. multi plexdx. com/). The sequences of primers and probes used in this study are listed in Supplementary Table S2 . This test is also available as an No. COV000, Exact Diagnostics, TX, USA) containing genomic DNA at a concentration of 75,000 copies/ml, resulting in samples with concentrations of 8 copies/μl (= 40 copies/reaction), 2 copies/μl (= 10 copies/reaction), 0.8 copies/μl (= 4 copies/reaction), 0.4 copies/μl (= 2 copies/reaction) and 0.2 copies/μl (= 1 copy/reaction) that were used in the analytical sensitivity test. The assay was performed in 8 replicates for each prepared dilution.

Clinical performance evaluation. We evaluated the clinical utility of our SARS-CoV-2 S gene and B.1.1.7

primer/probe sets using a selected set of 106 clinical samples, which were collected during December 2020 and January 2021 and confirmed by an RT-qPCR reference method used for routine testing by regional public health authorities of the Slovak Republic. Further sequencing revealed 67 of these samples belonging to the B. To ensure our clinical validation was sufficiently powered given the sample size, we conducted a post-hoc sample size estimate using a previously described approach for assessing diagnostic sensitivity and specificity 54 . The calculation contained the following parameters: expected sensitivity = 0.96; expected specificity = 0.95; prevalence of disease (p) = 0.25; precision (± expected) = 0.10; confidence level 100(1 − α) = 95%; expected drop out rate = 5%. This sample size estimate calculated a sample size for sensitivity, n sen = 64 (with 5% dropout) and a sample size www.nature.com/scientificreports/ for specifity, n spec = 25, or a total of 89 clinical samples. Therefore, our clinical evaluation of 106 SARS-CoV-2 positive samples, of which 67 samples were classified as B.1.1.7 and 39 samples classified as other lineages, and 94 SARS-CoV-2 negative samples was sufficiently powered. We conducted a detailed assessment of potential biases and applicability judgements of the clinical validation using the QUADAS-2 tool 55 (http:// www. brist ol. ac. uk/ popul ation-health-scien ces/ proje cts/ quadas/ quadas-2/; Supplementary Table S3 ). Inclusion criteria for the selected set of samples consisted of a positive result from a reference standard RT-qPCR test as well as available sequencing information to confirm lineage.

The clinical validation was conducted by the Biomedical Research Center, Institute of Virology, Slovak Academy of Sciences (BMC SAS). BMC SAS received clinical samples previously identified as SARS-CoV-2 positive by local laboratories, then extracted RNA and performed the reference standard test (rTEST COVID-19 qPCR Allplex kit; MultiplexDX) and the B.1.1.7 test in parallel. All samples were processed and tested in a timely manner to minimize the effects of RNA degradation. For the clinical evaluation, researchers used a prespecified criterion to interpret test results (Supplementary Table S4 ) and were blind to the sequencing outcome. The results from the B.1.1.7 assay were compared with sequencing outcome to determine the clinical sensitivity, specificity, and positive and negative likelihood ratios, which were calculated using MedCalc's free online diagnostic test evaluation calculator (for definitions of formulas see, https:// www. medca lc. org/ calc/ diagn ostic_ test. php) 56 . Estimates of the prevalence of B.1.1.7 and B.1.258 were derived from sequencing data available in Slovakia during December 2020 and January 2021 and were used to calculate positive and negative predictive values and accuracy. The results of the index test were not made available to the experimenters evaluating the reference standard or sequencing results. Reporting of the clinical evaluation of our B.1.1.7 RT-qPCR test is in accordance with STARD reporting guidelines (Supplementary Table S5 ).

All experiments were reviewed and approved by the Ethics committee of the Biomedical research Center of the Slovak Academy of Sciences, Bratislava, Slovakia (Ethics committee statement No. EK/BmV-02/2020) and were performed according to their regulations and guidelines. Informed consent was obtained from all subjects.

Surveillance of lineage B.1.1.7 prevalance throughout the Slovak Republic. This test (rTEST COVID-19 B.1.1.7 qPCR kit) was used to assess the prevalence of B.1.1.7 throughout the Slovak Republic over a period of 1 month. Multiple laboratories of the Public Health Authority of the Slovak Republic retested 6,886 samples that were identified as SARS-CoV-2 positive by a standard RT-qPCR test. Retesting occurred over a one month period on February 2nd, 2021 (1,962 samples), February 17th, 2021 (2,382 samples), and March 3rd, 2021 (2,542 samples). To ensure sufficient sample size for smaller geographic regions, we grouped districts based on the jurisdictions of the regional public health offices (there are 36 public health offices each covering 1-8 districts) and distributed tests accordingly. Laboratory personnel determined lineage status using the predefined criteria outlined in Supplementary Table S4 . For comparison with our observed B.1.1.7 prevalences, we also used the prevalences obtained on the first screening round (February 2nd, 2021) to estimate historical and future prevalences using a range of 2 weeks and applying spread factors (3.3 and 5.4) that were derived from estimated reproduction numbers from the UK 6, 8 . Data from the surveillance screening can be found on GitHub: https:// github. com/ Insti tut-Zdrav otnych-Analyz/ covid 19-data/ tree/ main/ PCR_ Tests.

Raw data from the clinical validation as well as the names and accession ID of sequences uploaded to GISAID are provided in Supplementary Table S7 . All commands and scripts for bioinformatic analyses of GISAID sequences are available here: https:// github. com/ Multi plexDX/ B117-RT-qPCR-design/ blob/ main/ B1117. ipynb. Data from the surveillance screening can be found on GitHub: https:// github. com/ Insti tut-Zdrav otnych-Analyz/ covid 19data/ tree/ main/ PCR_ Tests and in Supplementary Table S8 . 

Preliminary Genomic Characterisation of an Emergent SARS-CoV-2 Lineage in the UK Defined by a Novel Set of Spike Mutations

GISAID, EpiCoV Data Curation Team. GISAID-UK Reports New Variant, Termed VUI 202012/01

Enhanced binding of the N501Y-mutated SARS-CoV-2 spike protein to the human ACE2 receptor: Insights from molecular dynamics simulations

Deep mutational scanning of SARS-CoV-2 receptor binding domain reveals constraints on folding and ACE2 binding

Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England

SARS-CoV-2 variant of concern 202012/01 has about twofold replicative advantage and acquires concerning mutations

Recurrent emergence of SARS-CoV-2 spike deletion H69/V70 and its role in the Alpha variant B.1.1.7. Cell Rep

Transmission of SARS-CoV-2 lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data

Early analysis of a potential link between viral load and the N501Y mutation in the SARS-COV-2 spike protein

S-variant SARS-CoV-2 lineage B.1.1.7 is associated with significantly higher viral load in samples tested by TaqPath polymerase chain reaction

NERVTAG Paper on COVID-19 Variant of

Use of dual TaqMan probes to increase the sensitivity of 1-step quantitative reverse transcription-PCR: Application to the detection of SARS coronavirus

Evaluation of TaqMan qPCR system integrating two identically labelled hydrolysis probes in single assay

Detection of single DNA base differences by competitive oligonucleotide priming

Determination of a T/G polymorphism at nucleotide 3206 of the apolipoprotein C III gene by amplification refractory mutation system

Single-nucleotide polymorphisms and other mismatches reduce performance of quantitative PCR assays

Centers for Disease Control and Prevention. Research Use Only 2019-Novel Coronavirus (2019-nCoV) Real-Time RT-PCR Primers and Probes

Two-step strategy for the identification of SARS-CoV-2 variant of concern 202012/01 and other variants with spike deletion H69-V70, France

Identification of a Novel SARS-CoV-2 Spike 69-70 Deletion Lineage Circulating in the United States-SARS-CoV-2 coronavirus/SARS-CoV-2 Molecular Evolution

Detection of Non-B.1.1.7 Spike ∆69/70 Sequences (B.1.375) in the United States-SARS-CoV-2 coronavirus/nCoV-2019 Genomic Epidemiology

A SARS-CoV-2 mutant from B.1.258 lineage with ∆H69/∆V70 deletion in the Spike protein circulating in Central Europe in the fall 2020

Early introductions and transmission of SARS-CoV-2 variant B117 in the United States

Early assessment of diffusion and possible expansion of SARS-CoV-2 Lineage 20I/501Y.V1 (B.1.1.7, variant of concern 202012/01) in France

S gene dropout patterns in SARS-CoV-2 tests suggest spread of the H69del/V70del mutation in the US

B.1.1.7 became the dominant variant in Lebanon

Investigation of novel SARS-CoV-2 variant: Variant of concern 202012/01. Technical briefing 4

Quantification of the spread of SARS-CoV-2 variant B117 in Switzerland

Rapid SARS-CoV-2 variants spread detected in France using specific RT-PCR testing

Fast and cost-effective screening for SARS-CoV-2 variants in a routine diagnostic setting

SARS-CoV-2 N501Y introductions and transmissions in Switzerland from beginning of

40 minutes RT-qPCR assay for screening spike N501Y and HV69-70del mutations

Rapid detection of SARS-CoV-2 variants of concern identifying a cluster of B.1.1.28/P1 variant in British Columbia

A rapid and low-cost protocol for the detection of B.1.1.7 lineage of SARS-CoV-2 by using SYBR green-based RT-qPCR

A simple RT-PCR Melting temperature assay to rapidly screen for widely circulating SARS-CoV-2 variants

Evaluation of a fully automated high-throughput SARS-CoV-2 multiplex qPCR assay with built-in screening functionality for del-HV69/70-and N501Y variants such as B.1.1.7

Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2

RT-qPCR assay for detection of British (B.1.1.7) and South Africa (B.1.351) variants of SARS-CoV-2

Monitoring the propagation of SARS CoV2 variants by tracking identified mutation in wastewater using specific RT-qPCR

Multiplexing primer/probe sets for detection of SARS-CoV-2 by qRT-PCR

Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity

The COVID-19 XPRIZE and the need for scalable, fast, and widespread testing

AccuPower SARS-CoV-2 Multiplex RT-PCR kit (Bioneer, South Korea)" and "Allplex 2019-nCoV assay (Seegene, South Korea)" for SARS-CoV-2 RT-PCR diagnosis: Korean CDC EUA as a quality control proxy for developing countries

Clinical performance and analytical sensitivity of three SARS-CoV-2 nucleic acid diagnostic tests

Comparative performance of SARS-CoV-2 detection assays using seven different primer-probe sets and one assay kit

Limits of detection of 6 approved RT-PCR kits for the novel SARS-coronavirus-2 (SARS-CoV-2)

Analytical sensitivity and efficiency comparisons of SARS-CoV-2 RT-qPCR primer-probe sets

Parallelization of MAFFT for large-scale multiple sequence alignments

SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building

Sequence-dependent thermodynamic parameters for locked nucleic acid (LNA)−DNA duplex formation

Design considerations and effects of LNA in PCR primers

Position-dependent effects of locked nucleic acid (LNA) on DNA sequencing and PCR primers

Mfold web server for nucleic acid folding and hybridization prediction

SARS-CoV-2 genomes recovered by long amplicon tiling multiplex approach using nanopore sequencing and applicable to other sequencing platforms

Statistical methodology: I. Incorporating the prevalence of disease into the sample size calculation for sensitivity and specificity

QUADAS-2: A revised tool for the quality assessment of diagnostic accuracy studies

Diagnostic Test Evaluation Calculator

performed genomic sequencing and bioinformatics to determine sample lineage status

We gratefully acknowledge the authors from the submitting and originating laboratories that shared genetic sequencing data with the GISAID initiative. We acknowledge the contribution of the oligonucleotide synthesis and production team at MultiplexDX for synthesizing the primers and probes in this study, the laboratory staff of the Slovak regional public health authorities and diagnostic laboratories for screening samples during surveillance of B.1.1.7, and the team at the Institute for Healthcare Analyses for processing the surveillance data. 

VK, EDP, MR, RH, and PC are employees of MultiplexDX, a biotechnology company which has commercialized a kit called rTEST COVID-19 B.1.1.7 qPCR kit (https:// www. multi plexdx. com/ produ cts/ rtest-covid-19-b-1-1-7qpcr-kit, MultiplexDX, Inc., Bratislava, Slovakia) that is based on this research. BK is a Head of the Department of Virus Ecology, Institute of Virology, Biomedical Research Center of the Slovak Academy of Sciences (BMC SAS). MM is a Head of the Insititute for healthcare analyses at the Slovak Ministy of Health. BMC SAS has entered into collaboration with MultiplexDX, Inc. for development and validation of RT-qPCR tests for routine detection of SARS-CoV-2 and the test for detection of B.1.1.7 variant described in this study. The Ministry of Health procured said tests for the B.1.1.7 lineage surveillance. All other authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-99661-7.Correspondence and requests for materials should be addressed to P.Č. or B.K.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.