Appropriate sampling and long follow-up are required to rigorously evaluate longevity of humoral memory after vaccination ========================================================================================================================= * Vitaly V. Ganusov ## Abstract One of the goals of vaccination is to induce long-term immunity against the infection and/or disease. However, evaluating the duration of protection following vaccination often requires long-term follow-ups that can conflict with the desire to rapidly publish results. Arunachalam et al. JCI 2023 followed individuals receiving third or fourth dose of mRNA COVID19 vaccines for up to 6 months and in finding that the levels of SARS-CoV2-specific antibodies (Abs) declined with similar rates for the two groups came to the conclusion that additional boosting is unnecessary to prolong immunity to SARS-CoV-2. However, this may be premature conclusion to make. Accordingly, we demonstrate that measuring Ab levels at 3 time points and only for a short (up to 6 month) duration does not allow to accurately and rigorously evaluate the long-term half-life of vaccine-induced Abs. By using the data from a cohort of blood donors followed for several years, we show that after re-vaccination with vaccinia virus (VV), VV-specific Abs decay bi-phasically and even the late decay rate exceeds the true slow loss rate of humoral memory observed years prior to the boosting. We argue that mathematical modeling should be used to better optimize sampling schedules to provide more reliable advice about the duration of humoral immunity after repeated vaccinations. “To consult the statistician after an experiment is finished is often merely to ask him to conduct a post mortem examination. He can perhaps say what the experiment died of. “ Ronald Fisher. ## Results and discussion The holy grail for vaccinologists is to induce long term effective immunity against infections, such as against the SARS-CoV-2, the cause of the recent pandemic1. It is well known that specific antibodies (**Abs**) to the spike protein of SARS-CoV-2 account for disease protection2, and it is of great interest to discover how long such antibodies persist after infection or various forms and regimens of vaccinations3,4. In one such study, Arunachalam *et al*. 4, compared the durability of neutralizing antibody (**nAb**) levels in persons that received either 3 or 4 vaccinations with Pfizer or Moderna mRNA vaccines. By using a simple linear regression, Arunachalam *et al*. 4 showed that the decay of Abs between 1 and 6 months was similar in individuals receiving 3rd or 4th vaccination allowing them to conclude that “The durability of serum antibody responses improves only marginally following booster immunizations with the Pfizer-BioNTech or Moderna mRNA vaccines.” inferring that extra boosters are unnecessary to confer long-term immunity. From multiple studies in mice and humans it is clear that the kinetics of Ab decay are generally more complex than a simple exponentially decaying function would suggest5–7. Specifically, after the peak immune response, Ab titers initially decline rapidly and approach their long-term decay kinetics over months. To investigate how the dynamics of Ab titers after revaccination relates to their long-term maintenance, we re-analyzed unique data from a cohort of long-term blood donors that recorded Ab titers to multiple vaccines over a long time period8. For this analysis, we selected 4 individuals illustrated in Figure 2 of Amanna *et al*. 8 that were followed for 20+ years and subsequently been revaccinated with vaccinia virus (**Figure 1** and **Supplemental Figure S1**). The kinetics of Ab decay prior to revaccination showed remarkable stability of VV-specific Abs with half-life times of 31-91 years (Figure 1A). By fitting a modification of the mathematical model of Ab response reported previously7 to the data on Ab expansion and contraction after revaccination, we found much smaller half-life times of VV-specific Abs 1.4-2.6 years as compared to the half-life times prior to revaccination (except of one volunteer (ID=514) showing no long-term decay after revaccination, Figure 1B). Given that these individuals were followed for nearly 2 years after VV revaccination their VV-specific Abs levels still decayed more rapidly than during years prior to boosting. This analysis indicates that long follow up studies are required before any conclusions can be made about the true longevity of vaccine-induced Abs. Additional analysis of the kinetics of Ab decay after the peak following revaccination with VV or tetanus vaccine suggested that these decays are not well described by a simple exponential curve indicating that different sub-populations of antibody-secreting cells (**ASCs**) with discordant lifespans may be present (**Supplemental Figure S1**). Interestingly, the long-term loss of tetanus-specific Abs occurred at a slower rate than the loss of VV-specific Abs, which may be surprising since VV is a live vaccine and tetanus vaccine is a protein-based vaccine. The difference here likely lies in the different times of follow up allowing more rigorous estimates of the half-life time of tetanus-specific Abs (**Supplemental Figure S1**). ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/06/29/2023.06.28.23291950/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2023/06/29/2023.06.28.23291950/F1) Figure 1: Stable levels of antibodies are reached only months after revaccination with vaccinia virus. **A**: We analyzed kinetics of Ab titers in four long-term blood donors from subjects 1-4 shown in Figure 2 of Amanna *et al.*8). These individuals were followed up for 20+ years during which they had been revaccinated with vaccinia virus. **B**: We fitted a mathematical model of humoral immune response (**eqn. (S.1)**)) to subsets of the data that include revaccination and estimated the rate Ab expansion (*ρ*), the proportion of Ab conversion into long-lived population (*f**m*), and the half-life of the humoral immunity (*T*1/2, panel B and see **Supplemental Table S1** for estimated model parameters). **C**: Kinetics of Ab response following VV revaccination in one volunteer (Donor 97) suggests infinite half-life of the long-term memory (see Main text for best fit parameters)7. **D**: Sparse measurements of Ab titers after vaccination allow for alternative mathematical models (**eqn. (S.1)** with different parameter sets) with drastically different predicted longevities of Abs. Here the markers are average Ab titers from Figure 1Ciii of Arunachalam *et al.*4, and lines are predictions of two alternative mathematical models with different assumed sub-populations of ASCs. In our recent work, we formulated mathematical models aimed to describe kinetics of ASCs and Ab titer following vaccination of humans7. The example of the immune response following VV immunization showed that to capture Ab loss, one needs to measure Ab titers frequently around the peak of the immune response (**Figure 1**C; VV-specific response for Donor 97 in the model given in **eqns. (S.2)–(S.3)** are *P* = 0.022, *T*on= 9.5 day, *ρ* = 1.02 day, *T*off = 14.79 day, *δ**A* = 0.15/*day, T*mem = 28 day, *A* = 6.7, *p* = 261.3/*day*) 7. We then used the model that described well the Ab response to VV (**eqn. (S.1)**) and selected parameters that would allow matching the average Ab titers measured at three time points chosen in the report of Arunachalam *et al*. 4. Importantly, we could choose two extreme sets of parameters that predict either relatively short-lived Ab response with half-life *T*1/2 = 0.3 years (model parameters: *A* = 20000, *ρ* = 0.091/*day, T*off = 30 day, *f**m* = 0, *δ*1 = 0.0058/*day, δ*2 = 0.13/*day*) or extremely long Ab response (*T*1/2 *>*100 years, model parameters: *A* = 20000, *ρ* = 0.48/day, *T*off = 7 day, *f**m* = 0.23, *δ*1 = 0.04/*day, δ*2 = 0.000013/day, **Figure 1**). These results further demonstrate that measuring Ab titers at only 3 time sparsely spaced points does not allow to rigorously evaluate the duration of humoral memory following revaccination. What are the possible solutions? Obviously, following revaccinated individuals for longer times may provide some useful information, but inapparent boostings of immunity following reinfection from currently circulating SARS-CoV-2 cannot be excluded. A better approach would be to improve the experimental design by using mathematical modeling-assisted power analyses that should indicate the most informative time points that samples should be collected at to provide the most reliable estimate of the parameter of immunity that is being advocated such as the Ab half-life times9,10. Following Ab response longitudinally in individual volunteers may also allow to use a more powerful approach of mixed effect modeling that may help to better define average and variability in Ab half-life times11. Ultimately, better collaborations between experimentalists and mathematical modelers may help to design more reliable experiments that provide more rigorous estimates of the longevity of humoral immunity afforded by vaccination, and (to paraphrase Ronald Fisher) “not let postmortem analysis to identify reasons experiment died of”. ## Supporting information Data on Ab dynamics following vaccinia virus or tetanus revacciation [[supplements/291950_file03.zip]](pending:yes) ## Data Availability as supplement to the paper ## Supplemental information ### Materials and Methods #### Experimental data Experimental data on kinetics of Ab titers in blood donor volunteers have been provided by Drs. Mark Slifka and Ian Amanna and are from their previous publication8. Here we specifically analyzed the data on the kinetics of VV and tetanus-specific Abs from four subjects shown in Figure 2 of Amanna *et al*. 8. We focused specifically on VV and tetanus-specific Abs because revaccination with these antigens induced robust recall responses in all individuals. To fit mathematical models to the Ab kinetics after revaccination we only included the first measurement before the peak of response and all the available data after revaccination. More thorough analysis of the recall responses in other subjects and to other vaccines will be provided elsewhere. The analyzed data (VV- and tetanus-specific Ab titers) are provided with the submission. Data on the dynamics ASCs and Ab titers following boost immunization with vaccinia virus (Donor 97) was published previously7. As far as we are aware the listed donor IDs are not known to anyone outside the research group that performed these measurements. #### Mathematical model for the kinetics of humoral immune response In experiments involving immunization, it is common to only measure the Ab titers (either antigen-binding titers or netraulizating titers). To describe the kinetics of Ab response for those data we adapt a simple *T*on *T*off model proposed to describe kinetics of the T cell response to viruses7,12,13. In this model, immune response gets activated at *T*on time and grows exponentially at a per capita rate *ρ*. The immune response peaks at time *T*off and then decays at two rates, *δ*1 and *δ*2, that biologically corresponds to ASCs with different longevity (given at frequency 1 *− f**m* and *f**m*, respectively): ![Formula][1] where *A* is the initial Ab titer (prior to boosting), *T*on is the time of the boosting (typically, *T*on = 0), *ρ* is the expansion rate of the Ab response, *T*off is the time of the peak of the immune response, *δ*1 and *δ*2 are the rates of Ab decay in the short- and long-lived plasma cells, and *f**m* is the proportion of the long-lived plasma cells in the total ASC populations. Note that the half-life of humoral memory is determined by the lowest of the two decay rates *δ*1 and *δ*2, *T*1/2 = log(2)/ min(*δ*1, *δ*2). To describe the dynamics of both ASC and Ab titer response following revaccination we use a model proposed previously7. In model, ASC number (*P*) expands exponentially after a delay *T*on, peaks at *T*off, and then ASCs die and convert to long-lived plasma cell phenotype at a rate *δ**A*. At *T*mem long-lived plasma cells are formed; they decay at rate *δ**M*. ASCs produce Abs at a rate *p* and Abs decay at a rate *δ**a*: ![Formula][2] ![Formula][3] where *p* is the per capita rate of production of antibodies by virus-specific plasma cells, and *δ**a* is the natural decay rate of antibodies assumed to be *δ**a* = 0.0495 day*−*1 corresponding to a half-life of 2 weeks14,15. In our example, *δ**M* = 0 as this provided the best fit of the model to data7. To describe the kinetics of Ab decay after the peak following VV or tetanus revaccination we used several alternative mathematical models that included a different number *n* of subpopulations of ASCs with different lifespans: ![Formula][4] where *A* is the Ab titer from all sub-populations, *A**i* is the size of the *i**th* sub-population and *δ**i* is the decay rate of ASCs in the *i**th* sub-population. The ultimate half-life of Abs is determined by the slowest decay rate *δ* = min(*δ**i*) with the corresponding half-life time *T*1/2 = ln(2)/*δ*. To determine how many sub-populations are needed to best describe the kinetics of Ab loss after the peak we vary *n* = 1 … 3 and compare the quality of the model fits using F-test for nested models16. #### Statistics The models were fitted to data using nonlinear least squares (NLS) by log-transforming the model predictions and the data. We used routine FindMinimum in Mathematica for NLS. ## ADDITIONAL FIGURES AND TABLES ![Supplemental Figure S1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/06/29/2023.06.28.23291950/F2.medium.gif) [Supplemental Figure S1:](http://medrxiv.org/content/early/2023/06/29/2023.06.28.23291950/F2) Supplemental Figure S1: Approach to stable Ab levels takes months after revaccination with vaccinia virus or tetanus vaccine. We fitted several alternative mathematical models for Ab decay (**eqn. (S.4)**) to data by varying the number of sub-populations for four volunteers revaccinated with VV (**A**) or tetanus (**B**) vaccines. In the analysis we included only the data from the peak Ab response at the first instance after an apparent revaccination (noted by the horizontal dashed line). In both cases the models with 2 subpopulations fitted the data better than the model with a single sub-population but a larger model (*n* = 3) did not improve the model fit quality but did provide a longer half-life times for persisting Abs. For each fit we provide the estimated half-life time (*T*1/2)based on the slowest decay rate. Parameters of the best fit models are show in **Supplemental Table S2**. Note the different scales of the follow up for boosting with two different vaccines. View this table: [Supplemental Table S1:](http://medrxiv.org/content/early/2023/06/29/2023.06.28.23291950/T1) Supplemental Table S1: Estimates of parameters of the Ab expansion and decay kinetics during revaccination with VV. We fitted mathematical models (**eqn. (S.1)**) to the data on Ab expansion and decay kinetics in four volunteers following revaccination with VV (Figure 1B). View this table: [Supplemental Table S2:](http://medrxiv.org/content/early/2023/06/29/2023.06.28.23291950/T2) Supplemental Table S2: Estimates of parameters of the Ab decay mathematical model fitted to data on Ab kinetics in four volunteers vaccinated with VV or tetanus vaccine. We fitted several mathematical models (**eqn. (S.4)**) to the data on Ab decay in four volunteers following revaccination with VV or tetanus vaccine (**Supplemental Figure S1**). ## Acknowledgement We would like to thank Mark Slifka and Ian Amanna for providing experimental data for this analysis. This was work supposed by the NIH grant (R01AI158963) to VVG. ## Abbreviations VV : vaccinia virus EU : ELISA units NLS : nonlinear least squares ASCs : antibody-secreting cells NLME : nonlinear mixed effects. * Received June 28, 2023. * Revision received June 28, 2023. * Accepted June 29, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Plotkin, S. A. 2008 Vaccines: correlates of vaccine-induced immunity. Clin Infect Dis, 47(3), 401–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/589862&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18558875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F29%2F2023.06.28.23291950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257400300015&link_type=ISI) 2. 2.Gilbert, P. B., Montefiori, D. C., McDermott, A. B., Fong, Y., Benkeser, D., Deng, W., Zhou, H., Houchens, C. R., Martins, K. et al. 2022 Immune correlates analysis of the mRNA-1273 COVID-19 vaccine efficacy clinical trial. Science (New York, N.Y.), 375, 43–50. doi:10.1126/science.abm3425. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abm3425&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34812653&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F29%2F2023.06.28.23291950.atom) 3. 3.Goel, R. R., Painter, M. M., Apostolidis, S. A., Mathew, D., Meng, W., Rosenfeld, A. M., Lund-green, K. A., Reynaldi, A., Khoury, D. S. et al. 2021 mRNA vaccines induce durable immune memory to SARS-CoV-2 and variants of concern. Science (New York, N.Y.), 374, abm0829. doi: 10.1126/science.abm0829. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abm0829&link_type=DOI) 4. 4.Arunachalam, P. S., Lai, L., Samaha, H., Feng, Y., Hu, M., Hui, H. S.-Y., Wali, B., Ellis, M., Davis-Gardner, M. E. et al. 2023 Durability of immune responses to mRNA booster vaccination against COVID-19. The Journal of clinical investigation, 133. doi:10.1172/JCI167955. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI167955&link_type=DOI) 5. 5.Murali-Krishna, K., Altman, J., Suresh, M., Sourdive, D., Zajac, A., Miller, J., Slansky, J. & Ahmed, R. 1998 Counting antigen-specific CD8+ T cells: A re-evaluation of bystander actiation during viral infection. Immunity, 8, 177–187. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1074-7613(00)80470-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9491999&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F29%2F2023.06.28.23291950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000072201800005&link_type=ISI) 6. 6.Miller, J., van der Most, R., Akondy, R., Glidewell, J., Albott, S., Masopust, D., Murali-Krishna, K., Mahar, P., Edupuganti, S. et al. 2008 Human effector and memory CD8+ T cell responses to smallpox and yellow fever vaccines. Immunity, 28(5), 710–22. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2008.02.020.&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18468462.&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F29%2F2023.06.28.23291950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256003900017&link_type=ISI) 7. 7.Le, D., Miller, J. D. & Ganusov, V. V. 2014 Mathematical modeling provides kinetic details of the human immune response to vaccination. Front Cell Infect Microbiol, 4, 177. doi:10.3389/fcimb.2014.00177. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fcimb.2014.00177&link_type=DOI) 8. 8.Amanna, I. J., Carlson, N. E. & Slifka, M. K. 2007 Duration of humoral immunity to common viral and vaccine antigens. N Engl J Med, 357(19), 1903–1915. doi:10.1056/NEJMoa066092. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa066092&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17989383&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F29%2F2023.06.28.23291950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250732500005&link_type=ISI) 9. 9.Rajakaruna, H. & Ganusov, V. V. 2022 Mathematical modeling to guide experimental design: T cell clustering as a case study. Bulletin of mathematical biology, 84, 103. doi:10.1007/s11538-022-01063-x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11538-022-01063-x&link_type=DOI) 10. 10.Majumder, B., Budhu, S. & Ganusov, V. V. 2023 Mathematical modeling suggests cytotoxic T lympho-cytes control growth of B16 tumor cells in collagin-fibrin gels by cytolytic and non-lytic mechanisms. bioRxiv. doi:[https://doi.org/10.1101/2023.03.28.534600](https://doi.org/10.1101/2023.03.28.534600). [Https://doi.org/10.1101/2023.03.28.534600](https://doi.org/10.1101/2023.03.28.534600). 11. 11.Pinheiro, J. C. & Bates, D. M. 2000 Mixed Effects Models in S and S-Plus. Springer. 12. 12.De Boer, R. J., Oprea, M., Antia, R., Murali-Krishna, K., Ahmed, R. & Perelson, A. S. 2001 Recruitment times, proliferation, and apoptosis rates during the CD8(+) T-cell response to lymphocytic choriomeningitis virus. J Virol, 75(22), 10 663–10 669. 13. 13.Ganusov, V. V., Milutinovic, D. & De Boer, R. J. 2007 IL-2 regulates expansion of CD4+ T cell populations by affecting cell death: insights from modeling CFSE data. J Immunol, 179(2), 950–957. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6OToiMTc5LzIvOTUwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMjkvMjAyMy4wNi4yOC4yMzI5MTk1MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. 14.Slifka, M., Antia, R., Whitmire, J. & Ahmed, R. 1998 Humoral immunity due to long-lived plasma cells. Immunity, 8(3), 363–72. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1074-7613(00)80541-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9529153&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F29%2F2023.06.28.23291950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000072705300009&link_type=ISI) 15. 15.Kindt, T. J., Osborne, B. & Goldsby, R. A. 2007 Kuby Immunology, 6th edition. New York: W. H. Freeman. 16. 16.Bates, D. M. & Watts, D. G. 1988 Nonlinear regression analysis and its applications. Hoboken, NJ: John Wiles & Sons, Inc. [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif