Next Article in Journal
Changes in Attitudes towards Influenza and Pneumococcal Vaccination during the Subsiding COVID-19 Pandemic—Results of a Longitudinal Survey Study among Risk Groups in Germany between 2021 and 2023
Previous Article in Journal
Pseudorabies Virus Glycoproteins E and B Application in Vaccine and Diagnosis Kit Development
Previous Article in Special Issue
Intranasal Immunization of Pneumococcal pep27 Mutant Attenuates Allergic and Inflammatory Diseases by Upregulating Skin and Mucosal Tregs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antibody Persistence and Risk of COVID-19 Infection: Insights from Modeling

by
Laurent Coudeville
1,†,
Eleine Konate
2,†,
Tabassome Simon
3,
Xavier de Lamballerie
4,
Scott Patterson
5,
Clotilde El Guerche-Séblain
1,* and
Odile Launay
2,6,* on behalf of AP-PH COVIBOOST Vaccine Trial Group
1
Global Medical, Vaccines, 69007 Sanofi Lyon, France
2
Assistance Publique Hôpitaux de Paris (APHP), Hôpital Cochin, CIC Cochin Pasteur, Inserm, 75014 Paris, France
3
Sorbonne Université, Department of Clinical Pharmacology and Clinical Research Platform (URC-CRC-CRB), 75012 Paris, France
4
Aix-Marseille University, IRD-Inserm-IRBA, 13005 Marseille, France
5
Biostatistics Department, Sanofi Vaccines, Swiftwater, PA 18370, USA
6
Université Paris Cité; Inserm, I-REIVAC, French Clinical Research Infrastructure Network (F-CRIN), 75014 Paris, France
*
Authors to whom correspondence should be addressed.
These authors have contributed equally to this work.
Vaccines 2024, 12(9), 1079; https://doi.org/10.3390/vaccines12091079
Submission received: 5 August 2024 / Revised: 4 September 2024 / Accepted: 12 September 2024 / Published: 21 September 2024
(This article belongs to the Special Issue Immune Correlates of Protection in Vaccines)

Abstract

:
Background: In this post hoc exploratory study of the APHP-COVIBOOST trial (NCT05124171), we used statistical modeling to describe the evolution of neutralizing antibody (nAb) titers over time, asses its impact on SARS-CoV-2 infection, and explore potential differences between three booster vaccine formulations (D614, B.1.351, and BNT162b2). Methods: Antibody titers were measured for 208 adult participants at day 28, 3 months, and 6 months using a microneutralization assay against different Omicron subvariants. We developed four specific Bayesian statistical models based on a core model, accounting for vaccine-specific antibody decline, boosting of nAb for breakthrough infection, and risk of infection according to nAb levels. The model findings were cross-verified using another validated microneutralization assay. Results: The decrease in nAb titers was significantly lower for the B.1.351 vaccine than for the other booster formulations. An inverse relationship was found between risk of infection upon exposure and nAb levels. With Omicron BA.1 data, these results translated into a positive relative vaccine efficacy against any infection over 6 months for the B.1.351 vaccine compared to the BNT162b2 (31%) and D614 (21%) vaccines. Conclusions: Risk of infection decreased with increasing nAb titers for all vaccines. Statistical models predicted significantly better antibody persistence for the B.1.351 booster formulation compared to the other evaluated vaccines.

1. Introduction

As of December 2023, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused over 6.9 million deaths with more than 772 million confirmed cases across the globe [1,2]. Although the World Health Organization (WHO) ended the public health emergency of international concern for coronavirus disease (COVID-19) in May 2023, governments are recommended to remain vigilant and maintain the established COVID-19 infrastructure [3]. The major challenges during the COVID-19 pandemic were the rapid emergence of novel SARS-CoV-2 variants with immune escape capacity and waning of vaccine-derived protection [4,5].
Furthermore, another related challenge has been identifying correlates of protection (CoPs), which are immune markers that can be reliably deployed to assess the efficacy of a vaccine [6]. An established CoP is often used as the primary endpoint in the approval of vaccines. However, establishing a CoP is a rigorous process involving evidence generation from multiple sources such as natural history studies, efficacy trials, and mechanistic studies, which in this case is further complicated by the rapid emergence of novel SARS-CoV-2 variants [6,7]. Therefore, although a number of recent studies have accumulated strong evidence that the anti-SARS-CoV-2 neutralizing antibody (nAb) titer could be a potential CoP for COVID-19 vaccines [6,8,9,10,11], wide acceptance of this titer as a CoP for COVID-19 vaccines is yet to be achieved [6].
We have recently published the results of a multicenter, randomized, single-blind trial comparing safety and immunogenicity 15 days after administration of a monovalent ancestral strain recombinant protein vaccine (Sanofi-GSK MVD614, Lyon, France), a monovalent beta variant recombinant protein vaccine (Sanofi-GSK MVB.1.351, VidPrevtyn Beta, Lyon, France), and a monovalent ancestral strain mRNA vaccine (BNT162b2, Pfizer-BioNTech, Mainz, Germany). The results showed that the B.1.351 vaccine produced higher nAb against a wide range of variants (ancestral strain, Beta, Delta, Omicron BA.1) [12]. Subsequently, 3 months (M3) and 6 months (M6) post-booster data also showed that the beta variant recombinant protein vaccine induced higher and durable cross-neutralizing antibodies against Omicron subvariants [13,14]. In the present post hoc exploratory analysis, we utilized these immunogenicity data along with observations of the number of infections to construct a Bayesian model with three objectives: (i) to investigate the evolution of nAb titers over time, (ii) to assess the relationship between nAb and infection, and (iii) to explore the potential differences between different vaccine booster formulations [6]. The outcomes of this analysis will provide insights into the interplay between antibody persistence and protection against infection over time and will contribute to the existing pool of knowledge on CoPs for COVID-19.

