Characterizing Longitudinal Antibody Responses in Recovered Individuals Following COVID-19 Infection and Single-Dose Vaccination in British Columbia, Canada: a Prospective Cohort Study

2.1. Study Design

A prospective observational cohort termed CARE (Characterizing the Antibody Response to Emerging COVID-19) was established from individuals who recovered from SARS-CoV-2 infection, for the purposes of investigating antibody responses against several SARS-CoV-2 epitopes (full spike (S), receptor binding domain (RBD) and nucleocapsid (N)), as well as against the S protein of endemic HCoVs (OC43, HKU1, NL63, 229E), at least 2 weeks post natural SARS-CoV-2 infection, with or without subsequent SARS-CoV-2 vaccination. Information on age, biological sex and date of PCR-diagnosis was collected through medical records, while information on the duration of COVID-19 symptoms, hospitalization, and vaccination was collected through an online self-reporting survey. The date of confirmed PCR diagnosis of SARS-CoV-2 infection was used to estimate ‘days post-infection’. Survey data and participant informed consent were collected and managed using REDCap electronic data capture tools hosted at BC Children’s Hospital (Vancouver, BC). REDCap (Research Electronic Data Capture) is a secure, web-based application designed to support data capture for research studies (23). Participants were enrolled in the cohort from November 19^th, 2020, to September 7^th, 2021. During this time the most prevalent SARS-CoV-2 variant in British Columbia transitioned between the Alpha, Beta, Gamma and Delta genotypes. The Beta variant was detected in late 2020 and January 2021, increase prevalence of the Alpha variant shortly followed and it remained dominant until June 2020. The Gamma variant was first detected in late February 2021, it’s incidence surpassed Alpha in July 2021. Public health surveillance first recorded the Delta variant in March, and it was responsible for most sequenced cases over the summer from July to September 2021(24). All data analysis was performed in R version 4.0.4 using the packages: ‘DataExplorer’, ‘survival’, ‘survminer’, ‘dplyr’, ‘ggfortify’, ‘tableone’, ‘naniar’, ‘RColorBrewer’, ‘lme4’, ‘mgcv’, ‘gam.check’ and ‘readr’ (25).

2.2. Recruitment Criteria

Adults 18 years of age and older from the greater Vancouver metropolitan area were recruited if they had a confirmed PCR-positive SARS-CoV-2 infection and if they were no longer required to self-isolate per the BC provincial public health guidelines (i.e., tested positive for SARS-CoV-2 at least 14 days prior). The study protocol was approved by the University of British Columbia (UBC) Clinical Research Ethics Board (H20-01089).

2.3. Sample Collection and Processing

Participants were asked to donate blood samples collected by venipuncture (for serological testing) and concurrent self-collected saline gargle samples (for SARS-CoV-2 real-time-polymerase chain reaction (PCR) testing), every two weeks for 3 months post-recruitment (up to 7 collections total). Blood was drawn in gold-top serum separator tube with polymer gel (BD, cat# 367989); after at least 30 minutes of clotting at room temperature, the blood sample was then centrifuged at 1,400 G by staff at the collection site and sent to the British Columbia Centre for Disease Control Public Health Laboratory (BCCDC PHL). At the BCCDC PHL the samples were divided into serum aliquots that were frozen at −80°C within four hours of receipt. Blood collections occurred at four sites in the Greater Vancouver Area, British Columbia Canada: BC Children’s Hospital, St. Paul’s Hospital, Abbotsford Regional Hospital and Surrey Memorial Hospital. Saline gargle samples were self-collected by the participants at home, in accordance with well-validated instructions (26), on the day of blood collection and transported to the BCCDC PHL by the blood collection site. Self-collected saline gargle samples were tested for SARS-CoV-2 by PCR. SARS-CoV-2 whole genome sequencing was done on a subset of participants’ SARS-CoV-2 positive diagnostic clinical specimens (27). All molecular, genomic, and serological testing (described below) for participant specimens was conducted centrally at the BCCDC PHL.

2.4. Measurements of Humoral Immunity

All serum samples were initially tested using a combination of three Health Canada approved chemiluminescent immunoassays: 1) total antibodies to SARS-CoV-2 RBD (Siemens SARS-CoV-2 Total Assay [COV2T]), 2) total antibodies to SARS-CoV-2 S (Ortho VITROS™ Anti-SARS-CoV-2 Total) and IgG anti-N antibodies (Abbott ARCHITECT™ SARS-CoV-2 IgG), as per manufacturer guidelines, with results interpreted as reactive or non-reactive using the manufacturer-recommended signal to cut-off ratios (28, 29). All available samples were then tested using the V-PLEX COVID-19 Coronavirus Panel 2 (IgG) (Mesoscale Diagnostics LLC (MSD): #K15369U), the diagnostic accuracy of the MSD assay was previously validated through comparison with alternative Health Canada approved tests at the BCCDC PHL (30). The MSD assay provides quantitative measures of IgG antibodies against RBD, S and N SARS-CoV-2 epitopes, as well as IgG antibodies against S of the four seasonal HCoVs. Serological specimens were processed as previously reported (30). Quantitative antibody levels expressed as log10 antibody units (AU)/ml were recorded and evaluated for all tested samples. MSD results were interpreted as reactive or non-reactive using the MSD recommended signal thresholds for serum: SARS-CoV-2 anti-S IgG = 1960; anti-N IgG = 5000; anti-RBD IgG = 538. Cutoffs (derived at the BCCDC PHL) for seasonal HCoVs seropositive status are as follows: HCoV-229E anti-S IgG = 1700; HCoV-HKU1 anti-S IgG = 900; HCoV-NL63 anti-S IgG = 270; HCoV-OC43 anti-S IgG = 2000 (31). Samples were stratified by collection time to < 6 months and ≥ 6 months post infection and percent positivity was compared using a Chi-square test (χ²test).

2.5. Power Analysis for Investigating Association between IgG Concentration and Vaccination

A power calculation was conducted to determine the minimum number of paired participant samples needed to estimate at least a 70% association between COVID-19 vaccination and HCoV anti-IgG antibody concentration. Antibody concentrations were assumed to be normally distributed with a standard deviation of one. A significance level of 5% and two-sided alternative were used (32).

2.6. Analytic Data Selection

To analyze antibody dynamics, an analytic dataset was selected from the CARE COVID-19 cohort. At least k = 18 paired participant samples are required to estimate a 70% or greater association between COVID-19 vaccination and anti-HCoV IgG antibody concentration. Exclusion criteria were applied to select an analytic dataset from k = 57 participants with n = 341 observations. One participant had no follow up samples and was omitted from the analytic dataset (k = 1, n = 1). Six participants were excluded because they were vaccinated before collection of their baseline sample (k = 6, n = 37). Eight participants were removed from the analytic dataset because of missing data in their survey results (k = 8, n = 33). After applying the exclusion criteria, the analytic dataset contained k = 42 participants with n = 270 observations (Figure S1). There were k= 41 participants with >1 pre-vaccine sample (n=210 pre-vaccine observations) used for analysis of antibody waning pre-vaccination.

2.7. Bivariate Data Analysis

2.8. Primary Analysis

Primary analysis used a multivariable linear mixed-effects model to regress SARS-CoV-2 anti-IgG concentration on vaccine status adjusting for dependency within participant samples and covariates defined as potential confounders by the common cause criterion (33, 34). Separate models were fit for anti-S and anti-N IgG signals. Unconditional mean models were used to find the intraclass correlation coefficient (ICC) before covariates were added to build fixed effect models (35). Effect modification terms were assessed by the Akaike information criterion and included in the fixed effect models to understand if time from vaccination influences SARS-CoV-2 antibody concentrations (36).

2.9 Secondary Analysis

Secondary analysis employed a Kaplan-Meier curve to estimate the cumulative incidence of seroreactivity stratified by vaccine status. The survival function was transformed to cumulative incidence by 1-S(τ) (37). Seroreactivity was defined from the distribution of SARS-CoV-2 anti-S IgG concentration at first blood draw (baseline); the 95^th percentile was chosen as the threshold (5.5 log10 AU/ml). Participants were censored if they were not seroreactive before loss to follow-up (right censoring). A log-rank test was used to test the hypothesis that the cumulative incidence of seroreactivity between unvaccinated and vaccinated persons, who have been previously naturally infected with SARS-CoV-2, does not differ (38).

3.1. CARE COVID-19 Cohort

Fifty-seven individuals with COVID-19 were recruited through the CARE COVID-19 Cohort. Subjects (17 male, 40 female; 18 to 76 years old) represented a range of COVID-19 disease severity cases. Most subjects had a mild case of COVID-19, defined as not requiring hospitalization; 6 reported being asymptomatic and 12 reported experiencing fever. Only four subjects (7%) reported being hospitalized for COVID-19; one required intensive care. The observed case severity distribution was consistent with the general distribution of COVID-19 disease severity in BC (∼5% of diagnosed cases hospitalized as of April 2022) (39).

Participants were required to have recovered from COVID-19 (i.e., 10 days post-PCR diagnosis) before providing their first blood and saline gargle sample. Collection dates ranged from 18-490 days (median 152 days) since a positive PCR test (used as proxy for time since infection), with the baseline collection date ranging from 18-339 days (median 114 days). Participants submitted between 1 and 7 samples, with approximately 2 weeks (median 14 days; range 7–83 days) between each collection, with an average of 6 samples collected per participant and a total of 341 samples collected. No reinfections or persistent virus shedding were detected in self-collected saline gargle samples using RT-PCR (data not shown).

Virus whole genome sequencing was performed (27) on diagnostic samples from 14 participant specimens that were available for analysis, to determine the SARS-CoV-2 variant responsible for infection. SARS-CoV-2 variants were classified as ancestral (e.g., A.1) (n=4) and the D614G mutant (e.g., B.1) (n=10), which is consistent with variants circulating at the time of respective participants’ diagnoses (40). Whole genome sequencing data was missing for ∼71% of participants and; therefore, not included as a covariate in the analysis. Multiple studies corroborate no significant difference in neutralising antibodies between the alpha variant and the ancestral isolate post mRNA vaccination with BNT162b2 or mRNA-1273. Noteworthy reduction of post-vaccination neutralising sera was observed for the beta variant in persons vaccinated with mRNA-1273(41).

3.2 Comparison of anti-SARS-CoV IgG Antibody Responses Across Four Commercial Assays

All available samples (n=340; 1 missing) were initially tested using a combination of three commercial serology assays supplied by Siemens (COV2T), Abbott (ARCHITECT™), or Ortho (VITROS™) clinical diagnostics. Of 340 samples tested, 338 were classified as reactive using at least one assay (Table S1). All available samples (n=339) were subsequently tested using a highly sensitive and multiplex electro chemiluminescent assay offered by Meso Scale Diagnostics (MSD). Percent positivity differed across the platforms and by antigenic target. Overall detection of anti-S was more sensitive than anti-N SARS-CoV-2 IgG. Comparing anti-S results, the Ortho assay had the highest positivity rate (100%) followed by Siemens (95%) and MSD (89%) (Table S1). For anti-N results MSD (58%) outperformed Abbott (47%) with a 11% increase in positivity (χ²test, P=0.01). When samples were stratified by collection time to less than or greater than 6 months post-infection, the anti-N positivity rate decreased for both the Abbott (72% to 13%) and MSD, (76% to 33%) (P<0.001). A 7% decline in positivity was observed for anti-S (P=0.06) and 2% for anti-RBD (P=0.53) when tested by MSD (Table S1). Only antibody measurements from the MSD assay were used in the multivariable analysis as the anti-S IgG results compared well with Ortho and anti-N IgG results were superior to Abbott.

Waning of anti-S and anti-N IgG concentrations over time were measured between participants using the first baseline observation for each of the k=42 participants in the analytic dataset. Using linear regression analysis, overall waning was observed in both anti-N and anti-S and the slope did not differ significantly across the two measures (P = 0.46; Figure 1). On average SARS-CoV-2 antibodies wane at a rate of -0.0029 log10 AU/mL per day (P <0.001) or ∼ 4228 AU/mL per month. These results confirm waning of anti-SARS-CoV-2 antibodies over time in people who have recovered from natural SARS-CoV-2 infection before vaccination. Estimates of anti-ARS-CoV-2 IgG waning are calculated post-vaccination using a mixed-effects linear regression model and reported as the ‘primary analysis’.

• Download figure
• Open in new tab

Figure 1. Longitudinal decay of SARS-CoV-2 anti-N and anti-S IgG concentration over time in natural SARS-CoV-2 infected CARE participants prior to vaccination (k= 42, n= 42 samples).

Participant samples were restricted to the first collection date (baseline) and plotted independently. Linear regression was used to estimate the decrease in anti-S and anti-N titre over time since PCR test result.

3.3 Serological Response to SARS-CoV-2 Vaccination

Bivariate analysis was conducted on the analytic dataset to compare participant antibody responses pre- and post-vaccination for COVID-19. Participant’s serology results and survey responses are summarized and stratified at baseline by the exposure of interest, one dose of a COVID-19 vaccine (Table 1). No difference in the distribution of covariates between participants who received and did not receive a COVID-19 vaccine over the study period was observed for all variables except the number of participant visits. Though follow-up time did not significantly differ between the two groups, on average unvaccinated participants were observed 0.95 (approximately one) fewer times than those who received a COVID-19 vaccine (P=0.014) (Table 1). Importantly, age, biological sex, days from positive PCR test (diagnosis), symptom duration and endemic anti-coronavirus IgG signals did not differ by exposure at baseline; therefore, we expect limited confounding from these covariates when estimating the association between COVID-19 vaccination and anti-SARS-CoV-2 IgG signals. Covariates, which met the definition of a confounder by the common cause criteria, were adjusted for in the primary analysis using a linear mixed effects model.

In k=21 paired participants, SARS-CoV-2 anti-S and anti-RBD IgG antibody concentrations increased post vaccination by 1.63 (P≤0.001) and 1.82 (P≤0.001) log10 AU/ml (Figure 2A & 2B). Anti-N antibody concentration continued to decrease post vaccination by -0.3 (P=0.03) log10 AU/ml (Figure 2C), consistent with waning observed prior to vaccination. Most participants (>99%) were found to be seropositive for anti-S antibodies against the endemic HCoVs. Post vaccination, anti-S antibody concentrations for endemic human -coronaviruses HCoV-HKU1 and HCoV-OC43 slightly increased (P=0.11 and P=0.07) (Figure 3B & 3D). No increase in antibody concentration was observed for the endemic human -coronaviruses HCoV-229E and HCoV-NL63 (P=0.43 and P=0.67) (Figure 3A & 3C).

• Download figure
• Open in new tab

Figure 2. SARS-CoV-2 anti-IgG pre- and post-vaccination.

Antibody signals in k=21 paired participants, who received a COVID-19 vaccine during the study, at baseline and 14 to 56 days post-vaccination, presented by individual SARS-CoV-2 antigen (k=21): A) anti-S, B) anti-RBD, C) anti-N. Differences in antibody signals were examined with a paired t-test.

