medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Longitudinal variation in SARS-CoV-2 antibody levels and emergence of viral variants: implications for the ability of serological assays to predict immunity.

*Frauke Muecksch, Post-doctoral scientist1 *Helen Wise, Clinical Scientist2 Kate Templeton, Consultant Clinical Scientist2 Becky Batchelor, Clinical Scientist3 Maria Squires, Clinical Scientist2 Kirsty McCance, Chemical Pathology trainee2 Lisa Jarvis, Consultant Clinical Scientist4 Kristen Malloy, Clinical Scientist4 Elizabeth Furrie, Consultant Immunologist5, Claire Richardson, Clinical Scientist6 Jacqueline MacGuire, Consultant Clinical Scientist6 Ian Godber, Consultant Clinical Scientist7 Alana Burns, Clinical Scientist7 Sally Mavin, Director8 Fengwen Zhang, Post-doctoral scientist1 Fabian Schmidt, Post-doctoral scientist1 ^, Professor1,9 ^Sara Jenks, Consultant Biochemist2 ^Theodora Hatziioannou, Associate professor1

1Laboratory of Retrovirology, The , 1230 York Avenue, New York NY 10065 2Royal Infirmary of Edinburgh, NHS Lothian, 51 Little France Crescent, Edinburgh EH16 4SA 3Western General Hospital, Crewe Road South, Edinburgh, EH4 2XU 4SNBTS Microbiology Reference Laboratory, The Jack Copland Centre, 52 Research Avenue North, Heriot Watt Research Park, Edinburgh, EH14 4AP 5Ninewells Hospital and Medical School, NHS Tayside, Dundee, DD1 9SY 6University Hospital Monklands, NHS Lanarkshire, Airdrie, ML6 0JS 7Queen Elizabeth University Hospital, Govan Road, Glasgow, G51 4TF 8Scottish Microbiology Reference Laboratory, NHS Highland, Inverness, IV2 3UJ 9Howard Hughes Medical Institute, The Rockefeller University, 1230 York Avenue, New York NY 10065

* FM and HW contributed equally to this study

^Correspondence: Theodora Hatziioannou, Laboratory of Retrovirology, The Rockefeller University, 1230 York Avenue, NY10065, USA (+US code)-212-327-8040, [email protected], Sara Jenks, Royal Infirmary of Edinburgh, NHS Lothian, 51 Little France Crescent, Edinburgh EH16 4SA, UK, (+UK code)-01312426853, [email protected], Paul Bieniasz, Laboratory of Retrovirology, The Rockefeller University, 1230 York Avenue, NY10065, USA (+US code)-212-327-8848, [email protected].

NOTE: This preprint reports new research that has not been certified by peer review and should not be used to guide clinical practice. 1 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Abstract Background Serological assays are being deployed to monitor antibody responses in SARS- CoV-2 convalescents and vaccine recipients. There is a need to determine whether such assays can predict immunity, as antibody levels wane and viral variants emerge. Methods We measured antibodies in a cohort of SARS-CoV-2 infected patients using several high-throughput serological tests and functional neutralization assays. The effects of time and spike protein sequence variation on the performance and predictive value of the various assays was assessed. Findings Neutralizing antibody titers decreased over the first few months post-infection but stabilized thereafter, at about 30% of the level observed shortly after infection. Serological assays commonly used to measure antibodies against SARS-CoV-2 displayed a range of sensitivities that declined to varying extents over time. Quantitative measurements generated by serological assays based on the spike protein were better at predicting neutralizing antibody titers than assays based on nucleocapsid, but performance was variable and manufacturer positivity thresholds were not able to predict the presence or absence of detectable neutralizing activity. Even though there was some deterioration in correlation between serological measurements and functional neutralization activity, some assays maintained an ability to predict neutralizing titers, even against variants of concern. Interpretation The ability of high throughput serological assays to predict neutralizing antibody titers is likely crucial for evaluation of immunity at the population scale. These data will facilitate the selection of the most suitable assays as surrogates of functional neutralizing activity and suggest that such measurements may have utility in clinical practice.

Introduction The world has experienced an unprecedented pandemic following the emergence of severe acute respiratory coronavirus 2 (SARS-CoV-2). Millions have died and the repercussions have affected every aspect of life. The remarkable mobilization of the scientific community in response to the pandemic has led to the rapid development of safe and effective vaccines, as well as reagents and assays to aid in the detection and mitigation of virus spread. An early prominent issue in the pandemic was the accurate identification of infected individuals at a large scale. While PCR-based assays remain a reliable and sensitive test for infection they are not amenable to mass population screening. Thus, serological assays, despite limitations,1-3 have been instrumental for surveillance and providing selection criteria to recruit vaccine trial participants and convalescent plasma donors. Monitoring antibody titers is necessary to measure the magnitude and longevity of immune responses induced by natural infection or vaccination. As immune responses to SARS-CoV-2 antigens are increasingly elicited by infection and/or vaccination, the measurement of antibody titers and the ability of such measurements to predict protection from infection or disease will be of great importance. Whether simple serological tests will be able to predict neutralizing antibody titers and immunity to SARS-CoV-2 is yet to be determined. Moreover, as antibodies both mature, acquiring greater affinity, while total levels decline, 4-6 and new SARS-CoV-2 variants emerge, the predictive value of serological tests based on the prototype viral strain will need to be evaluated.

2 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

A number of high-throughput serological assays are routinely used to detect antibodies against nucleocapsid (N) or spike (S) viral antigens. These assays were initially designed to provide a positive or negative test result, but they also generate quantitative measurements of antibody levels. Prior studies of how these quantitative serological values correlate with neutralizing antibodies titers have yielded variable results7-9. Importantly, the sensitivity of the assays as diagnostic tools and their predictive value for immune parameters several months after infection and against variants of concern has not been assessed. We present a longitudinal study of COVID-19 recovered patients over 6 months post-infection, evaluating the diagnostic sensitivity of ten different serological assays and their ability to predict neutralizing antibody activity against SARS-CoV-2.

