medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Emergence in late 2020 of multiple lineages of SARS-CoV-2 Spike protein variants affecting amino acid position 677

Emma B. Hodcroft1,*, Daryl B. Domman2,*,a, Kasopefoluwa Oguntuyo5, Daniel J. Snyder3, Maarten Van Diest4, Kenneth H. Densmore4, Kurt C. Schwalm2, Jon Femling2, Jennifer L. Carroll4, Rona S. Scott4, Martha M. Whyte6, Michael D. Edwards7, Noah C. Hull8, Christopher G. Kevil4, John A. Vanchiere4, Benhur Lee5, Darrell L. Dinwiddie2,†, Vaughn S. Cooper9,†,a, Jeremy P. Kamil4†,a

*. these authors contributed equally a. corresponding authors †. senior authors

1 Institute of Social and Preventive Medicine, University of Bern, Switzerland: 2 University of New Mexico Health Sciences Center, Albuquerque, NM, USA. 3 Microbial Genome Sequencing Center, LLC, Pittsburgh, PA, USA 4 Louisiana State University Health Sciences Center, Shreveport, Shreveport, LA, USA. 5 Icahn Mt Sinai School of Medicine, New York, New York. USA 6 Louisiana Department of Health. Minden, LA, USA. 7 New Mexico Department of Health, Albuquerque, NM, USA 8 Wyoming Public Health Laboratory, Cheyenne, WY, USA 9 University of Pittsburgh, School of Medicine, Pittsburgh, PA, USA.

NOTE: This preprint reports new research that has not been certified by peer review and should not be used to guide clinical practice. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Abstract.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S) plays critical roles in host cell entry. Non-synonymous substitutions affecting S are not uncommon and have become fixed in a number of SARS-CoV-2 lineages. A subset of such enable escape from neutralizing antibodies or are thought to enhance transmission through mechanisms such as increased affinity for the cell entry receptor, ACE2. Independent genomic surveillance programs based in New Mexico and Louisiana contemporaneously detected the rapid rise of numerous clade 20G (lineage B.1.2) infections carrying a Q677P substitution in S. The variant was first detected in the US on October 23, yet between 01 Dec 2020 and 19 Jan 2021 it rose to represent 27.8% and 11.3% of all SARS-CoV-2 genomes sequenced from Louisiana and New Mexico, respectively. Q677P cases have been detected predominantly in the south central and southwest United States; as of 03 Feb 2021, GISAID data show 499 viral sequences of this variant from the USA. Phylogenetic analyses revealed the independent evolution and spread of at least six distinct Q677H sub-lineages, with first collection dates ranging from mid August to late November, 2020. Four 677H clades from clade 20G (B.1.2) , 20A (B.1.234), and 20B (B.1.1.220, and B.1.1.222) each contain roughly 100 or fewer sequenced cases, while a distinct pair of clade 20G clusters are represented by 754 and 298 cases, respectively. Although sampling bias and founder effects may have contributed to the rise of S:677 polymorphic variants, the proximity of this position to the polybasic cleavage site at the S1/S2 boundary are consistent with its potential functional relevance during cell entry, suggesting parallel evolution of a trait that may confer an advantage in spread or transmission. Taken together, our findings demonstrate simultaneous convergent evolution, thus providing an impetus to further evaluate S:677 polymorphisms for effects on proteolytic processing, cell tropism, and transmissibility. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Introduction. In mid December 2020, the reported a SARS-CoV-2 variant termed B.1.1.7 (20I/501Y.V1) that exhibited a rapid increase in its range and incidence following its initial detection in November (Andrew Rambaut, Nick Loman, Oliver Pybus, Wendy Barclay, Jeff Barrett, Alesandro Carabelli, Tom Connor, Tom Peacock, David L Robertson, Erik Volz, COVID-19 Consortium UK (CoG-UK), 2020; Volz, Mishra, et al., 2021). Since then, additional “variants of concern” have emerged, namely lineages B.1.351 (20H/501Y.V2) from S Africa (Tegally et al., 2020) and P.1 (20J/501Y.V3) and P.2 from Brazil (Voloch et al., 2020; Faria et al., 2021; Naveca, da Costa, et al., 2021; Sabino et al., 2021). A key concern is that certain polymorphisms may enhance SARS-CoV-2 infectivity or transmission, akin to what was seen for Spike D614G (Korber et al., 2020; Volz, Hill, et al., 2021), which has overtaken the original D614 form of the virus that dominated at the outset of the .

In areas where SARS-CoV-2 seroprevalence is high due to elevated rates of transmission during primary waves of the pandemic, selection pressures may have favored the emergence of variants that escape neutralizing antibodies. Such circumstances are thought to have contributed to the emergence of lineages B.1.351 and P.1 (501Y.V2 and 501Y.V3), which in addition to Spike N501Y, harbor at least two other non-synonymous substitutions, K417N/T and E484K, which have been found to confer escape from neutralizing antibodies (Weisblum, Schmidt, Zhang, DaSilva, Poston, J. C. C. Lorenzi, et al., 2020; Cele et al., 2021; Liu et al., 2021; Wang et al., 2021).

Despite accounting for roughly 25% of globally recorded COVID-19 cases, and approximately 20% of the available SARS-CoV-2 genome data, relatively few studies have detailed the introduction and emergence of SARS-CoV-2 lineages in the United States (Rochman et al., 2020; Worobey et al., 2020; Washington et al., 2021; Zeller et al., 2021). Furthermore, seroprevalence surveys indicate that between 1 in 10 and 1 in 3 people in the USA have already been infected with SARS-CoV-2, potentially high enough that selection for immune evasion may be present (Angulo, Finelli and Swerdlow, 2021; Bubar et al., 2021).

Variants affecting the Spike (S) protein are of great interest due to their potential to impact transmissibility and potential for effects on natural or -induced immunity. One early in Spike, D614G, quickly came to dominate the pandemic (currently accounting for >98% of sequences), at least in part because it promotes an S conformation that is more competent for binding to ACE2 (Yurkovetskiy et al., 2020) and medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