• Download figure
• Open in new tab

Figure 3. Endemic human coronavirus anti-S IgG antibody signals pre- and post-vaccination.

HCoV antibody signals in n=21 paired participants, who received a COVID-19 vaccine during the study, measured at baseline (before vaccination) and 14 to 56 days post-vaccination, presented by HCoV species: A) HCoV-229E anti-Spike (S), B) HCoV-HKU1 anti-S, C) HCoV-NL63 anti-S, D) HCoV-OC43 anti-S. Difference in antibody signal was examined with a paired t-test.

3.4 Primary Analysis

Linear mixed-effects regression models were used to estimate intraclass correlation within participant samples and the relationship between COVID-19 vaccination and SARS-CoV-2 anti-S or anti-N IgG antibody concentration. An unconditional mean model was fit to partition within participant variation from between participant variation (Table 2). The minority of variation in SARS-CoV-2 anti-S IgG concentration was attributable to differences between participants (ICC=0.43) (Table 2). On average, anti-S IgG concentration increased over time in participants who received one dose of a COVID-19 vaccine during the study by 2.06 log10 AU/mL (95%CI: 1.45-3.46) adjusting for age, biological sex, days from positive PCR test (time) and effect modification between COVID-19 vaccination and time (Table 2). In the adjusted model, the ICC increased to 0.89 indicating that between participant differences (e.g., COVID-19 vaccination) explains most of the variation in SASRS-CoV-2 anti-S IgG antibody concentration. COVID-19 vaccination has a positive association with SARS-CoV-2 anti-S IgG antibody concentration, which increases over time. Variation in anti-N IgG concentration was due to differences between participants in the unconditional mean model (ICC=0.88). The average, anti-N IgG concentration decreased in vaccinated participants over time (−0.243 log10 AU/mL, 95%CI: -1.2 – [=0.12]) adjusting for age, biological sex, days from positive PCR test (time) and effect modification between COVID-19 vaccination and time (Table 2). Variation in anti-N IgG concentration after fitting the adjusted model was explained by within participant variance (ICC=30). Overall, these results indicate that waning of SARS-CoV-2 anti-N IgG is unaffected by COVID-19 vaccination. Anti-S IgG titers increase post vaccination; therefore, vaccination of recovered individuals benefits the durability of their humoral immune response.

3.5 Secondary Analysis

Secondary analysis used the Kaplan-Meier method to estimate the cumulative incidence of seroreactivity above a defined threshold in vaccinated and unvaccinated participants over time. Seroreactive status was classified by the threshold of ≥ 5.5 log10 AU/mL SARS-CoV-2 anti-S IgG, as described in Methods. Participants with antibody measurements equal to or greater than the threshold were considered reactive. Over the 105 days follow up from baseline (first antibody measurement available for participants post-infection), 88% (95%CI: 42-98%) of vaccinated participants (n=16) were seroreactive compared to 5% (95%CI: 0-14%) of unvaccinated participants (n=1) (P= 0.03) (Figure 4). A single dose of COVID-19 vaccine increases the probability of a SARS-CoV-2 anti-S IgG antibody concentration ≥ 5.5 log10 AU/mL by 83%.

• Download figure
• Open in new tab

Figure 4. Cumulative incidence of seroreactivity (≥ 5.5 SARS-CoV-2 anti-S IgG Log10 AU/mL) days from participant’s first blood draw at baseline, stratified by vaccination status over the study period.

Vaccinated participants achieved antibody titres not possible from natural infection alone (unvaccinated participants). Within 105 days of follow up, 88% (95%CI: 42-98%) of vaccinated participants were seropositive, an increase of 83% in comparison to the unvaccinated group (P= 0.03). In previously naturally infected individuals, COVID-19 vaccination increases SARS-CoV-2 anti-S IgG concentration over time to levels which are not attained by natural infection alone. No re-infections were detected by qRT-PCR in the vaccinated or unvaccinated group during the study period, specimens were self-collected.

4.1. Summary of Results

A prospective cohort study was carried out in British Columbia to observe anti-SARS-CoV-2 and anti-endemic HCoV antibody dynamics in participants who were infected with SARS-CoV-2, a subset received the first dose of a Health Canada approved SARS-CoV-2 vaccine during the follow up. Several commercial serology assays were used to detect anti-Coronavirus antibodies; detection of anti-SARS-CoV-2 antibodies was confirmed in all available samples, although both anti-S and anti-N antibodies declines over time post-infection. Bivariate analysis found that vaccination significantly increased the titer of SARS-CoV-2 anti-S IgG antibodies 14-56 days post vaccination; no significant association was found between SARS-CoV-2 vaccination and endemic HCoV anti-S IgG antibodies, although antibodies against the S protein of beta-coronaviruses trended upwards, unlike those against the S protein of alpha-coronaviruses. Vaccination was not observed to boost SARS-CoV-2 anti-N IgG titers, which waned overtime in both vaccinated and unvaccinated individuals. The rate of anti-N waning was approximately double that of anti-S. Secondary analysis used a Kaplan-Meier model to estimate the cumulative incidence of anti-S antibody titers equal to or above 5.5 log10 AU/mL (‘seroreactivity threshold’) in those vaccinated and unvaccinated. In the vaccinated group, 88% (95%CI: 42-98%) of participants had SARS-CoV-2 anti-S IgG titers greater than or equal to this threshold, while this level was achieved in only one unvaccinated participant measured twenty-seven days post infection.

Despite overall antibody waning in unvaccinated participants, a few substantial increases in antibody levels were observed. No reinfections were confirmed using qRT-PCR in self-collected saline gargle samples throughout the study, however one participant had a large (>8-fold) average increase in mean antibody levels (anti-S, RBD and N IgG) seven months following initial SARS-CoV-2 diagnosis, which may be explained by a second exposure to SARS-CoV-2. A second participant had 6-fold increase of anti-S and anti-RBD levels, but not anti-N IgG levels, suggesting they may also have been reexposed. Other detected increases in antibody levels were of much smaller magnitude and might be secondary to rising titers early in convalescence or be explained by technical variations rather than a biological mechanism.

4.2. Comparison with Literature

Previous studies have measured changes in SARS-CoV-2 antibody titers over time. Repeated exposure to SARS-CoV-2 antigens increases IgG titer, while antibodies generated from a single exposure wane overtime (42, 43). Following infection, SARS-CoV-2 specific antibody waning has been observed to decrease from the 8^th to 9^th week post symptoms onset, with detectable levels observed up to the end of the 12^th week (44). In those with multiple SARS-CoV-2 exposures or a hybrid immune response from infection and vaccination, decrease of antibody titers stops shortly after the secondary antigen exposure when stimulation of the memory B cell response produces additional antibodies (45). A strong correlation between total lymphocyte count and SARS-CoV-2 anti-S IgG provides evidence that an ongoing/active immune response provides better protection than a dormant one (44). Waning of SARS-CoV-2 antibodies differs by their paratope, anti-N IgG antibodies wane faster than anti-S. The difference in reactivity between anti-N to anti-S IgG was observed at the population level, anti-N seroprevalence underestimated the number of confirmed infections by 9-31% (46). Vaccination post SARSCoV-2 infection prevents waning of anti-S but has no effect on anti-N IgG (47). Hybrid immunity also benefits the breadth of the antibody mediated response, increasing the probability that existing antibodies are effective against the novel variants. Persons who were infected prior to receiving one of two doses of a COVID-19 vaccine had more somatic mutations and antibody production from the IGHV2-5; IGHJ4-1 germline which was not active in the vaccinated but uninfected (48, 49). Additionally, hybrid immunity produces greater total and neutralizing anti-S titers than natural infection or vaccination alone (45). Our study both supports and builds upon prior findings, as we show that SARS-CoV-2 IgG antibodies wane in SARS-CoV-2 infected people over time, with the rate of decline being greater for anti-N IgG than anti-S; vaccination post-infection boosts anti-S IgG titers and participants with hybrid immunity possess anti-S antibody levels which are not common in those infected but unvaccinated. Our calculated rate of antibody decline may be used to help estimate infection timing in seroprevalence studies.

