Identification of broad, potent antibodies to functionally constrained regions of SARS-CoV-2 spike following a breakthrough infection

The antiviral benefit of antibodies can be compromised by viral escape especially for rapidly evolving viruses. Therefore, durable, effective antibodies must be both broad and potent to counter newly emerging, diverse strains. Discovery of such antibodies is critically important for SARS-CoV-2 as the global emergence of new variants of concern (VOC) has compromised the efficacy of therapeutic antibodies and vaccines. We describe a collection of broad and potent neutralizing monoclonal antibodies (mAbs) isolated from an individual who experienced a breakthrough infection with the Delta VOC. Four mAbs potently neutralize the Wuhan-Hu-1 vaccine strain, the Delta VOC, and also retain potency against the Omicron VOCs, including recently circulating BA.4/BA.5, in both pseudovirus-based and live virus assays, and one also potently neutralizes SARS-CoV-1. The potency of these mAbs was greater against Omicron VOCs than all but one of the mAbs that had been approved for therapeutic applications. The mAbs target distinct epitopes on the spike glycoprotein, three in the receptor binding domain (RBD) and one in an invariant region downstream of the RBD in subdomain 1 (SD1). The escape pathways we defined at single amino acid resolution with deep mutational scanning show they target conserved, functionally constrained regions of the glycoprotein, suggesting escape could incur a fitness cost. Overall, these mAbs are novel in their breadth across VOCs, their epitope specificity, and include a highly potent mAb targeting a rare epitope outside of the RBD in SD1.

Evolution of the SARS-CoV-2 mutational spectrum

SARS-CoV-2 evolves rapidly in part because of its high mutation rate. Here we examine whether this mutational process itself has changed during viral evolution. To do this, we quantify the relative rates of different types of single nucleotide mutations at four-fold degenerate sites in the viral genome across millions of human SARS-CoV-2 sequences. We find clear shifts in the relative rates of several types of mutations during SARS-CoV-2 evolution. The most striking trend is a roughly two-fold decrease in the relative rate of G[->]T mutations in Omicron versus early clades, as was recently noted by Ruis et al (2022). There is also a decrease in the relative rate of C[->]T mutations in Delta, and other subtle changes in the mutation spectrum along the phylogeny. We speculate that these changes in the mutation spectrum could arise from viral mutations that affect genome replication, packaging, and antagonization of host innate-immune factors--although environmental factors could also play a role. Interestingly, the mutation spectrum of Omicron is more similar than that of earlier SARS-CoV-2 clades to the spectrum that shaped the long-term evolution of sarbecoviruses. Overall, our work shows that the mutation process is itself a dynamic variable during SARS-CoV-2 evolution, and suggests that human SARS-CoV-2 may be trending towards a mutation spectrum more similar to that of other animal sarbecoviruses.

A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here we describe a new deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We demonstrate this new platform by making libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ~7000 distinct amino-acid mutations in the context of up to ~135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ~105 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

Evolution of antibody immunity following Omicron BA.1 breakthrough infection

Understanding the evolution of antibody immunity following heterologous SAR-CoV-2 breakthrough infection will inform the development of next-generation vaccines. Here, we tracked SARS-CoV-2 receptor binding domain (RBD)-specific antibody responses up to six months following Omicron BA.1 breakthrough infection in mRNA-vaccinated individuals. Cross-reactive serum neutralizing antibody and memory B cell (MBC) responses declined by two- to four-fold through the study period. Breakthrough infection elicited minimal de novo Omicron-specific B cell responses but drove affinity maturation of pre-existing cross-reactive MBCs toward BA.1. Public clones dominated the neutralizing antibody response at both early and late time points, and their escape mutation profiles predicted newly emergent Omicron sublineages. The results demonstrate that heterologous SARS-CoV-2 variant exposure drives the evolution of B cell memory and suggest that convergent neutralizing antibody responses continue to shape viral evolution.

Deep mutational scans for ACE2 binding, RBD expression, and antibody escape in the SARS-CoV-2 Omicron BA.1 and BA.2 receptor-binding domains

SARS-CoV-2 continues to acquire mutations in the spike receptor-binding domain (RBD) that impact ACE2 receptor binding, folding stability, and antibody recognition. Deep mutational scanning prospectively characterizes the impacts of mutations on these biochemical properties, enabling rapid assessment of new mutations seen during viral surveillance. However, the effects of mutations can change as the virus evolves, requiring updated deep mutational scans. We determined the impacts of all amino acid mutations in the Omicron BA.1 and BA.2 RBDs on ACE2-binding affinity, RBD folding, and escape from binding by the LY-CoV1404 (bebtelovimab) monoclonal antibody. The effects of some mutations in Omicron RBDs differ from those measured in the ancestral Wuhan-Hu-1 background. These epistatic shifts largely resemble those previously seen in the Beta variant due to the convergent epistatically modifying N501Y substitution. However, Omicron variants show additional lineage-specific shifts, including examples of the epistatic phenomenon of entrenchment that causes the Q498R and N501Y substitutions present in Omicron to be more favorable in that background than in earlier viral strains. In contrast, the Omicron substitution Q493R exhibits no sign of entrenchment, with the derived state, R493, being as unfavorable for ACE2 binding in Omicron RBDs as in Wuhan-Hu-1. Likely for this reason, the R493Q reversion has occurred in Omicron sub-variants including BA.4/BA.5 and BA.2.75, where the affinity buffer from R493Q reversion may potentiate concurrent antigenic change. Consistent with prior studies, we find that Omicron RBDs have reduced expression, and identify candidate stabilizing mutations that ameliorate this deficit. Last, our maps highlight a broadening of the sites of escape from LY-CoV1404 antibody binding in BA.1 and BA.2 compared to the ancestral Wuhan-Hu-1 background. These BA.1 and BA.2 deep mutational scanning datasets identify shifts in the RBD mutational landscape and inform ongoing efforts in viral surveillance. Author SummarySARS-CoV-2 evolves in part through mutations in its spike receptor-binding domain. As these mutations accumulate in evolved variants, they shape the future evolutionary potential of the virus through the phenomenon of epistasis. We characterized the functional impacts of mutations in the Omicron BA.1 and BA.2 receptor-binding domains on ACE2 receptor binding, protein folding, and recognition by the clinical LY-CoV1404 antibody. We then compared the measurements to prior data for earlier variants. These comparisons identify patterns of epistasis that may alter future patterns of Omicron evolution, such as turnover in the availability of specific affinity-enhancing mutations and an expansion in the number of paths of antibody escape from a key monoclonal antibody used for therapeutic treatment of COVID-19. This work informs continued efforts in viral surveillance and forecasting.

The landscape of antibody binding affinity in SARS-CoV-2 Omicron BA.1 evolution

The Omicron BA.1 variant of SARS-CoV-2 escapes convalescent sera and monoclonal antibodies that are effective against earlier strains of the virus. This immune evasion is largely a consequence of mutations in the BA.1 receptor binding domain (RBD), the major antigenic target of SARS-CoV-2. Previous studies have identified several key RBD mutations leading to escape from most antibodies. However, little is known about how these escape mutations interact with each other and with other mutations in the RBD. Here, we systematically map these interactions by measuring the binding affinity of all possible combinations of these 15 RBD mutations (215 = 32,768 genotypes) to four monoclonal antibodies (LY-CoV016, LY-CoV555, REGN10987, and S309) with distinct epitopes. We find that BA.1 can lose affinity to diverse antibodies by acquiring a few large-effect mutations and can reduce affinity to others through several small-effect mutations. However, our results also reveal alternative pathways to antibody escape that do not include every large-effect mutation. Moreover, epistatic interactions are shown to constrain affinity decline in S309 but only modestly shape the affinity landscapes of other antibodies. Together with previous work on the ACE2 affinity landscape, our results suggest that escape of each antibody is mediated by distinct groups of mutations, whose deleterious effects on ACE2 affinity are compensated by another distinct group of mutations (most notably Q498R and N501Y).

Molecular fate-mapping of serum antibodies reveals the effects of antigenic imprinting on repeated immunization

The ability of serum antibody to protect against pathogens arises from the interplay of antigen-specific B cell clones of different affinities and fine specificities. These cellular dynamics are ultimately responsible for serum-level phenomena such as antibody imprinting or "Original Antigenic Sin" (OAS), a proposed propensity of the immune system to rely repeatedly on the first cohort of B cells that responded to a stimulus upon exposure to related antigens. Imprinting/OAS is thought to pose a barrier to vaccination against rapidly evolving viruses such as influenza and SARS-CoV-2. Precise measurement of the extent to which imprinting/OAS inhibits the recruitment of new B cell clones by boosting is challenging because cellular and temporal origins cannot readily be assigned to antibodies in circulation. Thus, the extent to which imprinting/OAS impacts the induction of new responses in various settings remains unclear. To address this, we developed a "molecular fate-mapping" approach in which serum antibodies derived from specific cohorts of B cells can be differentially detected. We show that, upon sequential homologous boosting, the serum antibody response strongly favors reuse of the first cohort of B cell clones over the recruitment of new, naIve-derived B cells. This "primary addiction" decreases as a function of antigenic distance, allowing secondary immunization with divergent influenza virus or SARS-CoV-2 glycoproteins to overcome imprinting/OAS by targeting novel epitopes absent from the priming variant. Our findings have implications for the understanding of imprinting/OAS, and for the design and testing of vaccines aimed at eliciting antibodies to evolving antigens.

Receptor binding domain (RBD) antibodies contribute more to SARS-CoV-2 neutralization when target cells express high levels of ACE2

Neutralization assays are experimental surrogates for the effectiveness of infection- or vaccine-elicited polyclonal antibodies and therapeutic monoclonal antibodies targeting SARS-CoV-2. However, the measured neutralization can depend on details of the experimental assay. Here we systematically assess how ACE2 expression in target cells affects neutralization by antibodies to different spike epitopes in lentivirus pseudovirus neutralization assays. For high ACE2-expressing target cells, receptor binding domain (RBD) antibodies account for nearly all neutralizing activity in polyclonal human sera. But for lower ACE2-expressing target cells, antibodies targeting regions outside the RBD make a larger (although still modest) contribution to serum neutralization. These serum-level results are mirrored for monoclonal antibodies: N-terminal domain (NTD) antibodies and RBD antibodies that do not compete for ACE2 binding incompletely neutralize on high ACE2-expressing target cells, but completely neutralize on cells with lower ACE2 expression. Our results show that ACE2 expression level in the target cells is an important experimental variable, and that high ACE2 expression emphasizes the role of a subset of RBD-directed antibodies.

Compensatory epistasis maintains ACE2 affinity in SARS-CoV-2 Omicron BA.1

The Omicron BA.1 variant emerged in late 2021 and quickly spread across the world. Compared to the ancestral Wuhan Hu-1 strain and other pre-Omicron SARS-CoV-2 variants, BA.1 has many mutations, a number of which are known to enable antibody escape 1-3. Many of these antibody-escape mutations individually decrease the spike receptor-binding domain (RBD) affinity for ACE2 in the background of early SARS-CoV-2 variants 4, but BA.1 still binds ACE2 with high affinity 5,6. The fitness and evolution of the BA.1 lineage is therefore driven by the combined effects of numerous mutations. Here, we systematically map the epistatic interactions between the 15 mutations in the RBD of BA.1 relative to the Wuhan Hu-1 strain. Specifically, we measure the ACE2 affinity of all possible combinations of these 15 mutations (2 15 = 32,768 genotypes), spanning all possible evolutionary intermediates from the ancestral Wuhan Hu-1 strain to BA.1. We find that immune escape mutations in BA.1 individually reduce ACE2 affinity but are compensated by epistatic interactions with other affinity-enhancing mutations, including Q498R and N501Y. Thus, the ability of BA.1 to evade immunity while maintaining ACE2 affinity is contingent on acquiring multiple interacting mutations. Our results implicate compensatory epistasis as a key factor driving substantial evolutionary change for SARS-CoV-2 and are consistent with Omicron BA.1 arising from a chronic infection.

Mosaic RBD nanoparticles protect against multiple sarbecovirus challenges in animal models

To combat future SARS-CoV-2 variants and spillovers of SARS-like betacoronaviruses (sarbecoviruses) threatening global health, we designed mosaic nanoparticles presenting randomly-arranged sarbecovirus spike receptor-binding domains (RBDs) to elicit antibodies against conserved/relatively-occluded, rather than variable/immunodominant/exposed, epitopes. We compared immune responses elicited by mosaic-8 (SARS-CoV-2 and seven animal sarbecoviruses) and homotypic (only SARS-CoV-2) RBD-nanoparticles in mice and macaques, observing stronger responses elicited by mosaic-8 to mismatched (not on nanoparticles) strains including SARS-CoV and animal sarbecoviruses. Mosaic-8 immunization showed equivalent neutralization of SARS-CoV-2 variants including Omicron and protected from SARS-CoV-2 and SARS-CoV challenges, whereas homotypic SARS-CoV-2 immunization protected only from SARS-CoV-2 challenge. Epitope mapping demonstrated increased targeting of conserved epitopes after mosaic-8 immunization. Together, these results suggest mosaic-8 RBD-nanoparticles could protect against SARS-CoV-2 variants and future sarbecovirus spillovers.

