
==== Front
iScience
iScience
iScience
2589-0042
Elsevier

S2589-0042(24)01399-3
10.1016/j.isci.2024.110174
110174
Article
Compartment-specific antibody correlates of protection to SARS-CoV-2 Omicron in macaques
Tong Xin 13
Wang Qixin 13
Jung Wonyeong 1
Chicz Taras M. 1
Blanc Ross 1
Parker Lily J. 1
Barouch Dan H. dbarouch@bidmc.harvard.edu
12∗
McNamara Ryan P. rpmcnamara@mgh.harvard.edu
14∗∗
1 Ragon Institute of Mass General, MIT, and Harvard, Cambridge, MA 02139, USA
2 Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA 02215, USA
∗ Corresponding author dbarouch@bidmc.harvard.edu
∗∗ Corresponding author rpmcnamara@mgh.harvard.edu
3 These authors contributed equally

4 Lead contact

04 6 2024
20 9 2024
04 6 2024
27 9 11017429 3 2024
9 5 2024
31 5 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Summary

Antibodies represent a primary mediator of protection against respiratory viruses. Serum neutralizing antibodies (NAbs) are often considered a primary correlate of protection. However, detailed antibody profiles including characterization of antibody functions in different anatomic compartments are poorly understood. Here we show that antibody correlates of protection against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) challenge are different in systemic versus mucosal compartments in rhesus macaques. In serum, NAbs were the strongest correlate of protection and linked to spike-specific binding antibodies and other extra-NAb functions that create a larger protective network. In bronchiolar lavage (BAL), antibody-dependent cellular phagocytosis (ADCP) proved the strongest correlate of protection rather than NAbs. Within BAL, ADCP was linked to mucosal spike-specific immunoglobulin (Ig)G, IgA/secretory IgA, and Fcγ-receptor binding antibodies. Our results support a model in which antibodies with different functions mediate protection at different anatomic sites.

Graphical abstract

Highlights

• Correlates of protection to SARS-CoV-2 Omicron are highly compartment specific

• Antibody effector functions are primary correlates of protection at infection site

• Mucosal boosting enhances IgA and functionally levered IgG in lower respiratory tract

Health sciences; Biological sciences

Subject areas

Health sciences
Biological sciences
Published: June 4, 2024
==== Body
pmcIntroduction

COVID-19 vaccines, which generate antibodies to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein, have shown remarkable success at attenuating severe disease. Neutralizing antibodies, which most commonly target the receptor binding domain (RBD) of spike, were identified as a correlate of protection against ancestral strains of SARS-CoV-2.1 However, as Omicron-lineage SARS-CoV-2 variants emerged, vaccine-/infection-acquired antibody neutralization was largely lost due to the high degree of antigenic shift within the RBD.2,3,4,5,6 Yet protection from disease in vaccinated individuals did not see a concomitant drop,7,8,9,10,11 signifying that immune mediators of protection other than neutralizing antibodies existed.

Beyond their capacity to neutralize, antibodies exert several non-neutralizing functions such as antibody-dependent opsinophagocytosis, antibody-dependent cellular cytotoxicity, and complement deposition.12 These functions are largely modulated by post-translational modifications to the crystallizable fragment (Fc) of antibodies which dictate their binding to Fc receptors (FcγR for immunoglobulin [Ig]G subclasses, FcαR for IgA subclasses, etc.) on the surface of immune cells. Previous reports have demonstrated that FcγR-binding antibodies can recognize highly diverged SARS-CoV-2 spikes and confer protection even when neutralization is lost.13,14,15 To that end, antibodies mediate protection against pathogens such as SARS-CoV-2 through a variety of functions.

It is unclear how antibody correlates of protection for COVID-19 are shaped in different anatomic compartments. In this study, we show the unexpected results that neutralizing antibodies are a strong correlate of protection in serum but are not a clear correlate of protection in mucosa. Instead, extra-neutralizing functions such as antibody-dependent cellular phagocytosis (ADCP) were the strongest correlate of protection in bronchiolar lavage (BAL) against SARS-CoV-2 Omicron challenge. This extra-neutralizing role was conserved across SARS-CoV-2 variant spikes, including the challenge strain. Our results support a model in which antibody correlates of protection against SARS-CoV-2 Omicron are different in different anatomic compartments.

Results

Identification of humoral features inversely correlated with SARS-CoV-2 viral loads by anatomic compartment

Binding antibodies, Fcγ-receptor (FcγR) binding antibodies, neutralization titers, and Fc effector functions from serum- and lower respiratory tract-resident antibodies were quantified through systems serology (Figure 1, S1, and S2). A composite multivariate probable least-squares regression (PLSR) model was built to identify antibody features correlated with protection against viral loads across treatment groups. This was done for antibodies in serum (Figure 2A) and BAL (Figure 2B). Viral loads scattered on the latent variable 1 axis (LV1), which accounted for 66% and 68% of the total variance explained in the serum and BAL humoral profiles, respectively.Figure 1 Full systems serology data array

Serum (left) and BAL (right) antibody samples were assayed at week 0, week 4, and week 15 post-bivalent booster for binding to the antigens listed on the top left. Antibody subclass, isotype, Fc-binding, and functional outputs are shown on the right column. Route of administration and value legend is shown on the bottom right (IM = intramuscual and IN = intranasal). Callout boxes for ordering of antigens are shown on the left. All results were Z scored to account for differences in outputs. See also Figures S1 and S2.

Figure 2 Defining antibody correlates of protection within serum and bronchiolar lavage

(A) Partial least-squares regression (PLSR) model of serum antibody features of vaccinated and boosted non-human primates (NHPs) inversely correlated with viral loads within the lower respiratory tract. Heatmap gradient of viral loads is shown on the right.

(B) PLSR model of BAL antibody features of vaccinated and boosted NHPs inversely correlated with viral loads within the lower respiratory tract. Heatmap gradient of viral loads is shown on the right.

(C) Top serum antibody correlates of protection in the PLSR model.

(D) Top BAL antibody correlates of protection in the PLSR model. Neutralizing antibody titer was not selected as a bona fide correlate in the BAL and was manually plotted in purple for comparison.

(E) Validation of PLSR-selected serum correlates of protection. Viral loads were inversely correlated with (left) neutralizing antibody titer and (right) antibody-dependent natural killer cell activation (ADNKA) as measured by macrophage inflammatory protein 1 beta (MIP1β) production. Spearman’s R values and multiple comparisons adjusted p-values are shown. Trendline is shown with shaded areas being the 95% confidence interval.

(F) Validation of PLSR-selected BAL correlates of protection. Viral loads were inversely correlated with (left) antibody-dependent cellular phagocytosis (ADCP) to BQ.1.1 spike (challenge strain), but not to neutralizing antibody titers. Spearman’s R values and multiple comparisons adjusted p values are shown only for the statistically significant ADCP. Trendline is shown with shaded areas being the 95% confidence interval. See also Figures S3 and S4.

Distinguishing antibody features driving protection by compartment were identified. In the serum, neutralizing antibody titers, FcγR-binding antibodies, antibody-dependent natural killer cell activations (ADNKAs), and ADCP were all significant correlates of protection (Figure 2C, light blue bars). Interestingly, features driving protection in the BAL were focused on ADCP, IgG binding antibodies to Omicron XBB.1.5 spike, and antibody-dependent complement deposition (ADCD) (Figure 2D, light blue bars). Since neutralizing antibodies were not selected as a driver of protection in BAL, we separately added it for comparison relative to the other features (Figure 2D, purple bar).

To confirm these results, we plotted correlations of selected features with viral loads based on compartments. As expected, serum selected features at 4 and 15 weeks post-boost correlated with protection as defined by viral loads (Figures 2E and S3A). Likewise, BAL selected features correlated with protection at 4 and 15 weeks post-boost (Figures 2F left and S3B). A correlation could not be evaluated with BAL neutralizing antibodies as the majority of animals did not have neutralizing antibodies, and thus more data are required to evaluate BAL neutralizing correlates of protection (Figure 2F, right).

