key: cord-0761198-6yio5mjl
authors: Sheikh-Mohamed, Salma; Isho, Baweleta; Chao, Gary Y. C.; Zuo, Michelle; Cohen, Carmit; Lustig, Yaniv; Nahass, George R.; Salomon-Shulman, Rachel E.; Blacker, Grace; Fazel-Zarandi, Mahya; Rathod, Bhavisha; Colwill, Karen; Jamal, Alainna; Li, Zhijie; de Launay, Keelia Quinn; Takaoka, Alyson; Garnham-Takaoka, Julia; Patel, Anjali; Fahim, Christine; Paterson, Aimee; Li, Angel Xinliu; Haq, Nazrana; Barati, Shiva; Gilbert, Lois; Green, Karen; Mozafarihashjin, Mohammad; Samaan, Philip; Budylowski, Patrick; Siqueira, Walter L.; Mubareka, Samira; Ostrowski, Mario; Rini, James M.; Rojas, Olga L.; Weissman, Irving L.; Tal, Michal Caspi; McGeer, Allison; Regev-Yochay, Gili; Straus, Sharon; Gingras, Anne-Claude; Gommerman, Jennifer L.
title: Systemic and mucosal IgA responses are variably induced in response to SARS-CoV-2 mRNA vaccination and are associated with protection against subsequent infection
date: 2022-04-25
journal: Mucosal Immunol
DOI: 10.1038/s41385-022-00511-0
sha: c67d0959612d2647284a25a54b82c8903cee1c1e
doc_id: 761198
cord_uid: 6yio5mjl

Although SARS-CoV-2 infects the upper respiratory tract, we know little about the amount, type, and kinetics of antibodies (Ab) generated in the oral cavity in response to COVID-19 vaccination. We collected serum and saliva samples from participants receiving two doses of mRNA COVID-19 vaccines and measured the level of anti-SARS-CoV-2 Ab. We detected anti-Spike and anti-Receptor Binding Domain (RBD) IgG and IgA, as well as anti-Spike/RBD associated secretory component in the saliva of most participants after dose 1. Administration of a second dose of mRNA boosted the IgG but not the IgA response, with only 30% of participants remaining positive for IgA at this timepoint. At 6 months post-dose 2, these participants exhibited diminished anti-Spike/RBD IgG levels, although secretory component-associated anti-Spike Ab were more stable. Examining two prospective cohorts we found that participants who experienced breakthrough infections with SARS-CoV-2 variants had lower levels of vaccine-induced serum anti-Spike/RBD IgA at 2–4 weeks post-dose 2 compared to participants who did not experience an infection, whereas IgG levels were comparable between groups. These data suggest that COVID-19 vaccines that elicit a durable IgA response may have utility in preventing infection.

INTRODUCTION SARS-CoV-2 is a novel and highly contagious respiratory virus that has quickly spread across the globe. The virus uses a protein called Spike and its associated receptor binding domain (RBD) to interact with angiotensin converting enzyme 2 (ACE2) on host cells 1 . Interaction between viral Spike/RBD and ACE2 on the cell surface is the first essential step in SARS-CoV-2 infection, and expression of ACE2 on epithelial cells of the upper respiratory tract (URT) renders them susceptible to aerosolized virus. Thus, immunity in the oral and nasal mucosa is an important first line of defense against the development of COVID-19 2 .

Saliva is an important biofluid that can provide information about the mucosal antibody (Ab) response to SARS-CoV-2 3 . Indeed, salivary gland epithelial cells express ACE2 and harbor a significant population of IgA-producing plasma cells 4 . Secretory IgA (SIgA) in the saliva exists as IgA dimers that are associated with the secretory component, a proteolytic cleavage product which remains bound to IgA after it is transported across epithelial cells via the polymeric Ig receptor (pIgR) 5 . Secretory polymeric IgA has been shown to have potent neutralizing activity against SARS-CoV-2 in vitro 6 .

We and others have shown that IgM, IgG and IgA Ab against the SARS-CoV-2 Spike and RBD proteins are readily detected in the saliva of COVID-19 acute and convalescent patients 3, 7 . Whether COVID-19 vaccines delivered through the parenteral intramuscular route (i.m.) generate a similar salivary antibody response is unclear, and the nature and kinetics of this response are illcharacterized. Given the importance of mucosal immunity as a first line defense against SARS-CoV-2 infection we measured Spike/RBD-specific Ab in saliva samples from participants who had received either BNT162b2 (Pfizer/BioNTech) or mRNA-1273 (Moderna) vaccinations. We also determined whether levels of vaccine-induced anti-Spike/RBD IgG or IgA differed in people who subsequently experienced a SARS-CoV-2 infection. Collectively, our data show that a SIgA response is induced in~30% of participants who received 2 doses of a SARS-CoV-2 mRNA vaccine, and that IgA may play an important role in protection against infection.

Detection of anti-Spike and anti-RBD antibodies in saliva from participants receiving COVID-19 mRNA vaccines We first compared saliva from long-term care home (LTCH) workers that received either BNT162b2 or mRNA-1273 (Supplementary Fig. 1 Table 2b ). We expressed the data as a percentage relative to a pooled sample of saliva from COVID-19 acute and convalescent patients-the same pooled sample was present in each plate. We previously found that this method provided excellent plate-to-plate consistency and produced similar results as what we found when we normalized to total IgG/IgA 3 . Only 11% and 22% of vaccinated participants had a detectable IgM response to Spike and RBD respectively (data not shown). Focusing therefore on IgG and IgA responses, after two doses of mRNA vaccine 94% and 41% of participants were positive for anti-Spike IgG and IgA, and 93% and 20% of participants were positive for anti-RBD IgG and IgA Ab. While levels of anti-Spike/ RBD IgG were similar to or exceeded that of COVID-19 convalescent patients (Fig. 1A, B) , IgA levels were significantly lower (Fig. 1C, D) . Furthermore, as we observed before in COVID-19 recovered patients 3 , levels of salivary anti-Spike/RBD Ab positively correlated with anti-Spike/RBD Ab in the serum ( Supplementary Fig. 2 ). In multivariable analysis, age and prior SARS-CoV-2 infection were independently associated with the salivary anti-Spike IgA response (Supplementary Table 3a ). In contrast, male sex had a negative independent association with the salivary anti-Spike IgG response (Supplementary Table 3b ) as has been observed before for COVID-19 and other vaccines 8, 9 . Lastly, prior SARS-CoV-2 infection and time since vaccination were independently associated with higher and lower serum anti-RBD IgA levels, respectively (Supplementary Table 3c ).

Longitudinal assessment of anti-Spike and anti-RBD antibodies in saliva from participants receiving 2 doses of COVID-19 mRNA vaccines at a 3 month dose interval As of June 2021, although most LTCH workers had been fully vaccinated, significant sectors of the Canadian population had only been administered a single dose of a COVID-19 vaccine because the interval between dose 1 vs. dose 2 was extended as a dose sparing measure. Thus, we wished to ascertain if salivary Ab could be detected after a single dose of a COVID-19 mRNA vaccine, and how long these Ab would persist. We therefore collected samples from a second cohort of healthy adults that were followed longitudinally (Medical Sciences Building cohort-MSB; Supplementary Fig. 1 and Supplementary Table 4 ). These participants received 1 dose of BNT1162b2 and a second dose of BNT1162b2 3 months later, with samples taken at baseline, 2 weeks post-dose 1, 3 months post-dose 1, and 2 weeks post-dose 2. We observed that 97% and 93% of participants were positive for anti-Spike IgG and IgA, and 52% and 41% were positive for anti-RBD IgG and IgA Ab in their saliva 2 weeks postdose 1 (Fig. 1E-H) . Of note, 3 months after dose 1, the median level of salivary anti-Spike/RBD Ab had diminished nearly to baseline. Following administration of a second dose of mRNA, while antigen-specific IgG levels recovered upon administration of dose 2 as expected (Fig. 1E, F) , a second dose did not further augment antigen-specific salivary IgA levels in most subjects, and only approximately 30% of participants remained positive for IgA after dose 2 (Fig. 1G, H) . Therefore, the IgG and IgA response to COVID-19 mRNA vaccination differs upon administration of dose 2.

Detection of secretory component associated with anti-Spike and anti-RBD antibodies in saliva from participants receiving COVID-19 mRNA vaccines The i.m. route induces an immune response in the axillary draining lymph node that is biased towards class switch to IgG rather than IgA. Thus, we were curious as to why we were able to detect an IgA anti-Spike/RBD response in the saliva, particularly after one dose of mRNA vaccine. We hypothesized that some IgA may be produced locally. To test this, we designed an ELISA to detect any secretory component associated with anti-Spike/RBD Ab. We determined that the secretory component signal associated with anti-SARS-CoV-2 salivary Ab could be out-competed with recombinant secretory component, and no anti-Spike/RBD secretory component signal was detected in pre-pandemic colostrum, demonstrating assay specificity ( Supplementary Fig. 3 ). We then measured secretory component associated anti-SARS-CoV-2 Ab in the saliva of vaccinated LTCH participants who had received two doses of either BNT162b2 or mRNA-1273. We found that secretory component associated anti-Spike and anti-RBD Ab could be detected in 30% and 58% of participants, respectively, although the levels were significantly lower than what was observed in COVID-19 patients ( Fig. 2A, B) . The anti-SARS-CoV-2 associated secretory component signal was independent of prior SARS-CoV-2 exposure as we observed no significant difference in this signal comparing participants who were positive vs. negative for serum anti-nucleocapsid protein Ab (Fig. 2C, D) . Of note, if we divided the LTCH cohort into those participants who were positive versus negative for anti-Spike/RBD IgA, we observed that the secretory component signal was only detected in the IgA + participants (Fig. 2E, F) . Moreover, anti-Spike/RBD IgA and the level of secretory component positively correlated with each other ( Supplementary  Fig. 4) . Combined with the finding that most participants did not produce IgM Ab to anti-Spike/RBD (data not shown), we conclude that secretory component is associating with anti-Spike/RBD IgA (SIgA). Therefore, a local SIgA response to Spike/RBD is produced in response to mRNA vaccination in some participants.

Decay kinetics of anti-Spike/RBD IgG and IgA and neutralizing activity in the saliva of COVID-19 mRNA vaccinated participants Since the LTCH workers were the earliest recipients of mRNA vaccines during the Canadian roll-out, we had the opportunity to examine the stability of anti-Spike/RBD IgG and IgA in the saliva of a subset of the LTCH cohort who had reached 6 months post-dose 2. We examined the level of salivary Ab only in those participants who had remained anti-Spike IgA positive at 2 weeks post-dose 2 (n = 32,~30% of participants). Using a paired analysis, we observed a significant decline in antigen-specific IgG and IgA levels in the saliva at this time point compared to 2 weeks postdose 2 ( Fig. 3A-D) .

Although low, anti-Spike Ab associated with the secretory component remained stable in both BNT162b2 and mRNA-1273 vaccinated subjects, anti-RBD secretory component-associated antibodies were significantly reduced at 6 months post-dose 2 (Fig. 3E, F) .

