Epistasis lowers the genetic barrier to SARS-CoV-2 neutralizing antibody escape

Waves of SARS-CoV-2 infection have resulted from the emergence of viral variants with neutralizing antibody resistance mutations. Simultaneously, repeated antigen exposure has generated affinity matured B cells, producing broadly neutralizing receptor binding domain (RBD)-specific antibodies with activity against emergent variants. To determine how SARS-CoV-2 might escape these antibodies, we subjected chimeric viruses encoding spike proteins from ancestral, BA.1 or BA.2 variants to selection by 40 broadly neutralizing antibodies. We identify numerous examples of epistasis, whereby in vitro selected and naturally occurring substitutions in RBD epitopes that do not confer antibody resistance in the Wuhan-Hu-1 spike, do so in BA.1 or BA.2 spikes. As few as 2 or 3 of these substitutions in the BA.5 spike, confer resistance to nearly all of the 40 broadly neutralizing antibodies, and substantial resistance to plasma from most individuals. Thus, epistasis facilitates the acquisition of resistance to antibodies that remained effective against early omicron variants.

contributes to the residual effectiveness of ancestral variant-based vaccines against infection and particularly against serious disease. A key question, therefore, is whether and how SARS-CoV-2 might evolve to evade these evolved, broadly neutralizing antibodies.
In this work, we applied selection pressure, using a collection of 40 such antibodies to ancestral Wuhan-Hu-1, as well as BA.1 or BA.2 spike sequences, to generate resistant variants. These analyses reveal frequent epistasis between pre-existing substitutions in BA.1/ BA.2 and resistance mutations acquired during selection experiments, some of which correspond to substitutions appearing in BA.2 daughter lineages. Resistance to broadly neutralizing antibodies is more readily acquired by BA.1/BA.2 than ancestral lineages -single substitutions that confer resistance to more than half of the antibodies in the context of BA.1 or BA.2 fail to do so in the context of Wuhan-Hu-1. Moreover, only two or three of these substitutions introduced into the context of BA.5 are sufficient to confer resistance to nearly all broadly neutralizing antibodies tested. Our findings suggest that epistatic interaction between pre-existing and emergent substitutions in neutralizing epitopes could substantially affect the direction of future SARS-CoV-2 evolution, with implications for vaccine immunogen effectiveness and design.

Selection of SARS-CoV-2 spike variants with antibodies
We chose 40 antibodies, isolated from the memory B cells of multiple volunteer cohorts with varying exposures to SARS-CoV-2 antigens, using an ancestral Wuhan-Hu-1 RBD (Supplementary Data 1 and 2). Specifically, the antibodies were from (i) participants infected with Wuhan-Hu-1 at 6.2 or 12 months previously, some of whom had also received 1 or 2 mRNA vaccine doses prior to antibody isolation 13,17 (ii) recipients of 3 mRNA vaccine doses who had not been infected (antibodies isolated at~1 m after the third dose) 16 or (iii) individuals who had experienced an omicron BA.1 breakthrough infection after 3 vaccine doses (antibodies isolated at~1 m after infection) 21,22 . We chose 'broadly' neutralizing antibodies, defined as those that could neutralize both Wuhan-Hu-1 and omicron (BA.1) pseudotyped viruses, with IC 50 values of~100 ng/ml or less. Most of these antibodies (35 of 40) also neutralized BA.2. With the exception of two class 1/4 broadly neutralizing antibodies, isolated from infected individuals at 1.3 m post-infection 9,23 , the antibodies were chosen without regard to any property (e.g., V H , V L or RBD epitope class, Supplementary Data 2) other than potency and breadth. Based on competition experiments with prototype class 1 through 4 antibodies, the broadly neutralizing collection targeted diverse RBD epitopes that included all classes, but exhibited some skew toward class 3.
To select antibody escape mutants, we employed replicationcompetent recombinant vesicular stomatitis viruses (rVSV/SARS-CoV-2) 24 that encoded the spike proteins from ancestral Wuhan-Hu-1, BA.1 or BA.2 variants in place of the VSV G protein. Diversified rVSV/SARS-CoV-2 populations containing 10 6 infectious units were incubated with antibodies at 1 µg/ml before infection of target cells. Following 2 passages in the presence of each antibody, the acquisition of RBD substitutions in the selected progeny viral populations was evaluated by nextgeneration sequencing (Supplementary Data 3). There was no apparent difference in the ability of antibodies from different classes to select mutations. Indeed, in a total of 120 selection experiments, 39 out of the 40 antibodies tested yielded substitutions, present at a frequency of >10% of sequences in the selected viral populations. Overall, these substitutions were at 34 different positions in the RBD (Fig. 1a, b, Supplementary Fig. 1). Substituted positions in the selected viral populations were spatially proximal to, or within, the ACE2 binding site and were largely consistent with the selecting antibody class designation, determined by competition with prototype class 1-4 antibodies (Supplementary Data 2, Supplementary Fig. 1). With only 3 exceptions, two of which were reversion of BA.1 specific substitution (F375S and S496G), the antibody selection experiments in the BA.1 and BA.2 context generated substitutions at positions that were not already changed in BA.1 or BA.2 compared to Wuhan-Hu-1 (Fig. 1a, b, Supplementary Fig. 1).