To exclusively look at the humoral profile immediately before SARS-CoV-2 BQ.1.1 challenge, the top PLSR-selected features were correlated against viral loads at week 15 post-boost. Similar to our previous observations, serum neutralizing antibody levels were significantly correlated with a reduction in viral loads within the lower respiratory tract at this time point (Figure S3C). Within the BAL, ADCP against both wild-type (WT) and the challenge strain BQ.1.1 remained significantly correlated with reduced viral loads (Figure S3D). Therefore, whether combining weeks 4 + 15 as a total “area under the curve” humoral features or exclusively looking humoral profiles immediately before virus challenge, neutralizing antibodies remained strong predictor of protection within serum while ADCP was the strongest predictor of protection at the mucosa. PLSR models for both serum and BAL-resident selected features as driving protection were validated against permutated labels and random features (Figure S4).

BAL correlates of protection to SARS-CoV-2 are defined by binding, FcγR binding, and effector functions

To investigate how the humoral landscape operates at the systems level, we created constellation linkage profiles by compartment. To create these constellation linkage profiles by the compartment, humoral features driving protection in the specific compartment selected by PLSR were used as centering features, or nodes. Humoral features significantly correlated with these compartment-specific drivers of protection were then linked to the PLSR-selected feature (see STAR methods). Serum-resident antibody features driving protection were part of a large constellation of humoral features that included binding antibodies, FcγR-binding antibodies, and Fc-effector-mediated functions (Figure S5).

ADCP to BQ.1.1 spike was identified as a protective feature within the BAL in this unbiased initial analysis. Humoral features significantly correlating with ADCP to BQ.1.1 were identified which included IgG subclasses, IgA, secIgA, and FcγR-binding antibodies, as well as to other effector functions (Figure 3A). A separate correlate of protection in the BAL was total IgG, and a constellation of co-correlating features was mapped showing other binding IgG and FcγR-binding antibodies across variants of concern (VOC) (Figure 3B). The last key correlate of protection identified was ADCD to XBB.1.5 spike. As expected, this node correlated with ADCD to other spike variants including Omicron sublineages (Figure 3C). The full constellation for BAL antibody correlates of protection is shown in Figure S6. We chose three constellations based on the function of the PLSR-selected feature (ADCP, ADCD, and Total IgG). We elected to use ADCP to BQ.1.1 as it was the challenge strain and thus represented the most relevant ADCP output. PLSR-selected features in the serum were not highlighted as they were decompartmentalized; however, we did note that several humoral features identified as protective features were present in both compartments.Figure 3 Antibody network of strongest correlates of protection in the BAL

(A) ADCP to the challenge strain BQ. 1.1 spike (purple node), which was identified through PLSR as a correlate of protection against viral loads, was centered and co-correlating antibody features (R > 0.6, FDR p value <0.05) taken from the BAL are shown.

(B) Same as (A), but for total IgG to XBB. 1.5 spike.

(C) Same as (A), but for ADCD to BA.5 spike. Blue circles correspond to the selected features in Figure 1. Purple circles correspond to BQ.1.1 spike, which was the challenge strain. Legend is shown at the bottom right. Although features such as ADCD to XBB.1.16 were identified as protectors in both compartments, it was not highlighted in BAL correlates of protection as it was not a PLSR-selected feature. See also Figures S5 and S6.

For the BAL, neutralizing antibodies were neither selected as a driving feature of protection by the PLSR nor linked to features that were as shown in the constellation network. This is in stark contrast to the serum correlates of protection where neutralizing antibodies were selected by the PLSR and were thus a node in the serum constellation network. Within the serum, neutralizing antibodies were part of a complex network of co-correlating features, including to other PLSR-selected drivers of protection within the compartment.

Mucosal boosting increases serum- and lower respiratory tract-resident humoral responses to divergent spikes

To define if correlates of protection were influenced by route of vaccine booster, we analyzed the identified correlates of protection using the delivery site as a variable. For ADNKA responses, non-boosted non-human primates (NHPs) showed no changes to any spikes at any time points. Intramuscular (IM)-boosted (purple) NHPs showed significantly enhanced serum ADNKA to BQ.1.1 and XBB.1.5 spikes. Intranasal (IN)-boosted (orange) NHPs likewise showed enhanced ADNKA responses but were significant for BA.5 and XBB.1.5 spike. Mucosal-boosted (blue) NHPs showed significant ADNKA responses to BA.5, BQ.1.1, and XBB.1.5 spike (Figure 4A). Significant expansions of ADNKA to WT spike were not observed for any groups, likely indicating that existing profiles to ancestral spike from previous vaccinations were still present. Binding IgG within the serum to the spikes strongly correlated with ADNKA. We thus assayed for total IgG binding to these spikes from the serum of the boosted NHPs. We found that total IgG to spike variants significantly increased for IM- and mucosal-boosted NHPs, but not for IN-boosted NHPs, after multiple comparisons adjustments (Figure 4B).Figure 4 Mucosal boosting enhances serum and lower respiratory tract humoral responses to various SARS-CoV2 VOCs

(A) Post-booster fold enhancements of serum-resident antibody-dependent natural killer cell activation (ANDKA) to the indicated spike variants as quantified by macrophage inflammatory protein 1 beta (MIP1β) production. ADNKA was selected as a key correlate of protection in Figure 1. Individual data points as well as moving averages are shown; the solid line is the mean and the shaded regions are the 95% confidence intervals.

(B) Same as (A), but for the serum-resident networked feature of total IgG to the indicated spikes variants. Color scheme legend is shown at the bottom for the serum responses.

(C) Post-booster fold enhancements of BAL-resident antibody-dependent cellular phagocytosis (ADCP) by monocytes to the indicated spike variants. ADCP was selected as a key correlate of protection in Figure 1.

(D) Same as (C), but for the BAL-resident networked feature of secretory IgA to the indicated spike variants. Color scheme legend is shown at the bottom for the BAL responses. For all plots, fold enhancements were relative to the mean value of the group at the pre-boost time point, week 0. Shown here are the best-fit models of matched responses within 95% CI in regions shaded in corresponding colors. ∗ = p < 0.05, ∗∗ = p < 0.01; Wilcoxon test followed by FDR correction. See also Figures S7 and S8.

ADCP to several spikes was the strongest correlate of protection in the BAL. Similar to serum responses, we analyzed how ADCP responses were shaped by booster delivery site. Interestingly, mucosal-boosted (blue) NHPs showed the strongest and most consistent ADCP to spikes tested from BAL samples. IM-boosted (purple) NHPs also showed some significant ADCP expansions, but the magnitude of responses was much lower than mucosal-boosted NHPs. While some ADCP was observed from BAL samples from IN-boosted NHPs, only one spike showed significance after multiple test corrections (Figure 4C).

Within the mucosa, secIgA was identified as a significant co-correlate of protection along with ADCP of BQ.1.1 spike. Therefore, to validate our constellation analysis, we plotted responses of secIgA in the boosted NHPs. Similar to ADCP, secIgA was significantly boosted in NHPs that received a mucosal boost at 4 weeks post-boost. These secIgA levels waned but remained elevated for the mucosal-boosted NHPs at 15 weeks post-boost. IM-boosted NHPs showed low secIgA induction to some spikes, while IN-boosted NHPs did not show any significant secIgA increases for any spike assayed (Figure 4C). Therefore the correlation between secIgA and the identified BAL correlates of protection was strongly driven by booster site.

Other serum and BAL humoral profiles showed similar results, including functional assays (Figures S7 and S8). Within the BAL, mucosal-delivered boosts consistently yielded higher antibody levels and functional outputs. Collectively, these results confirm our machine learning approaches of classifying correlates of protection by compartment. These results are strongly influenced by the route of boosting.

Mucosal boosting results in tight correlations between FcγR binding and effector functions