In addition, we assessed if COVID-19 mRNA vaccination provokes SARS-CoV-2 neutralizing activity in the oral cavity, and if this would also be subject to decay. Specifically, saliva at twofold dilutions was added to hACE2-mCherry expressing HEK293 cells that were co-incubated with recombinant Vesicular Stomatitis Virus (rVSV)-eGFP-SARS-CoV-2-Spike. Infection of HEK293 cells was measured by fluorescence over the course of 72 h, with prevention of infection by added saliva read out as a green fluorescence reduction in neutralization titer (FRNT). We observed that although saliva from LTCH participants taken at 2-4 weeks post-dose 2 had variable capacity to prevent viral entry into hACE2 + HEK293 cells, this was significantly reduced at 6 months post-dose 2 ( Supplementary Fig. 5 ).

Taken together, with the exception of Spike-specfic SIgA, most Ab against SARS-CoV-2 in the saliva, as well as neutralizing capacity, significantly decline over a 6-month period.

Correlates of protection against SARS-CoV-2 breakthrough infection are ill-described and have only been examined for IgG 10 . Given that we have detected anti-Spike/RBD IgA in the serum and saliva of mRNA vaccinated participants, we examined whether IgA levels may be associated with protection against breakthrough infection. On 20th April 2021, an outbreak of P.1 lineage SARS-CoV-2 (Gamma variant) was declared in a Toronto LTCH home 11 . The residents, who were all doubly vaccinated with mRNA-1273, were part of a separate serum antibody study following LTCH residents that included sampling of blood, but not saliva, at 2-4 weeks post-dose 2 (received on February 2, 2021). Since anti-Spike and anti-RBD IgA levels correlate with each other at 2-4 weeks post-dose 2 ( Supplementary Fig. 2 ), we also measured serum anti-Spike IgA at this time point as a proxy of the salivary convalescent controls (n = 11) and individually run negative controls (n = 9). The positive cutoff (dotted line) was calculated as 2 standard deviations above the mean of a pool of negative control samples (n = 51) for each individual assay. E-H Anti-Spike/RBD IgG and IgA were also analyzed in the saliva of vaccinated participants that were followed at sequential timepoints before and after dose 1 and dose 2 of mRNA vaccination (MSB-1, n = 27, with n = 13 participants at 3 months post-dose 1, and n = 26 at 2 weeks post-dose 2). All data are expressed as a percentage of a pooled positive plate control of n = 11 COVID-19 convalescent saliva (see Methods). Solid black bars denote the median for each cohort. Mann-Whitney U test was used to calculate significance for (A-D). For (E-H) a Friedman pair-wise multiple comparison test was used for pair-wise comparisons for multiple timepoints, considering only data points from participants who had completed all timepoints (n = 11). *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. † = we were only able to recall n = 13 of the participants at 3 months post-dose 1.

IgA response. In the context of this isolated outbreak where n = 5 residents were infected we did not see a significant difference between exposed infected vs. exposed uninfected participants in terms of levels of anti-Spike and anti-RBD IgG at 2-4 weeks postdose 2 (Fig. 4A ). We noted a trend of reduced anti-Spike and anti-RBD IgA levels in exposed infected vs. exposed uninfected participants that reached significance but did not survive a post hoc multiple test correction (Fig. 4B ). Given the number of breakthrough infections in this isolated outbreak was small, we tested a replication cohort. We therefore measured serum antibody levels at 2-4 weeks post-dose 2 in a larger case-control cohort of healthcare workers (HCW) who had received 2 doses of BNT162b2 at the Sheba Medical Center in Ramat Gan, Israel, with the majority of cases infected with the Alpha SARS-CoV-2 variant 10 . We found that anti-Spike/RBD IgG levels were modestly but not significantly lower in the serum of cases versus controls (Fig. 4C ). In contrast, anti-Spike and anti-RBD IgA levels were both significantly lower in cases vs. controls, and this held up to a multiple comparison test (Fig. 4D ). Using the more well-powered Sheba cohort, if we convert these values to the BAU/ml using the 20/136 WHO serum standard (Supplementary Table 7 ), the median level of IgA that is associated with breakthrough cases is 152.78 BAU/ml and 162.31 BAU/ml for anti-Spike and anti-RBD respectively compared to uninfected controls whose median IgA levels were 471.36 BAU/ml and 495.68 BAU/ml for anti-Spike and anti-RBD respectively.

Thus, despite different vaccine regimes (mRNA-1273 vs. BNT162b2), geographical locations (Canada vs. Israel) and viral exposures (Gamma vs. Alpha), in both cohorts anti-Spike/RBD serum IgA but not IgG levels are lower at 2-4 weeks post-dose 2 in participants who subsequently are infected with SARS-CoV-2.

In this study we observed robust IgG levels to Spike/RBD in the saliva of participants immunized with BNT162b2 or mRNA-1273 that correlate with the systemic IgG response, and these IgG Ab are significantly diminished at 6 months post-dose 2. A different picture emerges however, for salivary anti-Spike/RBD IgA. In nearly all participants, we observed a modest SIgA response to Spike/RBD in the saliva after a single dose of mRNA vaccine which is maintained in~30% of participants after dose 2. Moreover, in these subjects, although the SIgA response to Spike is lower than what is observed with SARS-CoV-2 infection, it is more resilient to decay. This preserved SIgA response observed in a minority of vaccinated participants may be very important for preventing breakthrough infections. Indeed, vaccinated participants who subsequently experienced a SARS-CoV-2 infection had We show evidence of local SIgA production in the saliva by measuring secretory component levels associated with antigenspecific Ab. This is consistent with recent studies measuring anti-Spike IgA in breastmilk of BNT162b2 and mRNA-1273 vaccinated lactating mothers 16 . It is unclear how anti-Spike/RBD SIgA is generated in the saliva following i.m. immunization. Of note, Spike protein can be detected in the plasma, increasing one to 5 days after mRNA-1273 vaccination using an ultra-sensitive detection technique 17 . Thus, one possibility is that plasma-associated Spike antigen may reach the salivary glands (which are surrounded by capillaries 18 ,) provoking a local SIgA response. Another possibility is that a mucosal IgA response to mRNA vaccination takes place in the gut, and plasma cells generated at this location leave the gut (as we have shown before 19 ,) disseminating to other mucosal surfaces such as the oral cavity. Indeed, expression of the gut homing integrin (α4β7 on circulating immune cells has been observed following administration of yellow fever 20 and cholera toxin vaccines 21 via the systemic route in humans. Recirculation between mucosal compartments could explain the sustained levels of anti-Spike SIgA in the saliva at 6 months post-boost which is not observed for the systemic IgG response. Animal models are ideal for obtaining further insights into mechanisms that explain how SARS-CoV-2 mRNA vaccines provoke antigenspecific SIgA production at mucosal surfaces.

Dimeric secretory IgA has been shown to have potent neutralization activity against SARS-CoV-2 in vitro 6 . Using a flow cytometry-based pseudovirus entry assay we found that saliva from COVID-19 mRNA vaccinated subjects can variably prevent pseudovirus infection of ACE2 + cells. This is consistent with work by Nahass et al. who have examined neutralizing capacity of saliva in response to different vaccine regimes 16 . It is unclear whether The Wilcoxon signed-rank test was used to calculate significance between groups. ns not significant, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

IgG or IgA or both are the main contributors to the neutralizing capacity of saliva from vaccinated subjects. Indeed, although the neutralization activity of saliva at 6 months post-dose 2 was greatly diminished, anti-Spike SIgA levels had not significantly decayed at this time-point. It will be of interest to examine structural interactions between SIgA dimers with SARS-CoV-2 in the context of the mucosa. There are some limitations to our study: we did not have saliva available from the LTCH residents nor the Sheba HCW for confirmation that the association of high serum anti-Spike IgA levels with protection against breakthrough infection was also observed in a mucosal secretion. We also recognize potential age and health-related differences between these cohorts, as Abe et al. reported LTCH residents had lower nAb and IgG titers in serum compared to LTCH staff 22 . However, our observation that anti-Spike/RBD IgA in the serum and saliva at 2 weeks post-dose 2 correlate ( Supplementary Fig. 2) , suggests that what we are detecting in the serum at least partially mirrors the oral compartment. Furthermore, we are in the process of collecting saliva samples from the MSB and LTCH cohorts (Fig. 1) that have Serum samples from vaccinated Toronto LTCH residents were taken at 2-4 weeks post-dose 2. Serum anti-Spike/RBD IgG (A) and IgA (B) levels were compared in participants who were subsequently exposed to P.1 gamma SARS-CoV-2 and either infected (exposed infected, n = 5) or not infected (exposed uninfected, n = 12). In a separate cohort of double vaccinated healthcare workers from the Sheba Medical Center in Ramat Gan, Israel, serum samples were taken at 2-4 weeks post-dose 2. Serum anti-Spike/RBD IgG (C) and IgA (D) levels were compared in participants who experienced a breakthrough infection (cases, n = 11) vs. controls who did not (control, n = 56). Solid black lines denote the median for each cohort. Kruskal-Wallis test with a correction for multiple comparisons was used to conduct statistical analysis between groups. NS not significant, *p < 0.5; **p < 0.01. experienced Omicron breakthroughs, providing a potential means for assessing correlates of protection against breakthrough infection in the oral compartment.

In summary, we provide evidence that anti-Spike SIgA are induced and maintained in the saliva of approximately 30% of mRNA vaccinated participants and that high anti-Spike/RBD serum IgA levels are associated with protection against subsequent breakthrough infection. Intranasal (i.n.) adenoviral-based vaccines that robustly induce anti-Spike/RBD IgA have been shown to prevent transmission in golden hamsters 23 . Moreover, a SARS-CoV-2 mRNA vaccine i.m. prime followed by an adenoviral i.n. boost provides superior protection against infection compared to i.m. only regimes, even against divergent variants such as B.1.1.351 (beta) 24 . Our data suggest that an i.m. prime followed by i.n. boost strategy is worth considering for preventing personto-person transmission of SARS-CoV-2 and for broad protection against emerging variants.

The Mount Sinai Hospital Research Ethics Board (REB) granted approval for recruiting staff in long-term care facilities located in the Greater Toronto Area for blood and saliva collection and for conducting serum ELISAs at the Lunenfeld-Tanenbaum Research Institute (study number: 20-0339-E). The University of Toronto REB granted approval for participant recruitment to collect blood and saliva samples and for conducting saliva ELISAs (study number: 23901). The University of Saskatchewan REB granted approval for saliva sample collection during the pre-COVID era (study number: BIO-USask-1579). The Sheba Medical Center REB granted approval for healthcare worker recruitment (study number: 8008-20-SMC). Table 1 ) Staff working in Long-Term Care Homes (LTCHs) in Ontario were eligible to participate in the study if: (1) they were over the age of 18, (2) they were comfortable (a) reading and writing in English, (b) providing blood samples, and (c) sharing their COVID-19 diagnostic status with the study team. To invite LTCH staff to participate in the study, study staff approached the administrators and/or directors of LTCHs to assess if they were interested in being a participating site in the study. If they were interested, the administrators and/or directors shared information about the study with their staff and provided a deadline by which staff should opt-out if they did not want to be contacted about the study by a member of the study team. The LTCH administration then shared with the study team the contact information of all LTCH staff who did not opt-out of further communication about the study. The study team then contacted these staff by phone to invite them to participate in the study. This active recruitment strategy was paired with passive recruitment strategies, including having the participating LTCHs and/or relevant staff organizations share recruitment advertisements about the study, as well as having participants approach study staff when they were on-site for sample collection. A total of 12 LTCHs participated as a site in this study. Across these 12 sites, 647 individuals were invited to participate. Of these 647, a subset of individuals was not able to be reached by phone (n = 242, 37%), were not eligible to participate (n = 105, 16%), or refused to participate (n = 131, 20%). Common reasons for participant ineligibility included being beyond 6 weeks after their 2nd dose of their COVID-19 vaccine or not being interested in receiving a COVID-19 vaccine. The most common reason for refusal to participate was lack of interest in being part of research. Finally, a subset of individuals (n = 13) withdrew prior to their first sample collection, resulting in a final sample size of 156 participants.

