Longitudinal Study on Seroprevalence and Immune Response to SARS-CoV-2 in a Population of Food and Retail Workers Through Transformation of ELISA Datasets

1. Cohort composition and sample collection

All recruitment, data collection interviews and sample collection activities were carried out at the Centre Hospitalier Universitaire (CHU) de Québec-Université Laval and have been described in detail [39]. Briefly, the cohort included 304 food and retail workers from the Capitale Nationale (76%) and Chaudière-Appalaches (24%) administrative regions of the Canadian province of Quebec. Individuals provided written informed consent (approved by the “Comité d’éthique de la recherche du CHU de Québec-Université Laval”, registration number 2021-5744). The Canadian Immunity Task Force (CITF) questionnaire [60] was used for data collection and data collected during interviews have been reported [39]. With an average age of 41.3 years, 58 % were female and 42% were male; nearly half (49%) worked in restaurants or bars, 37% in grocery stores and 14% in hardware stores. None had a history of hospitalization linked to SARS-CoV-2 [39].

The observation period spanned 1.5 years, from April 2021 to October 2022. A total of 121 participants reported at least one positive test during the study period. Participants attended five evenly spaced visits between April 2021 and October 2022. A blood sample and information on SARS-CoV-2 risk factors, symptoms, antigen or PCR test results, and vaccination status were collected at every visit. In addition, information regarding potential risk factors (i.e., demographic, socioeconomic, behavioral, clinical, and occupational) was obtained, allowing multifactorial analysis of the cohort.

At each visit, a 30 mL blood sample was collected in 6 mL tubes (BD Vacutainer 367815). Sera collected between June 4^th and 9^th, and between July 2^nd and 8^th, 2020, from eight individuals (aged between 20-55 years old; 7 females and 1 male) who had never tested positive for SARS-CoV-2 were pooled and used as the negative control. Sample tubes were gently inverted, held at room temperature for 15–30 min to allow clotting and spun at 1600 g for 15 min. Serum samples (1 mL aliquots) were transferred into cryovials (Sarstedt Inc., product 72.694.006), frozen in an upright position at -20°C and stored at -80°C until shipment on dry ice to the assay site, where they were maintained at -80°C until use.

2. SARS-CoV-2 viral antigens

The National Research Council of Canada (NRC) provided trimeric recombinant SARS-CoV-2 spike proteins (Wuhan-Hu-1 (SmT1, lot PRO1-429), Delta (SmT1v3 (B.1.617.2)), and Omicron BA.1 (SmT1(BA.1), lot PRO7911-2]), as well as the spike receptor-binding domain (319–541 RBD) produced as described previously.[61,62] The NRC also supplied the SARS-CoV-2 nucleocapsid NCAP-1.[63]

3. In-house ELISA (Enzyme-linked Immunosorbent assays)

The semi-automated in-house colorimetric ELISA was performed as described in [45] and [37] with some modifications by adapting a semi-quantitative ELISA protocol from previous work [64–66]. Briefly, 96-well immunoassay plates were coated with 100 μL of the relevant SARS-CoV-2 antigen and incubated overnight at 4 °C. They were then washed and blocked with 100 μL of 3% (w/v) skimmed milk powder for 1h at room temperature (RT) and then washed again. The human serum samples were inactivated by incubating them in a block heater at 56 °C for one hour and diluted as specified (1: 15,000 for anti-spike IgGs, 1: 200 for anti-nucleocapsid IgGs, and 1: 100 for anti-spike and anti-nucleocapsid IgMs). Aliquots (100 μL) of the inactivated, diluted human serum samples were then added in the prepared immunoassay plates and incubated at RT for one hour. The plates were washed, and 100 μL of the appropriately diluted host-specific secondary antibody was added (1: 30,000 Goat anti-human IgG HRP Life Technologies (Invitrogen) Catalog #31413 or 1: 10,000 Goat anti-human IgM (µ-chain specific) HRP Sigma-Aldrich Catalog # A6907), covered to block out light, incubated at RT for 1h, and then washed again. 100 μL of 3,3ʹ,5,5ʹ-tetramethylbenzidine (TMB) was then added and incubated at RT for 1h, followed by the addition of 100 μL of 2 M HCl for color development. Absorbance was measured at 450 nm.

Quality control was performed using a calibration standard to ensure inter-day reproducibility throughout the testing period. The calibration standard consisted of pooled clinical samples from participants who were confirmed COVID-19-positive by PCR test. Independent replicates (n = 45, performed on different assay plates over several months) performed for each antigen and immunoglobulin type determined the standard signal range for each antigen. The standard was included on each ELISA plate (3 wells on each 96-well plate). When the value obtained for the standard fell outside of the acceptable range (+/- 25% of the median value), assay results for that plate were excluded. We note that the median value was determined retrospectively, in that the 45 independent replicates were performed alongside the query samples. According to that median, we then retrospectively evaluated which 96-well plates should be excluded. Those assays were repeated. The COVID-19-negative control (described above) was also included on each ELISA plate.

Independent replicates performed for each specific antigen and immunoglobulin type (n = 45, performed on different assay plates over several months) established the negative control range for each target. Positivity thresholds were established as values above the mean of negative controls for a given target plus one standard deviation. Again, positivity thresholds were determined retrospectively.

4. ELISA (Enzyme-linked Immunosorbent assays) at the uOttawa High Throughput Serology and Diagnostic Facility

The automated chemiluminescent ELISA, referred to as the ‘centralized assay’, was described in [36]. Briefly, automated chemiluminescent ELISAs were performed using MicroLab Star robotic liquid handlers (Hamilton, USA) and a 405 TS/LS LHC2 plate washer (Biotek Instruments; all wash steps included four washes with 100 μL PBST). All incubations were done at room temperature with shaking at 500 rpm. Antigens (Wuhan-Hu-1 ancestral spike, 319–541 RBD, and NCAP-1 nucleocapsid) were diluted in PBS and dispensed into the wells of a 384-well high-binding polystyrene Nunc plate (Thermo Fisher Scientific, #460372) at a final amount of 50 ng/well. The plates were centrifuged at 216 × g for 1 min to ensure even coating, incubated overnight, rocking at 4°C, and washed. Wells were blocked with 80 μL of 3% w/v skim milk powder dissolved in PBST for 1 h and then washed. Samples and controls were diluted as indicated to a final concentration of 1% w/v skim milk powder in PBST, and 10 μL was added to each well from a 96-well source plate. Plates were incubated for 2 h, and wells were washed. Secondary antibodies (as in the colorimetric assay) were diluted as indicated in 1% w/v skim milk powder in PBST, and 10 μL was added to each well. After incubation for 1 h, the wells were washed, and 10 μL of ELISA Pico Chemiluminescent Substrate (Thermo Scientific; diluted 1:2 in MilliQ H₂O) was dispensed into each well. After a 5 min incubation with shaking, plates were read on a Neo2 plate reader (BioTek Instruments) at 20 ms/well and a read height of 1.0 mm [36].

5. Data transformation and statistical analysis

An overview of seropositivity and IgG levels in food service and retail workers from April 2021 to October 2022

Higher binding antibody titer in infected or vaccinated individuals has been correlated with increased protection efficiency [6–12], allowing the evaluation of a population’s protection level and predicting the longevity of protection [7]. During the 1.5-year study period, blood samples from 304 food and retail workers were collected at five visits, each spanning a period of 4 to 7 months (Fig. 1). In total, 1299 samples were collected. The cohort composition is summarized under Methods, and complete details of the recruitment method and data collected by interview have been reported [39]. The percentage of fully vaccinated participants (defined by the Canadian government as having received at least two doses of one or a combination of approved vaccines [43]) was already 48% by the first visit and reached 97% at visit 5 (Table S2). Samples collected during the study period were analyzed with the following objectives: (i) to determine the trend of humoral immunological response over time and in response to critical pandemic events, such as the emergence of new variants; (ii) to identify any relevant statistical differences related to the workers’ occupations, their vaccination status, protective measures at their workplace as well as their underlying health conditions and lifestyle.

• Download figure
• Open in new tab

Figure 1. Overview of ELISA from all samples assayed over 18 months (n = 1290).

A) In-house ELISA, reporting the mean monthly O.D.₄₅₀. IgG levels for the ancestral, delta and Omicron spike protein variants and the nucleocapsid were measured (see colored legend). Positivity threshold: ancestral spike = 0.031; delta spike = 0.033; Omicron spike = 0.040; nucleocapsid = 0.048. B) Centralized ELISA, reporting the mean antibody units per mL (BAU mL^-1). IgG levels for the ancestral spike protein and the nucleocapsid were measured (see colored legend). Visits at which samples were taken are indicated in different colors under each graph. The approximate availability of vaccine doses to the general population in the province of Québec is indicated by arrows [47]. Note the scales on the y-axes of (B): left axis for ancestral spike, right axis for nucleocapsid.

It is advisable to use the spike protein as a reference in population-based seroprevalence studies, as the humoral immune responses it elicits are strong [44] and it is the antigen used for nearly all approved COVID-19 vaccines. IgG and IgM levels specific to the spike protein ectodomain, as well as the SARS-CoV-2 nucleocapsid protein, were determined for the ancestral strain and the Delta and Omicron VOCs, using a previously established in-house ELISA [37] [45]. The use of a calibration standard made up of pooled clinical samples from PCR-diagnosed COVID-positive participants ensured reproducibility. The standard was included on each ELISA plate, enabling the retrospective establishment of a standard signal range (median +/- 25%) for each antigen. Only plates where the standard signal was within this range were considered valid; invalid assays were repeated as needed. A pool of eight pre-COVID-19 serum samples was used as a negative control and was tested with each specific antigen and immunoglobulin type.

Our objective was first to provide a YES/NO answer to SARS-CoV-2 vaccine seroconversion or infection and then to report on statistical differences in antibody response levels between populations of workers in our cohort. Positivity thresholds (YES) were established as values above the mean of negative controls for a given target plus one standard deviation. The positivity thresholds were determined retrospectively after most samples had been tested (see Methods). The ELISA data were either used as is (YES/NO response) or further transformed into relative concentrations (R_c) by fitting them to a logistic function obtained using a calibration curve for each antigen for statistical analysis. [38]

Validation of the YES/NO results obtained with the in-house ELISA was achieved by analyzing the same clinical samples at the centralized testing site, using antigens from the same source. Seropositivity was determined at the centralized testing site by ELISA to IgG against both the ancestral Wuhan-Hu-1 strain spike protein ectodomain and receptor-binding domain (RBD) and against the nucleocapsid protein [24,36]. A World Health Organization (WHO) reference standard was used to calibrate the assay [36]. The scaled luminescence values were measured and converted to titer in the form of antibody units per mL [24].

For statistical analysis, calibration curves were obtained by serial dilution of the standard [37,38]. As the standard solution was of unknown concentration, the most diluted sample was arbitrarily set at a concentration of 1, hence the term “relative concentration”. This transformation made it possible to normalize and statistically analyze the data using Welch’s ANOVA, which is only applicable to normally distributed data, yet allows for variance inequality (see Methods) [38].

Seroprevalence in the population of workers over the timeline of the study

The monthly means of anti-spike (ancestral Wuhan-Hu-I strain) and anti-nucleocapsid IgG levels obtained with the in-house ELISA (Fig. 1A, reported as mean O.D. values) were compared with the centralized ELISA (Fig. 1B, reported as mean BAU/mL). Concordant trends of seroprevalence were observed by both methods for ELISA to the ancestral variant: anti-spike IgG levels peaked around July/August 2021 during Alpha variant dominance, coincident with the 2^nd mass vaccination campaign in the region, and in February to July 2022 during Omicron dominance, coincident with the 3^rd mass vaccination campaigns in the region; these observations are consistent with the trends observed more broadly in the Canadian population over the same period [46]. Concordant trends with the anti-nucleocapsid IgG levels (Fig. 1A and B) were also observed: nucleocapsid signal remained low because relatively few people were infected (Table 1) until the highly infectious Omicron variant was reported in Québec on November 29^th, 2021 [5]. The individual ELISA results are provided in Fig. S1 (in-house ELISA) and Fig. S2 (centralized ELISA).

We further investigated seroprevalence in this population by monitoring the IgM directed against the Wuhan-Hu-1 ancestral spike and nucleocapsid using the in-house ELISA assay. The IgM ELISA signals observed using ancestral spike and nucleocapsid antigens were consistently lower than the corresponding IgG signals, as reported in prior studies (Fig. S3)[37,48]. The pattern of rapid production and decline of anti-SARS-CoV-2 IgM is consistent with other reports [37,48–50]. This confirms the utility of reporting on IgM to monitor populational seroprevalence shortly after vaccination or infection events, while illustrating the greater overall practicality of reporting on IgG over longer observation periods.

Seropositivity in the population of workers over the timeline of the study

The ELISA data allowed determination of seropositivity for each individual, at each visit (Table 1; Fig. S1 and Fig. S2). In this standard YES/NO analysis, concordant patterns of vaccine seroprevalence (anti-spike IgG, ancestral strain) was observed between the in-house and centralized ELISA methods. A strong correspondence was also observed for naturally acquired seroprevalence (anti-nucleocapsid IgG) (Table 1). Furthermore, the anti-nucleocapsid IgG data is fully consistent with infection data, including reports of PCR-positive test or positive self-test, where infection among this cohort was low (fewer than 5 cases per 100 person-years) during the first waves, rapidly increasing with the emergence of Omicron (peaking at 80 cases per 100 person-years during the fifth wave) [39]. These observations reinforce the validity of using in-house ELISA as a reliable instrument for populational studies during the early phases of a pandemic [37].

The scope of the study was broadened by monitoring Delta and Omicron anti-spike IgG levels using the in-house ELISA assay, highlighting the adaptability of this method (Fig. 1A). Before December 2021, during the period where serological response resulted predominantly from vaccination using ancestral spike antigen and before the appearance of the Omicron VOC, the IgG ELISA response using the Delta spike antigen was indistinguishable from that using the ancestral Wuhan-Hu-1 spike antigen. This is consistent with prior reports of effective cross-reactivity between these strains [37,51,52]. As the study progressed, the Omicron variant became dominant. In contrast with the Delta variant, and also consistent with a prior report, cross-reactivity with the Omicron spike antigen was weaker but nonetheless consistent with immune protection even as the viral strains continued to evolve (Fig. 1A and 1B) [53].

Statistical correlation between vaccination and antigen levels

COVID-19 vaccines were made available to the general adult population in Québec beginning in March 2021 [47], coinciding with the beginning of this study. The Astra Zeneca COVISHIELD (Vaxzevria) was rapidly followed by the Pfizer-BioNTech and Moderna Spikevax RNA-based vaccines; all vaccines initially approved for use in the Québec region were based on the ancestral spike antigen.

Over the course of this 1.5-year longitudinal study, vaccination status increased even as long-term participation in the study decreased: at the first visit (V1), 48% (146 of 304 the study participants) were fully vaccinated (had received at least two doses of one or a combination of the approved vaccines [43]), V3 saw 65% (188/291 participants) fully vaccinated and V5 saw 97% (188/194) fully vaccinated (Table S2). We note that the increase in vaccination status from 65% at V3 to 97% at V5 resulted entirely from the drop from 291 to 194 study participants. Interestingly, all who withdrew from the study between V3 and V5 were non-vaccinated. These numbers reveal that vaccination status was strongly correlated with long-term engagement in the study and caution against assigning a significance of the near-perfect vaccination status at V5 to any criterion investigated here.

The in-house ELISA datasets were statistically analyzed to determine correlations of antibody levels with vaccination. Specifically, ELISA reporting on IgG using the Wuhan-Hu-1 ancestral spike, Omicron spike and nucleocapsid antigens were analyzed to determine the correlation between vaccination events and seroconversion, and the impact of viral evolution to the Omicron VOC on cross-reactivity. We note that, following seroconversion for anti-nucleocapsid IgG in any individual, further samples from that individual were excluded to eliminate the impact of infection on the immune response.

We performed statistical analysis of the in-house ELISA dataset by Welch’s ANOVA test. First, mathematical data curation and transformation was performed, fulfilling the requirement that datasets adopt a normal distribution; the unequal variances of the resulting dataset precluded application of ANOVA and indicated Welch’s ANOVA as appropriate. Then, the relationship between ELISA absorbance and relative IgG concentration was modeled, using a COVID-19-positive standard and applying a regression model, as previously reported by others; the data fit a five-parameter logistic function, detailed in Methods [38]. Comparison of ELISA standard curves to the model for each antigen shows excellent concordance (Fig. 2). Having confirmed the accuracy of each model, the in-house ELISA datasets were transformed using their respective antigen model and statistically analyzed by Welch’s ANOVA test.

• Download figure
• Open in new tab

Figure 2. Comparison of serial dilution ELISA data with model prediction.

ELISA data for IgG were compared with model predictions for A) ancestral spike, B) nucleocapsid, C) delta spike and D) Omicron spike antigens. X-axis: dual logarithm of the relative concentration of serial dilutions (experimental data points in black, model predicted points in red). Y-axes: O.D. of serial dilutions.

We observed clear statistical differences in anti-spike IgG levels between unvaccinated individuals and those who have received a greater number of doses, up to a maximum of four doses (Fig. 3). Analysis of cross-reactivity with the Omicron variant spike antigen gave similar results, although the statistical significance was generally greater. When the analysis was repeated with inclusion of samples following seroconversion for anti-nucleocapsid IgG, the trend was conserved with an increase statistical difference (Fig. S4). Although this may be due to the greater number of datapoints analyzed, it is more likely due to the additional immunity-provoking event of the infection in the nucleocapsid-seroconverted individuals. In addition to providing an additional event of antigen exposure, vaccination and infection together have been shown to provide a greater response than either alone [54].

• Download figure
• Open in new tab

Figure 3. Comparative assessment of humoral immune response in unvaccinated and vaccinated, nucleocapsid-negative individuals.

ELISA assays measuring ancestral or Omicron anti-spike IgG levels were analyzed collectively (all visits together). Each grouping includes all data collected (n is given) for that number of vaccine doses (0 to 4, x-axis); ‘n’ excludes samples of individuals following determination of seropositivity for anti-nucleocapsid IgG and includes any sample taken at different visits for the same individual whose vaccination status did not change. X-axis: number of vaccine doses received (minimum 7 days post-vaccination). Y-axis: data points after logarithmic transformation. The median (solid line) and quartiles (dashed lines) are shown in the violin plots; no outliers were removed. Statistical significance: *, p < 0.05; **, p < 0.01; ****, p < 0.0001.

It is important to note that the samples were collected without regard for the dates of vaccination of the individuals. Therefore, the observed serological response is averaged over many months post-vaccination for this population. This is seen in the modest increase (below statistical significance) in signal resulting from two doses of vaccine (‘primovaccination’) compared to non-vaccinated. Since sampling continued for months prior to the 3^rd vaccination campaign, it should not be interpreted as poor serological response as it is subject to the expected waning of humoral response over time following primovaccination [55]. A marked increase in median anti-spike IgG levels in samples following 3 or 4 vaccine doses was observed, consistent with other reports [56]. We note that the dataset included few participants having received a single vaccine dose, as many participants had already been vaccinated twice at their first sample collection visit.