Given that effector-mediated functions such as ADCP, antibody-dependent neutrophil phagocytosis (ADNP), and ADNKA were disproportionally observed in mucosal-boosted NHPs, we asked if these functions were indeed correlated with FcγR-binding antibodies, particularly within the BAL. We conducted this analysis using week 15 time points only from BAL as this captured the humoral landscape and linkages at the time and site of challenge. Non-boosted NHPs showed non-specific and variable degrees of correlations between FcγR-binding antibodies and effector functions across VOC spikes, including the challenge strain BQ.1.1 (Figure 5A). Some scattered and loose correlations were identified for NHPs that were boosted through IN or IM within the BAL, but the overall architectures were uncoordinated (Figures 5B and 5C).Figure 5 Mucosal boosting induces a highly correlated network of FcγR-binding antibodies and effector functions in the BAL

(A) Correlation heatmaps were generated for FcγR-binding antibodies and effector-mediated functions at week 15 for the antigens used in this study for NHPs who received no boost. Callout boxes for the antigens used for the specific assays (FcγR-binding, ADCP, ADNP, ADCD, and ADNKA) and their orders are shown.

(B) Same as (A), but for NHPs at week 15 that received an IM boost.

(C) Same as (A), but for NHPs at week 15 that received an IN boost.

(D) Same as (A), but for NHPs at week 15 that received a mucosal boost. Heatmap legend for the correlation coefficient is shown in the bottom right.

In stark contrast to the non-boosted, IM-, and IN-boosted NHPs, those that were mucosal boosted showed tight correlations across FcγRs and with functional outputs within the BAL immediately before challenge. This tight degree of correlation was almost completely ambivalent to the spike variant (Figure 5D). The sole exception was some natural killer (NK) cell readouts. However, overall, the FcγR-binding and functional output profile within the BAL were highly coordinated in mucosal-boosted NHPs. This provides further support that the humoral architecture within the BAL is highly leveraged for non-neutralizing functions that are linked with protection against viral loads.

Taken together, FcγR-binding antibodies and non-neutralizing functions are highly influenced by the route of booster administration. Moreover, humoral features driving protection within the BAL can act as a well-coordinated network with broad breadth of binding and functionality.

Discussion

In this study, we report that antibodies with distinct functions in different anatomic compartments are correlates of protection against SARS-CoV-2 Omicron challenge in rhesus macaques. Our data support a model whereby mucosal antibody correlates of protection in the lower respiratory tract are primarily functional antibodies that mediate opsinophagocytosis, while serum antibody correlates of protection are primarily neutralization, binding, and non-neutralizing functions. Our data suggest that antibodies with Fc effector functions may be more important than currently appreciated at mucosal surfaces.

Previous work demonstrated that neutralizing antibodies were a correlate of protection against ancestral SARS-CoV-2 across several clinical trials.1 As SARS-CoV-2 adapted itself to the human population, VOCs emerged and neutralizing antibody capacity was progressively lost.16,17,18 This did not translate to a loss of clinical protection,9,10,11 indicating that correlates beyond neutralization existed. Previous reports have shown that FcγR-binding antibodies and their corresponding effector functions were required for protection against antigenically diverged spikes.13,19

We employed a systems serology approach to arrive at the conclusion that humoral profiles that correlate with protection are different in different anatomic compartments. Deep antibody profiling including a comprehensive analysis of antibody functions revealed that neutralizing antibody titers are the strongest correlate of protection in serum. This is in agreement with several reports showing that serum-resident neutralizing antibodies are strongly induced after vaccination and/or infection and are a correlate of protection.20,21,22,23,24,25,26,27,28 Our approach extends these findings to show that antibody correlates of protection are different at distinct anatomic sites. We show that antibody correlates of protection at the mucosa are enriched for effector functions such as opsinophagocytosis. These extra-neutralizing functions were most tightly linked to mucosa IgG and IgA/secIgA. Stimulation of IgA by different COVID-19 vaccine formulations has been shown.29 Our study shows that IgA can be stimulated by distinct delivery sites of the same vaccine formulation. IgA is known to be a potent neutralizer30,31,32 and a strong driver of opsinophagocytic function.33 Due to its polyfunctionality and mucosal localization, vaccine formulations and platforms have sought to enhance IgA responses.29,34,35,36,37 Our study supports the notion that vaccine delivery route can greatly impact mucosal IgA induction.

Much recent work has focused on characterizing mucosal protection against COVID-19 and other respiratory pathogens. This work extends on previous observations by McMahan et al.,28 Ying et al.,38 Hassan et al.,39 and other groups that have shown that targeted boosting at the mucosa can elicit superior protection against COVID-19.34,35,36,37,40 In our study, we show that this enhanced protection against COVID-19 is largely attributable to effector functions, particularly ADCP, at mucosal sites. This is in agreement with previous work by others that have shown that FcγR-binding antibodies were key determinants of protection against antigenically diverged spikes such as the Omicron sublineages and other sarbecoviruses.13,14,15,19 Also, recently Tong et al. showed that intranasal-delivered influenza vaccinations yielded a stronger effector function per antibody unit than intramuscular vaccination.41 It is well established that ADCC plays a key role in protection against influenza.42,43,44

Previous work has demonstrated that antibody recall responses to challenge are skewed to effector functions.45 These non-neutralizing responses were correlated with viral clearance shortly after infection.46 In the present study, we identify that protection in the BAL is linked with non-neutralizing functions post-boost. How long these effector functions persist after challenge and/or vaccination within the BAL is unclear. Future work at characterizing compartment-specific waning of effector functions similar to those published on waning of neutralization47,48,49,50 is needed.

In conclusion, our study demonstrates that antibody correlates of protection may be different in different anatomic compartments. In the lower respiratory tract, antibody Fc effector functions leveraged by IgG and IgA drive protection against SARS-CoV-2, whereas, in the serum, neutralizing antibodies drive protection. Further work characterizing how antibodies’ roles are influenced by their compartment can lead to vaccination strategies conferring multiple layers of protection. This is particularly important for emerging infectious diseases such as SARS-related coronaviruses.

Limitations of the study

Our analysis to conclude that distinct compartments have unique humoral architecture relied on antibody binding, cellular response, and neutralization assays. Each of these has intrinsic limits of detection/sensitivity. For example, several animals had detectable neutralizing antibodies within the BAL and no viral loads; however, most did not have detectable neutralizing antibodies, while IgG and secIgA were readily detectable. We therefore cannot discount that the IgG and IgA present within the BAL may play a stronger role in neutralization. ADNKA is limited in its sensitivity and is a readout of low-affinity FcγR-signaling, namely through FcγIIIA (CD16a). How circulating antibody concentrations in the serum and BAL can activate ADNKA at different concentrations is unclear. Further work on harmonizing these distinct antibody profiling approaches is needed.

The use of human cells to characterize effector functions from NHPs has been done previously by others and us.29,51,52,53,54 How compartmentalization can affect the expression and/or induction of Fc receptors and their signaling is still being explored. It is known that copy numbers of FcγRIIIA and FcγRIIIB in humans can be linked to autoimmune diseases.55,56 Thus there exists a steady-state equilibrium to their expression. Future work into characterizing how stimuli affect the surface expression of these Fc receptors across compartments is warranted.

STAR★Methods

Key resources table

REAGENT or RESOURCE	SOURCE	IDENTIFIER	
Antibodies	
	
Anti-CD107a	BD Biosciences	RRID:AB_396136	
Anti-CD3	BD Biosciences	RRID:AB_397038	
Anti-CD16	BD Biosciences	RRID:AB_396864	
Anti-CD56	BD Biosciences	RRID:AB_396853	
Anti MIP-1β	BD Biosciences	RRID:AB_393549	
Anti-IFNγ	BD Biosciences	RRID:AB_400425	
Anti-guinea pig complement C3 goat IgG fraction, FITC	MP Biomedicals	RRID:AB_2334913	
Goat anti-Mouse IgG Fc Cross-Adsorbed Secondary Antibody, PE	Thermo Fischer	RRID:AB_429715	
Mouse Anti-rhesus IgG-PE (SB108a)	Southern Biotech	RRID:AB_2796071	
Anti-rhesus IgG1 [7H11]	NHP Reagent Resource	RRID:AB_2819310	
Anti-rhesus IgG2 [3C10]	NHP Reagent Resource	RRID:AB_2895607	
Anti-rhesus IgG3 [2G11]	NHP Reagent Resource	RRID:AB_2819316	
Anti-rhesus IgG4 [7A8]	NHP Reagent Resource	RRID:AB_2819322	
Anti-rhesus IgA [9B9]	NHP Reagent Resource	RRID:AB_2819303	
Anti-Secretory IgA	Invitrogen	RRID:AB_931432	
Anti-rhesus IgM	Life Diagnostics	2C11-1-5	
	