Broadly neutralizing antibody resistance mutations
We chose 23 substitutions that were enriched during the passage of rVSV/SARS-CoV-2 variants in the presence of the 40 antibodies and constructed spike expression plasmids based on the context in which the substitutions were selected. If a substitution was enriched by antibody selection in the ancestral rVSV/SARS-CoV-2 Wuhan-Hu-1 context, a corresponding mutant Wuhan-Hu-1 spike expression plasmid was constructed. In instances where a position was substituted in either of the antibody selected rVSV/SARS-CoV-2 BA.1 or rVSV/SARS-CoV-2 BA.2 populations, the corresponding position was substituted in all three parental spike expression plasmids. Thus, we generated 56 pseudotyped HIV-1 viruses with mutant spike sequences and tested their ability to resist neutralization by each of the 40 broadly neutralizing antibodies, at a concentration of 1 µg/ml (Fig. 2, Supplementary Data 4). Infection was quantified relative to uninhibited virus (absence of antibody), and antibody 'escape' was defined as a fivefold increase in infection relative to the antibody-inhibited parental pseudotype, provided infection reached >10% of the uninhibited control. Most of the substitutions conferred specific resistance to one or more antibodies of the same class, consistent with their position on the RBD surface. For example, class 4 or 1/4 broadly neutralizing antibodies, which bind epitopes that are well conserved in sarbecoviruses 9 were escaped by mutations at P384, R408, G413, or G504 (Fig. 3, Supplementary Data 4, Supplementary Fig. 1). Notably, some class 1/4 or 4 broadly neutralizing antibodies were ineffective or less potent against BA.2, which has substitutions that encroach on class 4 epitopes ( Fig. 2 and 3). While most substitutions resulted in class-specific antibody resistance, two substitutions positioned at the base of the RBD either gave modest class-independent reductions in sensitivity (F375S), or no reductions in antibody sensitivity (Y369F). We speculate that these represent fitnessenhancing substitutions that arose during the rVSV/SARS-CoV-2 passage. In the case of Y369F, which arose in selection experiments with diverse class 1, 1/2, 2/3, or 3 antibodies, there was an increase in sensitivity to some antibodies (Fig. 2), especially those in class 4, that is ordinarily concealed when the RBD is in the 'down' configuration. Y369 falls within the epitope of class 4 antibodies, and it is possible that this substitution sensitizes the spike to these antibodies by directly improving antibody binding affinity, or by increasing exposure to the class 4 antibody epitope.
Similar epistatic effects were evident for 14 of the 15 antibody resistance substitutions that were tested in all three pseudotype backgrounds (Fig. 2, Supplementary Figs.

3-5, Supplementary Data 4).
Substitutions at positions N354, F375, D420, L441, N460, and F490 conferred resistance to one or more antibodies when introduced into BA.1 or BA.2 but did not do so when introduced into Wuhan-Hu-1. For 11 antibodies, single-amino-acid substitutions (at positions T345, R346, K444, Y449, N450, L452, F486 or N487) conferred resistance in a context-independent manner. However, the very same substitutions conferred resistance to 10 other antibodies when introduced into BA.1 and/or BA.2 but not when introduced into the Wuhan-Hu-1. Overall, there was extensive epistatic interaction between antibody resistance substitutions and variation in the BA.1 or BA.2 spike proteins (Fig. 2, Supplementary Data 4).
Included in the set of broadly neutralizing antibodies was C099 (class 1), for which a structure of the Fab:RBD interface has been determined 14 . While we failed to generate single-amino-acid resistance mutations in the context of Wuhan-Hu-1, C099 resistance was readily generated by single-amino-acid substitutions in the BA.1 and BA.2 contexts (Fig. 2, Supplementary Fig. 3, Supplementary Data 4). Indeed, two substitutions (at D420 and N460) within or proximal to the target epitope ( Fig. 5a) were necessary to confer resistance to C099 in the Wuhan-Hu-1 context 14 , but each of these substitutions, as well as F486S or N487D, could individually lead to full C099 resistance in BA.1 and/or BA.2 (Fig. 2, Fig. 5a, Supplementary Fig. 3). Similarly, single target epitope substitutions at positions R346, K444 and L441 conferred resistance to C032 14 , a clonal ancestor of the affinity-matured class 3 broadly neutralizing antibody, C080, but none conferred resistance to C080 in the context of Wuhan-Hu-1 (Fig. 2, Fig. 5b, Supplementary  Fig. 4). Again, however, these single substitutions could generate C080 resistance in the BA.1 or BA.2 background. Overall, we conclude that BA.1 and BA.2 variants more readily escape broadly neutralizing antibodies such as C099 and C080 because they already possess a subset of the multiple target epitope substitutions that are required for antibody resistance (Fig. 5a, b).

Epistatic interactions for mutations selected in vitro and naturally in BA.2 daughter lineages
Recently, SARS-CoV-2 sublineages derived from BA.2 have displaced BA.1 and BA.  (Fig. 6a) 25,26 . These substitutions are identical to, or at the same positions as, substitutions that we found to exhibit epistatic interactions with BA.2 substitutions to confer resistance to multiple broadly neutralizing antibodies (Fig. 2, Fig. 4, Supplementary Figs. [3][4][5]. In particular, we found that RBD substitutions at positions that changed during the BA.2 to BA.5 transition (L452 and F486) conferred resistance to antibodies of multiple classes (class 2/3 and 3 for L452R and classes 1, 1/2, and 1/4 for F486S). Notably, for 7 of 14 antibodies, resistance conferred by L452 or F486 substitutions was dependent on the prior acquisition of substitutions in BA.2 (Fig. 6b). Similarly, the N460K substitution that arose during the BA.2 to BA. transition conferred resistance to 6 class 1 or 1/4 antibodies with VH 3-53 and VH 3-66 heavy chains (Fig. 6c), that have previously been shown to generate antibodies that can neutralize BA.1 through BA.5 omicron variants 28 . Notably, the N460K substitution gave resistance to all 6 of the VH 3-53/VH 3-66 antibodies in the BA.2 background, but did not affect sensitivity to the same antibodies in the Wuhan-Hu-1 background (Fig. 6c). Cumulatively, substitutions at positions R346, G446, and N460 that arose during the BA.2 to BA.2.75 to BA.2.75.2 transitions, conferred resistance to 20 broadly neutralizing antibodies. However, substitutions at these positions individually conferred resistance to only four antibodies in the Wuhan-Hu-1 background (Fig. 6c). Overall, therefore, substitutions that occurred naturally in emergent SARS-CoV-2 variants exhibit substantial overlap with those selected by antibodies in vitro and exhibit clear context-dependent effects on antibody sensitivity.

Resistance of synthetic SARS-CoV-2 variants to broadly neutralizing antibodies
Synthetic SARS-CoV-2 'polymutant' spike proteins, encoding combinations of changes arising during in vitro selection experiments, can have emergent escape characteristics similar to those of natural SARS-CoV-2 variants 3,7 . We generated synthetic RBD variants with combinations of mutations selected in vitro to assess the potential effect of selective pressure by broadly neutralizing antibodies. Specifically, we compiled RBD mutations that conferred resistance to non-overlapping subsets of antibodies (Fig. 2). For instance, the R346W substitution conferred resistance to several class 3 antibodies (e.g., C3102, C3108, C1068), the D420N substitution enabled escape from several class 1/4 antibodies (e.g., C099, C1002, C3095, C3111), while the K444N substitution enabled escape from various antibodies, particularly in class 2/3 (e.g., B1011, B1012, B1034). In some cases, two or more alternative substitutions at spatially proximal positions (e.g., T345A and R346W or D420N and N460K) conferred resistance to similar sets of antibodies (Fig. 2). In these instances, only one substitution (R346W and D420N) was chosen. We tested 12 synthetic RBD variants, with 2 to 8 substitutions introduced in the context of a BA.5 spike (BA.5(+2) through BA.5(+8), Fig. 7a) for neutralization by the antibody panel. The synthetic variants all generated infectious pseudotypes, and all resisted neutralization by 32 or more of the 40 antibodies, across all classes (Fig. 7b). Notably, the BA.5(+2) and BA.5(+3a/b/c) variants, with only two (R346W, D420N) and three (R346W and D420N with K444N, V445E or P499S) substitutions, were resistant to 36 to 40 of the 40 broadly neutralizing antibodies. For one of the synthetic variants (BA.5(+3a)) we introduced the same three substitutions (R346W, D420N, and K44N) in the Wuhan-Hu-1 as well as the BA.5 contexts (Fig. 7b.) In the BA.5 context these three substitutions conferred resistance to 10 additional antibodies that they did not confer resistance to in Wuhan-Hu-1 context. Overall, the 40 broadly neutralizing antibodies inhibited   Wuhan-Hu-1, BA.1, and BA.2 pseudotype infection by a median of 170-, 40-and 118-fold at 1 μg/ml, respectively (Fig. 7c). For BA.5, the median inhibition was 6-fold, while for each of the synthetic variants, median inhibition was only 1.0-to 1.2-fold. We conclude that compared to early omicron variants, none of the broadly neutralizing RBD antibodies examined herein remain effective against potential future SARS-CoV-2 variants derived from BA.5 with a limited number of additional substitutions.

Plasma neutralization of synthetic SARS-CoV-2 variants
Populations of antibodies in the memory B-cell compartment may differ from those represented in plasma 29 and in contrast to the monoclonal antibody panel used herein, are not selected for particular properties. We, therefore, selected three synthetic variants (BA.5(+3a), BA.5(+7b), and BA.5(+8)) and determined their neutralization sensitivity to plasma antibodies from individuals with heterogeneous antigen exposures (Fig. 8, Supplementary Data 1, 5). Individuals who were infected with Wuhan-Hu-1-like viruses (but not vaccinated) exhibited low plasma NT 50 against all omicron variants, as did those whose only antigen exposure was two doses of an mRNA vaccine. These groups had median NT 50 titers of 96-257 for BA.1 and BA.2 variants and close to background NT 50 values for BA.5 and the synthetic BA.5-based variants (Fig. 8).
Plasma from three groups who had received multiple antigen exposures, (i) 3x vaccinated, (ii) infected then vaccinated, or (ii) 3x vaccinated and breakthrough infection-the same groups from which the broadly neutralizing antibodies were obtained (Supplementary Data 1)-had higher neutralizing titers against BA.1, BA.2 and BA.5 (NT 50 values of 2231-18584). The synthetic variants were less sensitive to neutralization by these plasmas than either of the BA.2 or BA.5 parental variants (Fig. 8). For example, plasmas from 3x vaccinated individuals, had median NT 50 values against BA.5(+3a) that were sixfold and fourfold reduced compared to BA.2 and BA.5 respectively. Notably, the group with BA.1 breakthrough infection following 3x vaccination had higher NT 50 values against BA.1 (3.4-fold) and BA.2 (5.7-fold) than those without breakthrough infection (Fig. 8). However, these apparent gains in neutralizing activity after BA.1 breakthrough infection were absent against the synthetic variants. Consequently, the BA.5 (+3a) variant was 23-fold and sevenfold less sensitive than BA.2 and BA.5 to plasmas from the BA.1 breakthrough infection group (Fig. 8). We conclude that a significant portion of the residual BA.5 plasma neutralizing activity in these groups can be evaded by the acquisition of a few additional substitutions that epistatically interact with pre-existing substitutions.
Notably, plasma from individuals who had been infected early in the pandemic and subsequently vaccinated (and were thus exposed to ancestral Wuhan-Hu-1-like antigens) 17 had broader neutralization properties than the 3x vaccinated and BA.1 breakthrough groups (Fig. 8). Indeed, for the infected then vaccinated group, the decrease in median NT 50 for BA.5 compared to BA.2 was not statistically significant, and the synthetic variants were only twofold to 3.6-fold less sensitive than BA.5 to these plasmas. These decrements reached significance only for the BA.

Discussion
Prolonged or repeated exposure to antigen drives somatic mutation of antibody lineages resulting in greater diversity and affinity and, in the context of SARS-CoV-2 infection or vaccination, greater neutralization breadth and resilience to escape mutations [13][14][15][16][17][18] . A key consequence of antibody maturation is that SARS-CoV-2 escape from neutralization by individual RBD-specific antibodies often requires more than one substitution in the target epitope 14 . The mechanism underlying this effect is likely based on antibody affinity: as SARS-CoV-2 neutralizing antibodies evolve, they bind with higher affinity to their target epitope 14 . Thereafter, single-amino-acid substitutions may cause a reduction in binding affinity, which is insufficient to cause loss of neutralization activity for antibodies that bind with exceptionally high affinity 14 . A corollary, described herein, is that epistatic interactions between newly acquired and pre-existing viral substitutions enable escape from neutralizing antibodies that have a sufficiently high binding affinity that they tolerate single but not multiple target epitope substitutions. This feature is evident in antibody escape by BA.1 and BA.2 that encode many affinity-reducing mutations in class 1-4 RBD epitopes, which provide a background whereby a small number of additional RBD substitutions can result in additional affinity reductions that cause loss of neutralization for most of the broadly neutralizing antibodies tested.
Based on this premise, we posit that the rapid sequential displacement of early omicron variants by subsequent derivatives that are derived from BA.2, rather than from ancestral lineages, maybe in part because individual substitutions have greater impact on neutralization escape in omicron than in ancestral contexts, due to epistasis. Such a scenario is consistent with a transition to a SARS-CoV-2 phylogeny in which contemporary variants are derived from recent ancestors rather than a more distant ancestral root, which would typically characterize viral phylogeny in immunologically naïve populations 30 . Previously, studies with influenza virus have shown that one of two H3 hemaglutinin substitutions analyzed yielded context-dependent changes in neutralization sensitivity 31 , suggesting that epistatic interactions between existing and emergent substitutions may generally drive neutralizing antibody escape in immunologically drifting viral populations.
Limitations of this study include the fact that the broadly neutralizing antibodies were isolated using a Wuhan-Hu-1 RBD bait to sort memory B cells, and selected for breadth using a BA.1-based neutralization assay. It is possible that antibodies selected using a different bait or selected for different neutralizing properties would have   4401  257  172  89  50  50  96  14510  5478  2850  2231  668  505  1035  27105  8444  10810  8641  4274  4035  2419  36319  18584  16351  5164  648  719  1623  1086  96  159  67  53  51  yielded different results. Moreover, neutralizing antibodies that target other domains of the SARS-CoV-2 spike protein were not examined. However, the large decrement in polyclonal plasma neutralizing activity caused by the installation of a small number of RBD substitutions suggests that the activity of a large fraction of individual, affinitymatured, neutralizing antibody lineages are compromised by epistatic interactions between pre-existing and emergent RBD substitutions described herein. Analysis of antibodies in plasma and in the Wuhan-Hu-1-selected memory B-cell compartment suggests that BA.1 breakthrough infection selectively boosts the subset of vaccine-elicited neutralizing antibodies residing in memory that neutralize BA.1 21,22,32,33 . While some of these memory antibodies also neutralized BA.5, most did not neutralize synthetic variants with only three additional substitutions. BA.1 infection after three vaccine doses thus provides no increase in plasma neutralization over three vaccine doses alone against such variants. In contrast, boosting of a diverse set of evolved antibodies, elicited by ancestral Wuhan-Hu-1 infection, using a Wuhan-Hu-1-based mRNA vaccine, increases and further matures, a broader set of antibodies 17 . Subsets of these antibodies may cross-react with distinct variants, as suggested by the ability of plasma from infected and then vaccinated individuals to neutralize the synthetic variants tested herein. Our findings thus suggest that differences in antigen exposure may generate significant heterogeneity in antibody-neutralizing breadth and potency. Differences in plasma-neutralizing potency against contemporary variants between groups with different antigen exposures are increasingly evident as the antigenic distance between currently circulating and ancestral SARS-CoV-2 variants grows. Overall, our studies indicate that epistatic interactions between pre-existing and new viral mutations, along with population heterogeneity in neutralizing antibody responses, will likely continue to impact the future direction and pace of SARS-CoV-2 variant emergence, with implications for vaccine effectiveness and selection of viral sequences for inclusion in updated vaccines.

Cell culture
All cell lines were cultured in Dulbecco's Modified Eagle Medium (DMEM) (ThermoFisher Scientific, 21013024) supplemented with 10% fetal calf serum and 10 µg/ml gentamicin at 37°C and 5% CO 2 . Cells were periodically checked for mycoplasma and retrovirus contamination by DAPI staining and reverse transcriptase assays, respectively. The cell lines used have been previously described 24 . Derivatives of 293 T and HT1080 cells expressing ACE2 were generated by transducing 293 T cells with CSIB(ACE2). Single-cell clones were derived by limiting dilution from the bulk populations and clones designated 293 T/ACE2cl.22, and HT1080/ACE2cl.14 were used in this study.

SARS-CoV-2 variant spike expression plasmids
Generation of an omicron BA.1 spike expression plasmid lacking the C-terminal 19 amino acids and encoding the R683G furin cleavage site substitution (pCR3.1 BA.1Δ19) has been previously described 3 . The Omicron spike coding sequence was derived from sequence ID EPI_ISL_6640919. It was codon-optimized and synthesized as a C-terminally truncated Δ19 form in nine fragments (IDT). We also introduced a furin cleavage site mutation (R683G) that does not change the neutralization properties of the SASR-CoV-2 spike protein but enables higher titer pseudotyped viral stocks to be generated from transfected cells5. These synthetic DNA fragments, ranging in size from 444-599 bp and a NheI/XbaI-linearized pCR3.1 plasmid were Gibson assembled via 40 bps overlapping sequences. Individual plasmid clones were completely sequenced (Illumina MiSeq) and a single correct clone was used in these studies.

Monoclonal antibodies and plasma samples
Monoclonal antibodies were cloned from the memory B cells of individuals with varying exposure to SARS-CoV-2 antigen as a result of infection and/or vaccination.
Plasma samples were from five groups (i) Individuals who were infected with Wuhan-Hu-1-like viruses~1 y prior to blood collection, but not vaccinated 17 . (ii) Individuals who received two doses of mRNA vaccine (2nd dosẽ 1 m prior to plasma collection) 34 . (iii) Individuals who received three doses of mRNA vaccine (3rd dosẽ 1 m prior to plasma collection) 16 . (iv) Individuals from the same Wuhan-Hu-1-like infection cohort who were infected~1 y previously and also received two mRNA vaccine doses prior to blood collection 17 . (v) Individuals who received three mRNA vaccine doses and were then infected by omicron BA.1~1 m prior to blood collection 21,22 .
The study visits and blood draws were obtained with informed consent from all participants under a protocol that was reviewed and approved by the Institutional Review Board of the Rockefeller University (IRB no. DRO-1006, 'Peripheral Blood of Coronavirus Survivors to Identify Virus-Neutralizing Antibodies').

rVSV/SARS-CoV-2/GFP construction and rescue
The generation of rVSV/SARS-CoV-2/GFP chimeric viruses encoding SARS-2 spike proteins has been previously described 24 . A plaquepurified variant designated rVSV/SARS-CoV-2/GFP 2E1 encoding D215G/ R683G substitutions was used in this study. To introduce the spike proteins of VOCs BA.1.1 and BA.2, the respective spike sequences were amplified from the pCR3.1 constructs using PCR and primers specific for the VSV backbone. The rVSV backbone was digested with MluI/XhoI and the SARS-CoV-2 spike insert was introduced by Gibson assembly.
To rescue the recombinant viruses, HEK-293T/ACE2cl.22 cells were seeded 1 × 10 6 /well in poly-D-lysine-coated six-well plates 1 day prior to transfection. The following day, cells were rinsed with serum-free DMEM and infected with a recombinant vaccinia virus expressing T7 polymerase at an MOI of 5 for 45 min at 37°C, gently rocking every 10-15 min. The inoculum was removed and 1.5 ml DMEM supplemented with 10% FCS was added. Next, cells were transfected with a plasmid mixture containing the rVSV/SARS-CoV-2 genome (500 ng) and the helper plasmids pBS-N (300 ng) (Kerafast, EH1013), pBS-P (500 ng) (Kerafast, EH1014), pBS-L (100 ng) (Kerafast, EH1014) and pBS-G (800 ng) (Kerafast, EH1015) using the lipofectamine LTX (9 μl) (Ther-moFisher, 15338100) and PLUS (5.5 μl) (ThermoFisher, 15338100) transfection reagents. The DNA and PLUS reagent were mixed in 100 μl OptiMEM (ThermoFisher, 31985070) to which LTX was added in 105 μl OptiMEM. The reaction was incubated for 20 min at RT and then added to the cells. The transfected cells were monitored by microscopy and the supernatant of cells expressing GFP was collected~48 h post transfection, filtered using a 0.1 µm filter to remove residual vaccinia virus, and transferred to HEK-293T/ACE2cl.22 cells in six-well plates to amplify the rescued virus. After~48 h the supernatant of the infected cells was collected and filtered using a 0.22 µm filter and transferred onto HEK-293T/ACE2cl.22 cells were seeded in T175 flasks. The viruses were passaged several times in T175 flasks to create diversity in the viral population and generate virus stocks of >10 8 IU/ml. The titer of the viral stocks were assessed by serial dilution, infection of HEK-293T/ACE2cl.22 cells, and subsequent flow cytometry (see Supplementary Fig. 6) Antibody sequencing, cloning, and expression Antibodies were identified and sequenced 23 . In brief, RNA from single cells was reverse-transcribed (SuperScript III Reverse Transcriptase, Invitrogen, 18080-044) and the cDNA was stored at −20°C or used for subsequent amplification of the variable IGH, IGL, and IGK genes by nested PCR and Sanger sequencing. Sequence analysis was performed using MacVector. Amplicons from the first PCR reaction were used as templates for sequence-and ligation-independent cloning into antibody expression vectors. Recombinant monoclonal antibodies were subsequently produced and purified 23 . Paired heavy and light chain expression constructs were transfected into 293-6E cells (NRC) using branched polyethylenimine (PEI) 25 kDA (Sigma). After 7 days of culture, cells were centrifuged at 4200 x g for 40 mins at 4°C, and supernatants were filtered through 0.22 μM aPES (Thermo Nalgene Rapid-Flow). Antibodies were then purified from filtered supernatants using Protein G Sepharose 4 Fast Flow (GE Healthcare) according to standard protocols. Antibodies were buffer exchanged and concentrated into PBS using Amicon Ultra centrifugal filter (Millipore) with either a 30 or 50 kDA molecular weight cutoff.

Biolayer interferometry
Biolayer interferometry assays were performed as previously described 23 . In brief, we used the Octet Red instrument (ForteBio) at 30°C with shaking at 1000 r.p.m. Epitope binding assays were performed with protein A biosensor (ForteBio 18-5010), following the manufacturer's protocol "classical sandwich assay" as follows:

Selection of antibody-resistant rVSV/SARS-CoV-2 variants
To identify antibody escape mutations in rVSV/SARS-CoV-2/GFP 2E1 , rVSV/SARS-CoV-2/GFP BA1.1 , and rVSV/SARS-CoV-2/GFP BA2 , viral populations containing 1 × 10 6 infectious units were incubated with monoclonal antibodies (mAbs) at a concentration of 1 µg/mL for 1 h at 37°C. The virus/antibody mixture was then used to inoculate HEK-293T/ ACE2cl.22 cells in 6-well plates. The next day, the medium was replaced with fresh medium containing the monoclonal antibody at 1 µg/ml. After a further 24 h, the virus-containing supernatant was harvested and passed through a 0.22 μm 96-well filter plate. The filtered supernatant (100 μL) was added to medium containing mAbs to achieve a final mAb concentration of 1 µg/ml in a total volume of 1 ml. The virus:antibody mixture was incubated for 1 h at 37°C and used to inoculate 2 × 10 5 HEK-293T/ACE2cl.22 cells for a second passage (p2) in the presence of antibodies. Medium was again replaced with fresh antibody-containing medium after 24 h, and the putatively selected p2 virus population was harvested after 48 h. RNA was then extracted from 100 μL of filtered p2 supernatant and reverse-transcribed using the SuperScript VILO cDNA Synthesis Kit (Thermo Fisher Scientific). Sequences encoding the extracellular domain of the spike protein were initially amplified using KOD Xtreme Hot Start Polymerase (Sigma-Aldrich, 719753). In a subsequent PCR reaction, RBD-specific primers were used to amplify the receptor-binding domain (RBD), and the resulting products were purified and sequenced.

Receptor-binding domain sequencing
PCR products encoding the RBD were subjected to tagmentation reactions, that included 0.25 μl of Nextera TDE1 Tagment enzyme (Illumina, 20060059), 1.25 μl of TD Tagment buffer (Illumina, 20060059), and 10 ng of RBD PCR product. The tagmentation reaction was performed at 55°C for 5 min. Next, 3.5 μl of KAPA HiFi HotStart ReadyMix (Roche, 07958935001) and 1.25 μl of i5/i7 barcoded primers were added to the tagmented DNA to incorporate a unique indexing sequence. Following purification with AmpureXP beads (Beckman Coulter, A63882), the barcoded DNA from several reactions was pooled together and denatured with 0.2 N NaOH. The barcoded library was then diluted to 12 pM and sequenced using Illumina MiSeq Nano 300 V2 cycle kits (Illumina, MS-103-1001). Illumina sequencing reads were trimmed using the BBDuk function of Geneious Prime to remove adapter sequences and low-quality reads. The forward and reverse sequencing reads of each reaction were then aligned to the corresponding reference sequence and annotated for the presence of mutations. To detect RBD escape mutations, the minimum variant frequency (i.e.-the minimum fraction of reads that contained a substitution for it to be denoted as such) was set at 10%.
The generation of (HIV/NanoLuc) SARS-CoV-2 pseudotyped particles was described previously 24 . In brief, HEK-293T seeded in 10 cm dishes were transfected with 2.5 µg of a pCR3.1-SARS-CoV-2 spike expression plasmid and 7.5 µg of the pHIV-1 NL4-3 ΔEnv-NanoLuc reporter virus plasmid using 44 µl PEI (1 mg/ml) in 500 ml serum-free DMEM medium. The medium was changed 24 h post transfection. The virus-containing supernatant was harvested at 48 h post transfection, passed through a 0.22 µm pore-size polyvinylidene fluoride syringe filter, aliquoted, and stored at −80°C. The titers were assessed on HT1080/ACE2cl.14 cells seeded 10 4 cells/well in 100 µl medium in black flat-bottom 96-well plates. A fivefold serial dilution was performed and 100 µl was transferred to the target plate. The plates were incubated for 48 h and NanoLuc luciferase activity was determined as described below.

Neutralization assays
To assess the neutralization of pseudotyped HIV-1 viruses encoding the RBD point mutations that we had identified during our selection experiments by mAbs, neutralization assays were performed. First, HT1080 cells were seeded at 10 4 cells/well in 100 µl of medium. The next day, 60 µl of diluted virus (adjusted to give~10 7 relative light units in the infection assay) were incubated with 60 µl of mAb solution for 1 h at 37°C in a 96-well plate. Each mutant virus was initially screened for neutralization by the entire panel of monoclonal antibodies at a single antibody concentration of 1 µg/ml infection relative to no-antibody controls determined. For viruses deemed to be resistant to antibody neutralization, defined as having a relative infection value >0.4 relatives to the uninhibited virus control, IC 50 values were determined for parental virus and escape mutant. In this case, fivefold dilutions of the monoclonal antibody, starting at an initial concentration of 1 µg/ml, were incubated with the pseudotyped viruses. Thereafter, 100 µl of each antibody/virus mixture was added to the HT1080 cells and incubated at 37°C for 48 h. NanoLuc luciferase activity was then determined.

Reporter gene assays, and curve fitting
For NanoLuc luciferase assays, cells were washed twice with PBS and lysed with 30 µl of Cell Culture Lysis reagent (Promega, E1531). Nano-Luc luciferase activity was determined using the NanoGlo Luciferase Assay System (Promega, N1150). The substrate in NanoGlo buffer (30 µl) was added to the cell lysate. Luciferase activity was measured using a CLARIOstar Plus plate reader using 0.1 s integration time. The RLU readings were expressed as a decimal fraction of those derived from cells infected with pseudovirus in absence of antibodies.
To determine titers of rVSV/SARS-CoV-2/GFP, infected cells were trypsinized, fixed with 2% paraformaldehyde, resuspended in PBS/F-68, and the percentage of GFP-positive cells determined using an Attune NxT flow cytometer equipped with a 96-well CytKick Max autosampler.
Data were analyzed using Microsoft Excel for Mac (Version 16.65) and Prism V9.4.1 (GraphPad). NT 50 and IC 50 values were determined using a four-parameter nonlinear regression curve fit infectivity data measured as RLUs. The bottom values were set to zero, and the top values to one.

Statistics and reproducibility
Statistical significance between groups was assessed by a two-tailed t test using Welch's correction to account for unequal variances. All experiments were repeated successfully at least twice. No statistical method was used to predetermine the sample size. No data were excluded from the analyses; The experiments were not randomized. The Investigators were not blinded to allocation during experiments and outcome assessment, but the complexity of the experiments rendered the identification of single samples during experiments unlikely.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability
All source data for antibody escape mutation frequencies and for antibody neutralization information shown in the heatmaps are included in the supplementary data files 3 and 4.