2 didier.raoult@gmail.com 26
5 variant involved in almost all SARS-CoV-2 infections at that time (Table 1). Respiratory 78 samples collected from seven other SARS-CoV-2-positive patients living in the same 79 geographical area exhibited the same combination of mutations screened by qPCR. They were 80 two adults and five children (<15 years of age) (Table 1). The respiratory samples from these 81 eight patients were sent to university hospital institute Méditerranée Infection for SARS-CoV-82 2 genome sequencing as recommended by French public health authorities. A rapid NGS 83 procedure was launched overnight. It allowed obtaining SARS-CoV-2 genotype identification 84 in
6 genomes were classified into Nextclade and Pangolin lineages using web applications 103 (https://clades.nextstrain.org/;https://cov-lineages.org/pangolin.html) [10,11,13]. They were 104 deposited in the GISAID sequence database (https://www.gisaid.org/) [14] (Table 1). 105 Phylogenies were reconstructed with the nextstrain/ncov tool 106 (https://github.com/nextstrain/ncov) then visualized with Auspice 107 (https://docs.nextstrain.org/projects/auspice/en/stable/). 108 The analysis of viral genomes revealed the presence of 46 nucleotide substitutions and 109 37 deletions, resulting in 30 amino acid substitutions and 12 deletions (Figure 1a; 110 Supplementary Tables S1 and S2). Fourteen amino acid substitutions and 9 amino acid 111 deletions are located in the spike protein. Substitutions N501Y and E484K are combined as in 112 the Beta, Gamma, Theta and Omicron variants [5,15]. Substitution F490S is present as in the 113 Lambda variant, and substitution P681H is present as in the Lambda and Omicron variants. In 114 other structural proteins than the spike, amino acid changes include two substitutions in the 115 nucleocapsid protein and one in the membrane protein. In non-structural proteins, amino acid 116 changes include one substitution in proteins Nsp2, Nsp3, Nsp4, Nsp6, Nsp12 (RNA-117 dependent RNA polymerase), and Nsp13 (helicase); two substitutions in Nsp14 (
7 https://genome.ucsc.edu/cgi-bin/hgPhyloPlace) option, which showed the phylogenetic 128 placement of the genomes we obtained as an outgroup of the B.1.640 lineage and their 129 clustering with a genome obtained late October in France (Ile-de-France) 130 (EPI_ISL_5926666). The B.1.640 lineage corresponds to a variant first identified in France in 131 April 2021, in Indonesia in August 2021, and in Republic of the Congo (Brazzaville) in 132 September 2021, and it was involved in a cluster of cases in Brittany, France around mid-133 October 2021 [17]. As of 09/12/2021, 157 genomes were available from the GISAID database 134 including 92 from France and 36 from the Republic of the Congo. The sets of spike mutations 135 of the B.1.640 lineage and of genomes obtained here are similar, with 11 common nucleotide 136 substitutions and 1 common deletion of 9 codons (Supplementary Figure S1, Tables S1-2). 137 However, spike genes of both lineages differ by 7 mutations. In addition, 25 nucleotide 138 substitutions and 33 nucleotide deletions located elsewhere in the genomes differ between the 139 two genotypes. The pattern of mutations of present genomes hence indicates a new
8 stay in Cameroon, we sought for this variant in GISAID among genomes from this country 153 but as of 09/12/2021 none of the 556 available genomes belong to the B.1.640.1 or B.1.640.2 154 lineages. 155 We analyzed a complete structure of the spike protein of the IHU variant generated by 156 incorporating its specific mutational profile to the original 20B SARS-CoV-2 (Wuhan-Hu-1 157 isolate with D614G substitution) [18] and fixing all gaps in the pdb file by incorporating the 158 missing amino acids with the Robetta protein structure prediction tool 159 [https://robetta.bakerlab.org/], followed by energy minimization with the Polak-Ribière 160 algorithm as previously reported (Figure 1c) [19]. In the N-terminal domain (NTD), the 134-161 145 amino acid deletion is predicted to significantly affect the neutralizing epitope. Other 162 changes involve amino acids at positions 96 and 190: in Wuhan-Hu-1 isolate, E96 and R190 163 induce a turn in NTD secondary structure through electrostatic interactions between each 164 other. This interaction is conserved between substituted amino acids 96Q and 190S, which 165 suggests the co-evolution of these changes. In the receptor binding domain (RBD), aside the 166 well-known substitutions N501Y and E484K, several changes were predicted to significantly 167 affect the neutralizing epitopes. Particularly, P681H is located in the cleavage site of S1-S2 168 subunits of the spike and is observed in other variants including the recently emerging 169 Omicron [15]. Besides, D1139H substitution implies an amino acid involved in the fusion 170 between the virus and the infected cell.Also, D614G is combined with T859N in the IHU 171 variant. Interestingly, in the Wuhan-Hu-1 isolate, amino acids D614 and T859 from two 172 subunits of the trimeric spike are face to face and lock the trimer in a closed conformation. 173 Substitution D614G allows unlocking the trimer conformation, but this is predicted to be still 174 easier in case of additional presence of substitution T859N. 175 Respiratory samples collected until end of November 2021 from four other SARS-176 CoV-2 positive patients living in the same city or borough than the index case could be 177
9 identified as containing the IHU variant by NGS within 24 hours after their reception (Table 178 1). All 12 IHU variant-positive samples showed the same combination of spike mutations as 179 screened by real-time qPCR techniques: negativity for 452R and 484Q; positivity for 484K, 180 501Y [20], and 681H [3]. We also used the TaqPath COVID-19 kit (Thermo Fisher Scientific, 181 Waltham, USA) that provided positive signals for all three genes targeted (ORF1, S, and N). 182 Thus, the IHU variant can be distinguished by screening with qPCR assays from the Delta 183 (L452R-positive) and Omicron (L452R-negative and negative for S gene detection by the 184 TaqPath COVID-19 assay) variants that currently co-circulate in our geographical area. 185 Finally, scanning electron microscopy using a SUV 5000 microscope (Hitachi High-186 Technologies Corporation, Tokyo, Japan) [21] allowed a quick visualization of the virus from 187 a respiratory sample (Figure 1d). 188 189 Overall, these observations show once again the unpredictability of the emergence of 190 new SARS-CoV-2 variants and their introduction from abroad, and they exemplify the 191 difficulty to control such introduction and subsequent spread. They also warrant the 192 implementation of genomic surveillance of SARS-CoV-2 that we started from the very 193 beginning of the pandemic in our geographical area as soon as we diagnosed the first SARS-194 CoV-2 infection [21] and that we expanded during summer 2020 [2,3]. This surveillance has 195 been implemented at the country scale in 2021 through the French Emergen consortium 196 (https://www.santepubliquefrance.fr/dossiers/coronavirus-covid-19/consortium-emergen). It is 197 too early to speculate on virological, epidemiological or clinical features of this IHU variant 198 based on these 12 cases. For this purpose, respiratory samples from infected patients were 199 inoculated on Vero E6 cells as previously described [22] to be able assessing the sensibility to 200 neutralization by anti-spike antibodies elicited by vaccine immunization, or by prior infection 201 [23]. 202
13 resistant epitopes targeted by SARS-CoV-2 antibodies: A global consortium study. 265 Science 374: 472-478. 266 5. Harvey WT, Carabelli AM, Jackson B, Gupta RK, Thomson EC, Harrison EM, Ludden 267 C, Reeve R, Rambaut A (2021) COVID-19 Genomics UK (COG-UK) Consortium, 268 Peacock SJ, Robertson DL. SARS-CoV-2 variants, spike mutations and immune escape. 269 Nat Rev Microbiol 19: 409-424. 270 6. Tao K, Tzou PL, Nouhin J, Gupta RK, de Oliveira T, Kosakovsky Pond SL, Fera D, 271 Shafer RW (2021) The biological and clinical significance of emerging SARS-CoV-2 272 variants. Nat Rev Genet 22: 757-773. 273 7. Wilder-Smith A (2021) What is the vaccine effect on reducing transmission in the 274 context of the SARS-CoV-2 delta variant? Lancet Infect Dis S1473-3099(21)00690-3. 275 doi: 10.1016/S1473-3099(21)00690-3. Epub ahead of print. PMID: 34756187; PMCID: 276 PMC8554481. 277 8. Li H (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 278 34: 3094-3100. 279 9. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, 280 Durbin R, 1000 Genome Project Data Processing Subgroup (2009). The Sequence 281 Alignment/Map format and SAMtools. Bioinformatics 25: 2078-2079. 282 10. Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C, Sagulenko P, 283 Bedford T, Neher RA (2018). Nextstrain: real-time tracking of pathogen evolution. 284 Bioinformatics 34: 4121-4123. 285 11. Aksamentov I, Roemer C, Hodcroft EB, Neher RA (2021). Nextclade: clade 286 assignment, mutation calling and quality control for viral genomes. Zenodo 287 https://doi.org/10.5281/zenodo.5607694. 288
14 12. Garrison E, Marth G (2012). Haplotype-based variant detection from short-read 289 sequencing. arXiv.org. https://arxiv.org/abs/1207.3907 (accessed 10 December 2021). 290 13. Rambaut A, Holmes EC, O’Toole Ã, Hill V, McCrone JT, Ruis C, du Plessis L, Pybus 291 OG (2020). A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist 292 genomic epidemiology. Nat. Microbiol 5: 1403-1407. 293 14. Alm E, Broberg EK, Connor T, Hodcroft EB, Komissarov AB, Maurer-Stroh S, 294 Melidou A, Neher RA, O’Toole A, Pereyaslov D, WHO European Region sequencing 295 laboratories and GISAID EpiCoV group; WHO European Region sequencing 296 laboratories and GISAID EpiCoV group (2020) Geographical and temporal distribution 297 of SARS-CoV-2 clades in the WHO European Region, January to June 2020. Euro. 298 Surveill 25: 2001410. 299 15. Karim SSA, Karim QA (2021) Omicron SARS-CoV-2 variant: a new chapter in the 300 COVID-19 pandemic. Lancet 398: 2126-2128. 301 16. Rambaut A, Loman N, Pybus O, Barclay W, Barrett J, Carabelli A, Connor T, Peacock 302 T, Robertson DL, Volz E, on behalf of COVID-19 Genomics Consortium UK (CoG-303 UK) (2020) Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage 304 in the UK defined by a novel set of spike mutations. Virological Pre-print. 305 https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-306 2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563 307 17. Le Page M (2021) New variant gains ground. New Sci 252: 8. 308 18. Benton DJ, Wrobel AG, Roustan C, Borg A, Xu P, Martin SR, Rosenthal PB, Skehel JJ, 309 Gamblin SJ (2021). The effect of the D614G substitution on the structure of the spike 310 glycoprotein of SARS-CoV-2. Proc. Natl. Acad. Sci. U S A. 118:e2022586118. 311
10 203 204
12
13 resistant epitopes targeted by SARS-CoV-2 antibodies: A global consortium study. 265 Science 374: 472-478. 266 5. Harvey WT, Carabelli AM, Jackson B, Gupta RK, Thomson EC, Harrison EM, Ludden 267 C, Reeve R, Rambaut A (2021) COVID-19 Genomics UK (COG-UK) Consortium, 268 Peacock SJ, Robertson DL. SARS-CoV-2 variants, spike mutations and immune escape. 269 Nat Rev Microbiol 19: 409-424. 270 6. Tao K, Tzou PL, Nouhin J, Gupta RK, de Oliveira T, Kosakovsky Pond SL, Fera D, 271 Shafer RW (2021) The biological and clinical significance of emerging SARS-CoV-2 272 variants. Nat Rev Genet 22: 757-773. 273 7. Wilder-Smith A (2021) What is the vaccine effect on reducing transmission in the 274 context of the SARS-CoV-2 delta variant? Lancet Infect Dis S1473-3099(21)00690-3. 275 doi: 10.1016/S1473-3099(21)00690-3. Epub ahead of print. PMID: 34756187; PMCID: 276 PMC8554481. 277 8. Li H (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 278 34: 3094-3100. 279 9. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, 280 Durbin R, 1000 Genome Project Data Processing Subgroup (2009). The Sequence 281 Alignment/Map format and SAMtools. Bioinformatics 25: 2078-2079. 282 10. Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C, Sagulenko P, 283 Bedford T, Neher RA (2018). Nextstrain: real-time tracking of pathogen evolution. 284 Bioinformatics 34: 4121-4123. 285 11. Aksamentov I, Roemer C, Hodcroft EB, Neher RA (2021). Nextclade: clade 286 assignment, mutation calling and quality control for viral genomes. Zenodo 287 https://doi.org/10.5281/zenodo.5607694. 288
Πρόσφατα Σχόλια