SUMMARY
While the standard regimen of the BNT162b2 mRNA vaccine includes two doses administered three weeks apart, some public health authorities decided to space them, raising concerns about vaccine efficacy. Here, we analyzed longitudinal humoral responses including antibody binding, Fc-mediated effector functions and neutralizing activity against the D614G strain but also variants of concern and SARS-CoV-1 in a cohort of SARS-CoV-2 naïve and previously infected individuals, with an interval of sixteen weeks between the two doses. While the administration of a second dose to previously infected individuals did not significantly improve humoral responses, we observed a significant increase of humoral responses in naïve individuals after the 16-weeks delayed second shot, achieving similar levels as in previously infected individuals. Our results highlight strong vaccine-elicited humoral responses with an extended interval BNT162b2 vaccination for naïve individuals.
INTRODUCTION
Since the end of 2019, the etiological agent of the Coronavirus disease 2019 (COVID-19), the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has spread worldwide causing the current pandemic (Dong et al., 2020; World Health Organization). In the last months, several vaccines against SARS-CoV-2 have been approved in many countries, including the Pfizer/BioNtech BNT162b2 mRNA vaccine. This vaccine targets the highly immunogenic trimeric Spike (S) glycoprotein that facilitates SARS-CoV-2 entry into host cells via its receptor-binding domain (RBD) that interacts with angiotensin-converting enzyme 2 (ACE-2) (Hoffmann et al., 2020; Walls et al., 2020) and has shown an important vaccine efficacy (Polack et al., 2020; Skowronski and De Serres, 2021).
The approved BNT162b2 mRNA vaccine regimen comprises two doses administered 3-4 weeks apart (WHO, 2021). However, at the beginning of the vaccination campaign (Winter/Spring 2021) vaccine scarcity prompted some public health agencies to extend the interval between doses in order to maximize the number of immunized individuals. This strategy was supported by results indicating that a single dose affords ∼90% protection starting two weeks post vaccination, concomitant with the detection of some vaccine-elicited immune responses (Baden et al., 2021; Pilishvili, 2021; Polack et al., 2020; Skowronski and De Serres, 2021; Tauzin et al., 2021).
The rapid emergence of several variants of concerns (VOCs) and variants of interest (VOIs), which are more transmissible and in some cases more virulent (Allen et al., 2021; Brown et al., 2021; Davies et al., 2021; Fisman and Tuite, 2021; Pearson et al., 2021) remains a major public health preoccupation as the vaccine campaign advances worldwide. For example, the mutation D614G in the S glycoprotein which appeared very early in the pandemic is now present in almost all circulating strains (Isabel et al., 2020). The B.1.1.7 (Alpha) variant emerged in late 2020 in the United Kingdom and due to its increased affinity for the ACE2 receptor that leads to increased transmissibility (Davies et al., 2021), it became in just a few months a predominant strain worldwide (Davies et al., 2021; Prévost et al., 2021; Rambaut et al., 2020). The B.1.351 (Beta) and P.1 (Gamma) variants that first emerged in South Africa and Brazil respectively have largely spread and are now circulating in many countries (ECDC, 2021; Tang et al., 2021). The B.1.526 (Iota) variant first identified in New York in early 2021 (Annavajhala et al., 2021a) is in an upward trajectory in the United States (Annavajhala et al., 2021b). More recently, the B.1.617.2 (Delta) variant which emerged in India and has a high transmissibility is now the dominant strain in several countries (Allen et al., 2021; Dagpunar, 2021). Although several studies have shown that mRNA vaccines protect against severe disease caused by these variants, it has also been shown that some of them present resistance to some vaccine-elicited immune responses, notably against neutralizing antibodies (Annavajhala et al., 2021a; Goel et al., 2021a; Planas et al., 2021a; Puranik et al., 2021; Wall et al., 2021; Wang et al., 2021a). Most of these studies were based on the analysis of plasma samples collected from vaccinees following a short (3-4 weeks) interval between doses. Little is known about vaccine-elicited immune responses with longer dose intervals. Here, we characterized vaccine-elicited humoral responses in a cohort of SARS-CoV-2 naïve and previously infected individuals that received the two doses with an extended interval of sixteen weeks.
RESULTS
We analyzed the longitudinal humoral responses after vaccination with the BNT162b2 mRNA vaccine in blood samples, with an interval of around 16 weeks between the two doses (median [range]: 111 days [90–134 days]). The cohort included 22 SARS-CoV-2 naïve and 21 previously infected (PI) donors tested SARS-CoV-2 positive by nasopharyngeal swab PCR around 9 months before their first dose (median [range]: 280 days [116-326 days]). In the cohort of PI individuals, 10 donors did not receive the second injection, leaving 11 PI donors with two doses. The blood samples were collected at different time points: prior the first dose of vaccine (V0), three weeks (V1, median [range]: 21 days [13–28 days]) and three months (V2, median [range]: 85 days [73–104 days]) after the first dose of vaccine and three weeks after the second vaccine injection (V3, median [range]: 21 days [13–42 days]). Data collected at V0 and V1 have been previously described (Tauzin et al., 2021). Basic demographic characteristics of the cohorts and detailed vaccination timepoints are summarized in Table 1 and Figure 1A.
(A) SARS-CoV-2 vaccine cohort design. (B-E) Indirect ELISA was performed by incubating plasma samples from naïve and PI donors collected at V0, V1, V2 and V3 with recombinant SARS-CoV-2 RBD protein. Anti-RBD Ab binding was detected using HRP-conjugated (B) anti- human IgM+IgG+IgA (C) anti-human IgM, (D) anti-human IgG, or (E) anti-human IgA. Relative light unit (RLU) values obtained with BSA (negative control) were subtracted and further normalized to the signal obtained with the anti-RBD CR3022 mAb present in each plate. (F-I) Cell- based ELISA was performed by incubating plasma samples from naïve and PI donors collected at V0, V1, V2 and V3 with HOS cells expressing full-length SARS-CoV-2 S. Anti-S Ab binding was detected using HRP-conjugated (F) anti-human IgM+IgG+IgA (G) anti-human IgM, (H) anti- human IgG, or (I) anti-human IgA. RLU values obtained with parental HOS (negative control) were subtracted and further normalized to the signal obtained with the CR3022 mAb present in each plate. Naïve and PI donors with a long interval between the two doses are represented by red and black points respectively and PI donors who received just one dose by blue points. (Left panels) Each curve represents the normalized RLUs obtained with the plasma of one donor at every time point. Mean of each group is represented by a bold line. The time of vaccine dose injections is indicated by black triangles. (Right panels) Plasma samples were grouped in different time points (V0, V1, V2 and V3). Undetectable measures are represented as white symbols, and limits of detection are plotted. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
Elicitation of SARS-CoV-2 antibodies against the full Spike and its receptor-binding domain
To evaluate vaccine responses in SARS-CoV-2 naïve and PI individuals, we first measured the presence of SARS-CoV-2-specific antibodies (Abs) (IgG, IgM, IgA) recognizing the receptor-binding domain (Figure 1B-E) using an ELISA RBD assay or the native full-length S glycoprotein expressed at the cell surface (Figure 1F-I) using a cell-based ELISA assay. Both assays have been previously described (Anand et al., 2021; Beaudoin-Bussières et al., 2020; Prévost et al., 2020). Prior to vaccination (V0), no SARS-CoV-2 specific Abs were detectable in SARS-CoV-2 naïve individuals, except for anti-Spike IgM (31.8% seropositivity) which are likely to be cross-reactive antibodies against the S2 subunit (Fraley et al., 2021; Hicks et al., 2020; Ng et al., 2020). SARS-CoV-2 PI individuals still had detectable Abs several months post-symptoms onset, especially IgG, in agreement with previous observations (Anand et al., 2021; Dan et al., 2021; Tauzin et al., 2021; Wang et al., 2021b). For both groups, the first dose of vaccine induced a significant increase of total immunoglobulins (Igs) recognizing the RBD or the Spike protein three weeks post-vaccine (V1), with a significantly higher response for the PI group (Figure 1B- I). At V2 (i.e., 12 weeks post vaccination), while anti-Spike total Ig levels remained stable, we observed a decrease in anti-RBD total Ig levels in both groups, with the exception of some naïve donors where we observed an increase. We did not detect Abs recognizing the N protein for these donors (not shown), suggesting that they had not been infected between the two doses. This increase could therefore be linked to a delayed response or affinity maturation of the antibodies in the germinal center between V1 and V2. The second dose, which was administered ∼16 weeks after the first one, strongly boosted the induction of anti-RBD Igs in the SARS-CoV-2 naïve group, particularly IgG and IgA which reached higher levels, albeit not statistically significant, than those measured three weeks after the first dose (Figure 1D and E). For the PI group, the second dose also led to an increase in the level of total anti-RBD Igs that reached similar levels than after the first dose. Of note, the second dose in the naïve group elicited anti- RBD IgG levels that reached the same levels than in the PI group receiving two doses and significantly higher than PI receiving only one (Figure 1D). Similar patterns of responses were observed when we measured the level of Abs recognizing the full-length S glycoprotein (Figure 1F-I).
Recognition of SARS-CoV-2 Spike variants and other Betacoronaviruses
The BNT162b2 mRNA vaccine has been developed against the original Wuhan strain. However, SARS-CoV-2 is evolving, and many variants have emerged and spread rapidly worldwide. Some harbor specific mutations in S that are associated with increased transmissibility and/or immune evasion (Davies et al., 2020; Sabino et al., 2021; Tegally et al., 2020; Volz et al., 2021). Here, we evaluated the ability of Abs elicited by the Pfizer/BioNTech vaccine to recognize different S proteins of VOCs (B.1.1.7, B.1.351, P.1 and B.1.612.2) and the VOI B.1.526 expressed at the cell surface of 293T cells by flow cytometry, using a method we have previously described (Figure 2, S1) (Gong et al., 2021; Prévost et al., 2020; Tauzin et al., 2021).
293T cells were transfected with the indicated full-length S from different SARS-CoV-2 variants and other human Betacoronavirus Spike and stained with the CV3-25 Ab or with plasma from naïve or PI donors collected at V0, V1, V2 and V3 and analyzed by flow cytometry. The values represent the median fluorescence intensities (MFI) (G,H and I) or the MFI normalized by CV3- 25 Ab binding (A-F, J). Naïve and PI donors with a long interval between the two doses are represented by red and black points respectively and PI donors who received just one dose by blue points. (Left panels) Each curve represents the MFI or the normalized MFIs obtained with the plasma of one donor at every time point. Mean of each group is represented by a bold line. The time of vaccine dose injections is indicated by black triangles. (Right panels) Plasma samples were grouped in different time points (V0, V1, V2 and V3). Undetectable measures are represented as white symbols, and limits of detection are plotted. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
As expected, none of the SARS-CoV-2 naïve plasma samples collected at V0 were able to recognize the SARS-CoV-2 S (D614G) or any of the variants tested here (B.1.1.7, B.1.351, B.1.617.2, P.1, B.1.526) (Figure 2A-F). In contrast, plasma from PI individuals recognized all tested SARS-CoV-2 variants at V0 (Figure 2A-F, S1). The first dose of vaccine strongly enhanced the recognition of the full D614G S and all the tested variants in both groups (Figure 2A-F). Three months after the first dose, the recognition slightly decreased but not significantly. As expected, the second dose strongly increased recognition of all VOC Spikes in the naïve group and reached levels that where significantly higher than after the first dose. In contrast, for the PI group, the second dose did not result in a better recognition than after the first dose. Of note, we observed no significant differences at V3 between PI individuals who received one or two doses, despite a shorter period since the last dose for PI individuals who received two doses. The recognition of all VOCs was slightly lower at V3 by the naïve group compared to the PI that received two doses (Figure 2A-F). When we compared Spike recognition between the SARS-CoV-2 variants, we observed that plasma from PI individuals before vaccination recognized less efficiently the different S variants compared to the D614G S, with the exception of B.1.1.7 S (Figure S1A). After the first and second dose, only B.1.351 and B.1.617.2 S were less efficiently recognized by plasmas from PI individuals (Figure S1B-D). For naïve individuals, even if the vaccination strongly increased the recognition of every VOC Spike tested, we observed that plasmas recognized the different SARS-CoV-2 variants less efficiently compared to D614G S except for the B.1.1.7 S (Figure S1).
We also evaluated whether vaccination elicited Abs that were able to recognize S glycoproteins from endemic human Betacoronaviruses, (HCoV-OC43 and HCoV-HKU1). Interestingly, we observed that the first but not the second dose enhanced the recognition of HCoV-HKU1 S in the naïve group (Figure 2H). Moreover, we observed that plasma from PI donors better recognized HCoV-HKU1 S than plasma from naïve donors at every time point studied, suggesting that natural infection induced cross reactive Abs more efficiently than vaccination. In contrast, both doses did not significantly improve the recognition of HCoV-OC43 S (Figure 2G).
We then evaluated the capacity of the different plasma samples to bind S from highly pathogenic human coronaviruses (SARS-CoV-1 and MERS-CoV). We observed that PI individuals did not have Abs able to recognize MERS-CoV S before vaccination, in contrast to SARS-CoV-1 S (Figure 2I-J). This is likely related the closer genetic relationship between SARS- CoV-2 and SARS-CoV-1 than MERS-CoV (Rabaan et al., 2020; Sarkar et al., 2021). As previously observed (Tauzin et al., 2021), the first dose of vaccine significantly increased the recognition of the MERS-CoV S (Figure 2I) in both groups. However, the second dose did not enhance this recognition. On the other hand, both vaccine doses significantly increased the level of recognition of the SARS-CoV-1 Spike in the naïve group (Figure 2J). In the PI group, only the first dose significantly improved the recognition. We note that the long interval between doses brings SARS-CoV-2 naïve individuals to recognize the different variant Spikes and related HCoV to the same extent than previously-infected individuals.
Functional activities of vaccine-elicited antibodies
We (Tauzin et al., 2021) and others (Collier et al., 2021; Goel et al., 2021b; Planas et al., 2021b; Sahin et al., 2020) reported that three weeks post first Pfizer/BioNTech dose, SARS-CoV- 2 S specific antibodies with weak neutralizing properties are elicited. Nevertheless, these Abs present robust Fc-mediated effector functions as measured by their capacity to mediate antibody- dependent cellular cytotoxicity (ADCC) (Tauzin et al., 2021). To obtain a better understanding of this functional property over time, we tested all plasma samples with our previously reported ADCC assay (Anand et al., 2021; Beaudoin-Bussieres et al., 2021; Tauzin et al., 2021; Ullah et al., 2021). As expected, and in agreement with the absence of SARS-CoV-2 S specific antibodies at baseline, no ADCC activity was observed for the naïve group before vaccination (Figure 3A). Plasma from the PI group maintained some levels of ADCC activity before vaccination, in agreement with a longitudinal study following immune responses in convalescent donors (Anand et al., 2021). Three weeks after the first dose, ADCC activity was elicited in both groups, but significantly higher in the PI group. A decline in ADCC responses was observed in both groups nine weeks after V1 (V2, i.e., 12 weeks post vaccination). The second dose strongly boosted ADCC activity in the naïve group but remained stable for the PI groups. We noted that the levels of ADCC activity were significantly higher in the PI group at all timepoints (Figure 3A).
(A) CEM.NKr parental cells were mixed at a 1:1 ratio with CEM.NKr-Spike cells and were used as target cells. PBMCs from uninfected donors were used as effector cells in a FACS-based ADCC assay. (B) Neutralizing activity was measured by incubating pseudoviruses bearing SARS- CoV-2 S glycoproteins, with serial dilutions of plasma for 1 h at 37°C before infecting 293T-ACE2 cells. Neutralization half maximal inhibitory serum dilution (ID50) values were determined using a normalized non-linear regression using GraphPad Prism software. Naïve and PI donors with a long interval between the two doses are represented by red and black points respectively and PI donors who received just one dose by blue points. (Left panels) Each curve represents the values obtained with the plasma of one donor at every time point. Mean of each group is represented by a bold line. The time of vaccine dose injections is indicated by black triangles. (Right panels) Plasma samples were grouped in different time points (V0, V1, V2 and V3). Undetectable measures are represented as white symbols, and limits of detection are plotted. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
Neutralizing activity in plasma is thought to play an important role in vaccine efficacy (Jackson et al., 2020; Muruato et al., 2020; Polack et al., 2020). Accordingly, it has been recently identified as an immune-correlate of protection in the mRNA-1273 COVID-19 vaccine efficacy trial (Gilbert et al., 2021). To evaluate the vaccine neutralizing response over time, we measured the capacity of plasma samples to neutralize pseudoviral particles carrying the SARS-CoV-2 S D614G glycoprotein (Figure 3B). We did not detect a significant increase in neutralization in plasma isolated three weeks post vaccination of the naïve group, as previously described (Tauzin et al., 2021). Interestingly, nine weeks later (V2, i.e., 12 weeks post vaccination), we observed increased neutralizing activity in a few donors (Figure 3B). All donors presented a significant increase in neutralizing activity three weeks after the second dose. Importantly, the level of neutralizing activity of double vaccinated naïve individuals reached the same levels than in the PI group after one or two doses. In this latter group (PI), we measured low neutralizing activity before vaccination, consistent with remaining neutralizing activity in convalescent donors after several months post symptoms onset (Anand et al., 2021; Gaebler et al., 2021; Tauzin et al., 2021). As previously described, the first dose strongly increased neutralization activity (Stamatatos et al., 2021; Tauzin et al., 2021), but this activity significantly decreased a few weeks after (V2, i.e., 12 weeks post vaccination). The second dose boosted the neutralizing activity to the levels reached three weeks after the first dose. No difference in neutralization was observed between V1 and V3 for PI individuals. In contrast, in naïve individuals we observed a significantly higher neutralizing activity after the second dose compared to the first one (Figure 3B). Thus, while one dose is required to reach maximum neutralization activity in PI individuals, this activity decays over time and a second dose is required to bring back its maximum potential. On the other hand, naïve individuals requires both doses to achieve the same level of PI vaccinated individuals.
Neutralizing activity against variants of concern
SARS-CoV-2 is evolving, and variants of concern are emerging globally (Davies et al., 2020; Prévost and Finzi, 2021; Sabino et al., 2021; Tegally et al., 2020; Volz et al., 2021). To evaluate whether the long interval between the two doses impacted the capacity of vaccine- elicited antibodies to neutralize VOCs and VOI, we measured the neutralizing activity against pseudoviral particles bearing selected variant Spikes (Figure 4, S2). For all the variants tested, we observed a similar pattern than for the D614G S, with neutralizing Abs mainly induced after the second dose in the naïve group (Figure 4A-E). Previously-infected individuals followed a different pattern. While their plasma had some levels of neutralizing activity at baseline, it gained potency and breadth after the first dose. A second dose did not further enhance this activity. Accordingly, we found no significant differences in neutralizing activity in plasma from PIs that received one or two doses measured at V3 (i.e., week 19).
Neutralizing activity was measured by incubating pseudoviruses bearing the indicated SARS- CoV-2 VOCs or SARS-CoV-1 S glycoproteins, with serial dilutions of plasma for 1 h at 37°C before infecting 293T-ACE2 cells. Neutralization half maximal inhibitory serum dilution (ID50) values were determined using a normalized non-linear regression using GraphPad Prism software. Naïve and PI donors with a long interval between the two doses are represented by red and black points respectively and PI donors who received just one dose by blue points. (Left panels) Each curve represents the values obtained with the plasma of one donor at every time point. Mean of each group is represented by a bold line. The time of vaccine dose injections is indicated by black triangles. (Right panels) Plasma samples were grouped in different time points (V0, V1, V2 and V3). Undetectable measures are represented as white symbols, and limits of detection are plotted. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
We also noted that, with the exception of B.1.1.7, plasma from the PI group prior to vaccination (V0) neutralized less efficiently all pseudoviral particles bearing variant Spikes compared to the D614G (Figure S2A). Importantly, both doses boosted the neutralizing activity against all variants and SARS-CoV-1 Spike at V3 (Figure S2D).
Vaccination of PI individuals was shown to increase neutralization against pseudoviral particles bearing the SARS-CoV-1 Spike (Stamatatos et al., 2021; Tauzin et al., 2021). This Spike is used as a representative variant that is even more dissimilar to the vaccine, which was based on the ancestral Wuhan strain. While only one dose was sufficient to provide SARS-CoV-1 neutralizing capacity in PI individuals, two were required in naïve individuals. Remarkably, plasma from naïve individuals reached the same level of neutralizing activity against pseudoviral particles bearing the SARS-CoV-1 Spike than PI. Thus, suggesting that the delayed boosting in naïve individuals allows antibody maturation resulting in enhanced breath (Figure 4F).
Humoral responses in individuals receiving a short dose interval regimen
We also analyzed the humoral responses of 11 SARS-CoV-2 naïve donors from a separate cohort who received their two doses of Pfizer/BioNTech mRNA vaccine three weeks apart (median [range]: 21 days [21–21 days]) (Table 1 and Figure 5A). Unfortunately, for these donors, blood samples were only collected at baseline (V0, prior to vaccination) and 22 weeks (median [range]: 152 days [147–158 days]) after the first dose, thus precluding a direct side-by- side comparison of humoral responses with our cohort of naïve individuals that received the two doses 16 weeks apart. In other words, the V3 of individuals receiving the second dose following a short interval regimen was collected ∼19 weeks post boost whereas the V3 of those with a long interval corresponds to 3 weeks post boost (Figure 1A and 5A). While not a perfect match, we decided to compare the humoral responses of short interval vaccinated individuals with those from PI that received a single dose since their V3 was collected ∼19 weeks post vaccination. While we observed no significant differences in total Ig recognizing the RBD or the full S (Figure 5B and F), single dose PI individuals had more anti-RBD IgG and IgA than naïve donors who received their two doses 3 weeks apart (Figure 5D and E). Despite minor differences in the overall amount of anti-RBD or anti-Spike Abs, we observed major differences related to their capacity to recognize the full Spike. Plasma from short interval vaccinated individuals was significantly less efficient at recognizing the D614G S and all other S variants tested (Figure 5J). Major differences were also observed, particularly related to their capacity to mediate ADCC (Figure 5K) or neutralize pseudoviral particles bearing D614G or any of the variant Spike tested (Figure 5L). Indeed, the naïve donors with a short interval between doses had a very low ID50 against all the variant tested ∼19 weeks after the vaccine series. No neutralizing activity against SARS-CoV-1 was observed (Figure 5L). In contrast, plasmas from PI individual who received just one dose presented neutralizing activity against all the SARS-CoV-2 variants but also the SARS-CoV-1 pseudoviruses (Figure 5L).
(A) SARS-CoV-2 vaccine cohort design. (B-E) Indirect ELISA was performed by incubating plasma samples from naïve and PI donors collected at V3 with recombinant SARS-CoV-2 RBD protein. Anti-RBD Ab binding was detected using HRP-conjugated (B) anti-human IgM+IgG+IgA (C) anti-human IgM, (D) anti-human IgG, or (E) anti-human IgA. Relative light unit (RLU) values obtained with BSA (negative control) were subtracted and further normalized to the signal obtained with the anti-RBD CR3022 mAb present in each plate. (F-I) Cell-based ELISA was performed by incubating plasma samples from naïve and PI donors collected at V3 with HOS cells expressing full-length SARS-CoV-2 S. Anti-S Ab binding was detected using HRP- conjugated (F) anti-human IgM+IgG+IgA (G) anti-human IgM, (H) anti-human IgG, or (I) anti- human IgA. RLU values obtained with parental HOS (negative control) were subtracted and further normalized to the signal obtained with the CR3022 mAb present in each plate. (J) 293T cells were transfected with the indicated full-length S and stained with the CV3-25 Ab or with plasma from naïve or PI donors collected at V3 and analyzed by flow cytometry. The values represent the MFI normalized by CV3-25 Ab binding. (K) CEM.NKr parental cells were mixed at a 1:1 ratio with CEM.NKr-Spike cells and were used as target cells. PBMCs from uninfected donors were used as effector cells in a FACS-based ADCC assay. (L) Neutralizing activity was measured by incubating pseudoviruses bearing SARS-CoV-2 S glycoproteins or SARS-CoV-1 S glycoprotein, with serial dilutions of plasma for 1 h at 37°C before infecting 293T-ACE2 cells. Neutralization half maximal inhibitory serum dilution (ID50) values were determined using a normalized non-linear regression using GraphPad Prism software. PI donors who received one dose are represented by blue points and naïve donors with a short interval between the two doses by green points. Plasma samples were grouped at V3. Undetectable measures are represented as white symbols, and limits of detection are plotted. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
Integrated analysis of vaccine responses elicited with a sixteen-weeks interval between doses
When studying the network of pairwise correlations among all studied immune variables in SARS-CoV-2 naïve individuals (Figure 6), we observed a sparsely interconnected network after the first vaccine dose with focused clusters among binding and neutralization responses, respectively. Over time, the network induced upon 1st vaccination slightly collapsed until the delayed 2nd vaccination triggered a dense network of positive correlations involving binding neutralization responses against several SARS-CoV-2 variants and SARS-CoV-1 (but hardly OC43, HKU1, or MERS spike), ADCC, and memory B cell responses.
Edge bundling correlation plots where red and blue edges represent positive and negative correlations between connected parameters, respectively. Only significant correlations (p < 0.05, Spearman rank test) are displayed. Nodes are color coded based on the grouping of parameters according to the legend. Node size corresponds to the degree of relatedness of correlations. Edge bundling plots are shown for correlation analyses using four different datasets; i.e., SARS-CoV-2 naive individuals at V0, V1, V2 and V3 respectively.
DISCUSSION
The approved regimen of the BNT162b2 mRNA vaccine is the administration of two doses within a short interval of 3-4 weeks. Despite the rapid approval of different vaccine platforms, generating the required doses to immunize the world population represents a daunting task (Moore and Klasse, 2020). Confronted to vaccine scarcity, some jurisdictions decided to increase the interval between doses in order to increase the number of immunized individuals. This decision led to concerns about vaccine efficacy, notably against emergent variants rapidly spreading worldwide and more resistant notably against neutralizing Abs induced by vaccination (Annavajhala et al., 2021a; Planas et al., 2021a; Puranik et al., 2021; Wall et al., 2021; Wang et al., 2021a). Here, we measured the humoral responses of SARS-CoV-2 naïve and SARS-CoV-2 PI individuals who received their two doses sixteen weeks apart.
We observed that in the SARS-CoV-2 naïve group the BNT162b2 mRNA vaccine elicited antibodies with weak neutralizing activity but strong Fc-mediated functions three weeks after the first dose (Tauzin et al., 2021). These responses declined in the following weeks in the absence of a boost. However, administration of the second dose sixteen weeks later strongly enhanced these responses, notably neutralization against some VOCs/VOIs and even the divergent SARS- CoV-1. Therefore, despite initial concerns, the long interval between the doses did not result in poor immune responses, in agreement with recent findings (Parry et al., 2021). The idea behind the strategy of delaying the second dose was to provide some level of immunity to a larger number of individuals than if the second dose would have been saved to administer them three weeks later. However, despite the immunological benefits of increasing the interval between the two doses, this also increases the probability of being infected before the boost.
Several studies have shown that vaccination of previously-infected individuals elicits strong cellular and humoral responses (Efrati et al., 2021; Lozano-Ojalvo et al., 2021; Stamatatos et al., 2021; Tauzin et al., 2021; Urbanowicz et al., 2021). In agreement with these studies, we found that vaccination of these individuals resulted in the induction of strong humoral responses. These responses remained relatively stable over time. We noticed that the second dose did not result in a significant enhancement of these responses, even with a long interval of 16 weeks between doses. Our results demonstrate that, while the second dose boosts the humoral response, PI individuals reach their peak of immunity after the first dose. Altogether, these results suggest that a second dose for PI individuals might be delayed beyond sixteen weeks after the first dose. These observations are in agreement with recent studies showing that PI individuals had maximal humoral and CD4+ and CD8+ T cell responses after the first dose of an mRNA vaccine; the second did not strongly boost these responses (Goel et al., 2021a; Lozano-Ojalvo et al., 2021; Painter et al., 2021).
In contrast, here we show that a delayed second vaccine boost in naïve individuals significantly enhances several immune responses and tightens the network of linear correlations among those. The involved immune variables were humoral and cellular responses directed against SARS-CoV-2, including diverse variants, and SARS-CoV-1, but not or marginally against OC43, HKU1, or MERS. Thus, the potency, quality, and concerted triggering of immune responses appear enhanced in naïve individuals vaccinated with a prolonged interval of 16 weeks between first and second shot, comparable to those obtained after vaccination of previously infected individuals.
We also analyzed humoral responses in a cohort of naïve donors who received their two doses according to the approved three-week interval. Plasma collected ∼19 weeks post second dose, had poor humoral activities with low neutralizing activity against D614G strain and even weaker activity against some VOCs/VOIs including B.1.617.2. These results are in agreement with recent studies showing that circulating antibody levels and neutralizing activity decline over time (Doria-Rose et al., 2021; Goel et al., 2021a). In contrast, we observed that PI individuals who received one dose had better responses 19 weeks after their dose.
Field effectiveness studies in Israël and the USA, where a short interval between doses is recommended, suggest waning protection of the BNT162b2 mRNA vaccine series against non- severe disease after a period of approximately 5 months (CDC, 2021; Goldberg et al., 2021; JCVI, 2021; Tartof et al., 2021). However, SARS-CoV-2 specific memory B cell and CD4+ T cell responses remains stable for the following 6 months, likely protecting from severe disease (Goel et al., 2021a). Whether the same kind of waning will be observed with an extended schedule remains to be evaluated. It will be of critical importance to monitor immune responses and vaccine effectiveness of these extended schedules over time. If the strong humoral response seen with this extended schedule is longer-lasting than immune responses following the authorized schedule, the need of a third dose might be delayed and this could have significant implications regarding control of COVID-19.
To end this pandemic, it will be necessary to rapidly vaccinate the world’s population, including in countries where vaccines are poorly available. The research community around the globe rapidly generated a wealth of data related to vaccine-elicited immune responses and vaccine efficacy. Globally, these results suggest that the current vaccine strategy that was initially deployed could be improved. Our results suggest that modifying the interval at which the two doses are administered might be an important factor to take into account. It will be important to keep in mind that a fine balance needs to be achieved in order to avoid infection between the two doses and at the same time provide sufficient time to elicit optimal humoral responses.
Data Availability
This article includes all datasets generated and analyzed for this study. Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact Author (andres.finzi@umontreal.ca).
AUTHOR CONTRIBUTIONS
A.T. and A.F. conceived the study. A.T., G.B.B., R.G., J.Prévost, M.N., J.R., D.E.K., and A.F. designed experimental approaches. A.T., S.Y.G., G.B.B., D.V., R.G., L.N., L.M., M.Benlarbi, D.C., M.N., A.L., J.Prévost, M.Boutin, G.S., A.N., C.B., Y.B., M.D., D.E.K., and A.F. performed, analyzed, and interpreted the experiments. A.T. and R.D. performed statistical analysis. S.Y.G., G.B.B., A.L., J.Prévost, G.G.L., H.M., G.G., Y.B., J.R., M.C and A.F. contributed unique reagents. J.Perreault, L.G., C.M., P.A., R.B., R.R., G.C.M., C.T. and V.M.-L. collected and provided clinical samples. G.D.S., and N.B. provided scientific input related to VOC and vaccine efficacy. A.T., and A.F. wrote the manuscript with inputs from others. Every author has read, edited, and approved the final manuscript.
DECLARATION OF INTERESTS
The authors declare no competing interests.
STAR METHODS
RESOURCE AVAILABILITY
Lead contact
Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Andrés Finzi (andres.finzi{at}umontreal.ca)
Materials availability
All unique reagents generated during this study are available from the Lead contact without restriction.
Data and code availability
The published article includes all datasets generated and analyzed for this study. Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact Author (andres.finzi{at}umontreal.ca).
EXPERIMENTAL MODEL AND SUBJECT DETAILS
Ethics Statement
All work was conducted in accordance with the Declaration of Helsinki in terms of informed consent and approval by an appropriate institutional board. Blood samples were obtained from donors who consented to participate in this research project at Centre Hospitalier de l’Université de Montréal (CHUM) and at Azienda Socio-Sanitaria Territoriale (ASST) Santi Paolo e Carlo, Milan. The study was approved by the respective Institutional Review Boards (no. 19.381 at CHUM and no. 2020/ST/049 at ASST Santi Paolo et Carlo). Plasma and PBMCs were isolated by centrifugation and Ficoll gradient, and samples stored at -80°C and in liquid nitrogen, respectively, until use.
Human subjects
No specific criteria such as number of patients (sample size), clinical or demographic were used for inclusion, beyond PCR confirmed SARS-CoV-2 infection in adults.
Plasma and antibodies
Plasma from SARS-CoV-2 naïve and PI donors were collected, heat-inactivated for 1 hour at 56°C and stored at -80°C until ready to use in subsequent experiments. Plasma from uninfected donors collected before the pandemic were used as negative controls and used to calculate the seropositivity threshold in our ELISA, cell-based ELISA, ADCC and flow cytometry assays (see below). The RBD-specific monoclonal antibody CR3022 was used as a positive control in our ELISA, cell-based ELISA, and flow cytometry assays and was previously described (Anand et al., 2020; Beaudoin-Bussières et al., 2020; Meulen et al., 2006; Prévost et al., 2020). Horseradish peroxidase (HRP)-conjugated Abs able to detect all Ig isotypes (anti-human IgM+IgG+IgA; Jackson ImmunoResearch Laboratories) or specific for the Fc region of human IgG (Invitrogen), the Fc region of human IgM (Jackson ImmunoResearch Laboratories) or the Fc region of human IgA (Jackson ImmunoResearch Laboratories) were used as secondary Abs to detect Ab binding in ELISA and cell-based ELISA experiments. Alexa Fluor-647-conjugated goat anti-human Abs able to detect all Ig isotypes (anti-human IgM+IgG+IgA; Jackson ImmunoResearch Laboratories) were used as secondary Ab to detect plasma binding in flow cytometry experiments.
Cell lines
293T human embryonic kidney and HOS cells (obtained from ATCC) were maintained at 37°C under 5% CO2 in Dulbecco’s modified Eagle’s medium (DMEM) (Wisent) containing 5% fetal bovine serum (FBS) (VWR) and 100 μg/ml of penicillin-streptomycin (Wisent). CEM.NKr CCR5+ cells (NIH AIDS reagent program) were maintained at 37°C under 5% CO2 in Roswell Park Memorial Institute (RPMI) 1640 medium (Gibco) containing 10% FBS and 100 μg/ml of penicillin- streptomycin. 293T-ACE2 cell line was previously reported (Prévost et al., 2020). HOS and CEM.NKr CCR5+ cells stably expressing the SARS-CoV-2 S glycoproteins were previously reported (Anand et al., 2021).
METHOD DETAILS
Plasmids
The plasmids expressing the human coronavirus Spike glycoproteins of SARS-CoV-2, SARS- CoV-1 (Hoffmann et al., 2013, 2020), HCoV-OC43 (Prévost et al., 2020) and MERS-CoV (Park et al., 2016) were previously reported. The HCoV-HKU1 S expressing plasmid was purchased from Sino Biological. The plasmids encoding the different SARS-CoV-2 Spike variants (D614G, B.1.1.7, B.1.351, P.1, B.1.526 and B.1.617.2) were previously described (Beaudoin-Bussières et al., 2020; Gong et al., 2021; Li et al., 2021; Tauzin et al., 2021).
Protein expression and purification
FreeStyle 293F cells (Invitrogen) were grown in FreeStyle 293F medium (Invitrogen) to a density of 1 x 106 cells/mL at 37°C with 8 % CO2 with regular agitation (150 rpm). Cells were transfected with a plasmid coding for SARS-CoV-2 S RBD (Beaudoin-Bussières et al., 2020) using ExpiFectamine 293 transfection reagent, as directed by the manufacturer (Invitrogen). One week later, cells were pelleted and discarded. Supernatants were filtered using a 0.22 µm filter (Thermo Fisher Scientific). The recombinant RBD proteins were purified by nickel affinity columns, as directed by the manufacturer (Invitrogen). The RBD preparations were dialyzed against phosphate-buffered saline (PBS) and stored in aliquots at -80°C until further use. To assess purity, recombinant proteins were loaded on SDS-PAGE gels and stained with Coomassie Blue.
Enzyme-Linked Immunosorbent Assay (ELISA)
The SARS-CoV-2 RBD ELISA assay used was previously described (Beaudoin-Bussières et al., 2020; Prévost et al., 2020). Briefly, recombinant SARS-CoV-2 S RBD proteins (2.5 μg/ml), or bovine serum albumin (BSA) (2.5 μg/ml) as a negative control, were prepared in PBS and were adsorbed to plates (MaxiSorp Nunc) overnight at 4°C. Coated wells were subsequently blocked with blocking buffer (Tris-buffered saline [TBS] containing 0.1% Tween20 and 2% BSA) for 1h at room temperature. Wells were then washed four times with washing buffer (Tris-buffered saline [TBS] containing 0.1% Tween20). CR3022 mAb (50 ng/ml) or a 1/250 dilution of plasma were prepared in a diluted solution of blocking buffer (0.1 % BSA) and incubated with the RBD-coated wells for 90 minutes at room temperature. Plates were washed four times with washing buffer followed by incubation with secondary Abs (diluted in a diluted solution of blocking buffer (0.4% BSA)) for 1h at room temperature, followed by four washes. HRP enzyme activity was determined after the addition of a 1:1 mix of Western Lightning oxidizing and luminol reagents (Perkin Elmer Life Sciences). Light emission was measured with a LB942 TriStar luminometer (Berthold Technologies). Signal obtained with BSA was subtracted for each plasma and was then normalized to the signal obtained with CR3022 present in each plate. The seropositivity threshold was established using the following formula: mean of pre-pandemic SARS-CoV-2 negative plasma + (3 standard deviation of the mean of pre-pandemic SARS-CoV-2 negative plasma).
Cell-Based ELISA
Detection of the trimeric SARS-CoV-2 S at the surface of HOS cells was performed by a previously-described cell-based enzyme-linked immunosorbent assay (ELISA) (Anand et al., 2021). Briefly, parental HOS cells or HOS-Spike cells were seeded in 96-well plates (4×104 cells per well) overnight. Cells were blocked with blocking buffer (10 mg/ml nonfat dry milk, 1.8 mM CaCl2, 1 mM MgCl2, 25 mM Tris [pH 7.5], and 140 mM NaCl) for 30 min. CR3022 mAb (1 μg/ml) or plasma (at a dilution of 1/250) were prepared in blocking buffer and incubated with the cells for 1h at room temperature. Respective HRP-conjugated Abs were then incubated with the samples for 45 min at room temperature. For all conditions, cells were washed 6 times with blocking buffer and 6 times with washing buffer (1.8 mM CaCl2, 1 mM MgCl2, 25 mM Tris [pH 7.5], and 140 mM NaCl). HRP enzyme activity was determined after the addition of a 1:1 mix of Western Lightning oxidizing and luminol reagents (PerkinElmer Life Sciences). Light emission was measured with an LB942 TriStar luminometer (Berthold Technologies). Signal obtained with parental HOS was subtracted for each plasma and was then normalized to the signal obtained with CR3022 mAb present in each plate. The seropositivity threshold was established using the following formula: mean of pre-pandemic SARS-CoV-2 negative plasma + (3 standard deviation of the mean of pre- pandemic SARS-CoV-2 negative plasma).
Cell surface staining and flow cytometry analysis
293T cells were co-transfected with a GFP expressor (pIRES2-GFP, Clontech) in combination with plasmids encoding the full-length Spikes of SARS-CoV-2 variants or Spikes from different Betacoronaviruses. 48h post-transfection, S-expressing cells were stained with the CV3-25 Ab (Jennewein et al., 2021) or plasma (1/250 dilution). AlexaFluor-647-conjugated goat anti-human IgM+IgG+IgA Abs (1/800 dilution) were used as secondary Abs. The percentage of transfected cells (GFP+ cells) was determined by gating the living cell population based on viability dye staining (Aqua Vivid, Invitrogen). Samples were acquired on a LSRII cytometer (BD Biosciences) and data analysis was performed using FlowJo v10.7.1 (Tree Star). The seropositivity threshold was established using the following formula: (mean of pre-pandemic SARS-CoV-2 negative plasma + (3 standard deviation of the mean of pre-pandemic SARS-CoV-2 negative plasma). The conformational-independent S2-targeting mAb CV3-25 was used to normalize Spike expression. CV3-25 was shown to effectively recognize all SARS-CoV-2 Spike variants (Li et al., 2021).
ADCC assay
This assay was previously described (Anand et al., 2021). For evaluation of anti-SARS-CoV-2 antibody-dependent cellular cytotoxicity (ADCC), parental CEM.NKr CCR5+ cells were mixed at a 1:1 ratio with CEM.NKr cells stably expressing a GFP-tagged full length SARS-CoV-2 Spike (CEM.NKr.SARS-CoV-2.Spike cells). These cells were stained for viability (AquaVivid; Thermo Fisher Scientific, Waltham, MA, USA) and cellular dyes (cell proliferation dye eFluor670; Thermo Fisher Scientific) to be used as target cells. Overnight rested PBMCs were stained with another cellular marker (cell proliferation dye eFluor450; Thermo Fisher Scientific) and used as effector cells. Stained target and effector cells were mixed at a ratio of 1:10 in 96-well V-bottom plates. Plasma (1/500 dilution) or monoclonal antibody CR3022 (1 µg/mL) were added to the appropriate wells. The plates were subsequently centrifuged for 1 min at 300g, and incubated at 37°C, 5% CO2 for 5 hours before being fixed in a 2% PBS-formaldehyde solution. ADCC activity was calculated using the formula: [(% of GFP+ cells in Targets plus Effectors) - (% of GFP+ cells in Targets plus Effectors plus plasma/antibody)]/(% of GFP+ cells in Targets) x 100 by gating on transduced live target cells. All samples were acquired on an LSRII cytometer (BD Biosciences) and data analysis was performed using FlowJo v10.7.1 (Tree Star). The specificity threshold was established using the following formula: (mean of pre-pandemic SARS-CoV-2 negative plasma + (3 standard deviation of the mean of pre-pandemic SARS-CoV-2 negative plasma).
Virus neutralization assay
To produce the pseudoviruses, 293T cells were transfected with the lentiviral vector pNL4.3 R-E- Luc (NIH AIDS Reagent Program) and a plasmid encoding for the indicated S glycoprotein (D614G, B.1.1.7, P.1, B.1.351, B.1.617.2, B.1.526 and SARS-CoV) at a ratio of 10:1. Two days post-transfection, cell supernatants were harvested and stored at -80°C until use. For the neutralization assay, 293T-ACE2 target cells were seeded at a density of 1×104 cells/well in 96- well luminometer-compatible tissue culture plates (Perkin Elmer) 24h before infection. Pseudoviral particles were incubated with several plasma dilutions (1/50; 1/250; 1/1250; 1/6250; 1/31250) for 1h at 37°C and were then added to the target cells followed by incubation for 48h at 37°C. Then, cells were lysed by the addition of 30 µL of passive lysis buffer (Promega) followed by one freeze-thaw cycle. An LB942 TriStar luminometer (Berthold Technologies) was used to measure the luciferase activity of each well after the addition of 100 µL of luciferin buffer (15mM MgSO4, 15mM KPO4 [pH 7.8], 1mM ATP, and 1mM dithiothreitol) and 50 µL of 1mM d-luciferin potassium salt (Prolume). The neutralization half-maximal inhibitory dilution (ID50) represents the plasma dilution to inhibit 50% of the infection of 293T-ACE2 cells by pseudoviruses.
SARS-CoV-2-specific B cells characterization
To detect SARS-CoV-2-specific B cells, we conjugated recombinant RBD proteins with Alexa Fluor 488 or Alexa Fluor 594 (Thermo Fisher Scientific) according to the manufacturer’s protocol. Approximately 2×106 frozen PBMCs from SARS-CoV-2 naïve and prior infection donors were prepared in Falcon® 5ml-round bottom polystyrene tubes at a final concentration of 4×106 cells/mL in RPMI 1640 medium (GIBCO) supplemented with 10% of fetal bovine serum (Seradigm), Penicillin- Streptomycin (GIBCO) and HEPES (GIBCO). After a rest of 2h at 37°C and 5% CO2, cells were stained using Aquavivid viability marker (Biosciences) in DPBS (GIBCO) at 4°C for 20 min. The detection of SARS-CoV-2-antigen specific B cells was done by adding the RBD probes to the antibody cocktail (IgM BUV737, CD24 BUV805, IgG BV421, CD3 BV480, CD56 BV480, CD14 BV480, CD16 BV480, CD20 BV711, CD21 BV786, HLA DR BB700, CD27 APC R700; CD19 BV650, CD38 BB790, CD138 BUV661, CCR10 BUV395, IgD BUV563 and IgA PE). Staining was performed at 4°C for 30 min and cells were fixed using 1% paraformaldehyde (Sigma-Aldrich) at 4°C for 15 min. Stained PBMC samples were acquired on FACSymphony™ A5 Cell Analyzer (BD Biosciences) and analyzed using FlowJo v10.7.1 software.
QUANTIFICATION AND STATISTICAL ANALYSIS
Statistical analysis
Symbols represent biologically independent samples from SARS-CoV-2 naïve individuals or SARS-CoV-2 PI individuals. Lines connect data from the same donor. Statistics were analyzed using GraphPad Prism version 8.0.1 (GraphPad, San Diego, CA). Every dataset was tested for statistical normality and this information was used to apply the appropriate (parametric or nonparametric) statistical test. Differences in responses for the same patient before and after vaccination were performed using Kruskal-Wallis tests. Differences in responses between naïve and PI individuals at each time point were measured by Mann-Whitney (V0, V1 and V2) or Kruskal-Wallis (V3) tests. Differences in responses against the different Spikes for the same patient were measured by Friedman tests. P values < 0.05 were considered significant; significance values are indicated as ∗ p < 0.05, ∗∗ p < 0.01, ∗∗∗ p < 0.001, ∗∗∗∗ p < 0.0001. Spearman’s R correlation coefficient was applied for correlations. Statistical tests were two-sided and p < 0.05 was considered significant.
Software scripts and visualization
Edge bundling graphs were generated in undirected mode in R and RStudio using ggraph, igraph, tidyverse,and RColorBrewer packages (R; R studio). Edges are only shown if p < 0.05, and nodes are sized according to the connecting edges’ r values. Nodes are color-coded according to groups of parameters.
SUPPLEMENTAL INFORMATION
Supplemental information can be found online at …
293T cells were transfected with the indicated Betacoronavirus Spike and stained with the CV3-25 Ab or with plasma collected at V0 (A), V1 (B), V2 (C) and V3 (D) and analyzed by flow cytometry. Plasma recognitions are normalized with CV3-25 binding. Naïve and PI donors with a long interval between the two doses are represented by red and black points respectively, PI donors who received just one dose by blue points and naïve donors with a short interval between the two doses by green points. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
Neutralizing activity was measured by incubating pseudoviruses bearing SARS-CoV-2 S variant or SARS-CoV-1 S glycoproteins, with serial dilutions of plasma collected at V0 (A), V1 (B), V2 (C) and V3 (D) for 1 h at 37°C before infecting 293T-ACE2 cells. Neutralization half maximal inhibitory serum dilution (ID50) values were determined using a normalized non-linear regression using GraphPad Prism software. Naïve and PI donors with a long interval between the two doses are represented by red and black points respectively, PI donors who received just one dose by blue points and naïve donors with a short interval between the two doses by green points. Undetectable measures are represented as white symbols, and limits of detection are plotted. Error bars indicate means ± SEM. (* P < 0.05; ** P < 0.01; *** P < 0.001; **** P < 0.0001; ns, non-significant).
ACKNOWLEDGMENTS
The authors are grateful to the donors who participated in this study. The authors thank the CRCHUM BSL3 and Flow Cytometry Platforms for technical assistance. We thank Dr. Stefan Pöhlmann (Georg-August University, Germany) for the plasmid coding for SARS-CoV-2 and SARS-CoV-1 S glycoproteins and Dr. M. Gordon Joyce (U.S. MHRP) for the monoclonal antibody CR3022. This work was supported by le Ministère de l’Économie et de l’Innovation du Québec, Programme de soutien aux organismes de recherche et d’innovation to A.F. and by the Fondation du CHUM. This work was also supported by a CIHR foundation grant #352417, by a CIHR operating Pandemic and Health Emergencies Research grant #177958, a CIHR stream 1 and 2 for SARS-CoV-2 Variant Research to A.F., and by an Exceptional Fund COVID-19 from the Canada Foundation for Innovation (CFI) #41027 to A.F. and D.E.K. Work on variants presented was also supported by the Sentinelle COVID Quebec network led by the LSPQ in collaboration with Fonds de Recherche du Québec Santé (FRQS) to A.F. This work was also partially supported by a CIHR COVID-19 rapid response grant (OV3 170632) and CIHR stream 1 SARS-CoV-2 Variant Research to MC. A.F. is the recipient of Canada Research Chair on Retroviral Entry no. RCHS0235 950-232424. MC is a Tier II Canada Research Chair in Molecular Virology and Antiviral Therapeutics. V.M.L. is supported by a FRQS Junior 1 salary award. D.E.K. is a FRQS Merit Research Scholar. G.B.B. is the recipient of a FRQS PhD fellowship and J.P. is the recipient of a CIHR PhD fellowship. G.S. is supported by a scholarship from the Department of Microbiologie, Infectiologie et Immunology of Université de Montréal. R.G. and A.L. were supported by MITACS Accélération postdoctoral fellowships. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We declare no competing interests.
Footnotes
↵10 Lead contact
REFERENCES