Recruitment and participants-COVID-19 acute and convalescent participants (Supplementary Table 2b ) Acute and convalescent serum and saliva samples were obtained from patients identified by surveillance of COVID-19 (confirmed by PCR; in-and out-patients) by the Toronto Invasive Bacterial Diseases Network in metropolitan Toronto and the regional municipality of Peel in south-central Ontario, Canada (REB studies #20-044 Unity Health Network, #02-0118-U/05-0016-C, Mount Sinai Hospital). Consecutive consenting patients admitted to four Toronto Invasive Bacterial Disease Network (TIBDN) hospitals were enrolled: these patients had serum and saliva collected at hospital admission, and survivors were asked to submit repeat samples at 4-12 weeks post-symptom onset (PSO). Consecutive out-patients diagnosed at the same 4 hospitals prior to March 15th and on a convenience sample of later days were approached for consent to collect serum and saliva at 4-12 weeks PSO. Patients were interviewed and patient charts reviewed to determine age, sex, symptom onset date, and disease severity, as described in 3 . Table 4 Table 5 ) Residents in this LTCH were offered mRNA-1273 on 4 January and 2 February 2021. The facility has resident rooms on the second through fifth floors and administrative/service areas in the basement and ground floors. By 28 March, 100 of 124 residents (81%) and 120 of 224 staff (54%) had received 2 vaccine doses, and 7 residents (6%) and 29 staff (13%) had received 1 dose. Ongoing pandemic measures included no visitors except essential caregivers, symptom surveillance twice daily for residents and staff/caregivers, staff/caregiver nasopharyngeal swab SARS-CoV-2 testing 2-3 times weekly by reverse-transcriptase polymerase chain reaction (rt-PCR) or daily rapid antigen testing (Panbio), cohorting of staff by floor to limit any exposure to a single floor, and universal masking and eye protection for staff/caregivers. All positive rapid antigen tests were confirmed by SARS-CoV-2 PCR. All specimens PCR positive for SARS-CoV-2 were tested for N501Y and E484K mutations using a multiplex real-time PCR assay. Public Health Ontario Laboratory and/or Sunnybrook Health Sciences Center performed whole-genome sequencing and assessed predicted Pango lineage, as described elsewhere 25 . On the 2 outbreak units, positive specimens with concentrations too low for mutation detection were assumed to be lineage P.1. Table 6 ) Beginning January 20, 2021, 11 days after the first staff members had received a second dose of the BNT162b2 vaccine, a study was initiated to identify every breakthrough infection, including asymptomatic infections. Serum analysed herein was collected at~2-4 weeks post-dose 2. Data were collected until April 28, 2021. Efforts were extended to identify new cases with the use of daily health questionnaires, a telephone hotline, extensive epidemiologic investigations of exposure events, and contact tracing of infected patients and personnel. Testing for the presence of SARS-CoV-2 by RT-PCR was performed on fully vaccinated staff members who were symptomatic or had been exposed to an infected person, regardless of symptoms. Antigendetecting rapid diagnostic testing (Ag-RDT) was available as an initial screening tool in the personnel clinic in combination with RT-PCR testing. A breakthrough infection was defined as the detection of SARS-CoV-2 on RT-PCR assay performed 11 or more days after receipt of a second dose of BNT162b2 if no explicit exposure or symptoms had been reported during the first 6 days. A matched case-control analysis was used to identify possible correlates of breakthrough infection. For the case-control analysis, control serum samples were selected from a prospective cohort study to analyze vaccine-induced immune responses and dynamics at the Sheba Medical Center. Each breakthrough case was matched with control samples that had been obtained from uninfected controls according to the following variables: sex, age, immunosuppression status, and the interval between the second dose of BNT162b2 vaccine and serology at the peri-infection timepoint (see ref. 10 ) We compared antibody titers of cases versus controls obtained during the initial post-vaccination period. Breakthrough cases for which serologic samples were not available were excluded from this analysis.

Saliva collection LTCH, MSB-1 and MSB-2 participants were told not to eat, drink or smoke at least 30 min prior to collection. Subsequently, saliva were collected using Salivette ® tubes (Sarstedt, Numbrecht, Germany), a collection system which consists of a cotton ball which participants chew for exactly 3 min and place into a tube, which is then placed into a larger outer tube. The entire system is spun in a centrifuge at 1000 × g for 5 min at room temperature. The inner tube contains a hole at the bottom, which allows all the saliva absorbed by the cotton ball to filter into the larger outer tube. The total saliva volume from each participant was then separated into 300-500 µl aliquots and stored at −80°C until the time of testing. Given that these samples were collected from vaccinated participants who reported no symptoms of COVID-19 infection, we did not conduct any measures for viral inactivation.

Antigen production-Saliva assay The expression, purification and biotinylation of the SARS-CoV-2 RBD and Spike ectodomain were performed as recently described 3 . The human codon optimized cDNA of the SARS-CoV-2 Spike protein was purchased from Genscript (MC_0101081). The soluble RBD (residues 328-528, RFPN…CGPK) was expressed as a fusion protein containing a C-terminal 6 × His tag followed by an AviTag. The soluble trimeric Spike protein ectodomain (residues 1-1211, MFVF…QYIK) was expressed with a C-terminal phage foldon trimerization motif followed by a 6 × His tag and an AviTag. To help stabilize the Spike trimer in its prefusion conformation, residues 682-685 (RRAR) were mutated to SSAS to remove the furin cleavage site and residues 986 and 987 (KV) were each mutated to a proline residue (51). Stably transfected FreeStyle 293-F cells secreting the RBD and soluble Spike trimer were generated using a previously reported piggyBac transposon-based mammalian cell expression system (52). Protein production was scaled up in 1 L shake flasks containing 300 mL FreeStyle 293 medium. At a cell density of 10 6 cells/mL, 1 μg/mL doxycycline and 1 μg/mL Aprotinin were added. Every other day 150 mL of medium was removed and replaced by fresh medium. The collected medium was centrifuged at 10000 × g to remove the cells and debris and the His-tagged proteins were purified by Ni-NTA chromatography. The eluted protein was stored in PBS containing 300 mM imidazole, 0.1% (v/v) protease inhibitor cocktail (Sigma, P-8849) and 40% glycerol at −12°C. Shortly before use, the RBD and Spike proteins were further purified by size-exclusion chromatography on a Superdex 200 Increase (GE healthcare) or Superose 6 Increase (GE healthcare) column, respectively. Purity was confirmed by SDS-PAGE. For the Spike protein, negative stain electron microscopy was used show evidence of high-quality trimers. The Avi-tagged proteins, at a concentration of 100 μM or less, were biotinylated in reaction mixtures containing 200 μM biotin, 500 μM ATP, 500 μM MgCl2, 30 μg/mL BirA, 0.1% (v/v) protease inhibitor cocktail. The mixture was incubated at 30°C for 2 h followed by size-exclusion chromatography to remove unreacted biotin.

Enzyme-linked immunosorbent assays for detecting anti-Spike and anti-RBD IgA and IgG in saliva 96-well plates pre-coated with streptavidin were used for all saliva assays. We have previously determined that coating plates with 50 μl per well of 2 μg/ml of biotinylated RBD or 20 μg/ml of biotinylated Spike diluted in sterile phosphate-buffered saline (PBS) was the ideal coating solution. Control wells were coated with 50ul per well of sterile PBS. After coating with the antigen of interest and incubating overnight at 4°C, the coating solution was discarded and plates were blocked with 5% BLOTTO solution (5% w/vol skim milk powder (BioShop, CAT# SKI400.500)). Plates were incubated with the blocking solution at 37°C for 2 h, and the solution was discarded immediately prior to adding samples to each well. During the blocking incubation, frozen saliva samples were removed from −80°C storage, thawed and diluted using 2.5% BLOTTO at a range of 1:5-1:20. Sample dilutions were pre-incubated in a separate streptavidin-coated plate with no antigen to reduce anti-streptavidin activity in the saliva.

Samples were incubated in the pre-adsorption plate for 30 min at 37°C, after which 50 μl of each sample from the pre-adsorption plate was transferred to the corresponding wells of the antigencoated plates and incubated for 2 h at 37°C. Next, samples were discarded from the antigen-coated plates, and the plates were washed 3× with PBS + 0.05% Tween 20 (PBS-T (BioShop, CAT# TWN510)). 50 μl of Horse radish peroxidase (HRP)-conjugated goat anti-human-IgG and IgA secondary antibodies (Southern Biotech, IgG: 2044-05, IgA: 2053-05) were added to the appropriate wells at dilutions of 1:1000 and 1:2000 in 2.5% BLOTTO, respectively, and incubated for 1 h at 37°C. Plate development was done by adding 50 μL of 3,3′,5,5′tetramethylbenzidine (TMB) Substrate Solution (ThermoFisher, 00-4021-56) to each well. The reaction was then stopped by adding 50 µl/well of 1 N H2SO4, and optical density (OD) was read at a wavelength of 450 nm (OD450) on a spectrophotometer (Thermo Multiskan FC). For each sample, the raw OD450 for the PBS control well was subtracted from the raw antigen-specific OD450 value for each sample, at each sample dilution (1:5, 1:10, 1:20). The adjusted OD450 value for each sample dilution (1:5, 1:10, 1:20) was used to calculate the area under the curve (AUC) for each individual sample. The sample AUC was then normalized to the AUC of the positive control, which consisted of saliva collected from COVID-19 acute and convalescent participants. The normalized AUC was multiplied by 100 to give a final percentage, which we deemed the "% of positive control". Each plate also included 1-3 negative controls (with a total of n = 78 negative controls used between all experiments presented), which consisted of pre-COVID era saliva incubated in antigen-coated wells. Integrated scores were calculated for all negative control samples, using the same calculation method used for cohort 1 and 2 samples. "Positive" cutoff values for each antigen-specific isotype were calculated using the following formula: average integrated scores of negative samples + 2 (standard deviation of negative control integrated scores). Samples whose score was above the resulting cutoff for each antigen-specific isotype was used to determine which samples from cohorts 1 and 2 had detectable antibody levels in their saliva.

ELISA for detection of secretory component-associated anti-Spike/RBD antibodies Secretory component-associated antibodies were detected by modifying our saliva Spike/RBD ELISA by using an HRP-conjugated Goat anti-human secretory component detection reagent at a dilution of 1:750 from Nordic MUBio (Cat# GAHu/SC/PO).

Enzyme-linked immunosorbent assays for detecting anti-Spike and anti-RBD IgA and IgG in serum An automated chemiluminescent ELISA assay was used to analyze the levels of IgG and IgA antibodies to the Spike trimer, and its RBD with the following modifications: RBD (PRO1151, 20 ng/well) antigens were produced in CHO (Chinese Hamster Ovary) cells, and were a kind gift from Dr. Yves Durocher, National Research Council of Canada (NRC). The secondary antibody for IgG was an IgG-HRP fusion (PRO1146, 1:6700 or 0.9 ng/well), donated by the NRC. A standard curve of the VHH72 monoclonal antibody 26 fused to a human IgG1 Fc domain (PRO23, also from the NRC) was generated for calibrating the anti-Spike and anti-RBD IgG response. All other antigens, detection reagents and calibration reagents were as previously described 27 . The data analysis also proceeded as in 27 , with the following exceptions: Blanks were not subtracted from the chemiluminescence raw values of the samples, and the raw values were normalized to a blanksubtracted point in the linear range of the calibration standard curve (for Spike and RBD, the reference point was 0.0156 µg/ml and for N either 0.0625 µg/ml or a combination of 0.0156 and 0.0312 µg/ml was used to maintain consistent values across tests that used different anti-N antibody batches). The results are represented as a "relative ratio" to this reference point. To define the cutoff for positive antibody calls for each antigen for IgG when 0.0625 µl/well of sample was added, 3 standard deviations from the mean of the log negative control distribution from >20 different runs collected over 4 months was used. For IgA, negatives from 2 different runs over one month were used. In all cases, the selected cut offs correspond to <2% False Positive Rate (FPR) assessment, based on Receiver Operating Characteristic Curves.

Calibration to the WHO standard for serum IgA detection For anti-Spike and anti-RBD IgA, the relative ratio from the WHO International Standard (IS, National Institute for Biological Standards and Control (NIBSC, South Mimms, United Kingdom), Code 20/136, pooled convalescent plasma) at different sample dilutions (in binding antibody units (BAU)/ml) was represented in log-log scale. The response curve was modeled by an S-shaped sigmoid curve (y = a*x/(1 + b*x) + c, expressed in log(y), log(x) coordinates). The nls function (nonlinear least squares) from the R stats package was used to best match the response to the measured data. The interval where the log response was considered linear to the log of the BAU/ml was selected visually. The lm (linear model) function from the R stats package was used to obtain the parameters of the linear approximation of the curve. Estimated BAU/ml values were obtained from the above linear approximation using only sample dilutions that were within this linear range, and further adjusted by accounting for the dilution factor. To convert from relative ratios (RRs) to BAU/m1 for plasma or serum samples where only the reference curve was included in the same test, a conversion formula can be applied: log2 (sample BAU/ml at dilution fold d) = (log2 (sample RR) − a)/b + log2 (d) where a and b represent the y-intercept and slope of the linear interval of the IS curve. Please see Supplementary Table 7 for  conversion table and ref. 27 .

Flow cytometry method for detection of neutralizing activity Neutralizing activity was measured at two step saliva dilutions (1:12-1:384) following incubation with recombinant Vesicular Stomatitis Virus (rVSV)-eGFP-SARS-CoV-2-Spike in which the VSV-G protein was replaced with SARS-CoV-2-Spike protein was propagated on MA104 cells 28 15-140-163) . After visible cytopathic affect, supernatant was filtered, aliquoted and stored at −80°C. Supernatant was added to Human Embryonic Kidney (HEK) 293 cells were engineered to encode human angiotensin converting enzyme 2 (hACE2 in the pDEST-mCherry vector) as previously described 29 . HEK293-hACE2-mCherry cells were cultured in Gibco Dulbecco's Modified Eagle Medium (DMEM) formulation containing glucose, L-glutamine and sodium pyruvate (Gibco, Cat. No. 11995065) with 1× Penicillin/Streptomycin. Geneticin Selective Antibiotic (G418) (Gibco, Cat. No. 10131035) was added at a concentration of 500 µg/mL to maintain hACE2-mCherry expression. Cells were grown in 5% carbon dioxide (CO 2 ) at 37°C and passaged every 3 days using Versene solution (Gibco, Cat. No. 15040066). HEK293-hACE2-mCherry cells were seeded at a density of 25,000 cells per well in a 96-well, flat-bottom tissue culture coated plate. Outer rows were avoided to reduce assay variations resulting from edge effect in the IncuCyte. In a separate 96-well plate, samples were serially diluted and incubated with 50 μl of rVSV-eGFP-SARS-CoV-2-S for 2 h at 37°C in 5% CO 2 . Each sample plate included a dilution of anti-RBD antibody (Invitrogen, Cat. No. 703958) of 10, 5, 1, 0.5, 0.1, and 0.05 μg/mL as a positive control. After incubation, the mixture of sample and rVSV-eGFP-SARS-CoV-2-S was transferred to the plated HEK293-hACE2-mCherry cells at a 1:1 ratio of culture media to virus/sample suspension. Plates loaded in the IncuCyte were imaged every 1-4 h for a total of 72 h with 4 scans per sample well to visualize neutralization. The total integrated intensity of the fluorescent value of the lowest anti-RBD condition (0.05 μg/mL) with the rVSV-eGFP-SARS-CoV-2-S controls at 12 h intervals was used to normalize separate experiment runs. Each plate was normalized either to the mean of the rVSV-eGFP-SARS-CoV-2-S supernatant controls or to the 0.05 μg/mL anti-RBD antibody. Normalization to 0.05 μg/mL anti-RBD was performed only if division by the triplicate rVSV-eGFP-SARS-CoV-2 control conditions resulted in loss of a sigmoidal shape of the anti-RBD curve. To quantitatively determine assay sensitivity, fluorescent reduction of neutralization titers (FRNT) 50 and 70 were calculated to determine the amount of monoclonal anti-RBD needed to prevent 50% and 30%, respectively, of the maximum infection ( Supplementary Fig. 10B , C) and plotted in log10 scale for plasma and log2 scale for saliva due to variations in antibody titers in different tissue types. This assay is sufficiently sensitive down to 5 μg/mL of neutralizing antibodies.

Antibody levels from the cross-sectional LTCH cohort were compared between groups using Mann-Whitney U tests. For correlation analyses, a nonparametric Spearman correlation test was applied. A Friedman pair-wise multiple comparison test was used to compare participants who were assessed at multiple timepoints ( Fig. 1) , whereas the Wilcoxon pair-wise signed-rank test was used when examining 2 different timepoints (Fig. 3 , Supplementary Fig. 5 ). Kruskal-Wallis with Dunn's multiple comparison test was conducted to examine the differences in antibody levels to RBD and Spike between infected vs. non-infected participants in the LTCH and Sheba cohorts. Adjusted P value based on 2 comparisons per family (RBD case vs. control and Spike case vs. control) were used to determine statistical significance. Analyses were performed in Prism (GraphPad) Version 9.2.0.

The relationships between clinical predictors (age, sex, SARS-CoV-2 infection prior to vaccination, and time from vaccination to sample collection) and antibody levels were examined in bivariate and a priori multivariable linear regression models. For saliva, four multivariable linear regression models were constructed to examine potential independent associations between the four clinical predictorsand anti-RBD IgA/IgGand anti-SpikeIgA/IgG. Kaplan-Meier survival curve were plotted by dividing each group based on serum titer at above or below median. Log-rank (Mantel-Cox) test was done to compare survival curve of each isotype (IgG and IgA anti-Spike) above and below median cohort values.

Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein

SARS-CoV-2 viral load in upper respiratory specimens of infected patients

Persistence of serum and saliva antibody responses to SARS-CoV-2 spike antigens in COVID-19 patients

SARS-CoV-2 infection of the oral cavity and saliva

The structure and dynamics of secretory component and its interactions with polymeric immunoglobulins

Enhanced SARS-CoV-2 neutralization by dimeric IgA

IgA dominates the early neutralizing antibody response to SARS-CoV-2

Sex differences in vaccine-induced humoral immunity

COVID-19 mRNA vaccination generates greater IgG levels in women compared to men

Covid-19 breakthrough infections in vaccinated health care workers

COVID-19 outbreak associated with a SARS-CoV-2 P.1 Lineage in a long-term care home after implementation of a vaccination program-Ontario

Antibody responses to SARS-CoV-2 mRNA vaccines are detectable in saliva

Efficient mucosal antibody response to SARS-CoV-2 vaccination is induced in previously infected individuals

Mucosal immune response in BNT162b2 COVID-19 vaccine recipients

Cutting edge: serum but not mucosal antibody responses are associated with pre-existing SARS-CoV-2 spike cross-reactive CD4(+) T cells following BNT162b2 vaccination in the elderly

Intramuscular SARS-CoV-2 vaccines elicit varying degrees of plasma and salivary antibody responses as

Circulating SARS-CoV-2 vaccine antigen detected in the plasma of mRNA-1273 vaccine recipients

Saliva in the diagnosis of diseases

Recirculating Intestinal IgA-producing cells regulate neuroinflammation via IL-10

Dynamic T cell migration program provides resident memory within intestinal epithelium

Differential expression of tissue-specific adhesion molecules on human circulating antibody-forming cells after systemic, enteric, and nasal immunizations. A molecular basis for the compartmentalization of effector B cell responses

Neutralizing antibody responses to SARS-CoV-2 variants in vaccinated Ontario long-term care home residents and workers

Oral and intranasal Ad5 SARS-CoV-2 vaccines decrease disease and viral transmission in a golden hamster model

Protective mucosal immunity against SARS-CoV-2 after heterologous systemic prime-mucosal boost immunization

Detection of the novel SARS-CoV-2 European lineage B.1.177 in Ontario

Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation

A scalable serology solution for profiling humoral immune responses to SARS-CoV-2 infection and vaccination

Neutralizing antibody and soluble ACE2 inhibition of a replicationcompetent VSV-SARS-CoV-2 and a clinical isolate of SARS-CoV-2

Cholesterol 25-hydroxylase suppresses SARS-CoV-2 replication by blocking membrane fusion

performed saliva Ab measurements, prepared figures and were supervised by J.G. G.C. coordinated intake, database entry and maintenance for all samples and also consented and performed phlebotomy on MSB participants and was supervised by

performed serum antibody measurements and were supervised by A.C.G. A.J. performed multivariate analyses and was supervised by A.M. Z.L. purified and prepared biotinylated Spike and RBD proteins and was supervised by J

recruited, consented and collected saliva from COVID-19 convalescent patients and were supervised by M.O. WS recruited, consented and collected saliva from pre-COVID-19 participants. O.L.R. assisted in data analysis

Dr. Yves Durocher at the National Research Council of Canada (NRC) kindly donated reagents for the serum assays. We thank Elizabeth Yue, Antonio Estacio, Serena Loklam Chau, Ryan Law, and Eric Yixiao Cao for LTCH sample processing and sustaining the infrastructure of sample processing in the Ostrowski lab. We thank Jeff Browning for critical reading of this paper. We would like to thank Florian Krammer for advice on secretory component detection by ELISA.

The authors declare no competing interests.

Supplementary information The online version contains supplementary material available at https://doi.org/10.1038/s41385-022-00511-0.Correspondence and requests for materials should be addressed to Sharon Straus, Anne-Claude Gingras or Jennifer L. Gommerman.Reprints and permission information is available at http://www.nature.com/ reprintsPublisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.