A Novel SARS-CoV-2 Variant of Concern, B.1.526, Identified in New York

Manuscript Text

While evolution of SARS-CoV-2 was deemed to be slow at the beginning of the global pandemic (1), at least three major variants of concern have emerged over the past two months (2-4). These lineages are each characterized by numerous mutations in the spike protein, raising concerns that they may escape from therapeutic monoclonals and vaccine-induced antibodies. The hallmark mutation of B.1.1.7, the first SARS-CoV-2 variant of concern that emerged in the UK, is N501Y located in the receptor-binding domain (RBD) of spike (2). This variant is seemingly more transmissible and possibly more virulent (5-7). The other variants of concern, B.1.351 (first detected in South Africa) (3) and P.1 (first described in Brazilian travelers) (4), share the N501Y mutation with B.1.1.7 but contain an E484K substitution in RBD (3, 4). Epidemiological evidence suggests that P.1 emerged as part of a second surge in Manaus, Brazil despite a high pre-existing seroprevalence to SARS-CoV-2 in the population. Reinfections with P.1, as well as with another related Brazilian variant P.2 that also harbors E484K, have been documented (8, 9).

Our previous study on B.1.351 demonstrated that this variant is refractory to neutralization by a number of monoclonal antibodies directed to the top of RBD, including several that have received emergency use authorization (10). Moreover, this variant was markedly more resistant to neutralization by convalescent plasma and vaccinee sera. Importantly, these effects were largely mediated by the E484K mutation. These finding are worrisome in light of recent reports that two vaccine trials showed a substantial drop in efficacy in South Africa (11, 12). We therefore began an effort to survey our patient population at the Columbia University Irving Medical Center in New York City for B.1.351 and other E484K variants such as P.1 and P.2.

We first developed rapid PCR-based, single-nucleotide-polymorphism assays to search for N501Y and E484K mutations (see schematic in Fig. S1 and Methods in Supplement) in clinical samples known to be positive for SARS-CoV-2 and stored in the Columbia University Biobank, and patient information was extracted from the COVID-Care database (13). Between November 1, 2020 and February 15, 2021, a total of 60,539 nasopharyngeal swabs underwent clinical testing for SARS-CoV-2 at our medical center, with 4,358 positive samples identified. We screened 1,142 samples randomly chosen from this period for the two signature mutations. A total of 927 samples yielded a signal in our genotyping assays. We found that 83 (9.0%) were positive for E484K and 17 (1.8%) were positive for N501Y. Only one sample contained both mutations. The earliest case with E484K was collected in mid-November 2020. Subsequently, there was a substantial increase in E484K-positive cases over time (Fig. 1A), from 1.3% in early November to 5.3% by mid-January, and ultimately to 12.3% between February 8^th and 15^th. Viruses harboring N501Y also increased over time, from the earliest detection in mid-January to 2.6% of screened isolates by mid-February.

• Download figure
• Open in new tab

Figure 1. Increasing prevalence of E484K from November 2020 to February 2021 including E484K-harboring isolates within the B.1.526 lineage. (A)

The earliest detected E484K variant from our study period was collected in mid-November 2020, and a total of 83 isolates with this substitution were identified between November 2020 and February 15, 2021, out of 927 total samples that showed a signal in our screening assay (upper panel). The prevalence of E484K (n_E484K/(n_WT + n_E484K)) subsequently increased over time, from 1.3% in early November (11/09/2020-11/22/2020) to 11.5% in early February (2/1/2021-2/15/2021), as shown in the lower panel. (B) To place our E484K isolates within the context of globally circulating SARS-CoV-2 strains, we downloaded 140 genomes from GISAID randomly sampled across 50 distinct lineages. Phylogenetic analysis demonstrated that the majority of our B.1.526 isolates are tightly clustered (bottom; magenta), with the remaining E484K-harboring isolates interspersed within the parent B.1 lineage. The B.1.1.7 (UK; purple), B.1.351 (South Africa; green), and P.2 (Brazil; teal) samples identified fall within these branches, respectively, as expected, as shown. Samples from our collection which harbor E484K or N501Y but do not fall within one of these lineages are shown in magenta and purple, respectively.

We then performed whole genome nanopore sequencing on samples flagged as potential N501Y- or E484K-harboring strains (n=65). We also sequenced samples negative for these signature mutations obtained during the same time period, all with Ct values below 35 (n=65). Sequencing results verified the E484K and N501Y substitutions in all samples identified by our screening PCR assays. Based on phylogenetic analyses including publicly available genomes (Fig. 1B), six cases with N501Y were identified as belonging to the B.1.1.7 lineage, two cases with E484K as P.2, and one sample as B.1.351, which harbored both N501Y and E484K based on our screening assay. However, quite unexpectedly, the large majority (n=49) of the remaining cases with E484K fell within a single lineage, B.1.526 (14). It is this novel variant that is surging, alarmingly, in our patient population over the past few weeks.

