Neutralizing antibodies against SARS-CoV-2 variants of concern elicited by the comirnaty COVID-19 vaccine in nursing home residents.

Immunosenescence may impact the functionality and breadth of vaccine-elicited humoral immune responses. The ability of sera to neutralize the SARS-CoV-2 spike protein (S) from Beta, Gamma, Delta, and Epsilon variants of concern (VOCs) relative to the ancestral Wuhan-Hu-1 strain was compared in Comirnaty COVID-19-vaccinated elderly nursing home residents, either SARS-CoV-2 naïve (n = 22) or experienced (n = 8), or SARS-CoV-2 naïve younger individuals (n = 18) and non-vaccinated individuals who recovered from severe COVID-19 (n = 19). In all groups, except that including SARS-CoV-2-experienced nursing home residents, some participants lacked NtAb against one or more VOCs, mainly the Beta variant (15-20%). Serum NtAb titers were lowest against the Beta variant followed by Gamma, Delta and Epsilon variants. Overall, fold change reduction in NtAb titers relative to the ancestral strain was greatest for the Beta variant (6.7-19.4) followed by Gamma (4.8-16.0), Epsilon (2.9-13.4), and Delta (3.5-6.5) variants, although subtle differences were observed for Beta, Epsilon and Delta variants across comparison groups. In summary, older age, frailty, and concurrence of co-morbidities had no major impact on the serum NtAb activity profile against SARS-CoV-2 VOCs.

Introduction

Coronaviruses encode a proof-reading mechanism to minimize sequence variation in the long RNA genomes that may compromise their survival in nature[1]. Nevertheless, generation of variants during coronavirus replication in humans unfailingly ensues[2]. Many SARS-CoV-2 variants have emerged since its appearance late in 2019, some of which may pose a threat to pandemic control due to their relative resistance to neutralizing activity of antibodies elicited following natural infection or vaccination with the original SARS-CoV-2 isolate (Wuhan/Hu-1); these variants, termed variants of concern (VOC), incorporate amino acid substitutions or deletions in critical locations within the virus Spike (S) protein[3]. The Comirnaty mRNA COVID-19 vaccine encodes the ancestral strain S protein in a pre-fusion conformation[4]. Several studies have shown a decreased ability of serum antibodies elicited by this vaccine to neutralize a number of VOCs compared to the original strain, most notably Beta (B.1.351), Gamma (P.1), and more recently Delta (B.1.617.2)[5-12]. Nevertheless, elderly people with frailty and co-morbidities, who are at the highest risk to develop severe forms of COVID-19, were either not included or underrepresented in these studies. Since immunosenescence may be detrimental to the functionality and breadth of vaccine-elicited humoral immune responses[13], profiling of cross-neutralizing activity against VOCs of sera from fully vaccinated elderly people is of major public health interest. In the current work, we address this issue by comparing the neutralizing capacity of sera from Comirnaty COVID-19-vaccinated nursing home residents with that of sera from younger and seemingly healthy vaccinated individuals and non-vaccinated individuals who recovered from severe COVID-19.

Material and methods

Study groups

In this observational retrospective study, three population groups were considered. Group 1 included 30 nursing home (NH) residents with frailty and one or more co-morbidities (23 females), randomly selected from one NH affiliated to the Clínico-Malvarrosa Health Department, Valencia (Spain), aged a median of 90 years (range 64–100). Twenty-two were SARS-CoV-2 naïve, as judged by the absence of an historical records of having tested positive for SARS-CoV-2 RNA) or serological evidence of prior SARS-CoV-2 infection was available from any of these participants (all NHR tested negative for anti-N SARS-CoV-2 at the time of sampling), and 8 had contracted COVID-19 requiring hospitalization prior to receiving the first dose of the Comirnaty COVID-19 vaccine, as determined by previous records of positive RT-PCR results in nasopharyngeal specimens and the detection of SARS-CoV-2 N antibodies (at the time of sample collection). In all cases, the variant Wuhan/D614G was involved, as determined by whole-genome sequencing (not shown). All participants in this group were vaccinated on the same day and the sera were drawn 35 days after the second vaccine dose.

Group 2 comprised 18 healthy individuals (14 females) aged a median of 47.5 years (range 26–62), whose sera were drawn at a median of 15 days since the completion of the vaccination schedule (range 13–18). No molecular (historical records of having tested positive for SARS-CoV-2 RNA) or serological evidence of prior SARS-CoV-2 infection was available from any of these participants (all participants tested negative for anti-N SARS-CoV-2 at the time of sampling).

