S:D614G and S:H655Y are gateway mutations that act epistatically to promote SARS-CoV-2 variant fitness

SARS-CoV-2 variants bearing complex combinations of mutations that confer increased transmissibility, COVID-19 severity, and immune escape, were first detected after S:D614G had gone to fixation, and likely originated during persistent infection of immunocompromised hosts. To test the hypothesis that S:D614G facilitated emergence of such variants, S:D614G was reverted to the ancestral sequence in the context of sequential Spike sequences from an immunocompromised individual, and within each of the major SARS-CoV-2 variants of concern. In all cases, infectivity of the S:D614G revertants was severely compromised. The infectivity of atypical SARS-CoV-2 lineages that propagated in the absence of S:D614G was found to be dependent upon either S:Q613H or S:H655Y. Notably, Gamma and Omicron variants possess both S:D614G and S:H655Y, each of which contributed to infectivity of these variants. Among sarbecoviruses, S:Q613H, S:D614G, and S:H655Y are only detected in SARS-CoV-2, which is also distinguished by a polybasic S1/S2 cleavage site. Genetic and biochemical experiments here showed that S:Q613H, S:D614G, and S:H655Y each stabilize Spike on virions, and that they are dispensable in the absence of S1/S2 cleavage, consistent with selection of these mutations by the S1/S2 cleavage site. CryoEM revealed that either S:D614G or S:H655Y shift the Spike receptor binding domain (RBD) towards the open conformation required for ACE2-binding and therefore on pathway for infection. Consistent with this, an smFRET reporter for RBD conformation showed that both S:D614G and S:H655Y spontaneously adopt the conformation that ACE2 induces in the parental Spike. Data from these orthogonal experiments demonstrate that S:D614G and S:H655Y are convergent adaptations to the polybasic S1/S2 cleavage site which stabilize S1 on the virion in the open RBD conformation and act epistatically to promote the fitness of variants bearing complex combinations of clinically significant mutations. Highlights S:D614G is ubiquitous among SARS-CoV-2 B-lineage Spikes and is required for infectivity of the main Variants of Concern In an example of convergent evolution, SARS-CoV-2 A lineage viruses maintained transmission chains in the absence of S:D614G, but were instead dependent upon S:Q613H or S:H655Y S:D614G and S:H655Y are both adaptations to the polybasic S1/S2 cleavage site Increased infectivity of S:D614G and S:H655Y is associated with a more open RBD conformation and increased steady-state levels of virion-associated S1

[1]  A. Finzi,et al.  Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission , 2022, Nature Microbiology.

[2]  C. Davis,et al.  Cross-neutralization and viral fitness of SARS-CoV-2 Omicron sublineages , 2022, bioRxiv.

[3]  T. de Oliveira,et al.  Persistent SARS-CoV-2 infection with accumulation of mutations in a patient with poorly controlled HIV infection , 2022, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[4]  Shuntai Zhou,et al.  Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection , 2023, Cell Reports Medicine.

[5]  Pardis C Sabeti,et al.  Analysis of 6.4 million SARS-CoV-2 genomes identifies mutations associated with fitness , 2022, Science.

[6]  N. Loman,et al.  The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK , 2022, bioRxiv.

[7]  G. Ye,et al.  Cryo-EM structure of a SARS-CoV-2 omicron spike protein ectodomain , 2022, Nature Communications.

[8]  L. Poon,et al.  Transmission of SARS-CoV-2 delta variant (AY.127) from pet hamsters to humans, leading to onward human-to-human transmission: a case study , 2022, The Lancet.

[9]  H. Mo,et al.  Remdesivir and GS-441524 Retain Antiviral Activity against Delta, Omicron, and Other Emergent SARS-CoV-2 Variants , 2022, bioRxiv.

[10]  Liyuan Liu,et al.  Cryo-EM structure of the SARS-CoV-2 Omicron spike , 2022, Cell Reports.

[11]  N. Krogan,et al.  Mutations in SARS-CoV-2 variants of concern link to increased spike cleavage and virus transmission , 2022, Cell Host & Microbe.

[12]  T. Ndung’u,et al.  SARS-CoV-2 prolonged infection during advanced HIV disease evolves extensive immune escape , 2022, Cell Host & Microbe.

[13]  A. Darby,et al.  Rapid selection of P323L in the SARS-CoV-2 polymerase (NSP12) in humans and non-human primate models and confers a large plaque phenotype , 2021, bioRxiv.