Results A previously reported8 cohort of participants that developed mild symptoms following SARS- CoV-2 infection, was repeatedly sampled, up to 7.2 months post-infection, to evaluate how neutralizing antibody levels correlate over time with antibodies measured using ten high throughput serological assays. Neutralizing antibody titers were measured at up to five visits for each participant using a pseudovirus neutralization assay that correlates well with neutralization against authentic SARS-CoV-210. Consistent with prior studies4,8,11,12, the half-maximal neutralizing titers (NT50) declined over time in the majority of patients (Figure 1A,B). The most significant rate of decrease, approximately 25%, was observed between the early visits reaching a 45% decrease in NT50 by visit 3, approximately 70 days post-infection (Figure 1C). Thereafter the rate of decrease became less pronounced and NT50 values at visits 4 and 5, at approximately 3-7 months post-infection, appeared to stabilize at ~30% of the levels observed within the first 2 months post-infection (Figure 1B,C). NT50 values were higher in male than female participants at visit 1, but the overall rate of NT50 decline from the first to the last visit was greater for male participants (supplementary Figure 1). Thus, the sex difference gradually diminished and was not discernable by visit 4. No correlation between NT50 values and age were observed at any time point (supplementary Figure 1). The same sera were analyzed using nine different serological assays that detect antibodies against either the viral nucleocapsid (N) protein, or various forms of the spike (S) protein that included trimeric S, S1/S2 subunits or the receptor binding domain (RBD). First, the sensitivity of each assay was determined for three time windows over the course of the study. All assays were sensitive at the first window, 21-80 days post-infection; Abbott IgGII Quant, Roche S and Roche N had the highest sensitivities at 100%, Siemens COV2T and Diasorin Trimeric S were 95% sensitive and Diasorin S1/S2, Euroimmun and Abbott (N) ranged from 85-90%. While Abbott IgGII Quant, Roche N and Diasorin S1/S2 maintained their sensitivity over time, the sensitivity of other assays declined to varying degrees, ranging from 45% to 85% at >140 days post-infection (Figure 2A). Thus, the performance of the assays for serosurveillance applications at >140 days after infection was extremely variable. In addition to indicating whether a serum sample is negative or positive for antibodies against viral antigens, each assay indicates quantitative antibody levels within assay-specific scales. Analysis of antibody levels over time showed assay-dependent differences in trajectory that were not dependent on whether nucleocapsid or spike antigens were used (Figure 2A-J). Median antibody levels measured by the Roche S and Siemens COV2T assays increased slightly over time, those measured by the Roche N and Diasorin S1/S2 assays remained approximately

3 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

constant, while levels measured by the Siemens sCOVG, Diasorin Trimeric S, Euroimmun and both Abbott assays decreased over time. The deviation of individual participant antibody levels from the mean increased over time in all assays except Abbott IgGII Quant which exhibited high deviations from the first time point. We determined the ability of each serological assay to predict pseudotype virus neutralizing antibody titers over time (Figure 3, supplementary Figure 2). For each assay, the correlation with NT50 values was closest at early time points and deteriorated over time with the poorest correlation observed at visit 5 in each case. An additional assay, that measures antibodies that block the interaction between the RBD and the virus receptor (cPass), was included in this analysis using samples only from visits 1 and 5. At visit 1, antibody levels measured using the Diasorin S1/2 and trimeric S assays had the highest correlation with NT50 followed by Abbott IgGII, Siemens sCOVG, Euroimmun and cPass (Figure 3A, supplementary Figure 2). Antibody levels measured with the remaining spike antigen-based assays (Roche S and Siemens COV2T) which are designed to detect total antibody levels against the spike antigen regardless of antibody class had a lower correlation with NT50 titers, whereas the nucleocapsid based assays had the poorest correlation. For all serological assays, the correlation with neutralizing antibody titers decreased over time, but for all of the spike antigen based IgG assays (Diasorin S1/2 and trimeric S assays, Abbott IgGII, Siemens sCOVG) and cPass the correlation coefficient r, remained >0.75 even at visit 5 (Figure 3A, supplementary Figure 2). In general, the decrease in neutralizing antibody titers over time was proportionately greater than the corresponding decrease in levels measured using serological assays (Figure 3B), particularly between early time points (supplementary Figure 2). Thus, declines in antibody levels over time measured using serological assays did not, in some cases, accurately reflect the decrease in neutralization activity. This was particularly the case for the nucleocapsid based assays, where the magnitude of the decrease in antibody measurements did not correlate with the decrease in NT50. Nevertheless, for assays that correlated best with NT50 titers at early time points, declining levels of antibodies measured in the serological assays predicted declining NT50 quite well (Figure 3B, supplementary Figure 2). The ability of serological assays to qualitatively identify sera that did or did not have detectable neutralizing activity, was estimated (Table 1, supplementary Figure 3). This comparison was performed on sera collected across all time points and selected ‘cut off’ values for each serological assay scale were evaluated for sensitivity/specificity and predictive value (Table 1). None of the assays were effective in qualitatively discriminating neutralizing vs non-neutralizing sera using manufacturer recommended cut-offs, with specificity ranging from 3% for the Roche S and N total antibody assays to 72% for the EuroImmun assay. By selecting different cut off values the sensitivity, specificity and predictive values could be improved for some assays (Table 1). The occurrence of SARS-CoV-2 variants compromises the ability of first wave convalescent or vaccinee antibodies to neutralize contemporaneous viruses and might further erode the ability of serological assays based on proteins derived from the prototype (Wuhan-hu-1) virus to predict neutralizing antibody titers. We determined the ability of 58 serum samples obtained at visits 1 and 5 to neutralize selected variants of concern of the B.1.1.7, B.1.351, B.1.617.2 lineages that have been associated with partial resistance to neutralization13-16. Sera from both visits were able to neutralize viruses carrying the B.1.1.7 spike with potencies comparable to those observed with the Wuhan-hu-1 spike protein (Figure 4A, supplementary Figure 4). In contrast, titers against viruses with the B.1.1.7 (E484K) amino acid substitution, B.1.617.2 or either of the two B.1.351

4 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

spike variants were decreased by approximately 2.5- to-5 fold at visit 1 (Figure 4B-E, supplementary Figure 4). Differences between the ability to neutralize Wuhan-hu-1 and variant spike bearing viruses became somewhat less pronounced at visit 5 (Figure 4A-). Thus, the inclusion of amino acid substitution E484K in the context of B.1.1.7 or substitutions found in the spikes from other variants of concern appeared to significantly reduce neutralization titers, with the largest effect seen with substitutions found in the B1.351 variants, consistent with prior reports13,14,16-20. Overall, the correlation of antibody levels, as measured by serological assays, with neutralizing activity was marginally weaker for variant pseudotyped viruses than for the Wuhan-hu-1 variant (Figure 4, supplementary Figure 5). The weakest correlations were observed for the B.1.617 variant. These data suggest that some degree of variation in viruses circulating in a given population does not abolish the ability of serological assays to predict neutralizing potency against a given variant. Nevertheless, the calibration of predictions made using serological assays needs to account for which particular SARS-CoV-2 variants are in circulation in a given locale, and their susceptibility to neutralization.