Biological samples	
	
LowTox Guinea Pig Complement	CedarLane Labs	Cat # CL4051	
	
Chemicals, peptides, and recombinant proteins	
	
SARS-CoV-2 WT Spike	Sino Biological	40589-V08H4	
SARS-CoV-2 BA1 RBD	Sino Biological	40592-V08H129	
SARS-CoV-2 WT S2	Sino Biological	40590-V08B	
SARS-CoV-2 WT NTD	Sino Biological	40591-V49H	
SARS-CoV-2 Alpha Spike	Sino Biological	40589-V08B6	
SARS-CoV-2 Beta Spike	Sino Biological	40589-V08B7	
SARS-CoV-2 WT N	Sino Biological	40588-V08B	
SARS-CoV-2 Delta Spike	Sino Biological	40589-V08B16	
SARS-CoV-2 Gamma Spike	Sino Biological	40589-V08B10	
SARS-CoV-2 XBB1.5 Spike	Sino Biological	40589-V08H45	
SARS-CoV-2 BA1 Spike	Sino Biological	40589-V08H26	
SARS-CoV-2 HCMV gB	Sino Biological	10202-V08H1	
SARS-CoV-2 XBB1.16 Spike	Sino Biological	40589-V08H48	
SARS-CoV-2 Alpha RBD	Sino Biological	40592-V08H82	
SARS-CoV-2 SARS Spike	Sino Biological	40634-V08B	
SARS-CoV-2 Beta RBD	Sino Biological	40592-V08H59	
SARS-CoV-2 EG5 Spike	Sino Biological	40589-V08H55	
SARS-CoV-2 BA2.86 Spike	Sino Biological	40589-V08H58	
SARS-CoV-2 EG5 RBD	Sino Biological	40592-V08H151	
SARS-CoV-2 Delta RBD	Sino Biological	40592-V08H91	
SARS-CoV-2 Ebola Glycoprotein	Sino Biological	40459-V08H	
SARS-CoV-2 Gamma RBD	Sino Biological	40592-V08H86	
SARS-CoV-2 BA2 Spike	Sino Biological	40589-V08H28	
SARS-CoV-2 BA2 RBD	Sino Biological	40592-V08H123	
SARS-CoV-2 BA5 Spike	Sino Biological	40589-V08H32	
SARS-CoV-2 BA5 RBD	Sino Biological	40592-V08H131	
SARS-CoV-2 BQ1.1 Spike	Sino Biological	40589-V08H41	
SARS-CoV-2 BQ1.1 RBD	Sino Biological	40592-V08H143	
SARS-CoV-2 XBB1.5 RBD	Sino Biological	40592-V08H146	
SARS-CoV-2 XBB1.16 RBD	Sino Biological	40592-V08H136	
Rhesus soluble Fcγ2A-1	Duke University	Custom Order	
Rhesus soluble Fcγ2A-2	Duke University	Custom Order	
Rhesus soluble Fcy2A-3	Duke University	Custom Order	
Rhesus soluble Fcy2A-4	Duke University	Custom Order	
Rhesus soluble Fcγ3A-1	Duke University	Custom Order	
LC-LC-Sulfo-NHS Biotin	ThermoFisher	Cat # A35358	
Brefeldin A	Sigma Aldrich	Cat #B7651	
GolgiStop	BD Biosciences	Cat # 554724	
Streptavidin-R-Phycoerythrin	Prozyme	Cat # PJ31S	
	
Critical commercial assays	
	
Fix & Perm Cell Permeabilization Kit (Medium B)	ThermoFisher	GAS002S100	
Fix & Perm Cell Permeabilization Kit (Medium A)	ThermoFisher	GAS001S100	
EasySep™ Direct Human Neutrophil Isolation Kit	Stemcell Technologies	19666	
EasySep™ Human NK Cell Isolation Kit	Stemcell Technologies	17955	
NHS-Sulfo-LC-LC Kit	ThermoFisher	21435	
Zebra-Spin Desalting and Chromatography Columns	ThermoFisher	89882	
	
Experimental models: Cell lines	
	
THP-1 monocytes	ATCC	RRID: CVCL_0006	
Human Primary Natural Killer Cell Leukopack	StemCell	200–0092	
	
Software and algorithms	
	
GraphPad Prism 8	GraphPad Software, Inc.	RRID:SCR_002798	
iQue Forecyt 9.1	Sartorius	60028	
R Studio V 6.0	R Project for Statistical Computing	RRID:SCR_000432	
Flow Jo	BD Bioscience	RRID:SCR_008520	
MATLAB ver. R2019a	MathWorks	RRID:SCR_001622	
	
Other	
	
384-well HydroSpeed Plate Washer	Tecan	30190112	
iQue Screener Plus	Intellicyt/Sartorius	11811	
MagPlex Microspheres	Luminex MFG	MC12001-01 (Cataloged by region)	
Green Fluorescent Neutravidin Microspheres	ThermoFisher	Custom Synthesis	
Red Fluorescent Neutravidin Microspheres	ThermoFisher	Custom Synthesis	
Scarlet Fluorescent Neutravidin Microspheres	ThermoFisher	Custom Synthesis	
384-well HydroSpeed Plate Washer	Tecan	30190112	
MagPlex Microspheres	Luminex MFG	MC12001-01 (Cataloged by region)	
Green Fluorescent Neutravidin Microspheres	ThermoFisher	Custom Synthesis	
Red Fluorescent Neutravidin Microspheres	ThermoFisher	Custom Synthesis	
Scarlet Fluorescent Neutravidin Microspheres	ThermoFisher	Custom Synthesis	
Systems Serology Dataset	GitHub	BDRPM_20240422	
Existing dataset including neutralizing antibody titers and viral loads	McMahan et al.28		

Resource availability

Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Ryan P. McNamara (rpmcnamara@mgh.harvard.edu).

Materials availability

This study did not generate new unique reagents.

Data and code availability

• Raw systems serology data have been deposited at the Ragon Systems Serology homepage on GitHub under the accession: GitHub.com/RagonSystemSerology/DBRPM_20240422. This paper also analyzes existing, publicly available data. These accession numbers for the datasets are listed in the key resources table.

• This paper does not report original code. All scripts used for analysis can be found on the LoosC/systemsseRology on GitHub or can be obtained from the lead contact upon request

• Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

Experimental model and study participant details

The experimental model for this study in rhesus macaques (Macaca mulatta) has been published previously and serum and BAL were collected for secondary use.28 Briefly, 4–8 year old rhesus macaques were administered two doses of Ad26.COV2.S vaccine IM and one IM boost with either Ad26.COV2.S or Ad26.COV2.351 (Beta Spike). The NHPs were given boosters with a bivalent Ad26.COV2.S + Ad26.COV2.S.529 (BA.1 Spike) by IM, IN, or mucosal (intratracheal) route (6–8 animals/group). Serum and BAL collections were done at week 0 (pre-boost), week 4, and week 15 post-boost. Animals were subsequently challenged at week 16 with SARS-CoV-2 Omicron BQ.1.1 with 2E+6 PFU through intratracheal delivery. Viral loads were monitored and quantified in the lower respiratory tract (Figure S1A). As described in the previous study,28 all animal study protocols were designed and conducted in compliance with all relevant local, state, and federal regulations and were approved by the Bioqual Institutional Animal Care and Use Committee (IACUC).

Method DETAILS

Ig Subclassing/Isotyping and FcγR binding