[14]  M. Nelson,et al.  SARS-CoV-2 infection in free-ranging white-tailed deer , 2021, Nature.

[15]  E. Ghedin,et al.  Yearlong COVID-19 Infection Reveals Within-Host Evolution of SARS-CoV-2 in a Patient With B-Cell Depletion , 2021, The Journal of infectious diseases.

[16]  M. Kraemer,et al.  Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa , 2021, Nature.

[17]  Graham W. Taylor,et al.  SARS-CoV-2 lineage dynamics in England from September to November 2021: high diversity of Delta sub-lineages and increased transmissibility of AY.4.2 , 2021, BMC Infectious Diseases.

[18]  G. Wallau,et al.  Identification of a novel SARS-CoV-2 P.1 sub-lineage in Brazil provides new insights about the mechanisms of emergence of variants of concern , 2021, Virus evolution.

[19]  Philip L. Tzou,et al.  Coronavirus Resistance Database (CoV-RDB): SARS-CoV-2 susceptibility to monoclonal antibodies, convalescent plasma, and plasma from vaccinated persons , 2021, bioRxiv.

[20]  Vineet D. Menachery,et al.  The N501Y spike substitution enhances SARS-CoV-2 infection and transmission , 2021, Nature.

[21]  N. French,et al.  SARS-CoV-2 exposure in Malawian blood donors: an analysis of seroprevalence and variant dynamics between January 2020 and July 2021 , 2021, BMC Medicine.

[22]  M. C. Muenker,et al.  De novo emergence of a remdesivir resistance mutation during treatment of persistent SARS-CoV-2 infection in an immunocompromised patient: a case report , 2021, Nature Communications.

[23]  James B. Munro,et al.  Conformational dynamics and allosteric modulation of the SARS-CoV-2 spike , 2021, bioRxiv.

[24]  G. Kochs,et al.  Antibody escape and global spread of SARS-CoV-2 lineage A.27 , 2021, Nature communications.

[25]  Graham W. Taylor,et al.  Exponential growth, high prevalence of SARS-CoV-2, and vaccine effectiveness associated with the Delta variant , 2021, Science.

[26]  E. Ghedin,et al.  Yearlong COVID-19 Infection Reveals Within-Host Evolution of SARS-CoV-2 in a Patient With B-Cell Depletion , 2021, medRxiv.

[27]  S. Scheres,et al.  New tools for automated cryo-EM single-particle analysis in RELION-4.0 , 2021, bioRxiv.

[28]  M. Nussenzweig,et al.  High genetic barrier to SARS-CoV-2 polyclonal neutralizing antibody escape , 2021, Nature.

[29]  Pardis C Sabeti,et al.  Analysis of 6.4 million SARS-CoV-2 genomes identifies mutations associated with fitness , 2021, medRxiv.

[30]  James B. Munro,et al.  Stability and conformation of the dimeric HIV-1 genomic RNA 5′UTR , 2021, Biophysical journal.

[31]  Nuno R. Faria,et al.  Global disparities in SARS-CoV-2 genomic surveillance , 2021, medRxiv.

[32]  O. Tsang,et al.  The impact of spike N501Y mutation on neutralizing activity and RBD binding of SARS-CoV-2 convalescent serum , 2021, EBioMedicine.

[33]  D. Wesemann,et al.  Membrane fusion and immune evasion by the spike protein of SARS-CoV-2 Delta variant , 2021, bioRxiv.

[34]  Vineet D. Menachery,et al.  Delta spike P681R mutation enhances SARS-CoV-2 fitness over Alpha variant , 2021, bioRxiv.

[35]  D. O’Connor,et al.  Acute SARS-CoV-2 infections harbor limited within-host diversity and transmit via tight transmission bottlenecks , 2021, PLoS pathogens.

[36]  N. Loman,et al.  Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic , 2021, Cell.

[37]  Nuno R. Faria,et al.  Neutralisation of SARS-CoV-2 lineage P.1 by antibodies elicited through natural SARS-CoV-2 infection or vaccination with an inactivated SARS-CoV-2 vaccine: an immunological study , 2021, The Lancet Microbe.

[38]  J. Dushoff,et al.  The origins and potential future of SARS-CoV-2 variants of concern in the evolving COVID-19 pandemic , 2021, Current Biology.

[39]  P. Kaleebu,et al.  Emergence and spread of a SARS-CoV-2 lineage A variant (A.23.1) with altered spike protein in Uganda , 2021, Nature Microbiology.

[40]  Manoj Kumar Singh,et al.  Genomic characterization and epidemiology of an emerging SARS-CoV-2 variant in Delhi, India , 2021, Science.

[41]  William T. Harvey,et al.  SARS-CoV-2 variants, spike mutations and immune escape , 2021, Nature Reviews Microbiology.

[42]  S. Whelan,et al.  SARS-CoV-2 Spreads through Cell-to-Cell Transmission , 2021, bioRxiv.

[43]  G. Barton,et al.  Effects of common mutations in the SARS-CoV-2 Spike RBD domain and its ligand the human ACE2 receptor on binding affinity and kinetics , 2021, bioRxiv.

[44]  G. Kochs,et al.  Within-host evolution of SARS-CoV-2 in an immunosuppressed COVID-19 patient as a source of immune escape variants , 2021, Nature Communications.

[45]  M. Giacca,et al.  The furin cleavage site in the SARS-CoV-2 spike protein is required for transmission in ferrets , 2021, Nature Microbiology.

[46]  M. Suchard,et al.  Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil , 2021, Science.

[47]  Rommie E. Amaro,et al.  Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN , 2021, bioRxiv.

[48]  T. Ndung’u,et al.  Escape of SARS-CoV-2 501Y.V2 from neutralization by convalescent plasma , 2021, Nature.

[49]  Graham W. Taylor,et al.  Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England , 2021, Nature.

[50]  J. Todd,et al.  SARS-CoV-2 within-host diversity and transmission , 2021, Science.

[51]  A. Godzik,et al.  Detection of a SARS-CoV-2 variant of concern in South Africa , 2021, Nature.

[52]  D. Fremont,et al.  Resistance of SARS-CoV-2 variants to neutralization by monoclonal and serum-derived polyclonal antibodies , 2021, Nature Medicine.

[53]  W. Hanage,et al.  SARS-CoV-2 Sequence Characteristics of COVID-19 Persistence and Reinfection , 2021, medRxiv.

[54]  M. Beer,et al.  SARS-CoV-2 spike D614G change enhances replication and transmission , 2021, Nature.

[55]  Nikolaos S. Alachiotis,et al.  Population Genomics Insights into the First Wave of COVID-19 , 2021, Life.

[56]  Graham W. Taylor,et al.  SARS-CoV-2 evolution during treatment of chronic infection , 2021, Nature.

[57]  D. Ho,et al.  Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7 , 2021, bioRxiv.

[58]  Vineet D. Menachery,et al.  Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis , 2021, Nature.

[59]  M. Nussenzweig,et al.  mRNA vaccine-elicited antibodies to SARS-CoV-2 and circulating variants , 2021, bioRxiv.

[60]  L. Morris,et al.  SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma , 2021, bioRxiv.

[61]  Haiyong Peng,et al.  SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity , 2020, Nature Communications.

[62]  T. Halasa,et al.  SARS-CoV-2 Transmission between Mink (Neovison vison) and Humans, Denmark , 2020, Emerging infectious diseases.

[63]  Lisa E. Gralinski,et al.  SARS-CoV-2 D614G variant exhibits efficient replication ex vivo and transmission in vivo , 2020, Science.

[64]  Gaurav D. Gaiha,et al.  Persistence and Evolution of SARS-CoV-2 in an Immunocompromised Host , 2020, The New England journal of medicine.

[65]  M. Koopmans,et al.  Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans , 2020, Science.

[66]  E. Fischer,et al.  Case Study: Prolonged Infectious SARS-CoV-2 Shedding from an Asymptomatic Immunocompromised Individual with Cancer , 2020, Cell.

[67]  Vineet D. Menachery,et al.  Spike mutation D614G alters SARS-CoV-2 fitness , 2020, Nature.

[68]  E. Koonin,et al.  Ongoing Global and Regional Adaptive Evolution of SARS-CoV-2 , 2020, bioRxiv.

[69]  Yongfei Cai,et al.  Structural impact on SARS-CoV-2 spike protein by D614G substitution , 2020, bioRxiv.

[70]  D. A. Jackson,et al.  Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity , 2020, Cell.

[71]  Edward C. Holmes,et al.  A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology , 2020, Nature Microbiology.