The SARS-CoV-2 Delta variant induces an antibody response largely focused on class 1 and 2 antibody epitopes

Exposure histories to SARS-CoV-2 variants and vaccinations will shape the specificity of antibody responses. To understand the specificity of Delta-elicited antibody immunity, we characterize the polyclonal antibody response elicited by primary or mRNA vaccine-breakthrough Delta infections. Both types of infection elicit a neutralizing antibody response focused heavily on the receptor-binding domain (RBD). We use deep mutational scanning to show that mutations to the RBDs class 1 and class 2 epitopes, including sites 417, 478, and 484-486 often reduce binding of these Delta-elicited antibodies. The anti-Delta antibody response is more similar to that elicited by early 2020 viruses than the Beta variant, with mutations to the class 1 and 2, but not class 3 epitopes, having the largest effects on polyclonal antibody binding. In addition, mutations to the class 1 epitope (e.g., K417N) tend to have larger effects on antibody binding and neutralization in the Delta spike than in the D614G spike, both for vaccine- and Delta-infection-elicited antibodies. These results help elucidate how the antigenic impacts of SARS-CoV-2 mutations depend on exposure history.

Shifting mutational constraints in the SARS-CoV-2 receptor-binding domain during viral evolution

SARS-CoV-2 has evolved variants with substitutions in the spike receptor-binding domain (RBD) that impact its affinity for ACE2 receptor and recognition by antibodies. These substitutions could also shape future evolution by modulating the effects of mutations at other sites--a phenomenon called epistasis. To investigate this possibility, we performed deep mutational scans to measure the effects on ACE2 binding of all single amino-acid mutations in the Wuhan-Hu-1, Alpha, Beta, Delta, and Eta variant RBDs. Some substitutions, most prominently N501Y, cause epistatic shifts in the effects of mutations at other sites, thereby shaping subsequent evolutionary change. These epistatic shifts occur despite high conservation of the overall RBD structure. Our data shed light on RBD sequence-function relationships and facilitate interpretation of ongoing SARS-CoV-2 evolution.

Narrow transmission bottlenecks and limited within-host viral diversity during a SARS-CoV-2 outbreak on a fishing boat

The long-term evolution of viruses is ultimately due to viral mutants that arise within infected individuals and transmit to other individuals. Here we use deep sequencing to investigate the transmission of viral genetic variation among individuals during a SARS-CoV-2 outbreak that infected the vast majority of crew members on a fishing boat. We deep-sequenced nasal swabs to characterize the within-host viral population of infected crew members, using experimental duplicates and strict computational filters to ensure accurate variant calling. We find that within-host viral diversity is low in infected crew members. The mutations that did fix in some crew members during the outbreak are not observed at detectable frequencies in any of the sampled crew members in which they are not fixed, suggesting viral evolution involves occasional fixation of low-frequency mutations during transmission rather than persistent maintenance of within-host viral diversity. Overall, our results show that strong transmission bottlenecks dominate viral evolution even during a superspreading event with a very high attack rate.

Dynamics of infection-elicited SARS-CoV-2 antibodies in children over time

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection elicits an antibody response that targets several viral proteins including spike (S) and nucleocapsid (N); S is the major target of neutralizing antibodies. Here, we assess levels of anti-N binding antibodies and anti-S neutralizing antibodies in unvaccinated children compared with unvaccinated older adults following infection. Specifically, we examine neutralization and anti-N binding by sera collected up to 52 weeks following SARS-CoV-2 infection in children and compare these to a cohort of adults, including older adults, most of whom had mild infections that did not require hospitalization. Neutralizing antibody titers were lower in children than adults early after infection, but by 6 months titers were similar between age groups. The neutralizing activity of the childrens sera decreased modestly from one to six months; a pattern that was not significantly different from that observed in adults. However, infection of children induced much lower levels of anti-N antibodies than in adults, and levels of these anti-N antibodies decreased more rapidly in children than in adults, including older adults. These results highlight age-related differences in the antibody responses to SARS-CoV-2 proteins and, as vaccines for children are introduced, may provide comparator data for the longevity of infection-elicited and vaccination-induced neutralizing antibody responses.

Structural changes in the SARS-CoV-2 spike E406W mutant escaping a clinical monoclonal antibody cocktail

