key: cord-0732141-5f9fgj05 authors: Whitcombe, A. L.; McGregor, R.; Craigie, A.; James, A.; Charlewood, R.; Lorenz, N.; Dickson, J. M. J.; Sheen, C. R.; Koch, B.; Fox-Lewis, S.; McAuliffe, G.; Roberts, S. A.; Morpeth, S. C.; Taylor, S.; Webb, R. H.; Jack, S.; Upton, A.; Ussher, J.; Moreland, N. J. title: Comprehensive analysis of SARS-CoV-2 antibody dynamics in New Zealand date: 2020-12-11 journal: nan DOI: 10.1101/2020.12.10.20246751 sha: 2fee2ddb1935c7cb2d2d6ef808631f163a7526c6 doc_id: 732141 cord_uid: 5f9fgj05 Objectives Circulating antibodies are important markers of previous infection and immunity. Questions remain with respect to the durability and functionality of SARS CoV 2 antibodies. This study explored antibody responses in recovered COVID 19 patients in a setting where the probability of re-exposure is effectively nil, owing to New Zealand's successful elimination strategy. Methods A triplex bead-based assay that detects antibody isotype (IgG, IgM and IgA) and subclass (IgG1, IgG2, IgG3, IgG4) responses against Nucleocapsid (N) protein, Receptor Binding Domain (RBD) and Spike (S) protein of SARS-CoV-2 was developed. After establishing baseline levels with pre-pandemic control sera (n=113), samples from PCR-confirmed COVID-19 patients with mild-moderate disease (n=189) collected up to eight months post infection were examined. The relationship between antigen-specific antibodies and neutralising antibodies (NAbs) was explored with a surrogate neutralisation assay that quantifies inhibition of the RBD/hACE 2 interaction. Results While most individuals had broad isotype and subclass responses to each antigen shortly after infection, only RBD and S protein IgG, as well as NAbs, were stable over the study period, with 99%, 96% and 90% of samples, respectively, having responses over baseline 4-8 months post-infection. Anti-RBD antibodies were strongly correlated with NAbs at all time points (Pearsons r > 0.87) and feasibility of using finger prick sampling to accurately measure anti-RBD IgG was demonstrated. Conclusion Antibodies to SARS-CoV-2 persist for up to eight months following mild to moderate infection. This robust response can be attributed to the initial exposure without immune boosting given the lack of community transmission in our setting. It is now well established that antibody responses against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) are activated promptly after infection. 1 The virus, and causative agent of the Coronavirus Disease 2019 (COVID-19) pandemic, contains four structural proteins, the most immunodominant being the Nucleocapsid (N) protein, Spike (S) protein and the receptor binding domain (RBD) of the S protein. Measuring antibody responses to these antigens using serological assays has been critical for determining previous viral exposure in individuals, studying community transmission, and conducting population sero-surveys. 1 However, much is still to be learnt about the long-term duration and protective capacity of these responses. Antibody responses comprise different isotypes and subclasses, each associated with unique immune functions and dynamics over time. 2 Following SARS-CoV-2 infection there is an almost concurrent rise in Immunoglobulin M (IgM), IgA and IgG, with IgM then beginning to decline approximately three weeks after symptom onset. [3] [4] [5] [6] There have been conflicting reports with respect to IgG duration, ranging from a relatively short three months 7 , to six months or longer 8-10 , in part due to SARS-CoV-2 antibody dynamics being highly antigen-dependant. Anti-N antibodies are now known to wane faster than anti-RBD and anti-S, and may be more suited as a marker of recent COVID-19 infection, particularly since the N protein is associated with RNA packaging and anti-N antibodies are non-neutralising. 8, 11 The RBD of the Spike protein on the other hand binds to the Angiotensin-Converting Enzyme-2 (hACE2) on human host cells to facilitate viral entry and infection. Anti-S and anti-RBD antibodies can block this interaction leading to viral neutralisation, and as such, are better markers of functional immune responses. 12 Of note, the presence of neutralising antibodies (NAbs) protected a small number of individuals from re-infection during a SARS-CoV-2 outbreak on a fishing vessel 13 and anti-S protein IgG was associated with reduced re-infection in a recent study of health care workers . CC-BY-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020. 12.10.20246751 doi: medRxiv preprint in the United Kingdom. 14 This suggests NAbs and S protein antibodies are associated with protective immunity and may, in turn, underpin a correlate of protection for COVID-19 vaccine development. Levels of NAbs have proven relatively stable in recent reports 9,10,15 but further investigation is needed, particularly in non-severe cases of COVID-19. Bead-based serological assays capable of simultaneously measuring antibodies to N and S proteins, together with RBD, enable a comprehensive view of the SARS-CoV-2 antibody response. [16] [17] [18] In this study, a triplex Luminex-based assay that detects isotype and subclass responses to the major SARS-CoV-2 antigens was developed and utilised to interrogate the composition and duration of virus-specific antibodies up to eight months post-infection in COVID-19 cases in New Zealand. In parallel, levels of NAbs were measured using a surrogate Virus Neutralization Test (sVNT) previously shown to strongly correlate with neutralisation measured using live SARS-CoV-2 virus, and better suited to high throughput analyses due to improved speed and reproducibility. 12, 19 By exploring the relationship between NAbs and antigen-specific antibody features, this study adds to growing knowledge of the SARS-CoV-2 humoral immune response from a setting where the probability of re-exposure is effectively nil, owing to New Zealand's successful elimination strategy. 20, 21 . CC-BY-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint A triplex bead-based immunoassay was developed, with N protein, trimeric S protein and RBD coupled to spectrally unique beads. Compatibilty of the beads in a multiplex format was confirmed, as was comparability with previously described ELISA for S protein and RBD 22 , and N protein titres determined using the Abbott Architect SARS-CoV-2 IgG assay 23 , with highly significant correlations for all three antigens (Supplementary Figure 1) . The level of SARS-CoV-2 antigen specific isotypes (IgG, IgM, IgA) and subclasses (IgG1, 2, 3 and 4) were determined in a cohort of 112 PCR-confirmed COVID-19 participants, 50 of whom had multiple time points (n=189 samples), and the majority of whom had mild disease ( Table 1) . The days post-infection ranged up to 246 days, enabling the kinetics of 21 different antibody features (3 SARS-CoV-2 antigens and 7 secondary detectors) to be temporally investigated. The cut-off for positivity for each antibody type was determined using a previously described panel of pre-pandemic control samples (n=113) that includes donors with respiratory viruses and bacterial pneumonia that have symptom overlap with COVID-19. 22 Isotype responses for each antigen were plotted against time ( Figure 1 ). Participants were grouped by days/months post-infection (7-62 days (0-2 months, n=27), 63-124 days (2-4 months, n=79) and 125-250 days (4-8 months, n=80)) and the rate of positivity compared (Supplementary Table 1 ). With the exception of anti-N protein IgM, all antigen-specific isotypes were elevated above baseline levels in the acute stage through to early convalescence is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint days), respectively. Consistent with recent literature 8,11 , anti-N protein IgG waned faster, with only 54% (43/80) of patients having detectable levels above baseline in late convalescence. To explore the contribution of IgG subclasses to total IgG (IgGtot), the proportion of samples with antibody subclass levels above baseline in the three time groups were compared ( Figure 2 , Supplementary Table 1 ). The four IgG subclasses were elevated against each antigen in the majority of samples collected 7-63 days post-symptom onset. IgG3 responses were particularly strong in this earliest time group, with at least 96% of patients having a response above baseline against both RBD and S protein (Supplementary Table 1 ). In samples collected 63-124 days post symptom onset, IgG1 and IgG3 responses dominated for all three antigens, with IgG2 uncommonly, and IgG4 rarely, above baseline. In the 125-250 day group IgG1 antibodies against RBD and S protein persisted, with 64% (51/80) and 59% (47/80) of participants having antibodies above baseline, respectively, as did IgG3 against S protein with 60% (48/80) above baseline. As with IgGtot, the anti-N subclass responses waned faster than RBD and S protein with <30% participants having detectable levels above baseline of any subclass in the 125-250 day group. NAbs were measured using a surrogate assay that quantifies inhibition of the RBD/ACE-2 interaction 12 and found to be relatively stable over the eight month study period ( Figure 3 ). While there was a modest downward trend, 89% (70/79) in the 63-124 day group and 90% (72/80) in the late convalescent group (>125 day) had detectable NAbs above the cut-off. All bar one of the late convalescent samples (79/80) were from non-hospitalised individuals suggesting mild disease can induce long-lived NAbs. Correlation analysis found IgGtot-RBD most significantly correlated with NAbs in the early, mid-and late-convalescent groups (Pearson's r ≥ 0.87 across all time groups) while IgA-RBD was moderately associated, and IgM-RBD not at all ( Figure 3B ). To further explore NAb persistence, a comparative analysis was performed on participants for which two or more samples were collected across the study . CC-BY-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint period (Supplementary Figure 2) . There was no significant decrease in NAbs between the first and second time point, but a significant decline at the third time point (median 219 days post-infection, P < 0.05). Individual heterogeneity was evident in terms of both the level of NAbs, and their waning (or otherwise) over time. Indeed, approximately 1/3 of individuals demonstrated increasing levels, consistent with affinity maturation as recently proposed. 24 To estimate the rate of decay for NAbs a series of models were applied (see methods), and the expected time for NAbs to halve (half-life (t1/2)) was calculated ( Figure 3C , Supplementary Table 2 ). Exponential decay assumes NAbs decline immediately after infection, and estimated the longest t1/2 of 625 days and the largest margin of error (95% CI, 319-13465 days). A growth and decay model that assumes NAbs increase initially and then decline, estimated a shorter t1/2 of 425 days (95% CI, 253-1316). Using only data from the 50 participants for which two or more samples were available, excluding those with increasing NAbs (i.e. still in the growth phase), and modelling individual decay, estimated a much shorter t1/2 of 146 days (95% CI, 100-199). Lastly, since IgGtot persist, and the bead-based assay format is compatible with small quantities of sera (2µL), the feasibility of using dried blood samples to measure IgGtot was explored. Dried blood samples could expand the acceptability of serological assays and enable testing in settings where it maybe logistically challenging to collect, process and store venous blood. 25 Simulated dried blood eluents collected using Mitra samplers were compared to matched sera in a subset of patients (n=19). IgGtot measurements from Mitra eluents correlated strongly with serum samples, against all three antigens (r 2 = 0.9928, 0.9964 and 0.9931 for N, RBD and Spike, respectively) (Figure 4) . is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint The duration of antibody responses following SARS-CoV-2 infection is, in essence, being studied in real-time globally. Initial studies that purported early antibody waning were based on the N protein, now known to induce antibodies with a relatively short half-life 8 , or linear extrapolations of S protein data collected 2-3 months post-infection. 7 More recent analyses (5-8 months post-infection) show that S protein, RBD and NAbs are likely more stable than originally predicted. 6, 8, 10, 15, 18, 26 Here, further evidence is provided of the durability of the SARS-CoV-2 humoral immune response over an eight month period. Importantly, this study was performed in a setting where there is next to no likelihood of immune boosting by reexposure, given the successful public health response to COVID-19 in New Zealand. 20,21 That 99% of sera had anti-RBD IgG and 96% had anti-S protein IgG above baseline levels 4-8 months after infection in a setting with no circulating SARS-CoV-2 is notable, particularly given the majority of participants had non-severe infections. Utilising a multiplex assay enabled not only IgG, but other isotypes and subclass responses to be investigated. As expected, IgM had a shorter duration than IgG, with >75% of samples returned to baseline in late convalescence for all three antigens. Of the IgG subclasses, RBD and S protein specific IgG1 and IgG3 were dominant, in keeping with the association of these subclasses with convalescence and survival in COVID-19 systems serology studies. 16 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint strategies. 15 As such there is a pressing need to understand antibody persistence and time to seroreversion following SARS-CoV-2 infection. Using three different models the estimated t1/2 of 5-20 months for NAbs in this study is somewhat longer than 2-4 months calculated in other cohorts 18, 26 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint Plasma and serum samples were collected from multiple sources with appropriate ethical board approval. Samples collected before the pandemic were used as negative controls ("prepandemic," n=113), details of which have been described previously. 22 Samples from PCRconfirmed COVID-19 individuals were obtained from hospital in-patients in Auckland (n=18) (ethics HDEC 20NTB76) and convalescent participants in the Southern Region of New Zealand (ethics HDEC 20NTB101) as previously described. 22, 23 Samples were also obtained from donors with a prior COVID-19 diagnosis, collected at the New Zealand Blood service as part of a Medsafe approved process for convalescent plasma preparation. The final PCRconfirmed cohort in this study comprised 189 samples from 112 participants, 50 of whom had samples collected at multiple time points (Table 1) . Participant samples were grouped by days/months post symptom onset (infection) based on timeframes in which the samplers were obtained (acute and early convalescence, 7--62 days (0-2 months, n=27); mid-convalescence, 63-124 days (2-4 months, n=79); late convalescence, 125-250 days (4-8 months, n=80)). SARS-CoV-2 RBD and the trimeric S protein antigens were recombinantly expressed and purified from HEK293T or HEK293F cells as previously described. 22 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint temperature. Antigen-coupled beads were washed, enumerated using a haemocytometer, and stored at 4⁰C protected from light until further use. All serum and plasma samples were heat-inactivated at 56°C for 30 min prior to use. Samples is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint days post-infection (n=18), and to maximise sensitivity, the resulting 98% specificity cut-off for each antibody was taken as the baseline. Surrogate neutralisation assays were performed using the SARS-CoV-2 sVNT (GenScript) according to the manufacturers' instructions as previously described. 12, 22 Briefly, samples were diluted 1:20 and incubated with an equal volume of peroxidase-conjugated RBD for 30 min at 37°C, then added to wells pre-coated with human ACE-2 receptor protein and incubated for 15 min at 37°C prior to washing and development. Inhibition was calculated as (1-OD sample/OD of negative control) × 100. Any sample with a percentage inhibition ≥20% was deemed positive for NAbs, with this manufacturers cut-off previously shown to result in 100% specificity in our laboratory. 22 To simulate remote blood sampling from a fingerprick using volumetric absorptive microsampling, Mitra devices (Neoteryx) were placed onto the surface of whole blood collected from an EDTA tube and filled. The devices were left to dry at room temperature and stored for a maximum of 3 weeks or frozen at -20°C until use. To elute the dried blood, the tip of the Mitra device was placed in 400ul of PBS +1% BSA + 0.05% Tween 20 and incubated overnight at 4°C, with shaking (300rpm). Eluents were stored at 4°C for up to one week or at -80°C until use. The 'serum equivalent' dilution of eluents was calculated to be 1:40, based on the assumption that serum constitutes 50% of the whole blood volume. To compare antigenspecific IgG levels in eluents with levels measured in serum, eluents were diluted appropriately and measured alongside matched serum samples in the multiplex assay. Overall differences in raw values between eluents and sera were low for each antigen (average CV= 8.6%, 11.9% and 9.2% for N, RBD and Spike). is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint Data were analysed using Prism 8 (GraphPad) or R (version 4.0.2) within R Studio (version 1.2.5042) using packages rstatix (v 0.6.0; Kassambara, 2020) and the tidyverse (v1.3.0; Wickham, 2019). Locally estimated scatterplot smoothing (LOESS) with a span of 0.95 were used to visualise trends in antibody levels over time. Kruskal-Wallis followed by pairwise Dunn tests were performed to calculate statistical significance. Bonferroni method was used to correct for multiple comparisons. A P-value of ≤ 0.05 was considered statistically significant. The following models were applied to NAb data. Exponential decay in which an exponential decay curve was fitted to the full dataset where is measured in days. This method will likely underestimate the decay rate if the period in which antibodies increase after infection (growth period) for some individuals is long. It is also strongly affected by potential differences in both the size and timing of the peak level of NAbs between individuals. The best fit coefficients with 95% confidence interval are = 66 (56.7, 75.3) and = 0.0011 (0.00005, 0.00217). Growth and decay, which uses a two-part model fitted to the full dataset This allows an initial growth period to account for the underestimate that results from fitting a decay curve to measurements still in the growth phase. However, strong differences between individuals, particular in the timing of the peak level of NAbs will still affect estimations. This model predicts a slightly faster decay speed than the exponential decay model = 0.0016 (0.00053, 0.00274) that starts after 11 days ( = 11 (5.65, 16.35)). . CC-BY-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint 3. Of the 189 data points for NAbs, 127 are multiple measurements from 50 individuals, of which 33 (88 measurements) are in the exponential decay stage, i.e. inhibition is decreasing with time. The individual decay model incorporates data from these 33 individuals only. An exponential decay curve is fitted to each individual separately. An exponential distribution is fitted to the individual decay rates and the best fit parameter and 95% confidence interval reported. The mean decay rate much is faster than predicted by the other models with =0.00476 (95% confidence interval 0.00348, 0.00691). This work was funded by the School of Medicine Foundation (University of Auckland) and the Linfa Wang at Duke-NUS and Genscript for providing sVNT testing kits and technical advice. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint A systematic review of antibody mediated immunity to coronaviruses: kinetics, correlates of protection, and association with severity Prospects from systems serology research Profile of specific antibodies to SARS-CoV-2: The first report Persistence and decay of human antibody responses to the receptor binding domain of SARS-CoV-2 spike protein in COVID-19 patients Humoral Responses against SARS-CoV-2 Spike Defining the features and duration of antibody responses to SARS-CoV-2 infection associated with disease severity and outcome Antibodies in Persons with Mild Covid-19 Seroprevalence of anti-SARS CoV-2 antibodies in COVID-19 patients and healthy volunteers up to 6 months post disease onset Analysis of humoral immune responses in SARS-CoV-2 infected patients A SARS-CoV-2 surrogate virus neutralization test based on antibody-mediated blockage of ACE2-spike protein-protein interaction Neutralizing Antibodies Correlate with Protection from SARS-CoV-2 in Humans during a Fishery Vessel Outbreak with a High Attack Rate Antibodies to SARS-CoV-2 are associated with protection against reinfection Distinct systems serology features in children, elderly and COVID patients Distinct Early Serological Signatures Track with SARS-CoV-2 Survival Evolution of immunity to SARS-CoV-2 Successful Elimination of Covid-19 Transmission in New Zealand COVID-19 in New Zealand and the impact of the national response: a descriptive epidemiological study Collaborative networks enable the rapid establishment of serological assays for SARS-CoV-2 during nationwide lockdown in New Zealand SARS-CoV-2 antibodies in the Southern Region of New Zealand Mapping neutralizing and immunodominant sites on the SARS-CoV-2 spike receptor-binding domain by structure-guided high-resolution serology Standardization of enzyme-linked immunosorbent assays for serosurveys of the SARS-CoV-2 pandemic using clinical and athome blood sampling Immunological memory to SARS-CoV-2 assessed for greater than six months after infection Cross-reactive memory T cells and herd immunity to SARS-CoV-2 DNA vaccine protection against SARS-CoV-2 in rhesus macaques Toward superhuman SARS-CoV-2 immunity is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprintThe copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint Available from: http://www.nature.com/articles/s41591-020-01180-x . CC-BY-ND 4.0 International license It is made available under a perpetuity.is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprintThe copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprintThe copyright holder for this this version posted December 11, 2020. ; https://doi.org/10.1101/2020.12.10.20246751 doi: medRxiv preprint