Methods Participants 112 Participants with a history of SARS-CoV-2 infection, diagnosed by RT-PCR were recruited8. Recruits were surveyed to determine the date of the positive PCR test, the date of onset of symptoms, and if their symptoms required hospitalization. Serum samples were taken at a baseline visit (~3.5 to ~8.5 weeks post PCR test), and 2 weeks (visit 2), 4 weeks (visit 3), 8 weeks (visit 4) and 22 weeks (visit 5) later. In total, 101 participants, completed at least 3 visits, 58 participants completed the fifth visit and were included in the neutralization assays. Four patients out of 101 were hospitalized but none required intensive care. The mean age of the participants was 44.5 years (21 – 65 y), with 72 female (71% of cohort) participants. At visit 1 (baseline), the average number of days between PCR test and visit 1 (baseline) was 41.05 days (24 – 64 days); at visit 2 (2 weeks post-baseline), the average number of days post-PCR test was 55.1 days (40 – 79 days); at visit 3 (4 weeks post-baseline), the average number of days post- PCR test was 70.06 days (55 – 95 days); at visit 4 (8 weeks post-baseline), the average number of days post-PCR test was 97.81 days (85 – 110 days); at visit 5 (28 weeks post-baseline), the average number of days post-PCR test was 194 days (160-216 days). Ethical approval was obtained for this study to be carried out through the NHS Lothian BioResource (SR1407) and London-Brent Research Ethics Committee (REC ref: 20/HRA/3764 IRAS:28653). All recruits gave written and informed consent for serial blood sample collection. De-identified samples were shipped to the Rockefeller University whose IRB reviewed and approved the study.

Serological assays The Abbott SARS-CoV-2 IgG (NHS Lothian) and SARS-CoV-2 IgGII (NHS Greater Glasgow and Clyde) assays are two-step chemiluminescent microparticle immunoassays (CMIA) designed to detect IgG antibodies against nucleoprotein and RBD respectively. The DiaSorin LIASON SARS-CoV-2 S1/2 IgG (NHS Lothian), DiaSorin LIASON trimericS IgG (NHS Highland) and Siemens sCOVG RBD IgG (NHS Tayside) assays are also two-step CMIA designed to detect IgG antibodies. The Euroimmun SARS-CoV-2 IgG assay (Scottish National Blood Transfusion

5 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Service) is an indirect ELISA using the S1 domain of Spike as the antigen. The Roche Anti- SARS-CoV-2 N and S assays (NHS Lanarkshire) and the Siemens COV2T assay (NHS Tayside) are two-step bridging electrochemiluminesent immunoassays (ECLIA) using nucleocapsid or the RBD of the Spike protein as antigens. The cPass assay detects antibodies that block binding of a soluble RBD to an immobilized cellular receptor protein.

SARS-CoV-2 pseudotyped reporter virus SARS-CoV-2 pseudotyped particles were generated as previously described10. Briefly, 293Tcells were transfected with pNL4-3ΔEnv-nanoluc and pSARS-CoV-2-SΔ19. 48 hours later particles were harvested, filtered and stored at -80°C. The amino acid deletions and/or substitutions corresponding to variants of concern were incorporated into a spike expression plasmid using synthetic gene fragments (IDT) or overlap extension PCR mediated mutagenesis and Gibson assembly. Specifically, the variant-specific deletions and substitutions introduced were: B.1.1.7: ΔH69/V70, ΔY144, N501Y, A570D, D614G, P681H, T761I, S982A, D1118H B.1.1.7 E484K: ΔH69/V70, ΔY144, N501Y, A570D, E484K, D614G, P681H, T761I, S982A, D1118H B.1.351 #1: L18F, D80A, D215G, ∆242-4, K417N, E484K, N501Y, D614G, A701V B.1.351 #2: D80A, D215G, L242H, R246I, K417N, E484K, N501Y, D614G, A701V B.1.617.2: T19R, ∆156-8, L452R, T478K, D614G, P681R, D950N These spike proteins included the R683G substitution, which disrupts the furin cleavage site and increases particle infectivity and neutralization sensitivity. Therefore, in these neutralization assays a wildtype SARS-CoV-2 spike (NC_045512), carrying R683G was used for comparative purposes.

Pseudotyped virus neutralization assay Fivefold serially diluted serum from COVID-19-convalescent individuals were incubated with SARS-CoV-2 pseudotyped virus for 1 h at 37 °C. The mixture was subsequently added to 293TAce2 cl22 cells (for analyses using SARS-CoV-2 Wuhan-Hu-1 pseudovirus) or HT1080Ace2 cl14 cells (for analyses involving variant pseudovirus panels and respective Wuhan-Hu-1 R683G controls).10 The starting serum dilution on cells was 1:50. Nanoluc Luciferase activity in lysates was measured 48 hours post-inoculation using the Nano-Glo Luciferase Assay System (Promega) with the Glomax Navigator (Promega). Relative luminescence units were normalized to those derived from cells infected with SARS-CoV-2 pseudotyped virus in the absence of serum. The half-maximal neutralization titers for sera (NT50) were determined using four-parameter nonlinear regression (least squares regression method without weighting; constraints: top=1, bottom=0) (GraphPad Prism).

Role of funding source The funders of the study had no role in the design, collection, analysis and interpretation of the data presented.

Discussion Tracking transmission dynamics, spread and prevalence of viral infections, is critical in mitigating viral epidemics, particularly when a large number of cases remain asymptomatic

6 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

during infection as is the case with SARS-CoV-2. Moreover, the wide-spread use of vaccines necessitates the accurate determination of the vaccination and immune status of individuals. High-throughput serological assays address these needs but their usefulness obviously depends on their accuracy and reliability. In this study we compared the results provided by a number of SARS-CoV-2 serological assays with an emphasis on their ability to predict neutralization activity as antibodies both wane and evolve, and SARS-CoV-2 variants emerge4,5. Serological assays are generally optimized for increased sensitivity so they can reliably diagnose the presence or absence of antibodies against viral antigens1,21. The majority of the assays used herein accomplish this goal with sera obtained shortly after infection, however, their sensitivity was not always maintained over time and certain assays exhibited a sharp decline in sensitivity at later time points after infection. This loss of sensitivity was not related to the antigen on which the assays were based. Serological assays also provide a quantitative result that could potentially enable their use for estimation of antibody levels and prediction of immunity, especially if antibody levels correlate with functional neutralizing antibody titers. Assays that detect spike-binding antibodies can use various protein subdomains or conformations (e.g. isolated RBD, S1 or a stabilized trimeric spike) as their antigens. Moreover, some assays detect only specific antibody classes such as IgG, or those that directly interfere with RBD-receptor binding. In contrast, neutralization assays detect all antibodies capable of inhibiting spike-mediated virus entry into cells. While neutralizing antibodies are sometimes dominated by those targeting RBD22,23, including those that block ACE2 binding, antibodies targeting the N-terminal domain of S1 can significantly contribute to the overall serum neutralization activity in plasma24-26. Weak neutralizing activity has also been ascribed to antibodies targeting the region of S2 involved in fusion27. Multiple studies have shown that neutralizing antibody titers following natural infection or vaccination wane over time4,8,11,12,28, a decline that is not always accurately reflected by serological assays. Nevertheless, certain assays used herein that detect spike-specific antibodies maintained good levels of correlation with neutralizing titers over time. Assays measuring spike- specific IgG antibodies predicted neutralizing antibody titers more accurately than those measuring total antibodies against spike or those against the nucleocapsid protein. The Diasorin assays, Abbott IgGII, Siemens sCOVG, Euroimmun and cPass had the highest correlation with neutralizing antibody titers across all comparisons and changes in quantitative values over time for these assays were most closely correlated with changes in neutralizing antibody levels within individuals. The quantitative results from these assays are therefore best suited for estimating neutralizing antibody levels at a population level. In contrast, the qualitative assay results which, if based solely on the manufacturer recommended cut-offs, are poorly specific for detecting the presence of neutralizing antibodies and may lead to a significant over-estimation of antibody related immunity. For some assays it may be possible to improve specificity for the presence of neutralizing antibodies through selecting a higher quantitative cut-off value, thus improving the positive predictive value of these assays for neutralizing antibody detection. The majority of the naturally infected population studied thus far were infected with viral variants closely related to a prototype variant (Wuhan hu-1). The antigens used by all serological assays rely on protein sequences derived from that prototype. However, over the last several months new variants have emerged that encode multiple amino acid substitutions in their spike proteins, some of which impact neutralization by convalescent or vaccinee antibodies13,14,16,17,29,30. Our data indicate that antibody levels measured using several serological assays maintain good correlation with neutralization titers against some of the most important

