Abstract
Pre-existing antibodies to endemic coronaviruses (CoV) that cross-react with SARS-CoV-2 have the potential to influence the antibody response to COVID-19 vaccination and infection for better or worse. In this observational study of mucosal and systemic humoral immunity in acutely infected, convalescent, and vaccinated subjects, we tested for cross reactivity against endemic CoV spike (S) protein at subdomain resolution. Elevated responses, particularly to the β-CoV OC43, were observed in all natural infection cohorts tested and were correlated with the response to SARS-CoV-2. The kinetics of this response and isotypes involved suggest that infection boosts preexisting antibody lineages raised against prior endemic CoV exposure that cross react. While further research is needed to discern whether this recalled response is desirable or detrimental, the boosted antibodies principally targeted the better conserved S2 subdomain of the viral spike and were not associated with neutralization activity. In contrast, vaccination with a stabilized spike mRNA vaccine did not robustly boost cross-reactive antibodies, suggesting differing antigenicity and immunogenicity. In sum, this study provides evidence that antibodies targeting endemic CoV are robustly boosted in response to SARS-CoV-2 infection but not to vaccination with stabilized S, and that depending on conformation or other factors, the S2 subdomain of the spike protein triggers a rapidly recalled, IgG-dominated response that lacks neutralization activity.
Introduction
The ongoing pandemic of SARS-CoV-2 represents the third time in just two decades that a novel coronavirus (CoV) with significant morbidity and mortality has begun to circulate among humans1. Given the alarming frequency of these occurrences, the COVID-19 pandemic serves as a call to action to safeguard against continued emergence of novel human CoV. Unlike in the earlier outbreaks of SARS-CoV and MERS-CoV, however, SARS-CoV-2 has proven to be highly transmissible, infecting a substantial portion of the global population – a conservative estimate of the prevalence of infection based on only confirmed cases corresponds to over 2%2, while estimates of the total figure in the United States exceed 16% of the population3. The extent to which the immune responses generated by these and other endemic CoV exposures might serve as an effective deterrent against the emergence of novel CoV strains remains an open question.
Insights into the effect of exposure to endemic CoVs in the years leading up to the COVID-19 pandemic could suggest whether preexisting antibodies are likely to contribute beneficially or detrimentally to outcomes of infection by a novel strain. Under the hypothesis of original antigenic sin, preexisting humoral memory imprinted from prior exposure to related antigens partially predestines the antibody repertoire to focus on the epitopes of a new threat that closely resemble those for which there is an existing solution. Because neutralizing epitopes are subject to additional selective pressure, this recall and re-diversification of existing antibody lineages may hinder the immune system’s ability to generate effective neutralizing antibodies4. This phenomenon has been most convincingly established in the context of influenza virus infection, with support from both well-controlled animal model experiments4,5 and observational studies of human natural infection histories6-10. Similar observations have been made in the context of diverse dengue virus serotypes11. These and other studies have shown despite being a key hallmark of effective long-term immune defense, anamnestic responses are not without potential downsides.
Antibody-dependent enhancement (ADE), a phenomenon by which cross-reactive antibodies induced by a prior exposure to a related pathogen promote infection of cell types bearing antibody Fc receptors and potentially elevate morbidity and mortality (reviewed in12,13), suggests further potential advantages of a humoral blank slate. Though observed in vitro14,15, there has not been abundant evidence for the biological relevance of classical ADE in the context of SARS-CoV-2 in vivo15,16. Nonetheless, the role of non-neutralizing antibodies resulting from cross-strain challenges in promoting ADE in cases of dengue fever17 has motivated concern about COVID-19 enhancement18-20.
The extent to which preexisting responses to prior endemic CoV exposures may influence responses to SARS-CoV-2 infection and vaccination is not yet clear, but has been suggested from studies relating recent endemic CoV infection with reduced severity of COVID-1921,22. Numerous studies have observed elevated responses to endemic CoV following SARS-CoV-2 infection23-27, and more recent work has shown that the magnitude of this “back-boosting” effect is inversely correlated with the induction of IgG and IgM against to SARS-CoV-2 spike (S) protein28. Given differential degrees of homology between the receptor binding domain (RBD) in S1 that is the target of the majority of neutralizing antibodies, and the better conserved S2 domain that may be the target of antibodies with diverse effector functions but which are rarely neutralizing, the original antigenic sin hypothesis suggests that lower titers of neutralizing antibodies against SARS-CoV-2 may result from boosting of pre-existing cross-reactive lineages. Here, this apparent boosting effect is evaluated in diverse cohorts with the goal of beginning to understand its implications.
Results
Structural Analysis of SARS-CoV-2
Sequence conservation across the human CoV spike protein is not uniformly distributed: the N-terminal domain (NTD) of SARS-CoV-2 that contains the RBD responsible for interacting with human angiotensin-converting enzyme 2 (ACE2) and initiating viral entry has lower homology to the corresponding subdomains of other CoVs than do the S2 subdomains where the fusion peptide, heptad repeat, and central helix required for fusion are located (Figure 1A). In contrast to the S2 domain, structural comparison of various subdomains of the SARS-CoV-2 S protein show that both the NTD and RBD of the spike protein are poorly conserved across CoV (Figure 1B). Superimposition of the SARS CoV-2 S1 and S2 domains with the most well conserved widely circulating endemic human CoV, OC43, show high structural conservation in S2 and the NTD, and a complete lack of homology in the RBD that is consistent with the differing entry receptors used by these β-CoV (Figure 1C, Supplemental Figure 1). Based on both structural and sequence homology, it stands to reason that preexisting antibodies raised against endemic human CoV are more likely to target the better conserved regions of SARS-CoV-2 such as S2, and less likely to recognize the RBD29,30.
Subject Cohorts
A diversity of cohorts were evaluated in this observational study (Table 1). Subjects naturally infected with SARS-CoV-2 comprised a small cohort of convalescent subjects for which serum and mucosal samples were available, a larger cohort of convalescent plasma donors, acutely infected individuals, pregnant women infected in their third trimester, and a small set of subjects for whom pre- and post-infection samples were available. Subjects vaccinated against SARS-CoV-2 with mRNA included cohorts of healthy adults, and pregnant women vaccinated in their third trimester. Samples were analyzed alongside samples from naïve (n = 15) and historical (n = 38) negative controls.
Elevated Responses to Endemic CoV in Serum, Nasal Wash, and Stool Among SARS-CoV-2 Infected Subjects
The magnitude and specificity of IgM, IgA, and IgG responses were determined across a panel of SARS-CoV-2 and endemic CoV antigens (Supplemental Table 2) in two previously described cohorts of convalescent subjects31,32 (Table 1). In the smaller (DHMC, n = 26) cohort, mucosal samples were available, and responses in nasal wash and stool were also defined. Relative to SARS-CoV-2 naïve controls, elevated serum IgA and IgG, but not IgM responses to whole S of diverse endemic CoV were frequently observed (Figure 2). Among endemic CoV, elevated levels of OC43-specific responses in serum were most pronounced, but were largely restricted to the S2 domain and whole unstabilized S. Similarly, elevated IgG responses to other endemic CoV were observed in serum from convalescent subjects to the spike proteins of endemic CoV, but not to the S1 domain alone. Indeed, elevated serum responses specific to the S1 domain of these CoV were not observed for any antibody isotype or subclass (Figure 2B, Supplemental Figure 2).
These observations were validated in a larger (JHMI, n = 126) cohort of convalescent plasma donors: elevated levels of IgG and IgA, but not IgM, to diverse endemic CoV were observed in plasma (Figure 2B, Supplemental Figure 3). Elevated responses were most pronounced for β-CoV (OC43 and HKU1) but were present for α-CoV as well. The absence of elevated IgM responses among convalescent subjects is consistent with recall of class-switched antibodies as opposed to the de novo elicitation of cross-reactive antibodies. Additionally, differential reactivity profiles were observed across three different OC43 antigens: proline-stabilized OC43 S (S-2P), OC43 S, and OC43 S2. Whereas elevated responses to OC43 S-2P among convalescent subjects in either cohort were not observed compared to naïve controls, responses to both OC43 S and OC43 S2 (available only for the smaller cohort) were observed for IgG, but not IgM (Supplemental Figure 4).
This effect was not limited to serum; elevation of endemic CoV antibody responses was also observed in nasopharyngeal wash samples (Figure 2A,B, Supplemental Figure 2) and stool (Supplemental Figure 2). Again, these elevated responses were restricted to IgG and/or IgA isotypes, consistent with cross-reactive antibodies being boosted in a recall response during SARS-CoV-2 infection, and were observed for whole S or the S2 domain, but not S1.
OC43-Specific Antibody Responses in Acute infection
To better explore the kinetics of these responses, antibodies were measured in serum collected from a cohort of acutely infected subjects (n = 10, Table 1) two weeks after a positive PCR diagnosis afforded by twice weekly surveillance. At this early timepoint, subjects appeared to have robustly elevated levels of CoV-2 unstabilized S-specific IgG compared to naïve controls (Figure 3A). Again, these responses were observed for whole unstabilized S and the S2 domain, but not for the pre-fusion conformation stabilized S-2P form of OC43 spike. In contrast, only a minor increase in IgM recognizing unstabilized OC43 S was observed (<2-fold) (Figure 3A). Though cross-sectional in nature and reliant on a small number of subjects, these isotype profiles in acute infection nonetheless further suggest a recalled rather than novel response against the relevant cross-reactive epitopes in S2 and unstabilized spike, presumably in the post-fusion conformation.
In an effort to resolve whether elevated responses to endemic β-CoV were driven by boosting of pre-existing cross-reactive antibody lineages, a small cohort of subjects (n = 3) for which pre- and post-SARS-CoV-2 infection serum samples were available was analyzed (Table 1). The median increase in OC43 S2-specific IgA and IgG were 3.7- and 20-fold, respectively, between the pre- and post-infection timepoints (Figure 3B). In contrast, the median change in OC43 S2-specific IgM was within 2-fold, providing further support that the elevated responses to endemic β-CoV are due to boosting of pre-existing cross-reactive antibodies originally raised against homologous epitopes found in endemic CoV. Like in other convalescent subject cohorts, this boosting effect was observed in response to OC43 S and S2, but not the pre-fusion conformation stabilized S-2P form of OC43, which showed a median change less than 1.1-fold across isotypes. Boosting of IgG responses specific to the full-length spike of another endemic β-CoV, HKU1, was also observed. IgG and IgA but not IgM responses to whole unstabilized S, but not the S1 domain, of HKU1 were elevated following SARS-CoV-2 infection (Figure 3B, Supplemental Figure 5). Collectively, these results further suggest that recalled cross-reactive antibodies may be more likely to recognize both the better conserved S2 domain and post-rather than pre-fusion conformations of spike.
Correlations Between Antibody Responses to OC43 and SARS-CoV-2
We next examined correlative relationships between SARS-CoV-2- and OC43-specific IgG responses in convalescent cohorts. Whereas the magnitude of IgG responses to CoV-2 S2 were well correlated to those binding unstabilized OC43 S (RP = 0.61) and its S2 domain (RP = 0.45), they were less well correlated to responses to stabilized OC43 S-2P (RP = 0.29) (Figure 4A). Expanding this analysis to include additional CoV-2 specificities and other isotypes showed a hierarchy of correlative relationships (Figure 4B, Supplemental Figure 6). CoV-2-specific responses were better correlated to OC43 responses specific for unstabilized rather than stabilized spike. IgG responses showed stronger relationships than did IgA, which were in turn stronger than IgM. Consistent with the lack of elevated IgM to endemic CoV, these measures only rarely showed a statistically significant relationship with CoV-2-specific IgM responses. Lastly, among CoV-2 antigens tested, correlations with OC43 responses were strongest for the S2 domain and whole spike (S-2P), and weaker or absent for S1 and the RBD. Collectively, these correlative relationships are consistent with recall of class-switched antibodies recognizing shared epitopes in the S2 domain from prior endemic CoV exposure induced by SARS-CoV-2 infection.
Direct Evidence of Molecular Cross-Reactivity of SARS-CoV-2 and OC43-Specific Antibodies
To date, antibody cross-reactivity has been inferred from indirect evidence in the form of boosted responses to endemic CoV23-27, and more conclusively observed for select monoclonal antibodies that have been cloned and cross tested33,34. To better generalize the more definitive monoclonal studies, we sought to directly define the cross-reactivity of polyclonal antibodies raised following SARS-CoV-2 infection. Antibodies specific to stabilized CoV-2 S-2P, RBD, and S2 and unstabilized OC43 S were selectively purified from serum samples from 30 SARS-CoV-2-infected subjects with a range of disease severity and humoral response profiles. Unfractionated and antigen-specific antibodies eluted from affinity purification matrices presenting various epitopes of the CoV-2 spike protein and OC43 were then characterized to determine their cross-reactivity and isotype profiles. Successful affinity purification was confirmed by comparison of antigen-specific binding signal relative to total Ig levels for each isotype (Figure 5A, Supplemental Figures 7-18). While relative binding signal was elevated for each targeted antigen, it was not observed for control antigens such as influenza hemagglutinin (HA) (Figure 5A) or tetanus toxoid (Supplemental Figures 7-18).
A number of interesting differences in cross-reactivity for antibodies with differing isotypes and antigen-specificities were apparent. For example, OC43 S-specific fractions showed elevated recognition of CoV-2 S2 and S-2P, but not RBD (Figure 5A), demonstrating molecular cross-reactivity as a general feature of polyclonal IgG responses in convalescent subjects. Additionally, IgG in the CoV-2 S2-specific fractions cross-reacted robustly with OC43 S, but no cross-reactivity in the IgM fraction was observed, despite a robust IgM response to CoV-2 S2 (Figure 5B), suggesting that while SARS-CoV-2 infection elicits IgM responses to the S2 domain, these presumed de novo responses are not cross-reactive to OC43.
To begin to define the functional significance of these cross-reactive antibodies, eluted fractions were pooled and tested for neutralization potency. Whereas purified CoV-2 RBD and S-2P-specific fractions induced robust neutralization, neither purified OC43 S-specific nor SARS-CoV-2 S2 domain-specific serum antibodies had detectable neutralization activity (Figure 5C).
To quantify the extent of cross-reactivity across antigens and isotypes, the unfractionated and affinity purified fractions from each subject were titrated, and enrichment of antigen-specific antibodies was computed by calculating the difference between the area under the curve (AUC) for the best fit lines calculated from data points across subjects for each antibody specificity and each isotype in eluates versus unfractionated samples (Figure 6).
IgM, IgA, and IgG antibodies purified based on binding to SARS-CoV-2 RBD cross reacted with proline-stabilized spike (S-2P and S-6P) and, to a lesser extent, unstabilized S, but showed no evidence of cross-reactivity to SARS-CoV-2 S2 or to endemic CoV (Figure 6A). Because the overwhelming majority of neutralizing antibodies target the RBD, these results suggest that RBD-specific antibody responses raised by natural infection are unlikely to exhibit exceptionally broad CoV neutralization activity. Indeed, SARS-CoV-2 RBD-specific antibodies showed low, albeit detectable cross-reactivity toward even SARS-CoV S and S1.
Similarly, across all isotypes, antibody pools purified based on binding to stabilized SARS-CoV-2 S-2P showed the expected cross-reactivity to SARS-CoV-2 RBD, S2, unstabilized S, and stabilized S with six proline substitutions (S-6P) (Figure 6B). However, while both IgA and IgG components of SARS-CoV-2 S-2P-specific serum antibody fractions showed binding to OC43, cross-reactive IgM was not evident. Similarly, IgA and IgG, but not IgM, specific for SARS-CoV-2 S2 bound as expected to SARS-CoV-2 proteins as well as OC43, SARS-CoV, MERS, and other endemic CoV S2 domains (Figure 6C).
Lastly, OC43 S-specific antibodies were similarly purified and profiled for cross-reactivity. Despite the presumption that subjects had not experienced OC43 infection recently, IgM, IgA, and IgG to OC43 S were all successfully enriched (Figure 6D). IgA and IgM fractions showed cross-reactivity to a broad array of endemic CoV S but not S1 proteins. Likewise, OC43-specific IgA and IgG antibodies showed cross-reactivity to SARS-CoV-2 proteins containing the S2 domain. In contrast, OC43 S-specific IgM antibodies showed elevated recognition of diverse pandemic, pathogenic, and endemic S proteins in unstabilized forms, but a lack of recognition of stabilized prefusion conformations of CoV-2 spike, or the S1 domains of most other CoV tested. Given binding to whole S but not the S1 domain, this recognition is presumably principally driven via recognition of S2.
Immunization With Stabilized Spike Changes Cross-Reactivity Profiles
Based on the differential ability of affinity-purified antibodies to recognize stabilized and unstabilized forms of CoV spike proteins, we next sought to determine whether responses to natural infection differ from those that result from mRNA vaccination with stabilized SARS-CoV-2 spike (S-2P). To investigate this possibility, antibody responses were analyzed in two cohorts of pregnant women (Table 1) who were either infected (n = 38) or vaccinated (n = 50) during their third trimester, as well as a validation cohort of healthy vaccinated subjects (n = 37). Fascinatingly, despite inducing considerably greater levels of SARS-CoV-2-specific antibodies, immunization with stabilized spike (S-2P) failed to result in elevated IgG responses to endemic CoV among pregnant women (Figure 7). Similarly, elevated levels of OC43-reactive IgA responses were not observed among vaccinated mothers, and neither infected nor vaccinated mothers showed elevated levels of OC43-reactive IgM (Supplemental Figure 19). These differing immunogenicity profiles suggest that the presentation of native spike in the context of natural infection and stabilized spike in the context of mRNA vaccination are distinct. Lastly, to confirm this observation, responses, a validation cohort of mRNA-vaccinated healthy adults was evaluated (Table 1). In this cohort, while IgG responses to OC43 were statistically significantly elevated for S and S2 antigens, the effect was small (<2.5-fold) (Figure 7); elevated IgA or IgM responses were not observed (Supplemental Figure 19).
In sum, elevated IgA and IgG but not IgM responses to the endemic CoV OC43 were observed in five distinct cohorts, including acutely infected, cross-sectional, and longitudinal convalescent cohorts. In contrast, two cohorts of subjects immunized with stabilized spike in the form of mRNA-based vaccines showed no (IgA and IgM) or only a small (IgG only, and in only one cohort) elevation of binding to endemic CoV OC43. In combination with the lack of boosting of responses that cross-react with endemic CoV, the high levels of neutralization activity observed to result from vaccination35 suggest favorable antigenicity of the pre-fusion conformation of S, and that the costs of original antigenic sin might be avoided by immunogen design.
Discussion
Cross-reactivity of endemic CoV antibodies against SARS-CoV-2 is of potential clinical relevance. The original antigenic sin hypothesis suggests that pre-existing immunity results in the reactivation of a response to an earlier strain, as opposed to the formation of an unbiased response against the current strain. While this recall could offer temporal advantages, it may reduce the formation of neutralizing antibodies thereby dampening effective clearance of the novel virus. Seminal studies of responses to influenza, based on epidemiology, modeling, and repertoire profiling suggest that the antibodies generated from childhood exposure to influenza are ‘imprinted’ and exert a major influence on the nature of the antibody response elicited upon subsequent exposures in humans6-10. It has been suggested that this phenomenon may exist in SARS-CoV-2 infection28.
To better understand the role that prior exposure to CoV plays in antibody responses to SARS-CoV-2, we examined antigen-specific antibody responses against endemic and pandemic human CoV in subjects one month after SARS-CoV-2 infection. We measured heightened OC43-specific responses in SARS-CoV-2 convalescent subjects when compared to naïve controls and observed that the magnitude of these responses correlated with responses to SARS-CoV-2 spike protein and the S2 subdomain. Separately, we observed OC43-specific IgG, but not IgM, in acutely infected subjects two weeks post positive COVID-19 diagnosis, indicating that these elevated responses were not likely to result from newly-induced antibody lineages with cross-reactivity to OC43. Testing of serum samples pre- and post-SARS-CoV-2 infection reinforced the hypothesis that pre-existing clonal families were boosted by the SARS-CoV-2 infection.
By isolating antigen-specific antibodies from convalescent serum, we determined that antibodies targeting SARS-CoV-2 S-2P and S2, but not the RBD, are cross-reactive with the spike of OC43. Similarly, antibodies isolated using OC43 were cross-reactive with SARS-CoV-2 S-2P and S2 but not the RBD. These results indicate that the cross reactivity between these viral spike proteins is principally linked to the better-conserved S2 domain. Consistent with the observation that neutralizing antibodies tend to target the RBD of SARS-CoV-236, we determined that both OC43 S-specific and SARS-CoV-2 S2-specific antibodies were non-neutralizing.
Over the course of this work, differing degrees of cross-reactivity were observed depending on whether pre-fusion-stabilized or unstabilized forms of S were used. In a comparison of naturally-infected and vaccinated subjects representing exposure to unstabilized and proline-stabilized S antigens, respectively, the immunogens starkly contrasted. The immune response that resulted from vaccination with SARS-CoV-2 S protein stabilized in the pre-fusion conformation did not share the degree of cross-reactivity observed in those recovering from COVID-19. Reciprocally, elevated responses to stabilized OC43 were not observed in either infection or vaccination cohorts. Collectively, this data suggests both the critical importance of the conformational state of S to resulting humoral responses, but also the prevalence of epitopes that are not shared between pre and post-fusion or other conformations.
Whereas numerous studies have focused on the antibody responses elicited from vaccination towards the SARS-CoV-2 S protein and specifically the RBD, little is known about the effect previous exposure to endemic CoV may have on vaccination. However, our observations of limited boosting agree with two recent reports37,38. Evaluation of a cohort vaccinated with a non-stabilized form of S will be required to begin to elucidate whether this is a feature of the stabilization itself, the mRNA-based vaccine modality, or if some other aspect of the vaccine formulation is responsible. Finally, this evidence of cross-reactivity has implications for the design of vaccines targeting SARS-CoV-2 variants. There is a risk that subsequent antibody responses might disregard the novel epitopes in favor of those already more familiar, although it should be emphasized that cross-reactivity was not a deterrent to the development of multiple highly efficacious SARS-CoV-2 vaccines.
Limitations of this study include use of several small cohorts, some collected from distinct geographic locations, from subjects of varying age, who experienced differing disease severity, and who were not confirmed to lack a recent exposure to endemic CoV. While all natural infection cohorts showed evidence of boosting toward the endemic CoV OC43, it would be beneficial to survey boosting in cohorts that better span a range of disease severity and ages to examine whether those variables impact the observed boosting effect, and outcomes of infection. To the extent that question has been investigated, it appears that boosting of responses to endemic CoV may be associated with poorer responses to SARS-CoV-228.
Likewise, while both vaccination cohorts showed an absence of or reduction in boosting of cross-reactive responses, whether this absence is beneficial or detrimental cannot be determined from the data presented here. While cross-reactive antibodies were non-neutralizing, neutralization is not the sole mechanism by which antibodies confer protection. Cross-reactive antibodies could interact with receptors for the Fc region of antibodies found on the surface of innate immune cells and promote protective effector function activities including antibody-dependent cellular phagocytosis and antibody-dependent cellular cytotoxicity39-41. To this end, a subset of monoclonal antibodies isolated from SARS-CoV patients that were able to cross-react with, but not neutralize SARS-CoV-2 were able to confer protection in a mouse model42. Indeed, both passive transfer experiments of monoclonal antibodies43-47 and evaluation of polyclonal antibodies raised in the context of vaccination or infection48-50 have shown that effector mechanisms contribute to antiviral activity in vivo. This evidence extends beyond correlative observations45,48 to include mechanistic1 evidence of in vivo contributions via Fc sequence engineering to knockout or enhance these activities44,46, as well as on the basis of depletion of effector cells 46.
In sum, this study provides evidence that antibodies targeting OC43 are robustly boosted in response to SARS-CoV-2 infection but not vaccination with stabilized S, and that the S2 subdomain of the spike protein is likely responsible for triggering a recalled, IgG-dominated response. While non-neutralizing, the role of these cross-reactive antibodies in the context of infection is not yet known. Further work aimed at characterizing the effector function potential of these antibodies and understanding their role in vaccine-mediated protection could provide a more complete picture of the relative risks and benefits of this recall response relevant to vaccine design, particularly in the context of viral variants of concern.
Methods
Structure Visualization and Manipulation
The sequence alignments were performed using Geneious 2021.1.1. Sequences and coordinates for coronavirus spike proteins were retrieved from the Protein Data Bank (PDB) entries SARS-CoV-2 (PDB 6XKL), 229E (PDB 6U7H), OC43 (PDB 6OHW), NL63 (PDB 5SZS), HKU1 (PDB 5I08), and SARS-CoV-2 Closed (PDB 6×6P) (Supplemental Table 1). SARS-CoV-2 was structurally aligned to the other models by domain using the MatchMaker function with default parameters and visualized using Chimera version 1.1551. For structural characterization of conservation among the 60 complete genomes of the Coronaviridae suborder (https://www.ncbi.nlm.nih.gov/genomes/GenomesGroup.cgi?taxid=11118), Batch Entrez was used to find 585 associated proteins, which were then further down selected to spike proteins (N=56) and aligned using Clustal Omega. This alignment was used to render by conservation and visualized using ChimeraX version 1.252.
Human Subjects
Initial study cohorts comprised 126 adult subjects interested in donating COVID-19 convalescent plasma. All were diagnosed with SARS-CoV-2 infection by PCR-based assays of nasopharyngeal swab, met the standard eligibility criteria for blood donation and were collected in the Baltimore, MD and Washington DC area (Johns Hopkins Medical Institutions, JHMI cohort), as previously described32 and partially reported in Natarajan, et al.53, and 26 SARS-CoV-2 convalescent individuals from the Lebanon, New Hampshire area (Dartmouth-Hitchcock Medical Center, DHMC cohort) and partially reported previously in Butler, et al.31 SARS-CoV-2 infection status was confirmed in all subjects by nasopharyngeal swab PCR. Plasma (JHMI) or serum (DHMC) was collected approximately one month after symptom onset or first positive PCR test.
The acute infection cohort was comprised of ten subjects under the age of thirty with mostly asymptomatic infections diagnosed by positive results of PCR swabs in twice-weekly screening tests. Serum samples were collected two weeks following a positive test result. Three subsequently infected subjects over the age of sixty with mild symptoms for which pre-pandemic serum samples had been banked were evaluated pre-infection and at two- and ten-weeks post-infection were also evaluated.
To support comparisons of natural infection and vaccination, two cohorts of pregnant women were evaluated. The first was a cohort of 38 subjects collected in Belgium who tested positive for SARS-CoV-2 in the third trimester of pregnancy. The second was a cohort of 50 subjects who were vaccinated with the Pfizer/BioNTech (BNT162) SARS-CoV-2 vaccine in the third trimester of pregnancy in Israel. Maternal serum samples were collected at time of delivery. A final cohort of 37 healthy adults vaccinated with Pfizer/BioNTech (BNT162) (n = 35) or Moderna (n = 2) vaccines whose serum was collected approximately one week after the second vaccine dose was evaluated to confirm observations of the differing immunogenicity observed among pregnant vaccine recipients.
Negative controls included samples from 16 naïve subjects collected from the Hanover, New Hampshire area, and from 38 subjects collected pre-pandemic by a commercial vendor (BioIVT).
Human subject research was approved by the Johns Hopkins University School of Medicine’s Institutional Review Board, the Dartmouth-Hitchcock Medical Center, CHU St. Pierre, Hadassah Medical Center, and BioIVT clinical site Committees for the Protection of Human Subjects as described in Table 1. Participants provided informed written consent. Table 1 provides basic clinical and demographic information for each cohort.
Fc Array Assay
SARS-CoV-2 antigens (Supplemental Table 2), including spike protein in its trimeric and subdomain forms (i.e., S1, S2, RBD), endemic CoV, and the control antigens influenza hemagglutinin and tetanus toxoid were covalently coupled to Luminex Magplex magnetic microspheres by two-step carbodiimide chemistry as previously described54. Anti-isotype and subclass primary antibodies were used to quantify the total amount of each immunoglobulin isotype in a sample before (load) and after (eluate) antigen-specific antibody purification. The load was diluted in 1x PBS by 1:100 followed by seven five-fold serial dilutions. The eluate was diluted beginning from 1:10 with six five-fold serial dilutions. Antibody isotypes and subclasses were detected using R-phycoerythrin (PE) conjugated secondary Abs as previously described55. A FlexMap 3D array instrument was then used to measure the median fluorescence intensity (MFI) of the bead sets.
Antigen-Specific Antibody Purification
Human CoV antigens were covalently attached to magnetic Dynabeads (ThermoFisher, 65011) using carbodiimide chemistry and the manufacturer’s instructions. Briefly, 1 nanomole of antigen was coupled to 300 μL of beads using a solution of 10 mg/mL EDC and 10 mg/mL sulfo-NHS. Bead activation and antigen coupling took place at room temperature with end-over-end mixing for 30 minutes and overnight, respectively. Following washes, the beads were reconstituted to 150 μL in PBS-TBN (Teknova, P0220).
The initial DHMC and JHMI cohorts were downselected to a total of 30 subjects, 15 from each cohort, representing a range of responses. Briefly, the self-reported case severities provided by the subjects of the DHMC cohort were used to select 4 mild, 5 moderate, and 6 severe subjects’ samples. From the JHMI cohort, the subjects were ranked by their anti-spike titers and the five subjects closest to the 25th, 50th, and 75th percentiles were selected for purification.
In a non-binding, clear-bottom 96-well plate, 5 μL of beads were diluted with 5 μL of 1X PBS before adding 50 μL of serum or plasma to the well. The plate was covered and shook at 800 rpm for 2.5 hours at room temperature. Using a magnetic base insert from a plate washer, the beads were pulled down from suspension for 1 minute before decanting the waste. The beads were washed 3 times using 100 μL of PBS-TBN and 3 minutes of shaking for each wash. Following the third wash, the magnetic separation and decanting step was followed up by pipetting all residual buffer out of the wells while the plate sat on the magnetic base. Antibodies specific to the antigen found on the beads were eluted using 20 μL of 1% formic acid (pH 2.9). After incubating with shaking for 10 minutes at room temperature and separating the beads, the eluate was pipetted from each well and transferred to a plate with wells containing 8 μL of 0.5 M sodium phosphate.
Neutralization Assay
Neutralization was performed using a VSV-SARS-CoV-2 pseudovirus assay as previously described31,56. Briefly, serum or plasma samples were serially diluted two-fold starting from a 1:25 dilution and incubated at 37°C for one hour with VSV-SARS-CoV-2 pseudovirus. Virus-serum/plasma mixtures were then added to pre-plated 293T-ACE2-expressing target cells in white 96-well plates at a final volume of 100 μL per well and incubated at 37°C for 24 hrs. To test for neutralization by antibody eluates, samples were first concentrated 10-fold using Amicon Ultra 0.5mL centrifugal filter devices (molecular weight cut-off of 100K Da) followed by two-fold serial dilution from 1:30, incubation with VSV-SARS-CoV-2 S pseudovirus, and addition of the mixtures to 293T-ACE2 target cells as described. Luciferase activity was measured using the Bright-Glo system and percent neutralization determined relative to control wells consisting of 293T-ACE2 cells infected with the pseudovirus alone.
Data Analysis
Basic statistical data analysis and visualization or raw Fc Array data was performed using GraphPad Prism, with statistical tests described in each figure legend. Heatmaps were visualized using the ‘pheatmap’57 package with hierarchical cluster analysis58 defined using Manhattan distance59 in R version 3.6.1. Fc Array features were log transformed then scaled and centered by their standard deviation from the mean (z-score).
In order to quantify the enrichment of antigen-specific antibodies, individual features (antigen-detection pair) were plotted relative to the total immunoglobulin isotype in serum across a titration range for each sample. A generalized additive model (GAM)60 with a cubic spline basis was used to fit a smoothed curve to load and eluate sample data using the ‘ggplot2’ package in R61. To quantify the difference between the load and eluate curves, we fit a separate GAM curve to the differences of the predicted values from the load and eluate GAM curves over total antibody as previously suggested62. We approximated the area under the difference curve (AUC) using the trapezoidal rule63,64 in the ‘pracma’65 package in R.
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Author Contributions
J.L., E.M.B., A.A.R.T., A.D.R., J.N.B., D.W., T.G., A.M., and P.F.W. contributed samples and/or other materials. H.N., A.R.C., A.P.H. W.W-A., and R.I.C. collected experimental data. H.N., A.R.C., A.P.B., C.A.B., W.W-A., R.I.C., S.X., and J.A.W. performed data analysis. H.N., A.R.C., and A.P.H. drafted the manuscript. All authors reviewed and edited the manuscript. M.E.A. and P.F.W. conceived of and supervised the research.
Acknowledgements
The authors acknowledge the Immune Monitoring and Flow Cytometry Resource (IMFCSR) at the Norris Cotton Cancer Center at Dartmouth supported by NCI Cancer Center Support Grant 5P30 CA023108-41. We thank all participants who enrolled and laboratory staff who helped collect and process the samples. VSV pseudovirus expression plasmids were provided by Michael Letko (Rocky Mountain Laboratories), CoV S-2P, S-6P, and RBD-Fc expression constructs were provided by Jason McLellan (University of Texas at Austin), and fusion peptide was provided by Laura Walker and Mrunal Sakharkar (Adimab). The following reagent was produced under HHSN272201400008C and obtained through BEI Resources, NIAID, NIH: Spike Glycoprotein Receptor-binding domain (RBD) from SARS-Related Coronavirus 2, Wuhan-Hu-1 with C-Terminal Histidine Tag, Recombinant from Baculovirus, NR-52307. The following reagent was deposited by the Centers for Disease Control and Prevention and obtained through BEI Resources, NIAID, NIH: SARS Related Coronavirus 2, Isolate USA-WA1/2020, NR-52281. This work was supported in part by the Division of Intramural Research, National Institute of Allergy and Infectious Diseases, as well as extramural support from the National Institute of Allergy and Infectious Diseases (U19AI145825 to M.E.A, R01AI120938, R01AI120938S1, and R01AI128779 to A.A.R.T.), National Heart Lung and Blood Institute (K23HL151826 to E.M.B.), National Institute of General Medical Sciences (P20-GM113132 BioMT Molecular Tools Core).