Patient Recruitment
Individuals confirmed to have COVID-19 by SARS-CoV-2 real-time reverse transcriptase-polymerase chain reaction (RT-PCR), admitted to the National Centre for Infectious Diseases (NCID), Singapore, and fully vaccinated were eligible for inclusion in the study. Fully vaccinated was defined as infection onset ≥14 days after completion of any primary COVID-19 vaccine series with a WHO Emergency Use Listing (EUL). Illness onset was defined as whichever was early from date of symptom onset or first PCR with SARS-CoV-2 detected.
Two age-matched (+/-10 years) cohorts were selected: Omicron (B.1.1529) or Delta (B.1.617.2) variant infections as assigned by direct whole-genome sequencing (WGS) of RNA extracted from nasopharyngeal swabs by the National Public Health Laboratory (NPHL) at NCID. Pangolin COVID-19 Lineage Assigner and CoVsurver were used to assign sequence lineages. Individuals were included if they were admitted to NCID between 1-18 Dec 2021 for Omicron and 27 April to 11 Aug 2021 for Delta. During these two periods, all individuals with Omicron or Delta infection respectively were admitted to hospital for isolation and evaluation regardless of disease severity, and NPHL attempted WGS for all infections with RT-PCR Cycle threshold (Ct) value <30.
Clinical Chart Review and Sample Collection
Clinical information was extracted from the medical record using a standardized data collection form adapted from the International Severe Acute Respiratory and Emerging Infection Consortium (ISARIC) case record form as previously described.14 Laboratory data including baseline investigations, Ct values from RT-PCR and serologic tests were recorded. Serial blood and respiratory samples were collected during hospitalization as part of clinical care. Serological results from Elecsys® (Roche, Basel, Switzerland) Anti-SARS-CoV-2 chemiluminescent immunoassays [anti-nucleocapsid (anti-N) and anti-spike protein receptor binding domain (anti-S)]. SARS-CoV-2 PCR was performed with a range of commercially available assays used by the clinical laboratory. Repeat RT-PCR testing was conducted at intervals based on institutional protocol and as decided by the managing physicians. At the time of the study, patients were de-isolated based on discharge criteria of at least 10 days from illness onset with Ct values ≥ 30 on RT-PCR of nasopharyngeal swabs.
Ethics Statement and Data Availability
Waiver of informed consent for collection of clinical data from individuals infected with the Omicron variant was granted by the Ministry of Health (MOH), Singapore, under the Infectious Diseases Act as part of the COVID-19 outbreak investigation. Retrospective data collection from individuals with Delta infection was approved by the institutional ethics committee (REF: 2020/01122). Written informed consent was obtained from study participants for collection of biological samples after review (REF DSRB: 2012/00917). All data sharing requests should be addressed to the corresponding authors.
Virus Isolation
Four hundred μl of the processed and filtered supernatant of nasopharyngeal swab obtained from SARS-CoV-2-positive individuals was inoculated into Vero E6 TMPRSS2 cells15 in a biosafety level 3 laboratory. Cells were cultured at 37 °C under 5% CO2 for two rounds of passages (with the second passage performed at 3 days post infection of the first passage).16,17 The cytopathic effect (CPE) was observed on day 6 post infection for each passage under an inverted microscope. Positive isolation was confirmed by observation of CPE and positive RT-PCR of SARS-CoV-2 N2 gene (CDC N2). At the end of second passage, positive CPE cultures were extracted for viral RNA to perform PCR to confirm if isolated virus was from Omicron infection using Omicron specific S-gene primers (F: 5’-TACAAACTTGTGCCCTTTTGA-3’; R: 5’-AAAAGTGAAAAATGGTGCGAG-3’).
Multiplex microbead-based immunoassay
Plasma samples were inactivated with Triton™ X-100 (Thermo Fisher Scientific) to a final concentration of 1% for two hours in the dark.18,19 Measurement of immune mediators was done using the Cytokine/Chemokine/Growth Factor 45-plex Human ProcartaPlex™ (Thermo Fisher Scientific) with the Luminex™ assay (Supplementary Appendix 1). Immune mediator levels were measured in 17 vaccinated healthy donor plasma as baseline controls. The value of Limit of Quantification (LOQ) was assigned to samples with concentrations out of the measurement range.
Whole blood immunophenotyping
Twenty five µL of whole blood was stained with antibodies as stated in Supplementary Table 1 for 20 min in the dark at room temperature (RT) and followed by lysis of erythrocytes with 1 mL of 1.2X BD FACS Lysing solution (BD cat#349202) for 10 min at RT. 300 μL of PBS was added after the incubation and centrifugated at 800 x g for 5 min. Samples were then washed with 1 mL of PBS and transferred to polystyrene FACS tubes containing 5 mL (2.4 x 103 beads) of CountBright Absolute Counting Beads (Invitrogen Cat#C36950, Lot 2324193). Cells were acquired with the CytekTM Aurora cytometer running SpectroFlo® Version 2.2.0.3 with automated unmixing and analyzed using FlowJo v10.8.0.
SARS-CoV-2-specific T cells by Intracellular Cytokine Staining
To profile the SARS-CoV-2 specific T effector subsets in COVID-19 patients, intracellular cytokine staining was carried out as previously described20 with modification. Briefly, PBMCs collected from whole blood of vaccine breakthrough cases and healthy vaccinated controls were stimulated with phorbol 12-myristate 13-acetate (PMA 100ng/mL, Sigma Aldrich) and ionomycin (1µg/mL, Sigma Aldrich), or pooled SARS-CoV-2 PepTivator S, S1, M and N peptides (0.6nmol/mL each) (Miltenyi Biotec) for six hours. Brefeldin A and Monesin (1× , ThermoFisher Scientific) were added at two hours post-stimulation. Cells were stained for surface markers in the dark at room temperature for 30 minutes (Supplemental Table 2, Row 1-21), followed by fixation and permeabilization for 30 mins with Foxp3/Transcription Factor Staining Buffer Set (ThermoFisher Scientific). Permeabilized cells were then stained for intracellular cytokines in the dark at room temperature for 30 minutes (Supplemental Table 2, Row 22-29). Cells were then washed with PBS and centrifuged at 800 x g for 5 minutes before transferring to respective polystyrene FACS tubes containing 5 mL (4.8 x 103 beads) of CountBright Absolute Counting Beads (Invitrogen Cat#C36950, Lot 2324193). Cells were acquired with the CytekTM Aurora cytometer running SpectroFlo® Version 2.2.0.3 with automated unmixing and analyzed using FlowJo v10.8.0.
Spike protein flow cytometry-based assay (SFB assay) for antibody detection
The SFB assay was performed as previously described.21,22 S protein-expressing cells were seeded at 1.5 x 105 cells per well in 96 well V-bottom plates. Cells were incubated with human serum (diluted 1:100 in 10% FBS) followed by a secondary incubation with a double stain, comprising Alexa Fluor 647-conjugated anti-human IgG (1:500 dilution) and propidium iodide (PI; 1:2500 dilution). Cells were acquired using a BD Biosciences LSR4 laser and analyzed using FlowJo (Tree Star). The assay was performed as two independent experiments with technical duplicates each time.
Memory B cell ELISpot
SARS-CoV-2 RBD-specific memory B cell numbers were counted using ELISpot. MultiScreenHTS IP Filter Plate, 0.45 µm plates (Merck Millipore) were coated overnight at 4°C with purified anti-human-IgG (MT91/145, Mabtech) prepared at 15μg/ml in PBS 1X. Plates were washed 4 times with PBS 1X and blocked at least 30 minutes at room temperature with RPMI 10% FBS. Fresh PBMCs were plated directly into ELISpot plate (100,000 and 400,000) to determine RBD-specific plasmablast numbers. In parallel, 1 x 106 PBMCs were resuspended in 1ml RPMI + 10% FBS + 1 μg/ml R848 + 10 ng/ml IL-2, and incubated at 37°C, 5% CO2 for 4-5 days to differentiate memory B cells into antibody-secreting cells. After incubation, cells were counted, and 100,000 or 400,000 live cells were taken for ELISpot plating to determine RBD-specific memory B cell numbers. Total IgG secreting cells were determined by plating 1,500 or 3,000 live cells. In both cases (plasmablasts or memory B cells), cells were incubated 18-22h at 37°C, 5% CO2 in the ELISpot plate before detection. A combination of RBD-WASP/anti-WASP-ALP or anti-IgG-biotinylated/streptavidin-ALP (Mabtech) were used to detect RBD-specific or total IgG secreting cells respectively. Plates were then read on an IRIS ELISpot reader (Mabtech). Spots were calculated based on the average of two wells.
Statistical Analysis
For descriptive analysis, data are presented as median (interquartile range [IQR]) for continuous parameters and frequency (percentage) for categorical variables. Fisher’s exact test was used to compare categorical variables, while for continuous variables, t-test was used for normal data and Mann-Whitney U test for non-normal data. All statistical tests were 2-sided and P values <0.05 were considered statistically significant. For asymptomatic patients, day one of illness was assigned to the day of confirmatory diagnosis by RT-PCR. For symptomatic patients, day of symptom onset was denoted as day one of illness.
For serial Ct values, we fitted a generalized additive mixed model (GAMM) with a random intercept by patient. A Ct value of 45 was imputed where the PCR result was not detected. We compared this with previously reported serial Ct values of Delta infected patients.23 We plotted Ct values with marginal effect of day of illness by variant type and 95% confidence intervals (CI) from the GAMM.
For immunological analysis, to compare between multiple groups, Kruskal-Wallis tests with post hoc tests using Dunn’s multiple comparison tests were used to identify significant differences. For whole blood immunophenotyping, gated cells were manually exported using FlowJo v10.8.0. Samples were then used for Uniform Manifold Approximation and Projection for Dimension Reduction (UMAP) analysis using cytofkit2 R Packages with RStudio v3.5.2.24 Thirteen healthy vaccinated controls, 15 Delta vaccine breakthrough patients and 17 Omicron breakthrough patients were each concatenated to their respective groups and 150,000 cells were analyzed using the ceil method. Principal Component Analysis (PCA) was performed on the systemic cytokine data using Singular Value Decomposition (SVD) method in ClustVis.25