2. Methods

2.1. Study DATA

2.1.1. Study Population

The present study describes a post hoc exploratory analysis based on the results of APHP-COVIBOOST trial, a randomized, single-blinded multicenter clinical trial (ClinicalTrials.gov number, NCT05124171; EudraCT number, 2021-004550-33) [12]. Data from 208 adult participants previously primed with two doses of BNT162b2 and corresponding to the original per-protocol population at day 28 (D28) post vaccination were included in this analysis [12]. Participants were randomized into three groups based on the booster vaccine formulation received at baseline: Sanofi-GSK MVB.1.351 (B.1.351) (n = 65), Sanofi-GSK MVD614 (D614) (n = 73), and Pfizer-BioNTech BNT162b2 (BNT162b2) (n = 70). All data were collected between January and July 2022.

2.1.2. Immunological Assays

A virus microneutralization assay, previously tested for both Zika and SARS-CoV-2 viruses [15,16], was used for the detection of nAb at Aix-Marseille University Laboratories (Marseille, France). Although nAb titers were measured for different subvariants of Omicron in the COVIBOOST trial [12], only data for the Omicron BA.1 strain that was in circulation during the study period were used for modeling. In addition, the nAb titers against Omicron BA.4/5 strain data from the validated pseudovirus microneutralization assay (Monogram Biosciences, South San Francisco, CA, USA) were used as a part of the post hoc analysis, in order to cross-validate the model findings [17].

2.1.3. COVID-19 Detection

All participants were screened for COVID-19 on day 0 (D0) to meet the inclusion criteria. In this post hoc analysis, both symptomatic and asymptomatic COVID-19 infections were identified between D28 and M6 by testing the participants at D28, M3, and M6, or whenever they developed clinical symptoms. Positive cases were defined as having a positive test (RT-PCR or serology), anti-nucleocapsid protein serology, or a 4-fold or higher increase in nAb titers against the corresponding strains.

2.2. Statistical Analyses

2.2.1. Statistical Model

A specific Bayesian statistical model accounting for vaccine-specific antibody decline over time, boosting of nAb in case of breakthrough infection, and risk of infection according to nAb levels was developed for this study. The model is defined below:
X i ,   j + 1 g   ~   N f g X i ,   j g ,   t i ,   j + 1 g t i , j g ,   σ ,                                   i f   I n f i , j g = 0 N ( X 1 , j g + b ) ,   σ ,                                                             i f   I n f i , j g = 1 I n f i , j g   ~   B e r n o u l l i π i , j g                                                                                                                                   π i , j g = 1 exp [ t i , j g t i , j + 1 g λ j R i s k X t d t ]                                                         X i , j g = X ~ i , j g + ϵ i , j                                                                                                                                                               ϵ i , j   ~   U n i f o r m 1 ,   1                                                                                                                                                                                    
where:
t i , j g : Time at which data point j for subject i in group g is observed.
X i , j g : Actual log-transformed titer for subject i in group g at data point j.
X ~ i , j g : Observed log-transformed titer for subject i in group g at data point j.
I n f i , j g : Infection observed between time tgi,j and time tgi,j+1.
ϵ i , j : Titer measurement error.
b : Boosting effect in case of breakthrough infection.
σ : Standard deviation of the antibody persistence model.
π i , j g : Probability of infection between time tgi,j and time tgi,j+1.
λ j : Force of infection during the period between data point j and data point j + 1, i.e., rate at which susceptible individuals acquire infection per unit time [18].
R i s k X : Level of risk of infection for given titer.
f g X , t : Group-specific antibody decline as a function of the previously measured titer and time between measurements.
Although antibody persistence was specific to each study arm in this model, the relationship between antibody titers and protection against infection was common to all age groups. All estimations were performed based on the observed data on antibody titers (X) and occurrence of infection ( I n f i , j g ) between data points. The model used log-transformed antibody titers, i.e., an increase in X by 1 corresponds to an increase in titer dilution and an increase by 0 corresponds to subjects with no detectable titers.
Based on the general model described above, four specific models were considered with varying antibody persistence functions (fg(X,t)) and levels of risk (Risk(X)) for a given titer. Model 1 consisted of a linear antibody decline function along with a logistic curve for the evolution of protection with titers. For this model:
f g X , t = X d g t
R X = 1 1 1 + e ( X m ) / s
R i s k X = R X / R ( 0 )
where m and s are the location parameter and scale parameter of the logistic function, respectively, and dg is the antibody waning rate for the group g. For Model 1 and all other models, the level of risk obtained is normalized to 1 for participants with no detectable titers, i.e., Risk(0) = 1.
Model 2 is based on an exponential antibody decline function and the same logistic curve for the evolution of protection with titer as Model 1.
f g X , t = X e d g t
R X = 1 1 1 + e ( X m ) / s
R i s k X = R X / R ( 0 )
Model 3 is based on the same linear antibody decline function as Model 1 and on piecewise linear evolution of protection. Model 3 is defined as
f g X , t = X d g t
R X = 1 ,                                                                                         i f   X ( m s ) m X s ,                                                   i f   m s < X < m 0 ,                                                                                                               i f X   m
R i s k X = R X / R ( 0 )
where m is maximum titer for risk of infection (i.e., location parameter), and s is the rate of decrease (i.e., scale parameter) in risk according to titer (Risk = 1 if X = m − s).
Finally, Model 4 is based on the exponential antibody decline function and piecewise linear evolution of protection:
X , t = X e d g t
R X = 1 ,                                                                                         i f   X ( m s ) m X s ,                                                   i f   m s < X < m 0 ,                                                                                                               i f X   m
R i s k X = R X / R ( 0 )