Levels of antigen-specific antibody subclasses/isotypes and Fc-gamma receptor (FcγR) interaction were evaluated via the multiplexing Luminex microsphere-based assay, as previously described.57 Antigens of target were covalently linked to carboxyl group-labeled MagPlex microspheres (Luminex) through NHS-ester linkages using Sulfo-NHS and EDC (Thermo Fisher). Serum and BAL samples were diluted (serum isotypes/subclasses and FcγR binding: 1:250, BAL isotypes/subclasses and FcγR binding: 1:25) and added to the antigen-coupled microspheres to form the immune complexes in 384-well plates, and subsequently incubated at 4°C overnight, shaking at 750 rpm. After incubation, plates were washed with the washing buffer containing 0.1% BSA and 0.02% Tween 20 in PBS. Following the wash step, antibody isotype/subclass-specific mouse anti-rhesus antibodies (NHP Reagent Resource) were added to the immune complexes and incubated at room temperature for 1 h. Following a second wash step, the anti-mouse IgG Fc cross-adsorbed secondary antibody (PE, Thermo Fisher) was added to detect the anti-rhesus antibodies with fluorescence. For measurement of FcγR binding activities, Avi-tagged Rhesus macaque FcγRs (Duke Human Vaccine Institute) were biotinylated using BirA500 kit (Avidity) per manufacturer’s instructions and tagged with streptavidin-PE. The PE-labeled FcγR was subsequently incubated with the immune complexes for 2 h at room temperature. The plates were then washed and subject to Flow Cytometry measurements (iQue, Intellicyt) to determine the median fluorescence intensity (MFI). All Luminex experiments were conducted in duplicate, and the final results reported show the average values of the duplicates. The reagents and materials used are listed in the key resources table.

Antibody-dependent cellular phagocytosis and neutrophil phagocytosis

ADCP and ADNP experiments were performed as previously described.58,59 Briefly, antigen proteins of the target were biotinylated using the EZ-linkSulfo-NHS-LC-LC-Biotin kit (Thermo Fisher), then coupled to the fluorescent neutravidin beads (Thermo Fisher, F8776). The bead-antigen conjugates were incubated with diluted serum and BAL samples (serum: 1:100, BAL: 1:10) for 2 h at 37°C. The unbound antibody was removed by washing buffer. The immune complexes were then incubated overnight with cultured THP-1 cells (ADCP), or for 1 h with primary neutrophils isolated from human whole blood (ADNP) using negative selection (Stemcell). Treated THP-1 cells were subsequently washed and fixed in 4% paraformaldehyde (PFA), while the treated neutrophils were washed, stained for CD66b+ marker (Biolegend), and fixed in 4% (PFA) prior to flow cytometry analysis. A phagocytosis score for THP-1 or neutrophil was eventually determined as (% cells positive × Median Fluorescent Intensity of positive cells). Flow cytometry was performed with an iQue (IntelliCyt) instrument and population measurements were conducted using IntelliCyt ForeCyt (v8.1). The reagents and materials used are listed in the key resources table.

Antibody-dependent complement deposition (ADCD)

ADCD assays were designed and performed as previously described.60 Antigens of target were covalently linked to the carboxyl group-labeled MagPlex microspheres (Luminex) through NHS-ester linkages using Sulfo-NHS and EDC (Thermo Fisher) as described for Luminex. Diluted serum and BAL samples (serum: 1:50, BAL: 1:10) were incubated with coupled antigens for 2 h at 37°C to form immune complexes in 384-well plates. Plates were washed and incubated with lyophilized guinea pig complement (Cedarlane) diluted in gelatin veronal buffer with calcium and magnesium (Sigma Aldrich) for 20 min at 37°C. The deposition of C3 complement component was evaluated by an anti-guinea pig C3 FITC detection antibody (MpBio). Fluorescent intensity was acquired using an iQue Flow Cytometer (Intellicyt). The antibody-specific complement C3 deposition is calculated as the median fluorescence intensity of FITC. All ADCD experiments were conducted in duplicate, and final values were reported as average of the duplicates. The reagents and materials used are listed in the key resources table.

Antibody-dependent natural killer cell (NK) activation (ADNKA)

