SARS-CoV-2 specific memory B-cells from individuals with diverse disease severities recognize SARS-CoV-2 variants of concern

Human Research Ethics

The study was reviewed and approved by the Oregon Health & Science University Institutional Review Board (IRB# 21230). Informed consent was obtained from subjects on initiation of their participation in the study.

Human subjects

Samples (n=24) were collected from subjects with confirmed COVID-19 infection (OHSU IRB# 21230). Forty mL of whole blood was collected for PBMCs and plasma (BD Vacutainer® Lavender Top EDTA Tubes). Whole blood was spun @ 1,000 x g and the plasma fraction was collected and stored at -80 until use. PBMCs were isolated using Sepmate tubes (Stemcell) and stored in liquid nitrogen until needed.

Human MBC limiting dilution analysis

PBMCs were thawed and resuspended in LDA media, RPMI 1640 (Gibco), 1×Anti-Anti (Corning), 1X non-essential amino acids (HyClone), 20 mM HEPES (Thermo Scientific), 50μMβ-ME, and 10% heat-inactivated FBS (VWR). Cells were cultured in serial 2-fold diluted doses (10 wells per dose), starting with 3-5 x 10⁵ PBMCs per well at the highest dose. In 96-well round-bottom plates in a final volume of 200μl per well. Cells were incubated with IL-2 (Prospec) 1000U/ml and R848 (InvivoGen) 2.5ug/mL.[8] To determine background absorbance values, supernatants were used from 8 wells unstimulated PBMCs only. Plates were incubated at 37 °C and 5% CO² for 7 days. Stimulation was determined by running total IgG ELISAs. Antigen-specific MBC frequencies were calculated by assaying LDA supernatants by antigen-specific ELISAs [9].

Antigen-specific ELISA for plasma and LDA

Ninety-six well (plasma, Corning-3590) or half-well ELISA plates (LDA, Greiner Bio-one) were coated with 100µL (plasma) or 50µL (LDA) of 1ug/mL antigen in PBS, recombinant RBD (provided by Dr. David Johnson), recombinant spike subunit 1 (40591-V08H, Sino Biological Inc.), and recombinant VoC-RBD. Plates were incubated overnight at 4°C, washed with PBS-T (0.05% Tween) and blocked for 1 hour with 5% milk prepared in PBS-T. For plasma, samples were serially 3-fold diluted in dilution buffer and 100uL was added to wells. Plates were incubated at room temperature (RT) for 1 hour. For LDA, 20μL of LDA supernatants were added to each well and incubated at RT for 1 hour. Plates were washed 4 times with wash buffer, and 100µL (plasma) or 50 μL (LDA) of 1:3000 dilution of anti-human IgG-HRP (BD Pharmingen, 555788) detection antibody was added and incubated at RT for 1 hour. Plates were washed 4 times with wash buffer, 100µL (plasma) or 50 μL (LDA) of colorimetric detection reagent containing 0.4□mg/ml o-phenylenediamine and 0.01% hydrogen peroxide in 0.05□M citrate buffer (pH 5) were added and the reaction was stopped after 20□minutes by the addition of 1□M HCl. Optical density (OD) at 492□nm was measured using a CLARIOstar ELISA plate reader. Plasma Ab endpoint titers were set as the lowest dilution with an OD 4-fold above PBS only wells. LDA wells were scored positive at ODs at least 2-fold above background (unstimulated PBMC wells).

MBC frequency

MBC precursor frequencies were calculated by the semi-logarithmic plot of the percent of negative cultures versus the cell dose per culture, as previously described [9] and frequencies were calculated as the reciprocal of the cell dilution at which 37% of the cultures were negative for antigen-specific IgG production. Rows which yielded 0% negative wells were excluded, since this typically resides outside of the linear range of the curve and artificially reduced the MBC precursor frequency. For subjects with low frequency of antigen-specific antibody secreting cells frequency was determined by number of positive wells divided by the total number of IgG positive secreting wells, multiplied by one million, giving a frequency per million PBMCs stimulated.

Generation of VoC-RBD and protein expression

To generate the triple mutant VoC-RBD we started with the WT RBD gene and performed site directed mutagenesis to make the following changes: N501Y, E484K, K417N using Q5 site-directed mutagenesis kit (NEB). Purified SARS-CoV-2 VoC-RBD protein was prepared as described previously [40-41]. Briefly, His-tagged VoC-RBD bearing lentivirus was produced in HEK 293-T cells and used to infect HEK 293-F suspension cells. The suspension cells were allowed to grow for 3 days, shaking at 37°C, 8% CO2. Cells were spun down, supernatant collected, sterile filtered, and purified by Ni-NTA chromatography. The purified protein was then buffer exchanged into PBS and concentrated.

Statistical methods

For values below the limit of detection a value of 49 was assigned for ELISA titers, and a value of <0.0000001 was assigned for MBC frequencies (corresponding to <1 SARS-CoV-2 antibody secreting cell per total PBMCs stimulated). Subjects who received vaccination prior to second draw were included in the longitudinal data, but excluded from all statistical analysis. Paired t-test was done to compare ELISA titers between RBD-WA-1 and RBD-VoC as well as used to compare log transformed specific MBC frequencies against WA-1 and VoC RBD. Spearman’s rank-order correlation was performed for all correlation analysis. All analysis was done in Prism version 9.1.1.

Univariate and multivariate analyses

We conducted univariate analyses with each pre-determined clinical covariate of interest and MBC frequency was log transformed. Each univariate analysis was conducted as follows: age (continuous), age (categorical, <50 vs. ≥50), disease severity (ordinal), disease severity (categorical, <5 vs. ≥5), and hospitalization status (yes vs. no). For multivariable models our first model included age (continuous), disease severity score (ordinal), and the interaction between age and disease severity. Contrast estimates were evaluated for each disease severity score interaction with age to determine if any effects within each score and mean age existed. Our second model included age, hospitalization status, and the interaction between hospitalization status and age. All statistical analyses were performed in SAS version 9.4.

COVID-19 cohort population

SARS-CoV-2 PCR positive subjects were recruited in Oregon under an OHSU IRB approved protocol (IRB# 21230). On enrollment, a comprehensive medical history was taken in addition to serum, plasma, and peripheral blood mononuclear cells (PBMC). When possible, one or more follow-up serial sample visits were arranged. Samples from 24 subjects were taken during the convalescent period (>28 days post-infection) (Table 1). Time post-infection, defined by months post-PCR diagnosis, ranged from 1-11.75 months. WHO Disease severity scores [42] (ranged from 1 (asymptomatic) to 6 (severe, requiring intubation and mechanical ventilation), with 2 asymptomatic subjects, 16 symptomatic/not hospitalized subjects, and 6 hospitalized subjects. Subjects were 46% female, with a median age of 60 (Table 1). Of the 24 subjects, two were vaccinated before their final time-point, providing two comparator samples for vaccine elicited boosting vs. natural decay in un-boosted subjects.

The magnitude, durability, and breadth of LLPC-derived antibodies

Forty plasma samples were collected from 24 subjects, with 15 subjects providing multiple timepoints. The VoC-RBD contained the canonical K417N, E484K, and Y501N mutations present in the currently circulating VoC. Y501N is present in the UK, SA, and P.1 strains and has been shown to increase transmission via enhanced binding to ACE-2 receptor [10]. E484K is present in 15% of strains circulating in the US (CDC reporting 5/21/2021). K417N and E484K are present in the B1.351 and P.1 strains [9, 25-26]. Antibody ELISA titers were below our lower-limit of detection for two asymptomatic study subjects at two separate timepoints (Fig.1 A, B) against full-length S1, RBD and VoC RBD. All but two symptomatic non-hospitalized cases had detectable antibody titers by ELISA against parental WA-1 RBD, but five symptomatic subjects had S1 ELISA titers at or below the limit of detection (at least one timepoint), and four subjects fell below the limit of detection against VoC RBD (at least one timepoint) (Figure 1A). All of the hospitalized subjects ELISA titers were above the limit of detection against S1, WA-1 RBD and VoC RBD. Overall GMT titers across antigens were highest for the hospitalized subjects (S1 858, RBD 1480, VoC 643), with symptomatic, non-hospitalized patients having lower overall GMT ELISA titers vs S1 (171), WA-1-RBD (365) and VoC-RBD (152) compared to hospitalized (Fig 1A).

• Download figure
• Open in new tab

Figure 1:

Plasma antibody titers overtime and stratified by disease severity A) Plasma ELISA titers for pre-2020 plasma (black n=6), asymptomatic (blue n=2), non-hospitalized green (n=16), and hospitalized red (n=6). Subjects who were vaccinated between sample collection are boxed out in both panels and were excluded from all statistical and data analysis. Geometric mean titers (GMT) are listed above each group. B) Comparison between ELISA titers against RBD-WA1 and VoC-RBD. Samples from subjects who received a vaccination prior to the second draw were excluded from the data analysis. Paired t-test was performed indicating significant differences between the values (p-value 0.0004) significance denoted by ^***.

Two vaccinated subjects are included in the plotted data, with both showing a boost in ELISA titer post vaccination and both achieving the highest ELISA titers against all both RBDs in this study. The vaccinees were excluded from statistical analyses. There was a significant 2.1 -fold decrease in VoC ELISA geometric mean titers when compared to WA-1 (Fig. 1B) across all samples analyzed.

SARS-CoV-2 specific MBC frequency remains detectable and out to 11.75 months post-infection

Thirty-five PBMC samples were collected from 23 subjects with 12 subjects providing multiple timepoints. In contrast to plasma Ab titers, SARS-CoV-2 S1 and RBD specific MBCs were present above the limit of detection in all subjects in this analysis, even asymptomatic subjects (Figure 2A). Moreover, SARS-CoV-2 S1 and RBD-specific MBCs remained detectable for as long as 11.75 months post-infection (Figure 2A). One subject at time-point 1 had a VoC-RBD MBC frequency below the limit of detection, however at time-point 2 this subject had detectable MBCs against RBD-VoC (Figure 2A). This is consistent with previous reports of continued germinal centers, leading to further affinity maturation over time in response to SARS CoV-2 infection [5,8,27]. When stratified by disease severity, SARS-CoV specific MBC-frequencies were closely grouped within asymptomatic and hospitalized subjects compared to non-hospitalized subjects whose frequencies that spanned the entire range from asymptomatic to hospitalized (Figure 2A). Again, two vaccinated subjects were included in the figure but excluded from analyses. Both subjects show evidence of post vaccination boost and have high overall SARS-CoV-2 specific MBC frequencies. In contrast to LLPC-derived antibodies RBD-WA-1 and RBD-VoC, MBC– derived antibodies showed little variability in their ability to recognize RBD-WA-1 vs RBD-VoC (Figure 2B) with no significant difference detected between groups.

• Download figure
• Open in new tab

Figure 2:

SARS-CoV-2 specific MBC frequency overtime and stratified by disease severity. A) MBC SARS CoV-2 specific frequency for pre-2020 plasma (black n=6), asymptomatic (blue n=2), non-hospitalized green (n=15), and hospitalized red (n=6). Subjects who were vaccinated between sample collection are boxed out in both panels and were excluded from all statistical and data analysis. Geometric mean titers (GMT) are listed above each group. B) Comparison between RBD-specific MBC frequency against WA1-RBD and VoC-RBD. T-test of log-transformed data, no statistically significant difference between WA-1 and VoC MBC Ab binding. Samples from subjects who had been vaccinated prior to blood draw were excluded from this analysis.

Relationship between WA-1 and VoC ELISA titer and MBC frequency

We next examined the relationship between RBD-WA-1 specific plasma antibodies and RBD-VoC specific plasma antibodies. We observed a strong linear correlation between WA-1 and VoC RBD specific antibodies, with a Spearman correlation coefficient R²=0.877 p-value <0.0001, with ELISA titers against RBD-VoC approximately 3-fold lower than RBD-WA-1 across all subjects (Figure 3A), WA-1 and VoC RBD specific MBC frequencies were strongly correlated as well with a Spearman correlation coefficient R²=0.885 p-value <0.0001, with some samples exhibiting lower binding to VoC-RBD, but many binding equally well to both WA-1 and VoC RBDs (Figure 3B). Finally, we noted a much weaker relationship between WA-1 or VoC RBD specific ELISA titers and WA-1 (R²=0.346, p-value <0.0001) and VoC (R²=0.331, p-value 0.0004) specific MBC frequency (Figure 3C).

• Download figure
• Open in new tab

Figure 3:

The relationship between WA-1 RBD and VoC RBD antibody binding by ELISA, WA-1 RBD and VoC RBD MBC frequency as well as the relationship between LLPC derived antibody titers and MBC frequency for each antigen. A) Correlation between RBD-VoC specific ELISA titer and RBD-WA-1 specific ELISA titer. Dotted line indicates line of identity. Spearman correlation coefficient R²=0.877 p-value <0.0001. B) Correlation between RBD-VoC specific MBC frequency and RBD-WA-1 specific MBC frequency. Samples from subjects who received vaccination were excluded from this analysis. Dotted line indicates line of identity. Spearman correlation coefficient R²=0.885 p-value <0.0001. C) Correlation between WA-1 RBD MBC frequency and ELISA titer (R²=.346) and between VoC-RBD MBC frequency and ELISA titer (R²=.331).

The relationship between clinical/disease severity and MBC frequency

Subjects reported a variety of symptoms during acute illness (Supplemental Table 1). We found no relationship between reported symptoms and hospitalization status except for diarrhea (Fisher’s exact test, p-value <0.01). In order to test whether there was a relationship between clinical data and MBC frequencies, we conducted univariate and multivariate analyses of the relationship between age, ordinal disease severity score, and MBC frequencies. Univariate analyses found no relationship between MBC frequency, ordinal disease severity score or age (Supplemental Table 2). To complete this evaluation, we then conducted two multivariate analyses testing whether there was a relationship between age, disease score and MBC frequency. When including an interaction term between age and disease score, the model found no significant relationship between age, disease score and MBC frequency (P=0.310, Supplemental Table 3). Because the distribution of ordinal disease severity scores was largely determined by hospitalization status, we revised the model to test the relationship between age, hospitalization status, and MBC frequency, again including an interaction term for age and hospitalization status and found no significant relationship between age, hospitalization status, and MBC frequency (P=0.747).