Group 3 included 19 unvaccinated individuals (7 females; median age, 63 years; range 40–80) who were infected with the Wuhan/D614G strain (determined by whole-genome sequencing; data not shown), required hospitalization either in clinical wards or intensive care units, and eventually recovered from severe COVID-19. Sera from these individuals were obtained at a median of 105 days (range 12–147) after molecular diagnosis of SARS-CoV-2 infection.

The study was approved by INCLIVA Research Ethics Committee, ethical statement number (2021/165). Given the retrospective nature of the analysis, the INCLIVA Research Ethics Committee exempted us from obtaining informed consent from patients. All methods were performed in accordance with Declaration of Helsinki and the Belmont Report guidelines and regulations.

Generation of SARS-CoV-2 with variant spikes

We introduced mutations in a mammalian expression vector encoding a codon-optimized SARS-CoV-2 S sequence from the Wuhan reference strain[14]. First, the D614G mutation was introduced by site-directed mutagenesis. Subsequently, additional site-directed mutagenesis and/or the cloning of synthetic fragments harboring the mutations (Gblocks, IDT) was performed to introduce all mutations using the NEBuilder HiFi DNA assembly mix (NEB). For the Beta variant (PANGO Lineage: B.1.351; mutations D80A, L241Del, L242Del, A243Del, E484K, N501Y, D614G, and A701V), a synthetic fragment encoding mutations D80A, L241Del, L242Del, A243Del, E484K, and N501Y was cloned, and the A701V mutation introduced by site-directed mutagenesis. For the Gamma variant (PANGO Lineage: P1 or B.1.1.28.1; mutations L18F, T20N, P26S, D138Y, R190S, K417T, E484K, N501Y, D614G, H655Y, T1027I, and V1176F), a synthetic fragment encoding mutations (L18F, T20N, P26S, D138Y, R190S, K417T, E484K, and N501Y) was cloned, and additional mutations H655Y, T1027I, and V1176F were inserted by site-directed mutagenesis. For the Delta variant (PANGO lineage: B.1.617.2; mutations T19R, del157/158, L452R, T478K, D614G, P681R, and D950N) mutations were introduced by site-directed mutagenesis. Finally, for the Epsilon variant (PANGO lineages B.1.427/B.1.429; mutations S13I, W152C, L452R, D614G), a synthetic fragment encoding mutations S13I, W152C, and L452R was cloned. All plasmids were verified by Sanger sequencing.

Virus neutralization assay

The neutralization capacity of circulating antibodies (NtAb) against the SARS-CoV-2 S protein was assessed using a GFP-expressing vesicular stomatitis virus pseudotyped with different Spike variants as previously described[15], but using A549-ACE2-TMPRSS2 cells (InvivoGen catalog code a549-hace2tps). All tests were done in duplicate using fivefold serum dilutions ranging from 1:20 to 1:62,500, with ~ 1,000 focus forming units per well. Following 16 h of infection, the GFP signal in each well was quantified using a live-cell microscope (Incucyte S3, Sartorius). Background fluorescence from uninfected wells was subtracted from all infected well, and the GFP fluorescence in each antibody-treated dilution was standardized to the average fluorescence observed in mock-treated wells. Any value resulting in a relative GFP signal of < 0.001 versus was assigned a value of 0.001 to eliminate negative values. Finally, the reciprocal antibody dilution resulting in 50% virus neutralization was calculated using the drc package (version 3.0‐1) in R via a three‐parameter log‐logistic regression model (LL.3 model). For the Beta, Gamma, and Epsilon variants, all assays for the same serum were performed simultaneously with the Wuhan reference strain unless neutralization was repeated due to high variability. As the Delta variant was performed subsequently, we retested each serum was against both Delta and the Wuhan strains unless no neutralization of Wuhan was observed in the previous assays. Sera testing negative (undetectable) were arbitrarily ascribed a titer of 1/20.

Quantitation of SARS-CoV-2 receptor binding domain-reactive antibodies

Total antibodies (IgG and IgM) against the SARS-CoV-2 S protein receptor-binding domain (RBD) were measured by Roche Elecsys electrochemiluminescence sandwich immunoassays (Roche Diagnostics, Pleasanton, CA, USA). The limit of detection of the assay is 0.4 IU/ml and its quantification range is between 0.8 and 250 IU/ml. Plasma specimens were further diluted for antibody quantitation when appropriate. The assay is calibrated with the first WHO International Standard and Reference Panel for anti-SARS-CoV-2 antibody[16]. Cryopreserved plasma (− 20 °C) specimens were thawed and assayed in singlets.

Statistical methods

Geometric mean titers (GMTs) are reported throughout the study. Comparison of GMTs and medians across groups was carried out by the Student's t-test or the Mann–whitney U-test, as appropriate. Correlations between variables of interest were calculated by the Spearman Rank test. Frequency comparison across groups was performed using Fisher’s exact test. Two-sided exact P values were reported. A P value < 0.05 was considered statistically significant. The analyses were performed using SPSS version 20.0 (SPSS, Chicago, IL, USA).

Results

Detection of SARS-CoV-2-RBD and neutralizing antibodies in participants

Antibodies targeting the SARS-CoV-2 RBD were detectable in all but one (belonging to group 3) out of the 67 participants. Antibody levels across the study groups, shown in Fig. 1 (panel A), were found to differ significantly, with SARS-CoV-2-experienced nursing home residents displaying significantly higher levels as compared to that in SARS-CoV-2 naïve and participants from the other two study groups. In contrast, anti-RBD antibody levels were comparable across SARS-CoV-2 naïve, and individuals in groups 2 and 3 (P ≥ 0.5). Similarly, using a pseudotyped vesicular stomatitis virus system, NtAb against the ancestral variant were detected in all but one participants (belonging group 3, who also lacked anti-RBD antibodies). NtAb levels against the Wuhan-Hu-1 variant were significantly higher in SARS-CoV-2-experienced than in SARS-CoV-2 naïve, and participants in groups 2 and 3 participants (Fig. 1, panel B). Measurable NtAb against the different SARS-CoV-2 VOCs included in the study were found at variable frequencies depending upon the study group and the virus variant considered (Table 1). In all groups except that corresponding to SARS-CoV-2-experienced nursing home residents, in which all individuals had detectable NtAb against all variants, some participants lacked NtAb against one or more VOCs, most notably the Beta variant (15–20%). As for vaccinated individuals (groups 1 and 2), no significant differences in the detection rate of NtAb for any SARS-CoV-2 VOCs were found (P ≥ 0.5) (Table 1).

Figure 1

Box-Whisker plots depicting initial (Upper pannel) anti-SARS-CoV-2 RBD total antibody levels, as measured by the Roche Elecsys electrochemiluminescence sandwich immunoassays (Roche Diagnostics, Pleasanton, CA, USA) and neutralizing antibodies against SARS-CoV-2 Wuhan-Hu-1 variant (Lower pannel), in nursing home residents (Group 1), healthy individuals fully vaccinated with the Comirnaty vaccine (Group 2) and unvaccinated, COVID-19 recovered individuals that required hospitalization in internal medicine or intensive care units (Group 3). P value for comparison across the three groups is shown.

Table 1

Detectable neutralizing antibodies against different SARS-CoV-2 variants across study groups.

SARS-CoV-2 variant	Number of participants with detectable neutralizing antibodies/total participants (%)a
	Group 1		Group 2	Group 3
	SARS-CoV-2 naïve	SARS-CoV-2 experienced
Ancestral (Wuhan-Hu-1)	22/22 (100)	8/8 (100)	18/18 (100)	18/19 (94.7)
Beta	16/22 (72.7)	8/8 (100)	16/18 (88.9)	16/19 (84.2)
Gamma	19/22 (86.3)	8/8 (100)	17/18 (94.4)	17/19 (89.5)
Delta	20/22 (90.9)	8/8 (100)	18/18 (100)	17/19 (89.5)
Epsilon	20/22 (90.9)	8/8 (100)	17/18 (94.4)	17/19 (89.5)

aGroup 1, nursing home residents fully vaccinated with the Comirnaty vaccine; Group 2, healthy individuals fully vaccinated with the Comirnaty vaccine; Group 3, unvaccinated, COVID-19 recovered individuals that required hospitalization in internal medicine or intensive care units.

Antibody neutralizing titers against SARS-CoV-2 S variants across study groups