The SARS-CoV-2 receptor-binding domain (RBD) E406W mutation abrogates neutralization mediated by the REGEN-CoV therapeutic monoclonal antibody (mAb) COVID-19 cocktail and the cilgavimab (AZD1061) mAb. Here, we show that this residue substitution remodels the ACE2-binding site allosterically, thereby dampening receptor recognition severely and altering the epitopes recognized by these three mAbs. Although vaccine-elicited neutralizing antibody titers are decreased similarly against the E406 mutant and the Delta or Epsilon variants, broadly neutralizing sarbecovirus mAbs, including a clinical mAb, inhibit the E406W spike mutant.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any genomes within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

An antibody-escape calculator for mutations to the SARS-CoV-2 receptor-binding domain

A key goal of SARS-CoV-2 surveillance is to rapidly identify viral variants with mutations that reduce neutralization by polyclonal antibodies elicited by vaccination or infection. Unfortunately, direct experimental characterization of new viral variants lags their sequence-based identification. Here we help address this challenge by aggregating deep mutational scanning data into an "escape calculator" that estimates the antigenic effects of arbitrary combinations of mutations to the viruss spike receptor-binding domain (RBD). The calculator can be used to intuitively visualize how mutations impact polyclonal antibody recognition, and score the expected antigenic effect of combinations of mutations. These scores correlate with neutralization assays performed on SARS-CoV-2 variants, and emphasize the ominous antigenic properties of the recently described Omicron variant. An interactive version of the calculator is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/escape-calc/, and we provide a Python module for batch processing.

Structure, receptor recognition and antigenicity of the human coronavirus CCoV-HuPn-2018 spike glycoprotein

The recent isolation of CCoV-HuPn-2018 from a child respiratory swab indicates that more coronaviruses are spilling over to humans than previously appreciated. Here, we determined cryo-electron microscopy structures of the CCoV-HuPn-2018 spike glycoprotein trimer in two distinct conformational states and identified that it binds canine, feline and porcine aminopeptidase N (APN encoded by ANPEP) orthologs which serve as entry receptors. Introduction of an oligosaccharide at position N739 of human APN renders cells susceptible to CCoV-HuPn-2018 spike-mediated entry, suggesting that single nucleotide polymorphisms could account for the detection of this virus in some individuals. Human polyclonal plasma antibodies elicited by HCoV-229E infection and a porcine coronavirus monoclonal antibody inhibit CCoV-HuPn-2018 S-mediated entry, indicating elicitation of cross-neutralizing activity among -coronaviruses. These data provide a blueprint of the CCoV-HuPn-2018 infection machinery, unveil the viral entry receptor and pave the way for vaccine and therapeutic development targeting this zoonotic pathogen.

Antibody-mediated broad sarbecovirus neutralization through ACE2 molecular mimicry

Understanding broadly neutralizing sarbecovirus antibody responses is key to developing countermeasures effective against SARS-CoV-2 variants and future spillovers of other sarbecoviruses. Here we describe the isolation and characterization of a human monoclonal antibody, designated S2K146, broadly neutralizing viruses belonging to all three sarbecovirus clades known to utilize ACE2 as entry receptor and protecting therapeutically against SARS-CoV-2 beta challenge in hamsters. Structural and functional studies show that most of the S2K146 epitope residues are shared with the ACE2 binding site and that the antibody inhibits receptor attachment competitively. Viral passaging experiments underscore an unusually high barrier for emergence of escape mutants making it an ideal candidate for clinical development. These findings unveil a key site of vulnerability for the development of a next generation of vaccines eliciting broad sarbecovirus immunity.

A SARS-CoV-2 variant elicits an antibody response with a shifted immunodominance hierarchy

Many SARS-CoV-2 variants have mutations at key sites targeted by antibodies. However, it is unknown if antibodies elicited by infection with these variants target the same or different regions of the viral spike as antibodies elicited by earlier viral isolates. Here we compare the specificities of polyclonal antibodies produced by humans infected with early 2020 isolates versus the B.1.351 variant of concern (also known as Beta or 20H/501Y.V2), which contains mutations in multiple key spike epitopes. The serum neutralizing activity of antibodies elicited by infection with both early 2020 viruses and B.1.351 is heavily focused on the spike receptor-binding domain (RBD). However, within the RBD, B.1.351-elicited antibodies are more focused on the "class 3" epitope spanning sites 443 to 452, and neutralization by these antibodies is notably less affected by mutations at residue 484. Our results show that SARS-CoV-2 variants can elicit polyclonal antibodies with different immunodominance hierarchies.