reduces shedding of the S1 subunit that contains the receptor binding domain (Zhang et al., 2020). Missense mutations at other positions, for example S477N and N439K, have appeared multiple times in large infection clusters in Australia and Europe, and are associated with resistance to certain antibodies and/or increased affinity for ACE2 (Hodcroft et al., 2020; Liu et al., 2020; Weisblum, Schmidt, Zhang, DaSilva, Poston, J. C. Lorenzi, et al., 2020; Gaebler et al., 2021; Thomson et al., 2021).

Here, we describe evidence from two independent SARS-CoV-2 genomic surveillance programs in Louisiana and New Mexico that each detected the rise of S variants affecting position 677 in the later months of 2020. We further provide phylogenetic analyses that identify six independent Q677H sub-lineages and one Q677P sub-lineage that all appear to have emerged within the United States. These variants were not detected until mid-August 2020, but as of 03 Feb 2021 already account for over 2,327 of the 102,462 genomes deposited to GISAID from the USA. Given the broad detection of the lineages across multiple states and the apparent increase in frequency of detection, these novel emergent Q677H and Q677P lineages merit further study for potential differences in transmissibility.

Results. In late January of 2021, our two independent SARS-CoV-2 genomic surveillance programs, based at the University of New Mexico Health Sciences (UNM HSC) in Albuquerque, New Mexico and the Louisiana State University Health Sciences Center (LSUHS) in Shreveport, Louisiana, each noticed increasing numbers of PANGO lineage B.1.2 (Rambaut et al., 2020) / Nextstrain clade 20G viruses carrying an S:Q677P mutation, and that this variant had increased in frequency in samples collected in late 2020 to mid January (FIG. 1A). We also noted broad but uneven geographical distribution of these S: 677 polymorphic viruses across the United States, with certain states approaching 15% of total viral sequences submitted since the beginning of the pandemic. Given the proximity of this residue to the polybasic cleavage site (681-685), and our observations prompted communication between our two surveillance programs, and motivated us to carry out phylogenetic analyses of variation at S position 677.

Site frequency dynamics In addition to the Q677P sub-lineage, our analyses indicated that SARS-CoV-2 variants carrying non-synonymous mutations affecting S codon 677 have arisen at least six other times in the United States (FIG. 2A, FIG. 3). Four of these occurred within clade 20G, a large clade that has accounted for over 30% of US sequences since October 2020, and as of the beginning of February 2021 represents approximately 50% of all SARS-CoV-2 sequence data in the USA (Nextstrain Build, North America, Feb 04, 2021, 2021). The amino acid Q changes to H due to a mutation at nucleotide position medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

23593. Notably, in four of these six lineages, the mutation changes from G to U, whereas in the other two, it changes from G to C (FIG 2B). In contrast, the S:Q677P variant occurs by virtue of an A to C change at position 23592. All mutations leading to Q677H or Q677P involve transversions. Hence, their spontaneous occurrence is generally disfavored relative to transitions. In SARS-CoV-2 samples from human infections, A to C and G to C transversions occur at only ~10% the frequency of C to U transitions, while G to U mutations are more common, occurring half as frequently as C to U. (Wright, Lakdawala and Cooper, 2020) (Ratcliff and Simmonds, 2021). For ease of discussion and to avoid using geographically-associated names or nicknames, we have named each of the seven S: 677 mutant lineages identified here after American bird species.

Six U.S. lineages and a provisional naming system. The largest of the 677 variant sub-lineages (“Robin 1”) is a B.1.2 / 20G clade virus carrying Q677H that first appears in GISAID data from a sample with a August 17, 2020 collection date. As of Feb 4, 2021, this sub-lineage contained 754 sequences (Table 1, FIG. 3). Robin 1 is found in over 30 US states, but predominates in the Midwest. A second Q677H clade, distinguished from Robin 1 by an N2361K substitution in orf1a, first appeared from a Oct 6, 2020 sample from Alabama and is named “Robin 2” owing to its similarity to the parental Robin 1 sub-lineage. This cluster contains 303 sequences, and is found mostly in the Southeast. The next largest cluster is the Q677P variant of 20G (B.1.2) (“Pelican”), which was first detected in Oregon from a sample with collection date of Oct 23, 2020 and as of Feb 3, 2021 contains 504 sequences. The Q677P variant has been detected in LA, NM, NC, WY, MA, ID, MI AZ, CA, TX, WI, and MD, and five international sequences (Australia (2), Denmark, Switzerland, India)(Emma B. Hodcroft, 2021). The remaining Q677H sub-lineages each contain around 100 or fewer sequences, and are named: Yellowhammer, detected mostly in the southeast US; Bluebird, mostly in the northeast United States; Quail, mainly in the Southwest and Northeast; and Mockingbird, mainly in the South-central and East coast states (Table 1, Fig. 3). A schematic summarizing the key lineage-specific and shared protein polymorphisms of the US S: Q677P and S:Q677H variants is shown in FIG 4.

Epidemiological details Broad-scale genomic surveillance efforts conducted by the New Mexico Department of Health (NM DOH) and the UNM HSC in December and January, 2021 revealed that 83 of 733 SARS-CoV-2 genomes sequenced in New Mexico between December 1st and January 19th contained the S:Q677P variant at a general frequency of 11.3% in the SARS-CoV-2 positive individuals. In New Mexico, the S:Q677P substitution was first observed from a December 12, 2020 sample, and the frequency of the variant in medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

genomic sequences increased through January. In December 2020, our New Mexico- based surveillance effort detected 23 genomes harboring the Q677P mutation, with an additional 59 detected in January 2021. However, limited genomic sequencing of positive cases in New Mexico in October- November, hinder the ability to accurately determine the true rate of increase. Similarly, genomic surveillance efforts conducted by LSUHS together with the Louisiana Department of Health detected a first occurrence of the S:Q677P variant on Dec 1, 2020, which over the month rose to 187 complete genomes, many from a single, large congregate facility outbreak. In part due to deep sampling of the large outbreak, the Q677P virus amounted to 37.2% of all viral genomes sequenced from Dec 2020 samples collected in Louisiana (LSUHS submissions account for 481 out of the 503, or 95.6% of the total sequenced viruses available for the state in Dec on GISAID). However, this Q677P variant also amounted to 14% of all LSUHS samples collected in Jan 2021, and 11.4% of all Louisiana samples collected between Jan 01 - Jan 19 (LSUHS sequence data for 2021 collection dates is currently available up to Jan 19). An S:Q677H variant first detected from Louisiana in the summer of 2021, appearing in LSUHS samples collected on July 21, 2020 and Aug 11, 2020. Although only 2 additional S:Q667H sequences were collected from the entire state of Louisiana in Nov, both on the 27th. Strikingly, however, the total for Dec 2020 rose to 53, amounting to 10.5% of statewide sequence data from Dec 2021. According to the most recent LSUHS data, which covers collection dates up to Jan 19, 2021, the Q677H polymorphism occurs in 5.7% of total SARS-CoV-2 genome data. When considered in unison, S: Q677P and S: Q677H samples together comprise 47.5% of all Dec 2020 viral genomes collected from Louisiana, and 17.1% of the genome sequences collected in the state from Jan 01 - 19, 2021. Albeit that founder effects lead to stochastic fluctuations in the abundance of even neutral genotypes, they also contribute to the expansion of the fittest viruses. Thus, the observation of a sudden and contemporaneous increase in the abundance of S: Q677H and Q677P variants by surveillance programs in two different states along the southeast / southwest corridor is remarkable.

Modeled structure of variants The S1/S2 cleavage site contains the multibasic cleavage site and is found in a disordered region within Spike (S)(Wrapp et al., 2020). We used SWISS MODEL to model this inherently flexible region and spotlight the Q677P residue (FIG. 5). Although the Q677P position is outside the furin binding pocket (Tian, Huajun and Wu, 2012), we speculate that the presence of a proline at this site may introduce a favorable kink that promotes the dynamic, conformational changes necessary for cleavage at the S1/S2 junction, which is governed not only by furin-like activities, but also by trypsin-like proteases and cathepsins(Jaimes, Millet and Whittaker, 2020). Moreover, the medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

introduction of a proline in this model appears to be 3.7 angstroms away from the carbon backbone of S689 (relative to 4.9 angstroms for the native glutamine), which may promote atomic interactions that encourage conformations favorable for proteolytic cleavage. In the case of the S: Q677H substitution, histidine protonation could similarly act as a conformational switch affecting accessibility to proteases. Cleavage at the S1/S2 boundary promotes a more ‘open’ S conformation that is more competent to bind ACE2(Wrobel et al., 2020), these putative mechanisms for enhanced cleavage at the S1/S2 junction may promote more efficient viral entry (Hoffmann, Kleine-Weber and Pöhlmann, 2020).

Discussion. Caveats. Selectively neutral mutations can become fixed in a lineage purely by chance and human behaviour. For instance, the 20E (EU1) lineage characterized by an S: A222V polymorphism emerged suddenly in Europe over the summer, but has not been found to show any evidence for increased transmissibility (Hodcroft et al., 2020) and instead is thought to have been spread via holiday travel and relaxing summertime restrictions. Additional S variants such as N439K and S477N also rapidly increased in frequency in Europe over the summer and into the fall of 2020. Although S477N reportedly increases affinity for the entry receptor, ACE2 (Zahradník et al., 2021), and both mutations may impact antibody neutralization to some degree (Starr et al., 2020; Weisblum, Schmidt, Zhang, DaSilva, Poston, J. C. Lorenzi, et al., 2020; Liu et al., 2021; Thomson et al., 2021), neither shows any signature of increased transmissibility over the S: D614G background from which they emerged, and neither have become prominent in the United States. Therefore, SARS-CoV-2 variants can emerge and increase greatly in number over time in the absence of any clear or sustained selective advantage.

Convergent evolution is a hallmark of positive selection. The repeated evolution of a trait in independent populations provides strong evidence of adaptation. Between August and November, 2020, seven independent lineages of SARS-CoV-2 with S:Q677H or S:Q677P mutations arose and gained in frequency. This coincidental rise and spread on independent genetic backgrounds is remarkable and suggests some fitness advantage. Observed frequencies undoubtedly incorporate some sampling biases, which may over-estimate the relative amount of S variants affecting position 677. However, sampling bias cuts both ways. U.S. states with fewer deposited sequences may simply have missed detection of these variants. To date, all 677H/P mutants collected from the US stem from the S:614G lineage that now predominates worldwide, but alongside varied polymorphisms in S, N, ORF1a, ORF1b, and other genes, suggesting any fitness advantage of S:677 mutations is largely medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

independent of these other mutations (Table 1, FIG. 2-4). Nonetheless, the relatively slow rise of lineages with S:677 substitutions may suggest that any fitness benefits are modest relative to other circulating variants, which may have also independently gained adaptive mutations. Given their relatively recent emergence, however, Q677P/H lineages may continue to rise as a percentage of total cases. Additional laboratory and genetic surveillance data will be needed to define whether S: 677 polymorphisms are biologically relevant. Although we focus here on the appearance and expansion of S: 677H/P mutants in the United States, global analyses reveal that 677H mutants have arisen multiple times elsewhere in the world as well, including large clusters of 677H mutants in Egypt and Denmark, and multiple clusters in India (Emma B. Hodcroft, 2021). Furthermore, a newly designated, emergent PANGO lineage, B.1.525, carries S: Q677H in addition to several mutations seen in B.1.1.7 (501Y.V1), such as S: del 69-70 and S: del 144, and also, S: E484K (Rambaut et al., 2020; Ainé O’toole Verity Hill, 2021). Remarkably, a 19B cluster harboring the ostensibly less fit, ‘ancestral’ D614 Spike, which has been circulating at ≤2% of global frequency since August 2020, recently resurfaced as a newly re-emergent lineage carrying N501Y together with 677H (Wagner, 2021). N501Y is notably found in all three of the ‘variants of concern’(Tegally et al., 2020; Faria et al., 2021; Lauring and Hodcroft, 2021; Naveca, Nascimento, et al., 2021; Volz, Mishra, et al., 2021). Albeit that this observation is circumstantial, it further suggests that the 677H mutation may confer an evolutionary advantage to SARS-CoV-2.

The importance of sequencing for viral surveillance and genetic epidemiology. Global surveillance of genomic changes in SARS-CoV-2 varies widely, with leading countries such as Australia, New Zealand, the United Kingdom, and Denmark sequencing viruses from 5-50% of all cases and lagging countries such as the United States, France, Spain, and Brazil sequencing less than 1% of all cases. It is notable that these Q677 variants were detected in the undersampled US population, suggesting that these lineages may actually be more prevalent. The finding of lineages containing 677H mutations in better sampled countries like Denmark indicates that they repeatedly emerge but may be outcompeted by lineages with larger gains in transmission like B.1.1.7 (Davies et al., 2020). Collectively these findings demonstrate the value of greater genomic sequencing and the importance of tracking the emergence and spread of lineages that combine multiple mutations which could enhance transmissibility or evade immunity from prior infection or .

At least two emergent lineages of concern, B.1.1.7 (501Y.V1), and a newer variant whose prevalence is on the rise in Uganda both contain amino acid substitutions medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

affecting the first position of the polybasic cleavage site, S:P681H and S:P681R, respectively (Andrew Rambaut, Nick Loman, Oliver Pybus, Wendy Barclay, Jeff Barrett, Alesandro Carabelli, Tom Connor, Tom Peacock, David L Robertson, Erik Volz, COVID- 19 Genomics Consortium UK (CoG-UK), 2020; Lule Bugembe et al., 2021). The polybasic site strongly impacts viral replication in culture and as well as pathogenesis in animal models (Hoffmann, Kleine-Weber and Pöhlmann, 2020; Johnson et al., 2021). Although it is too early to predict whether any particular S: 677 polymorphic lineages will persist, given these observations, the recurrent parallelism affecting S: 677 suggests that this position will continue to surface in variants that show signs of increased transmissibility or fitness. It will thus be critical to not only to continue genomic surveillance of SARS-CoV-2 to monitor the prevalence of such variants over time, but also to formally define any biological characteristics of these polymorphisms in cell culture and small animal model systems.

Methods.

Sequencing. Genome sequencing of SARS-CoV-2 at the UNM HSC was conducted from RNA isolated from residual clinical specimens that had previously been determined to be positive for SARS-CoV-2 by FDA-approved diagnostic testing. SARS-CoV-2 genomic sequences were amplified by PCR using the widely adopted ARTIC primer set (v3) and were sequenced either on an Illumina MiSeq or Oxford Nanopore GridION system using in-house, slightly modified versions of protocols (https://www.protocols.io/view/sars-cov- 2-illumina-miseq-protocol-v-1-bjd9ki96 & https://www.protocols.io/view/ncov-2019- sequencing-protocol-v3-locost-bh42j8ye), respectively.

For LSUHS samples, SARS-CoV-2 genomic RNA sequencing was done as follows. 13 µL of de-identified total RNA from patient anterior turbinate nasal swabs previously determined as SARS-CoV-2 positive CDC N1 / N2 RT-PCR assay results of Ct <26, was subjected to hybridization capture enrichment sequencing. For each sample, 13µL of extracted RNA was reverse transcribed using Maxima H-minus ds cDNA kits (ThermoFisher Scientific). Libraries were enriched using a Nextera Flex for Enrichment Library Preparation kit with a Respiratory Virus Oligo Set v1, with samples being pooled in 12-plex enrichment reactions, essentially as described elsewhere (O’Flaherty et al., 2018). The resulting pools were quantified and grouped in sets of no more than 48 samples and run on a NextSeq 550 using a 150cyc High Output Flow Cell. We used breseq (Deatherage and Barrick, 2014) (v.0.34.1) to map reads to Wuhan-Hu-1 SARS- CoV-2 (NC_045512) or 2019-nCoV WIV04 (GISAID EPI_ISL_402124, NCBI Genbank MN996528)(Zhou et al., 2020) and call the consensus sequence. All predicted mutations were reported for isolates exceeding mean 40x coverage. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Phylogenetic analyses. A ‘focal’ Nextstrain build was prepared as described in Hodcroft et al, 2020, initially selecting for sequences with any mutation at position 23592 or 23593. These sequences are used as a ‘focal set’ around which background sequences are selected, capturing both the sequences most genetically similar to the focal set, as well as a selection of sequences distributed across time and by country. Both of these are then processed through the open-source Nextstrain `ncov` pipeline to produce a time- resolved phylogenetic analysis. A time-stamped version of the Nextstrain build used in this manuscript can be found here: https://nextstrain.org/groups/neherlab/ncov/S.Q677/2020-02-04, and the most recently updated version of this build can be viewed here: https://nextstrain.org/groups/neherlab/ncov/S.Q677. The final designations of the clusters and counts of sequences within them were determined from the resulting phylogenetic tree rather than the ‘raw counts’ from on the nucleotide mutations given above, as the phylogenetic structure can better account for gaps and reversions in sequences which might prevent a sequence being picked up by mutations alone. The resulting JSON file from the Nextstrain pipeline was used to visualize the phylogenetic trees using the baltic package (https://github.com/evogytis/baltic). A full table of acknowledgements for the data included in this analysis is available in the supplement.

Structural analyses. SWISS-MODEL (Bienert et al., 2017; Waterhouse et al., 2018) was used to model the full length, wild type and Q677P mutant spike glycoprotein. The template utilized for the model was 7BBH(Wrobel, A.G., Benton, D.J., Rosenthal, P.B., Gamblin, S.J., 2020), which has all receptor binding domains in the down conformation.the entire furin cleavage site for the D614G ("WT") and D614G+Q677P (mutant) version of the SARS-CoV-2 Spike, based on PDB . Then PyMol 2.4 (Schödinger, Inc.) was subsequently used to superimpose the resulting structures and highlight the side chains of each individual structure.

Acknowledgements. We gratefully acknowledge all submitting authors and collecting authors on whose work this research is based, and to all researchers, clinicians, and public health authorities who make SARS-CoV-2 sequence data available in a timely manner via the GISAID initiative (Elbe and Buckland-Merrett, 2017; Shu and McCauley, 2017). This work was supported by a COVID-19 Fast Grants award from Emergent Ventures, an initiative of the Mercatus Center at George Mason University (J.P.K.), and by an intramural grant and other funding from the Office of the Vice Chancellor for Research at LSU Health Sciences Center Shreveport (J.P.K., R.S.S., J.A.V.); an medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Institutional Development Award from the National Institutes of General Medical Sciences of the NIH under grant number P20 GM121307 (C.G.K.); by the Swiss National Science Foundation through grant number 31CA30196046 (E.B.H.), by the U.S. National Center for Research Resources and the National Center for Advancing Translational Sciences of the National Institutes of Health through Grant Number UL1TR001449 (D.L.D., D.B.D), a KL2 Mentored Career Development Award KL2R001448 (D.B.D), and a Translational and Clinical Pilot Project Award CTSC008-11 (D.B.D). We would like to thank the UNM Center for Advanced Research Computing, supported in part by the National Science Foundation, for providing the high performance computing resources used in this work.

IRB Approvals. SARS-CoV-2 genome sequences were generated under IRB approved protocols STUDY00001445 (LSU Health Sciences Center), and 14-039 and 20-151 (University of New Mexico Health Sciences Center).

Conflicts of Interest. V.S.C. and D.J.S. are co-founders of Microbial Genome Sequencing Center, LLC. All other authors declare no conflicts of interest.

Author contributions. Collected samples: J.A.V., M.D.E., N.C.H.; Sample Preparation: R.S.S., J.L.C.; Generation and processing of sequence data: D.L.D, D.J.S., V.S.C.; Phylogenetic Analyses: E.B.H., D.B.D., V.S.C.; Structural modeling: K.O.; B.L.; Data Curation: K.H.D., C.G.K., M.v.D., D.L.D.; Analysis and Interpretation of Data: E.B.H., D.B.D., D.L.D., M.v.D., K.H.D., C.G.K, K.O., B.L., V.S.C., J.P.K.. Study Design: E.B.H., D.L.D., D.B.D., V.S.C., J.P.K.; Wrote the paper E.B.H., V.S.C., J.P.K., with comments from all authors. Obtained Funding: D.L.D., D.B.D., C.G.K., J.P.K. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

FIGURE LEGENDS.

Figure 1. Rising prevalence of SARS-CoV-2 S:Q677P and S:Q677H variants in the United States. (A) Monthly number of Q677H and Q677P variant viruses in the USA and worldwide found in GISAID data. Note that data for samples collected in December and January is incomplete for many surveillance facilities. (B) State-by-state prevalence of United States Q677H and Q677P variants (percent frequency shaded from 0%- 15% out of total GISAID data up to Feb 03, 2021).

Figure 2. Newly emergent SARS-CoV-2 spike position 677 variants. (A) A simplified Nextstrain phylogeny of 677H and 677P lineages. (B) A nucleotide and protein alignment of representative S genes from Q677P 92C, and Q677H 93T, Q677H 93C lineages, set against the reference sequence: Dec 2019 Wuhan, China, WIV04 (GISAID: EPI_ISL_402124, GenBank: MN996528). The polybasic or “furin” cleavage site is labeled, as are the portions of the S1 and S2 subunits pictured in the alignment.

Figure 3. Time-scaled phylogeny of recently expanded lineages within the USA harboring mutations at spike position 677. A ‘focal’ Nextstrain build was prepared initially selecting for sequences with any mutation at position 23592 or 23593. These sequences were used as a ‘focal set’ around which background sequences are selected, capturing both the sequences most genetically similar to the focal set, as well as a selection of sequences distributed across time and by country. Both of these were processed through the Nextstrain `ncov` pipeline to produce a time-resolved phylogenetic analysis. The resulting JSON file from the Nextstrain pipeline was used to visualize the phylogenetic trees using the baltic package (https://github.com/evogytis/baltic). A full table of acknowledgements for the data included in this analysis is available in the supplement.

Figure 4. Key lineage-specific and shared polymorphisms of US S:Q677P and S:677H variants. The U.S. S:Q677P and S:Q677H variants defined in Figs.1-2 are illustrated for the major defining protein variants that are shared among most or all of the viruses (pink), or restricted to one (green) or two (yellow). Q677P and Q677H polymorphisms are shown in bold and in red. A full table of acknowledgements for the data included in this analysis is available in the supplement.

Figure 5. Structure of SARS-CoV-2 Spike protein denoting location of Q677P within a disordered loop adjacent to the polybasic (furin) cleavage site. Structure was modeled from PDB: 7BBH using SWISS-MODEL and visualized using PyMol. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Table 1. Genetic attributes and relative abundance of six different SARS-CoV-2 lineages identified by S:Q677 mutations.

Supplementary Tables 1, 2, 3, 4, 5, and 6. Acknowledgement of the Authors and the Originating laboratories where the clinical specimens or virus isolates were first obtained and the Submitting laboratories, where sequence data have been generated and submitted to GISAID. We gratefully acknowledge the Authors from the Originating laboratories responsible for obtaining the specimens and the Submitting laboratories where genetic sequence data were generated and shared via the GISAID Initiative, on which this research is based.

medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

REFERENCES.

Ainé O’toole Verity Hill (2021) Proposal for new lineage within B.1 #4: B.1.525, cov- lineages / pango-designation. Available at: https://github.com/cov-lineages/pango- designation/issues/4 (Accessed: 12 February 2021).

Andrew Rambaut, Nick Loman, Oliver Pybus, Wendy Barclay, Jeff Barrett, Alesandro Carabelli, Tom Connor, Tom Peacock, David L Robertson, Erik Volz, COVID-19 Genomics Consortium UK (CoG-UK) (2020) Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations, Virological.org. Available at: 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 (Accessed: 12 February 2021).

Angulo, F. J., Finelli, L. and Swerdlow, D. L. (2021) ‘Estimation of US SARS-CoV-2 Infections, Symptomatic Infections, Hospitalizations, and Deaths Using Seroprevalence Surveys’, JAMA network open, 4(1), p. e2033706. doi: 10.1001/jamanetworkopen.2020.33706.

Bienert, S. et al. (2017) ‘The SWISS-MODEL Repository-new features and functionality’, Nucleic acids research, 45(D1), pp. D313–D319. doi: 10.1093/nar/gkw1132.

Bubar, K. M. et al. (2021) ‘Model-informed COVID-19 vaccine prioritization strategies by age and serostatus’, Science. doi: 10.1126/science.abe6959.

Cele, S. et al. (2021) ‘Escape of SARS-CoV-2 501Y.V2 variants from neutralization by convalescent plasma’, bioRxiv. medRxiv. doi: 10.1101/2021.01.26.21250224.

Davies, N. G. et al. (2020) ‘Estimated transmissibility and severity of novel SARS-CoV-2 Variant of Concern 202012/01 in England’, bioRxiv. medRxiv. doi: 10.1101/2020.12.24.20248822.

Deatherage, D. E. and Barrick, J. E. (2014) ‘Identification of mutations in laboratory- evolved microbes from next-generation sequencing data using breseq’, Methods in molecular biology , 1151, pp. 165–188. doi: 10.1007/978-1-4939-0554-6_12.

Elbe, S. and Buckland-Merrett, G. (2017) ‘Data, disease and diplomacy: GISAID’s innovative contribution to global health’, Global challenges (Hoboken, NJ), 1(1), pp. 33– 46. doi: 10.1002/gch2.1018.

Emma B. Hodcroft, R. N. (2021) Phylogenetic analysis of SARS-CoV-2 clusters in their international context - cluster S.Q677. Available at: =S.677P,677H&m=div&p=grid&r=country&tl=clade_membership">https://nextstrain.org/ groups/neherlab/ncov/S.Q677?c=gt- medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

nuc_23593>=S.677P,677H&m=div&p=grid&r=country&tl=clade_membership.

Faria, N. R. et al. (2021) Genomic characterisation of an emergent SARS-CoV-2 lineage in Manaus: preliminary findings. Available at: https://virological.org/t/genomic- characterisation-of-an-emergent-sars-cov-2-lineage-in-manaus-preliminary- findings/586.

Gaebler, C. et al. (2021) ‘Evolution of antibody immunity to SARS-CoV-2’, Nature, pp. 1–10. doi: 10.1038/s41586-021-03207-w.

Hodcroft, E. B. et al. (2020) ‘Emergence and spread of a SARS-CoV-2 variant through Europe in the summer of 2020’, medRxiv : the preprint server for health sciences. doi: 10.1101/2020.10.25.20219063.

Hoffmann, M., Kleine-Weber, H. and Pöhlmann, S. (2020) ‘A Multibasic Cleavage Site in the Spike Protein of SARS-CoV-2 Is Essential for Infection of Human Lung Cells’, Molecular cell. doi: 10.1016/j.molcel.2020.04.022.

Jaimes, J., Millet, J. and Whittaker, G. (2020) ‘Proteolytic Cleavage of the SARS-CoV-2 Spike Protein and the Role of the Novel S1/S2 Site’, SSRN, p. 3581359. doi: 10.2139/ssrn.3581359.

Johnson, B. A. et al. (2021) ‘Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis’, Nature. doi: 10.1038/s41586-021-03237-4.

Korber, B. et al. (2020) ‘Spike mutation pipeline reveals the emergence of a more transmissible form of SARS-CoV-2’, bioRxiv. doi: 10.1101/2020.04.29.069054.

Lauring, A. S. and Hodcroft, E. B. (2021) ‘Genetic Variants of SARS-CoV-2—What Do They Mean?’, JAMA: the journal of the American Medical Association, 325(6), pp. 529– 531. doi: 10.1001/jama.2020.27124.

Liu, Z. et al. (2020) ‘Landscape Analysis of Escape Variants Identifies SARS-CoV-2 Spike Mutations that Attenuate Monoclonal and Serum Antibody Neutralization’, Cell host & microbe, In Press. doi: 10.2139/ssrn.3725763.

Liu, Z. et al. (2021) ‘Identification of SARS-CoV-2 spike mutations that attenuate monoclonal and serum antibody neutralization’, Cell host & microbe, 0(0). doi: 10.1016/j.chom.2021.01.014.

Lule Bugembe, D. et al. (2021) ‘A SARS-CoV-2 lineage A variant (A.23.1) with altered spike has emerged and is dominating the current Uganda epidemic’, medRxiv. medRxiv. doi: 10.1101/2021.02.08.21251393.

Naveca, F., Nascimento, V., et al. (2021) ‘Phylogenetic relationship of SARS-CoV-2 sequences from Amazonas with emerging Brazilian variants harboring mutations E484K and N501Y in the Spike protein’, Virological. org. Available at: https://virological.org/t/phylogenetic-relationship-of-sars-cov-2-sequences-from- medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

amazonas-with-emerging-brazilian-variants-harboring-mutations-e484k-and-n501y-in- the-spike-protein/585.

Naveca, F., da Costa, C., et al. (2021) ‘SARS-CoV-2 reinfection by the new Variant of Concern (VOC) P. 1 in Amazonas, Brazil’, virological. org. Preprint available at: https://virological. org/t/sars-cov-2-reinfection-by-thenew-variant-of-concern-voc-p-1-in-a mazonas-brazil/596. Available at: https://virological.org/t/sars-cov-2-reinfection-by-the- new-variant-of-concern-voc-p-1-in-amazonas-brazil/596.

Nextstrain Build, North America, Feb 04, 2021 (2021). doi: https://nextstrain.org/ncov/north-america/2021-02- 04?d=frequencies&f_country=USA&p=full.

O’Flaherty, B. M. et al. (2018) ‘Comprehensive viral enrichment enables sensitive respiratory virus genomic identification and analysis by next generation sequencing’, Genome research, 28(6), pp. 869–877. doi: 10.1101/gr.226316.117.

Rambaut, A. et al. (2020) ‘A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology’, Nature microbiology, 5(11), pp. 1403–1407. doi: 10.1038/s41564-020-0770-5.

Ratcliff, J. and Simmonds, P. (2021) ‘Potential APOBEC-mediated RNA editing of the genomes of SARS-CoV-2 and other coronaviruses and its impact on their longer term evolution’, . doi: 10.1016/j.virol.2020.12.018.

Rochman, N. D. et al. (2020) ‘Ongoing Global and Regional Adaptive Evolution of SARS-CoV-2’, Cold Spring Harbor Laboratory. doi: 10.1101/2020.10.12.336644.

Sabino, E. C. et al. (2021) ‘Resurgence of COVID-19 in Manaus, Brazil, despite high seroprevalence’, The Lancet. doi: 10.1016/S0140-6736(21)00183-5.

Shu, Y. and McCauley, J. (2017) ‘GISAID: Global initiative on sharing all data - from vision to reality’, Euro surveillance: bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin, 22(13). doi: 10.2807/1560- 7917.ES.2017.22.13.30494.

Starr, T. N. et al. (2020) ‘Prospective mapping of viral mutations that escape antibodies used to treat COVID-19’, Cold Spring Harbor Laboratory. doi: 10.1101/2020.11.30.405472.

Tegally, H. et al. (2020) ‘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’, bioRxiv. medRxiv. doi: 10.1101/2020.12.21.20248640.

Thomson, E. C. et al. (2021) ‘Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity’, Cell. doi: 10.1016/j.cell.2021.01.037.

Tian, S., Huajun, W. and Wu, J. (2012) ‘Computational prediction of furin cleavage sites medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

by a hybrid method and understanding mechanism underlying diseases’, Scientific reports, 2, p. 261. doi: 10.1038/srep00261.

Voloch, C. M. et al. (2020) ‘Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil’, medRxiv. Available at: https://www.medrxiv.org/content/10.1101/2020.12.23.20248598v1.abstract.

Volz, E., Hill, V., et al. (2021) ‘Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity’, Cell, 184(1), pp. 64–75.e11. doi: 10.1016/j.cell.2020.11.020.

Volz, E., Mishra, S., et al. (2021) ‘Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data’, bioRxiv. medRxiv. doi: 10.1101/2020.12.30.20249034.

Wagner, C. (2021) SARS-CoV-2 genomic analysis: 19B. Available at: https://nextstrain.org/groups/blab/ncov/19B?c=gt-S_677&m=div.

Wang, Z. et al. (2021) ‘mRNA vaccine-elicited antibodies to SARS-CoV-2 and circulating variants’, Cold Spring Harbor Laboratory. bioRxiv. doi: 10.1101/2021.01.15.426911.

Washington, N. L. et al. (2021) ‘Genomic epidemiology identifies emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States’, medRxiv, p. 2021.02.06.21251159. doi: 10.1101/2021.02.06.21251159.

Waterhouse, A. et al. (2018) ‘SWISS-MODEL: homology modelling of protein structures and complexes’, Nucleic acids research, 46(W1), pp. W296–W303. doi: 10.1093/nar/gky427.

Weisblum, Y., Schmidt, F., Zhang, F., DaSilva, J., Poston, D., Lorenzi, J. C. C., et al. (2020) ‘Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants’, bioRxiv : the preprint server for biology, p. 2020.07.21.214759. doi: 10.1101/2020.07.21.214759.

Weisblum, Y., Schmidt, F., Zhang, F., DaSilva, J., Poston, D., Lorenzi, J. C., et al. (2020) ‘Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants’, eLife, 9. doi: 10.7554/eLife.61312.

Worobey, M. et al. (2020) ‘The emergence of SARS-CoV-2 in Europe and North America’, Science. doi: 10.1126/science.abc8169.

Wrapp, D. et al. (2020) ‘Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation’, Science, 367(6483), pp. 1260–1263. doi: 10.1126/science.abb2507.

Wrobel, A. G. et al. (2020) ‘SARS-CoV-2 and bat RaTG13 spike glycoprotein structures inform on virus evolution and furin-cleavage effects’, Nature structural & molecular biology, 27(8), pp. 763–767. doi: 10.1038/s41594-020-0468-7. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Wrobel, A.G., Benton, D.J., Rosenthal, P.B., Gamblin, S.J. (2020) Structure of Coronavirus Spike from Smuggled Guangdong Pangolin, PDB: 7bbh. Available at: https://www.wwpdb.org/pdb?id=pdb_00007bbh.

Yurkovetskiy, L. et al. (2020) ‘Structural and functional analysis of the D614G SARS- CoV-2 spike protein variant’, Cell, 183(3), pp. 739–751.e8. doi: 10.1016/j.cell.2020.09.032.

Zahradník, J. et al. (2021) ‘SARS-CoV-2 RBD in vitro evolution follows contagious mutation spread, yet generates an able infection inhibitor’, Cold Spring Harbor Laboratory. doi: 10.1101/2021.01.06.425392.

Zeller, M. et al. (2021) ‘Emergence of an early SARS-CoV-2 epidemic in the United States’, medRxiv, p. 2021.02.05.21251235. doi: 10.1101/2021.02.05.21251235.

Zhang, L. et al. (2020) ‘SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity’, Nature communications, 11(1), p. 6013. doi: 10.1038/s41467-020-19808-4.

Zhou, P. et al. (2020) ‘A pneumonia outbreak associated with a new coronavirus of probable bat origin’, Nature, 579(7798), pp. 270–273. doi: 10.1038/s41586-020-2012-7.

medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Figure 1. Rising prevalence of SARS-CoV-2 S:Q677P and Q677H variants in the United States. (A) Monthly number of Q677H and Q677P variant viruses in the USA and worldwide found in GISAID data. Note that data for samples collected in December and January is incomplete for many surveillance facilities. (B) State by state prevalence of United States Q677H and Q677P variant frequency (out of total GISAID data up to Feb 03, 2021). A full table of acknowledgements for the data included in this analysis is available in the supplement.

medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Figure 2. Newly emergent SARS-CoV-2 spike position 677 variants. (A) A simplified Nextstrain phylogeny of 677H and 677P lineages. (B) A nucleotide and protein alignment of representative S genes from Q677P 92C, and Q677H 93T, Q677H 93C lineages, set against the reference sequence: Dec 2019 Wuhan, China, WIV04 (GISAID: EPI_ISL_402124, GenBank: MN996528). The polybasic or “furin” cleavage site is labeled, as are the portions of the S1 and S2 subunits pictured in the alignment. A full table of acknowledgements for the data included in this analysis is available in the supplement.

medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Figure 3. Time-scaled phylogeny of recently expanded lineages within the USA harboring mutations at spike position 677. A ‘focal’ Nextstrain build was prepared initially selecting for sequences with any mutation at position 23592 or 23593. These sequences were used as a ‘focal set’ around which background sequences are selected, capturing both the sequences most genetically similar to the focal set, as well as a selection of sequences distributed across time and by country. Both of these were processed through the Nextstrain `ncov` pipeline to produce a time-resolved phylogenetic analysis. The resulting JSON file from the Nextstrain pipeline was used to visualize the phylogenetic trees using the baltic package (https://github.com/evogytis/baltic). A full table of acknowledgements for the data included in this analysis is available in the supplement.

medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Figure 4. Lineage-specific and shared polymorphisms of US S:Q677P and S:677H variants. The U.S. S:Q677P and S:Q677H variants defined in Figs.1-2 are illustrated for the major defining protein variants that are shared among most or all of the viruses (pink), or restricted to one (green) or two (yellow). Q677P and Q677H polymorphisms are shown in bold and in red. A full table of acknowledgements for the data included in this analysis is available in the supplement. medRxiv preprint doi: https://doi.org/10.1101/2021.02.12.21251658; this version posted February 14, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Figure 5. Structure of SARS-CoV-2 Spike protein denoting location of Q677P within a disordered loop adjacent to the polybasic (furin) cleavage site. Structure was modeled from PDB: 7BBH using SWISS-MODEL and visualized using PyMol.

medRxiv preprint (which wasnotcertifiedbypeerreview)

667H 677H 677P 677H 677H 677H 677H “Robin1” “Robin2” “Pelican” “Yellowhammer “Bluebird” “Quail” “Mockingbird” ” doi:

Defining G23593T G23593T A23592C G23593C G23593T G23593C G23593T https://doi.org/10.1101/2021.02.12.21251658 mutation Lineage 20G / B.1.2 20G / B.1.2 20G / B.1.2 20G / B.1.2 20B / 20B / B.1.1.222 20A / B.1.234 (Nextstrain / B.1.1.220

pangolin) It ismadeavailableundera Approximate 2020-08-17 2020-10-06 2020-10-23 2020-11-24 2020-08-17 2020-10-07 2020-11-26 date of origin istheauthor/funder,whohasgrantedmedRxivalicensetodisplaypreprintinperpetuity. Deposited 754 298 504 125 122 123 52 sequences S D614G D614G D614G L5F T29I D614G G142S Q677H Q677H Q677P D614G T572I Q677H E180V Q677H D614G T732S** D614G CC-BY-NC-ND 4.0Internationallicense

Q677H Q677H ; this versionpostedFebruary14,2021. D936N ORF3a Q57H Q57H Q57H Q57H S58N G172V G172V G172V G172V H78Y ORF1a T265I T265I T265I T265I V665I M2606I N2361K L3352F M2606I K2059R L3352F M2606I Q3729R L3352F D2980G L3352F A3454V N P67S P67S P67S P67S R203K R203K S194L P199L* P199L* P199L* P199L* G204R G204R T391I D377Y N:M210I A376V M234I .

(many) E378Q The copyrightholderforthispreprint ORF1b P314L P314L P314L P314L P314L P314L P314L N1653D T1555I (many) N1653D N1653D R2613C N1653D R2613C P1666L R2613C R2613C Example USA/WI-UW- USA/AL- USA/TX-HMH- USA/AL-HGSC- USA/CT-Yale- USA/NY- USA/TX-HMH- sequence near 2237/2020 Cullman- MCoV- JFBB/2020 352/2020 MSHSPSP- MCoV- clade root JEMT/2020 19011/2020 PV19649/2020 18986/2020 medRxiv preprint (which wasnotcertifiedbypeerreview) * shared on recent common ancestor within 20G / B.1.2 ** T732A occurs on a branch prior to the MRCA of Quail, then A732S occurs on a branch inside of Quail, with only 2 sequences inside the cluster not carrying this mutation. doi: https://doi.org/10.1101/2021.02.12.21251658 It ismadeavailableundera istheauthor/funder,whohasgrantedmedRxivalicensetodisplaypreprintinperpetuity. CC-BY-NC-ND 4.0Internationallicense ; this versionpostedFebruary14,2021. . The copyrightholderforthispreprint