NATIONAL SCALE REAL-TIME SURVEILLANCE OF SARS-COV-2 VARIANTS DYNAMICS BY WASTEWATER MONITORING IN ISRAEL

In this report, we describe a national-scale monitoring of the SARS-COV-2 (SC-2) variant dynamics in Israel, using multiple-time sampling of twelve wastewater treatment plants. We used a combination of inclusive and selective quantitative PCR assays that specifically identify variants A19 or B.1.1.7 and tested each sample for the presence and relative viral RNA load of each variant. We show that between December-2020 and March-2021, a complete shift in the SC-2 variant circulation was observed, where the B.1.1.7 replaced the A19 in all examined test points. We further show that the normalized viral load (NVL) values and the average new cases per week reached a peak in January 2021, and then decreased gradually in almost all test points, in parallel with the progression of the national vaccination campaign, during February-March 2021. This study demonstrates the importance of monitoring SC-2 variant dynamics on a national scale through wastewater sampling. It also provides a proof-of-concept methodology for continuous surveillance by using a combination of inclusive and selective PCR tests, which is far more amendable for high throughput monitoring compared with sequencing. This approach may be useful for real-time dynamics surveillance of current and future variants, such as the Omicron (BA.1) variant. SynopsisThis study describes the continuous monitoring of the SARS CoV-2 variant B.1.1.7 circulation in wastewater in Israel using a positive/negative quantitative PCR assay.

SPECIFIC DETECTION OF SARS-COV-2 VARIANTS B.1.1.7 (ALPHA) AND B.1.617.2 (DELTA) USING A ONE-STEP QUANTITATIVE PCR ASSAY

In this report, we describe the development of an RT-qPCR assay, termed Alpha Delta assay, which can detect SARS-COV-2 (SC-2) and distinguish between the Alpha (B.1.1.7) and Delta (B.1.617.2) variants. The Alpha- and Delta-specific reactions in the assay target mutations that are strongly linked to the target variant. The Alpha reaction targets the D3L substitution in N gene, and the Delta reaction targets the spike gene 156-158 mutations. Additionally, we developed a second Delta-specific assay, used as a confirmatory test for the Alpha Delta assay that targets the 119-120 deletion in the Orf8 gene. Both reactions have similar sensitivities of 15-25 copies per reaction, similar to the sensitivity of commercial SC-2 detection tests. The Alpha Delta assay and the Orf8-119del assay were successfully used to classify clinical samples that were subsequently analyzed by whole genome sequencing. Lastly, we show that the Alpha Delta and Orf8-119del assays correctly identified the presence of Alpha and Delta lineages RNA in wastewater samples. This study provides a rapid, sensitive and cost-effective tool for detecting and classifying two worldwide dominant SC-2 variants. It also highlights the importance of a timely diagnostic response to the emergence of new SC-2 variants with significant consequences on global health.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing provides critical information to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments and vaccines, and guide public health responses. To investigate the spatiotemporal heterogeneity in the global SARS-CoV-2 genomic surveillance, we estimated the impact of sequencing intensity and turnaround times (TAT) on variant detection in 167 countries. Most countries submit genomes >21 days after sample collection, and 77% of low and middle income countries sequenced <0.5% of their cases. We found that sequencing at least 0.5% of the cases, with a TAT <21 days, could be a benchmark for SARS-CoV-2 genomic surveillance efforts. Socioeconomic inequalities substantially impact our ability to quickly detect SARS-CoV-2 variants, and undermine the global pandemic preparedness. One-Sentence SummarySocioeconomic inequalities impacted the SARS-CoV-2 genomic surveillance, and undermined the global pandemic preparedness.

The unique evolutionary dynamics of the SARS-CoV-2 Delta variant

The SARS-Coronavirus-2 (SARS-CoV-2) driven pandemic was first recognized in late 2019, and the first few months of its evolution were relatively clock-like, dominated mostly by neutral substitutions. In contrast, the second year of the pandemic was punctuated by the emergence of several variants that bore evidence of dramatic evolution. Here, we compare and contrast evolutionary patterns of various variants, with a focus on the recent Delta variant. Most variants are characterized by long branches leading to their emergence, with an excess of non-synonymous substitutions occurring particularly in the Spike and Nucleocapsid proteins. In contrast, the Delta variant that is now becoming globally dominant, lacks the signature long branch, and is characterized by a step-wise evolutionary process that is ongoing. Contrary to the "star-like" topologies of other variants, we note the formation of several distinct clades within Delta that we denote as clades A-E. We find that sequences from the Delta D clade are dramatically increasing in frequency across different regions of the globe. Delta D is characterized by an excess of non-synonymous mutations, mostly occurring in ORF1a/b, some of which occurred in parallel in other notable variants. We conclude that the Delta surge these days is composed almost exclusively of Delta D, and discuss whether selection or random genetic drift has driven the emergence of Delta D.