Overall, all VOCs were neutralized to a lesser extent than the ancestral strain (Fig. 2) by sera belonging to participants in the different study groups. Specifically, serum NtAb titers were lowest against the Beta variant followed by Gamma, Delta and Epsilon variants, although wide ranges were observed for all variants. This pattern was shared by all study groups except for SARS-CoV-2-experienced participants in group 1, who displayed a remarkable NtAb activity against the Delta variant, and group 2 individuals who overall displayed an slightly decresased neutralizing activity against Epsilon as compared to the Delta variant (Table 2 and Supplementary Fig. 1). Fold change reduction in NtAb titers relative to the ancestral variant was overall greatest for the Beta variant (6.7–19.4) followed by Gamma (4.8–16.0), Epsilon (2.9–13.4), and Delta (3.5–6.5) variants. The decrease in NtAb titers for SARS-CoV-2 Gamma variant was not significantly different across comparison groups (Fig. 3); nevertheless, there were differences across comparison groups for the remaining VOCs; specifically, overall, group 3 participants exhibited a less prominent fold-change decrease in NtAb titers against the Beta variant, SARS-CoV-2-experienced nursing home residents did so against the Delta variant, whereas the latter individuals also displayed a more marked decrease against the Epsilon variant.

Figure 2

Box-whisker plots depicting neutralizing antibody titers against the SARS-CoV-2 ancestral strain (Wuhan-Hu-1) and variants of concern Beta, Gamma, Delta and Epsilon in sera form participants. P values for pairwise comparisons are shown.

Table 2

Serum neutralizing antibody titers against SARS-CoV-2 variants across study groups.

SARS-CoV-2 variant	Parameter	Study groupsa
		Group 1		Group 2	Group 3
		SARS-CoV-2 Naïve	SARS-CoV-2 experienced
Beta	GMT (SD)	73.2 (612)	6,172 (5,078)	124 (316)	122 (999)
	Mean fold change reduction of NtAb titer (CI 95%)b	19.4 (147.8)	17.4 (27.5)	13.1 (6.9)	6.7 (84.9)
Gamma	GMT (SD)	241.7 (2,013)	14,641 (6,777)	256 (435)	225 (1607)
	Mean fold change reduction of NtAb titer (CI 95%)b	5.88 (10.4)	7.35 (16.0)	6.3 (6.6)	3.6 (4.8)
Delta	GMT (SD)	216.6 (3,274)	77,157 (79,273)	543 (1,138)	231 (1,162)
	Mean fold change reduction of NtAb titer (CI 95%)b	6.57 (14.7)	1.39 (1.62)	3.0 (0.6)	3.5 (4.1)
Epsilon	GMT (SD)	330 (4,309)	8,047 (7,189)	261 (324)	283 (1,869)
	Mean fold change reduction of NtAb titer (CI 95%)b	4.3 (8.1)	13.4 (15.4)	6.2 (6.6)	2.9 (2.2)
(Wuhan-Hu-1)	GMT (SD)	1,423 (5,173)	107,541 (20,934)	1,622 (2,259)	442 (2,083)

CI confidence interval, GMT geometric mean titer, NtAb neutralizing antibodies, SD standard deviation.

aGroup 1, nursing home residents fully vaccinated with the Comirnaty vaccine; Group 2, healthy individuals fully vaccinated with the Comirnaty vaccine; Group 3, unvaccinated, COVID-19 recovered individuals that required hospitalization in internal medicine or intensive care units.

bRelative to the NtAb titer against WA1/2020.

Figure 3

Box-whisker plots depicting fold change reduction in neutralizing antibody titers against SARS-CoV-2 variants of concern as compared to the ancestral strain as a function of the study group considered. P value for comparison across the three groups is shown.

Correlation between SARS-CoV-2 neutralizing and anti-RBD antibodies

The degree of correlation between levels of SARS-CoV-2 antibodies neutralizing SARS-CoV-2 VOCs and total anti-RBD antibodies across the different population groups was examined next. Overall, the correlation between these two parameters was rather similar for the ancestral strain and and all VOCs and moderate in all cases (Fig. 4). Of interest, this pattern was rather similar for both the vaccinated participants from nursing home residents, irrespective of their SARS-CoV-2 infection status prior to vaccination, and healthy controls; nevertheless, a slightly lower correration between anti-RBD antibodies and NtAb levels was noticed for SARS-CoV-2-experienced nursing home residents (Supplementary Table 1).

Figure 4

Overall correlation between neutralizing antibody titers and anti-SARS-CoV-2 RBD total antibody levels, as measured by the Roche Elecsys electrochemiluminescence sandwich immunoassays (Roche Diagnostics, Pleasanton, CA, USA), according to the SARS-CoV-2 S-variant. Rho and P values are shown.

Discussion

