At Least Three Doses of Leading Vaccines Essential for Neutralisation of SARS-CoV-2 Omicron Variant

SARS-CoV-2 Stock Generation and Characterisation

Three SARS-CoV-2 isolates, VIC31 (hCoV-19/Australia/VIC31/2020 containing D614G mutation), Delta variant of concern (hCoV-19/Australia/VIC18440/2021), and Omicron BA.1.1 variant of concern (hCoV-19/Australia/VIC28585/2021) were provided by the Victorian Infectious Diseases Reference Laboratory (VIDRL; Melbourne, Australia). Virus stocks were propagated and titrated in Vero E6 cells (American Type Culture Collection (ATCC), Manassas, VA, USA) prior to use as previously outlined in [5]. Briefly, VeroE6 cells were grown in 150 cm² flasks in Dulbecco’s Modified Eagle Medium (DMEM) containing 10% heat-inactivated fetal bovine serum (FBS), 2mM GlutaMAX supplement, 100U/mL penicillin, and 100 μg/mL streptomycin (all components from ThermoFisher Scientific; Scoresby, VIC, Australia) until >80% confluent. Virus isolates were diluted 1:100 in DMEM (containing 2% FBS, 2mM GlutaMAX supplement, 100U/mL penicillin, and 100 μg/mL streptomycin), and 5 mL was used to inoculate Vero E6 cells for 1 hr at 37°C/5% CO₂ before additional media was added to the flask. The flasks were incubated at 37°C/5% CO₂ for 48 h (for VIC31 and Delta) or 72 hr (for Omicron) before supernatant was clarified at 2,000 x g for 10 min, and harvested and stored in 1 mL aliquots at -80 °C.

Identity of virus stocks were confirmed by next-generation sequencing using a MiniSeq platform (Illumina, Inc; San Diego, CA, USA). In brief, 100 µL cell culture supernatant from infected Vero E6 cells was combined with 300 µL TRIzol reagent (Thermo Fisher Scientific) and RNA was purified using a Direct-zol RNA Miniprep kit (Zymo Research; Irvine, CA, USA). Purified RNA was further concentrated using an RNA Clean-and-Concentrator kit (Zymo Research), followed by quantification on a DeNovix DS-11 FX Fluorometer. RNA was converted to double-stranded cDNA, ligated then isothermally amplified using a QIAseq FX single cell RNA library kit (Qiagen, Hilden, Germany). Fragmentation and dual-index library preparation was conducted with an Illumina DNA Prep, Tagmentation Library Preparation kit. Average library size was determined using a Bioanalyser (Agilent Technologies; San Diego, CA, USA) and quantified with a Qubit 3.0 Fluorometer (Invitrogen; Carlsbad, CA, USA). Denatured libraries were sequenced on an Illumina MiniSeq using a 300-cycle Mid-Output Reagent kit as per the manufacturer’s protocol. Paired-end Fastq reads were trimmed for quality and mapped to the published sequence for the SARS-CoV-2 reference isolate Wuhan-Hu-1 (RefSeq: NC_045512.2) using CLC Genomics Workbench version 21 from which consensus sequences were generated. Stocks were confirmed to be free from contamination by adventitious agents by analysis of reads that did not map to SARS-CoV-2 or cell-derived sequences.

Human Plasma Samples

Blood samples in EDTA blood collection tubes (BD Biosciences, Australia) were collected from healthy volunteers, aged between 25-70 (both female and male), who received COVID-19 vaccine manufactured by either AstraZeneca (Vaxzevria/AZD1222; University of Oxford-AstraZeneca, Oxford/Cambridge, UK), Pfizer (Comirnaty/BNT162b2; BioNTech-Pfizer, Brooklyn, New York, USA), or Moderna (Spikevax/mRNA-1273; Moderna Inc, Cambridge, MA, USA) (Table 1). AstraZeneca vaccine was advised for the age group >60 [6]. Plasma was separated and stored at -70°C until use. The study protocols were approved by institutional ethics committees of the Peter Doherty Institute, Melbourne Australia [University of Melbourne Central Human Research Ethics Committee (2021-21198-15398-3)] and the CSIRO (CSIRO Human Research Ethics Committee ID 2021_123_RR). Biosafety protocols for handling the human samples and the infectious agents were approved by the Institutional Biosafety Committee of ACDP, Geelong before infectious work commenced. Blood samples were collected on the day of vaccination (baseline; pre-1^st dose) and 2 weeks post second dose (2wk-2^nd dose) for all the vaccine groups and after six months of 2^nd dose (6mo-2^nd dose) and two weeks post third dose (2wk-3^rd dose) for the Pfizer group.