BNT162b2 Vaccination efficacy is marginally affected by the SARS-CoV-2 B.1.351 variant in fully vaccinated individuals.

BackgroundIsraeli has vaccinated over 80% of its adult population, with two doses of the Pfizer BNT162b2 vaccine. This intervention has been highly successful in curtailing the coronavirus 2 outbreak. One major concern is the ability of the virus to mutate which potentially can cause SARS-CoV-2 to partially escape from the immune system. Here we evaluate the efficacy of the Pfizer vaccine against the B.1.351 variant. MethodsThe Ministry of Health, initiated sequencing of selected positive swab samples identified as being of interest. We used logistic regression, with variant type as the dependent variable, vaccination status as the main explanatory variable, controlling for age, sex, subpopulation, place of residence and time of sample, to estimate the odds ratio for a vaccinated case to have the B. 1.351 versus the B.1.1.7 variant, within vaccinated and unvaccinated persons who tested positive. FindingsThere were 19 cases of B.1.351 variant (3.2%) among those vaccinated more than 14 days before the positive sample and 88 (3.5%) among the unvaccinated. The estimated odds ratio was 1.29 [95% CI: 0.66-2.50]. From this result, assuming the efficacy against the B.1.1.7 variant to be 95%, the estimated efficacy against the B.1.351 variant was 94% [95% CI: 87-97%]. InterpretationDespite the concerns caused by the B.1.351 variant, the BNT162b2 vaccine seems to provide substantial immunity against both that variant and the B.1.1.7. Our results suggest that from 14 days following the second vaccine dose the efficacy of BNT162b2 vaccine is at most marginally affected by the B.1.351 variant. FundingNo funding

The rise and fall of an emerging SARS-CoV-2 variant with the spike protein mutation L452R

Emerging SARS-CoV-2 variants may threaten global vaccination efforts and awaited reduction in outbreak burden. In this study, we report a novel variant carrying the L452R mutation that emerged from a local B.1.362 lineage, B.1.362+L452R. The L452R mutation is associated with the Delta and Epsilon variants and was shown to cause increased infection and reduction in neutralization in pseudoviruses. Indeed, the B.1.362+L452R variant demonstrated a X4-fold reduction in neutralization capacity of sera from BNT162b2-vaccinated individuals compared to a wild-type strain. The variant infected 270 individuals in Israel between December 2020 and March 2021, until diminishing due to the gain in dominance of the Alpha variant in February 2021. This study demonstrates an independent, local emergence of a variant carrying a critical mutation, L452R, which may have the potential of becoming a variant of concern and emphasizes the importance of routine surveillance and detection of novel variants among efforts undertaken to prevent further disease spread.

Rapid And high throughput RT-qPCR assay for identification and differentiation between SARS-CoV-2 variants B.1.1.7 and B.1.351

Emerging SARS-CoV-2 (SC-2) variants with increased infectivity and vaccine resistance are of major concern. Rapid identification of such variants is important for the public health activities and provide valuable data for epidemiological and policy decision making. We developed a multiplex quantitative RT-qPCR (qPCR) assay that can specifically identify and differentiate between the emerging B.1.1.7 and B.1.351 SC-2 variants. In a single assay, we combined four reactions: one that detects SC-2 RNA independently of the strain, one that detects the D3L mutation, which is specific to variant B.1.1.7, and one that detects the 242-244 deletion, which is specific to variant B.1.351. The fourth reaction identifies human RNAseP gene, serving as an endogenous control for RNA extraction integrity. We show that the strain-specific reactions target mutations that are strongly associated with the target variants, and not with other major known variants. The assays specificity was tested against a panel of respiratory pathogens (n=16), showing high specificity towards SC-2 RNA. The assays sensitivity was assessed using both In-vitro transcribed RNA and clinical samples, and was determined to be between 20 and 40 viral RNA copies per reaction. The assay performance was corroborated with Sanger and whole genome sequencing, showing complete agreement with the sequencing results. The new assay is currently implemented in the routine diagnostic work at the Central Virology Laboratory, and may be used in other laboratories to facilitate the diagnosis of these major worldwide circulating SC-2 variants.