4.3. Clinical and Epidemiological Interpretation

Our findings have important implications for clinical practice and public health guidelines as the pandemic progresses into its third year, novel viral variants continue to emerge, and vaccine doses are more widely distributed globally. Humoral immunity from natural infection wanes and vaccination with at least one dose of COVID-19 vaccine increases SARS-CoV-2 anti-S IgG titers immediately and over time. Therefore, we recommend that naturally infected individuals receive COVID-19 vaccination to increase protection from re-infection and severe disease and the duration of their humoral immune response against SARS-CoV-2. We demonstrate that a single dose of SARS-CoV-2 vaccine is effective in boosting anti-S antibody titers to high levels, which has implications in distribution of vaccine supplies in those countries with scarce access and low vaccination levels in the setting of high numbers of natural infection.

4.4 Strengths and Limitations

The strength of the described study stems from the prospective design, use of multiple serological tests, including the quantitative MSD option, and thorough analysis. A prospective cohort design offers several benefits, which allowed us to observe SARS-CoV-2 antibody dynamics over time with minimal bias. Recruiting participants post-infection but prior to vaccination delineated the sequence of temporal events, limiting the probability that any changes in antibody titers observed post-vaccination were due to causes other than the vaccine. Selection bias was minimized as the participants exposure and outcome status were not known when they were recruited into the study. At the beginning of the study, the measured covariates were exchangeable between participants who were unvaccinated or vaccinated during follow-up. Balance of the covariates allowed for estimation of the relationship between vaccination and anti-SARS-CoV-2 IgG antibody titer with minimal bias from confounding. Utilizing multiple serological tests strengthened our observations, as well as allows generalizability to study populations in many different laboratories. Statistical power was optimized by analysis with a mixed effects linear regression model, which accommodated multiple repeated measures per participant.

Limitations of the work include differential loss to follow up in the vaccinated and unvaccinated groups, a small sample size, and incomplete/missing survey responses. Unvaccinated participants were observed to have approximately one fewer visit than those who received a COVID-19 vaccine. Vaccines were not an originally planned intervention in the study and were made available in British Columbia on a stage roll-out basis about half-way through the study period. The difference in visit numbers between the vaccinated and unvaccinated groups is likely related to surveillance biasthose who receive a medical intervention are more open to clinical follow up than those who do not. Obtaining a larger sample size initially planned for the study was difficult due to low enrollment uptake, likely related to the social and economic stress of the pandemic on the public and geographic limitations on recruitment related to the availability of sample collection sites.

While the reached sample size was adequate for the primary analysis-which included repeated measures-enrolling additional participants would have allowed for a more precise estimate of the relationship between COVID-19 vaccination and endemic HCoVs anti-IgG titers. Although not significant, a trend was observed for increased average antibody titers for HCoV-HKU1 and HCoV-OC43 following vaccination which, like SARS-CoV-2, are both β-coronaviruses in contrast to decreases in the averages observed for both HCoV-NL63 and HCoV-229E which are α-coronaviruses. The study design may have also underestimated any association between existing endemic coronavirus IgG titers and COVID-19 vaccination as the sample was restricted to persons previously infected with SARS-CoV-2. COVID-19 has been shown to affect endemic coronavirus antibody levels and as such, the effect of vaccination should be observed in a cohort of SARS-CoV-2 naive persons prior to vaccination (50–52). The overall effect of SARS-CoV-2 vaccination and/or infection on the circulating antibodies against endemic HCoVs in the population may have implications for their seasonal epidemiology.