Live-Virus Neutralisation Assays

Virus neutralization assays (VNT) were carried out using Vero E6 cells as described previously [5]. Briefly, each plasma sample for the three vaccines was diluted 1:10 in Dulbecco’s Phosphate Buffered Saline (DPBS; Thermo Fisher Scientific, Waltham, MA, USA) in a deep-well plate, followed by a two-fold serial dilution up to 1:1,280. The dilution series for each plasma sample was dispensed into triplicate columns of a 96-well plate (one plate per isolate), for a total volume of 50 µl per well. For the plasma-containing wells, 50 µl virus diluted in DMEM (Thermo Fisher Scientific, Waltham, MA, USA) containing 2% FBS, 2mM GlutaMAX supplement, 100U/mL penicillin, and 100 μg/mL streptomycin to contain approximately 100 TCID₅₀ (checked by back-titration) was added to each well. The plates were incubated at 37°C/5% CO₂ for 1 h to allow neutralization complexes to form between the antibodies present in the plasma and the virus. A positive control serum was included to confirm the reproducibility of the assay. At the end of the incubation, 100 µl VeroE6 cells were added to each well and the plates were returned to the incubator for 4 days. Each well was scored for the presence of viral CPE, readily discernible on Day 4 post-infection. NT₅₀ neutralization titres calculated using the Spearman–Karber formula [7], and transformed to log2 values for analysis. Replicates that did not show neutralisation at 1:10 were scored <1:10 dilution and assigned a value 1:5 for statistical analysis in line with the current best practice [8].

In silico Modelling

Molecular models were made of the trimeric Omicron Spike protein from residues 13 to 1160 (omitting the transmembrane domain) consistent with variant BA.1.1 and including glycosylation to visualize structural changes from the original strain. ‘AlphaFold’ [9] was used to reconstruct the N-terminal domain due to insertions and deletions. Models were built to include the ‘up’ and ‘down’ conformations of the Receptor Binding Domain (RBD, residues 330-530) as well as the human angiotensin converting enzyme 2 (ACE2) receptor. Building upon our previous work, [10] models were solvated and simulated using ‘NAMDv2.14’ software]11] for over 200 nanoseconds before visual assessment of changes to surface epitopes and receptor binding interfaces. More modelling details included in Supplementary Methods.

Statistical Analysis

The VNT data were grouped based on vaccines, sex (female and male), and age groups (<35 – young; 35-60 – middle and >60 – senior). Both one-way and two-way Analysis of Variance (ANOVA) were used to test the statistical differences between vaccines, sex, age, day post vaccination/booster, variants, and their interactions. Levene’s test for Homogeneity of Variance was performed before two-way ANOVA. If the ANOVA results returned a p-value of <0.05, a post hoc test with a Tukey’s range test (Tukey’s Honestly Significant Difference) was performed to measure the interactions between two variables. Plots were drawn using ‘ggplot2’ and all statistical procedures used the ‘car’ and ‘lme’ libraries in R [12]. Longitudinal data for continuous outcomes in multiple vaccine groups were compared using a linear mixed effects model. Subject identification number was added as a random variable while day post vaccination/booster (as a factor) and vaccine, sex and age groups and the interactions were considered as explanatory variables. In all models, explanatory variables were selected based on the lowest Akaike Information Criterion (AIC) using forward selection.

Neutralisation of Delta and Omicron Variants of Concern

Table 1 shows the median ages for AstraZeneca sample donors are 57.5 compared to 38.5 for Moderna, and 33 for Pfizer vaccine donors respectively, the reasons for which are as follows. The Australian Technical Advisory Group on Immunisation (ATAGI) changed its recommendations regarding the AstraZeneca vaccine following reports of thrombotic thrombocytopenia syndrome (TTS) emerging from vaccinated adults [13]. It is worth noting that some individuals younger than 60 have received AstraZeneca vaccine following an assessment by a qualified health processional and with a verbal or written consent as provisioned by the Department Health, Government of Australia [6], but there is no brand preference for people aged 60 years and above because the benefits outweigh the risk of TTS in this age group. As a result, samples predominantly originated from young and middle-aged individuals for the Moderna and Pfizer vaccines, and from middle-aged and senior adults for the AstraZeneca vaccine. Given the large age range for the Pfizer vaccine, for this manuscript we classified the age groups in to young (<35 age), middle (35-60) and senior (>60). We chose 35 as the cut-off for middle age as it generally corresponds to parents of teen aged children and is used in metrics such as the Peterson-KFF health System Tracker [14] which has identified COVD-19 as the second and third most important cause of death of people who are over and under 35, respectively. In Australia, Moderna vaccine was approved for use relatively recently (9 August 2021), while a third dose of the AstraZeneca is only approved in exceptional circumstances, therefore this study did not have access to third-dose samples for these vaccines.

Calculation of fold changes using the data presented in Table 2 and Figure 1, demonstrate that for the Pfizer vaccine there was an 8.2- and 5.5-fold reduction in NT₅₀ titres against VIC31 and Delta respectively 6 months after second dose administration (compared to two weeks), while the average titre against Omicron decreased from an estimated 1:6.3 (only 4/15 samples had quantifiable neutralization) to below assay detection limit for all samples at the same time points. Subsequently there was a 25.5- and 24.7-fold increase in titres against VIC31 and Delta respectively two weeks after a third dose administration (compared to 6 month post-second dose), while the average titre against Omicron increased from below assay detection limit to 1:48.5 at the same time points. Because of the aforementioned limitations, calculation of such fold changes for AstraZeneca and Moderna were not possible. However, it is worthwhile to compare their neutralisation titres two weeks post-second dose relative to the Pfizer vaccine samples. Across all three variants, we did not find any significant difference between Pfizer and Moderna vaccines by 2^nd week post second dose of the vaccines (Supplementary Table S1).

• Download figure
• Open in new tab

Figure 1: Box Plots Showing 50% Neutralisation Titres (NT₅₀) Against SARS-CoV-2 Variants of Concern with Plasma from Human Volunteers

Neutralisation of VIC31, Delta, and Omicron by human donor samples collected at different timepoints following vaccination with AstraZeneca (A), Moderna (B), or Pfizer (C) COVID-19 vaccines. The titres are shown in Log2 scale on the y-axis.

ANOVA and Linear Mixed Models (lm and lmer) analyses were used to investigate the interactions of age, sex, vaccine, and day post vaccination/booster with neutralising antibody titres to the three different virus variants (Supplementary Text S1). One-way and two-way ANOVA analyses showed that, for all vaccines (Pfizer, Moderna and AstraZeneca), day post vaccination/booster and virus variant were the main contributors to differences in neutralising antibody responses, while age and sex did not have significant effects. Based on ANOVA analysis there was no significant difference in neutralising titres based on age or sex for any of the vaccines (Supplementary Tables S2 & S3; Figure 2A-C). However, based on Linear Mixed Models (lm) analysis, younger individuals tended to have higher neutralising antibody titres (Supplementary Tables S4 & S5; Figure 3A-C). The titres were dependent on the vaccine, day post-vaccination and variant (p<0.001) and to an extent the age where antibody titres in young showed statistically significant difference when compared to the middle and senior age groups (p=0.022). The antibody titres against Omicron have a negative slope and indicating that there is poor neutralisation with the variant.

• Download figure
• Open in new tab

Figure 2: Box Plots Showing 50% Neutralization Titres (NT₅₀) Against SARS-CoV-2 Variants of Concern Assessed by Sex of Donors

Plasma samples from human volunteers vaccinated with either AstraZeneca (A), Moderna (B), or Pfizer (C) vaccines were tested in neutralisation assays against VIC31 (i), Delta (ii), and Omicron (iii), with titres assessed with respect to sex. The titres are shown in Log2 scale in the y-axis. The ANOVA results of pairwise comparison between pre and post vaccinated sera for the different sexes are presented as adj p-values using Bonferroni post hoc test.

• Download figure
• Open in new tab

Figure 3: Box Plots Showing NT₅₀ Values Against SARS-CoV-2 Variants of Concern Assessed by Age Group

Plasma samples from human volunteers vaccinated with either AstraZeneca (A), Moderna (B), or Pfizer (C) vaccines were tested in neutralisation assays against VIC31 (i), Delta (ii), and Omicron (iii), with titres assessed with respect to age group. The titres are shown in Log2 scale in the y-axis. The ANOVA results of pairwise comparison between pre and post vaccinated sera for the differen sexes are presented as adj p-values using Bonferroni post hoc test.