Nearly all of the newly identified B.1.526 variants have a set of common mutations in the spike protein: L5F, T95I, D253G, E484K, D614G, and A701V. Fig. 2A displays all of the spike mutations found in all variant viruses identified in the study, along with their phylogenetic relationship. Fig. 2B shows that D253G resides in the antigenic supersite within the N-terminal domain (15), which is a target for neutralizing antibodies (16), whereas the E484K is situated at the RBD interface with the cellular receptor ACE2. The A701V mutation near the furin cleavage site is also shared with variant B.1.351. The impact of the E484K mutation on antibody neutralization was assessed using 4 monoclonal antibodies with emergency use authorization, 10 convalescent plasma, and 10 vaccinee sera. As shown in Fig. S2, the neutralizing activity of REGN10987 against E484K pseudovirus is unaltered, but the activities of REGN10933, CB6, and LY-CoV555 are either impaired or abolished. Likewise, neutralizing activities of convalescent plasma or vaccinee sera are lower by 7.7-fold or 3.4-fold, respectively, against the E484K variant. These results signify an important antigenic drift in B.1.526 that could have clinical consequences, as have been noted for B.1.351 (10-12).

• Download figure
• Open in new tab

Figure 2. S-protein amino acid substitutions and structural changes represented in E484K- and N501Y-harboring isolates. (A)

Consensus sequences for 65 genomes were aligned against the Wuhan-Hu-1 reference genome (NC_045512). A phylogenetic reconstruction of these variant genomes is displayed on the left, generated using IQtree with 1,000 bootstrap replicates and maximum-likelihood tree search. Lineage assignments for major clades as determined by Pangolin are shown on the appropriate branches of the phylogenetic tree. On the right, residues at which at least one sample harbored a mutation compared to the reference genome are displayed above the S protein schematic. Amino acid substitutions identified in each sample are highlighted in color and wildtype residues are shown in black; an asterisk (*) indicates a silent nucleotide mutation. (B) Structural changes reflecting substitutions identified in the majority of B.1.526 isolates. Of the amino acid changes largely conserved in the B.1.526 lineage (L5F, T95I, D253G, E484K, D614G, A701V), D253G resides in the antigenic supersite within the N-terminal domain, a target for neutralizing antibodies, E484K at the RBD interface with the cellular receptor ACE2, and A701V near the furin cleavage site. Abbreviations – NTD – N-terminal domain; RBD – receptor binding domain; FP – fusion peptide.

Patients with E484K variant viruses were comparable in gender, race and ethnicity to those with wildtype SARS-CoV-2, but on average were older (58.1 vs 52.4 years, p=0.049) and more likely to present to the ED or be admitted to the hospital (85.9% vs 70.8%, p=0.007, Table S1). There were no differences in the frequency of ICU admissions and length of stay between groups. The majority of patients with E484K variants were geographically concentrated in two distinct neighborhoods in the catchment area of our hospital system, but many others were found scattered throughout the metropolitan area without evidence for a single outbreak (Fig. S3). An early case of our novel B.1.526 variant (NP-3581 in Fig. 2A) was detected in November 2020 in a patient with advanced AIDS who had first presented in August with infection by a SARS-CoV-2 strain (NP-3005 in Fig. 2A) that initially lacked the L5F, D253G, and E484K mutations. This case is reminiscent of reported examples of extensive intra-host evolution due to persistent SARS-CoV-2 infection in immunocompromised individuals (17, 18).

We also investigated SARS-CoV-2 sequences in public databases and found ∼140 genomes highly related to the newly identified B.1.526 variant (14) (Fig. 1B). These were predominantly from samples collected in the Northeastern US, suggesting that E484K in the B.1.526 lineage is now widespread in the region, the original epicenter of COVID-19 in the US (19). This supports concerns that novel variants may spread in regions with a relatively high sero-prevalence. Moreover, it appears that the E484K mutation has emerged in at least 59 different lineages of SARS-CoV-2 (20), a real testament to convergent evolution. In conclusion, we identified B.1.526 as a local lineage of concern due to E484K in particular, which could threaten the efficacy of current antibody therapies and vaccines. This discovery also highlights the need for a concerted national surveillance program to track and contain the spread of novel SARS-CoV-2 variants.