A unique SARS-CoV-2 spike protein P681H strain detected in Israel

Routine detection, surveillance and reporting of SARS-CoV-2 novel variants is important, as these threaten to hinder vaccination efforts. Herein we report a local novel strain that includes a non-synonymous mutation in the spike (S) protein - P681H and additional synonymous mutations. The P681H Israeli strain has not been associated with higher infection rates and was neutralized by sera from vaccinated individuals in comparable levels to the B.1.1.7 strain and a non-P681H strain from Israel.

HiSpike: A high-throughput cost effective sequencing method for the SARS-CoV-2 spike gene

The changing nature of the corona virus of the SARS-CoV-2 pandemic poses unprecedented challenges to the worlds health systems. New and virulent emerging spike gene variants, such as the UK 20I/501Y.V1 and South African 20H/501Y.V2, could jeopardize global efforts to produce immunity and reduce mortality. These challenges require effective real-time genomic surveillance solutions that the medical community can quickly adopt. The SARS-CoV-2 spike protein mediates host receptor recognition and entry into the cell and therefore, it is most susceptible to generation of variants with increased transmissibility and pathogenicity. The spike protein is also the primary target of neutralizing antibodies in COVID-19 patients and the most common antigen for induction of effective vaccine immunity. Therefore, tight monitoring of the spike protein gene variants is key to mitigating COVID-19 spread and vaccine escape mutants. Currently, the ARTIC method for SARS-CoV-2 whole genome sequencing is applied worldwide. However, this method commonly requires more than 96 hours (4-5 days) from start to finish and at present high sample sequence demands, sequencing resources are quickly exhausted. In this work, we present HiSpike, a method for high-throughput targeted next generation sequencing (NGS) of the spike gene. This simple three-step method can be completed in less than 30 hours and can sequence 10-fold more samples compared to the conventional ARTIC method and at a fraction of the cost. HiSpike was proven valid, and has identified, at high quality, multiple spike variants from real-time field samples, such as the UK and the South African variants. This method will certainly be effective in discovering future spike mutations. Therefore, running HiSpike for full sequencing of the spike gene of all positive SARS-CoV-2 samples could be considered for near real-time detection of known and emerging spike mutations as they evolve. HiSpike provides affordable sequencing options to help laboratories conserve resources, hence it provides a tool for widespread monitoring, that can support critical knowledge-based decisions.

Improved Sensitivity, Safety, and Rapidity of COVID-19 Tests by Replacing Viral Storage Solution with Lysis Buffer

Conducting numerous, rapid, and reliable PCR tests for SARS-CoV-2 is essential for our ability to monitor and control the current COVID-19 pandemic. Here, we tested the sensitivity and efficiency of SARS-CoV-2 detection in clinical samples collected directly into a mix of lysis buffer and RNA preservative, thus inactivating the virus immediately after sampling. We tested 79 COVID-19 patients and 20 healthy controls. We collected two samples (nasopharyngeal swabs) from each participant: one swab was inserted into a test tube with Viral Transport Medium (VTM), following the standard guideline used as the recommended method for sample collection; the other swab was inserted into a lysis buffer supplemented with nucleic acid stabilization mix (coined NSLB). We found that RT-qPCR tests of patients were significantly more sensitive with NSLB sampling, reaching detection threshold 2.1{+/-}0.6 (Mean{+/-}SE) PCR cycles earlier then VTM samples from the same patient. We show that this improvement is most likely since NSLB samples are not diluted in lysis buffer before RNA extraction. Re-extracting RNA from NSLB samples after 72 hours at room temperature did not affect the sensitivity of detection, demonstrating that NSLB allows for long periods of sample preservation without special cooling equipment. We also show that swirling the swab in NSLB and discarding it did not reduce sensitivity compared to retaining the swab in the tube, thus allowing improved automation of COVID-19 tests. Overall, we show that using NSLB instead of VTM can improve the sensitivity, safety, and rapidity of COVID-19 tests at a time most needed.