7 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

variants that have emerged thus far. However, this property will need to be monitored in the future if serological assays are to be used to predict immunity, particularly as antibodies diversify in response to variant virus infection and, potentially, variant booster vaccination. The need for serological assays in monitoring natural infection at the population level remains. Moreover, the introduction of vaccines raises new requirements for serological assays. These requirements include (i) distinction between vaccinated and naturally infected individuals, (ii) prognostication of levels of protection against infection and disease afforded by vaccines and (iii) identification of individuals where a boosting immunization or monoclonal antibody therapy is indicated. Furthermore, in instances where countries are considering deploying ‘immunity passports’ to allow, for example travel, the selection of assays and the diagnostic cut off used can have important ramifications. Finally, since measurement of antibody function (e.g., neutralization) at a population level is not currently practical, establishing the ability of serological assays to predict neutralization and immunity will help determine correlates of protection that can be applied at a large scale, and perhaps as part of routine clinical practice.

Contributors HW, SJ, TH and PDB conceived and designed the study. FM performed and analyzed the neutralization and cPASS assays and performed cross-assay data analyses and data visualization. HWKT and SJ acquired and analyzed data using the serological assay platforms. HW and SJ performed ROC and sensitivity analyses and respective visualization. HW, SJ and MS provided Abbott N data. EF provided Siemens data. SM provided DiaSorin data. IG and AB provided Abbott S IgG II data. LJ and KM provided EuroImmun data. CR and JM provided Roche data. BB and KM assisted with collection of research clinic samples. FZ and FS provided spike plasmids. TH wrote the manuscript with help from FM, PDB, SJ and input from all authors. All authors had access to all the data in the study and accept responsibility for publication.

Declaration of interests SJ & EF received honoraria from Siemens for an online webinar in October 2020 which was paid to their institutions. EF received sCOVG reagent from Siemens to support this study.

Acknowledgements We acknowledge the support of NRS BioResource and NHS Lothian Outpatients with the provision of this research study. We are grateful to Dr Linfa Wang and GenScript for providing the cPass assay kits. Funding sources NIH (R01AI78788 to TH, R01 AI50111 to PDB) and NRS BioResource (SR1407 to SJ).

Data sharing De-identified data will be available for share following manuscript publication.

References

1. GeurtsvanKessel CH, Okba NMA, Igloi Z, et al. An evaluation of COVID-19 serological assays informs future diagnostics and exposure assessment. Nat Commun 2020; 11(1): 3436.

8 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

2. La Marca A, Capuzzo M, Paglia T, Roli L, Trenti T, Nelson SM. Testing for SARS-CoV- 2 (COVID-19): a systematic review and clinical guide to molecular and serological in-vitro diagnostic assays. Reprod Biomed Online 2020; 41(3): 483-99. 3. Ozcurumez MK, Ambrosch A, Frey O, et al. SARS-CoV-2 antibody testing-questions to be asked. J Allergy Clin Immunol 2020; 146(1): 35-43. 4. Gaebler C, Wang Z, Lorenzi JCC, et al. Evolution of antibody immunity to SARS-CoV- 2. Nature 2021. 5. Muecksch F, Weisblum Y, Barnes CO, et al. Development of potency, breadth and resilience to viral escape mutations in SARS-CoV-2 neutralizing antibodies. bioRxiv : the preprint server for biology 2021. 6. Wang Z, Muecksch F, Schaefer-Babajew D, et al. Naturally enhanced neutralizing breadth to SARS-CoV-2 one year after infection. Nature 2021. 7. Luchsinger LL, Ransegnola BP, Jin DK, et al. Serological Assays Estimate Highly Variable SARS-CoV-2 Neutralizing Antibody Activity in Recovered COVID-19 Patients. Journal of clinical microbiology 2020; 58(12). 8. Muecksch F, Wise H, Batchelor B, et al. Longitudinal analysis of serology and neutralizing antibody levels in COVID19 convalescents. The Journal of infectious diseases 2020. 9. Tang MS, Case JB, Franks CE, et al. Association between SARS-CoV-2 neutralizing antibodies and commercial serological assays. Clin Chem 2020. 10. Schmidt F, Weisblum Y, Muecksch F, et al. Measuring SARS-CoV-2 neutralizing antibody activity using pseudotyped and chimeric viruses. The Journal of experimental medicine 2020; 217(11). 11. Crawford KHD, Dingens AS, Eguia R, et al. Dynamics of Neutralizing Antibody Titers in the Months After Severe Acute Respiratory Syndrome Coronavirus 2 Infection. The Journal of infectious diseases 2021; 223(2): 197-205. 12. Seow J, Graham C, Merrick B, et al. Longitudinal observation and decline of neutralizing antibody responses in the three months following SARS-CoV-2 infection in humans. Nature microbiology 2020; 5(12): 1598-607. 13. Wang Z, Schmidt F, Weisblum Y, et al. mRNA vaccine-elicited antibodies to SARS- CoV-2 and circulating variants. Nature 2021. 14. Wibmer CK, Ayres F, Hermanus T, et al. SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma. bioRxiv : the preprint server for biology 2021: 2021.01.18.427166. 15. Weisblum Y, Schmidt F, Zhang F, et al. Escape from neutralizing antibodies by SARS- CoV-2 spike protein variants. eLife 2020; 9. 16. Wibmer CK, Ayres F, Hermanus T, et al. SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma. Nature medicine 2021; 27(4): 622-5. 17. Wang Z, Schmidt F, Weisblum Y, et al. mRNA vaccine-elicited antibodies to SARS- CoV-2 and circulating variants. Nature 2021; 592(7855): 616-22. 18. Tada T, Dcosta BM, Samanovic-Golden M, et al. Neutralization of viruses with European, South African, and United States SARS-CoV-2 variant spike proteins by convalescent sera and BNT162b2 mRNA vaccine-elicited antibodies. bioRxiv : the preprint server for biology 2021. 19. Tan Y, Liu F, Xu X, et al. Durability of neutralizing antibodies and T-cell response post SARS-CoV-2 infection. Front Med 2020; 14(6): 746-51.

9 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

20. Beaudoin-Bussieres G, Laumaea A, Anand SP, et al. Decline of Humoral Responses against SARS-CoV-2 Spike in Convalescent Individuals. mBio 2020; 11(5). 21. Galipeau Y, Greig M, Liu G, Driedger M, Langlois MA. Humoral Responses and Serological Assays in SARS-CoV-2 Infections. Frontiers in immunology 2020; 11: 610688. 22. Greaney AJ, Loes AN, Crawford KHD, et al. Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies. Cell host & microbe 2021. 23. Rogers TF, Zhao F, Huang D, et al. Isolation of potent SARS-CoV-2 neutralizing antibodies and protection from disease in a small animal model. Science (New York, NY) 2020; 369: 956-63. 24. Liu L, Wang P, Nair MS, et al. Potent neutralizing antibodies against multiple epitopes on SARS-CoV-2 spike. Nature 2020; 584(7821): 450-6. 25. McCallum M, De Marco A, Lempp FA, et al. N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2. Cell 2021; 184(9): 2332-47 e16. 26. Cerutti G, Guo Y, Zhou T, et al. Potent SARS-CoV-2 neutralizing antibodies directed against spike N-terminal domain target a single supersite. Cell host & microbe 2021; 29(5): 819- 33 e7. 27. Sauer MM, Tortorici MA, Park YJ, et al. Structural basis for broad coronavirus neutralization. Nature structural & molecular biology 2021; 28(6): 478-86. 28. Widge AT, Rouphael NG, Jackson LA, et al. Durability of Responses after SARS-CoV-2 mRNA-1273 Vaccination. The New England journal of medicine 2021; 384(1): 80-2. 29. Harvey WT, Carabelli AM, Jackson B, et al. SARS-CoV-2 variants, spike mutations and immune escape. Nat Rev Microbiol 2021. 30. Abdool Karim SS, de Oliveira T. New SARS-CoV-2 Variants - Clinical, Public Health, and Vaccine Implications. The New England journal of medicine 2021; 384(19): 1866-8.

Figure Legends

Figure 1 – Neutralization activity in longitudinal COVID-19 patient sera. (A) Half-maximal neutralization titers (NT50s) for each sample collected at the visit indicated (v1-5). (B) Relative NT50 values in sera obtained at visit 1-5, normalized to visit 1. Colored horizontal bars indicate median values with 95% confidence interval. Statistical significance was determined with Kruskal-Wallis test and Dunn’s multiple comparisons test. (C) Relative decay of NT50 per day between visits. Red horizontal bars indicate median. Statistical significance was assessed with the Wilcoxon test.

Figure 2 - Serological analysis of longitudinal COVID-19 patient sera. (A) Sensitivity of the indicated serological assays in samples collected at three different time-intervals post PCR, as indicated. Mean with 95% confidence intervals are shown. (B-J) Relative serological results at visit 1-5 (v1-5), normalized to visit 1 (left panels) as well as serological results per participant over time with each line representing a single participant (right panels) for the assay indicated. Horizontal bars in B-J indicate median with 95% confidence interval and statistical significance was determined with the Kruskal-Wallis test and subsequent Dunn’s multiple comparisons test. Dotted lines indicate serological assay thresholds.

10 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Figure 3 – Correlation of neutralization titers and serology assays. (A) Correlation of NT50 (x-axis) with serological assay values (y-axis) obtained at visit 1 (black) and visit 5 (red) for each participant. Statistical significance was determined using the Spearman correlation test for samples obtained at visit 1 and visit 5 independently, as indicated. Dotted lines indicate serological assay thresholds. (B) Correlation of fold-change (visit1 to visit5) of NT50s with corresponding fold-change in serological assay values for indicated serology assays. Statistical significance was determined using the Spearman correlation test. Dotted lines at x=1 and y=1 indicate unchanged assay results over time.