Reduced neutralization of VOCs as compared to the ancestral Wuhan-Hu-1 strain has been observed in sera from individuals that were fully vaccinated against SARS-CoV-2 or have recovered from infection by either the ancestral Wuhan or the Wuhan/D164G SARS-CoV-2 strains[5-12,17,18]. This has been noticed in both the full SARS-CoV-2 virus or pseudotyped virus systems and is most notable for the Beta variant. As for recipients of the Comirnaty COVID-19 vaccine, this observation has been mostly made in vaccine clinical trials participants and seemingly healthy health workers[5-12]. Here, using S-pseudotyped SARS-CoV-2 variants, we extend this finding to include elderly nursing home residents, both in SARS-CoV-2 naïve participants and those infected by the Wuhan/D614G variant prior to vaccination. In effect, 1 month after the second vaccine dose, sera from fully vaccinated SARS-CoV-2 naïve and experienced nursing home residents retained the ability to neutralize VOCs in most and all cases, respectively, but did so to a lesser extent than the ancestral strain, especially in the case of the Beta variant. Importantly, reduction in NtAb activity against VOCs in this population group was not significantly different from that seen in healthy fully-vaccinated controls (in sera obtained a median of two weeks after vaccination); in fact, it was significantly less pronounced for the Delta variant in SARS-CoV-2-experienced individuals. Nevertheless, unvaccinated individuals who had recovered from severe COVID-19 diagnosed within a median time of 3 months had a less marked fold-change reduction in NtAb levels against the Beta and Epsilon variants. Leaving aside participants who remained unvaccinated and naturally contracted the infection, which are not strictly comparable to vaccinees, our data strongly support the idea that older age, frailty, and concurrence of co-morbidities in fully-vaccinated individuals had no major impact on the serum NtAb activity profile against VOCs, including Delta variant that has been predominant in many countries until the emergence of the Omicron variant. The reduction in NtAb activity for all VOCs relative to that observed for the Wuhan-Hu-1 reference strain (fold-change reduction in NtAb activity) was higher for all VOCs compared to other studies, especially the Beta variant. For example, using a pseudotyped virus system, Liu et al.[5-7] reported a 2.7- and 1.4-fold reduction in the NtAb activity against Beta and Delta variants, respectively, for sera from participants in the Comirnaty COVID-19 vaccine the trial, who were seemingly comparable to our healthy vaccinated controls (13.1 and 3.0 fold reduction, respectively). In turn, using live virus isolates in microneutralization assays, Lustig and colleagues[11] found fold change reduction in NtAb titers of 10.4, 2.3, and 2.1 against the Beta, Gamma, and Delta variants compared with the ancestral virus in health care workers fully vaccinated with the Comirnaty COVID-19 vaccine. Comparison of fold change reduction figures across studies is not straightforward, due not only to differences between population group characteristics, especially in terms of age, sex, comorbidities, and SARS-CoV-2 infection status prior to vaccination, but also to differences in the pseudotyped virus platforms, the chosen S sequences used to represent each variant or the use of live viruses. Despite these differences, a similar pattern of NtAb activity against VOCs can be derived from these studies as well as our current work, with the Beta variant being less efficiently neutralized than the ancestral strain, followed by the Gamma, Delta and Epsilon variants.

Levels of total antibodies binding SARS-CoV-2 RBD as quantitated by the Roche Elecsys electrochemiluminescence sandwich immunoassays strongly correlate with NtAb titers against the ancestral SARS-CoV-2 strain[19,20]. Here we extended this observation to those against Beta, Gamma, Delta, and Epsilon variants, suggesting that SARS-CoV-2 RBD antibody levels measured by this assay may be a reliable proxy for serum NtAb activity against these VOCs, regardless of the vaccinated population group considered.

The relatively small sample size could be considered a limitation of the current study; this particularly concerned the subgroup of SARS-CoV-2-experienced nursing home residents, so that conclusions drawn regarding this population group should be taken with caution and warrants further investigations. Also, a limitation inherent to the use of S-pseudotyped viruses is the inability to assess the impact of variations in other viral genes outside of the S protein that could affect sera neutralization.

In summary, herein we show that neutralizing activity of sera against SARS-CoV-2 variants carrying critical escape mutations in the S gene is decreased relative to the ancestral strain in both qualitative and quantitative terms in nursing home residents recently vaccinated with the Comirnaty COVID-19 vaccine, most notably in those who were SARS-CoV-2 naïve prior to vaccination. Nevertheless, the fold reduction in NtAb activity in this population group was not significantly different from that seen in vaccinated younger controls. Nevertheless, the degree to which reduced neutralization activity of the different VOCs, and especially the Delta variant, impact infection and/or disease in vaccinated nursing home residents translate into diminished effective protection is currently under investigation.

Supplementary Information 1.

Supplementary Information 2.