Two-step strategy for the identification of SARS-CoV-2 variant of concern 202012/01 and other variants with spike deletion H69-V70, France, August to December 2020 =================================================================================================================================================================== * Antonin Bal * Gregory Destras * Alexandre Gaymard * Karl Stefic * Julien Marlet * Sébastien Eymieux * Hadrien Regue * Quentin Semanas * Constance d’Aubarede * Geneviève Billaud * Frédéric Laurent * Claudia Gonzalez * Yahia Mekki * Martine Valette * Maude Bouscambert * Catherine Gaudy-Graffin * Bruno Lina * Florence Morfin * Laurence Josset * COVID-Diagnosis HCL Study Group ## Abstract We report the implementation of a two-step strategy for the identification of SARS-CoV-2 variants carrying the spike deletion H69-V70 (ΔH69/ΔV70). This spike deletion resulted in a S-gene target failure (SGTF) of a three-target RT-PCR assay (TaqPath kit). Whole genome sequencing performed on 37 samples with SGTF revealed several receptor-binding domain mutations co-occurring with ΔH69/ΔV70. More importantly, this strategy enabled the first detection of the variant of concern 202012/01 in France on December 21th 2020. Since September a SARS-CoV-2 spike (S) deletion H69-V70 (ΔH69/ΔV70) has attracted increasing attention. This deletion was detected in the cluster-5 variant identified both in minks and humans in Denmark. This cluster-5 variant carries a receptor binding domain (RBD) mutation Y453F and was associated with reduced susceptibility to neutralizing antibodies to sera from recovered COVID-19 patients [1–3]. The ΔH69/ΔV70 has also co-occurred with two other RBD mutations of increasing interest [4]: N439K that is currently spreading in Europe and might also have reduced susceptibility to SARS-CoV-2 antibodies [5]; and N501Y that is part of the SARS-CoV-2 variant of concern (VOC) 202012/01 recently detected in England [6]. Although the impact of ΔH69/ΔV70 on SARS-CoV-2 pathogenesis is not clear, enhanced surveillance is urgently needed. Herein we report the implementation of a two-step strategy enabling a rapid detection of VOC 202012/01 or other variants carrying ΔH69/ΔV70. ## ΔH69/ΔV70 associated with S-gene target failure of a three-target RT-PCR assay As part of routine SARS-CoV-2 genomic surveillance performed at national reference centre (NRC) for respiratory viruses (Lyon, France) [7], a 6-nucleotide deletion (21765-21770) within the S gene was identified in two nasopharyngeal samples collected on September 1st and 7th, respectively. The SARS-CoV-2 infection diagnosis had been performed with the Applied Biosystems TaqPath RT-PCR COVID-19 kit (Thermo Fisher Scientific, Waltham, USA) that includes the ORF1ab, S, and N gene targets. For these two samples, a S-gene target failure (SGTF) was reported while ORF1ab and N targets were positive with Ct values < 25 (Figure 1A). ![Figure 1A.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/11/2020.11.10.20228528/F1.medium.gif) [Figure 1A.](http://medrxiv.org/content/early/2021/01/11/2020.11.10.20228528/F1) Figure 1A. Amplification curves obtained with TaqPath COVID-19 RT-PCR kit for samples with the S deletion 21765-21770. The three targets included in the RT-PCR kit are represented by a different color. The amplification curve of the internal control is also represented (MS2, red curve). 1B. Pairwise sequence alignment from nucleotide position 21758 to 21775 of the spike gene using CoV-GLUE resource. Sequence with the deletion 21765-21770 is represented in green and the reference sequence in blue (Wuhan-Hu-1). The 21765-21770 deletion results in deletion of amino acid residues 69 and 70; ATC (21764-21771-21772) encoding for an isoleucine amino acid (I). The mean coverage for the whole genome sequences generated was 6903x and 6898x, respectively and the S deletion 21765-21770 was present in 100% of the reads. Using CoV-GLUE online resource [8], we found that the S deletion 21765-21770 led to the removal of 2 amino acids (ΔH69/ΔV70) in the N-terminal domain of the S1 subunit of the S protein (Figure 1B). The whole genome sequencing (WGS) method used was the amplicon-based ARTIC v3 protocol ([https://artic.network/ncov-2019](https://artic.network/ncov-2019)) combined with Nextera DNA Flex library and sequencing on NextSeq 500 platform (Illumina, San Diego, USA). To confirm the presence of the deletion, one sample was also sequenced with an untargeted metagenomic protocol that yielded the same sequence. Of note, this metagenomic approach could not be applied for the second sample due to low viral load [9]. Although the coordinates of the primer/probe binding regions were not available for the TaqPath kit, the manufacturer confirmed that the S deletion H69-V70 was in the area targeted by the test. ## ΔH69/ΔV70 screening with RT-PCR followed by WGS We then performed a retrospective analysis of RT-PCR results obtained using the TaqPath kit from August 3rd to December 20th. We selected only positive samples with a Ct value < 25 for the N target, the most sensitive target of the test. By doing so, we found that 59/9,266 (0.6%) of positive tests had no amplification of the S gene. No significant increase of the SGTF was noticed over time; the proportion ranging from 0% (week # 32, 33, 34, 42, 48-51) to 2.91% (week # 35; Figure 2). Among the 59 samples with SGTF, 36 were available for WGS. These 36 samples were collected from August 5th to November 11th (18/36 were collected after October 9th). A total of 11 samples that presented an amplification of the S target were also sequenced. The sequencing results were fully concordant with the RT-PCR profiles (100% of the samples with SGTF had the S deletion ΔH69/ΔV70, while 100% of the S-gene positive samples did not contain ΔH69/ΔV70). For samples with SGTF, other S mutations were detected and are summarized in Table 1. The most frequent S mutations co-occurring with ΔH69/ΔV70 deletion were S477N and D614G that were found in 21/36 samples (58.3%). The co-occurrence of N439K and D614G mutations was found in 10/36 samples; the first sample containing this combination of mutation was collected on August 5th. Of note the complete combination of S mutations detected in cluster-5 variant was not found. View this table: [Table 1.](http://medrxiv.org/content/early/2021/01/11/2020.11.10.20228528/T1) Table 1. Spike mutations co-occurring with ΔH69/ΔV70 deletion in 36 samples with S negative profiles (negative for S target and positive for N & ORF1ab targets) obtained with the RT-PCR TaqPath kit. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/11/2020.11.10.20228528/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/01/11/2020.11.10.20228528/F2) Figure 2. Prevalence of the S negative profile (negative for S target and positive for N & ORF1ab targets) with TaqPath COVID-19 RT-PCR kit from August 3rd (week 32) to December 20th (week 51). The 2-step strategy presented herein and based on a screening with TaqPath Kit followed by WGS for samples with SGTF has been implemented in France since December 20th. On December 21th, the virology laboratory of university hospital of Tours reported a SGTF on a nasopharyngeal sample from a patient with a recent travel history from England (London). The sample was addressed to NRC for WGS and the detection of VOC 202012/01 (lineage B.1.1.7) was confirmed on December 25th that corresponded to the first detection of this variant in France (GISAID accession number EPI_ISL_735391). ## Discussion and conclusion According to CoV-GLUE resource [8] (last update from GISAID: December 14th), the S deletion 21765-21770 has been identified in 4,632 sequences worldwide (>99% in Europe) Interestingly, only 16 sequences containing this deletion were sampled between March 15th and July 23th corresponding to the first wave of COVID-19 pandemic in Europe. Herein, using data obtained with TaqPath RT-PCR kit, we found an overall prevalence of 0.6%, suggesting a limited circulation of variants presenting the ΔH69/ΔV70 deletion during the second wave of the pandemic in Lyon, France. It should be underlined that N439K, Y453F, or N501Y RBD mutations that can co-occurred with ΔH69/ΔV70 deletion might be associated with an increased affinity to ACE2 or reduced sensitivity to SARS-CoV-2 antibodies [3, 5, 10–12]. It has been hypothesized that the ΔH69/ΔV70 deletion might compensate some RBD mutations and might be involved in the transmissibility of variant containing these mutations [4, 6]. In addition, it has been recently shown that the combined ΔH69/V70 and D796H mutant was less sensitive to neutralizing antibodies [13]. As the N-terminal domain may interact with lung receptors [14] and might be a target of neutralizing antibodies [15, 16], further studies are needed to understand the consequences of ΔH69/ΔV70 deletion on SARS-CoV-2 transmissibility and host-immune response. Importantly, the TaqPath kit used for this study did not lead to a false negative conclusion as the two other targets remain positive. The data presented herein emphasize that the TaqPath RT-PCR assay is a useful and cost-effective tool enabling a rapid, large-scale screening of SARS-CoV-2 variants with ΔH69/ΔV70. Samples with SGTF should be further addressed to national referral laboratories for SARS-CoV-2 WGS. This 2-step strategy can contribute to the early detection of SARS-CoV-2 VOC 202012/01 which has been found to be more transmissible than non-VOC lineage [17]. This strategy is currently being reinforced in France as national diagnostic platforms have mainly implemented the TaqPath RT-PCR kit. ## Data Availability SARS-CoV-2 genomes sequenced in this study were deposited in the GISAID database (EPI\_ISL\_582110, EPI\_ISL\_582111, EPI\_ISL\_582112, EPI\_ISL\_582113, EPI\_ISL\_582114, EPI\_ISL\_582115, EPI\_ISL\_582116, EPI\_ISL\_582117, EPI\_ISL\_582118, EPI\_ISL\_582119, EPI\_ISL\_582120, EPI\_ISL\_582508, EPI\_ISL\_623098, EPI\_ISL\_623099, EPI\_ISL\_623100, EPI\_ISL_623101, EPI_ISL_623102) ## Data availability SARS-CoV-2 whole genomes sequenced in this study were deposited in the GISAID database (EPI\_ISL\_582110, EPI\_ISL\_582111, EPI\_ISL\_582112, EPI\_ISL\_582113, EPI\_ISL\_582114, EPI\_ISL\_582115, EPI\_ISL\_582116, EPI\_ISL\_582117, EPI\_ISL\_582118, EPI\_ISL\_582119, EPI\_ISL\_582120, EPI\_ISL\_582508, EPI\_ISL\_623098, EPI\_ISL\_623099, EPI\_ISL\_623100, EPI\_ISL\_623101, EPI\_ISL_623102, EPI_ISL_735391) ## Ethics statement Samples used in this study were collected as part of an approved ongoing surveillance conducted by the national reference centre for respiratory viruses in Lyon, France (WHO reference laboratory providing confirmatory testing for COVID-19). The investigations were carried out in accordance with the General Data Protection Regulation (Regulation (EU) 2016/679 and Directive 95/46/EC) and the French data protection law (Law 78–17 on 06/01/1978 and Décret 2019–536 on 29/05/2019). Samples were collected for regular clinical management during hospital stay, with no additional samples for the purpose of this study. Patients were informed of the research and their non-objection approval was confirmed. This study was presented by the ethics committee of the Hospices Civils de Lyon (HCL), Lyon, France and registered on the HCL database of RIPHN studies (AGORA N°41). * Received November 10, 2020. * Revision received January 11, 2021. * Accepted January 11, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.ECDC to assess risk associated with spread of SARS-CoV-2 in mink farms. European Centre for Disease Prevention and Control. 2020. [https://www.ecdc.europa.eu/en/news-events/ecdc-assess-risk-associated-spread-sars-cov-2-mink-farms](https://www.ecdc.europa.eu/en/news-events/ecdc-assess-risk-associated-spread-sars-cov-2-mink-farms). Accessed 8 Nov 2020. 2. 2.WHO | SARS-CoV-2 mink-associated variant strain – Denmark. WHO. [http://www.who.int/csr/don/06-november-2020-mink-associated-sars-cov2-denmark/en/](http://www.who.int/csr/don/06-november-2020-mink-associated-sars-cov2-denmark/en/). Accessed 25 Nov 2020. 3. 3. Lassaunière Ria, Fonager Jannik, Rasmussen Morten, Frische Anders, Polacek Strandh Charlotta, Bruun Rasmussen Thomas, Bøtner Anette, Fomsgaard Anders. Working paper on SARS-CoV-2 spike mutations arising in Danish mink, their spread to humans and neutralization data.[https://files.ssi.dk/Mink-cluster-5-short-report\_AFO2](https://files.ssi.dk/Mink-cluster-5-short-report_AFO2). 4. 4.Kemp SA, Datir RP, Collier DA, Ferreira I, Carabelli A, Harvey W, et al. Recurrent emergence and transmission of a SARS-CoV-2 Spike deletion ΔH69/ΔV70. bioRxiv. 2020;:2020.12.14.422555. 5. 5.Thomson EC, Rosen LE, Shepherd JG, Spreafico R, Filipe A da S, Wojcechowskyj JA, et al. The circulating SARS-CoV-2 spike variant N439K maintains fitness while evading antibody-mediated immunity. bioRxiv. 2020;:2020.11.04.355842. 6. 6. Andrew Rambaut, Nick Loman, Oliver Pybus, Wendy Barclay, Jeff Barrett, Alesandro Carabelli, Tom Connor, Tom Peacock, David L Robertson, Erik Volz, on behalf of COVID-19 Genomics Consortium UK (CoG-UK). Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations. Virological. [https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563/1](https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563/1). Accessed 21 Dec 2020. 7. 7.Bal A, Destras G, Gaymard A, Bouscambert-Duchamp M, Valette M, Escuret V, et al. Molecular characterization of SARS-CoV-2 in the first COVID-19 cluster in France reveals an amino acid deletion in nsp2 (Asp268del). Clin Microbiol Infect. 2020;26:960–2. 8. 8.Singer J, Gifford R, Cotten M, Robertson D. CoV-GLUE: A Web Application for Tracking SARS-CoV-2 Genomic Variation. 2020. doi:10.20944/preprints202006.0225.v1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20944/preprints202006.0225.v1&link_type=DOI) 9. 9.Charre C, Ginevra C, Sabatier M, Regue H, Destras G, Brun S, et al. Evaluation of NGS-based approaches for SARS-CoV-2 whole genome characterisation. Virus Evol. 2020;6:veaa075. 10. 10.Weisblum Y, Schmidt F, Zhang F, DaSilva J, Poston D, Lorenzi JC, et al. Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants. Elife. 2020;9. 11. 11.Baum A, Fulton BO, Wloga E, Copin R, Pascal KE, Russo V, et al. Antibody cocktail to SARS-CoV-2 spike protein prevents rapid mutational escape seen with individual antibodies. Science. 2020;369:1014–8. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwNi8xMDE0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMTEvMjAyMC4xMS4xMC4yMDIyODUyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 12. 12.Starr TN, Greaney AJ, Hilton SK, Ellis D, Crawford KHD, Dingens AS, et al. Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding. Cell. 2020;182:1295-1310.e20. 13. 13.Kemp SA, Collier DA, Datir R, Gayed S, Jahun A, Hosmillo M, et al. Neutralising antibodies drive Spike mediated SARS-CoV-2 evasion. medRxiv. 2020;:2020.12.05.20241927. 14. 14.Soh WT, Liu Y, Nakayama EE, Ono C, Torii S, Nakagami H, et al. The N-terminal domain of spike glycoprotein mediates SARS-CoV-2 infection by associating with L-SIGN and DC-SIGN. bioRxiv. 2020;:2020.11.05.369264. 15. 15.Chi X, Yan R, Zhang J, Zhang G, Zhang Y, Hao M, et al. A neutralizing human antibody binds to the N-terminal domain of the Spike protein of SARS-CoV-2. Science. 2020;369:650– 5. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwNC82NTAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8xMS8yMDIwLjExLjEwLjIwMjI4NTI4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 16. 16.McCarthy KR, Rennick LJ, Nambulli S, Robinson-McCarthy LR, Bain WG, Haidar G, et al. Natural deletions in the SARS-CoV-2 spike glycoprotein drive antibody escape. bioRxiv. 2020;:2020.11.19.389916. 17. 17. Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data. Virological. 2020. [https://virological.org/t/transmission-of-sars-cov-2-lineage-b-1-1-7-in-england-insights-from-linking-epidemiological-and-genetic-data/576](https://virological.org/t/transmission-of-sars-cov-2-lineage-b-1-1-7-in-england-insights-from-linking-epidemiological-and-genetic-data/576). Accessed 5 Jan 2021.