Abstract
The SARS-CoV-2 antibody neutralization response and its evasion by emerging viral variants are unknown. Antibody immunoreactivity against SARS-CoV-2 antigens and Spike variants, inhibition of Spike-driven virus-cell fusion, and infectious SARS-CoV-2 neutralization were characterized in 807 serial samples from 233 RT-PCR-confirmed COVID-19 individuals with detailed demographics and followed up to seven months. A broad and sustained polyantigenic immunoreactivity against SARS-CoV-2 Spike, Membrane, and Nucleocapsid proteins, along with high viral neutralization were associated with COVID-19 severity. A subgroup of ‘high responders’ maintained high neutralizing responses over time, representing ideal convalescent plasma therapy donors. Antibodies generated against SARS-CoV-2 during the first COVID-19 wave had reduced immunoreactivity and neutralization potency to emerging Spike variants. Accurate monitoring of SARS-CoV-2 antibody responses would be essential for selection of optimal plasma donors and vaccine monitoring and design.
One Sentence Summary Neutralizing antibody responses to SARS-CoV-2 are sustained, associated with COVID19 severity, and evaded by emerging viral variants
Introduction
Control of the SARS-CoV-2 pandemic relies on population resistance to infection due to a post-infection and vaccination-induced immunity. Current questions relate to the level, breadth, and longevity of generated immunity, and whether mutation of the virus will compromise immunity. Previous studies reported varying results in longitudinal changes of the virus-specific antibody response. Some detected stable antibody titers 4-6 months after diagnosis (1, 2), while others reported waning of the antibody response 2-3 months after infection (3, 4). Differences in assay sensitivity and antigen targets may account for these discrepancies, with Spike and nucleocapsid being the main antigens investigated. Immunoreactivity to other abundant antigens, such as Membrane or Envelope, are unknown. Neutralization of SARS-CoV-2 has been reported for antibodies that bind to Spike, a large homo-trimeric glycoprotein studded across the viral surface (5, 6), whereas Membrane and Envelope proteins, although exposed on the viral surface, remain to be identified as neutralizing antibody targets. Rapid development of neutralizing antibody response to Spike correlates with viral immunity, and individuals who seroconvert may develop a lasting neutralization response (7).
The SARS-CoV-2 virus has accumulated many polymorphisms across its genome, especially within the Spike gene (8). Shortly after the introduction of SARS-CoV-2 into the human population, many early and dominant amino acid polymorphisms were associated with viral entry fitness, such as D614G (9, 10). However, the pressure of the neutralising antibody response might select for escape mutations in Spike that limit post-infectious immunity or vaccine protection (11). One example is the S477N/D614G Spike variant which appeared in Australia during July and August 2020 and was traced to a single event from Australian hotel quarantine (12). The S477N/D614G Spike variant currently represents greater than 5% of Spike variants worldwide, 15% in Europe, and 58% in Oceania (13).
Using the lessons learned from research of other viral pathogens and neuroimmunological autoantibodies (14, 15), we have developed a suite of novel high-content assays that sensitively assess antibody responses against the native oligomeric structure of Spike and its emerging variants (16). To measure the neutralizing capacity, we have also developed a Biosafety Level 2 surrogate Spike-driven virus-cell fusion assay that has been cross-validated with a novel high content, machine-scored, Biosafety Level 3 authentic SARS-CoV-2 neutralization assay.
Herein we characterize the longevity, polyantigenic breadth, and neutralization capacity of the SARS-CoV-2 antibody response in individuals and their responses to globally emerging SARS-CoV-2 variants. Using two longitudinal SARS-CoV-2 community- and hospital-based Australian cohorts representative of the broad spectrum of disease severity at acute infection, we showed that the polyantigenic and neutralizing responses to SARS-CoV-2 are sustained, associated with COVID19 severity, and are evaded by emerging viral variants.
This work provides a community snapshot of humoral immunity in those recovering from infection and sheds light on important considerations for vaccine design and selection of donors for convalescent plasma therapy. Additionally, the modular assays used herein can be adapted for novel viral pathogens to respond rapidly to emerging pathogens.
Results
SARS-CoV-2 antibody responses are sustained for up to seven months post-infection and are focused on Spike
SARS-CoV-2 antibodies were assessed in RT-PCR-confirmed COVID-19 convalescent adults in two Australian cohorts; ADAPT, a hospital-based cohort of patients recruited during the first and second wave of infection in Australia (n=83 and n=17), and LIFE, a cohort of plasma donors (n=159) (Table 1, Fig. 1A). Antibody immunoreactivity to SARS-CoV-2 antigens, inhibition of virus-cell fusion, live SARS-CoV-2 neutralization, and immunoreactivity to Spike emerging variants were assessed and antibody features were compared with demographic data (Fig. 1A). At first date of collection post-infection, 96% (81/83 ADAPT, median 71 days, mean 74 days, post first PCR positivity) and 98% (152/159 LIFE, median 59 days, mean 61 days) of infected patients were Spike IgG-positive, and 81% (66/83 ADAPT) and 91% (139/152 LIFE) were Spike IgM-positive (Table 2). A broad range of Spike IgG levels was observed. No differences in Spike IgG and IgM levels were observed between females and males, but higher IgG and IgM levels were associated with older age (P<0.0001) (Fig. S1). Detection of convalescent positive serostatus was more sensitive when Spike IgG were detected by live cell flow cytometry compared to Nucleocapsid IgG or Spike IgG using commercial assays (Table 2).
The longevity of antibody responses was assessed in 807 Spike IgG-positive serial samples from 233 individuals (n=162 ADAPT, n=645 LIFE), spanning up to 205 days post-PCR positivity (Fig. 1A and Table 1). There was a range of Spike IgG titers at first collection date, and among all Spike IgG-positive individuals, no individual seroreverted, even up to 205 days post-PCR positivity. The majority of ADAPT patients had stable IgG responses (85%), whereas most LIFE donors exhibited decreased IgG over time (59%), where a decrease was defined as >30% change from first collected sample (Fig. 1B) (14). A two-phase decay in those with decreasing responses characterized by an initial high rate of decay followed by stabilization, and the breakpoint between the two phases was estimated at 85 days post-PCR positivity (Fig. 1B). The level at which Spike IgG stabilized was dependent on intial antibody response. High Spike IgG levels decayed to mid-level reactivity and mid-low level reactivity to low level (Fig. 1B). In Spike IgM-positive patients, the majority had decreased IgM levels over time (68% ADAPT; 84% LIFE), in which levels initially decreased and then stabilised at lower levels, but did not serorevert up to 205 days (Fig. 1C). Only five ADAPT (6%) and fourteen LIFE (9%) individuals sero-reverted for Spike IgM at median 146 days post-PCR positivity (Fig. 1C). The breakpoint between the two phases of IgM decay was at 93 days post-PCR positivity.
The polyantigenic breadth of Spike IgG-positive individuals against the virus was examined by detecting IgG targeting the SARS-CoV-2 Membrane, Envelope, and Nucleocapsid proteins (Table 2). 54% (45/83 ADAPT) and 57% (87/152 LIFE) individuals harboured IgG targeting the SARS-CoV-2 Membrane protein, whereas 78% had antibody targeting the Nucleocapsid protein (65/83 ADAPT, 118/152 LIFE) (Fig. 1D and 1E). Antibody titers toward the Membrane protein remained stable over the period of observation in most individuals (91%, 41/45 ADAPT; 95%, 83/87, LIFE), whereas responses toward the Nucleocapsid protein differed between ADAPT and LIFE, and were reminiscent of the Spike IgG response, i.e. mostly stable in ADAPT, and mostly decreased in LIFE over time (Fig. 1D and 1E). Across both cohorts, reactivity to the Envelope protein was very limited with only two ADAPT patients (2%) positive for Envelope IgG (Table 2). Antibody responses to SARS-CoV-2 were highly focused on Spike, followed by the Nucleocapsid and Membrane proteins. Individuals with higher Spike IgG had also high levels of Nucleocapsid and Membrane IgG (Fig. S1).
The overall decay of SARS-CoV-2 antibodies between both cohorts behaved similarly for Spike IgM, but not for Spike IgG, Membrane IgG, and Nucleocapsid IgG, with LIFE donors exhibiting a higher proportion of decreased profiles (Fig. 1). The first collected sample in ADAPT started later post-PCR positivity, and the time duration between paired samples was shorter than for LIFE samples, therefore some ADAPT patients may have been captured during the 2nd, more stable, phase (Fig. 1 and Table 1). Furthermore, few ADAPT patients underwent plasmapheresis, whereas all LIFE donors underwent plasmapheresis as part of convalescent plasma donations (median 6 donations, IQR 3-9, max 14). However, donors with more than 10 donations (n=30) had decay profiles similar to the whole cohort, in which donors stabilized at mid-low level, and none of these highly recurrent donors became seronegative (Fig. S2).
Neutralization of SARS-CoV-2 is correlated with Spike antibody levels and is maintained over time
The neutralization capacity of these individual responses was assessed on a Spike-driven virus-cell fusion assay and a whole-virus neutralization assay (Fig. 1A). Most sera were capable of inhibiting virus-cell fusion (82%, 68/83 ADAPT; 68%, 104/152 LIFE) and mediating viral neutralization (88%, 73/83 ADAPT; 94%, 143/152 LIFE) (Fig. 2A, Table 2). In both cohorts, the virus-cell fusion assay was more stringent than the SARS-CoV-2 neutralization assay as a proportion of individual sera with lower titers in the SARS-CoV-2 neutralization assay were negative in the virus-cell fusion assay (7%, 6/83 ADAPT; 27%, 41/152 LIFE), and most individuals had higher titers in the neutralization assay (Fig. 2A). To understand the discrepancy between both viral assays, live SARS-CoV-2 viral particles were enumerated and directly compared to Spike-pseudotyped lentiviral particles. Cell permeable RNA-specific staining of live virions detected viral particles that were Nucleocapsid-positive (Fig. 2B). The particle to transduction ratios from the fusion assay were 1.03×105, consistent with the low specific infectivity of lentiviruses such as HIV-1(17). In contrast, the SARS-CoV-2 particle to infectivity ranged from 58 (HekAT14) to 578 (VeroE6), consistent with the ratio reported for influenza virus (18). However, the absolute viral particle number was 74-fold higher in Spike-pseudotyped particle preparation (1.64×108 particles per ml) compared to authentic SARS-CoV-2 (2.22×106 particles per ml). Thus, the specific infectivity of SARS-CoV-2 was higher than that of Spike-expressing lentiviral particles, which may account for the higher sensitivity of the SARS-CoV-2-based neutralization assay.
Most ADAPT patients had stable virus-cell fusion inhibition (99%) and neutralization (89%) titers over time (Fig. 2C and 2D). Most of LIFE donors had decreased virus-cell fusion inhibition (82%) and neutralization (56%) capacity over time, and the majority exhibited a single-phase decay in both assays (Fig. 2C and 2D). The greater number of samples per LIFE donor enabled finer characterization of the decay profile in 34 donors in the virus-cell fusion and 44 donors in the neutralization assay). Most donors had a single-phase decay, whilst a two-phase decay was observed in those with >1:320 titers at first collection. These rapidly dropped and then stabilized over time at 1:80 to 1:160 (28% and 25% of LIFE donors in the virus-cell fusion and neutralization assays, respectively). Individuals with two-phase decay had much higher starting titers than individuals with a single-phase decay (Fig. 2D). In the neutralization assay, LIFE donors with decreased profile had a similar median follow up as the stable profile (∼63 days and 56 days, respectively). In both cohorts, the neutralization and fusion profiles were similar to the Spike IgG profiles, in which ADAPT had more stable responses than LIFE. Indeed, Spike IgG and IgM titers were strongly correlated with virus-cell fusion inhibition and SARS-CoV-2 neutralization (Fig. 2E).
A broad antigenic repertoire and high neutralization capacity against SARS-CoV-2 is associated with COVID-19 severity
Approximately half of individuals (55% ADAPT and 49% LIFE) had broad polyantigenic immunoreactivity as defined by IgG responses against each of SARS-CoV-2 Spike, Membrane, and Nucleocapsid proteins (Fig. 3A). Around a third of individuals exhibited antibodies against only two proteins (27% and 30%, Nucleocapsid and Spike; 2% and 9% toward Membrane and Spike in ADAPT and LIFE respectively), and a smaller proportion had responses against Spike alone (12 and 17%) (Fig. 3A). Polyantigenic immunoreactivity did not change overtime in most individuals (82%, 15/81 ADAPT, 83%, 41/152 LIFE, data not shown). No individual developed IgG to new antigens at any point of follow up, but instead, lost immunoreactivity to one antigen, either Nucleocapsid or Membrane.
Patients had broader responses across the spectrum of severity in ADAPT (Fig. 3B). ADAPT and LIFE hospitalized patients with more severe symptoms were more likely to exhibit a broader antibody response to SARS-CoV-2, i.e. polyreactive toward the three antigens (Fig. 3B). Interestingly, two of seven hospitalized LIFE patients who had a short 24 hour hospitalization harboured non-broad responses, and Spike-only responses were exclusively observed in non-hospitalized, mild, and moderate individuals (Fig. 3B). Higher IgG titers against Membrane and Nucleocapsid proteins were also associated with disease severity in both cohorts (Fig. 3C). Patients with broader SARS-CoV-2 responses and higher disease severity had greater viral neutralization and virus-cell fusion inhibition (Fig. 3D, E). This polyreactive, high severity subgroup was populated almost exclusively by older males (Fig. 3F). Similarly, higher neutralization and virus-cell fusion inhibition titers were more enriched in older males with moderate disease and who were hospitalized (Fig. 3G and H).
High responders with strong and broad SARS-CoV-2 antibody responses are ideal plasma donors
A small subgroup of individuals were “high responders” characterized by high Spike IgG, Spike IgM-positive, broad polyantigenic immunoreactivity (binding to Nucleocapsid, Spike, and Membrane), virus-cell fusion inhibition (>1:160), and neutralization (>1:320). They maintained this high response over time (n=14, 17% ADAPT, n=19, 12% LIFE). High responders were more likely to be male, hospitalized, and were of older age (Fig. 4A). Further characterization was performed on a series of increasingly permissive cell lines, VeroE6, HekAT14, HekAT10, and HekAT24 (Fig. 4B and Fig. S3). Low, i.e. non-high responders, and high responders sera neutralized live SARS-CoV-2 in VeroE6, HekAT14, and HekAT10 cell lines, whereas limited neutralization was observed in the hyper-permissive HekAT24 cell line (Fig. 4B). Using the HekAT24 cell line, two elite responders were identified in LIFE (Fig. 4C), with high Spike IgG and IgM levels, and neutralization titers 30-to 4-fold greater than other individuals (Fig. 4D). Interestingly, elite responders had the highest detectable IgM levels, and early IgM decay coincided with a decrease in neutralization titers, whereas Spike IgG remained stable overtime (Fig. 4E). This association between decreased IgM and neutralization titers was observed in ∼10% of individuals in both cohorts (data not shown).
Spike IgG antibody binding and neutralizing capacity are dependent on Spike mutations in emerging new variants
Numerous Spike polymorphisms have evolved over the course of the pandemic (11) with the most attention given to the transmission fitness gain variants, such as D614G (10, 11). To test the breadth of the antibody response, Spike IgG immunoreactivity to several Spike variants implicated in the Receptor Binding Domain (RBD) and S1 was assessed (Fig. 1A). Expression of all Spike variants was similar across each transfected cell line used in the flow cytometry antibody assays (Fig. S4A). Compared to the Wuhan-1 D614 variant, most patients had similar binding and were able to recognise the Spike RBD variants G476F, V483A, and V367S (Fig. 5A). However, across both cohorts, there was an overall reduced binding to D614G, a prominent non-RBD S1 variant present during the Australian first wave (Fig. 5A). 65% of ADAPT and 91% of LIFE individuals, infected from the first world-wide wave, generated antibodies that bound broadly to G476F, V483A, V367S, and D614G Spike, whereas 35% of ADAPT and 9% of LIFE had more restricted Spike recognition; i.e. they recognized G476F, V483A, V367S, but had a decreased binding to D641G (Fig. 5A). Immunoreactivity toward all Spike variants was stable overtime in most patients (data not shown).
Importantly, sera with reduced D614G IgG binding also had lower neutralization and virus-cell fusion inhibition compared to those who recognised D614G Spike (Fig. 5B), suggesting implications for blocking infection in patients who cannot induce robust Spike antibody recognition. Furthermore, patients who bound D614G Spike had broad SARS-CoV-2 polyantigenic immunoreactivity, whereas patients who displayed reduced binding to D614G had more limited antigenic recognition, with 36% recognising Spike only (Fig. 5C). In a D614G virus-cell fusion assay, patients who maintained binding to D614G showed enhanced virus-cell fusion inhibition, compared to when parental Wuhan-1 D614 Spike was used (Fig. 5D). Individuals with lower IgG binding to D614G, i.e. restricted variant recognition had limited D614G Spike virus-cell fusion inhibition, and most (8/11) were unable to prevent Spike fusion (data not shown), emphasizing the need to maintain robust binding to Spike variants for efficient viral neutralization. Patients with restricted Spike variant recognition were not distinguished by age and severity, but were more likely to be female (Fig. S5).
Although D614G Spike remains a predominant variant globally (Fig. 5E), in the second wave of Australian infection between July to September, an isolate with additional polymorphisms, primarily S477N, and in some cases an additional V1068F, was identified. These variants were not detected during the first Australian wave (Fig. 5F), which included the original Wuhan-1 D614 or the D614G variant equally. To assess the antibody binding capacity between original and emerging variants, patients infected by two Spike variants, S477N/D614G and S477N/D614G/V1068F, were recruited during the second wave in Australia (n=17, from the ADAPT cohort, Table 1). All ADAPT patients from the first and second wave had detectable IgG against all Spike variants (Fig. S4B). Compared to the D614G variant, a strong decrease in immunoreactivity to S477N/D614G and S477N/D614G/V1068F was observed in all ADAPT patients from the second wave, whereas the third mutation within the Spike S2 domain V1068F did not have an additive effect (Fig. 5G). This decrease was also observed irrespective of the virus variant that had infected the ADAPT patients (Fig. 5G). Importantly ADAPT patients from the first wave, who had not encountered the new variants, had reduced binding to S477N/D614G and S477N/D614G/V1068F, suggesting a global decrease of immunoreactivity toward both new variants (Fig. 5G).
Discussion
The current study characterizes the breadth, longevity, and neutralizing capacity of SARS-CoV-2 antibody response in two Australian cohorts, encompassing a wide range of demographics and disease states, up to seven months after COVID-19 diagnosis. We show the development of broad and sustained immunoreactivity against SARS-CoV-2 antigens, and found high titers of Spike-binding and virus-neutralizing antibodies were associated with COVID-19 severity. A group of high responders were identified with high, broad, and sustained neutralizing responses, who may represent ideal donors for convalescent plasma donations. Most importantly, although most patients seroconverted, antibodies generated after early infection displayed a significantly reduced antibody binding and neutralization potency to emerging evasive variants. Our data has important implications on hyperimmune therapy, monoclonal antibody treatments, and vaccine development strategies against emerging viral variants.
The longevity of the immune response against SARS-CoV-2 is a fundamental yet currently unresolved question. Like others, we observed a strong correlation between Spike IgG levels and neutralization capacity (19, 20). Although reports on neutralization prevalence and average titers vary widely depending on sampling and detection assay strategies (21, 22), our results expand on previous findings by comparing neutralization levels with antigen-specific response over a longer follow-up period with more timepoints than most previous studies. The decline in IgG titers and neutralization often stabilized at different levels later into convalescence, addressing whether decreasing IgG levels eventually plateau. Especially in LIFE whose samples were collected later post-infection and with a longer follow-up period than ADAPT. Spike IgM levels decreased more rapidly than IgG, but were still detectable up to 205 days after diagnosis, much later than previously reported (1, 19), and consistent with mathematical modelling of decline of IgM titers in a smaller convalescent cohort (23). While our results reveal widely different magnitudes of initial responses and a decrease in neutralizing antibodies titres, most patients have detectable Spike IgG and neutralizing responses more than 5 months after diagnosis, suggesting extended humoral protection, even in those with mild manifestations of the disease.
IgG and IgM against conformational Spike antibody assays have been seldom used, and Spike IgM detection has been challenging. Although many serological assays have reported 100% sensitivity at ∼15 days post-infection (24), prevalence studies, vaccine efficacy, and assessment for convalescent COVID-19 plasma donors may not recruit so early post-infection or -vaccination. In this context, and future seroprevalence studies, more sensitive antibody assays will be essential. Flow cytometry assays are used in clinical diagnostics, mainly in the sensitive and specific detection of neuroimmunological autoantibodies in which antigen conformation and discrimination of seropositive patients from healthy controls are critical (14, 15). Within the follow up time, the detection of Nucleocapsid and Spike IgG by high capacity commercial assays was significantly less sensitive compared to the flow cytometry assay. Integration of the flow cytometry assay to detect Spike IgG would be valuable to include in the diagnostic pipeline in addition to resource-intensive whole virus neutralization. Given the sensitivity of the flow cytometry assay, this methodology would be ideally suited towards seroprevalence in populations to reveal the true rates of community infection.
The majority of individuals in both cohorts were treated in the community. COVID-19 severity, from mild to hospitalization, was associated with an antibody immune response against SARS-CoV-2 that was reactive toward an increasing number of SARS-CoV-2 antigens, as recently reported (25). As our cohorts included only convalescent individuals, the role of broad polyantigenic immunoreactivity in the acute response of hospitalized patients remains unknown. Indeed, reports of patients with absent humoral immune responses have hinted at the role of T cells and innate immune response during the acute disease. Nonetheless, the presence of a broad polyantigenic viral immunoreactivity can be useful to monitor the quality of the antibody response after vaccination.
Whilst the correlation of Spike IgG levels with viral neutralization was strong, high Spike IgM levels were also associated with high viral neutralization in some, especially during the early convalescent days. A lack of somatic mutations was observed in hundreds of cloned neutralizing human antibodies from convalescent patients (26). In addition, many antibody precursor sequences were observed in naïve B cells from pre-pandemic patient samples, highlighting the importance of pre-existing germline antibody sequences in the neutralization response. The lack of somatic mutations observed in IgG may be consistent with IgM being potent in a neutralization response as both isotypes could have similar affinity binding sites for Spike, but with multiple binding sites per molecule on IgM, the avidity for Spike would be higher.
Full virus neutralization and prevention of virus-cell fusion were associated. Whilst many assays aim to assess neutralization surrogates outside of level 3 biosafety laboratories, key differences were observed between Spike-driven virus-cell fusion and the authentic SARS-CoV-2 assay. In our study, the particle to transduction ratio in the virus-cell fusion assay was much higher than the SARS-CoV-2 neutralization assay. This is consistent with the respective infectivity of HIV-1 compared to respiratory viruses such as SARS, non-SARS coronaviruses, and influenza (17, 18). The virus-cell fusion assay involves a single round of infection, whereas the full virus in the neutralization assay is replication-competent and undergoes multiple rounds of replication over a three day culture. Therefore, the spread of the virus must be considered alongside the capacity of antibodies to inhibit the initial single particle entry and blocking of the virus spread between cells. Although the pseudotyping fusion assay had lower sensitivity, most individuals across both cohorts had titers in this assay with potency ranking similar to full virus neutralization.
Transfusion of convalescent COVID-19 plasma has been proposed as a therapy, with >70 ongoing randomized controlled trials. The few clinical trials to date have supported an acceptable safety profile, but evidence regarding efficacy is mixed (27). Most unsuccessful trials included donors with unknown neutralizing status or FDA-classified low titers, whereas successful trials either used high-titer or convalescent COVID-19 plasma delivered within three days of hospitalization (28). Indeed, utility of convalescent plasma could be improved by donations during early disease stages and by selecting donors with high neutralization antibody titers (29). Our findings that the immunological response to SARS-CoV-2 is widely heterogeneous, with large variations in SARS-CoV-2 antibodies and neutralization, polyantigenic immunoreactivity, and longitudinal responses complement these assertions. To take into account the first phase of decay observed during early convalescence, we propose an optimal window for plasmapheresis, up to 100 days post-diagnosis. Furthermore, the occurrence of a small group of individuals, termed “high and elite responders” with highly neutralising, broad, and sustained SARS-CoV-2 antibody responses over time, may be due to the rapid and lasting generation of memory B cells (30). These patients were likely to be hospitalized older males. Alongside appropriate serology screening programs, the targeted recruitment for plasma donations could help to identify optimal convalescent resources available within affected communities.
A clear advantage of the methodologies used in this study is the capacity of both level 2 biosafety pseudotyped fusion and flow cytometry assays to monitor the effects of viral polymorphisms in real time. Indeed with acceleration of global viral spread, we are now observing evolution of viral fitness and/or immune escape across millions of infected people. The SARS-CoV-2 fitness gain of D614G (9) appeared very early in the pandemic and still represents the majority of viral infections globally (>80%)(10). Zoonosis of a virus is often followed by finer tuning of replication, as observed in the 2014 Ebola outbreak, in which the variant A82V enabled more efficient receptor NPC1 usage (31). Although D614G is a polymorphism outside of the RBD, it significantly impacts the RBD positioning and Spike quaternary structure. The release of hydrogen bonds leading to structural changes are proposed to expose Spike to increase ACE2-dependent fusion (32). RBD exposure in the D614G variant may explain the association with great inhibition of virus-cell fusion in patients who recognized the D614G Spike variant. These results are consistent with recent studies in hamsters (9) and preliminary data on protection from the ongoing vaccine human trials in areas where the D614G Spike variant remains prevalent. However, our data also highlighted that a subgroup of patients who displayed limited antibody binding to D614G Spike also had reduced virus neutralization irrespective of the viral variant that had infected them. This could be a major concern for vaccine candidate design, especially given the emergence of the S477N/D614G polymorphism in the majority of patients that were infected in the Australian second wave and in Europe (12, 13). Seroconversion was observed in all patients from the first and second wave, and good antibody binding to Wuhan-1 D614 and D614G, but there also was a significant decrease in binding to S477N/D614G Spike independent of the variant that had infected individuals. Therefore the emergence of the additional polymorphism S477N/D614G could represent an immune evasive variant leading to less antibody binding and a resistance to virus neutralization, which could imply a need for periodic variation in vaccine design, as for the influenza vaccine (33). Indeed, the mapping of Spike monoclonal antibody escape in vitro has recently shown that S477N/D614G is broadly resistant to many neutralizing antibody clones (34). Whilst the mechanism behind these observations is unknown, the appearance of a N-glycosylation site within Spike RBD could lead to glycan shielding, as in HIV (35), and our evidence that S477N/D614G-infected patients have a similar binding to this variant, albeit reduced, compared to first wave patients, may suggest changing the Wuhan-1 D614 Spike to the S477N/D614G variant in vaccine generation may not overcome the resistance of this variant to the neutralising antibody response. As antibodies against Spike harness the majority of neutralizing activity, selecting the optimal Spike variants in monovalent or multivalent vaccine strategies may be critical.
Our study has important translatable implications to understand the natural history of COVID-19, and post-infection and vaccination-induced immunity. We have highlighted that molecular epidemiology and sero-surveillance will both be required to detect emerging polymorphisms. Furthermore, sensitive monitoring of antibody binding and neutralization capacity will be paramount in vaccine design strategy and convalescent plasma therapy, and in seroprevalence studies, and this would require involvement of more rapidly adaptive methodologies to characterize the magnitude of the neutralization antibody responses against emerging variants.
Materials and Methods
Subjects
This study investigated two cohorts of RT-PCR-confirmed convalescent individuals recruited from February to October 2020 in Australia (Table 1 and Fig. 1A). The Adapting to Pandemic Threats (ADAPT) cohort included 83 patients diagnosed at a community-based fever clinic whose sera was collected at two time points post PCR-positivity during the first wave (March-August, n=166 samples). The second wave included sera from 17 patients recruited between July to October. The Australian Red Cross Lifeblood (Lifeblood) cohort (LIFE) included 645 sera samples from 159 donors collected at multiple timepoints post-PCR-positivity (at least 28 days post-recovery) from volunteers presenting to Lifeblood for whole blood or plasma donation. The disease severity of ADAPT patients ranged from mildly symptomatic (mild), community-managed (moderate) to critically unwell and hospitalized (hosp), whereas the self-reported disease severity of LIFE donors included community-managed (non-hosp) and hospitalized (hosp) (Table 1). A healthy adult non-infected pre- pandemic cohort was collected in Australia and consisted of healthy and non-inflammatory neurological disorder donors (n=24). No re-exposure to SARS-CoV-2 and no re-infection was reported. Ethics approval for this study was granted by St Vincent’s Hospital (2020/ETH00964) and Lifeblood (30042020) Research Ethics Committees. Written consent was obtained from all ADAPT patients. In LIFE, the donor consent form included a statement that blood donation may be used in research.
Flow cytometry cell-based assay for detection of SARS-CoV-2 antibodies
A flow cytometry cell-based assay detected patient serum antibodies against SARS-CoV-2 antigens as for neuroimmunological autoantibodies (14, 15). SARS-CoV-2 full-length Spike (Wuhan-1 D614, V367F, G476S, V483A, D614G, S477N/D614G, and S477N/D614G/V1068F) (10, 11), Membrane, and Envelope proteins were expressed on transfected HEK293 cells. Serum (1:80) was added to live Spike-expressing cells, and Membrane-, and Envelope-expressing cells were treated with 4% paraformaldehyde and 0.2% saponin, followed by AlexaFluor 647-conjugated anti-human IgG (H+L) (ThermoFisher Scientific) or anti-human IgM (A21249, ThermoFisher Scientific). Cells were acquired on the LSRII flow cytometer (BD Biosciences). Patients were SARS-CoV-2 antibody-positive if their delta median fluorescence intensity (ΔMFI □=□ MFI transfected cells – MFI untransfected cells) was above the positive threshold (mean ΔMFI+4SD of 24 pre-pandemic controls) in at least two of three quality-controlled experiments (14). Binding to Spike variants was expressed as a percentage of reduced binding compared to Spike. Data was analysed using FlowJo 10.4.1 (TreeStar, USA), Excel (Microsoft, USA), and GraphPad Prism (GraphPad Software, USA).
Commercial SARS-CoV-2 ELISA
Nucleocapsid IgG assay on the ARCHITECT-I (Abbott Diagnostics, USA), quantitative Spike-1/Spike-2 (S1/S2) IgG on LIASON-155 XL (DiaSorin S.p.A, Italy), and Spike (S1) IgG immunoassay (EUROIMMUN, Germany) were performed. Samples were reported positive if the signal was greater than the published cut-off value (>1.4). Signal to cut-off ratios were used.
SARS-CoV-2 viral-cell fusion assay
The hACE2 ORF (Addgene# 1786) was cloned into a 3rd generation lentiviral expression vector and clonal stable ACE2-expressing Hek293T cells were generated by lentiviral transductions (36). Lentiviral particles pseudotyped with SARS-CoV-2 Spike envelope were produced by co-transfecting Hek293T cells with a GFP encoding lentiviral plasmid HRSIN- CSGW (37), psPAX2, and plasmid expressing C-terminal truncated Spike (pCG1-SARS-2-S Delta18) (38) including D614 or D614G (36). Neutralization activity of sera was measured using a single round infection of ACE2-HEK293T with Spike-pseudotyped lentiviral particles. Virus particles were incubated with serially diluted donor sera for 1 hour at 37°C. Virus-serum mix was then added onto ACE2-HEK293T cells (2.5×103/well) in a 384-well plate. Following spinoculation at 1200g for 1 hour at 18°C, the cells were moved to 37°C for 72 hours. Entry of Spike particles was imaged by GFP-positive cells (InCell Analyzer) followed by enumeration with InCarta software (Cytiva, USA). Neutralization was measured by reduction in GFP expression relative to control group infected with the virus particles without any serum treatment.
The virus entry pathway in VeroE6, used in live virus neutralization assays, is primarily endosomal (39). In contrast, cells derived from nasopharyngeal tissues, express ACE2 in addition to the surface serine protease TMPRSS2 which drives virus-cell membrane fusion and can signficantly enhance viral entry (38). To address viral neutralization in the presence of ACE2 and TMPRSS2, a portfolio of Hek293T expressing clonal cell lines with ACE2 and TMPRSS2 (HekAT) was generated. The coexpression of ACE2 and TMPRSS2 led to a series of increasingly permissive cell lines that were readily susceptible to SARS-CoV-2 cytopathic effects, VeroE6, HekAT14, HekAT10, HekAT24 (Fig. S3).
High content fluorescdent live SARS-CoV-2 neutralization assay
Sera were serially diluted and mixed in duplicate with an equal volume of virus solution at 1.5×103 TCID50/mL. After 1 hour of virus-serum coincubation at 37°C, 40μL were added to equal volume of freshly-trypsinized VeroE6 cells, and three clonal HekAT cells in 384-well plates (5×103/well) selected on SARS-CoV-2 permissiveness. After 72h, cells were stained with NucBlue (Invitrogen, USA) and the entire well was imaged with InCell Analyzer. Nuclei counts, proxy for resulting cytopathic effect, were compared between convalescent sera, mock controls (defined as 100% neutralisation), and infected controls (defined as 0% neutralization) using the formula; % viral neutralization = (D-(1-Q))x100/D, where Q = nuclei count normalized to mock controls, and D = 1-Q for average of infection controls (InCarta software).
Enumeration of SARS-CoV-2 particles
Live SARS-CoV-2 and lentiviral particles were stained using SYTO™ RNASelect™ Green Fluorescent cell Stain (Invitrogen, USA) at a final concentration of 10µM for 30 minutes at 37°C in freshly thawed unpurified viral particles. Particles were then diluted 1/10 and 1/100 in sterile PBS and then adhered to Poly-L-Lysine coated glass bottom 96-well Greiner Sensoplates (Sigma Aldrich, USA) through spinoculation at 1200g for 1 hour at 18°C. Particles were either imaged live or immune-fluorescently counter-stained using a rabbit polyclonal SARS-CoV-2 Nucleocapsid antibody, followed by Alexa647-conjugated goat anti-rabbit IgG (Novus Biologicals, USA). Viral particles were then imaged and quantified as previously described (40). Particle to infectivity ratios were determined by dividing the total particle count per ml with the calculated TCID50/ml. Particle to GFP transduction ratios were used for lentiviruses.
SARS-CoV-2 Spike sequencing and analysis
Clinical respiratory samples were sequenced using an existing amplicon-based Illumina sequencing approach. The raw sequence data were subjected to an in-house quality control procedure before further analysis as reported in (41). Non-synonymous SARS-CoV-2 Spike mutations (read frequency >0.8, minimum coverage 10x) were inferred from variant calling files during bioinformatic analysis using phylogenetic assignment of named global outbreak lineages (PANGOLIN)(11). All consensus SARS-CoV-2 genomes identified have been uploaded to GISAID (www.gisaid.org).
Statistics
Statistical analyses were performed in R v4.0.3. Loess curves were generated using ggplot2 v3.3.2. For categorical variables, a log-linear model was fitted and Pearson residuals plotted in a mosaic plot (MASS v7.3-51.6). Shapiro-Wilk test was used to test for normality in continuous variables and a Dwass-Steel-Critchlow-Fligner test was used to test for significance between continuous and categorical variables. Correlations were measured using the Spearman method (psych v2.0.8). Virus-cell fusion and neutralization data were fitted using an exponential decay curve (Origin Lab). Patient curves unable to be fitted, <3 collection dates or low viral fusion and neutralization, were undetermined. Statistical significance was determined as p<0.05.
Data Availability
Correspondence and requests for data should be addressed to FB and may be limited due to ethical considerations. Plasmids transfer should be obtained through a MTA.
Funding
This work was supported by Snow Medical (Australia), The University of New South Wales Rapid Response grant (Australia), the University of Sydney Research Excellence Initiative grant (Australia), and the MRRF NHRMC COVID-19 grant. The Australian Governments fund Australian Red Cross Lifeblood for the provision of blood, blood products, and services for the Australian community;
Author contributions
FB, SGT, AK, and GM designed the study. FT, DP, DV, AOS, AA, VM, FXL, GW, RR, VS, WR, PC conducted and analysed experiments. GD, GM, DRD, IG, VH, DI enrolled and managed patients. FT, AOS, AA, and DRD wrote the manuscript first draft and FT prepared Figs and tables. FB, SGT, AK, and GM designed and coordinated research and verified results. All authors reviewed the draft before submission;
Competing Interests
FB has received honoraria from Biogen Idec and Merck Serono as invited speaker. All other authors declare no competing interests;
Data and materials availability
Correspondence and requests for data should be addressed to FB and may be limited due to ethical considerations. Plasmids transfer should be obtained through a MTA.
Acknowledgements
We thank all the patients and donors who participated in this study. We thank Drs. Suat Dervish, Edwin Lau, and Maggie Wang for providing advices at the Flow Cytometry Core Facility of the Westmead Research Hub. We thank Ms Rebecca Rielly and the Operation Team at Kids Research for providing acces to the PC2 facility;