In line with recent findings in similar studies by other groups [15,16], Omicron titres were significantly lower than VIC31-D614G and Delta titres at all time points and regardless of vaccine, highlighting the concern of poor protection to Omicron following vaccination by any of the three leading vaccines. One shortcoming of our data is the lack of later samples and post-boost (third vaccination) from individuals vaccinated with AstraZeneca and Moderna vaccines. However, the data available from Pfizer vaccinated individuals at 6 months after receiving their second dose indicate a significant waning of immunity to all three variants over this period. Neutralising antibody levels to Omicron were low or undetectable in most individuals two weeks after receiving two Pfizer doses, and below the assay limit of detection in all samples tested at 6 months after the second dose. However, at two weeks after the booster vaccination (third Pfizer dose), a significant booster effect in neutralising titres was noted to all three variants. This highlights the importance of the need for a booster (third) vaccine dose for ongoing protection against D614G and Delta variants, but more importantly that two doses are insufficient to stimulate significant neutralising protection against Omicron BA.1, therefore making the third dose essential. The cause of better protection against Omicron BA.1 following three doses of vaccine incorporating Spike protein lacking Omicron specific mutations, is likely due to higher overall neutralising titres due to a booster effect despite neutralising antibodies to certain epitopes being circumvented.

Linear Mixed Model analysis of the variance between variables confirms the ANOVA analyses which showed that day post vaccination/booster and virus variant are the major contributors to titre differences, but further shows that there was a significant difference in neutralising titres elicited by the three different vaccines, mostly due to poor responses to the AstraZeneca vaccine. This analysis model further suggests that younger individuals (aged less than 35 years) are likely to have a stronger neutralising antibody response than older individuals, but we concede that this finding might be confounded by the distribution of the ages within the vaccine groups. This might also explain the poorer responses in the AstraZeneca group because of the vaccine advice in Australia (AstraZeneca mostly used in the older age group).

In silico modelling and analysis

From structural perspectives, there are significantly more mutations in the Spike protein of the Omicron variant shown in Figure 4A, (37 mutations in the BA.1 variant, 32 in the BA.2 variant including insertions and deletions), compared to the Delta variant (10 mutations), while only one in the VIC31 strain (D614G). The Omicron mutations contribute to its ability to circumvent vaccine or infection induced neutralising antibody responses [15,16]. Three of these mutations (G142D, T478K and D614G) are common to Delta and Omicron, however we have previously shown experimentally and through in silico modelling that D614G has no significant impact on vaccines [4]. In the receptor binding domain (RBD) of the Spike protein alone, which is most immunologically exposed in its ‘up’ configuration and critical for ACE2 binding, there are 15 and 17 mutations in the Omicron BA.1 and BA.2 variants compared to only 2 in Delta as shown in Figures 4A & C. Many of the observed Omicron mutations by our analysis appear to be immune evasive, as they do not appear to be at the ACE2 binding interface and represent a significant change in the residue characteristics; in the RBD, these include N440K, K417N, S477N, T484K, and E484A.

• Download figure
• Open in new tab

Figure 4: Structure of SARS-CoV-2 Spike Protein and Mutations of Interest

A. The Omicron spike protein with mutations shown in pink. A single S1 domain is highlighted in blue with the receptor binding domain (RBD) in the ‘up’ configuration is highlighted in cyan. Glycosylation is shown as green lines. B. & C. Comparison of the RBD in the Omicron and Delta variants, 15 mutations are observed in Omicron while only 2 in the Delta variant.

Interestingly our isolate VIC28585 belongs to the BA.1.1 sub-lineage as it also contains the R346K spike mutation, which possibly further decreases antibody neutralization based on the reduced binding of class 2 antibodies in theoretical studies with the Mu variant of interest [17]. Next-generation sequencing results also indicated the presence of additional Spike mutations S373P and S375F in our VIC28585 isolate at frequency exceeding 99.8%; these two mutations with Grantham scores of 74 (moderately conservative change) and 155 (radical change) respectively are near the S371L Spike mutation characteristic to BA.1 which has a Grantham score of 145 (moderately radical change). Thus, we have investigated the neutralising responses to Omicron BA.1.1 in individuals vaccinated with three leading COVID-19 vaccines and compared these neutralising titres to those against an early SARS-CoV-2 isolate with the D614G mutation (VIC31), and the Delta variant (B.1.617.2) that subsequently became the dominant variant worldwide prior to Omicron’s emergence. For example, on 11 December 2021 the ratio of Delta : BA.1 : BA.1.1 : BA.2 : Other in USA (where extensive data is available from the CDC; [18]) was 92.5 : 5.2 : 2.2 : 0 : 0.2, and this steadily changed over two months so that on 12 February 2021 the ratio was 0 : 22.9 : 73.2 : 3.9 : 0. This shows that BA.1.1, used in our study, could emerge as the dominant variant although we need to wait for more comprehensive data from around the world – as of 16 February 2022, 719,364 sequences of BA.1 were reported from 139 countries, while 410,015 sequences of BA.1.1 were reported from 127 countries.

Other Omicron mutations are likely to confer increased binding affinity to ACE2 by increasing molecular interactions such as salt bridges and pi-stacking [9]; these mutations include Q493R and Q498R (increasing salt bridge interactions with ACE2’s D38 and E35) and N501Y (also seen in Alpha, Beta and Gamma variants, pi-stacking with ACE2’s Y41, shown in Supplementary Figure S1. A cluster of mutations S371L, S373P and S375F, interestingly appear co-located with Y505H in adjacent RBD domains in the ‘down’ position. These serine mutations, as well as significantly altering epitope presentation of the RBD domain by creating a more rigid structural motif, may allow H505 to interact with the adjacent F375, potentially modulating ‘up’/’down’ conformations via the ionizable H505, (pKa of 6.0). This structural change may coincide with the observed Omicron preference of infection to the bronchus, where pH is usually lower at 5.5 to 6.5, compared to the lung [19,20] (Supplementary Figure S2).

With the large number of mutations, it is unsurprising that the Omicron variant displays the observed antibody evasion and replicative advantage over earlier variants [21], as experimentally observed in our study. Our biomolecular modelling shows sufficient distinction between BA.1 and BA.2, therefore future studies should experimentally compare the two.

Public Health Implications and Future Studies

The data we present here highlight the importance of booster vaccination for adequate protection against newly emerging variants such as Omicron BA.1. Further studies are needed to perform similar analyses of neutralising titres in individuals receiving “ heterologous” booster doses, e.g., initial AstraZeneca vaccination followed by Moderna or Pfizer boost, or Pfizer-Moderna / Moderna-Pfizer.

In the early days of the pandemic, the scientific prediction was that SARS-CoV-2 will have a mutation rate similar to SARS-CoV (0.8-2.4×10⁻³) thanks to the 3′-to-5′ exoribonuclease ‘proofreading mechanism’ in coronaviruses [22], which was thought to result in a reduction in the mutation rate (∼10⁻⁶ per site per cycle) compared to influenza (∼3×10⁻⁵ per site per cycle) [23,24]. Therefore, it was predicted that vaccine matching may not be required as frequently as seasonal influenza (once a year each for Northern and Southern hemispheres), although this was hard to say definitively because SARS-CoV-2’s genetic drift was uncertain with respect to the long-lasting efficacy of vaccine candidates then under development [25]. Early studies, for instance between 1 February to 1 May 2020 seemed to confirm the mutation rate prediction with a median estimate of 1.12×10⁻³ per site-year ([26]; we believe their upper 95% confidence interval should read 1.85×10⁻³).

Two years on, we have increasing evidence, including from this work, that the primary immunization schedule should consist of three rather than two doses, and indeed many countries are redefining ‘fully vaccinated’ status as three doses of these vaccines [27]. It is early to speculate on annual boosters beyond the third dose, although Israel has started to expand its fourth dose administration to groups beyond the immunocompromised [28,29]. These developments make it all the more important to ensure vaccine equity across the world so everyone is adequately protected, and highlight the potential need for vaccine matching to increase protection against infection. Vaccines using Spike as the primary antigen should consider including mutations common to BA.1, BA.1.1, and BA.2; development of vaccines with additional antigens, other than Spike, targeting more conserved viral factors should also be explored. This view is further supported by a recent preprint currently under peer review suggesting that, although a fourth dose restored anti-Omicron BA.1 neutralising antibody titres to a level observed post-third dose, breakthrough Omicron infections were common with high virus titres, albeit with very mild clinical signs [29].

In this context, it is important to explore vaccines that can better withstand mutations to the Spike protein better, and ideally have less dependence on cold chain. In addition to repeating this work with Omicron BA.2, it is also important to explore the benefits of T cell immunity and heterologous boosting (which are beyond the scope of this paper, but see for instance these studies [30-33]. Another aspect worth exploring is intranasal administration of selected vaccines. In the earliest study (conducted from March 2020) looking at intranasal administration of the AstraZeneca vaccine [34], we found that although neutralising antibody titres were 2-fold lower in ferrets receiving two doses of this vaccine intranasally compared to the intramuscular route, the former was more protected using virus shedding as an infection metric [35]. These findings were subsequently confirmed in hamsters and non-human primates [36], and also reportedly led to a small clinical trial in 30 humans [37]. We feel that the intranasal route, and the role of mucosal immunity, needs to be explored more comprehensively given that the AstraZeneca vaccine has the widest reach (183 countries worldwide), with over 2.5 billion doses administered [38].