Abstract
Background Antibodies to Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) have been shown to neutralize the virus in-vitro. Similarly, animal challenge models suggest that neutralizing antibodies isolated from SARS-CoV-2 infected individuals prevent against disease upon re-exposure to the virus. Understanding the nature and duration of the antibody response following SARS-CoV-2 infection is therefore critically important.
Methods Between April and October 2020 we undertook a prospective cohort study of 3555 healthcare workers in order to elucidate the duration and dynamics of antibody responses following infection with SARS-CoV-2. After a formal performance evaluation against 169 PCR confirmed cases and negative controls, the Meso-Scale Discovery assay was used to quantify in parallel, antibody titers to the SARS-CoV-2 nucleoprotein (N), spike (S) protein and the receptor-binding-domain (RBD) of the S-protein. All seropositive participants were followed up monthly for a maximum of 7 months; those participants that were symptomatic, with known dates of symptom-onset, seropositive by the MSD assay and who provided 2 or more monthly samples were included in the analysis. Survival analysis was used to determine the proportion of sero-reversion (switching from positive to negative) from the raw data. In order to predict long-term antibody dynamics, two hierarchical longitudinal Gamma models were implemented to provide predictions for the lower bound (continuous antibody decay to zero, “Gamma-decay”) and upper bound (decay-to-plateau due to long lived plasma cells, “Gamma-plateau”) long-term antibody titers.
Results A total of 1163 samples were provided from 349 of 3555 recruited participants who were symptomatic, seropositive by the MSD assay, and were followed up with 2 or more monthly samples. At 200 days post symptom onset, 99% of participants had detectable S-antibody whereas only 75% of participants had detectable N-antibody. Even under our most pessimistic assumption of persistent negative exponential decay, the S-antibody was predicted to remain detectable in 95% of participants until 465 days [95% CI 370-575] after symptom onset. Under the Gamma-plateau model, the entire posterior distribution of S-antibody titers at plateau remained above the threshold for detection indefinitely. Surrogate neutralization assays demonstrated a strong positive correlation between antibody titers to the S-protein and blocking of the ACE-2 receptor in-vitro [R2=0.72, p<0.001]. By contrast, the N-antibody waned rapidly with a half-life of 60 days [95% CI 52-68].
Discussion This study has demonstrated persistence of the spike antibody in 99% of participants at 200 days following SARS-CoV-2 symptoms and rapid decay of the nucleoprotein antibody. Diagnostic tests or studies that rely on the N-antibody as a measure of seroprevalence must be interpreted with caution. Our lowest bound prediction for duration of the spike antibody was 465 days and our upper bound predicted spike antibody to remain indefinitely in line with the long-term seropositivity reported for SARS-CoV infection. The long-term persistence of the S-antibody, together with the strong positive correlation between the S-antibody and viral surrogate neutralization in-vitro, has important implications for the duration of functional immunity following SARS-CoV-2 infection.
Introduction
Since appearing as a cluster of pneumonia cases in December 2019 in Wuhan, China, Coronavirus disease (COVID-19) has rapidly spread worldwide1. As of October 26th, there have been 43,187,134 cases, resulting in over 1.1 million deaths and a global health crisis, with significant social, economic and public health implications2. COVID-19 is caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), an enveloped RNA β-coronavirus3. Specific immunoglobulin (IgG) antibody responses to the SARS-CoV-2 trimeric spike (S) protein, nucleoprotein (N) protein and the receptor-binding domain (RBD) develop between 6-15 days following disease-onset4. The S-protein, which contains the RBD, binds to host cells via the angiotensin-converting-enzyme-2 (ACE-2) receptor, and membrane fusion occurs before viral entry5,6. The N-protein plays an important role in transcription enhancement and viral assembly7.
SARS-CoV-2-specific antibodies, particularly to the S- and RBD-antigens, have been shown to correlate with T-cell responses and viral neutralization in vitro as well as to protect against disease in animals, following passive transfer of convalescent serum or selected monoclonal antibodies8–12. It is unclear, however, whether re-infection can occur in humans who mount a humoral response following primary SARS-CoV-2 infection and achieve viral clearance. Recent case reports have emerged describing new respiratory samples positive for SARS-CoV-2 RNA after confirmed negativity, although these are few compared to the worldwide scale of infection, with several potential explanations proposed13,14. Furthermore, no recurrence of disease was reported in rhesus macaques or Syrian hamsters that were re-challenged in the presence of detectable endogenous antibodies (although protection against infection varied between studies).9,15,16. These findings highlight the importance of characterising humoral dynamics following SARS-CoV-2 infection.
SARS-CoV IgG and neutralizing antibodies have been shown to commonly persist up to 2-3 years post-infection, particularly in hospitalized patients,17,18 with recent reports demonstrating seropositivity as late as 12-17 years after infection19,20. Following severe disease caused by Middle Eastern Respiratory Syndrome (MERS), antibodies have been detected up to 34-months post-infection21,22. Existing longitudinal studies of SARS-CoV-2 are limited by inadequate modeling of antibody dynamics, short duration, low sampling density and frequency of longer-term follow-up23–32. Fitting Locally Estimated Scatterplot Smoothing (LOESS) or equivalent lines of best fit to the data23–25,33 also fails to provide a mathematical framework for evaluating long-term antibody responses.
In order to evaluate antibody kinetics and longevity following SARS-CoV-2 infection, we undertook the prospective Covid-19 Staff Testing of Antibody Responses Study (Co-STARS). Seropositive and symptomatic participants were followed up monthly with repeated antibody titer quantification. Detailed demographic, clinical and socioeconomic data were collected and mathematical models developed to characterize longitudinal humoral kinetics from initial antibody boosting to subsequent decay. To predict long-term antibody dynamics, we fitted two different models based on the gamma distribution: one which assumed persistent antibody decay34, and an alternate that allowed for an eventual plateau,35,36 to account for sustained antibody production by long-lived plasma cells (LLPCs).
Results
Participant Demographics
After providing informed consent, a total of 3555 individuals - all healthcare workers at Great Ormond Street Hospital - were enrolled in the study. Of this group, 349 were both symptomatic, seropositive by the MSD assay and provided 2 or more monthly samples for the primary outcome analysis of antibody dynamics. These 349 seropositive participants were followed up monthly for a maximum of 7 months and provided 1163 serial monthly serological samples. The median follow-up time per participant was 122 days (IQR 65 – 157 days) with a maximum follow up time of 262 days from symptom onset. The majority of participants 252/349 (72%) donated 3 or more samples with a maximum of 7 samples donated during follow up. Most seropositive participants were women (80%) with a mean age of 39 years representative of the underlying population structure of the hospital. The predominant symptoms reported were cough 225/349 (64%), myalgia 225/349 (64%), followed by ageusia and anosmia at 210/349 (60%) and 201/349 (58%) respectively.
Factors Associated with Increased Peak Antibodies and Rapid Decay
Multivariate analysis demonstrated that fever, rigors, ageusia, anosmia, a previous medical condition, high BMI and Black Asian Minority Ethnic (BAME) backgrounds were all associated with higher peak spike protein antibody titers (Table 1). No variables were identified to be independently associated with the rate of antibody decay.
Demographic details of study participants and variables associated with high peak antibody titers.
Observed Antibody Kinetics and Seroreversion
Serial monthly serological measurements from 349 participants who provided 2 or more samples following the onset of symptoms demonstrated a rapid rate of decay of the N-antibody relative to the S and RBD antibody (Figure 1). The spike antibody assay detected a total of 342/349 (98%) participants who were seropositive to any one of the S, RBD or N-antibodies. In comparison the RBD and N-assays detected 332/349 (95%) and 333/349 (95%) respectively. The sensitivity of the RBD and N-assays further declined with time relative to the S-antibody assay. At 200 days following the onset of symptoms, only 75% of 349 participants tested positive for N-antibody whereas 99% remained positive for S-antibody (Figure 2).
Samples from the same participant are linked with a thin black line and the red-dotted line is shown to indicate seroreversion. The gamma-plateau model is superimposed to show antibody trajectory: The predicted antibody trajectory (black line) is the median of the posterior distribution of the best model fit and 95% CI with and without individual effects (light blue and dark blue shading respectively) a) The spike (S) protein b) RBD antibody and c) the nucleoprotein (N) antibody with a relatively steep rate of decay.
Repeated serological measurements to the spike (S, Blue), RBD (Green) and the nucleoprotein (N, Red) demonstrating the time to a negative test from all possible starting positive tests.
Modeled Serological Reversion and Proportion of Positive Tests Over Time
Monte-Carlo Markov traces converged well for both the gamma distribution and decay-to-plateau curves demonstrating a stable model fit to the data (Supplementary Figure 1). The maximum R-hat for any parameter was 1.0035, while the minimum effective sample size (ESS) was 842.8 (Supplementary data, Table S1). Comparison of goodness of fit between models showed that for all antigens the decay-to-plateau model provided a better fit to the data than the gamma-decay model, although this difference was not statistically significant (Supplementary data, table S2). Even under the most pessimistic assumption of continuous gamma-decay we estimate that 95% of individuals following infection with SARS-CoV-2 will have measurable S-antibody until 465 days [95% CI 370-575] after the symptom start date. Under the gamma-plateau model S-antibody will remain detectable indefinitely, in line with recent reports on antibody kinetics following SARS-CoV-1 infection (Figure 3a)20. The most pessimistic gamma-decay model (lower bound) and most optimistic gamma-plateau model (upper bound) for each antibody are shown in Figure 3b. Under both models the N-antibody decayed to undetectable levels. Even after accounting for long-lived plasma cells under the gamma-plateau model, 75% of participants were predicted to have seroreverted N-antibody by 610 days [95% CI 420-530] whereas under the gamma-decay model 100% of participants had seroreverted N-antibody by 460 days [95% CI 420-530] following symptom onset.
Model based predictions of time to seronegativity. a, Comparison of the three tested antibodies against the S-protein (S, blue), receptor-binding domain (RBD, green) and nucleocapsid protein (N, red) for the gamma plateau model (top) and the gamma decay model (bottom). b, Differences between the two proposed models, the gamma plateau model (blue) and the gamma decay model (red) for the three tested antibodies (S; top right, RBD; bottom, and N; top left). Coloured lines represent the median estimates of the posterior density, while the shaded ribbons encompass the 95% confidence interval.
Antibody Peak, Half-Life and Plateau
Measured weekly average titer data for each of the S-, RBD-, and N-antibodies demonstrated that peak antibody response to infection was itself a plateau/slowly increasing line. Antibody titers rapidly increased during the first 3 weeks with prolonged high titers reached and maintained between week 4 and week 10 after the onset of symptoms. The peak antibody response for the S-antibody, RBD, and N-antibody from both raw weekly average serial titer and modeled data occurred at 40 [95% CI 30-63] days, 31 days [95% CI 26-38], and 35 [95% CI 31-42] days respectively. This was supported by the both the gamma-decay and gamma-plateau models which provided a similar close fit to this early stage of the humoral response (Figure 4 a, b, c).
a, b, c, d, e) Measured and Modeled Weekly Mean Antibody Titer. Real data (green), gamma-plateau model (red), gamma-decay model (blue) for a) the spike antibody, b) the RBD antibody and c) the N-antibody. Modeled half-lives of antibody decay for d) The gamma-plateau model and e) gamma-decay model. Colors represent the three different antibodies tested: those for the spike protein (blue), nucleocapsid protein (red), and receptor-binding domain (green).
The modeled half-life under the gamma-decay model and the gamma-plateau model were also very similar and both models showed a rapid decay of the N-relative to the RBD- and S-antibody. The half-life for the N-, RBD- and S-antibody was 60 days [95% CI 52-68], 102 days [95% CI 92-114] and 126 days [95% CI 112-146] respectively under the gamma decay model, while the half-lives under the gamma-plateau model were 60 days [95% CI 52-70], 110 days [95% CI 74-148], and 364 days [95% CI 212-997] respectively. The half-lives under the gamma-plateau model were widened as a consequence of being closer to the time to plateau where the half-lives are stretched by a flattening curve (Figure 4 d and e).
Under the gamma-plateau model, the S-antibody was characterized by a slow decay, with an eventual stabilized plateau at 1825 days [95% CI 250-3700] and none of the posterior probability distribution of the titers at the eventual plateau crossed the threshold for a negative test, whereas 75% of the posterior probability distribution for the N-antibody crossed the threshold for a negative test by 610 days.
Surrogate Neutralization Assay
There was a sigmoidal relationship between raw antibody titers and percentage binding/ACE-2 receptor blocking for both the S- and RBD-antibodies. The sigmoidal relationship demonstrated that above a threshold spike antibody titer of 8586 [95% CI 8160-9095] there was a dramatic increase in percentage binding/ACE-2 receptor blocking. When titer and blocking data were log transformed the relationship was linear with a strong positive correlation coefficient R2=0.72 and R2=0.77 for the spike antibody and the RBD respectively. In order to visualize the binding activity of spike and RBD antibodies at the plateau we mapped the maximum of the second derivative (the point at which the change in percentage binding increased the most) of the sigmoid curve to the final titers at the plateau (Figure 5). The point of maximum slope increase for the sigmoid curve was at 5430.11 [4553.50-6427.50] for the spike protein, and 5733.670 [4867.225-6423.525] for the receptor-binding domain. Whilst the full range of the distribution of spike antibodies were predicted to remain detectable indefinitely at plateau under the gamma-plateau model, only a small proportion of individuals were predicted to have titers sufficient to enable measurable functional binding under our surrogate neutralization assay.
a,b Percentage binding plotted against antibody titer. Red line represents the median amount of antibody titer at which the change in percentage binding is greatest, with the 95% CI indicated by dotted red lines. Black line is the median posterior distribution of the generalized logistic model, while the blue ribbon represents the 95% CI. c,d, Posterior distribution of the antibody titers at the long-term plateau as predicted by the gamma-plateau model. Shaded grey area corresponds to the threshold of seronegativity.
Discussion
This prospective cohort study of antibody responses following SARS-CoV-2 infection has demonstrated that 99% of 349 healthcare workers symptomatic with SARS-CoV-2 remained seropositive for the spike protein antibody 200 days after symptoms. Our study is the first to provide a mathematical modeling framework capable of predicting the long-term dynamics of the 3 key SARS-CoV-2 antibodies following natural infection. Even under our most pessimistic assumptions of continuous exponential decay, 95% of individuals were predicted to remain seropositive to S-antibody at 465 days [95% CI 370-575 days] while our more optimistic upper bound gamma-decay model predicted a permanent long-lasting plateau of detectable S-antibody.
These data contradict conclusions from studies that have reported rapid waning of antibodies after a few months29,31,32. Furthermore, our findings are in line with the duration of humoral responses observed following SARS-CoV and MERS infections17–20. Importantly, the long-lasting S- and RBD-antibodies also correlated well with a surrogate SARS-CoV-2 neutralization assay of ACE-2-receptor-blocking, strongly suggesting that long-term measurable S-antibody levels are functionally important. When the humoral correlates of protection against reinfection are known, our model of longitudinal S-antibody dynamics will therefore enable predictions to be made about the duration of long-lasting protective immunity following infection with SARS-CoV-2.
In contrast to the S-antibody, the N-antibody was observed to serorevert in 56/349 participants over the course of the study alone and had a modeled half-life of 60 days. This has important implications for diagnostic testing, epidemiological modeling and public health decision making that often rely on the N-antibody to estimate SARS-CoV-2 seroprevalence. This finding may also explain some unexpectedly low population level prevalence estimates in high burden countries37 and confound the finding that children with multi-system inflammation have a higher S:N ratio compared to adults30. Antibodies peaked at 30-40 days, this is significantly longer than other reports that are likely to have missed the prolonged peak/plateau due to inadequate sampling density38. It is notable that the delay between the SARS-CoV-2 epidemic curve and the mini-epidemic of paediatric multi-system inflammatory syndrome coincides temporally with peak antibodies39.
The persistence of detectable S- and/or RBD-antibody compared to the rapid decay of the N-antibody has also been observed in convalescent sera obtained from SARS survivors, seventeen years after infection19, although the exact underlying mechanisms warrant further investigation. Differences in the epitope structure40, immunogenicity and presentation to B-cells may distinctly impact the production, maturation and longevity of the plasma cells that secrete these antibodies41–44. Distinct T-helper cell interactions at the germinal centre may further determine B-cell and humoral dynamics, as previously observed in the context of the response to different HIV proteins45. Independently or in parallel, transcriptional programs and epigenetic imprinting may also selectively influence the kinetics and survival of both N-versus S-antibody-producing long-lived plasma cells43,44,46. Finally, cross-reactive memory, rather than naive B cells, may play a role in responses targeting nucleoprotein, given that it is more conserved across CoVs than RBD47.
To date, no studies have comprehensively modeled the nature and duration of antibody responses to different SARS-CoV-2 epitopes. Long et al26, Seow at al29, Ibarrondo et al27, demonstrated rapid decay of SARS-CoV-2 IgG and neutralizing antibodies within the first 3 months following infection, particularly in mildly symptomatic cases. In July, our earlier publicly available pre-print of serial antibody and modeled data demonstrated longer lasting spike antibody and rapidly decaying N-antibody48. In comparison, others have reported that the S-antibody and/or RBD-antibody correlate with neutralizing responses and decay slowly, persisting during the study period, up to at 90-150 days post-infection23–25,49. These studies, however, are limited by their shorter sampling time frame, lower sampling density and lack of appropriate modeling to predict antibody trajectory. Implementing LOESS lines of best fit to the data23–25,29 or comparing the variance of average antibody titers at different time intervals31,32,49 does not permit evaluation of long-term antibody trajectory.
Our study is strengthened by the density, frequency, and duration of longitudinal sampling collection. The parallel evaluation of absolute antibody titers by the chemiluminescent MSD assay to three major SARS-CoV-2 proteins also enabled us to demonstrate the decay of the N-antibody relative to the S- and RBD-antibodies. Importantly, this is the first study to provide a mathematical framework for long-term SARS-CoV-2 antibody responses, modeling both the peak and decay following infection and enabling realistic best-case and worst-case predictions of future antibody titers. Our work provides a detailed, shareable and reproducible model, with parameters that are useful for epidemiological purposes. Additionally, some of our parameters are fitted at both the population and individual level which is informative when inspecting risk factors and variability in the population. A third possible trajectory may be that the humoral response stabilizes but then continues to decline (‘plateau then decay’), albeit at a slower rate, as previously demonstrated in the context of vaccine-induced HPV responses and Hepatitis A infection35,36. Further serological measurements of seropositive recruits will take place in 6 months’ time to confirm which of these three models is superior in the longer term.
None of the seropositive healthcare workers identified in this study required hospitalization. This is important, given that the overwhelming majority of COVID-19 cases are not hospitalized. Our study population is, therefore, representative of most community SARS-CoV-2 infections50. Severe disease has however been associated with higher antibody titers and a longer duration of antibody response following both SARS and MERS17,21,51,52.
To date, no definitive quantitative or qualitative correlate-of-protection has been identified for SARS-CoV-2 infection, disease or onward transmission. Nevertheless, findings from animal studies support the role of neutralizing antibodies as a correlate of anti-viral immunity9,15,16. Mapping our surrogate neutralization assay to the final distribution of antibody titers under the gamma-plateau model suggested that some measurable functional binding (surrogate neutralization) of antibody would occur at plateau albeit in a small proportion. Whether this level of long-term detectable antibody is sufficient to induce sterilizing immunity and limit transmission, or primarily attenuates severity of disease, remains to be seen. Formal neutralization assays >1-year post infection are required to clarify this further. Finally, SARS-CoV-2-specific T- and B-memory cellular responses must also be characterised to accurately determine durability of immunity. Indeed, robust memory T-cell responses to specific SARS-CoV-2 peptides are elicited following infection and their magnitude correlated with antibody titers53–56; nevertheless, these responses were also observed in seronegative and/or pauci-symptomatic individuals53,56. Prospective evaluation of re-infection alongside humoral, T-cell, B-cell and mucosal IgA/IgG/cellular dynamics should therefore be an urgent priority, particularly during the evolving second wave of the pandemic.
One potential limitation of this study is the fact that only 38% of participants had an available confirmatory positive PCR result. In order to mitigate this concern, prior to the study taking place, a formal performance evaluation of the MSD assay was undertaken among 169 confirmed PCR positive SARS-CoV-2 participants that demonstrated a 97.9% sensitivity and 97.4% specificity at 21 days post infection57. This makes the proportion of false positive serological tests likely to be small and therefore have little impact on our findings.
Viral neutralization assays remain the gold-standard in vitro correlate of protection; as such, the lack of formal ‘authentic’ neutralization tests is another limitation of the study. However, ACE-2 receptor competition assays, such as the MSD competitive binding assay, have been shown to correlate well with formal viral neutralization assays, enabling their use as a suitable surrogate functional test58.
Whilst the severity of infection among our study participants is likely to be representative of community infection, our findings may be biased to healthcare workers. Recent studies have hypothesized that previous exposure to seasonal CoVs - to which healthcare workers may be disproportionately exposed - may confer some protection against SARS-CoV-212,19,53–56,58 and may need to be accounted for when modeling transmission or longevity dynamics59. Our estimates of the time-to-negativity are also dependent on the negative thresholds and lower limits of detection of the assay, respectively. However, our model fits, estimates of the rate of decay and the raw serial antibody titer trajectory (Figure 1) are not dependent on the threshold for a negative test.
In summary, this prospective cohort study has shown that the SARS-CoV-2 S-antibody, which correlated well with functional receptor blocking in-vitro, remained detectable in 99% of individuals up to 200 days post infection. In comparison, the N-antibody waned rapidly with a half-life of 60 days and 54/349 participants seroreverting over the course of the study. This study therefore has immediate consequences for diagnostic testing and public health decision making that often depend on the N-antibody as a reliable measure of past infection. Our most pessimistic continuous decay model, predicted that 95% of individuals would continue to have detectable spike antibody at 465 days while our gamma-plateau predicted that spike antibody would plateau at detectable levels indefinitely. The long-term presence of functional S (and RBD-antibody) has important implications for the duration of protective immunity following natural infection. It remains to be seen whether the SARS-CoV-2 vaccine candidates will replicate the long-lasting spike antibody duration observed and modeled here following natural infection.
Materials and methods
Study setting and design
Co-STARS is a 1 year single-centre, two-arm, prospective longitudinal cohort study of healthcare workers at a central London paediatric hospital Great Ormond Street Hospital for Children (GOSH). The study was approved to start by the United Kingdom NHS Health Research Authority on 29th April 2020 and registered on ClinicalTrials.gov (NCT04380896). Informed consent was obtained from all participants. The Study Protocol and Supplementary Materials submitted with this paper include detailed methods, power calculations and the data analysis approach.
Study participants
All hospital staff members ≥18 years of age were eligible for the study, provided they did not display symptoms consistent with SARS-CoV-2 infection at recruitment. Those significantly immunosuppressed or those who had previously received blood products (including immunoglobulins or convalescent sera) since September 2019 were excluded from the study.
Data Collection
After providing informed consent, participants undertook a detailed, standardised online questionnaire at study entry. This included socio-demographic factors, details of previous exposure to and symptomatic episodes consistent with COVID-19, any subsequent complications, previous SARS-CoV-2 diagnostic test results, past medical and contact history, and a comprehensive assessment of risk factors for exposure, susceptibility to infection and severe disease. Blood samples were also taken at baseline and each follow-up visit for determination of SARS-CoV-2 serology.
Measurement of SARS-CoV-2 serum antibody and viral RNA by PCR
Serum antibodies titers were measured by the Meso Scale Discovery (MSD) Chemiluminescent binding assay that simultaneously detects and quantifies anti-SARS-CoV-2 IgG specific for trimeric S-protein, RBD and N-protein. Assay qualification and performance were evaluated as described in our accompanying methods paper57. Briefly, IgG levels for the MSD assay were expressed as arbitrary units calibrated against a set of reference sera distributed by the National Institute of Biological Standards and Control (Potters Bar, UK) under the auspices of the World Health Organisation. SARS-CoV-2 real-time reverse-transcriptase polymerase chain reaction (RT-PCR) targeting the N-gene was performed following RNA extraction, as previously described by colleagues at our laboratory60.
Surrogate Neutralization Assay
An MSD® 96-well Custom Competition Assay designed to measure the inhibition of ACE-2 receptor binding to S or RBD by serum-derived antibody (MSD, Maryland) was run on 94 serial samples from 46 participants (two participants had 3 serial samples) in order to establish in vitro correlates of functional immunity.
Follow-Up Appointments
All seropositive participants were followed up monthly (ongoing) for repeat antibody testing. Seronegative participants will be followed up 6-monthly. At each follow-up visit, participants completed a shortened version of the baseline questionnaire, focussing on any recurrent COVID-19 exposure and/or symptoms.
Study outcomes
The primary outcome of the study was to establish humoral dynamics following SARS-CoV-2 infection. Nested sub-studies studies will explore the secondary outcome measures including the incidence of SARS-CoV-2 re-infection, the dynamics of the cellular response, IgA dynamics and the clinical and demographic factors that are associated with SARS-CoV-2 infection.
Statistical Analysis and Power Calculations
Power calculations were based on a negative exponential model of antibody decay from the peak using the pwr.f2.test function in R. We assumed a study power of 80% and explored a variety of hypothesized effect sizes (decreases in antibody titers over 1 year) and co-variates on the required study size with an alpha of 0.05 (Supplementary Materials, Study Protocol).
Statistical analysis and modeling of antibody dynamics
To calculate the proportion of individuals that serorevert over the course of the study, we performed a survival analysis to account for censoring using the survival package in R. An “event” was defined as a persistent negative test after the first positive test, while positive tests were counted as “censored” events.
The dynamics of antibody response following infection with SARS-CoV-2 were estimated by fitting two mixed effects models based on a gamma curve on every sample that had more than three antibody observations. Two gamma models were chosen (“gamma-plateau” and “gamma-decay”) to enable modeling of an optimistic upper bound estimate (eventual stabilization of decay to a plateau) to a pessimistic lower bound estimate (continuous exponential decay to zero). The gamma-decay model hypothesizes continuous antibody decay and does not account for long-lived humoral responses (although B-memory cells may still be induced to secrete antibody on challenge). On the other hand, the gamma-plateau model is based on the assumption that there are two phases of plasma cell production: ‘short-lived’ plasma cells are initially generated which secrete antibody with a short half-life, followed by a subsequent robust long-lived plasma cell (LLPC) response that maintains circulating, high-specificity, high-avidity antibody long-term61.
Gamma-decay and gamma-plateau models were fitted to the entire antibody response curve from day 0 (symptom onset) to the peak and then subsequent decay of each SARS-CoV-2 antibody. The confidence limits around the curves were derived by repeated sampling from the posterior distribution of the different model parameters, including individual effects.
The gamma-decay model assumed a uninterrupted continuous decay, given by the formula:
where f(t) is the log antibody titer at time t after symptom onset. The gamma function is described in terms of the shape (a) and scale (b), parameterized to reduce confounding of the parameters. β1 is the initial titer value at baseline; β2 determines the level of antibody rise; and u is an individual effect.
To account for the contribution of long lived plasma cells, a third term was added to the model, allowing for a long-term plateau expressed as:
Where β3 represents the long-term plateau, and k is the rate at which the latter term rises.
The relationship between ACE-2 receptor blocking and antibody titers was modeled with a 4 parameter generalized logistic curve, where the percentage binding at titer level t is given by:
with parameters a,b, and c that represent, respectively, the upper receptor blocking asymptote, the growth rate, and the titer at which maximum growth occurs. The parameter d is an asymmetry factor that affects the point of inflection on the y axis.
All the models were fitted using RSTAN in R (R: A language and environment for statistical computing. Foundation for Statistical Computing, Vienna, Austria). For each model, we ran 4 independent chains for 15,000 iterations for the gamma distributions, and 10,000 iterations for the sigmoid model. Model comparison was performed using Pareto-smoothed importance sampling leave-one-out cross-validation (PSIS-LOO) as implemented in the loo R package.
Ethics Statement
This study was approved by the UK Health Research Authority (www.hra.nhs.uk) and registered with www.clinical-trials.gov (NCT04380896). Written informed consent was obtained from all participants before recruitment to the study.
Data Availability
Data available upon request.
Funding statement
In line with UK government policy GOSH NHS Trust made the first diagnostic test available to all staff members in the hospital either as part of the study or outside the study. The GOSH charity provided the funding for follow up testing over the first 6 months of the study. LG declares funding from a Wellcome Trust Post-Doctoral Research Fellowship (201470/Z/16/Z www.wellcome.ac.uk). AS declares funding from Wellcome Trust Global Health Clinical Ph.D. Fellowship (220565/Z/20/Z www.wellcome.ac.uk). GOSH NHS Trust hosts a NIHR Funded Biomedical Research Centre (www.nihr.ac.uk) which provides infrastructure support permitting translational research. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests’ statement
The authors have declared that no competing interests exist.
Supplementary data
Acknowledgements
We would like to dedicate this article to the staff members who died of COVID-19 at Great Ormond Street Hospital during the first wave of the pandemic. We would also like to thank all the staff at Great Ormond Street Hospital who have taken part in the study. In addition, we are very grateful for all the hard work undertaken by the Great Ormond Street laboratory staff and the staff in the immunology laboratories both in the Camelia Botnar Laboratory and the Great Ormond Street Institute of Child Health who ensured that all the PCR tests and serological assays were completed in a timely manner. Finally, we would like to acknowledge the support of the Great Ormond Street Hospital Research & Development, Governance, Finance, Management, Estates, Operations and Communications departments.
Footnotes
↵† The Co-Stars Study Team: Dorcas Mirambe-Korsah Fernanda Torrente Jakub Wyszynski Victoria Gander Amy Leonard Louise Myers Aimee Vallot Camille Paillas Rose Fitzgerald Adam Twigg Rabia Manaf Lois Gibbons Hollie Powell Richard Nar-Dorh Ally Gray Elias Fernandez Aline Minja Emily Beech Waffa Girshab Pei Shi Chia Kate Webb Malti Nakrani Kim Gardiner Valerija Karaluka Karen Ryan Dorothy Lee Katie Groves Hamad Khan Shamime Nsubuga Olivia Rosie-Wilkinson Julia Spires Nuria Sanchez-Clemente Sapriya Kaur Natasha Carroll Jemma Efford Gabriel Bredin Celma Marisa Dos Santos Domingues Sophie Foxall Helen Ashton Abbey Afzal Sally Mainland Kate Crumpler Lucinda Dawson Claire Smith Maria Tabbu Laura Chiverton Jade Sugars Jordan Mooney Dorothy Chikusu Fariba Tahami Baratth Samy Shomona Begum Dhimple Patel Philippa Wiltshire Annie Susay Anna Ryan Luke Lancaster Kavita Thind Kate Speller Rachel Sterling Connor Tugulu Sandhya Ghurburrun Steffi Gray Joy Mugas Moe Kishma Kathleen Akpokomua Sophie White Eleana Pieri Sabina Shamsad Demi Alexandrou Odera Aguele Katherine Miles Anamika Jain Subishma Gautam Oliver Simms Rachel Goff Zarif Shams Tinya Chirinda Aaliya Nur Tarekur Rahman
The study was designed and implemented by: LG, AS, AT, TL, JH, RT, MH, TB, MJ, HW, ER, AM, CC, JB, KG, DG & the Co-Stars Team.
Testing for SARS-CoV-2 RNA and antibody was conducted by: LG, AS, AT, TL, JH, RT, MH, TB, MJ, HW, ER, AM, CC, JB, KG, DG
Funding: GOSH charity, Wellcome Trust (201470/Z/16/Z and 220565/Z/20/Z). GOSH NIHR Funded Biomedical Research Centre.
Trial registration number: NCT04380896.