[72]  Pardis C Sabeti,et al.  Structural and Functional Analysis of the D614G SARS-CoV-2 Spike Protein Variant , 2020, bioRxiv.

[73]  S. Rowland-Jones,et al.  Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus , 2020, Cell.

[74]  F. Balloux,et al.  No evidence for increased transmissibility from recurrent mutations in SARS-CoV-2 , 2020, Nature Communications.

[75]  M. Hoffmann,et al.  A Multibasic Cleavage Site in the Spike Protein of SARS-CoV-2 Is Essential for Infection of Human Lung Cells , 2020, Molecular Cell.

[76]  Shinji Makino,et al.  An Infectious cDNA Clone of SARS-CoV-2 , 2020, Cell Host & Microbe.

[77]  K. Shi,et al.  Structural basis of receptor recognition by SARS-CoV-2 , 2020, Nature.

[78]  A. Walls,et al.  Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein , 2020, Cell.

[79]  Qiang Zhou,et al.  Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2 , 2020, Science.

[80]  B. Graham,et al.  Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation , 2020, Science.

[81]  P. Niu,et al.  Genome Composition and Divergence of the Novel Coronavirus (2019-nCoV) Originating in China , 2020, Cell Host & Microbe.

[82]  E. Holmes,et al.  A new coronavirus associated with human respiratory disease in China , 2020, Nature.

[83]  Kai Zhao,et al.  A pneumonia outbreak associated with a new coronavirus of probable bat origin , 2020, Nature.

[84]  E. Holmes,et al.  Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding , 2020, The Lancet.

[85]  Y. Hu,et al.  Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China , 2020, The Lancet.

[86]  G. Gao,et al.  A Novel Coronavirus from Patients with Pneumonia in China, 2019 , 2020, The New England journal of medicine.

[87]  D. Veesler,et al.  Structural insights into coronavirus entry , 2019, Advances in Virus Research.

[88]  Timothy B Sackton,et al.  Convergent evolution in the genomics era: new insights and directions , 2019, Philosophical Transactions of the Royal Society B.

[89]  James B. Munro,et al.  Direct Visualization of the Conformational Dynamics of Single Influenza Hemagglutinin Trimers , 2018, Cell.

[90]  A. Gorbalenya,et al.  A planarian nidovirus expands the limits of RNA genome size , 2018, bioRxiv.

[91]  Simon R. Harris,et al.  SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments , 2016, bioRxiv.

[92]  Daniel S. Terry,et al.  Single-molecule imaging of non-equilibrium molecular ensembles on the millisecond timescale , 2016, Nature Methods.

[93]  Henry Pinkard,et al.  Advanced methods of microscope control using μManager software. , 2014, Journal of biological methods.

[94]  Marc A Suchard,et al.  Stability-mediated epistasis constrains the evolution of an influenza protein , 2013, eLife.

[95]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[96]  T. Nabeshima,et al.  Discovery of the First Insect Nidovirus, a Missing Evolutionary Link in the Emergence of the Largest RNA Virus Genomes , 2011, PLoS pathogens.

[97]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[98]  Xiaotao Lu,et al.  High Fidelity of Murine Hepatitis Virus Replication Is Decreased in nsp14 Exoribonuclease Mutants , 2007, Journal of Virology.

[99]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[100]  D. J. Kiviet,et al.  Empirical fitness landscapes reveal accessible evolutionary paths , 2007, Nature.

[101]  J. J. Bull,et al.  Theory of Lethal Mutagenesis for Viruses , 2007, Journal of Virology.

[102]  Alexander E Gorbalenya,et al.  Discovery of an RNA virus 3'->5' exoribonuclease that is critically involved in coronavirus RNA synthesis. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[103]  John L. Sullivan,et al.  Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus , 2003, Nature.

[104]  Manfred Eigen,et al.  Error catastrophe and antiviral strategy , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[105]  F Sachs,et al.  A direct optimization approach to hidden Markov modeling for single channel kinetics. , 2000, Biophysical journal.

[106]  J. Qin,et al.  Miscellanea. Combining parametric and empirical likelihoods , 2000 .

[107]  D. Myszka,et al.  Improving biosensor analysis , 1999, Journal of molecular recognition : JMR.

[108]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[109]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[110]  COVID-19 Genomics UK (COG-UK) Consortium Report Number 6 , 2020 .