SUMMARY
Vaccination elicits immune responses capable of potently neutralizing SARS-CoV-2. However, ongoing surveillance has revealed the emergence of variants harboring mutations in spike, the main target of neutralizing antibodies. To understand the impact of globally circulating variants, we evaluated the neutralization potency of 48 sera from BNT162b2 and mRNA-1273 vaccine recipients against pseudoviruses bearing spike proteins derived from 10 strains of SARS-CoV-2. While multiple strains exhibited vaccine-induced cross-neutralization comparable to wild-type pseudovirus, 5 strains harboring receptor-binding domain mutations, including K417N/T, E484K, and N501Y, were highly resistant to neutralization. Cross-neutralization of B.1.351 variants was weak and comparable to SARS-CoV and bat-derived WIV1-CoV, suggesting that a relatively small number of mutations can mediate potent escape from vaccine responses. While the clinical impact of neutralization resistance remains uncertain, these results highlight the potential for variants to escape from neutralizing humoral immunity and emphasize the need to develop broadly protective interventions against the evolving pandemic.
INTRODUCTION
Since the first described human infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in December of 2019, nine vaccines have been approved for use in humans (“COVID-19 Vaccine Tracker” n.d.). Two of the vaccines currently in use worldwide, BNT162b2 (manufactured by Pfizer) and mRNA-1273 (manufactured by Moderna), are based on lipid nanoparticle delivery of mRNA encoding a prefusion stabilized form of spike protein derived from SARS-CoV-2 isolated early in the epidemic from Wuhan, China. Both of these vaccines demonstrated >94% efficacy at preventing coronavirus disease 2019 (COVID-19) in phase III clinical studies performed in late 2020 in multiple countries (Polack et al. 2020; Baden et al. 2021). However, the recent emergence of novel circulating variants has raised significant concerns about geographic and temporal efficacy of these interventions. Indeed, more recently completed trials of two adenovirus-based vaccines (AZD1222 from Astrazeneca and JNJ-78436735 from Johnson & Johnson), a nanoparticle-based vaccine (NVX-CoV2373 from Novavax), and an inactivated protein vaccine (Coronavac) have demonstrated reduced overall efficacy (“Novavax COVID-19 Vaccine Demonstrates 89.3% Efficacy in UK Phase 3 Trial” n.d., “AZD1222 Vaccine Met Primary Efficacy Endpoint in Preventing COVID-19” 2020), and subset analyses suggest marked geographic variation with lower efficacy against mild-to-moderate disease in countries such as South Africa and Brazil where the epidemic is dominated by variant strains. Taken together, these data suggest that neutralization-resistant variants may have contributed to these outcomes (“South Africa Suspends Use of AstraZeneca’s COVID-19 Vaccine after It Fails to Clearly Stop Virus Variant” 2021).
One of the earliest variants that emerged and rapidly became globally dominant was D614G. While several studies demonstrated that this strain is more infectious (Korber et al. 2020; Yurkovetskiy et al. 2020a, [b] 2020; Plante et al. 2020; Zhou et al. 2020; Hou et al. 2020), we and others found that sera from convalescent individuals showed effective cross-neutralization of both wildtype and D614G variants (Garcia-Beltran et al. 2021; Legros et al. 2021; Hou et al. 2020). However, recent genomic surveillance in the United Kingdom has revealed rapid expansion of a novel lineage termed B.1.1.7 (also known as VOC-202012/01 or 501Y.V1). B.1.1.7 harbors 3 amino acid deletions and 7 missense mutations in spike, including D614G as well as N501Y in the ACE2 receptor-binding domain (RBD), and has been reported to be more infectious than D614G (Santos and Passos 2021; Galloway et al. 2021; Liu et al. 2021). Several studies have demonstrated that convalescent and vaccinee sera cross-neutralize B.1.1.7 variants with only slightly decreased potency, suggesting that prior infection or vaccination with wild-type SARS-CoV-2 may still provide protection against B.1.1.7 variants (Wu et al. 2021; Muik et al. 2021; Shen et al. 2021; Rees-Spear et al. 2021; P. Wang et al. 2021). There have also been reports of SARS-CoV-2 transmission between humans and minks in Denmark with a variant called mink cluster 5 or B.1.1.298, which harbors a 2-amino acid deletion and 4 missense mutations including Y453F in RBD. Concerns relating to ongoing interspecies transmission resulted in the culling of over 17 million Danish minks to prevent further viral spread and evolution (Oude Munnink et al. 2021; Oxner 2020). Another variant that recently emerged in California, United States, designated as B.1.429, contains 4 missense mutations in spike, one of which is a single L452R RBD mutation. The ability of B.1.1.298 and B.1.429 variants to evade neutralizing humoral immunity from prior infection or vaccination has yet to be determined.
Of particular concern is an E484K mutation in RBD, which was recently identified through deep mutational scanning as a variant with the potential to evade monoclonal and serum antibody responses (Greaney et al. 2020, 2021). Novel variants arising from the B.1.1.28 lineage first described in Brazil and Japan, termed P.2 (with 3 spike missense mutations) and P.1 (with 12 spike missense mutations), contain this E484K mutation, and P.1 in particular also contains K417T and N501Y mutations in RBD. These strains have been spreading rapidly, and both P.2 and P.1 were recently found in documented cases of SARS-CoV-2 re-infection (Paiva et al. 2020; Faria et al. 2021; Resende et al. 2021; Naveca et al. 2021; Nonaka et al. 2021).
Of greatest concern has been the emergence of multiple strains of the B.1.351 lineage (also known as 501Y.V2), which were first reported in South Africa and have since spread globally (Tegally et al. 2021). This lineage bears three RBD mutations, K417N, E484K, and N501Y, in addition to several mutations outside of RBD, and several reports have suggested that convalescent and vaccinee sera have decreased cross-neutralization of B.1.351 lineage variants (P. Wang et al. 2021; Wibmer et al. 2021; Wu et al. 2021; Hu et al. 2021). A key limitation of several of these reports has been the use of single mutations or combinations of mutations that do not naturally occur. Regardless, the emergence of novel variants that appear to escape immune responses has spurred vaccine manufacturers to develop boosters for these spike variants (“Moderna Developing Booster Shot for New Virus Variant B.1.351” n.d.).
Here, we systematically assessed the neutralization potential of sera from a primarily young cohort of individuals who received one or two doses of the BNT162b2 (Pfizer) or mRNA-1273 (Moderna) vaccine against SARS-CoV-2 pseudoviruses assembled to recapitulate circulating strains. We used our previously described high-throughput pseudovirus neutralization assay (Garcia-Beltran et al. 2021) to quantify neutralization against variants first arising in the United Kingdom (B.1.1.7), Denmark (B.1.1.298), United States (B.1.429), Brazil and Japan (P.2 and P.1), and South Africa (three variants of the B.1.351 lineage), as well as SARS-CoV from the 2002 Hong Kong outbreak and the pre-emergent bat coronavirus WIV1-CoV. We find that although neutralization is largely preserved against many variants, those containing the K417N/T, E484K, and N501Y RBD mutations, namely, P.1 and B.1.351 variants, have significantly decreased neutralization even in fully vaccinated individuals. Individuals that received only a single recent dose of vaccine had weaker neutralization titers overall and did not exhibit detectable neutralization of B.1.351 variants in our assays.
Taken together, our results highlight that BNT162b2 and mRNA-1273 vaccines achieve only partial cross-neutralization of novel variants and support the reformulation of existing vaccines to include diverse spike sequences. Ultimately, development of new vaccines capable of eliciting broadly neutralizing antibodies may be necessary to resolve the ongoing pandemic.
RESULTS
Emergence and global spread of SARS-CoV-2 variants of concern
Since the first described case of SARS-CoV-2 infection in Wuhan, China in December of 2019, over 107 million confirmed infections have been documented (Coronavirus.jhu.edu), enabling viral diversification and the emergence of six distinct major lineages with numerous variants. A subset of these variants have been denoted as variants of concern by the World Health Organization given the presence of mutations with potential to increase transmissibility, virulence, or evade immune responses. We focused on variants of concern first described in the United Kingdom (B.1.1.7), Denmark (B.1.1.298), United States (B.1.429), Brazil and Japan (P.2 and P.1), and South Africa (B.1.351), most of which arose in late 2020 (Figure 1A-B). Although the naturally arising mutations in these variants span the entire spike protein, they mainly occur in S1 and RBD, the main target of neutralizing antibodies (Figure 1C and S1). After analyzing SARS-CoV-2 spike sequences deposited for each lineage in GISAID, we identified consensus sequence mutations that represent the dominant circulating strain for each lineage (Figure 1D and S1). In the case of the B.1.351 lineage, we studied three of the dominant variants comprising the majority of deposited sequences, which we designated v1, v2, and v3. The three main RBD mutations of concern are: (i) N501Y, present in B.1.1.7, P.1, and B.1.351 variants; (ii) E484K, present in the P.2, P.1, and B.1.351 variants; and (iii) K417T for the P.1 variant and K417N for the B.1.351 variants. Separately, the B.1.1.298 variant found in Danish minks contained a Y453F mutation in RBD, and the California variant B.1.429 contained an L452R. Although distantly related, SARS-CoV from the 2002 Hong Kong coronavirus outbreak and pre-emergent bat-derived WIV1-CoV both share ∼76% spike homology to SARS-CoV-2 and were assessed to provide comparison to serologically distinct coronaviruses (Figure 1A and 1D).
Vaccine-elicited neutralizing responses to circulating SARS-CoV-2 strains
We accrued a cohort of individuals who received either one or two doses of the BNT162b2 or mRNA-1273 vaccines, and used a luminescence-based lentiviral pseudovirus neutralization assay that we and others have previously validated (C. Wang et al. 2020; Ju et al. 2020; Pinto et al. 2020; Yang et al. 2020; Moore et al. 2004; Crawford et al. 2020; Garcia-Beltran et al. 2021) to assess neutralization of SARS-CoV-2 variants in a high-throughput system (Figure 2A-B). This cohort was relatively young (median age: 33 years, range: 23 - 67 years) and was 56% (27/48) female and 44% (21/48) male. We found that individuals who received 2 full doses of BNT162b2 or mRNA-1273 (>7 days out from second dose) had robust neutralization of wild-type SARS-CoV-2, while those who received only one dose had lower but detectable neutralization (Figure 2C). The magnitude of neutralization correlated with total anti-RBD antibody levels (Figure 2D). Although sex did not appear to influence wild-type neutralization titers in our cohort, there appeared to be an age dependency, with younger individuals in their 20s and 30s having comparatively higher responses than older individuals (Figure 2E-F).
Interestingly, neutralization of D614G was slightly decreased in individuals who received two doses of vaccine (∼2-fold decrease for 2-dose BNT162b2) (Figure 3A-C and S2), which was in contrast to previous studies in convalescent sera that we and others conducted demonstrating slightly increased neutralization of D614G variant versus wild type (Garcia-Beltran et al. 2021) following natural infection. The effect was even more pronounced in individuals who received one vaccine (p < 0.05 for 1-dose mRNA-1273), some of which had undetectable neutralization of D614G despite detectable neutralization of wild-type SARS-CoV-2 (Figure 3B). This difference may be a consequence of the vaccine encoding the wild-type spike sequence, while many convalescent individuals in previous studies were likely infected with D614G variant SARS-CoV-2, given it had already become the globally dominant strain by the summer of 2020 (Lemieux et al. 2020; Korber et al. 2020).
When assessing variants containing one RBD mutation as part of their mutational landscape, the UK variant B.1.1.7 (N501Y), Danish mink variant B.1.1.298 (Y453F), and California variant B.1.429 (L452R) exhibited neutralization that was similar to that of wild-type and the parental D614G variant. Among people who received 2 doses of BNT162b2, the mean fold decrease in neutralization relative to wild type was 2.3-fold for B.1.1.7, 3.4-fold for B.1.1.298, and 2.9-fold for B.1.429 (Figure 3C and S2). However, neutralization of the Brazilian/Japanese P.2 variant, whose RBD contains a E484K mutation, was significantly decreased (13.4-fold, p < 0.001) (Figure 3C and S2). This is in line with previous studies suggesting that the E484K mutation can evade polyclonal antibody responses (Greaney et al. 2020; Jangra et al. 2021) and has been found in cases of SARS-CoV-2 re-infection (Paiva et al. 2020; Faria et al. 2021; Naveca et al. 2021; Resende et al. 2021; Nonaka et al. 2021). Similarly, neutralizing antibody responses were also significantly decreased for the Brazilian/Japanese P.1 strain (15.1-fold, p < 0.0001), which harbors three mutations in RBD (K417T, E484K, and N501Y) and has also been found in cases of re-infection.
Strikingly, neutralization of all three South African B.1.351 strains was substantially decreased (75.4-fold for v1, 97.2-fold for v2, and 93.6-fold for v3, p < 0.0001 for all three). These strains contain the same three RBD mutations as P.1 except for an asparagine versus threonine substitution at K417 (K417N) (Figure 3C and S2) and several additional mutations in non-RBD regions. Surprisingly, we found that the neutralization of B.1.351 v2 and v3 strains was worse than that of more distantly related coronaviruses SARS-CoV (84.1-fold decrease) and WIV1-CoV (72.3-fold decrease), suggesting that a relatively small number of mutations can medicate potent escape from vaccine-induced neutralizing antibody responses. These substantial differences in neutralization could not be explained by differences in spike protein expression (Figure 3D) nor in the amount of pseudovirus used in the assay (Figure S3). Notably, 36% (8/22) recipients of 2-dose BNT162b2 and 50% (2/4) recipients of 2-dose mRNA-1273 vaccines did not have detectable neutralization of B.1.351 v2 (Figure 3A-B). Of the individuals who received only 1 dose of BNT162b2 or mRNA-1273, 0% (none) had detectable neutralization of B.1.351 v2 and v3, except for two individuals who reported having prior COVID-19 infection or a significant exposure. This suggests that a single dose of existing mRNA vaccines may be insufficient to induce cross-neutralizing antibody responses in previously uninfected individuals. Interestingly, the three individuals who reported having prior COVID-19 infection or significant exposure had the highest neutralization titers for most variants and exhibited substantial cross-neutralization of the B.1.351, SARS-CoV, and WIV1-CoV (Figure 3A-B). This suggests that prior infection combined with vaccination may result in the greatest breadth of cross-reactive neutralizing antibody responses, even against distantly related coronaviruses.
Neutralization resistance of B.1.351 variants is primarily due to RBD mutations
To better characterize the mutational context responsible for B.1.351 neutralization resistance, we explored the contribution of mutations located both within and outside of the RBD region of spike. Of note, many of the substitutions and deletions found in the B.1.351 lineage are in close structural proximity in the S1 domain (Figure 4A). However, 86.2% of available GISAID sequences of this lineage harbor these mutations in three distinct patterns represented by B.1.351 v1 (46.6%), v2 (31.6%) and v3 (8.0%). We created pseudoviruses in which the RBD regions of B.1.351 v1 and v2 were reverted to their original wild-type Wuhan sequence (v1/wtRBD and v2/wtRBD) while retaining all other B.1.351 v1 or v2 mutations. Remarkably, neutralization assays conducted with sera from 22 individuals that received the 2-dose BNT162b2 vaccine revealed that neutralization of B.1.351 v1 and v2 in the absence of RBD mutations was comparable to that of D614G (Figure 4B-C). A pseudovirus bearing only the three RBD mutations (K417N, E484K, and N501Y), largely, but not entirely, recapitulated the escape phenotype (Figure 4C). In regression analyses, neutralization of the triple RBD mutant explained the majority of variance in neutralization of B.1.351 v1 (90.5%) and v2 (75.3%), while a minor, but significant, fraction was explained by mutations outside of RBD: v1/wtRBD (3.1%) and v2/wtRBD (17.6%). Moreover, there was significant statistical interaction between mutations within and outside of the RBD in mediating escape for both variants, consistent with synergistic (ie. non-additive) effects (p < 0.05). These results suggest that while RBD mutations contribute the majority of the observed escape from vaccine-induced neutralization, they are more effective when in the context of additional mutations, particularly those found in B.1.351 v2 variants.
DISCUSSION
Traditionally, polyclonal immune responses that arise in the context of infection and vaccination are thought to target multiple antigenic epitopes. Given this assumption, the expectation would be for small numbers of variations in antigen sequence to have only modest effects on recognition by the immune system. Here, we find that while many strains, such as B.1.1.7, B.1.1.298, or B.1.429, continue to be potently neutralized despite the presence of individual RBD mutations, other circulating SARS-CoV-2 variants escape vaccine-induced humoral immunity. The P.2. variant, which contains an E484K mutation within the RBD region, was capable of significantly reducing neutralization potency of fully vaccinated individuals, in line with what has been suggested by deep mutational scanning (Greaney et al. 2020; Jangra et al. 2021). Similarly, the P.1 strain, which has three RBD mutations, more effectively escaped neutralization, possibly explaining recently reported cases of re-infection with this variant (Paiva et al. 2020; Faria et al. 2021; Resende et al. 2021; Naveca et al. 2021; Nonaka et al. 2021). Finally, we found that B.1.351 variants exhibited remarkable resistance to neutralization, largely due to three mutations in RBD but with measurable contribution from non-RBD mutations. The magnitude of the effect is such that B.1.351 strains escape neutralizing vaccine responses as effectively as distantly related coronaviruses.
Given the loss of vaccine potency against a number of circulating variants, individuals receiving a single dose of vaccine did not raise sufficient antibody titers to provide any detectable cross neutralization against B.1.351 v2 or v3. While our studies are limited by the relatively short follow-up time after vaccination, our findings support the importance of 2-dose regimens to achieve titers, and perhaps breadth, to enhance protection against novel variants. These findings are important to consider in the context of proposals to administer a single dose of vaccine across a larger number of individuals instead of using doses to boost prior recipients
Importantly, our studies rely on pseudoviruses which are only capable of modeling the ACE2-dependent entry step of the SARS-CoV-2 lifecycle. While numerous studies have now demonstrated a close correlation between neutralization titers measured against pseudovirus and live SARS-CoV-2 cultures (Wang et al. 2020; Ju et al. 2020; Pinto et al. 2020; Yang et al. 2020; Moore et al. 2004; Crawford et al. 2020), it is unclear what impact additional mutations located outside of the spike may have on immunological escape, virulence, infectivity, or pathogenesis. Several recent reports and pre-prints, including studies conducted by Pfizer as well as Moderna, have produced similar findings in terms of vaccine potency against B.1.1.7 and B.1.1.298 variants, but substantially less neutralization resistance by B.1.351 than we measured. However, we would caution that each study was done with distinct serum samples measured using different conditions (e.g. replication-competent chimeric vesicular stomatitis virus versus single-entry lentiviral pseudovirus) and with different spike expression plasmids and combinations of mutations (e.g. triple RBD mutations, lineage-defining mutations, full complement of mutations not yet observed together in nature, etc.). In the present study, we mimicked the natural pattern of mutations found in circulating strains, in phase. Importantly, these include the spike variants found in two replication-competent live-virus strains of the B.1.351 lineage, which were recently reported to exhibit nearly complete neutralization resistance in response to convalescent plasma (Cele et al. 2021). The strains used in the aforementioned study were designated as 501Y.V2.HVdF002 and 501Y.V2.HV001, which have identical spike protein sequences as the B.1.351 v1 and B.1.351 v2 pseudoviruses reported in this study.
One important aspect of immunity not addressed by our work is cellular immunity contributed by cytotoxic lymphocytes, including T and NK cells. Even in the absence of neutralizing humoral immunity, previous studies have suggested that cellular immunity can mitigate severe or prolonged infection (Le Bert et al. 2020). In convalescent individuals, T-cell immunity would not be restricted to spike-derived epitopes, but also from other more abundant proteins such as nucleocapsid. As such, it would be reasonable to assume that T-cell-mediated immunity elicited by infection would remain largely intact for circulating variants including B.1.351. However, with the exception of killed whole virus vaccines, all currently available vaccine designs only provide spike protein as the target immunogen, thus limiting T-cell immunity to spike epitopes. It has yet to be determined if the circulating mutations in these epitopes will similarly impair T-cell responses in addition to escaping neutralizing antibodies.
In summary, our data highlights the challenges facing all vaccines whose designs were finalized early in the pandemic and based on the sequence of the first-reported virus from Wuhan, China. Given the global scale and magnitude of the ongoing pandemic, including case reports of re-infection, it is clear that viral evolution will continue. It is possible that currently available vaccines will still provide clinical benefit against variants that exhibit poor cross-neutralization, such as P.1 and B.1.351, by reducing COVID-19 disease severity, but this has yet to be determined. Ultimately, it will be important to develop interventions capable of preventing transmission of diverse SARS-CoV-2 variants, including vaccine boosters that target these variants or technologies capable of eliciting or delivering broadly neutralizing antibodies.
Data Availability
This study assessed the neutralization potency of sera obtained from vaccinated individuals and was not a clinical trial. De-identified raw data of the results of serological assays can be provided upon request.
AUTHOR CONTRIBUTIONS
W.F.G.B., E.C.L., K.S.D., and A.B.B. designed the experiments. W.F.G.B., E.C.L., K.S.D., Z.G., A.D.N., V.N., and A.B.B. carried out experiments and analyzed data. J.F., B.M.H., and A.G.S. provided key reagents and useful discussions and insights. D.J.G., M.C.P., and M.N.P. provided human samples for the study. A.D.N. and V.N. contributed to statistical and sequence analyses. A.J.I. and V.N. provided key discussions and input into experimental design. W.F.G.B., V.N., and A.B.B. wrote the paper with contributions from all authors.
DECLARATIONS OF INTEREST
There are no conflicts of interest for the authors to declare.
SUPPLEMENTAL FIGURE TITLES AND LEGENDS
STAR ★ METHODS
RESOURCE AVAILABILITY
Lead Contact
Further information and requests for resources and reagents should be directed to and will be fulfilled by Alejandro Balazs (abalazs{at}mgh.harvard.edu).
Materials Availability
Plasmids generated in this study will be available through Addgene. Recombinant proteins and antibodies are available from their respective sources.
Data and Code Availability
This study did not generate sequence data or code. Data generated in the current study (including ELISA, neutralization, and cytokine measurements) have not been deposited in a public repository but are available from the corresponding author upon request.
EXPERIMENTAL MODEL AND SUBJECT DETAILS
Human subjects
Use of human samples was approved by Partners Institutional Review Board (protocol 2020P002274). Serum samples from 48 vaccine recipients that received one or two doses of the BNT162b2 or mRNA-1273 vaccine were collected. For each individual, basic demographic information including age and sex as well as any relevant COVID-19 history was obtained.
Cell lines
HEK 293T cells (ATCC) were cultured in DMEM (Corning) containing 10% fetal bovine serum (VWR), and penicillin/streptomycin (Corning) at 37°C/5% CO2. 293T-ACE2 cells were a gift from Michael Farzan (Scripps Florida) and Nir Hacohen (Broad Institute) and were cultured under the same conditions. Confirmation of ACE2 expression in 293T-ACE2 cells was done via flow cytometry.
METHOD DETAILS
Construction of variant spike expression plasmids
To create variant spike expression plasmids, we performed multiple PCR fragment amplifications utilizing oligonucleotides containing each desired mutation (Integrated DNA Technology) and utilized overlapping fragment assembly to generate the full complement of mutations for each strain. Importantly we generate these mutations in the context of our previously described codon-optimized SARS-CoV-2 spike expression plasmid harboring a deletion of the C-terminal 18 amino acids that we previously demonstrated to result in higher pseudovirus titers. Assembled fragments were inserted into NotI/XbaI digested pTwist-CMV-BetaGlobin-WPRE-Neo vector utilizing the In-Fusion HD Cloning Kit (Takara). All resulting plasmid DNA utilized in the study was verified by whole-plasmid deep sequencing (Illumina) to confirm the presence of only the intended mutations.
SARS-CoV-2 pseudovirus neutralization assay
To compare the neutralizing activity of vaccinee sera against coronaviruses, we produced lentiviral particles pseudotyped with different spike proteins as previously described (Garcia- Beltran et al. 2021). Briefly, pseudoviruses were produced in 293T cells by PEI transfection of a lentiviral backbone encoding CMV-Luciferase-IRES-ZsGreen as well as lentiviral helper plasmids and each spike variant expression plasmid. Following collection and filtering, production was quantified by titering via flow cytometry on 293T-ACE2 cells. Neutralization assays and readout were performed on a Fluent Automated Workstation (Tecan) liquid handler using 384-well plates (Grenier). Three-fold serial dilutions ranging from 1:12 to 1:8,748 were performed for each serum sample before adding 50-250 infectious units of pseudovirus for 1 h. Subsequently, 293T-ACE2 cells containing polybrene were added to each well and incubated at 37°C/5% CO2 for 60-72 h. Following transduction, cells were lysed using a luciferin-containing buffer (Siebring-van Olst et al. 2013) and shaken for 5 min prior to quantitation of luciferase expression within 1 h of buffer addition using a Spectramax L luminometer (Molecular Devices). Percent neutralization was determined by subtracting background luminescence measured in cell control wells (cells only) from sample wells and dividing by virus control wells (virus and cells only). Data was analyzed using Graphpad Prism and NT50 values were calculated by taking the inverse of the 50% inhibitory concentration value for all samples with a neutralization value of 80% or higher at the highest concentration of serum.
Titering
To determine the infectious units of pseudotyped lentiviral vectors, we plated 400,000 293T-ACE2 cells per well of a 12-well plate. 24 h later, three ten-fold serial dilutions of neat pseudovirus supernatant were made in 100 µL, which was then used to replace 100 µL of media on the plated cells. Cells were incubated for 48 h at 37°C/5% CO2 to allow for expression of ZsGreen reporter gene and harvested with Trypsin-EDTA (Corning). Cells were resuspended in PBS supplemented with 2% FBS (PBS+), and analyzed on a Stratedigm S1300Exi Flow Cytometer to determine the percentage of ZsGreen-expressing cells. Infectious units were calculated by determining the percentage of infected cells in wells exhibiting linear decreases in transduction and multiplying by the average number of cells per well determined at the initiation of the assay. At low MOI, each transduced ZsGreen cell was assumed to represent a single infectious unit.
Spike Expression
To compare the relative surface expression of pseudovirus spike variant proteins, we plated 400,000 293T cells per well of a 12-well plate. 24 h later, 1 µg of variant spike expression plasmid was mixed with 250ng of vector encoding GFP and transfected using PEI. Cells were incubated for 48 h at 37°C and harvested into PBS+ Cells transfected with each vector were divided and stained with 5 µg/mL of either, B38, ADI-55689, or ADI-56046 for 30 minutes at room temperature. Cells were then washed with 1 mL PBS+, spun at 900 x g, and stained with 2 µg/mL of anti-human IgG-AF647 polyclonal antibody (Invitrogen) for 30 minutes at room temperature. Cells were washed with 1 mL of PBS+, spun at 900 x g, resuspended in 100 µL of PBS+, fixed with 100ul 4% PFA and analyzed on a Stratedigm S1300Exi Flow Cytometer.
SARS-CoV-2 receptor binding domain and spike IgG, IgM, and IgA ELISA
Quantitative detection of total antibodies to SARS-CoV-2 receptor binding domain (RBD) was performed as previously described (Garcia-Beltran et al. 2021). Briefly, we used an indirect ELISA with a standard consisting of anti-SARS-CoV and -CoV-2 monoclonal antibody (CR3022) (IgG1 isotype). 96-well ELISA plates were coated with purified wild-type SARS-CoV-2 RBD. Plates were blocked with BSA and washed. A seven-point standard curve was created using CR3022-IgG1 starting at 2 µg/mL by performing 1:3 serial dilutions with dilution buffer, and serum samples were diluted 1:100 with dilution buffer. Diluted samples and standards were added to corresponding wells and incubated for 1 h at 37°C, followed by washing. Total antibodies were detected with anti-human IgG/A/M(H+L)-HRP (Bethyl) diluted 1:25,000 for a 30 min incubation at room temperature. After washing, TMB substrate (Inova) was added to each well and incubated for 5-15 min before stopping with 1 M H2SO4. Optical density (O.D.) was measured at 450 nm with subtraction of the O.D. at 570 nm as a reference wavelength on a SpectraMax ABS microplate reader. Anti-RBD antibody levels were calculated by interpolating onto the standard curve and correcting for sample dilution; one unit per mL (U/mL) was defined as the equivalent reactivity seen by 1 µg/mL of CR3022.
QUANTIFICATION AND STATISTICAL ANALYSIS
Data and statistical analyses were performed using GraphPad Prism 9.0.1, JMP Pro 15.0.0 (SAS Institute), and R v4.0.2. Flow cytometry data was analyzed using FlowJo 10.7.1. Non-parametric multivariate ANOVAs were performed on the indicated figures where several cohorts were present; all p values were adjusted for multiple comparisons. Statistical significance was defined as p < 0.05. Error bars throughout all figures represent one standard deviation.
ACKNOWLEDGEMENTS, FUNDING SUPPORT
We wish to thank Michael Farzan, PhD for providing ACE2-expressing 293T cells. B.M.H. is supported by award Number T32GM007753 from the National Institute of General Medical Sciences. J.F. is supported by T32AI007245. D.J.G., M.C.P., and M.N.P. were supported by the VIC Innovation fund. A.S was supported by the Bill and Melinda Gates Investment INV-018944 (AS) and by the South African Medical Research Council and the Department of Science and Innovation (TdO), A.G.S. was supported by NIH R01 AI146779 and a Massachusetts Consortium on Pathogenesis Readiness (MassCPR) grant. A.J.I. is supported by the Lambertus Family Foundation. A.B.B. was supported by the National Institutes for Drug Abuse (NIDA) Avenir New Innovator Award DP2DA040254, the MGH Transformative Scholars Program as well as funding from the Charles H. Hood Foundation. This independent research was supported by the Gilead Sciences Research Scholars Program in HIV.