key: cord-0740165-mtme3nyx authors: Abdel Sater, Fadil; Younes, Mahmoud; Nassar, Hassan; Nguewa, Paul; Hamze, Kassem title: 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 date: 2021-10-06 journal: Mol Biol Rep DOI: 10.1007/s11033-021-06717-y sha: 78344c0b40d6bc74532c7d8e41997a68c60905fe doc_id: 740165 cord_uid: mtme3nyx BACKGROUND: The new SARS-CoV-2 variant VOC (202012/01), identified recently in the United Kingdom (UK), exhibits a higher transmissibility rate compared to other variants, and a reproductive number 0.4 higher. In the UK, scientists were able to identify the increase of this new variant through the rise of false negative results for the spike (S) target using a three-target RT-PCR assay (TaqPath kit). METHODS: To control and study the current coronavirus pandemic, it is important to develop a rapid and low-cost molecular test to identify the aforementioned variant. In this work, we designed primer sets specific to the VOC (202012/01) to be used by SYBR Green-based RT-PCR. These primers were specifically designed to confirm the deletion mutations Δ69/Δ70 in the spike and the Δ106/Δ107/Δ108 in the NSP6 gene. We studied 20 samples from positive patients, detected by using the Applied Biosystems TaqPath RT-PCR COVID-19 kit (Thermo Fisher Scientific, Waltham, USA) that included the ORF1ab, S, and N gene targets. 16 samples displayed an S-negative profile (negative for S target and positive for N and ORF1ab targets) and four samples with S, N and ORF1ab positive profile. RESULTS: Our results emphasized that all S-negative samples harbored the mutations Δ69/Δ70 and Δ106/Δ107/Δ108. This protocol could be used as a second test to confirm the diagnosis in patients who were already positive to COVID-19 but showed false negative results for S-gene. CONCLUSIONS: This technique may allow to identify patients carrying the VOC (202012/01) or a closely related variant, in case of shortage in sequencing. As the pandemic of SARS-CoV-2 continues to affect the planet, researchers around the world are monitoring the virus and detecting acquired mutations that may lead to higher threats of spreading COVID-19 [1] . A new variant was recently reported in the UK as a VOC 202012/01 or B.1.1.7 lineage [2] . Its rate of transmission is estimated to be >70%, and its reproductive number (Ro) seems to be up to 0.4 higher [3] . This variant harbors 14 non-synonymous mutations, 6 synonymous mutations and 3 deletions [2, 3] . A wide variety of diagnostic tests have been used by high-throughput national testing systems around the world, to monitor the SARS-CoV-2 infection [4] . The arising prevalence of new SARS-CoV-2 variants such as B.1.1.7 has become of great concern, as most of the existing RT-PCR tests will not be able to specifically distinguish these new variants because they were not designed for such a purpose. Therefore, public health officials most rely on their current testing systems and their sequencing results to draw conclusions on the prevalence of new variants in their territories [2, 5] . An example of such cases has been seen in the UK. In fact, scientists were able to identify the augmentation of the B.1.1.7 SARS-CoV-2 variant infection in the population through an increase in the S-gene target failure in their three target gene assay (N+, ORF1ab+, S−) when using the Applied Biosystems TaqPath RT-PCR COVID-19 kit (Thermo Fisher Scientific, Waltham, USA) that included the ORF1ab, S, and N gene targets [1, 6, 7] . In 24 December 2020, Thermo Fisher Scientific confirmed that the S deletion Δ69/Δ70 was in the area targeted by the TaqPath Kit. Whilst other variants with Δ69/70 are also circulating worldwide, the absence of detection of the S gene target increasingly appears to be a highly specific marker for the new variant [2, 7] . In addition, the European CDC recommended that multi-target RT-PCR assays that included an S gene target affected by the deletions could be used as a signal for the presence of the Δ69/Δ70 mutation for further investigations and could be helpful to keep tracking these mutant strains [8] . Genome sequencing is the gold method to confirm the new variant, but observational studies provide also stronger evidence if similar models are observed in multiple countries, especially when randomized studies are not possible. In Lebanon, surveillance data from Beirut Medical Center and Bahman Hospital showed a rapid and dramatic increase in S-negative profile in PCR testing for SARS-CoV-2 in the first twelve days in January, reaching approximately 60% of positive cases and 95% in February [9, 10] . In 15 January 2021, in GISAID we found 374,613 SARS-CoV-2 sequences, among them 17,473 sequences belonged to the new variant (VOC 202012/01). It is important to notice that the co-occurrence of mutations Δ69/Δ70 in the spike and Δ106/Δ107/Δ108 in the NSP6, was present in this new variant. We then decided to use these two deletions as targets for our rapid and low-cost protocol, performing SYBR Green-Based RT-PCR. In this work, we propose primer sets that can be applied as a second step to confirm the diagnosis in cases that were already detected as positive to SARS-CoV-2 when using the Applied Biosystems TaqPath RT-PCR COVID-19 kit (Thermo Fisher Scientific, Waltham, USA). It is well known that confirmatory diagnosis based on specific diagnostic biomarkers remains a great challenge allowing the further control and eradication of infections including COVID-19 [11] . Therefore, this novel method may lead to specifically identify individuals carrying the Δ69/ Δ70 and the Δ106/Δ107/Δ108 mutations. 20 clinical samples for SARS-CoV-2 positive patients, with Ct (cycle threshold) value < 32, were selected for this study. All samples were previously tested positive by the Applied Biosystems™ TaqPath™ COVID-19 assay which targeted the RdRP, N and Spike genes in the Molecular Laboratories at Bahman Hospital and Beirut Cardiac Institute in the city of Beirut. 16 of these samples were S-negative and 4 were S-positive. All patients provided written and signed informed consents. RNA was extracted from the clinical samples using Qiamp viral RNA mini kit (Qiagen). Total RNA was converted to cDNA using iScript cDNA synthesis kit (BioRad), following the manufacturer's recommended procedures. Primers were designed based on the recently available full sequence of the new variant VUI (202012/01) from the NCBI Reference Sequence Database (https:// www. ncbi. nlm. nih. gov/ nucco re/ NC_ 045512). We utilized NCBI-Primer BLAST (https:// www. ncbi. nlm. nih. gov/ tools/ primer-blast/) to design specific primers (Fig. 1 ). The specificity of the primers was verified using in silico prediction analyses with the online Basic Local Alignment Search Tool (BLAST) (https:// blast. ncbi. nlm. nih. gov/ Blast. cgi). None of our designed primers showed genomic cross-reactivity with other viruses, the human genome, or other probable interfering genomes in the BLAST database analysis. The primer sets were synthesized and delivered by Macrogen (Republic of Korea), they are listed in Table 1 . Genomic sequencing was performed for the 20 samples, collected between 03/01/2021 and 04/04/2021. Samples were sequenced using the ARTIC network methodology (https:// artic. netwo rk/) with the V3 amplicon scheme, V3-LoCost library prep method (Quick 2020) and consensus sequences were generated using the bioinformatics SOP (https:// artic. netwo rk/ ncov-2019/ ncov2 019-bioin forma tics-sop. html) in Nano polish mode. 8,9 Consensus sequences were assigned a lineage using pangolin v2.0 (github.com/cov-lineages/ pangolin). Quantitative RT-PCR was carried out using iTaq universal SYBR Green super mix (Bio Rad). In brief, each reaction consisted of a total volume of 20 µL containing 2 µL of each primer [10 pM/µL], 5 µL of cDNA, 10 µL SYBR Green super mix and 1 µL of Rnase free Water. Real-time PCR was performed using Bio Rad CFX96 Real-Time PCR Machine. The thermal cycling conditions used were as follows: 94 °C for 2 min, followed by 40 cycles of amplification at 94 °C for 10 s, and 60 °C for 1 min. The reaction was completed by determining the dissociation curve of all amplicons. In order to validate our method, we re-tested 20 samples that had been tested as positive for SARS-CoV-2 with Ct < 32, by the TaqPath kit. 16 samples with S-negative profile and four samples with S-positive profile ( Table 2) . N primer pairs were used as a control of cDNA synthesis and integrity of the sample, spike WT and NSP6 WT primer pairs were used to detect the variants not harboring the deletion mutations Δ69/Δ70 in the spike nor the Δ106/ Δ107/Δ108 in the NSP6, respectively. Spike del 69/70 primer pairs were designed to detect the deletions Δ69/ Δ70 in the spike protein (Fig. 1) . Our results showed that the N-gene amplicons were detected in both, S-negative and S-positive profiles ( Table 2 ; Fig. 2 ), whereas the spike del 69/70-gene amplicons were detected only in the S-negative profile confirming the absence of the amino acids 69 and 70 in the spike protein of these samples (Fig. 2B) . For the spike WT and The results of spike del 69/70 primers were fully concordant with spike WT and NSP6 WT primers, 100% of the S-negative profile had the Spike deletions Δ69/Δ70 (Fig. 3A) , while 100% of the S-positive profile did not contain the deletions Spike Δ69/Δ70 and NSP6 Δ106/ Δ107/Δ108 ( Fig. 2A) . To confirm the lineage present in these 20 patient samples, genomic sequencing was performed and the obtained sequences were deposited in Genbank under NCBI accession number MZ772870 to MZ772873, MZ772924, MZ772926, MZ772929, MZ772933, MZ772935, MZ772973 to MZ772975, MZ773014, MZ773900, MZ773902, MZ773903, MZ773924 and MZ773926 to MZ773928. Genomic sequencing revealed that all the 20 samples belonged to the B.1.1.7 lineage (Table 3) . The analytical sensitivity of the assays was determined with 5 serials tenfold dilutions using an isolated sample having Ct values 25 by TaqMan RT-PCR (TaqPath kit) and 28 by our SYBR Green-based protocol. All diluted samples were tested in triplicate by gold standard TaqMan RT-qPCR (TaqPath kit) and by our SYBR Green-based protocol. The unamplified dilution by SYBR Green-based RT-PCR has Ct value 37.55 by the TaqMan RT-qPCR. Thus, the tested RT-PCRs could detect positive samples with maximal Ct of 34.45 by TaqMan RT-PCR or with Ct value of 37.67 by SYBR Green RT-PCR. All results generated by the serials tenfold dilution are shown in Table 4 . The limit of detection of TaqPath kit is 250 copies/mL (according to the manufacturer). Starting with the same RNA concentration, the Ct values of both TaqPath kit and SYBR Green-based RT-PCR were similar to a difference in Ct of 2.1 on average ( Table 2) . This difference means a decrease of 4.2-folds in the limit of detection for the SYBR Green-based RT-PCR (~ 1050 copies/mL) compared to TaqPath kit. One drawback of the SYBR Green assay is that the dye is nonspecific, and this lack of specificity can generate falsepositive signals if nonspecific products or primer-dimers are present in the sample. However, including a melting curve analysis at the end of each PCR assay to determine the specificity and efficiency of each RT-PCR reaction will ensure the accuracy of the results when multiple peaks or primer-dimers are not observed. Our primer set "spike del 69/70" gave specific melt peak at 79.5 °C C in S-negative samples (Fig. 3B ), while primers "spike" and "NSP8" gave melt peaks at 79 and 79.5 °C, respectively, in S-positive samples (Fig. 2C, D) . Sequencing is the current gold standard diagnostic method to confirm the presence of any new variant, its efficient but a time consuming method and not accessible to all laboratories especially in developing countries based on the elevated costs and the present global demands for supplies and reagents. In this study, we have developed a rapid, low-cost, large-scale screening protocol for the detection of the deletions Δ69/Δ70 and Δ106/Δ107/Δ108. We adopted SYBR Green-based RT-PCR method over Taqman-based RT-PCR method. Taqman-based diagnosis is accepted in the field as more accurate method but also it's not accessible to all laboratories especially in developing countries. In the present work we are dealing with samples already confirmed SARS-CoV-2 positive by TaqPath Kit (ThermoFisher), so even if SYBR Green-based RT-PCR method is less accurate it will be very useful since we need just to confirm the presence of the two mutations Δ69/Δ70 in the spike and Δ106/ Δ107/Δ108 in the NSP6. We hope that our efforts will be helpful and can contribute to the early detection of the new variant (VOC 202012/01) belonging to the B.1.1.7 lineage, for the prevention of transmission and early intervention. This protocol should not be limited to this variant, but also for any other future variant to come, just by designing the appropriate primers. Author contributions KH analysed and interpreted the data and wrote the draft of the manuscript. FAS analysed and interpreted the data and was involved in revising the manuscript critically. MY, HN and PG revised the manuscript critically. All the authors have given the final approval of the version to be published. All authors read and approved the final manuscript. Funding No funding was provided for this work. Data availability We declared that materials described in the manuscript, including all relevant raw data, will be freely available to any Ethical approval and consent to participate All presented cases or their legal guardian provided consent to data collection and use according to institutional guidelines. All presented cases or their legal guardian provided consent to publish according to institutional guidelines. Design of specific primer set for detection of B. 1.1. 7 SARS-CoV-2 variant using deep learning Transmission of SARS-CoV-2 lineage B. 1.1. 7 in England: insights from linking epidemiological and genetic data Investigation of novel SARS-COV-2 variant: variant of concern 202012/01 Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations: COVID-19 genomics UK Consortium Molecular diagnostic technologies for COVID-19: limitations and challenges Two-step strategy for the identification of SARS-CoV-2 variants co-occurring with spike deletion H69-V70, Lyon S gene dropout patterns in SARS-CoV-2 tests suggest spread of the H69del/V70del mutation in the US Abdel Sater F (2021) Emergence and fast spread of B.1.1.7 lineage in Lebanon B. 1.1. 7 became the dominant variant in Lebanon. medRxiv Machine learning, artificial intelligence, and data science breaking into drug design and neglected diseases Acknowledgements This work was supported by the Lebanese University and Fundación La Caixa (LCF/PR/PR13/11080005), Fundación Caja Navarra, Fundación Roviralta, Ubesol, Inversiones Garcilaso de la Vega, COST Actions CA18217 and CA18218, and EU Project uncover (Grant/Award Number: 101016216). The authors declare that they have no competing interests.