Using the same mathematically transformed ELISA dataset, we monitored for any statistical differences in seroprevalence to the spike antigen as a result of the type of vaccine administered. Two vaccine doses (no matter the combination of vaccine products used) resulted in no statistical difference with the unvaccinated individuals (Fig. 4A). For three vaccine doses, a statistical difference between unvaccinated individuals and fully vaccinated participants was observed for the combination of three RNA-based vaccines (no matter the combination of Moderna/Pfizer-BioNTech products used); cross-reactivity with the Omicron spike antigen gave a higher statistical significance. Although the trend was maintained for the combination of Astra-Zeneca/2RNA vaccines, it was not statistically significant (Fig. 4B). However, upon including the nucleocapsid positive individuals in the analysis (Fig. S5 A and B), greater statistical differences were seen. As above, this could be due to the greater number of datapoints analyzed, and to the additional immunity-provoking event of the infection in the nucleocapsid-seroconverted individuals.

• Download figure
• Open in new tab

Figure 4. Comparative assessment of humoral immune response for different vaccine types.

ELISA assays measuring ancestral or Omicron anti-spike IgG levels were analyzed collectively (all visits together) following administration of A) 2 vaccine doses or B) 3 vaccine doses. Each grouping includes all data collected (n is given) for that type of vaccine; ‘n’ excludes samples of individuals following determination of seropositivity for anti-nucleocapsid IgG and includes any sample taken at different visits for the same individual whose vaccination status did not change. “RNA only”: either the Pfizer-BioNTech Comirnaty® mRNA COVID-19 vaccine, the Moderna Spikevax® mRNA COVID-19 vaccine, or a combination of both. “AZ/RNA” and “AZ/2RNA”: one dose of the AstraZeneca COVISHIELD® viral vector-based COVID-19 vaccine and one or two doses of an RNA vaccine, respectively. Individuals who had received the AstraZeneca vaccine only, or two doses of AstraZeneca and one of an RNA vaccine, were excluded as there was insufficient data to allow statistical analysis. Y-axes: data points after logarithmic transformation. The median (solid line) and quartiles (dashed lines) are shown in the violin plots; no outliers were removed. Statistical significance: *, p < 0.05; ****, p < 0.0001.

Our findings are consistent with the report that vaccine efficacy is roughly comparable, reaching around 90% seroconversion when the two doses administered are from an mRNA vaccine (i.e. Pfizer and Moderna), or when the first dose administered is AstraZeneca (viral vector vaccine) but followed by a second dose with an mRNA vaccine [57].

Statistical analysis of antigen levels in relation to the workplace in an under-reported yet highly exposed population

The past years have seen a worldwide, concerted effort to understand the human immunological response to SARS-CoV-2. Our focus on the humoral immunological response of the underrepresented yet highly exposed group of food service and retail workers has yielded a comprehensive dataset spanning the 18-month period during which vaccination was deployed and the virus evolved into several VOCs. We solicited the participation of individuals working in bars, restaurants, grocery stores and hardware stores because these jobs, by their very nature, require direct contact with customers and are therefore considered at greater risk than other types of work. Workers in these sectors in Québec often have below-average incomes and no health-related benefits [39]. A certain variability of risk exists in this cohort: while grocery stores and hardware stores remained open throughout the pandemic in Québec, bars and restaurants were periodically closed by health authorities. However, the use of protective equipment was better enforced in grocery and retail stores than in restaurants and bars, where clients stayed longer and removed masks to consume food and beverages [39].

We observed no overall statistical difference in humoral response between the different occupations according to ELISA when the data from all visits were analyzed collectively (all visits together) (Fig. 5A and B). However, certain trends emerged upon stratification by visit (Fig. 5C and D), specifically over the period of the two first visits, during which the cohort mostly became fully vaccinated (Fig. 1) but where the rate of infection-acquired immunity was low (Table 2). At visit 1, ancestral strain anti-spike IgG levels were statistically higher for the individuals working in restaurants than in any of the three other areas of occupation (Fig. 5C); at visit two, anti-spike IgG levels were statistically higher in those working in restaurants and hardware stores. Bars were closed in the region investigated for 3 of the 10 months spanning visit 1, coincident with evening curfews; this measure of health protection appears to have been successful in reducing exposure of bar workers relative to the other areas of occupation as nucleocapsids levels for this group are lower than for the others. In contrast, grocery and hardware stores were open throughout visit 1; restaurants were also open, although they were restricted to limited capacity by periods. Interestingly, while restaurant and hardware store workers had the highest level of anti-nucleocapsid IgG response at visit 1, bar and grocery store workers had higher levels from visit 2 on (although only statistically significant for visit 2 itself) (Fig. 5D).

• Download figure
• Open in new tab

Figure 5. Effect of different occupations on the humoral immune response.

ELISA measuring A) ancestral anti-spike or B) anti-nucleocapsid IgG levels analyzed collectively (all visits together), or C) ancestral anti-spike and D) anti-nucleocapsid IgG levels stratified by visit (V1 to V5) for serum samples collected during the study period. Results are shown in logarithmic form, post-transformation. The median (solid line) and quartiles (dashed lines) are shown in the violin plots; no outliers were removed. The number of data points (n) per plot is shown. Statistical significance: **, p < 0.01; ***, p < 0.001; ****, p < 0.0001.

An important and contentious subject throughout the pandemic has been whether wearing a mask contributes to protection against exposure to the virus. There is some evidence in the literature that masks provide useful protection [58]. We analyzed self-reported mask-wearing habits at the workplace, finding no statistical difference in the level of anti-nucleocapsid IgG (infection-induced immunity) for all visits, with the exception of visit 4 (Fig. 6). Data for visit 4 were collected during the peak of the highly infectious Omicron variant (Fig. 1). Here, more people (65) reported not wearing masks than in previous visits (5 in visit 1, 11 in visit 2 and 21 in visit 3), and their anti-nucleocapsid IgG levels were significantly higher than those of people wearing masks during the same period. We hypothesize that the incidence of the Omicron variant and more relaxed mask-wearing habits contributed to the observed difference.

• Download figure
• Open in new tab

Figure 6. Impact of wearing a mask in the workplace on humoral immune response.

ELISA measuring A) ancestral anti-spike and B) anti-nucleocapsid IgG levels analyzed by visit (V1 to V5). Results are shown in logarithmic form, post-transformation. The median (solid line) and quartiles (dashed lines) are shown in the violin plots; no outliers were removed. The total number of data points (n) per plot is shown. Statistical significance: ****, p < 0.0001.

We found no significant differences for participants who used public or private transport to get to and from the workplace (Fig. S6) or for those who reported being in contact with infected people (Fig. S7).

Statistical findings unrelated to work-related activities

While the main focus of this study was to analyze the effects of vaccination and exposure for this cohort, we had the opportunity to examine other factors thanks to the information provided by participants during the survey at each visit, including level of education (Fig. 7A), chronic illnesses (Fig. 7B), age (Fig. 7C), body mass index (BMI) (Fig. 7D and Fig. S8) and participants’ smoking habits (Fig. 7E).

• Download figure
• Open in new tab

Figure 7. Impact of non-occupational activities on humoral immune response.

ELISA measuring ancestral anti-spike or anti-nucleocapsid IgG levels are analyzed collectively (all visits together) and separated by A) level of education, B) number of chronic illnesses, C) age, D) BMI and E) smoking status. Results are shown in logarithmic form, post-transformation. The median (solid line) and quartiles (dashed lines) are shown in the violin plots; no outliers were removed. The total number of data points (n) per plot is shown. Statistical significance: *, p < 0.05. In (D), two data points were removed because the category (BMI < 18) was not sufficiently populated to allow statistical analysis.

No statistically significant differences were found for any of the factors analyzed, except a statistically higher ancestral anti-spike IgG response in non-smokers than in smokers (p < 0.05). This difference was not observed in the case of natural immunity (anti-nucleocapsid IgG). This is consistent with prior observations that smokers generally show lower antibody titers or a more rapid decline of vaccine-induced IgG levels than non-smokers [59]. The results shown above were analyzed collectively (all visits together); no significant difference was observed when statistical analysis was stratified by visit.