key: cord-0976519-aeqhdu9f authors: Bal, Antonin; Destras, Gregory; Gaymard, Alexandre; Stefic, Karl; Marlet, Julien; Eymieux, Sébastien; Regue, Hadrien; Semanas, Quentin; d’Aubarede, Constance; Billaud, Geneviève; Laurent, Frédéric; Gonzalez, Claudia; Mekki, Yahia; Valette, Martine; Bouscambert, Maude; Gaudy-Graffin, Catherine; Lina, Bruno; Morfin, Florence; Josset, Laurence; Casalegno, Jean-Sébastien; Frobert, Emilie; Escuret, Vanessa; Icard, Vinca; Jeannoel, Marion; Milon, Marie-Paule; Ramière, Christophe; Scholtès, Caroline; Tardy, Jean-Claude; Trabaud, Mary-Anne; Schuffenecker, Isabelle title: 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 date: 2021-01-21 journal: Euro Surveill DOI: 10.2807/1560-7917.es.2021.26.3.2100008 sha: b0eb926c3d36dfc305f80e588535adfddfe3304d doc_id: 976519 cord_uid: aeqhdu9f We report the strategy leading to the first detection of variant of concern 202012/01 (VOC) in France (21 December 2020). First, the spike (S) deletion H69–V70 (ΔH69/ΔV70), identified in certain SARS-CoV-2 variants including VOC, is screened for. This deletion is associated with a S-gene target failure (SGTF) in the three-target RT-PCR assay (TaqPath kit). Subsequently, SGTF samples are whole genome sequenced. This approach revealed mutations co-occurring with ΔH69/ΔV70 including S:N501Y in the VOC. We report the strategy leading to the first detection of variant of concern 202012/01 (VOC) in France (21 December 2020). First, the spike (S) deletion H69-V70 (ΔH69/ΔV70), identified in certain SARS-CoV-2 variants including VOC, is screened for. This deletion is associated with a S-gene target failure (SGTF) in the threetarget RT-PCR assay (TaqPath kit). Subsequently, SGTF samples are whole genome sequenced. This approach revealed mutations co-occurring with ΔH69/ΔV70 including S:N501Y in the VOC. Since September 2020 a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) deletion H69-V70 (ΔH69/ΔV70) in the spike (S) protein 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 neutralising antibodies of sera from recovered coronavirus disease patients [1] [2] [3] . The ΔH69/ΔV70 has also co-occurred with either one of two other noteworthy RBD mutations [4] : N439K that is currently spreading in Europe and might also be related to reduced susceptibility to SARS-CoV-2 antibodies [5] and N501Y that was identified for example in 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 two-step strategy that enabled the timely detection of VOC 202012/01 in France, as well as other variants carrying ΔH69/ΔV70. As part of routine SARS-CoV-2 genomic surveillance performed at the national reference centre (NRC) for respiratory viruses (Lyon, France) [7] , a 6-nt deletion (21765-21770) within the S gene was identified in two nasopharyngeal samples collected in Lyon, France on 1 and 7 September 2020, respectively. The SARS-CoV-2 infection diagnosis had been performed on 2 and 8 September 2020, respectively with the Applied Biosystems TaqPath RT-PCR COVID-19 kit (Thermo Fisher Scientific, Waltham, United States (US)) that includes the open reading frame (ORF) 1ab, S, and nucleocapsid (N) gene targets. For these two samples, a S-gene target failure (SGTF) was observed while ORF1ab and N targets was correctly amplified with quantification cycle (Cq) values < 25 ( Figure 1A ). The two samples were subjected to whole genome sequencing (WGS) and the sequences were deposited in the GISAID Initiative's EpiCOV database on 15 October 2020 (GISAID accession numbers: EPI_ ISL_582112, EPI_ISL_582120). The mean coverage was 6903× and 6898×, 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 two amino acids (ΔH69/ΔV70) in the N-terminal domain of the S1 subunit of the S protein ( Figure 1B ). The WGS method used was the amplicon-based ARTIC v3 protocol (https://artic.network/ncov-2019) combined with Nextera DNA Flex library and sequencing on NextSeq 500 platform (Illumina, San Diego, US). To confirm 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. We then performed a retrospective analysis of the TaqPath kit results obtained at a community-based testing platform hosted by the virology laboratory of Lyon university hospital, from 3 August to 20 December 2020. We selected only the samples with a Cq value < 25 for the N target, the most sensitive target of the test. Of note, a Cq value of 25 with the TaqPath kit corresponds to a cycle threshold (Ct) value of 30 with the real time RT-PCR IP4 (Institut Pasteur assay, Paris, France) [10] . By doing so, we found that 59 of 9,266 (0.6%) 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 numbers 32, 33, 34, 42, 48-51) to 2.84% (week number 37; Figure 2 ). Among the 59 samples with SGTF, 36 were available for WGS. These 36 samples were collected from 5 August to 11 November (18/36 were collected after 9 October). A total of 11 samples that presented an amplification of the S target were also sequenced. Week 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). According to CoV-GLUE resource [8] (update from GISAID: 14 December 2020), 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 15 March and 23 July 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 ΔH69/ΔV70 during the second wave of the pandemic in Lyon, France. Of note, at the time of the study no travel restriction measures were implemented in France. The two-step strategy presented herein enabled to identify SARS-CoV-2 variants carrying ΔH69/ ΔV70 and RBD mutations. More importantly, this strategy allowed the first detection of the VOC 202012/01 in France. It should be underlined that N439K, Y453F or N501Y RBD mutations that can co-occur with ΔH69/ΔV70 might be associated with an increased affinity to angiotensinconverting enzyme 2 (ACE2) or reduced sensitivity to SARS-CoV-2 antibodies [3, 5, [11] [12] [13] . The N501Y mutation of the VOC 202012/01, in particular, might be also responsible of the higher transmissibility reported for this lineage [14] . It has been hypothesised that ΔH69/ΔV70 might compensate some RBD mutations and might be involved in the transmissibility of the variants containing these mutations [4, 6] . In addition, it has been recently shown that the combined ΔH69/V70 and D796H mutant was less sensitive to neutralising antibodies [15] . As the N-terminal domain of the S1 subunit of the S protein may interact with lung receptors [16] and might be a target of neutralising antibodies [17, 18] , further studies are needed to understand the consequences of ΔH69/ ΔV70 on SARS-CoV-2 transmissibility and host-immune response. The present study has several limitations. First we selected only samples with a Cq value <25 for the N gene that corresponds to the limit of sensitivity of the WGS method used. The detection of VOC 202012/01 could therefore be underestimated, but not the prevalence. In addition, this strategy cannot identify other variants of concern that do not carry ΔH69/ΔV70 as the variant 501Y.V2 first detected in South Africa [19] . However, it allows rapid retrospective and prospective evaluation of the prevalence of VOC 202012/01, including in countries with limited sequencing capacity. Importantly, the TaqPath kit did not lead to false negative conclusions regarding the SARS-CoV-2 diagnosis as the two other targets remained positive. The data presented herein emphasise 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. In the meantime, a similar approach has been adopted in the United Kingdom where the frequency of SGTF is used as a proxy for estimating and monitoring the spread of the VOC 202012/01 [14, 20] . Samples with SGTF should be further addressed to national referral laboratories for SARS-CoV-2 WGS in order to confirm the presence of VOC 202012/01. This two-step strategy can contribute to the timely detection and isolation of VOC 202012/01 cases and has been reinforced in France as national diagnostic platforms have mainly implemented the TaqPath RT-PCR kit. ECDC to assess risk associated with spread of SARS-CoV-2 in mink farms. Stockholm: ECDC; 2020 ecdc-assess-risk-associated-spread-sars-cov-2-mink-farms 2. World Health Organization (WHO) Working paper on SARS-CoV-2 spike mutations arising in Danish mink, their spread to humans and neutralization data Recurrent emergence and transmission of a SARS-CoV-2 Spike deletion ΔH69/ΔV70 The circulating SARS-CoV-2 spike variant N439K maintains fitness while evading antibodymediated immunity Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations. Virological.org Molecular characterization of SARS-CoV-2 in the first COVID-19 cluster in France reveals an amino acid deletion in nsp2 (Asp268del) CoV-GLUE: A Web Application for Tracking SARS-CoV-2 Genomic Variation Evaluation of NGS-based approaches for SARS-CoV-2 whole genome characterisation Molecular assays to diagnose COVID-19: Summary table of available protocols. Geneva: WHO Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants. eLife Antibody cocktail to SARS-CoV-2 spike protein prevents rapid mutational escape seen with individual antibodies Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data Neutralising antibodies drive Spike mediated SARS-CoV-2 adaptation The N-terminal domain of spike glycoprotein mediates SARS-CoV-2 infection by associating with L-SIGN and DC-SIGN A neutralizing human antibody binds to the N-terminal domain of the Spike protein of SARS-CoV-2 Recurrent deletions in the SARS-CoV-2 spike glycoprotein drive antibody escape Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa Investigation of novel SARS-COV-2 variant. Variant of Concern 202012/01. London: PHE We would like to thank all microbiologist, clinicians, laboratory technicians who contributed to this investigation. We thank the authors, the originating and submitting laboratories for their sequence and metadata shared through GISAID on which this research is based. We gratefully acknowledge all the members and authors of CoV-GLUE, Nextstrain.org, and Virological.org for sharing their analysis in real-time. We thank Philip Robinson (DRCI, Hospices Civils de Lyon) for his help in manuscript preparation. This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY 4.0) Licence. You may share and adapt the material, but must give appropriate credit to the source, provide a link to the licence and indicate if changes were made.Any supplementary material referenced in the article can be found in the online version.