ADNKA assays were designed and performed as described previously.61 ELISA plates were coated 3 μg/mL of selected antigen and incubated at 4°C overnight. The coated plates were washed with PBS and blocked with 5% bovine serum albumin (BSA) for 2 h at 37°C. Natural Killer (NK) cells were isolated from Leukopaks (Stemcell Technologies) using EasySep Human NK Cell Isolation Kit (Stemcell Technologies). The isolated NK cells were incubated overnight at 37°C 5% CO2 in R10 (RPMI-1640 (Sigma Aldrich) media supplemented with 10% fetal bovine serum (FBS) (Sigma Aldrich), 5% penicillin/streptomycin (Corning, 50 μg/mL), 5% L-glutamine (Corning, 4 mM), 5% HEPES buffer (pH 7.2) (Corning, 50 mM) supplemented with 2 ng/mL IL-15. The ELISA plates were washed, and diluted serum and BAL samples (serum: 1:40, BAL: 1:10) were added to plates for 2 h at 37°C to form immune complexes. After wash, NK cells were added to plates at a concentration of 2.5E+5 cells/mL in R10 media supplemented with anti-CD107a–phycoerythrin (PE)–Cy5 (BD Biosciences, lot # 0149826, 1:1000 dilution), brefeldin A (10 μg/mL) (Sigma-Aldrich), and GolgiStop (BD Biosciences). The NK cells were incubated with immune complexes for 5 h at 37°C. The incubated NK cells were stained for cell surface markers with anti-CD3 Pacific Blue (BD Biosciences, clone G10F5)), anti-CD16 allophycocyanin (APC)-Cy5 (BD Biosciences, clone 3G8), and anti-CD56 PE-Cy7 (BD Biosciences, clone B159) for 15 min at room temperature. The washed NK cells were then fixed with PermA (Life Technologies), permeabilized with PermB (Life Technologies), and labeled with anti-MIP-1β PE (BD Biosciences) and anti-IFNγ FITC for 15 min at room temperature. Fluorescent intensity was measured using iQue Cytometer (Intellicyt). NK cells were gated as CD56+/CD16+/CD3-and the NK activation was evaluated as the percentage of NK cells positive for CD107a, IFNγ, or MIP-1b. All assays were performed with at least two healthy donors and the results shown here report the average of the donors. The reagents and materials used are listed in the key resources table. A gating strategy figure can be found in Figures S1B–S1H.

Pseudovirus neutralization assay

The Nab data represent lentivirus-based pseudovirus neutralization assays that have been previously published and incorporated in this study for comparison.28

Quantification and statistical analysis

Compartment-specific partial least-squares regression analysis (PLSR)

Partial Least Square Regression (PLSR) model was utilized to determine the best feature combination that describes peak viral load within the lower respiratory tract. Systems serology data including Ig isotype-, Ig subclass-, and FcγR-binding, as well as functional outputs ADCD, ADNP, and ADNKA, were regressed against peak viral loads within the lower respiratory tract of all animals. Neutralizing antibody titers from McMahon et al. were also incorporated into this analysis. Regressions were initially plotted as ambivalent to the route of booster delivery. The features that contributed most to the regression model were selected by LASSO (Least Absolute Shrinkage and Selection Operator). Features that were selected from more than 90% (Serum samples) or 80% (BAL samples) of 100 LASSO selections were finally selected for PLSR. PLSR model performance was determined by 5-fold cross-validations that were repeated 40 times. In addition, control models with random features or permuted output labels were built 25 times, whose accuracies were compared with the accuracy of the original model with 5-fold cross-validations. The LASSO-selected features were shown in the order of Variable Importance in Projection (VIP). For the BAL, neutralizing antibody levels regression against viral loads were manually plotted within the VIP. The compartment PLSR analysis was done for both time of challenge (week 15 post-booster only), as well as an “area under the curve” (AUC, 4 + 15 weeks post-booster).

Univariate correlations of PLSR

To further validate the PLSR, univariate correlations against viral loads for the selected humoral features were plotted. Viral loads are the log 10 of RNA copies/mL detected within the lower respiratory tract, and humoral features are plotted as the fold induction over baseline (week 0). This was done for both time of challenge (week 15 post-booster only) as well as an AUC (4 + 15 weeks post-booster). For the AUC univariate correlations, route of booster delivery is shown by the color of the symbol and time since booster is shown by the shape of the symbol. For the time of challenge univariate correlations, the route of booster delivery is shown by the color of the symbol. Spearman’s correlation coefficient (R) for each comparison is shown with an FDR adjusted p-value.

Compartment-specific constellation networks of protection

Humoral correlates with the PLSR-identified and validated correlates of protection were defined as having a Spearman’s correlation >0.6 with a false discovery rate (FDR) < 0.05 to the features. PLSR-selected humoral correlates of protection are shown in blue nodes, and features against the challenge strain BQ1.1 are shaded in purple. For these co-correlates networking (constellation networks), humoral features that were selected as significant drivers of protection by the PLSR were used as the centering features. Other humoral features that were correlated with that specific PLSR-selected feature (p < 0.05 after FDR correction, and an R value >0.6) were plotted as a linked constellation network. The length of the link does not contain any statistical information beyond achieving correlation with the above threshold.

Univariate comparisons

For univariate analysis, the “rstatix” R package was used, and two-sided Wilcoxon tests were performed to determine if data from different timepoints (week 0 vs. week 4, week 4 vs. week 15, and week 0 vs. week 15) significantly differed. p values were then corrected for multiple comparisons through FDR correction. Resulting p-values <0.05 were considered as significant. For all groups, baseline antibody levels were standardized to 1 to quantify fold changes. Moving averages from week 0, week 4, and week 15 post-booster were modeled and plotted showing the mean and 95% confidence intervals. Individual data points are also shown.

Supplemental information

Document S1. Figures S1–S8

Acknowledgments

We thank Terry and Susan Ragon and Mark and Lisa Schwartz for their contributions to the Ragon Institute of Mass General, MIT, and Harvard. This study was supported by the grants CA260476 (D.H.B.), 1P01AI165072-01 (R.P.M.), and 2U19AI135995-06 (R.P.M.) and the Bill and Melinda Gates Foundation Global Health and Vaccine Accelerator Program (GH-VAP) INV-001650.

Author contributions

X.T. and Q.W. contributed equally as co-first authors to this manuscript. X.T., Q.W., D.H.B., and R.P.M. conceptualized the study. X.T., Q.W., T.M.C., R.B., and L.J.P. performed the wet-bench systems serology experiments. Q.W., W.J., and R.P.M. performed bioinformatic systems serology analyses. X.T., Q.W., W.J., D.H.B., and R.P.M. analyzed all the data and wrote the manuscript. D.H.B. and R.P.M. provided supervision and funding acquisition. All authors reviewed and edited the manuscript.

Declaration of interests

The authors declare no competing interests.

Supplemental information can be found online at https://doi.org/10.1016/j.isci.2024.110174.
==== Refs
References

1 Gilbert P.B. Donis R.O. Koup R.A. Fong Y. Plotkin S.A. Follmann D. A Covid-19 Milestone Attained - A Correlate of Protection for Vaccines N. Engl. J. Med. 387 2022 2203 2206 10.1056/NEJMp2211314 36507702
2 Tegally H. Moir M. Everatt J. Giovanetti M. Scheepers C. Wilkinson E. Subramoney K. Makatini Z. Moyo S. Amoako D.G. Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa Nat. Med. 28 2022 1785 1790 10.1038/s41591-022-01911-2 35760080
3 Viana R. Moyo S. Amoako D.G. Tegally H. Scheepers C. Althaus C.L. Anyaneji U.J. Bester P.A. Boni M.F. Chand M. Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa Nature 603 2022 679 686 10.1038/s41586-022-04411-y 35042229
4 Cao Y. Wang J. Jian F. Xiao T. Song W. Yisimayi A. Huang W. Li Q. Wang P. An R. Omicron escapes the majority of existing SARS-CoV-2 neutralizing antibodies Nature 602 2022 657 663 10.1038/s41586-021-04385-3 35016194
5 Planas D. Saunders N. Maes P. Guivel-Benhassine F. Planchais C. Buchrieser J. Bolland W.H. Porrot F. Staropoli I. Lemoine F. Considerable escape of SARS-CoV-2 Omicron to antibody neutralization Nature 602 2022 671 675 10.1038/s41586-021-04389-z 35016199
6 Pajon R. Doria-Rose N.A. Shen X. Schmidt S.D. O'Dell S. McDanal C. Feng W. Tong J. Eaton A. Maglinao M. SARS-CoV-2 Omicron Variant Neutralization after mRNA-1273 Booster Vaccination N. Engl. J. Med. 386 2022 1088 1091 10.1056/NEJMc2119912 35081298
7 Arbel R. Peretz A. Sergienko R. Friger M. Beckenstein T. Duskin-Bitan H. Yaron S. Hammerman A. Bilenko N. Netzer D. Effectiveness of a bivalent mRNA vaccine booster dose to prevent severe COVID-19 outcomes: a retrospective cohort study Lancet Infect. Dis. 23 2023 914 921 10.1016/S1473-3099(23)00122-6 37062302
8 Altarawneh H.N. Chemaitelly H. Ayoub H.H. Tang P. Hasan M.R. Yassine H.M. Al-Khatib H.A. Smatti M.K. Coyle P. Al-Kanaani Z. Effects of Previous Infection and Vaccination on Symptomatic Omicron Infections N. Engl. J. Med. 387 2022 21 34 10.1056/NEJMoa2203965 35704396
9 Lauring A.S. Tenforde M.W. Chappell J.D. Gaglani M. Ginde A.A. McNeal T. Ghamande S. Douin D.J. Talbot H.K. Casey J.D. Clinical severity of, and effectiveness of mRNA vaccines against, covid-19 from omicron, delta, and alpha SARS-CoV-2 variants in the United States: prospective observational study BMJ 376 2022 e069761 10.1136/bmj-2021-069761
10 Nyberg T. Ferguson N.M. Nash S.G. Webster H.H. Flaxman S. Andrews N. Hinsley W. Bernal J.L. Kall M. Bhatt S. Comparative analysis of the risks of hospitalisation and death associated with SARS-CoV-2 omicron (B.1.1.529) and delta (B.1.617.2) variants in England: a cohort study Lancet 399 2022 1303 1312 10.1016/S0140-6736(22)00462-7 35305296
11 Shao W. Chen X. Zheng C. Liu H. Wang G. Zhang B. Li Z. Zhang W. Effectiveness of COVID-19 vaccines against SARS-CoV-2 variants of concern in real-world: a literature review and meta-analysis Emerg. Microbes Infect. 11 2022 2383 2392 10.1080/22221751.2022.2122582 36069511
12 Bowman K.A. Kaplonek P. McNamara R.P. Understanding Fc function for rational vaccine design against pathogens mBio 15 2024 e0303623 10.1128/mbio.03036-23
13 Mackin S.R. Desai P. Whitener B.M. Karl C.E. Liu M. Baric R.S. Edwards D.K. Chicz T.M. McNamara R.P. Alter G. Diamond M.S. Fc-γR-dependent antibody effector functions are required for vaccine-mediated protection against antigen-shifted variants of SARS-CoV-2 Nat. Microbiol. 8 2023 569 580 10.1038/s41564-023-01359-1 37012355
14 Bartsch Y.C. Tong X. Kang J. Avendaño M.J. Serrano E.F. García-Salum T. Pardo-Roa C. Riquelme A. Cai Y. Renzi I. Omicron variant Spike-specific antibody binding and Fc activity are preserved in recipients of mRNA or inactivated COVID-19 vaccines Sci. Transl. Med. 14 2022 eabn9243 10.1126/scitranslmed.abn9243 35289637
15 Kaplonek P. Fischinger S. Cizmeci D. Bartsch Y.C. Kang J. Burke J.S. Shin S.A. Dayal D. Martin P. Mann C. mRNA-1273 vaccine-induced antibodies maintain Fc effector functions across SARS-CoV-2 variants of concern Immunity 55 2022 355 365.e4 10.1016/j.immuni.2022.01.001 35090580
16 Miller J. Hachmann N.P. Collier A.R.Y. Lasrado N. Mazurek C.R. Patio R.C. Powers O. Surve N. Theiler J. Korber B. Barouch D.H. Substantial Neutralization Escape by SARS-CoV-2 Omicron Variants BQ.1.1 and XBB.1 N. Engl. J. Med. 388 2023 662 664 10.1056/NEJMc2214314 36652339
17 Wang Q. Iketani S. Li Z. Liu L. Guo Y. Huang Y. Bowen A.D. Liu M. Wang M. Yu J. Alarming antibody evasion properties of rising SARS-CoV-2 BQ and XBB subvariants Cell 186 2023 279 286.e8 10.1016/j.cell.2022.12.018 36580913
18 Hachmann N.P. Miller J. Collier A.R.Y. Ventura J.D. Yu J. Rowe M. Bondzie E.A. Powers O. Surve N. Hall K. Barouch D.H. Neutralization Escape by SARS-CoV-2 Omicron Subvariants BA.2.12.1, BA.4, and BA.5 N. Engl. J. Med. 387 2022 86 88 10.1056/NEJMc2206576 35731894
19 Adams L.E. Leist S.R. Dinnon K.H. West A. Gully K.L. Anderson E.J. Loome J.F. Madden E.A. Powers J.M. Schäfer A. Fc-mediated pan-sarbecovirus protection after alphavirus vector vaccination Cell Rep. 42 2023 112326 10.1016/j.celrep.2023.112326
20 Laing E.D. Weiss C.D. Samuels E.C. Coggins S.A. Wang W. Wang R. Vassell R. Sterling S.L. Tso M.S. Conner T. Durability of Antibody Response and Frequency of SARS-CoV-2 Infection 6 Months after COVID-19 Vaccination in Healthcare Workers Emerg. Infect. Dis. 28 2022 828 832 10.3201/eid2804.212037 35203111
21 Walls A.C. Sprouse K.R. Bowen J.E. Joshi A. Franko N. Navarro M.J. Stewart C. Cameroni E. McCallum M. Goecker E.A. SARS-CoV-2 breakthrough infections elicit potent, broad, and durable neutralizing antibody responses Cell 185 2022 872 880.e3 10.1016/j.cell.2022.01.011 35123650
22 Pérez-Then E. Lucas C. Monteiro V.S. Miric M. Brache V. Cochon L. Vogels C.B.F. Malik A.A. De la Cruz E. Jorge A. Neutralizing antibodies against the SARS-CoV-2 Delta and Omicron variants following heterologous CoronaVac plus BNT162b2 booster vaccination Nat. Med. 28 2022 481 485 10.1038/s41591-022-01705-6 35051990
23 Cohen A.A. van Doremalen N. Greaney A.J. Andersen H. Sharma A. Starr T.N. Keeffe J.R. Fan C. Schulz J.E. Gnanapragasam P.N.P. Mosaic RBD nanoparticles protect against challenge by diverse sarbecoviruses in animal models Science 377 2022 eabq0839 10.1126/science.abq0839
24 Goldblatt D. Alter G. Crotty S. Plotkin S.A. Correlates of protection against SARS-CoV-2 infection and COVID-19 disease Immunol. Rev. 310 2022 6 26 10.1111/imr.13091 35661178
25 Corbett K.S. Nason M.C. Flach B. Gagne M. O'Connell S. Johnston T.S. Shah S.N. Edara V.V. Floyd K. Lai L. Immune correlates of protection by mRNA-1273 vaccine against SARS-CoV-2 in nonhuman primates Science 373 2021 eabj0299 10.1126/science.abj0299
26 Garcia-Beltran W.F. Lam E.C. Astudillo M.G. Yang D. Miller T.E. Feldman J. Hauser B.M. Caradonna T.M. Clayton K.L. Nitido A.D. COVID-19-neutralizing antibodies predict disease severity and survival Cell 184 2021 476 488.e11 10.1016/j.cell.2020.12.015 33412089
27 Liu J. Liu Y. Xia H. Zou J. Weaver S.C. Swanson K.A. Cai H. Cutler M. Cooper D. Muik A. BNT162b2-elicited neutralization of B.1.617 and other SARS-CoV-2 variants Nature 596 2021 273 275 10.1038/s41586-021-03693-y 34111888
28 McMahan K. Wegmann F. Aid M. Sciacca M. Liu J. Hachmann N.P. Miller J. Jacob-Dolan C. Powers O. Hope D. Mucosal boosting enhances vaccine protection against SARS-CoV-2 in macaques Nature 626 2024 385 391 10.1038/s41586-023-06951-3 38096903
29 Deng Y. Atyeo C. Yuan D. Chicz T.M. Tibbitts T. Gorman M. Taylor S. Lecouturier V. Lauffenburger D.A. Chicz R.M. Beta-spike-containing boosters induce robust and functional antibody responses to SARS-CoV-2 in macaques primed with distinct vaccines Cell Rep. 42 2023 113292 10.1016/j.celrep.2023.113292
30 Maurer M.A. Meyer L. Bianchi M. Turner H.L. Le N.P.L. Steck M. Wyrzucki A. Orlowski V. Ward A.B. Crispin M. Hangartner L. Glycosylation of Human IgA Directly Inhibits Influenza A and Other Sialic-Acid-Binding Viruses Cell Rep. 23 2018 90 99 10.1016/j.celrep.2018.03.027 29617676
31 Wang Z. Lorenzi J.C.C. Muecksch F. Finkin S. Viant C. Gaebler C. Cipolla M. Hoffmann H.H. Oliveira T.Y. Oren D.A. Enhanced SARS-CoV-2 neutralization by dimeric IgA Sci. Transl. Med. 13 2021 eabf1555 10.1126/scitranslmed.abf1555
32 Sterlin D. Mathian A. Miyara M. Mohr A. Anna F. Claër L. Quentric P. Fadlallah J. Devilliers H. Ghillani P. IgA dominates the early neutralizing antibody response to SARS-CoV-2 Sci. Transl. Med. 13 2021 eabd2223 10.1126/scitranslmed.abd2223
33 Bartsch Y.C. Cizmeci D. Kang J. Zohar T. Periasamy S. Mehta N. Tolboom J. Van der Fits L. Sadoff J. Comeaux C. Antibody effector functions are associated with protection from respiratory syncytial virus Cell 185 2022 4873 4886.e10 10.1016/j.cell.2022.11.012 36513064
34 Hemmi T. Ainai A. Hashiguchi T. Tobiume M. Kanno T. Iwata-Yoshikawa N. Iida S. Sato Y. Miyamoto S. Ueno A. Intranasal vaccination induced cross-protective secretory IgA antibodies against SARS-CoV-2 variants with reducing the potential risk of lung eosinophilic immunopathology Vaccine 40 2022 5892 5903 10.1016/j.vaccine.2022.08.049 36064667
35 Diallo B.K. Ní Chasaide C. Wong T.Y. Schmitt P. Lee K.S. Weaver K. Miller O. Cooper M. Jazayeri S.D. Damron F.H. Mills K.H.G. Intranasal COVID-19 vaccine induces respiratory memory T cells and protects K18-hACE mice against SARS-CoV-2 infection NPJ Vaccines 8 2023 68 10.1038/s41541-023-00665-3 37179389
36 Lund F.E. Randall T.D. Scent of a vaccine Science 373 2021 397 399 10.1126/science.abg9857 34437109
37 Adler J.M. Martin Vidal R. Langner C. Vladimirova D. Abdelgawad A. Kunecova D. Lin X. Nouailles G. Voss A. Kunder S. An intranasal live-attenuated SARS-CoV-2 vaccine limits virus transmission Nat. Commun. 15 2024 995 10.1038/s41467-024-45348-2 38307868
38 Ying B. Darling T.L. Desai P. Liang C.Y. Dmitriev I.P. Soudani N. Bricker T. Kashentseva E.A. Harastani H. Raju S. Mucosal vaccine-induced cross-reactive CD8+ T cells protect against SARS-CoV-2 XBB.1.5 respiratory tract infection Nat. Immunol. 25 2024 537 551 10.1038/s41590-024-01743-x 38337035
39 Hassan A.O. Shrihari S. Gorman M.J. Ying B. Yuan D. Raju S. Chen R.E. Dmitriev I.P. Kashentseva E. Adams L.J. An intranasal vaccine durably protects against SARS-CoV-2 variants in mice Cell Rep. 36 2021 109452 10.1016/j.celrep.2021.109452
40 Li B. Jiang A.Y. Raji I. Atyeo C. Raimondo T.M. Gordon A.G.R. Rhym L.H. Samad T. MacIsaac C. Witten J. Enhancing the immunogenicity of lipid-nanoparticle mRNA vaccines by adjuvanting the ionizable lipid and the mRNA Nat. Biomed. Eng. 2023 10.1038/s41551-023-01082-6 Online ahead of print
41 Tong X. Deng Y. Cizmeci D. Fontana L. Carlock M.A. Hanley H.B. McNamara R.P. Lingwood D. Ross T.M. Alter G. Distinct Functional Humoral Immune Responses Are Induced after Live Attenuated and Inactivated Seasonal Influenza Vaccination J. Immunol. 212 2024 24 34 10.4049/jimmunol.2200956 37975667
42 Boudreau C.M. Alter G. Extra-Neutralizing FcR-Mediated Antibody Functions for a Universal Influenza Vaccine Front. Immunol. 10 2019 440 10.3389/fimmu.2019.00440 30949165
43 Jegaskanda S. Vanderven H.A. Tan H.X. Alcantara S. Wragg K.M. Parsons M.S. Chung A.W. Juno J.A. Kent S.J. Influenza Virus Infection Enhances Antibody-Mediated NK Cell Functions via Type I Interferon-Dependent Pathways J. Virol. 93 2019 e02090-18 10.1128/JVI.02090-18
44 Jegaskanda S. Weinfurter J.T. Friedrich T.C. Kent S.J. Antibody-dependent cellular cytotoxicity is associated with control of pandemic H1N1 influenza virus infection of macaques J. Virol. 87 2013 5512 5522 10.1128/JVI.03030-12 23468501
45 Dugan H.L. Stamper C.T. Li L. Changrob S. Asby N.W. Halfmann P.J. Zheng N.Y. Huang M. Shaw D.G. Cobb M.S. Profiling B cell immunodominance after SARS-CoV-2 infection reveals antibody evolution to non-neutralizing viral targets Immunity 54 2021 1290 1303.e7 10.1016/j.immuni.2021.05.001 34022127
46 McNamara R.P. Maron J.S. Boucau J. Roy V. Webb N.E. Bertera H.L. Barczak A.K. Positives Study Staff T. Franko N. Logue J.K. Anamnestic humoral correlates of immunity across SARS-CoV-2 variants of concern mBio 14 2023 e0090223 10.1128/mbio.00902-23
47 Tong X. McNamara R.P. Avendaño M.J. Serrano E.F. García-Salum T. Pardo-Roa C. Bertera H.L. Chicz T.M. Levican J. Poblete E. Waning and boosting of antibody Fc-effector functions upon SARS-CoV-2 vaccination Nat. Commun. 14 2023 4174 10.1038/s41467-023-39189-8 37443074
48 Moreira E.D. Kitchin N. Xu X. Dychter S.S. Lockhart S. Gurtman A. Perez J.L. Zerbini C. Dever M.E. Jennings T.W. Safety and Efficacy of a Third Dose of BNT162b2 Covid-19 Vaccine N. Engl. J. Med. 386 2022 1910 1921 10.1056/NEJMoa2200674 35320659
49 Barda N. Dagan N. Cohen C. Hernán M.A. Lipsitch M. Kohane I.S. Reis B.Y. Balicer R.D. Effectiveness of a third dose of the BNT162b2 mRNA COVID-19 vaccine for preventing severe outcomes in Israel: an observational study Lancet 398 2021 2093 2100 10.1016/S0140-6736(21)02249-2 34756184
50 Levin E.G. Lustig Y. Cohen C. Fluss R. Indenbaum V. Amit S. Doolman R. Asraf K. Mendelson E. Ziv A. Waning Immune Humoral Response to BNT162b2 Covid-19 Vaccine over 6 Months N. Engl. J. Med. 385 2021 e84 10.1056/NEJMoa2114583
51 Barouch D.H. Stephenson K.E. Borducchi E.N. Smith K. Stanley K. McNally A.G. Liu J. Abbink P. Maxfield L.F. Seaman M.S. Protective efficacy of a global HIV-1 mosaic vaccine against heterologous SHIV challenges in rhesus monkeys Cell 155 2013 531 539 10.1016/j.cell.2013.09.061 24243013
52 Tuero I. Mohanram V. Musich T. Miller L. Vargas-Inchaustegui D.A. Demberg T. Venzon D. Kalisz I. Kalyanaraman V.S. Pal R. Mucosal B Cells Are Associated with Delayed SIV Acquisition in Vaccinated Female but Not Male Rhesus Macaques Following SIVmac251 Rectal Challenge PLoS Pathog. 11 2015 e1005101 10.1371/journal.ppat.1005101
53 Routhu N.K. Stampfer S.D. Lai L. Akhtar A. Tong X. Yuan D. Chicz T.M. McNamara R.P. Jakkala K. Davis-Gardner M.E. Efficacy of mRNA-1273 and Novavax ancestral or BA.1 spike booster vaccines against SARS-CoV-2 BA.5 infection in nonhuman primates Sci. Immunol. 8 2023 eadg7015 10.1126/sciimmunol.adg7015
54 Francica J.R. Flynn B.J. Foulds K.E. Noe A.T. Werner A.P. Moore I.N. Gagne M. Johnston T.S. Tucker C. Davis R.L. Protective antibodies elicited by SARS-CoV-2 spike protein vaccination are boosted in the lung after challenge in nonhuman primates Sci. Transl. Med. 13 2021 eabi4547 10.1126/scitranslmed.abi4547
55 Aitman T.J. Dong R. Vyse T.J. Norsworthy P.J. Johnson M.D. Smith J. Mangion J. Roberton-Lowe C. Marshall A.J. Petretto E. Copy number polymorphism in Fcgr3 predisposes to glomerulonephritis in rats and humans Nature 439 2006 851 855 10.1038/nature04489 16482158
56 Fanciulli M. Norsworthy P.J. Petretto E. Dong R. Harper L. Kamesh L. Heward J.M. Gough S.C.L. de Smith A. Blakemore A.I.F. FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity Nat. Genet. 39 2007 721 723 10.1038/ng2046 17529978
57 Brown E.P. Dowell K.G. Boesch A.W. Normandin E. Mahan A.E. Chu T. Barouch D.H. Bailey-Kellogg C. Alter G. Ackerman M.E. Multiplexed Fc array for evaluation of antigen-specific antibody effector profiles J. Immunol. Methods 443 2017 33 44 10.1016/j.jim.2017.01.010 28163018
58 Butler A.L. Fallon J.K. Alter G. A Sample-Sparing Multiplexed ADCP Assay Front. Immunol. 10 2019 1851 10.3389/fimmu.2019.01851 31456799
59 Karsten C.B. Mehta N. Shin S.A. Diefenbach T.J. Slein M.D. Karpinski W. Irvine E.B. Broge T. Suscovich T.J. Alter G. A versatile high-throughput assay to characterize antibody-mediated neutrophil phagocytosis J. Immunol. Methods 471 2019 46 56 10.1016/j.jim.2019.05.006 31132351
60 Fischinger S. Fallon J.K. Michell A.R. Broge T. Suscovich T.J. Streeck H. Alter G. A high-throughput, bead-based, antigen-specific assay to assess the ability of antibodies to induce complement activation J. Immunol. Methods 473 2019 112630 10.1016/j.jim.2019.07.002
61 Ackerman M.E. Moldt B. Wyatt R.T. Dugast A.S. McAndrew E. Tsoukas S. Jost S. Berger C.T. Sciaranghella G. Liu Q. A robust, high-throughput assay to determine the phagocytic activity of clinical antibody samples J. Immunol. Methods 366 2011 8 19 10.1016/j.jim.2010.12.016 21192942