2.2.2. Model Outcomes

In addition to parameter estimation, statistical model analyses generated a series of outcomes to aid the interpretation of model results. The posterior distributions of model parameters were used to estimate the following outcomes:
Estimated risk of infection according to titer given by the risk function:
R i s k ( X )
Estimated evolution of antibody titers in the absence of infection by study group in terms of geometric mean titer (GMT), the standard metric for population-level immunogenicity data [12,14]:
G M T j g = 1 n g i = 1 n g X i , j g
where ng is the number of subjects in study group g.
Relative vaccine efficacy between study groups:
R V E a , b = 1 n b j = 1 T i = 1 n a π i , j a n a j = 1 T i = 1 n b π i , j b
where na and nb are the number of subjects in study groups a and b (relative efficacy of group a versus group b), respectively, and T is the number of periods of observation.
Monthly evolution of protection against infection over time by study arm:
P r o t m g = 1 i = 1 n g R i s k ( f g X i , 1 g , m ) n g
where m is number of months since the first observed titer.

2.2.3. Statistical Methods

The model was coded and analyzed using RStan and loo packages in R programming [19,20]. The microneutralization assay data from the Aix-Marseille University laboratories, analyzed using the statistical models, were based on titers against Omicron BA.1 strain, one of the circulating strains during the period of the study. However, relevant data (i.e., titers at D28, M3, and M6) from the Monogram Biosciences laboratory were only available for Omicron BA.4/5 strain and the same were used for analyses.
Non-informative priors were used for parameters related to antibody decline and weakly informative priors for parameters related to risk of infection for a given titer. R-hat convergence diagnostic, effective sample sizes, and Pareto smoothed importance sampling diagnostic plots were used to assess the accuracy and robustness of model fit. Subsequently, expected log pointwise predictive densities (ELPDs), measured using the loo package (elpd_loo), were used to compare the fit estimated by the different models considered [20].

3. Results

3.1. Immunological Data and Infection Events

nAb titers against Omicron BA.1 strain were measured with the Aix-Marseille assay for all 208 participants at D28, for 205 participants at M3, and for 142 participants at M6 (Table 1) [12,15]. For the Monogram assay (Omicron BA.4/5 strain), the number of data points available was, however, lower: 161, 152, and 118 participants at D28, M3, and M6, respectively.
For the Aix-Marseille assay, a total of 85 infection data points were available for the analysis. Some differences according to vaccine group in the proportion of patients infected were observed: an infection was identified for 23/65 = 35.4% of the subjects included in the B.1.351 group, 32/73 = 43.8% in the D614 group, and 30/70 = 42.8% in the BNT162b2 group. For the Monogram assay, the more limited number of participants resulted in a lower number of available infection data points (71) (Table 2).

3.2. Model Selection

Satisfactory convergence was achieved for all four models (Table 3). The force of infection was found to decrease between two consecutive data points (λ1 vs. λ2) across the models. Furthermore, all four models predicted the lowest waning of antibody (dg) for the B.1.351 vaccine group and highest for the BNT162b2 vaccine group.
Model 2 and Model 4 provided a similar level of fit for the Aix-Marseille data and outperformed the other two models, as evident from the ELPD (−1245.4 vs. −1245.2) [20]. The ELPD for Model 4 differed from those for Model 1 and Model 3 by −16.5 and −16.0 units, respectively (Table S1). Thus, Model 2 and Model 4 were considered as the preferred models for the Aix-Marseille assay. Interestingly, leave-one-out (LOO) cross-validation predicted that Model 3 and Model 4 provided the best fit for the Monogram assay data. However, the differences in the levels of fit were limited for the remaining models (Table S2).

3.3. Observed vs. Model-Estimated Parameters

The model-estimated infection rates and GMT were compared with the observed infection rates and GMT to assess the quality of the models. For the Aix-Marseille data, the observed infection rates in the two time intervals (D28–M3 and M3–M6) lay within the 95% credible interval (CI) of the infection rates estimated by Model 2 and Model 4 (i.e., the preferred models) (Figure 1A). Similarly, the observed GMT also fell between the 95% CI of the estimated GMT, except for the GMT increase between M3 and M6 for the D614 vaccine (Figure 1B). The observed infection rates for the Omicron BA.4/5 strain (Monogram data) lay between the 95% CI of the infection rate estimated by Model 3 (preferred model). The estimated GMT values also accorded with the observed data except for the D614 vaccine group (Figure S1).

3.4. Model Outcomes for Aix-Marseille University Data (Omicron BA.1)

Risk of infection vs. microneutralization assay titer (against Omicron BA.1).
The risk curve was obtained from the risk function Risk(X) and is presented as the relative risk compared to an individual with no detectable titers. Both Model 2 and Model 4 established a correlation between antibody level and the risk of SARS-CoV-2 infection upon exposure. There was a significant reduction in the risk level with increasing microneutralization titers. Compared to subjects with no detectable antibodies, the relative risk of developing infection upon exposure was <10% for a titer ≥1280 (Figure 2). A similar relationship was also predicted by Model 1 and Model 3 (Figure S2).

3.4.1. Evolution of Antibody Titers (GMT) in the Absence of Infection

In the absence of infection, GMT against Omicron BA.1 steadily declined for all three study groups (i.e., three vaccine groups). The mean GMT estimated by Model 2 decreased from 188.3 at month 1 (M1) to 75.6 at month 6 (M6) for the B.1.351 vaccine group. The changes in the mean GMT of the D614 and BNT162b2 vaccine groups were 120.9 to 35.2 and 110.3 to 18.0, respectively. Similar changes in GMT values were predicted by Model 4. The rate of decrease in GMT was significantly lower for the B.1.351 vaccine than for the other vaccines (Figure 3).

3.4.2. Evolution of Protection against Infection over Time by Study Arm (against Omicron BA.1)

The results from Model 2 indicate that in the absence of reinfection, for the B.1.351 vaccine, protection against infection is maintained for at least 6 months at a higher level (39.9% [9.0–65.2]) as compared to the D614 (22.1% [95% CI: 2.0–44.4]) and BNT162b2 (8.8% [95% CI: 0.1–22.9]) vaccines. (Figure 4). A similar trend was also predicted by Model 4. However, the differences between the vaccine groups were not significant considering the large CIs.

3.4.3. Relative Vaccine Efficacy between Study Groups (against Omicron BA.1)

The relative efficacy of the B.1.351 vaccine is presented against those of the D614 and BNT162b2 vaccines. This efficacy was calculated using the expected evolution of titers and corresponding probability of infection over the 6 months for which observed data were obtained. As predicted by the antibody persistence and evolution of protection curves, the better protection conferred by the B.1.351 vaccine translated into a significantly positive relative efficacy of the B.1.351 vaccine as compared to those of the BNT162b2 (31% [95% CI: 15–45]) and D614 (21% [95% CI: 10–32]) vaccines. Model 4 also predicted a similar relative vaccine efficacy of the B.1.351 vaccine against those of the BNT162b2 (30% [95% CI: 17–41]) and D614 (21% [95% CI: 11–31]) vaccines.

3.5. Model Outcomes for Monogram Biosciences Data (Omicron BA.4/5)

Analysis of the Monogram assay data (titers against Omicron BA.4/5) with the preferred model (Model 3) showed similar outcomes to those obtained by the Aix-Marseille assay with Omicron BA.1 data (Figure S3). The antibody waning rate was lowest in the B.1.351 vaccine group vs. the D614 and BNT162b2 vaccine groups (0.28 [95% CI: 0.14–0.41], 0.34 [95% CI: 0.21–0.48], and 0.45 [95% CI: 0.31–0.60], respectively). Predicted GMT in the absence of infection remained significantly higher for the B.1.351 vaccine group. At M1, the mean GMTs for the B.1.351, D614, and BNT162b2 groups were predicted to be 824.7, 448.1, and 378.0, which changed to 317.8, 135.9, and 78.3, respectively, at M6. The relative effectiveness of the B.1.351 vaccine was found to be 27% against both the D614 (95% CI: 17–35) and BNT162b2 (95% CI: 17–36) vaccines. Finally, Model 3 predicted that for the B.1.351 vaccine, protection against infection is maintained for at least 6 months at a higher level (37.7% [95% CI: 11.0–62.4]), as compared to the D614 (22.7% [95% CI: 3.5–47.3]) and BNT162b2 (13.1% [95% CI: 0.2–37.8]) vaccines.

4. Discussion

Statistical modeling using real-world clinical data has been a powerful tool in predicting long-term antibody persistence against various pathogens as well as in establishing a relationship between immune markers and vaccine-derived protection [21,22,23]. In the present study, we analyzed clinical data from the APHP-COVIBOOST trial using Bayesian statistical modeling to predict the evolution of nAb titers over time and its relationship with the risk of infection in adults. Adapting the general model, we constructed and compared four specific models to ensure the best possible outcomes for two SARS-CoV-2 subvariants, Omicron BA.1 and BA.4/5. The four models used different mathematical formulations for the antibody decline and risk functions. Despite the differences in the construction of the models and the values of measured parameters, two common tendencies were evident from all four models: (i) the risk of infection is significantly correlated to the microneutralization assay titer, and (ii) compared to the other two study groups, antibody waning was lowest in the recipients of the B.1.351 vaccine.
Of the four models tested, Model 2, which presented an exponential antibody decline function and a logistic curve for the evolution of protection, and Model 4, which presented a linear antibody decline function and piecewise linear evolution of protection function, outperformed the other two models for Omicron BA.1 data and provided the best fit. Despite the differences in the mathematical constitution of the two models, there was remarkable similarity in the relative risk vs. antibody titer curves. The relative risk of infection (with respect to individuals with no nAb titers) decreased with increasing microneutralization assay titers and became very low for individuals with high nAb titers (<10% for titers ≥1280). This finding was consistent with the results of recent studies evaluating the relationship of nAb titers with the risk of symptomatic or asymptomatic SARS-CoV-2 infection, which also predicted a corresponding inverse relationship [6,8,9,10,11]. The projected GMT decreased over time in all three study groups but remained significantly higher in the B.1.351 vaccine group than in the other two vaccine groups. Furthermore, the model showed significantly positive relative vaccine efficacy for the B.1.351 vaccine compared to those for the D614 and BNT162b2 vaccines. Although, the clinical significance of this finding has not yet been established, the result is aligned with the findings of similar studies showing that vaccines conferring higher antibody titers present higher efficacy [9,10].
Monogram assay analyses with Omicron BA.4/5 data presented similar trends in outcomes compared to the Aix-Marseille assay analyses with Omicron BA.1 data. These trends included an inversely proportional relationship between the risk level and microneutralization assay titer, a significantly lower decrease in antibody titers over time, and a significantly positive relative vaccine efficacy for the B.1.351 vaccine as compared to the other booster vaccines.
Establishing a CoP for COVID-19 vaccines remains a challenge owing to the rapid emergence of novel SARS-CoV-2 variants [7]. In particular, studies assessing Omicron-specific CoPs are not yet widely reported [7]. Nevertheless, the results of this study showcase the value of modeling data on antibody persistence and vaccine-derived protection and are expected to guide public health considerations and evaluations of new vaccines in the future [24].
However, there are several limitations to this study. First, as the study describes a post hoc exploratory analysis, the case definitions were not predetermined. Second, the results are not transposable to other assays (in terms of technology and variants), and the characteristics of the study period, such as the level of virus circulation and immunity profile of the population, may have impacted the results. In particular, the rapid evolution of subvariants precludes any extrapolation of the current study results. Third, no results were derived for symptomatic cases considering their limited number. Additional limitations include the relatively small sample size, lack of exhaustive model search, and absence of a placebo control group for efficacy evaluation.

5. Conclusions

Based on a general Bayesian model, we constructed and used four specific statistical models to predict the antibody persistence and estimate the risk of COVID-19 infection by Omicron BA.1 and BA.4/5 strains. Model analyses revealed an inverse relationship between antibody level and the risk of any asymptomatic or symptomatic SARS-CoV-2 infection upon exposure and estimated a significantly better antibody persistence for the beta-adjuvanted B.1.351 booster vaccine formulation as compared to the other evaluated booster vaccine formulations. These results translated into a positive relative vaccine efficacy against any infection for the B.1.351 booster vaccine formulation as compared to the other booster vaccine formulations. These results highlight the importance of continued data collection for better documentation of the correlation between immune response and protection against COVID-19 for different assays, variants, and outcomes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/vaccines12091079/s1, Figure S1: Observed vs model-estimated parameters for Omicron BA.4/5 data: (A) infection rates and (B) GMT by study arms; Figure S2: Risk of infection vs microneutralization assay titer (against Omicron BA.1) estimated by (A) Model 1 and (B) Model 3; Figure S3: Outcomes of Model 3 (Omicron BA.4/5) based on the Monogram assay data: (A) Risk of infection vs microneutralization assay titer. (B) Evolution of antibody titers (GMT) in the absence of infection. (C) Evolution of protection against infection over time (monthly) by study arm; Table S1: Summary of leave-one-out cross validation analysis and Pareto k diagnostics for the four models (Omicron BA.1 data); Table S2: Summary of leave-one-out cross validation analysis and Pareto k diagnostics for the four models (Omicron BA.4/5 data).

Author Contributions

Conceptualization: L.C., C.E.G.-S. and O.L.; methodology: L.C., S.P., C.E.G.-S. and O.L.; software: L.C., S.P., C.E.G.-S. and O.L.; validation: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; formal analysis: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; investigation: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; resources: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; data curation: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; writing—original draft preparation: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; writing—review and editing: L.C., E.K., T.S., X.d.L., S.P., C.E.G.-S. and O.L.; visualization: L.C., S.P. and C.E.G.-S.; supervision, C.E.G.-S. and O.L.; project administration, C.E.G.-S. and O.L.; funding acquisition, C.E.G.-S. and O.L. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Sanofi and the French Ministries of Solidarity and Health and of Higher Education, Research, and Innovation.

Institutional Review Board Statement

The COVIBOOST trial was conducted in accordance with the Declaration of Helsinki and French law for biomedical research. It was approved by the “CPP Ile de France III” Ethics Committee and the French Health Products Safety Agency (ANSM).

Informed Consent Statement

Informed consent was obtained from all participants involved in the study.

Data Availability Statement

The datasets generated and/or analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Acknowledgments

Medical writing assistance was provided by Sridhar Banerjee (Sanofi), while editorial assistance and guidance on coordination of manuscript development were provided by Hanson Geevarghese (Sanofi).

Conflicts of Interest

L.C., S.P., and C.E.G.-S. are Sanofi employees and may hold stocks/shares in the company. O.L. reports a grant from the French Ministry of Health and grants or contracts from Pfizer, Sanofi-Pasteur, GlaxoSmithKline, MSD, MD, AstraZeneca, and Johnson & Johnson. T.S. reports research grant support (through the institution) from Novartis, Sanofi, Merck, AstraZeneca, Lilly, Daichy Sankyo, GlaxoSmithKline, and French Ministry of Health and speaker or consulting honoraria (e.g., DSMB, Advisory Board, and Educational Events) from Ablative Solutions, Air Liquide, AstraZeneca, Sanofi, Servier, Novartis, and 4Living Biotech.

Abbreviation

CI, credible interval; CoP, correlate of protection; COVID-19, coronavirus disease; ELPD, expected log pointwise predictive density; GMT, geometric mean titer; nAb, neutralizing antibody; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2; SD, standard deviation; WHO, World Health Organization.

References

  1. World Health Organization. COVID-19 Weekly Epidemiological Update: 24 November 2023; WHO: Geneva, Switzerland, 2023; pp. 1–25.
  2. WHO Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int/ (accessed on 11 December 2023).
  3. World Health Organization. COVID-19 Weekly Epidemiological Update: 20 July 2023; WHO: Geneva, Switzerland, 2023; pp. 1–15.
  4. Fernandes, Q.; Inchakalody, V.P.; Merhi, M.; Mestiri, S.; Taib, N.; Moustafa Abo El-Ella, D.; Bedhiafi, T.; Raza, A.; Al-Zaidan, L.; Mohsen, M.O.; et al. Emerging COVID-19 variants and their impact on SARS-CoV-2 diagnosis, therapeutics and vaccines. Ann. Med. 2022, 54, 524–540. [Google Scholar] [CrossRef] [PubMed]
  5. Arashiro, T.; Arima, Y.; Muraoka, H.; Sato, A.; Oba, K.; Uehara, Y.; Arioka, H.; Yanai, H.; Kuramochi, J.; Ihara, G.; et al. Coronavirus Disease 19 (COVID-19) Vaccine Effectiveness Against Symptomatic Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infection During Delta-Dominant and Omicron-Dominant Periods in Japan: A Multicenter Prospective Case-control Study (Factors Associated with SARS-CoV-2 Infection and the Effectiveness of COVID-19 Vaccines Study). Clin. Infect. Dis. 2023, 76, e108–e115. [Google Scholar] [CrossRef] [PubMed]
  6. Gilbert, P.B.; Donis, R.O.; Koup, R.A.; Fong, Y.; Plotkin, S.A.; Follmann, D. A COVID-19 Milestone Attained—A Correlate of Protection for Vaccines. N. Engl. J. Med. 2022, 387, 2203–2206. [Google Scholar] [CrossRef] [PubMed]
  7. Perry, J.; Osman, S.; Wright, J.; Richard-Greenblatt, M.; Buchan, S.A.; Sadarangani, M.; Bolotin, S. Does a humoral correlate of protection exist for SARS-CoV-2? A systematic review. PLoS ONE 2022, 17, e0266852. [Google Scholar] [CrossRef] [PubMed]
  8. Feng, S.; Phillips, D.J.; White, T.; Sayal, H.; Aley, P.K.; Bibi, S.; Dold, C.; Fuskova, M.; Gilbert, S.C.; Hirsch, I.; et al. Correlates of protection against symptomatic and asymptomatic SARS-CoV-2 infection. Nat. Med. 2021, 27, 2032–2040. [Google Scholar] [CrossRef] [PubMed]
  9. Fong, Y.; Huang, Y.; Benkeser, D.; Carpp, L.N.; Áñez, G.; Woo, W.; McGarry, A.; Dunkle, L.M.; Cho, I.; Houchens, C.R.; et al. Immune correlates analysis of the PREVENT-19 COVID-19 vaccine efficacy clinical trial. Nat. Commun. 2023, 14, 331. [Google Scholar] [CrossRef] [PubMed]
  10. Fong, Y.; McDermott, A.B.; Benkeser, D.; Roels, S.; Stieh, D.J.; Vandebosch, A.; Le Gars, M.; Van Roey, G.A.; Houchens, C.R.; Martins, K.; et al. Immune correlates analysis of the ENSEMBLE single Ad26.COV2.S dose vaccine efficacy clinical trial. Nat. Microbiol. 2022, 7, 1996–2010. [Google Scholar] [CrossRef] [PubMed]
  11. Gilbert, P.B.; Montefiori, D.C.; McDermott, A.B.; Fong, Y.; Benkeser, D.; Deng, W.; Zhou, H.; Houchens, C.R.; Martins, K.; Jayashankar, L.; et al. Immune correlates analysis of the mRNA-1273 COVID-19 vaccine efficacy clinical trial. Science 2022, 375, 43–50. [Google Scholar] [CrossRef] [PubMed]
  12. Launay, O.; Cachanado, M.; Luong Nguyen, L.B.; Ninove, L.; Lachâtre, M.; Ben Ghezala, I.; Bardou, M.; Schmidt-Mutter, C.; Lacombe, K.; Laine, F.; et al. Immunogenicity and Safety of Beta-Adjuvanted Recombinant Booster Vaccine. N. Engl. J. Med. 2022, 387, 374–376. [Google Scholar] [CrossRef] [PubMed]
  13. Launay, T.; Konate, E.; Cachanado, M.; Lachatre, M.; Ghezala, I.B.; Lacombe, K.; Laine, F.; Schmidt-Mutter, C.; de Lamballerie, X.; Simon, T. Persistence of neutralizing antibodies induced by recombinant protein vaccine SARS-CoV-2 Omicron variants strain B1.531. In Proceedings of the Journées Nationales D’infectiologie (JNI) 2023, Grenoble, France, 7–9 June 2023. [Google Scholar]
  14. Launay, O.; Gupta, R.; Machabert, T.; Konate, E.; Rousseau, A.; Claire, V.; Beckers, F.; Chicz, R.; Botelho-Nevers, E.; Cachanado, M.; et al. Beta-variant recombinant SARS-CoV-2 vaccine induces durable cross-reactive antibodies against Omicron variants. Res. Sq. (Prepr. Version 1) 2023. [Google Scholar] [CrossRef]
  15. Gallian, P.; Pastorino, B.; Morel, P.; Chiaroni, J.; Ninove, L.; de Lamballerie, X. Lower prevalence of antibodies neutralizing SARS-CoV-2 in group O French blood donors. Antiviral Res. 2020, 181, 104880. [Google Scholar] [CrossRef] [PubMed]
  16. Nurtop, E.; Villarroel, P.M.S.; Pastorino, B.; Ninove, L.; Drexler, J.F.; Roca, Y.; Gake, B.; Dubot-Peres, A.; Grard, G.; Peyrefitte, C.; et al. Combination of ELISA screening and seroneutralisation tests to expedite Zika virus seroprevalence studies. Virol. J. 2018, 15, 192. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, Y.; Borisov, O.; Kee, J.J.; Carpp, L.N.; Wrin, T.; Cai, S.; Sarzotti-Kelsoe, M.; McDanal, C.; Eaton, A.; Pajon, R.; et al. Calibration of two validated SARS-CoV-2 pseudovirus neutralization assays for COVID-19 vaccine evaluation. Sci. Rep. 2021, 11, 23921. [Google Scholar] [CrossRef] [PubMed]
  18. Kaslow, D.C. Force of infection: A determinant of vaccine efficacy? npj Vaccines 2021, 6, 51. [Google Scholar] [CrossRef] [PubMed]
  19. Stan Development Team. RStan: The R Interface to Stan. R Package; Version 2.21.2; Stan Development Team, 202. Available online: http://mc-stan.org/ (accessed on 11 December 2023).
  20. Vehtari, A.; Gelman, A.; Gabry, J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat. Comput. 2017, 27, 1413–1432. [Google Scholar] [CrossRef]
  21. Coudeville, L.; Bailleux, F.; Riche, B.; Megas, F.; Andre, P.; Ecochard, R. Relationship between haemagglutination-inhibiting antibody titres and clinical protection against influenza: Development and application of a bayesian random-effects model. BMC Med. Res. Methodol. 2010, 10, 18. [Google Scholar] [CrossRef] [PubMed]
  22. Desai, K.; Coudeville, L.; Bailleux, F. Modelling the long-term persistence of neutralizing antibody in adults after one dose of live attenuated Japanese encephalitis chimeric virus vaccine. Vaccine 2012, 30, 2510–2515. [Google Scholar] [CrossRef] [PubMed]
  23. López, E.L.; Contrini, M.M.; Mistchenko, A.; Kieffer, A.; Baggaley, R.F.; Di Tanna, G.L.; Desai, K.; Rasuli, A.; Armoni, J. Modeling the long-term persistence of hepatitis A antibody after a two-dose vaccination schedule in Argentinean children. Pediatr. Infect. Dis. J. 2015, 34, 417–425. [Google Scholar] [CrossRef] [PubMed]
  24. Michos, A.; Tatsi, E.B.; Filippatos, F.; Dellis, C.; Koukou, D.; Efthymiou, V.; Kastrinelli, E.; Mantzou, A.; Syriopoulou, V. Association of total and neutralizing SARS-CoV-2 spike -receptor binding domain antibodies with epidemiological and clinical characteristics after immunization with the 1(st) and 2(nd) doses of the BNT162b2 vaccine. Vaccine 2021, 39, 5963–5967. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Observed vs. model-estimated parameters for Omicron BA.1 data: (A) infection rates and (B) geometric mean titers (GMTs) by study arms. The error bars (in Panel A) and shaded areas (in Panel B) represent 95% credible interval. D28, M3, and M6 denote day 28, month 3, and month 6, respectively.