Figure 4 – Neutralization of variants of concern. NT50, normalized to wt of pseudoviruses with spikes from variants(A) B.1.1.7, (B) .1.1.7 (E484K) (C) B.1.351 (#1) (D) B.1.351 (#2) or (E) B.1.617.2 measured for visits indicated. Statistical significance was determined using Wilcoxon test. Dotted line indicates 100% (equal NT50). (F) Correlation coefficients of NT50 values between original and variant viruses and each serological assay as indicated.

Supplementary Figure 1 – Correlation of neutralization activity and demographic parameters. (A) Correlation of NT50s in samples collected at different timepoints. (B) NT50 values per patient over time with each line representing a single participant. (C) NT50 values in female and male participants, sampled at the indicated visit. (D) Relative NT50 values at visit 5, normalized to visit 1 for female and male participants. Statistical significance was assessed with the Mann-Whitney test. (E) Correlation of NT50 and age at the indicated timepoints. Statistical significance in (A) and (E) was determined using the Spearman correlation test. Dotted lines in (A) and (B) indicate limit of detection.

Supplementary Figure 2 – Correlation of neutralization titers and serology assays per visit. (A-J) Correlation of NT50 (x-axis) and indicated serological assay measurements (y axis) at visit 1 through 5. Sampling timepoints are indicated by color and shown collectively (left) and individually. Statistical significance was determined using the Spearman correlation test and r and P-values are indicated for total samples (left panels) and for the individual visits. Dotted lines in left panels indicate serological assay thresholds.

Supplementary Figure 3 – ROC analysis. Receiver-operating characteristic curve (ROC) for prediction of NT50>50 based on results obtained with indicated serology assays. Shown are mean values.

Supplementary Figure 4 – Neutralization of variants of concern. NT50s for wt (A) B.1.351 (#1) (B), B.1.351 (#2) (C), B.1.1.7 (D) and B.1.1.7 (E484K) (E) pseudovirus at visit 1 and visit 5. Statistical significance was determined using Wilcoxon test. Dotted line indicates limit of detection. Deletions/substitutions present in VOCs, as well as respective wt control are in R683G background, as indicated.

Supplementary Figure 5 – Neutralization of variants of concern versus serology assays. (A- K) Correlation of, from left to right NT50s against the pseudoviruses indicated with the Abbott IgGII Quant (A), Roche S (B), Roche N (C), Siemens COV2T (D), Diasorin S1/2 (E), Diasorin Trimeric S (F), Siemens sCOVG (G), Euroimmun (H), Abbott NC (I) and the cPass assay (J) as well as with wt NT50 values. Samples obtained at visit 1 and visit 5 are included. Statistical

11 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

significance was determined using the Spearman correlation. Deletions/substitutions present in VOCs, as well as respective wt control are in R683G background.

12 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. figure 1 It is made available under a CC-BY-NC-ND 4.0 International license . A B C 0.0040 >0.9999

<0.0001 0.2142 4000 v1 5.0 v2 600 <0.0001 400 3000 200 v3 150 2.5 125 NT50 2000 v4 100 0.0 v5 % NT50 75 1000 50 -2.5 25 %decay NT50 per day -5 0 0 -10 Participants v1 v2 v3 v4 v5 v1to3 v3to5 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. figure 2 It is made available under a CC-BY-NC-ND 4.0 International license .

0.0090 0.5431

A B 0.0016 >0.9999 100 175 100000 90 150 80 125 10000 70 Abbott IgGII Quant (RBD IgG) 100 Roche S (RBD total) 1000 60 Roche N (NC total) 75 Siemens COV2T (RBD total) Sensitivity/% 50 Diasorin S1/2 (S1/2 IgG) 50 Diasorin Trimeric S (Sstablized native IgG) 100 (RBD IgG) (AU/ml) Abbott IgGII Quant Abbott IgGII Quant Siemens sCOVG (RBD IgG) (RBD IgG) (% of v1) 40 Euroimmun (S1 IgG) 25 Abbott (NC IgG) 30 0 10 v1 v2 v3 v4 v5 50 100 150 200 21-80 (n = 305) 81-140 (n = 52) >140 (n = 79) days post PCR time post PCR/days C D 0.5004 0.8135 >0.9999 0.0222

0.0155 >0.9999 <0.0001 >0.9999 6000 10000 500 1000 4000 2000 1000 400 500 100 400 100 300 10 300 10 200 200 1 1 100 100 Roche N (NC total) (COI) Roche S (RBD total) (U/ml) Roche N (NC total) (% of v1)

Roche S (RBD total) (% of v1) 0 0.1 0 0.1 v1 v2 v3 v4 v5 50 100 150 200 v1 v2 v3 v4 v5 50 100 150 200 >0.9999 >0.9999 days post PCR days post PCR E F 0.6113 >0.9999 >0.9999 0.5898 >0.9999 0.5429 6000 10000 600 1000 3000 400 1000 200 400 150 100 100 125 300 100 10 200 75 10 Siemens COV2T Siemens COV2T Siemens COV2T Siemens COV2T

(RBD total) (COI) 50 100 1 (RBD total) (% of v1) 25 0 0.1 0 1 Diasorin S1/2 (S1/2 IgG) (AU/ml) v1 v2 v3 v4 v5 50 100 150 200 Diasorin S1/2 (S1/2 IgG) (% of v1) v1 v2 v3 v4 v5 50 100 150 200 days post PCR days post PCR 0.0150 >0.9999 0.0010 >0.9999 G H <0.0001 0.0010 <0.0001 0.0121 1000 1000 175 1000 500 150 150 125 100 100 125 100 IgG) (AU/ml) IgG) (% of v1) 100 10 75 75 10 50

(RBD IgG) (COI) 1 Siemens sCOVG

50 Siemens sCOVG Diasorin Trimeric S Trimeric Diasorin Diasorin Trimeric S Trimeric Diasorin (RBD IgG) (% of v1) stabilized native

stabilized native 25 25 (S (S 0 1 0 0.1 v1 v2 v3 v4 v5 50 100 150 200 v1 v2 v3 v4 v5 50 100 150 200 days post PCR days post PCR I 0.0068 0.2442 J 0.1414 <0.0001 <0.0001 0.0395

200 175 <0.0001 0.1551 100 175 10 150 150 125 10 125 100 100 1 1 75 75 50 50 0.1

25 25 Abbott (NC IgG) (s/c) Abbott (NC IgG) (% of v1) Euroimmun (S1 IgG) (S/C)

Euroimmun (S1 IgG) (% of v1) 0 0.1 0 0.01 v1 v2 v3 v4 v5 50 100 150 200 v1 v2 v3 v4 v5 50 100 150 200 days post PCR days post PCR medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. figure 3 It is made available under a CC-BY-NC-ND 4.0 International license .

A r=0.8977 r=0.9073 r=0.8254 r=0.6106 r=0.8543 v1 P<0.0001 v1 v1 v1 v1 P<0.0001 5 4 P<0.0001 3 P<0.0001 P<0.0001 10 r=0.7609 10 10 4 3 r=0.7356 v5 r=0.6855 r=0.3971 10 r=0.6149 10 v5 P<0.0001 v5 P<0.0001 v5 P=0.002 v5 P<0.0001 3 P<0.0001 4 10 3 10 102 10 2 2 10 10 102 103 101 1 1 Roche S 10 Roche N 10 1 2 0 10 Diasorin S1/2 10 0 (NC total) (COI) 10 (RBD total) (U/ml) Abbott IgGII Quant (RBD IgG) (AU/ml) Siemens COV2T Siemens COV2T 10 (RBD total) (COI) 0 10 (S1/2 IgG) (AU/ml)

1 -1 -1 10 10 10 10-1 100 r=0.9030 r=0.8893 r=0.8960 r=0.8988 v1 P<0.0001 r=0.6574 v1 v1 P<0.0001 v1 P<0.0001 v1 P<0.0001 P<0.0001 103 103 r=0.6835 15 r=0.4254 103 r=0.7814 r=0.7689 r=0.7516 v5 v5 v5 v5 P<0.0001 v5 P<0.0001 10 P<0.0001 P=0.0008716 P<0.0001 102 102 10 102

IgG) (AU/ml) 1 10 1 1 1 10 Euroimmun 5 10

0 (S1 IgG) (S/C) cPass (sVNT50) (RBD IgG) (COI) 10 Siemens sCOVG Diasorin trimeric S Abbott (NC IgG) (S/C) stabilized native 0 -1 0 (S 10 10 0.1 0 10 101 102 103 104 101 102 103 104 101 102 103 104 101 102 103 104 101 102 103 104 NT50 NT50 NT50 NT50 NT50

B 64 64 64 64 64

16 16 ) 16 16 16

4 4 total 4 4 4

NC

1 1 ( 1 1 1 N 0.25 0.25 0.25 0.25 0.25 (fold-change v1-5) (fold-change v1-5) (fold-change v1-5) (fold-change v1-5) (fold-change v1-5) Roche 0.0625 r=0.7638 Roche S (RBD total) 0.0625 r=0.5576 0.0625 r=-0.2361 0.0625 r=0.5455 0.0625 r=0.8015 P<0.0001 Diasorin S1/2 (S1/2 IgG) P<0.0001 P<0.0001 P=0.07440 (RBD total) siemens COV2T P<0.0001 Abbott IgG II Quant (RBD IgG) 0.015625 0.015625 0.015625 0.015625 0.015625

64 64 64 64 64 IgG) 16 16 16 16 16

stab.native 4 4 4 4 4

1 1 1 1 1

0.25 0.25 0.25 0.25 0.25 Abbott (NC IgG) cPass (sVNT50) (fold-change v1-5) (fold-change v1-5) (fold-change v1-5) (fold-change v1-5) (fold-change v1-5) 0.0625 r=0.7477 0.0625 r=0.7906 Euroimmun (S1 IgG) 0.0625 r=0.4501 0.0625 r=-0.1648 0.0625 r=0.6993

P<0.0001 Siemens sCOVG (RBD IgG) P<0.0001 P<0.0001 P=0.2164 P<0.0001

Diasorin trimeric S (S 0.015625 0.015625 0.015625 0.015625 0.015625 1 2 4 8 1 2 4 8 1 2 4 8 1 2 4 8 1 2 4 8 0.5 0.5 0.5 0.5 0.5 0.25 0.25 0.25 0.25 0.25 0.125 0.125 0.125 0.125 0.125 0.0625 0.0625 0.0625 0.0625 0.0625 0.03125 0.03125 0.03125 0.03125 0.03125 NT50 (fold-change v1-5) NT50 (fold-change v1-5) NT50 (fold-change v1-5) NT50 (fold-change v1-5) NT50 (fold-change v1-5) medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. figure 4 It is made available under a CC-BY-NC-ND 4.0 International license .

A B C D E B.1.1.7 B.1.1.7 E484K B.1.351 #1 B.1.351 #2 B.1.617.2 0.0072 0.0017 0.0063 0.0103 0.7168 103 103 103 103 103

102 102 102 102 102

101 101 101 101 101 NT50 relative to wt (%) NT50 relative to wt (%) NT50 relative to wt (%) NT50 relative to wt (%) NT50 relative to wt (%) 100 100 100 100 100 v1 v5 v1 v5 v1 v5 v1 v5 v1 v5

IgG) F stabilized native

wt B.1.351 #1B.1.351 #2B.1.1.7 B.1.1.7 E484KB.1.617 Abbott IgGRoche II Quant S (RBDRoche total)IgG) N (NCSiemens total) COV2TDiasorin (RBD S1/2Diasorin (S1/2total) Trimeric SiemensIgG) S sCOVGEuroimmun (S (RBDAbbott (S1 IgG) IgG)(NCcPass IgG) (sVNT50)

wt 1.00 0.85 0.83 0.94 0.83 0.90 0.83 0.76 0.53 0.74 0.83 0.79 0.83 0.81 0.53 0.88

B.1.351 #1 0.85 1.00 0.95 0.82 0.88 0.79 0.77 0.77 0.58 0.76 0.82 0.71 0.72 0.72 0.44 0.82

B.1.351 #2 0.83 0.95 1.00 0.80 0.87 0.80 0.73 0.73 0.54 0.72 0.76 0.66 0.69 0.67 0.43 0.79

B.1.1.7 0.94 0.82 0.80 1.00 0.83 0.88 0.76 0.72 0.47 0.72 0.77 0.72 0.74 0.72 0.43 0.82

B.1.1.7 E484K 0.83 0.88 0.87 0.83 1.00 0.78 0.78 0.80 0.60 0.79 0.85 0.75 0.73 0.74 0.45 0.86

B.1.617 0.90 0.79 0.80 0.88 0.78 1.00 0.69 0.63 0.42 0.61 0.68 0.60 0.69 0.65 0.43 0.76 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license . supplementary figure 1

4 4 4 4 4 A 10 r=0.9588 10 r=0.9590 10 r= 0.8203 10 r=0.8294 10 r=0.5100 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 103 103 103 103 103

2 2 2 2 2 v3-NT50 v4-NT50 v5-NT50 v5-NT50 v2-NT50 10 10 10 10 10

101 101 101 101 101 101 102 103 104 101 102 103 104 101 102 103 104 101 102 103 104 101 102 103 104 v1-NT50 v2-NT50 v3-NT50 v4-NT50 v1-NT50

0.0008 B C p=0.01801 ns ns ns ns D 10000 10000 1000

1000 1000 100 NT50 NT50 100 100 10 NT50 v5 (% of v1)

10 10 1 50 100 150 200 F M F M F M F M F M F M days post PCR v1 v2 v3 v4 v5

E r=0.1305 80 r=0.2759 80 r=0.3136 80 80 r=0.3702 80 r=0.3025 P=0.3288 P=0.005222 P=0.001404 P=0.002107 P=0.02411 60 60 60 60 60

40 40 40 40 40 age (years) age (years) age (years) age (years) 20 20 age (years) 20 20 20

0 0 0 0 0 101 102 103 104 101 102 103 104 101 102 103 104 101 102 103 104 101 102 103 104 v1-NT50 v2-NT50 v3-NT50 v4-NT50 v5-NT50

1 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license . supplementary figure 2

103 r=0.8505 r=0.9073 r=0.8802 r=0.8654 r=0.8511 r=0.7356 A P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 102

101 Diasorin S1/2 (S1/2 IgG) (AU/ml) v1 v2 v3 v4 v5 100 103 r=0.8749 r=0.9030 r=0.8692 r=0.8719 r=0.7500 r=0.7689 B P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 102 IgG) (AU/ml)

101 Diasorin trimeric S

stabilized native v1 v2 v3 v4 v5 0 (S 10 104 r=0.7051 P<0.0001 r=0.8543 r=0.8015 r=0.7660 r=0.7598 r=0.6149 C 103 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 102

101 Siemens COV2T Siemens COV2T (RBD total) (COI) 100

10-1 v1 v2 v3 v4 v5 3 10 r=0.8525 r=0.8893 r=0.8434 r=0.8131 r=0.7738 r=0.7516 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P=1.057e-011 D 102

101

0

(RBD IgG) (COI) 10 Siemens sCOVG v1 v2 v3 v4 v5 10-1 15 r=0.6743 r=0.6574 r=0.6793 r=0.6603 r=0.6729 r=0.4254 E P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P=0.0008716 10 Abbott 5 (NC IgG) (S/C)

v1 v2 v3 v4 v5 0 105 r=0.8440 r=0.8977 r=0.7980 r=0.8578 r=0.7492 r=0.7609 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 F 104

103

102 Abbott IgGII Quant (RBD IgG) (AU/ml) v1 v2 v3 v4 v5 101 103 r=0.5259 r=0.6106 r=0.6045 r=0.5940 r=0.6620 r=0.3971 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P=0.002025 G 102

101 Roche N 0 (NC total) (COI) 10 v1 v2 v3 v4 v5 10-1 104 r=0.7591 r=0.6855 r=0.7068 r=0.8254 r=0.8098 r=0.8087 P<0.0001 P=2.897e-009 P<0.0001 P<0.0001 P<0.0001 P<0.0001 H 103

102

1

Roche S 10

0 (RBD total) (U/ml) 10 v1 v2 v3 v4 v5 10-1 r=0.8456 10 r=0.8484 r=0.8960 r=0.8578 r=0.6316 r=0.6835 I P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001

1 Euroimmun (S1 IgG) (S/C)

v1 v2 v3 v4 v5 0.1 10 100 1000 10000 10 100 1000 10000 10 100 1000 10000 10 100 1000 10000 10 100 1000 10000 10 100 1000 10000 NT50 NT50 NT50 NT50 NT50 NT50 3 10 r=0.8168 r=0.8988 r=0.7814 J P=5.344e-029 P=1.021e-021 P=4.555e-013 102

101 cPass (sVNT50) v1 v5 100 10 100 1000 10000 101 102 103 104 101 102 103 104 NT50 NT50 NT50

2 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license . supplementary figure 3

ROC NT50≥50 1.0

0.8

0.6

Abbott IgGII Quant (RBD IgG) Roche S (RBD total) 0.4 Roche N (NC total) Siemens COV2T (RBD total) True positive fraction Diasorin S1/2 (S1/2 IgG) Diasorin Trimeric S (Sstablized native IgG) 0.2 Siemens sCOVG (S-RBD IgG) Euroimmun (S1 IgG) Abbott (NC IgG) cPASS 0.0 0.0 0.2 0.4 0.6 0.8 1.0 False positive fraction (1-specificity)

3 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

supplementary figure 4 A B C D E F wt B.1.1.7 B.1.1.7 E484K B.1.351 #1 B.1.351 #2 B.1.617.2 0.0143 0.2700 0.1207 0.0004 105 0.0013 105 105 0.2209 105 105 105

104 104 104 104 104 104

103 103 103 103 103 103 NT50 NT50 NT50 NT50 NT50 NT50

102 102 102 102 102 102

101 101 101 101 101 101 v1 v5 v1 v5 v1 v5 v1 v5 v1 v5 v1 v5 3% 4% 26% 19% 31% 9%

v1 v1 v1 v1 v1 v1 NT50 <30 30-100 3% 7% 28% 9% 28% 14% 100-500 500-1000 >1000 v5 v5 v5 v5 v5 v5

4 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license . supplementary figure 5

5 A 10 104

103

102 Abbott IgGII Quant (RBD IgG) (AU/ml) r=0.8345 r=0.7584 r=0.7819 r=0.7673 r=0.7256 r0.6854 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 101 4 B 10 103

102 Roche S 101 (RBD total) (U/ml) r=0.7628 r=0.7248 r=0.7977 r=0.7676 r=0.7261 r=0.6290 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 100 C 103

102

Roche N 101 (NC total) (COI) r=0.5318 r=0.4730 r=0.6043 r=0.5758 r=0.5415 r=0.4249 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 100 D 104 103

102

101 Siemens COV2T Siemens COV2T

(RBD total) (COI) 0 10 r=0.7409 r=0.7233 r=0.7904 r=0.7591 r=0.7155 r=0.6122 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 10-1 E 103 102

101 Diasorin S1/2

(S1/2 IgG) (AU/ml) r=0.8328 r=0.7706 r=0.8507 r=0.8243 r=0.7611 r=0.6760 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 100 3 F 10 102 IgG) (AU/ml)

101

Diasorin trimeric S r=0.7869 r=0.7192 r=0.7496 r=0.7066 r=0.6618 r=0.6047 stabilized native P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 0 (S 10 G 103 102

101

0

(RBD IgG) (COI) 10 Siemens sCOVG r=0.8290 r=0.7378 r=0.7317 r=0.7218 r=0.6947 r=0.6889 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 10-1 H 101

100 Euroimmun (S1 IgG) (S/C) r=0.8136 r=0.7216 r=0.7376 r=0.7198 r=0.6708 r=0.6547 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 10-1 1010

I 108

106

104

102 r=0.5332 r=0.4303 r=0.4546 r=0.4439 r=0.4278 r=0.4336 Abbott (NC IgG) (S/C) P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 100 J 103

102

101 cPass (sVNT50) r=0.8770 r=0.8202 r=0.8599 r=0.8200 r=0.7932 r=0.7599 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 100 K 105 104

103 wt NT50 102 r=0.9040 r=1 r=0.9353 r=0.8285 r=0.8470 r=0.8275 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 P<0.0001 101 101 102 103 104 105 101 102 103 104 105 101 102 103 104 105 101 102 103 104 105 101 102 103 104 105 101 102 103 104 105 wt NT50 B.1.1.7 NT50 B.1.1.7 E484K NT50 B.1.351 #1 NT50 B.1.351 #2 NT50 B.1.617.2

5 medRxiv preprint doi: https://doi.org/10.1101/2021.07.02.21259939; this version posted July 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license .

Table 1. Ability of serological assays to qualitatively identify presence of neutralizing antibodies.

Cut off TP FN TN FP n %SENS %SPEC %PPV %NPV ≥1.4 307 43 42 29 421 88 59 91 49 Abbott N (NC IgG) ≥4.5 166 184 66 5 421 47 93 97 26 Diasorin S1/2 (S1/2 IgG) ≥15 331 15 34 37 417 96 48 90 69 Siemens COV2T ≥1 346 4 9 62 421 99 13 85 69 (RBD total) ≥4.8 311 39 36 35 421 89 51 90 48 Roche N (NC total) ≥1 346 4 2 69 421 99 3 83 33 ≥0.8 350 0 2 69 421 100 3 84 100 Roche S (RBD total) ≥132 188 162 63 8 421 54 89 96 28 ≥1.1 315 33 51 20 419 91 72 94 61 Euroimmun (S1 IgG) ≥3.5 193 155 69 2 419 55 97 99 31 Diasorin Trimeric S (S stabilised native ≥ 13 340 9 18 53 420 97 25 87 67 IgG) Abbott IgGII IgG ≥50 347 0 3 68 418 100 4 84 100 (RBD IgG) ≥500 226 121 67 4 418 65 94 98 36 Siemens sCOVG ≥1 328 21 36 35 420 94 51 90 63 (RBD IgG) ≥4.8 217 132 71 0 420 62 100 100 35 .. ≥2 88 0 0 28 116 100 0 76 cPASS (sVNT50) ≥10 69 19 26 2 116 78 93 97 58 NT50 relative to wt ≥50 350 0 71 0 421 100 100 100 100

TP, true positives are samples that score positive in each assay and have NT50 >50 (the limit of detection in the neutralization assay). FN, false negatives score negative in serological assays but have NT50>50. TN, true negatives score negative in serological assays and have NT50<50. FP, false positives score positive in serological assays but have NT50<50. %SENS, sensitivity measured as the ration of TP over samples with NT50>50. %SPEC, specificity measured as the ratio of TN over samples that with NT50<50. %PPV, positive predictive value is the ratio of TP over total positive serological assay results. %NPV, negative predictive valus is the ratio of TN over total negative serological assay results. Top row is manufacturer recommended cut off, bottom FDA recommended cut off when available.