Figure 1. Observed vs. model-estimated parameters for Omicron BA.1 data: (A) infection rates and (B) geometric mean titers (GMTs) by study arms. The error bars (in Panel A) and shaded areas (in Panel B) represent 95% credible interval. D28, M3, and M6 denote day 28, month 3, and month 6, respectively.
Vaccines 12 01079 g001
Figure 2. Risk of infection vs. microneutralization assay titer (against Omicron BA.1) estimated by (A) Model 2 and (B) Model 4. The shaded areas represent 95% credible interval. Logarithmic scale has been used for the X-axes of the graphs to better visualize a long range of data.
Figure 2. Risk of infection vs. microneutralization assay titer (against Omicron BA.1) estimated by (A) Model 2 and (B) Model 4. The shaded areas represent 95% credible interval. Logarithmic scale has been used for the X-axes of the graphs to better visualize a long range of data.
Vaccines 12 01079 g002
Figure 3. Evolution of antibody titers (GMT) in the absence of infection as estimated by (A) Model 2 and (B) Model 4 (Omicron BA.1 data). The shaded areas represent 95% credible interval.
Figure 3. Evolution of antibody titers (GMT) in the absence of infection as estimated by (A) Model 2 and (B) Model 4 (Omicron BA.1 data). The shaded areas represent 95% credible interval.
Vaccines 12 01079 g003
Figure 4. Evolution of protection against infection over time (monthly) by study arm (against Omicron BA.1) estimated by (A) Model 2 and (B) Model 4. The shaded areas represent 95% credible interval.
Figure 4. Evolution of protection against infection over time (monthly) by study arm (against Omicron BA.1) estimated by (A) Model 2 and (B) Model 4. The shaded areas represent 95% credible interval.
Vaccines 12 01079 g004
Table 1. Number of participants with available antibody data.
Table 1. Number of participants with available antibody data.
DatasetTimepointB.1.351D614BNT162b2
Aix-MarseilleD28657370
M3637270
M6 *485044
MonogramD28564758
M3534554
M6 *424234
* Antibody data at M6 for subjects with infection between D28 and M3 were considered as missing (accounting for possible delay in boosting effect that impacts interpretation of titer evolution between M3 and M6).
Table 2. Number of infection events during the D28–M3 and M3–M6 time intervals.
Table 2. Number of infection events during the D28–M3 and M3–M6 time intervals.
DatasetTime IntervalB.1.351D614BNT162b2
Aix-MarseilleD28–M3101318
M3–M6131912
MonogramD28–M310118
M3–M6141711
Table 3. Summary results for the four models (Omicron BA.1 data).
Table 3. Summary results for the four models (Omicron BA.1 data).
ModelParameter *MeanSD95% CIneffectiveRhat
Model 1Force of infection, D28–M3 (λ1) 0.220.070.120.3845871
Force of infection, M3–M6 (λ2) 0.130.030.080.1920921
Location parameter of the risk function (m)2.881.330.35.1824391
Scale parameter of the risk function (s)1.540.550.482.689321.01
Antibody waning rate for B.1.351 group (d1)0.270.060.150.3921181
Antibody waning rate for D614 group (d2)0.390.060.270.579321
Antibody waning rate for BNT162b2 (d3)0.490.060.370.649871
Standard deviation for Ab persistence model ( σ )1.190.051.11.2949891
Boost in titers post infection (b)2.080.151.782.3769111
Model 2Force of infection, D28–M3 (λ1) 0.220.070.120.3855281
Force of infection, M3–M6 (λ2) 0.130.030.080.1956541
Location parameter (m)2.881.340.335.2739481
Scale parameter (s)1.50.540.562.744571
Antibody waning rate for B.1.351 group (d1)0.080.020.040.1174221
Antibody waning rate for D614 group (d2)0.140.020.090.1872891
Antibody waning rate for BNT162b2 (d3)0.280.040.210.3764101
Standard deviation for Ab persistence model ( σ )1.190.041.11.2860381
Boost in titers post infection (b)2.080.151.782.3889321
Model 3Force of infection, D28–M3 (λ1) 0.210.050.120.3175291
Force of infection, M3–M6 (λ2) 0.120.020.080.1785351
Location parameter of the risk function (m)6.680.925.488.9255861
Scale parameter of the risk function (s)72.642.8913.353341
Antibody waning rate for B.1.351 group (d1)0.270.060.150.3887151
Antibody waning rate for D614 group (d2)0.390.060.270.599671
Antibody waning rate for BNT162b2 (d3)0.490.060.370.688581
Standard deviation for Ab persistence model ( σ )1.190.051.11.2960961
Boost in titers post infection (b)2.080.151.772.3712,2951
Model 4Force of infection, D28–M3 (λ1) 0.210.050.120.3167041
Force of infection, M3–M6 (λ2) 0.120.020.080.1778411
Location parameter of the risk function (m)6.520.975.269.0456161
Scale parameter of the risk function (s)6.862.692.7713.1155241
Antibody waning rate for B.1.351 group (d1)0.080.020.040.1178181
Antibody waning rate for D614 group (d2)0.140.020.090.1878671
Antibody waning rate for BNT162b2 (d3)0.280.040.210.3775741
Standard deviation for Ab persistence model ( σ )1.190.051.11.2868211
Boost in titers post infection (b)2.080.151.782.38105351
Ab, antibody; CI, credible interval; SD, standard deviation. * The measured parameters were common in all three models, but the physical interpretations differed for m and s. Force of infection refers to the rate at which susceptible individuals acquire infection per unit time.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Coudeville, L.; Konate, E.; Simon, T.; de Lamballerie, X.; Patterson, S.; El Guerche-Séblain, C.; Launay, O., on behalf of AP-PH COVIBOOST Vaccine Trial Group. Antibody Persistence and Risk of COVID-19 Infection: Insights from Modeling. Vaccines 2024, 12, 1079. https://doi.org/10.3390/vaccines12091079

AMA Style

Coudeville L, Konate E, Simon T, de Lamballerie X, Patterson S, El Guerche-Séblain C, Launay O on behalf of AP-PH COVIBOOST Vaccine Trial Group. Antibody Persistence and Risk of COVID-19 Infection: Insights from Modeling. Vaccines. 2024; 12(9):1079. https://doi.org/10.3390/vaccines12091079

Chicago/Turabian Style

Coudeville, Laurent, Eleine Konate, Tabassome Simon, Xavier de Lamballerie, Scott Patterson, Clotilde El Guerche-Séblain, and Odile Launay on behalf of AP-PH COVIBOOST Vaccine Trial Group. 2024. "Antibody Persistence and Risk of COVID-19 Infection: Insights from Modeling" Vaccines 12, no. 9: 1079. https://doi.org/10.3390/vaccines12091079

APA Style

Coudeville, L., Konate, E., Simon, T., de Lamballerie, X., Patterson, S., El Guerche-Séblain, C., & Launay, O., on behalf of AP-PH COVIBOOST Vaccine Trial Group. (2024). Antibody Persistence and Risk of COVID-19 Infection: Insights from Modeling. Vaccines, 12(9), 1079. https://doi.org/10.3390/vaccines12091079

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop