Research Article

Immunology and Inflammation

Modelling the response to vaccine in non-human primates to define SARS-CoV-2 mechanistic correlates of protection

University of Bordeaux, Department of Public Health, Inserm Bordeaux Population Health Research Centre, Inria SISTM, France
Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, France
Vaccine Research Institute, France
Inserm U955, France
Baylor Scott and White Research Institute, United States
Virus & Immunity Unit, Department of Virology, Institut Pasteur, France
CNRS UMR 3569, France
Department of Medical Microbiology, Amsterdam UMC, University of Amsterdam Amsterdam Infection & Immunity Institute, Netherlands
AP-HP, Hôpital Henri-Mondor Albert-Chenevier, Service d'Immunologie Clinique et Maladies Infectieuses, France
CHU Bordeaux, Department of Medical information, France

Jul 8, 2022

Open access
Copyright information

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Data availability
References
Article and author information
Metrics

Abstract

The definition of correlates of protection is critical for the development of next-generation SARS-CoV-2 vaccine platforms. Here, we propose a model-based approach for identifying mechanistic correlates of protection based on mathematical modelling of viral dynamics and data mining of immunological markers. The application to three different studies in non-human primates evaluating SARS-CoV-2 vaccines based on CD40-targeting, two-component spike nanoparticle and mRNA 1273 identifies and quantifies two main mechanisms that are a decrease of rate of cell infection and an increase in clearance of infected cells. Inhibition of RBD binding to ACE2 appears to be a robust mechanistic correlate of protection across the three vaccine platforms although not capturing the whole biological vaccine effect. The model shows that RBD/ACE2 binding inhibition represents a strong mechanism of protection which required significant reduction in blocking potency to effectively compromise the control of viral replication.

Editor's evaluation

This work should be of interest to a broad readership in infectious diseases, especially those people interested in modeling of infections. It combines statistical and mechanistic modeling to find assayable correlates of immunity for vaccines. This method could be relevant to many diseases or vaccines, although the particular markers identified here likely will be limited in their generalizability.

https://doi.org/10.7554/eLife.75427.sa0

Introduction

There is an unprecedented effort for SARS-CoV-2 vaccine development with 294 candidates currently evaluated (World Health Organization, 2021). However, variants of concern have emerged before the vaccine coverage was large enough to control the pandemics (Cobey et al., 2021). Despite a high rate of vaccine protection, these variants might compromise the efficacy of current vaccines (Kuzmina et al., 2021; Planas et al., 2021; Lustig et al., 2021; Zhou et al., 2021). Control of the epidemic by mass vaccination may also be compromised by unknown factors such as long-term protection and the need of booster injections in fragile, immuno-compromised, elderly populations, or even for any individual if protective antibody levels wane. Furthermore, the repeated use of some of the currently approved vaccine could be compromised by potential adverse events or by immunity against vaccine viral vectors (Greinacher et al., 2021). Finally, the necessity to produce the billions of doses required to vaccinate the world’s population also explains the need to develop additional vaccine candidates.

The identification of correlates of protection (CoPs) is essential to accelerate the development of new vaccines and vaccination strategies (Koch et al., 2021; Jin et al., 2021). Binding antibodies to SARS-CoV-2 and in vitro neutralization of virus infection are clearly associated with protection (Khoury et al., 2021; Yu et al., 2020; Earle et al., 2021; Feng et al., 2021). However, the respective contribution to virus control in vivo remains unclear (Zost et al., 2020), and many other immunological mechanisms may also be involved, including other antibody-mediated functions (antibody-dependent cellular cytotoxicity [ADCC], antibody-dependent complement deposition [ADCD], antibody-dependent cellular phagocytosis [ADCP]; Yu et al., 2020; Mercado et al., 2020; Tauzin et al., 2021), as well as T cell immunity (McMahan et al., 2021). Furthermore, CoP may vary between the vaccine platforms (Plotkin, 2013; Plotkin, 2020; Bradfute and Bavari, 2011; Dagotto et al., 2020).

Non-human primate (NHP) studies offer a unique opportunity to evaluate early markers of protective response (Muñoz-Fontela, 2020; Eyal and Lipsitch, 2021). Challenge studies in NHP allow the evaluation of vaccine impact on the viral dynamics in different tissue compartments (upper and lower respiratory tract) from day 1 of virus exposure (Yu et al., 2020; Mercado et al., 2020; Corbett et al., 2020). Such approaches in animal models may thus help to infer, for example, the relation between early viral events and disease or the capacity to control secondary transmissions.

Here, we propose to apply a model-based approach on NHP studies to evaluate (i) the immune mechanisms involved in the vaccine response and (ii) the markers capturing this/these effect(s) leading to identification of mechanisms of protection and definition of mechanistic CoP (Plotkin and Gilbert, 2012). First, we present a mechanistic approach based on ordinary differential equation (ODE) models reflecting the virus-host interaction inspired from models proposed for SARS-CoV-2 infection (Gonçalves et al., 2020; Kim et al., 2021; Gonçalves et al., 2021; Wang et al., 2020; Marc et al., 2021; Ke et al., 2021) and other viruses (Myers et al., 2021; Baccam et al., 2006; Goyal et al., 2019; Goyal et al., 2017). The proposed model includes several new aspects refining the modelling of viral dynamics in vivo, in addition to the integration of vaccine effect. A specific inoculum compartment allows distinguishing the virus coming from the challenge inoculum and the virus produced de novo, which is a key point in the context of efficacy provided by antigen-specific pre-existing immune effectors induced by the vaccine. Then, an original data mining approach is implemented to identify the immunological biomarkers associated with specific mechanisms of vaccine-induced protection.

We apply our approach to a recently published study (Marlin et al., 2021) testing a protein-based vaccine targeting the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein to CD40 (αCD40.RBD vaccine). Targeting vaccine antigens to dendritic cells via the surface receptor CD40 represents an appealing strategy to improve subunit-vaccine efficacy (Flamar et al., 2012; Zurawski et al., 2017; Cheng et al., 2018; Godot et al., 2020) and for boosting natural immunity in SARS-CoV-2 convalescent NHP.

We show that immunity induced by natural SARS-CoV-2 infection, as well as vaccine-elicited immune responses contribute to viral load control by (i) blocking new infection of target cells and (ii) by increasing the loss of infected cells. The modelling showed that antibodies inhibiting binding of RBD to ACE2 correlated with blockade of new infections and RBD-binding antibodies correlate with the loss of infected cells, reflecting importance of additional antibody functionalities. The role of RBD/ACE2-binding inhibition has been confirmed in two other vaccine platforms.

Results

A new mechanistic model fits the in vivo viral load dynamics in nasopharyngeal and tracheal compartments

The mechanistic model aims at capturing the viral dynamics following challenge with SARS-CoV-2 virus in NHP. For that purpose, we used data obtained from 18 cynomolgus macaques involved in the vaccine study reported by Marlin et al., 2021, and exposed to a high dose (1 × 10⁶ pfu) of SARS-CoV-2 administered via the combined intra-nasal and intra-tracheal route. The viral dynamics during the primary infections were characterized by a peak of genomic RNA (gRNA) production 3 days post-infection in both tracheal and nasopharyngeal compartments, followed by a decrease toward undetectable levels beyond day 15 (Figure 1—figure supplement 1). At the convalescent phase (median 24 weeks after the primary infection), 12 macaques were challenged with SARS-CoV-2 a second time, 4 weeks after being randomly selected to receive either a placebo (n=6) or a single injection of the αCD40.RBD vaccine (n=6) (Figure 1A). A third group of six naïve animals were infected at the same time. Compared to this naïve group, viral dynamics were blunted following the second challenge of convalescent animals with the lowest viral load observed in vaccinated animals (Figure 1B, Figure 1—figure supplement 2).

Figure 1 with 5 supplements see all

Download asset Open asset

Design of the study 1 and viral dynamics.

(A) *Study design*. Cynomolgus macaques (*Macaca fascicularis*), aged 37–58 months (8 females and 13 males). 24–26 weeks post-infection with SARS-CoV-2, 12 of these animals were randomly assigned in two experimental groups. The convalescent-vaccinated group (n=6) received 200 µg of αCD40.RBD vaccine. The other six convalescent animals were used as controls. Additional six age matched (43.7 months±6.76) cynomolgus macaques from same origin were included in the study as controls naïve from any exposure to SARS-CoV-2. Four weeks after immunization, all animals were exposed to a total dose of 10⁶ pfu of SARS-CoV-2 virus via the combination of intra-nasal and intra-tracheal routes. In this work, only data collected from the second exposure were considered. (B) Individual log₁₀ transformed genomic RNA (gRNA) viral load dynamics in nasopharyngeal swabs (top) and tracheal swabs (bottom) after the initial exposure to SARS-CoV-2 in naïve macaques (black, right) and after the second exposure in convalescent (blue, middle) and αCD40.RBD-vaccinated convalescent (green, left) groups. Horizontal red dashed lines indicate the limit of quantification.

Figure 1—source data 1 Genomic RNA (gRNA) viral load longitudinally measured in the trachea and nasopharynx after the second exposure in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data1-v2.xlsx
Download elife-75427-fig1-data1-v2.xlsx
Figure 1—source data 2 Genomic RNA (gRNA) viral load longitudinally measured in the trachea and nasopharynx after the first exposure for convalescent non-human primates (NHPs) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data2-v2.xlsx
Download elife-75427-fig1-data2-v2.xlsx
Figure 1—source data 3 Anti-spike IgG longitudinally measured post-immunization and quantified by Luminex in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data3-v2.xlsx
Download elife-75427-fig1-data3-v2.xlsx
Figure 1—source data 4 Quantification of the spike/ACE2-binding inhibition longitudinally measured post-immunization and quantified by Mesoscale Discovery (MSD) assay (in 1/ECL) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data4-v2.xlsx
Download elife-75427-fig1-data4-v2.xlsx
Figure 1—source data 5 Anti-N and anti-RBD binding antibodies longitudinally measured post-immunization and quantified by Mesoscale Discovery (MSD) assay (in AU mL^–1) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data5-v2.xlsx
Download elife-75427-fig1-data5-v2.xlsx
Figure 1—source data 6 Subgenomic RNA (sgRNA) viral load longitudinally measured in the trachea and nasopharynx after the second exposure in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data6-v2.xlsx
Download elife-75427-fig1-data6-v2.xlsx
Figure 1—source data 7 Antigen-specific T-cell response longitudinally measured post-exposure in % of CD4+ T cells measured by ICS in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data7-v2.xlsx
Download elife-75427-fig1-data7-v2.xlsx
Figure 1—source data 8 Antigen-specific T-cell response longitudinally measured post-exposure in % of CD8+ T cells measured by ICS in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data8-v2.xlsx
Download elife-75427-fig1-data8-v2.xlsx
Figure 1—source data 9 T-cell response expressing IFN-γ longitudinally measured post-exposure by ELISpot in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data9-v2.xlsx
Download elife-75427-fig1-data9-v2.xlsx
Figure 1—source data 10 Cytokine concentrations measured post-exposure in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data10-v2.xlsx
Download elife-75427-fig1-data10-v2.xlsx
Figure 1—source data 11 Quantification of the neutralization function of antibodies against three variants (B117, B1351, and D614G) longitudinally measured post-exposition (in ED50) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig1-data11-v2.xlsx
Download elife-75427-fig1-data11-v2.xlsx

We developed a mathematical model to better characterize the impact of the immune response on the viral gRNA and subgenomic RNA (sgRNA) dynamics, adapted from previously published work (Gonçalves et al., 2020; Kim et al., 2021; Baccam et al., 2006), which includes uninfected target cells (T) that can be infected (I₁) and produce virus after an eclipse phase (I₂). The virus generated can be infectious (V_i) or non-infectious (V_ni). Although a single compartment for de novo produced viruses (V) could be mathematically considered, two distinct ODE compartments were assumed for a better understanding of the model. We completed the model by a compartment for the inoculum to distinguish between the injected virus (V_s) and the virus produced de novo by the host (V_i and V_ni). In both compartments of the upper respiratory tract (URT), the trachea and nasopharynx, viral dynamics were distinctively described by this model (Figure 2A). Viral exchange between the two compartments was tested (either from the nasopharynx to the trachea or vice versa). However, as described in the literature (Gonçalves et al., 2021; Ke et al., 2020; Pinky et al., 2021) and demonstrated by the additional modelling work in Appendix 1 ‘Model building’, viral transport within the respiratory tract plays a negligible role in viral kinetics compared with viral clearance. Consequently, no exchange was considered in the model. Using the gRNA and sgRNA viral loads, we jointly estimated (i.e., shared random effects and covariates) the viral infectivity (β), the viral production rate (P), and the loss rate of infected cells (δ) in the two compartments. We assumed that gRNA and sgRNA were proportional to the free virus and the infected cells, respectively. This modelling choice relied on both biological and mathematical reasons (see section Materials and methods for more details). Due to identifiability issues, the duration of the eclipse phase (1/k), the clearance of free viruses from the inoculum (c_i) and produced de novo (c) were estimated separately by profile likelihood and assumed to be identical in the two compartments of the URT. In addition, infectious and non-infectious viruses were assumed to be cleared at the same rate. We estimated the viral infectivity at 0.95 × 10^–6 (CI_95% [0.18 × 10^–6; 4.94 × 10^–6]) (copies/mL)^–1 day^–1 in naïve animals, which is in the range of previously reported modelling results whether in the case of SARS-CoV-2 virus (Kim et al., 2021; Wang et al., 2020) or influenza (Myers et al., 2021; Baccam et al., 2006). We found estimates of the loss rates of infected cells of 1.04 (CI_95% [0.79; 1.37]) day^–1, corresponding to a mean half-life of 0.67 day. This estimation was consistent with previously published results obtained on SARS-CoV-2 virus showing the mean value of this parameter ranging from 0.60 to 2 day^–1 (i.e., half-life between 0.35 and 1.16 days) (Gonçalves et al., 2020; Kim et al., 2021; Gonçalves et al., 2021; Wang et al., 2020; Marc et al., 2021). The eclipse phase (3 day^–1) was found similar to the values commonly used in the literature (Gonçalves et al., 2020; Marc et al., 2021; Myers et al., 2021; Baccam et al., 2006). Here, we distinguished the clearance of the inoculum which was much higher (20 virions day^–1) as compared to the clearance of the virus produced de novo (3 virions day^–1). While the half-life of the virus de novo produced usually approximates 1.7 hr (i.e., c=10 day^–1) (Gonçalves et al., 2020; Gonçalves et al., 2021; Marc et al., 2021; Myers et al., 2021), because of this distinction, our model provided a higher estimation of 5.5 hr which remained in accordance with the estimations obtained by Baccam et al., 2006, on influenza A. Furthermore, the viral production by each infected cells was estimated to be higher in the nasopharyngeal compartment (12.1 × 10³ virions cell^–1 day^–1, CI_95% [3.15 × 10³; 46.5 × 10³]) as compared to the tracheal compartment (0.92 × 10³ virions cell^–1 day^–1, CI_95% [0.39 × 10³; 2.13 × 10³]). These estimations are in agreement with the observation of the intense production of viral particles by primary human bronchial epithelial cells in culture (Robinot et al., 2021). In particular, they are in the range of estimates obtained within the URT, either in NHP (Gonçalves et al., 2021) or in humans (Wang et al., 2020), with the product p × T₀ equals to 15.1 × 10⁸ (CI_95% [3.98 × 10⁸; 58.1 × 10⁸]) and 0.21 × 10⁸ (CI_95% [0.088 × 10⁸; 0.48 × 10⁸]) virions mL^–1 day^–1 in the nasopharynx and the trachea, respectively. By allowing parameters to differ between animals (through random effects), the variation of cell infectivity and of the loss rate of infected cells captured the observed variation of the dynamics of viral load. The variation of those parameters could be partly explained by the group to which the animals belong reducing the unexplained variability of the cell infectivity by 66% and of the loss rate of infected cells by 54% (Supplementary file 1). The model fitted well the observed dynamics of gRNA and sgRNA (Figure 2B).

Figure 2 with 2 supplements see all

Download asset Open asset

Mechanistic modelling.

(A) Description of the model in the two compartments: the nasopharynx and the trachea. (B) Model fit to the log₁₀ transformed observed genomic RNA (gRNA) viral loads in tracheal (top) and nasopharyngeal (bottom) compartments after the initial exposure to SARS-CoV-2 in naïve macaques (black, right) and after the second exposure in convalescent (blue, middle) and vaccinated (green, left) animals. Thick solid and dashed lines indicate mean viral load dynamics predicted and observed, respectively. Shaded areas indicate the 95% confidence intervals of the predictions. Dots represents observations. (C) Model predictions of unobserved quantities in the tracheal compartment for naïve (black, solid lines), convalescent (blue, dashed lines) and vaccinated (green, dotted lines) animals: target cells as percentage of the value at the challenge (top, left), infected cells (top, middle), productively infected cells (top, right), inoculum (bottom, right), infectious (bottom, left) and non-infectious virus (bottom, middle). Thick lines indicate mean values over time within each group. Shaded areas indicate the 95% confidence interval. Horizontal dashed red lines indicate the limit of quantification and horizontal solid red lines highlight the threshold of one infected cell.

Figure 2—source data 1 Volumes of the trachea and nasopharynx, and weights measured at the time of exposure in four non-human primates (NHPs) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig2-data1-v2.xlsx
Download elife-75427-fig2-data1-v2.xlsx
Figure 2—source data 2 Weights of the 18 non-human primates (NHPs) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig2-data2-v2.xlsx
Download elife-75427-fig2-data2-v2.xlsx
Figure 2—source data 3 Genomic RNA (gRNA) viral load measured in the trachea and nasopharynx in the two additional non-human primates (NHPs) receiving inoculum via intra-gastric and intra-nasal routes.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig2-data3-v2.xlsx
Download elife-75427-fig2-data3-v2.xlsx

Modelling of the dynamics of viral replication argues for the capacity of αCD40.RBD vaccine to block virus entry into host cells and to promote the destruction of infected cells

We distinguish the respective contribution of the vaccine effect and post-infection immunity on the reduction of the cell infection rate and the increase of the clearance of infected cells. Because blocking de novo infection and promoting the destruction of infected cells would lead to different viral dynamics profile (Figure 2—figure supplement 1), we were able to identify the contribution of each mechanism by estimating the influence of the vaccine compared to placebo or naïve animals on each model parameter. The αCD40.RBD vaccine reduced by 99.6% the infection of target cells in the trachea compared to the naïve group. The estimated clearance of infected cells was 1.04 day^–1 (95% CI 0.75; 1.45) in naïve macaques. It was increased by 80% (1.86 day^–1) in the convalescent macaques vaccinated by αCD40.RBD or not.

The mechanistic model allows predicting the dynamics of unobserved compartments. Hence, a very early decrease of the target cells (all cells expressing ACE2) as well as of the viral inoculum which fully disappeared from day 2 onward were predicted (Figure 2C). In the three groups, the number of infected cells as well as infectious viral particles increased up to day 2 and then decreased. We show that this viral dynamic was blunted in the vaccinated animals leading to a predicted maximum number of infectious viral particles in the nasopharynx and the trachea below the detection threshold (Figure 2C). The number of target cells would be decreased by the infection in the naïve and the convalescent groups, whereas it would be preserved in vaccinated animals.

The RBD-ACE2-binding inhibition is the main mechanistic CoP explaining the effect of the αCD40.RBD vaccine on new cell infection

In our study (Marlin et al., 2021), an extensive evaluation of the immunological response has been performed with quantification of spike-binding antibodies, antibodies inhibiting the attachment of RBD to ACE2, antibodies neutralizing infection, SARS-CoV-2-specific CD4⁺ and CD8⁺ T cells producing cytokines and serum cytokine levels (Figure 3, Figure 1—figure supplements 3–5). Therefore, based on our mechanistic model, we investigated if any of these markers could serve as a mechanistic CoP. Such a CoP should be able to capture the effect of the natural immunity following infection, associated or not to the vaccine (group effect) estimated on both the rate of cell infection and the rate of the loss of infected cells. To this aim, we performed a systematic screening by adjusting the model for each marker and we compared these new models with the model without covariates and with the model adjusted for the groups. In particular, our approach allowed us to benefit from all the information provided by the overall dynamics of the immunological markers after the exposure by integrating them as time-varying covariates (see the Materials and methods section for a detailed description of the algorithm). We demonstrate that the RBD-ACE2-binding inhibition measure is sufficient to capture most of the effect of the groups on the infection of target cells (Figure 4A and B). The integration of this marker in the model explains the variability of the cell infection rate with greater certainty than the group of intervention, reducing the unexplained variability by 87% compared to 66% (Supplementary file 1). The marker actually takes into account the variation between animals within the same group. Hence, it suggests that the levels of anti-RBD antibodies induced by the vaccine that block attachment to ACE2 are highly efficient at reflecting the neutralization of new infections in vivo. Furthermore, when taking into account the information provided by the RBD-ACE2-binding inhibition assay, the effect of the group of intervention was no longer significant (Supplementary file 1). Finally, we looked at the estimated viral infectivity according to the binding inhibition assay in each animal. A positive dependence was found between the viral infectivity and the RBD-ACE2-binding inhibition measure, linking an increase of 10³ AU of the marker, whether over time or between animals, with an increase of 1.8% (95CI% [1.2%; 2.3%]) of the viral infectivity (see Supplementary file 4). Accordingly, the values at the time of exposure were not overlapping at all, distinguishing clearly the vaccinated and unvaccinated animals (see Figure 4C).

Figure 3

Download asset Open asset

Harvest times and measurements.

Nasopharyngeal and tracheal fluids were collected at 0, 1, 2, 3, 4, 6, 9, 14, and 20 days post-exposure (d.p.e) while blood was taken at 0, 2, 4, 6, 9, 14, and 20 d.p.e. Genomic and subgenomic viral loads were measured by RT-qPCR. Anti-spike IgG sera were titrated by multiplex bead assay, anti-RBD, and anti-nucleocapside (N) IgG were titrated using a commercially available multiplexed immunoassay developed by Mesoscale Discovery (MSD, Rockville, MD). The MSD pseudo-neutralization assay was used to measure antibodies neutralizing the binding of the spike protein and receptor-binding domain (RBD) to the ACE2 receptor. Neutralizing antibodies against B.1.1.7, B.1.351, and D614G strains were measured by S-Fuse neutralization assay and expressed as ED50 (effective dose 50%). T-cell responses were characterized as the frequency of PBMC expressing cytokines (IL-2, IL-17a, IFN-γ, TNF-α, IL-13, CD137, and CD154) after stimulation with S or N sequence overlapping peptide pools. IFN-γ ELISpot assay of PBMCs were performed on PBMC stimulated with RBD or N sequence overlapping peptide pools and expressed as spot-forming cell (SFC) per 1.0 × 10⁶ PBMC.

Figure 4 with 2 supplements see all

Download asset Open asset

Immune markers.

(A) *Dynamics of biomarker selected as mechanistic correlate of protection (mCoP)*. Quantification of antibodies inhibiting RBD-ACE2 binding, measured by the Mesoscale Discovery (MSD) pseudo-neutralization assay (electro-chemiluminescence [ECL], in arbitrary unit [AU]) (top) and anti-RBD IgG titrated by ELISA assay (in IgG titer) (bottom). Thin lines represent individual values. Thick lines indicate medians of observations within naïve (black, solid line), convalescent (blue, dashed line), and αCD40.RBD-vaccinated convalescent (green, dotted line) animals. Shaded areas indicate 5th–95th confidence intervals of observations. (B) *Systematic screening of effect of the markers*. For every single marker, a model has been fitted to explore whether it explains the variation of the parameter of interest better or as well than the group indicator. Parameters of interest were β, the infection rate of ACE2+ target cells, and δ, the loss rate of infected cells. Models were compared according to the Bayesian information criterion (BIC), the lower being the better. The green line represents the reference model that includes the group effect (naïve/convalescent/vaccinated) without any adjustment for immunological marker (see Figure 3 for more details about measurement of immunological markers). (C) *Thresholds of inhibition of RBD-ACE2 binding*. Estimated infection rate (in (copies/mL)^–1 day^–1) of target cells according to the quantification of antibodies inhibiting RBD-ACE2 (in ECL) at exposure. Thin dotted lines and circles represent individual values of infection rates (right axis) and neutralizing antibodies (left axis). Shaded areas delimit the pseudo-neutralization/viral infectivity relationships within each group. (D) *Reproduction number over time*. Model predictions of the reproduction number over time in the trachea (right) and nasopharynx (left). The reproduction number is representing the number of infected cells from one infected cell if target cells are unlimited. Below one, the effective reproduction number indicates that the infection is going to be cured. Horizontal solid red lines highlight the threshold of one. Same legend than (A). (E) *Conditions for controlling the infection*. Basic reproduction number (R₀) at the time of the challenge according to the levels of antibodies inhibiting RBD-ACE2 binding (the lower the better) and of anti-RBD IgG-binding antibodies (the higher the better) assuming they are mechanistic correlates of blocking new cell infection and promoting infected cell death, respectively. The red area with R>1 describes a situation where the infection is spreading. The green area with R<1 describes a situation where the infection is controlled. The dotted red line delimitates the two areas. Black long dashed lines represent the values of neutralizing and binding antibodies measured at exposure. Observed values for three different animals belonging to the naïve (bottom, right), convalescent (bottom, left), and vaccinated (top, left) groups are represented. For each animal, individual values of R₀ were estimated considering their individual values of the model parameters (β and δ).

Figure 4—source data 1 Anti-N and anti-receptor-binding domain (RBD)-binding antibodies longitudinally measured post-immunization and quantified by ELISA in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig4-data1-v2.xlsx
Download elife-75427-fig4-data1-v2.xlsx
Figure 4—source data 2 Anti-receptor-binding domain (RBD) and anti-spike neutralizing antibodies longitudinally measured post-exposition and quantified by Mesoscale Discovery (MSD) assay (in electro-chemiluminescence [ECL]) in the study 1.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig4-data2-v2.xlsx
Download elife-75427-fig4-data2-v2.xlsx

In the next step, several markers (IgG-binding anti-RBD antibodies, CD8⁺ T cells producing IFN-γ) appeared to be associated to the rate of loss of infected cells (Figure 4—figure supplement 1A). Both specific antibodies and specific CD8⁺ T cells are mechanisms commonly considered important for killing infected cells. We retained the anti-RBD binding IgG Ab that were positively associated to the increase of the loss of infected cells. For unknown reason the IFN-γ response was high in unstimulated conditions in the naïve group. Thus, although this marker was associated with a decrease of the loss rate of infected cells, it appears essentially here as an indicator of the animal group. Further studies would be needed to fully confirm the place of IFN-γ response as a mechanistic marker.

A large part of the variation of the infection rate (71%) and loss rate of infected cells (60%) were captured by the two markers of CoP: the RBD-ACE2-binding inhibition and the anti-RBD-binding Ab concentration. Using the estimated parameters, the effective reproduction number could be calculated (R) which is representing the number of cells secondarily infected by virus from one infected cell (Figure 4D). When looking at this effective reproduction number according to the groups, the vaccinated animal presented from the first day of challenge an effective R below 1 meaning that no propagation of the infection started within the host. These results were consistent when taking the value of RBD-ACE2-binding inhibition at the time of the challenge without considering the evolution of the inhibition capacity over time (Figure 4—figure supplement 1B). This means that the dynamics of the viral replication is impacted very early during the infection process in immunized (i.e., both convalescent and vaccinated) animals and that vaccinated animals were protected from the beginning by the humoral response. Then, we looked at the threshold of the markers of interest leading to the control of the within-host infection (as defined by R<1) which was around 30,000 AU for the RBD-ACE2-binding inhibition assay. For the animals in the naïve and the convalescent groups, the observed values of binding inhibition measured by ECL RBD (the lower the better) and of IgG anti-RBD-binding antibodies (the higher the better) led to R>1, whereas in vaccinated animals, the value of ECL RBD led to R<1. Therefore, our modelling study shows that the inhibition of binding of RBD to ACE2 by antibodies is sufficient to control initial infection of the host (Figure 4E). According to the observed value of ECL RBD in vaccinated animals (e.g., 66 AU in Figure 4E), a decrease of more than 2 log₁₀ of the inhibition capacity (to reach 81,000 AU), due to variant of concern (VoC) or waning of immunity, would have been necessary to impair the control of the within-host infection. Moreover, a decrease of the neutralizing activity (i.e., increased ECL) could be compensated by an increase of cell death as measured by an increase of binding IgG anti-RBD as a surrogate. As an example, increasing IgG anti-RBD from 2.5 to 10 in the animal MF7 of the convalescent group would lead to a control of the infection.

In conclusion, the αCD40.RBD vaccine-elicited humoral response leads to the blockade of new cell infection that is well captured by measure of the inhibition of attachment of the virus to ACE2 through the RBD of the spike protein. Hence, the inhibition of binding of RBD to ACE2 is a promising mechanistic CoP. Indeed, this CoP fulfills the three criteria of leading to the best fit (lower BIC), the best explanation of interindividual variability, and fully captured the effect of the group of intervention.

The model revealed the same CoP related to another protein-based vaccine but not with mRNA-1273 vaccine

We took the opportunity of another study testing a two-component spike nanoparticle protein-based vaccine performed in the same laboratory and using the same immune and virological assays (Brouwer et al., 2021), measured only at the time of exposure, for applying the proposed model and methodology. In this study, six animals were vaccinated and compared to four naïve animals (Figure 5A and B). The good fit of the data (Figure 5C and D) allows for estimating the effect of the vaccine that appeared here also to decrease the infectivity rate (by 99%) and increase the clearance of the infected cells by 79%. Looking at the best mechanistic CoP following the previously described strategy, we ended here again with the inhibition of RBD binding to ACE2 as measured by ECL RBD. In fact, this marker measured at baseline before challenge fulfilled the three criteria: (i) it led to the best model in front of a model adjusted for group effect, (ii) it rendered the group effect non-significant, and (iii) it explained around 71% of the infectivity rate variability, compared to 65% of variability explained by the groups. Interestingly, here again, the inhibition assay led to a clear separation of the estimated rate of infectivity between vaccinees and the placebo group (Figure 5E).

Figure 5

Download asset Open asset

Study design and modeling results for the second study testing two-component spike nanoparticle vaccine.

(A) Study design. Cynomolgus macaques were randomly assigned in two experimental groups. Twelve, eight, and two weeks post-infection with SARS-CoV-2 virus, six of them were successively immunized with 50 µg of SARS-CoV-2 S-I53-50NP vaccine. The four other animals received no vaccination. Two weeks after the final immunization, all monkeys were exposed to a total dose of 10⁶ pfu of SARS-CoV-2 virus via intra-nasal and intra-tracheal routes. (B) Harvest times and measurements. Nasopharyngeal and tracheal fluids were collected at 0, 1, 2, 3, 4, 5, 6, 10, 14, and 21 days post-exposure (d.p.e.) while blood was taken at 0, 2, 4, 6, 10, 14, and 21 d.p.e. Genomic and subgenomic viral loads were measured by RT-qPCR. Anti-spike, anti-RBD, and anti-nucleocapside (N) IgG were titrated using a multiplexed immunoassay developed by Mesoscale Discovery (MSD, Rockville, MD) and expressed in AU mL^–1. The MSD pseudo-neutralization assay was used to quantify antibodies neutralizing the binding of the spike protein and RBD to the ACE2 receptor and results were expressed in electro-chemiluminescence (ECL). (C) Genomic viral load dynamics in nasopharyngeal and tracheal swabs after the exposure to SARS-Cov-2 in naïve (black, solid line) and vaccinated (green, dashed line) animals. Thin lines represent individual values. Thick lines indicate medians within each group. (D) Model fit to the log₁₀-transformed observed genomic RNA (gRNA) viral load in nasopharynx and trachea after the exposure to SARS-CoV-2 in naïve and vaccinated macaques. Solid thin lines indicate individual dynamics predicted by the model adjusted for groups. Thick dashed lines indicate mean viral load over time. (E) Thresholds of inhibition of RBD-ACE2 binding. Estimated infection rate of target cells ((copies/mL)^–1 day^–1) according to the quantification of antibodies inhibiting RBD-ACE2 binding (ECL) at exposure for naïve (black) and vaccinated (green) animals. Thin dotted lines and circles represent individual infection rates (right axis) and neutralizing antibodies (left axis). Thick dashed lines and dashed areas delimit the pseudo-neutralization/viral infectivity relationships within each group. (**C,D**) Horizontal red dashed lines represent the limit of quantification and shaded areas the 95% confidence intervals.The second study testing two-component spike nanoparticle vaccine.

Figure 5—source data 1 Anti-spike, anti-receptor-binding domain (RBD), and anti-N-binding antibodies quantified by Mesoscale Discovery (MSD) assay (AU mL^–1), and quantification of the spike/ACE2-binding inhibition by MSD assay (in 1/ECL), at the time of exposure in the study 2.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig5-data1-v2.xlsx
Download elife-75427-fig5-data1-v2.xlsx
Figure 5—source data 2 Genomic RNA (gRNA) and subgenomic RNA (sgRNA) viral loads longitudinally measured in the trachea and nasopharynx in the study 2.: https://cdn.elifesciences.org/articles/75427/elife-75427-fig5-data2-v2.xlsx
Download elife-75427-fig5-data2-v2.xlsx

Finally, we applied our approach to a published NHP study performed to evaluate several doses of mRNA-1273 vaccine (Corbett et al., 2020). Using available data, we compared the viral dynamics in the 100 µg, 10 µg, and placebo groups, enrolling a total of 12 rhesus macaques in a 1:1:1 ratio. Similar to the previous study, only immune markers measured at the time of exposure were available in this study, in addition to viral dynamics. We started from the same model as defined previously. We estimated a reduction of the infection rate by 97% but we did not find any additional effect. Looking at potential mechanistic CoP, we retained neutralization as measured on live cells with Luciferase marker. Although this marker led to the best fit and replaced the group effect (which was non-significant after adjustment for the marker), it explained only 15% of the variability of estimated viral infectivity, while 19% were explained by the groups.

In conclusion, we demonstrated, based upon challenge studies in NHP vaccinated with two different protein-based vaccine platforms, that both vaccines lead to the blockade of new cell infection. Neutralizing antibodies likely represent a consistent mechanistic correlate of protection (mCoP). This could change across vaccine platforms especially because mechanisms of action are different.

Discussion

We explored the mechanistic effects of three SARS-CoV-2 vaccines and assessed the quality of markers as mCoP. This model showed that neutralizing and binding antibodies elicited by a non-adjuvanted protein-based vaccine targeting the RBD of spike to the CD40 receptor of antigen presenting cells are reliable mCoP. Interestingly, we found the simpler and easier to standardize and implement binding inhibition assay may be more relevant to use as a CoP than cell-culture neutralization assays. This result has been replicated in another study testing a nanoparticle spike vaccine. The model was able to capture the effect of the vaccines on the reduction of the rate of infection of target cells and identified additional effects of vaccines beyond neutralizing antibodies. This latter consisted of increasing the loss rate of infected cells which was better reflected by the IgG-binding antibodies and CD8⁺ T-cell responses in the case of the CD40-targeting vaccine. One limitation of our study is that the prediction potential of our model relies on the range of the immune markers measured. However, our approach would allow a full exploitation of the data generated as in systems serology where non-neutralizing Ab functions, such as ADCC, ADCP, ADCD, and Ab-dependent respiratory burst (ADRB) are explored (Chung et al., 2015). The role of ADCC in natural infection has been previously shown (Dufloo et al., 2021), ADCD in DNA vaccine recipients (Yu et al., 2020) and with Ad26 vaccine (Alter et al., 2021). Here, we extended significantly these data by modelling the viral dynamic, showing that two other protein-based vaccines exert an additional effect on infected cell death which relied on the level of IgG anti-RBD-binding antibodies especially for the CD40.RBD-targeting vaccine. Measurements of other non-neutralizing Ab functions would probably also capture this additional effect.

The next question after determining which marker is a valid mCoP is to define the concentration that leads to protection, looking for a threshold effect that will help to define an objective (Khoury et al., 2021; Jin et al., 2021). In the context of SARS-CoV-2 virus, several emerged variants are leading to a significant reduction of viral neutralization as measured by various approaches. However, a 20-fold reduction of viral neutralization might not translate in 20-fold reduction of vaccine efficacy (Emary et al., 2021). First, there are many steps between viral neutralization and the reduction of viral infectivity or the improvement of clinical symptoms. Second, the consequences of a reduction of viral neutralization could be alleviated by other immunological mechanisms not compromised by the variant. In the context of natural immunity, when the level of neutralizing antibodies was below a protective threshold, the cellular immune response appeared to be critical (McMahan et al., 2021; Chandrashekar et al., 2020). We showed with our model that an improvement of infected cell destruction could help to control the within-host infection and is quantitatively feasible.

The control of viral replication is the key for reducing infectivity (Leung et al., 2020; Marks et al., 2021) as well as disease severity (Néant, 2021; Gutmann et al., 2021). According to our non-linear model linking the neutralization to the viral replication, a decrease of 4- to 20-fold in neutralization as described for the variants of concern (Planas et al., 2021; Zhou et al., 2021) is not enough, especially in the context of the response to CD40.RBD-targeting vaccine, to compromise the control of viral replication. The results showing a conserved effectiveness of mRNA vaccines in humans infected by the alpha or beta variants (Charmet et al., 2021), although a decrease of neutralization has been reported (Planas et al., 2021), are consistent with this hypothesis. However, this is highly dependent upon the mode of action of currently used vaccines and upon the VoC that may much more compromise the neutralization but being also intrinsically less pathogenic such as Omicron (Nyberg et al., 2022).

The analysis performed extended significantly the observation of associations between markers as previously reported for SARS-CoV-2 vaccine (Yu et al., 2020) and other vaccines (Kester et al., 2009) because it allows a more causal interpretation of the effect of immune markers. However, our modelling approach requires the in vivo identification of the biological parameters under specific experimentations. On the other hand, the estimation of parameters included in our model also provided information on some aspect of the virus pathophysiology. Notably, we found an increased capacity of virion production in nasopharynx compared to the trachea which could be explained by the difference in target cells according to the compartment (Travaglini et al., 2020). This result needs to be confirmed as it may also be the consequence of a different local immune response (Pizzorno et al., 2020). The choice of the structural model defining the host-pathogen interaction is a fundamental step in the presented approach. Here, it was well guided by the biological knowledge, the existing models for viral dynamics (Goyal et al., 2019; Gonçalves et al., 2021; Smith et al., 2018), and the statistical inference allowing the selection of the model that best fit the data. As the number of observations was relatively small in regard to the number of model parameters, we investigated overfitting issues. This was done using a bootstrap approach to evaluate the stability of confidence intervals of the estimated parameters. Results are provided in Appendix 2 ‘BICc as selection criteria and multiple testing adjustment’. Many modelling choices for the statistical model were made in this approach and more theoretical work evaluating the robustness of the results in their regards may be relevant for future works. In particular, we could relax the constraint of linear interpolation of marker dynamics by using simple regression models, allowing in the same time the integration of error model to account for measurement error for time-varying covariates (Dafni and Tsiatis, 1998; Carroll et al., 2006; Wu, 2009). Moreover, by construction, we assumed similar interindividual variability and effects of covariates within the two URT compartments as well as similar values for the viral infectivity and the loss rate of infected cells. Viral load dynamics measured in lungs being different from those in the URT (Lui et al., 2020; Goyal et al., 2020), the relaxation of this hypothesis of homogeneous physiological behavior in the URT may be pertinent to extend the model to the LRT. Finally, it should be underlined that the dynamics of the immune response has not been modelled as suggested for instance for B-cell response (Balelli et al., 2020). This clearly constitutes the next step after the selection of the markers of interest as done in the present work.

In conclusion, the modelling of the response to two new promising SARS-CoV-2 vaccines in NHP revealed a combination of effects with a blockade of new cell infections and the destruction of infected cells. For these two vaccines, the antibody inhibiting the attachment of RBD to ACE2 appeared to be a very good surrogate of the vaccine effect on the rate of infection of new cells and therefore could be used as a mechanistic CoP. This modelling framework contributes to the improvement of the understanding of the immunological concepts by adding a quantitative evaluation of the contributions of different mechanisms of control of viral infection. In terms of acceleration of vaccine development, our results may help to develop vaccines for ‘hard-to-target pathogens’, or to predict their efficacy in aging and particular populations (Pollard and Bijker, 2021). It should also help in choosing vaccine dose, for instance at early development (Rhodes et al., 2018) as well as deciding if and when boosting vaccination is needed in the face of waning protective antibody levels (Gaebler et al., 2021; Vanshylla et al., 2021), at least in NHP studies although the framework could be extended to human studies using mixed approaches of within and between hosts modelling (Goyal et al., 2022) providing that enough information is collected.

Materials and methods

Experimental model and subjects details

Request a detailed protocol

Cynomolgus macaques (Macaca fascicularis), aged 37–66 months (18 females and 13 males) and originating from Mauritian AAALAC certified breeding centers were used in this study. All animals were housed in IDMIT facilities (CEA, Fontenay-aux-roses), under BSL2 and BSL-3 containment when necessary (Animal facility authorization #D92-032-02, Préfecture des Hauts de Seine, France) and in compliance with European Directive 2010/63/EU, the French regulations and the Standards for Human Care and Use of Laboratory Animals, of the Office for Laboratory Animal Welfare (OLAW, assurance number #A5826-01, US). The protocols were approved by the institutional ethical committee ‘Comité d’Ethique en Expérimentation Animale du Commissariat à l’Energie Atomique et aux Energies Alternatives’ (CEtEA #44) under statement number A20-011. The study was authorized by the ‘Research, Innovation and Education Ministry’ under registration number APAFIS#24434-2020030216532863v1.

Evaluation of anti-spike, anti-RBD, and neutralizing IgG antibodies

Request a detailed protocol

Anti-spike IgG were titrated by multiplex bead assay. Briefly, Luminex beads were coupled to the spike protein as previously described (Fenwick et al., 2021) and added to a Bio-Plex plate (Bio-Rad). Beads were washed with PBS 0.05% tween using a magnetic plate washer (MAG2x program) and incubated for 1 hr with serial diluted individual serum. Beads were then washed and anti-NHP IgG-PE secondary antibody (Southern Biotech, clone SB108a) was added at a 1:500 dilution for 45 min at room temperature (RT). After washing, beads were resuspended in a reading buffer 5 min under agitation (800 rpm) on the plate shaker then read directly on a Luminex Bioplex 200 plate reader (Bio-Rad). Average MFI from the baseline samples were used as reference value for the negative control. Amount of anti-spike IgG was reported as the MFI signal divided by the mean signal for the negative controls.

Anti-RBD and anti-nucleocapside (N) IgG were titrated using a commercially available multiplexed immunoassay developed by Mesoscale Discovery (MSD, Rockville, MD) as previously described (Johnson et al., 2020). Briefly, antigens were spotted at 200–400 μg mL^–1 in a proprietary buffer, washed, dried, and packaged for further use (MSD Coronavirus Plate 2). Then, plates were blocked with MSD Blocker A following which reference standard, controls, and samples diluted 1:500 and 1:5000 in diluent buffer were added. After incubation, detection antibody was added (MSD SULFO-TAGTM Anti-Human IgG Antibody) and then MSD GOLDTM Read Buffer B was added and plates read using a MESO QuickPlex SQ 120 MM Reader. Results were expressed as arbitrary unit (AU) mL^–1.

Anti-RBD and anti-N IgG were titrated by ELISA. The nucleocapsid and the spike RBD (Genbank # NC_045512.2) were cloned and produced in Escherichia coli and CHO cells, respectively, as previously described (Flamar et al., 2012). Antigens were purified on C-tag column (Thermo Fisher) and quality-controlled by SDS-PAGE and for their level of endotoxin. Antigens were coated in a 96-well plates Nunc-immuno Maxisorp (Thermo Fisher) at 1 μg mL^–1 in carbonate buffer at 4°C overnight. Plates were washed in TBS Tween 0.05% (Thermo Fisher) and blocked with PBS 3% BSA for 2 hr at RT. Samples were then added, in duplicate, in serial dilution for 1 hr at RT. Non-infected NHP sera were used as negative controls. After washing, anti-NHP IgG coupled with HRP (Thermo Fisher) was added at 1:20,000 for 45 min at RT. After washing, TMB substrate (Thermo Fisher) was added for 15 min at RT and the reaction was stopped with 1 M sulfuric acid. Absorbance of each well was measured at 450 nm (reference 570 nm) using a Tristar2 reader (Berthold Technologies). The EC50 value of each sample was determined using GraphPad Prism 8 and antibody titer was calculated as log (1/EC50).

The MSD pseudo-neutralization assay was used to measure antibodies neutralizing the binding of the spike protein to the ACE2 receptor. Plates were blocked and washed as above, assay calibrator (COVID- 19 neutralizing antibody; monoclonal antibody against S protein; 200 μg mL^–1), control sera, and test sera samples diluted 1:10 and 1:100 in assay diluent were added to the plates. Following incubation of the plates, an 0.25 μg mL^–1 solution of MSD SULFO-TAGTM-conjugated ACE2 was added after which plates were read as above. Electro-chemiluminescence (ECL) signal was recorded.

Viral dynamics modelling

Request a detailed protocol

The mechanistic approach we developed to characterize the impact of the immune response on the viral gRNA and sgRNA dynamics relies on a mechanistic model divided in three layers: first, we used a mathematical model based on ODEs to describe the dynamics in the two compartments, the nasopharynx and the trachea. Then, we used a statistical model to take into account both the interindividual variability and the effects of covariates on parameters. Finally, we considered an observation model to describe the observed log₁₀ viral loads in the two compartments.

For the mathematical model, we started from previously published models (Gonçalves et al., 2020; Kim et al., 2021; Baccam et al., 2006) where the nasopharynx and trachea were respectively described by a target cell limited model, with an eclipse phase, as model of acute viral infection assuming target-cell limitation (Baccam et al., 2006). We completed the model by adding a compartment for the inoculum that distinguishes the injected virus (V_s) from the virus produced de novo (V_i and V_ni). To our knowledge, this distinction has not been proposed in any previous work. Two main reasons led us to make this choice. First, it allowed us to study the dynamics of the inoculum, in particular during the early phase of viral RNA load dynamics. Second, as described in more detail below, it gave us the opportunity to use all the information provided by the preclinical studies, such as the known number of inoculated virions, to define the initial conditions of the ODE model rather than estimating or randomly fixing them for V_i and V_ni, as is usually done. Consequently, for each of the two compartments, the model included uninfected target cells (T) that can be infected (I₁) either by infectious viruses (V_i) or inoculum (V_s) at an infection rate β. After an eclipse phase, infected cells become productively infected cells (I₂) and can produce virions at rate P and be lost at a per capita rate δ. The virions generated can be infectious (V_i) with proportion µ while the (1−µ) remaining proportion of virions is non-infectious (V_ni). Mathematically, a single compartment (V) for de novo produced virions could be considered in the model, with µV and (1−µ)V representing the respective contributions of infectious and non-infectious viruses to the biological mechanisms. However, to have a better visual understanding of the distinction between the two types of viruses, we wrote the model with distinct compartments, V_i and V_ni.

Finally, virions produced de novo and those from the inoculum are cleared at a rate c and c_i, respectively. Distinct clearances were considered to account for the effects of experimental conditions on viral dynamics. In particular, it is hypothesized that, animals being locally infected with large numbers of virions, a large proportion of it is assumed to be rapidly eliminated by swallowing and natural downstream influx, in contrast to the de novo-produced virions. However, it is important to keep in mind that this distinction was possible because of the controlled experimental conditions performed in animals, (i.e., exact timing and amount of inoculated virus known, and frequent monitoring during the early phase of the viral dynamics). Because of identifiability issues, similar clearances for infectious and non-infectious viruses were used. Accordingly, the model can be written as the following set of differential equations, where the superscript X denotes the compartment of interest (N, nasopharynx or T, trachea):

{\begin{cases} \frac{d T^{X}}{d t} = - β^{X} V_{i}^{X} T^{X} - μ β^{X} V_{s}^{X} T^{X} \frac{d I_{1}^{X}}{d t} = β^{X} V_{i}^{X} T^{X} + μ β^{X} V_{s}^{X} T^{X} - k I_{1}^{X} \frac{d I_{2}^{X}}{d t} = k I_{1}^{X} - δ^{X} I_{2}^{X} \frac{d V_{i}^{X}}{d t} = μ P^{X} I_{2}^{X} - c V_{i}^{X} - β^{X} V_{i}^{X} T^{X} \frac{d V_{n i}^{X}}{d t} = (1 - μ) P^{X} I_{2}^{X} - c V_{n i}^{X} \frac{d V_{s}^{X}}{d t} = - c_{i} V_{s}^{X} - μ β^{X} V_{s}^{X} T^{X} T^{X} (t = 0) = T_{0}^{X}; I_{1}^{X} (t = 0) = 0; I_{2}^{X} (t = 0) = 0 V_{i}^{X} (t = 0) = 0; V_{n i}^{X} (t = 0) = 0; V_{s}^{X} (t = 0) = V_{S, 0}^{X} \end{cases}

where $T^{X} (t = 0)$ , $I_{1}^{X} (t = 0)$ , $I_{2}^{X} (t = 0)$ , $V_{i}^{X} (t = 0)$ , $V_{n i}^{X} (t = 0)$ , and $V_{s}^{X} (t = 0)$ are the initial conditions at the time of exposure. The initial concentration of target cells, that are the epithelial cells expressing the ACE2 receptor, is expressed as $T_{0}^{X} = \frac{T_{0}^{X, n b c}}{W^{X}}$ , where $T_{0}^{X, n b c}$ is the initial number of cells and W^X is the volume of distribution of the compartment of interest (see the subsection ‘Consideration of the volume of distribution’). Each animal was exposed to 1 × 10⁶ pfu of SARS-CoV-2 representing a total of 2.19 × 10¹⁰ virions. Over the total inoculum injected (5 mL), 10% (0.5 mL) and 90% (4.5 mL) of virions were respectively injected by the intra-nasal route and the intra-tracheal route leading to the following initial concentrations of the inoculum within each compartment: $V_{S, 0}^{N} = \frac{0.10 \times I n o c_{0}}{W^{N}}$ and $V_{S, 0}^{T} = \frac{0.90 \times {I n o c}_{0}}{W^{T}}$ , with Inoc₀ the number of virions injected via the inoculum.

Using the gRNA and sgRNA viral loads, we estimated the viral infectivity, the viral production rate, and the loss rate of infected cells within each of the two compartments of the URT (Supplementary file 2). To account for interindividual variability and covariates, each of those three parameters was described by a mixed-effect model and jointly estimated between the two compartments as follows:

\begin{array}{ll} {\begin{cases} \log_{10} (β_{i}^{N}) = β_{0} + ϕ_{c o n v}^{β} \times I_{g r o u p = c o n v} + ϕ_{C D 40}^{β} \times I_{g r o u p = C D 40} + u_{i}^{β} \\ β_{i}^{T} = β_{i}^{N} \times \exp (f_{β}^{T}) \\ \log (δ_{i}^{N}) = \log (δ_{0}) + ϕ_{c o n v}^{δ} \times I_{g r o u p = c o n v} + ϕ_{C D 40}^{δ} \times I_{g r o u p = C D 40} + u_{i}^{δ} \\ δ_{i}^{T} = δ_{i}^{N} \times \exp (f_{δ}^{T}) \\ \log (P_{i}^{N}) = \log (P_{0}) + ϕ_{c o n v}^{P} \times I_{g r o u p = c o n v} + ϕ_{C D 40}^{P} \times I_{g r o u p = C D 40} + u_{i}^{P} \\ P_{i}^{T} = P_{i}^{N} \times \exp (f_{P}^{T}) \end{cases} \end{array}

where $β_{0}, \log (δ_{0})$ , and $\log (P_{0})$ are the fixed effects, ${ϕ_{c o n v}^{θ} ∣ θ \in {β, δ, P}}$ and ${ϕ_{C D 40}^{θ} ∣ θ \in {β, δ, P}}$ are respectively the regression coefficients related to the effects of the group of convalescent and αCD40.RBD-vaccinated animals for the parameters β, δ, and P, and $u_{i}^{θ}$ is the individual random effect for the parameter θ, which is assumed to be normally distributed with variance $ω_{θ}^{2}$ . A log-transformation was adopted for the parameters δ and P to ensure their positivity while a log₁₀-transformation was chosen for viral infectivity to also improve the convergence of the estimation. Because of the scale difference between the parameter β and the other parameters (see Supplementary file 2), the mere use of the log-transformation for this parameter led to convergence issues. The use of a log₁₀-transformation allowed to overcome this problem. Moreover, as shown in Equation 2, a joint estimation of the parameters β, δ, and P between the two compartments of the URT was considered. In this regard, a homogeneous interindividual variability within the URT was assumed as well as a similar contribution of the covariates to the value of the parameters. Parameters in the trachea were then either equal or proportional to those in the nasopharynx. This modelling choice, resulting in a smaller number of parameters to be estimated, was made mainly to address identifiability issues and to increase the power of the estimation. All other parameters included in the target-cell limited models were assumed to be fixed (see the subsection ‘Parameter estimation’ for more details).

In practice, after the selection of the optimal statistical model (see Appendix 1 ‘Model building’), random effects were added only to the parameters β and δ (i.e., ω_β ≠0, ω_δ ≠0, and ω_P=0), and the estimation of multiple models identified the viral production rate P as the only parameter taking different values between the trachea and nasopharynx. (i.e., β^N=β^T with f_β^T=0, δ^N= δ^T with f_δ^T=0, while P^N≠P^T). Finally, the adjustment of the model for the categorical covariates of groups of treatment, natural infection, and/or vaccination identified β and δ as the parameters with a statistically significant effect of these covariates (i.e., $ϕ_{c o n v}^{P} = 0$ and $ϕ_{C D 40}^{P} = 0$ ).

For the observation model, we jointly described genomic and subgenomic viral loads in the two compartments of the URT. We defined genomic viral load, which characterizes the total viral load observed in a compartment (nasopharynx or trachea), as the sum of inoculated virions (V_s), infectious (V_i), and non-infectious virions (V_ni). The sgRNA was described as proportional to the infected cells (I₁ + I₂). This choice was driven by two main reasons. First, sgRNA is only transcribed in infected cells (Sawicki et al., 2007). Second, as described by Miao et al., 2011, to overcome identifiability issues between the parameters β and P typically observed in target-cell limited models. The comparison of the two observation models describing sgRNA as either proportional to virions produced de novo (V_i +V_ni) or proportional to infected cells (I₁ + I₂) confirmed this conclusion. In addition to a better BICc value (–25 points) compared with the first model, the second one allowed the estimation of both β and P by counteracting identifiability problems faced with the first model (results not shown). Accordingly, the log₁₀-transformed gRNA and sgRNA of the ith animal at the jth time point in compartment X (nasopharynx or trachea), denoted ${g R N A}_{i j}^{X}$ and ${s g R N A}_{i j}^{X}$ , respectively, were described by the following equations:

{\begin{cases} g R N A_{i j}^{X} = \log_{10} [(V_{i}^{X} + V_{n i}^{X} + V_{s}^{X}) (Θ_{i}^{X}, t_{i j})] + ε_{i j, g}^{X} ε_{i j, g}^{X} \sim N (0, σ_{g X}^{2}) \\ s g R N A_{i j}^{X} = α_{s g R N A} \times \log_{10} [(I_{1}^{X} + I_{2}^{X}) (Θ_{i}^{X}, t_{i j})] + ε_{i j, s g}^{X} ε_{i j, s g}^{X} \sim N (0, σ_{s g X}^{2}) \end{cases}

where $Θ_{i}^{X}$ is the set of parameters of the subject i for the compartment X and ε are the additive normally distributed measurement errors.

Consideration of the volume of distribution

Request a detailed protocol

To define the concentration of inoculum within each compartment after injection, nasopharyngeal and tracheal volumes of distribution, labelled W^N and W^T, respectively, were needed. Given the estimated volumes of the trachea and the nasal cavities in four monkeys similar to our 18 macaques (Figure 2—figure supplement 2A–C) and the well-documented relationship between the volume of respiratory tract and animal weights (Asgharian et al., 2012), the volume of distribution of each compartment was defined as a step function of NHP weights:

\begin{array}{ll} W_{i}^{N} = {\begin{matrix} 4 & {i f w e i g h t}_{i} \leq 4.5 \\ 5.5 & o t h e r w i s e \end{matrix} \\ W_{i}^{N} = {\begin{matrix} 2 & {i f w e i g h t}_{i} \leq 4.5 \\ 3 & o t h e r w i s e \end{matrix} \end{array}

where weight_i is the weight of the monkey i in kg. Using Equation 4 and weights of our 18 NHPs (mean = 4.08; [Q1; Q3] = [3.26; 4.77]), we estimated WT = 2 and WN = 4 mL for a third of them (n=12) (Figure 2—figure supplement 2D), leading to the initial concentration of target cells $T_{0}^{X}$ (see ‘Viral dynamics modelling’ for equation) fixed at 3.13 × 10⁴ cells mL^–1 and 1.13 × 10⁴ cells mL^–1 in nasopharynx and trachea, respectively. Similarly, their initial concentrations of challenge inoculum $V_{S, 0}^{X}$ were fixed at 5.48 × 10⁸ copies,mL^–1 and 9.86 × 10⁹ copies,mL^–1 in nasopharynx and trachea respectively. For the last third of NHPs (n=6), WT = 3 and WN = 5.5 mL leading to $T_{0}^{X}$ fixed at 2.27 × 10⁴ cells mL^–1 in nasopharynx and 7.50 × 10³ cells mL^–1 in trachea while $V_{S, 0}^{X}$ was fixed at 3.98 × 10⁸ copies mL^–1 in nasopharynx and 6.57 × 10⁹ copies mL^–1 in trachea. Through this modelling, we assumed a homogenous distribution of injected virions and target cells within nasopharyngeal and tracheal compartments. In addition, the natural downward flow of inoculum toward lungs, at the moment of injection, was indirectly taken into account by the parameter of inoculum clearance, c_i.

Parameter estimation

Request a detailed protocol

Among all parameters involved in the three layers of the mechanistic model, some of them have been fixed based on experimental settings and/or literature. That is the case of the proportion of infectious virus (µ) that has been fixed at 1/1000 according to previous work (Gonçalves et al., 2021) and additional work (results not shown) evaluating the stability of the model estimation according to the value of this parameter. The initial number of target cells, that are the epithelial cells expressing the ACE2 receptor, $T_{0}^{X, n b c}$ was fixed at 1.25 × 10⁵ cells in the nasopharynx and 2.25 × 10⁴ cells in trachea (Gonçalves et al., 2021; Supplementary file 2). The duration of the eclipse phase (1/k), the clearance of the inoculum ( $c_{i}$ ) and the clearance of the virus produced de novo (c) were estimated by profile likelihood. The profile likelihood consists in defining a grid of values for the parameters to be evaluated and sequentially fixing these parameters to one of these combinations of values. The model and all the parameters that are not fixed are then estimated by maximizing the log-likelihood. In this process, all parameters that are assumed to be fixed in the model (i.e., μ and the initial conditions) are held fixed. Finally, the optimal set of parameters is chosen as the one optimizing the log-likelihood. Although the available data did not allow the direct estimation of these three parameters, the use profile likelihood enabled the exploration of various potential values for k, c, and $c_{i}$ . In a first step, we explored the 18 models resulting from the combination of three values of k∈{1, 3, 6} day^–1 and six values for c∈{1, 5, 10, 15, 20, 30} day^–1, assuming that the two parameters of virus clearance were equal, as first approximation. As shown in Supplementary file 3a, an eclipse phase of 8 hr (k=3) and virus clearance higher than 15 virions per day led to lowest values of –2log-likelihood (–2LL, the lower the better). In a second step, we fixed the parameter k at 3 day^–1 and estimated the 70 models resulting from the combination of 10 values for c∈{1, 2, 3, 4, 5, 10, 15, 20, 25, 30} day^–1 and 7 values for∈{1, 5, 10, 15, 20, 25, 30} day^–1 (Supplementary file 3b). The distinction of the two parameters of free virus clearance enabled to find much lower half-life of inoculum (~50 min) than half-life of virus produced de novo (~5.55 hr), with c=3 day^–1 compared to c_i = 20 day^–1.

Once all these parameters have been fixed, the estimation problem was restricted to the determination of the viral infectivity β, the viral production rate P, the loss rate of infected cells δ for each compartment, the parameter $α_{v l s g}$ in the observation model, regression coefficients for groups of intervention $(ϕ_{c o n v}, ϕ_{C D 40})$ , and standard deviations for both random effects ( $ω$ ) and error model (σ). The estimation was performed by maximum likelihood estimation using a stochastic approximation EM algorithm implemented in the software Monolix (http://www.lixoft.com). The Fisher information matrix was calculated by stochastic approximation, providing for each estimated parameter its variance, from which we were able to derive its 95% confidence interval. Selection of the compartment effect on parameters (β, δ, P) as well as random effects and covariates on the statistical model (Equation 2) was performed by the estimation of several models that were successively compared according to the corrected Bayesian information criterion (BICc) (to be minimized). After the removal of random effect on the viral production ( $ω_{P} = 0$ ) allowing the reduction of the variance on the two other random effects, all combinations of compartment effects were evaluated, leading to the final selection of a single effect on P $(f_{β}^{T} = f_{δ}^{T} = 0)$ . Then, the effect of group intervention was independently added on model parameters among β, δ, P, and c. Once the group effect on the viral infectivity identified as the best one, the addition of a second effect on the remaining parameters was tested, resulting in the selection of the loss rate of infected cells. Finally, the irrelevance of the addition of a third effect was verified.

The possibility of migration of free plasma virus between the nasopharynx and the trachea was tested. However, as widely described in the literature, the transport of viral particles within the respiratory tract is negligible in the viral dynamics and is difficult to estimate. The reader can refer to Appendix 1 ‘Model building’ for an additional modelling work conducted to estimate this exchange and provided the same conclusion. Accordingly, the two compartments of the URT were assumed are distinct in our model.

Algorithm for automatic selection of biomarkers as CoP

Request a detailed protocol

After identifying the effect of the group of intervention on both the viral infectivity (β) and the loss rate of infected cells (δ), we aimed at determining whether some immunological markers quantified in the study could capture this effect. Nowadays, many methods for selecting constant covariates already exist (Chowdhury and Turin, 2020) and are implemented in software like Monolix. However, these latter do not allow time-varying covariates. In this section, we present the algorithm we implemented to select time-varying covariates. We proposed a classical stepwise data-driven automatic covariate modelling method (Figure 4—figure supplement 2). However, initially implemented to select covariates from more than 50 biomarkers, computational time restricted us to consider only a forward selection procedure. Nevertheless, the method can be easily extended to classical stepwise selection in which both forward selection and backward elimination are performed sequentially. Although the method was developed for time-varying covariates, it can also be applied to constant covariates.

At the initialization step (k=0) (see Figure 4—figure supplement 2), the algorithm requests three inputs: (World Health Organization, 2021) a set of potential $M$ covariates, labelled Marker m for $m \in {1, \dots, M}$ (e.g., immunological markers); (Cobey et al., 2021) a set of P parameters on which covariates could be added, labelled θ_p for $p \in {1, \dots, P}$ (e.g., β and δ); and (Kuzmina et al., 2021) an initial model (e.g., the model without covariates), labelled M⁰, with $θ_{p}^{0}$ being the definition of the parameter θ_p. At each step k>0, we note M^k⁻¹ the current model resulting in the model built in the step k−1. Then, each combination of markers and parameters that have not already been added in M^k⁻¹, labelled $r (r \in {M a r k e r m ⨂$ $θ_{p} \notin M^{k - 1} ∣$ $m \in {1, \dots M}, p \in {1, \dots P}})$ , are considered and tested in an univariate manner (each relation r is independently added in M^k⁻¹ and ran). To this end, the parameter θ_p involved in this relationship r is modified as $θ_{p}^{k} (t) = θ_{p}^{k - 1} (t) \times e x p (ϕ_{m}^{p} \times M a r k e r_{m} (t))$ , where $ϕ_{m}^{p}$ is the regression coefficient related to the marker and $M a r k e r_{m} (t)$ being the trajectory of the marker over time, while other parameters remain unchanged $(\forall θ_{q} \notin r, θ_{q}^{k} (t) = θ_{q}^{k - 1} (t))$ . Once all these models evaluated, the one with the optimal value of a given selection criterion defining the quality of the fits (e.g., the lowest BICc value) is selected and compared to the model M^k⁻¹. If the value of the criterion is better than the one found for M^k⁻¹, then this model is defined as the new current model, M^k, and the algorithm moves to the step k+1. Otherwise, the algorithm stops. The algorithm can also be stopped at the end of a fixed number of step k.

The objective of this algorithm being to identify mechanistic CoP, at each step, the selected model should respect, in addition to the best fits criterion, the two other criteria defining mCoP meaning the ability to capture the effect of the group of intervention and the ability to better explain the variability on individual parameters than the model adjusted for the group effect. To this end, we verify that in the selected model additionally adjusted for the group of intervention, the group effect appears as non-significantly different from 0 using a Wald test. Then, we check that the variances of random effects in the selected model are lower or equal to the ones obtained in the model adjusted only for the group effect.

Modelling hypothesis for time-dependent covariates in our application

Request a detailed protocol

Using a population-based approach to estimate our mechanistic model and similar to the adjustment of the model for constant covariates (e.g., groups of intervention), time-varying covariates are incorporated into the statistical model as individual-specific explanatory variables in the mixed-effects models. To implement the algorithm for selecting the time-varying covariates, many modelling choices were made. First, targeting covariates able to fully replace the group of intervention, we kept a similar mathematical relationship between parameters and immune markers than the one used with the constant covariate (see Equation 2). Accordingly, we adjusted the model parameters additively in logarithmic scale. In this regard, at each step k (k>0), the parameter θ_p was defined as $\log (θ_{p}^{k} (t))$ $= \log (θ_{p}^{k - 1} (t))$ $+ ϕ_{m}^{p} \times M a r k e r_{m} (t)$ . However, this choice may affect the results and other choices may be more relevant under different conditions. Second, because immune markers are observed only at discrete time points, whereas the estimation of the model is performed in a continuous way, we introduced immune markers as time-varying covariates using linear interpolation. Lets denote Marker_i,j the value of the marker observed for the ith animal at the jth time point, with $i \in {1, \dots, n}$ and $j \in {1, \dots, J}$ . By linear interpolation, the time-continuous marker was defined as, $\forall t > 0$ ,

\begin{array}{ll} M a r k e r_{i}^{i n t} (t) = & \sum_{j = 1}^{J - 1} I_{[t_{j}; t_{j + 1}]} (t) [\frac{M a r k e r_{i, j + 1} - M a r k e r_{i, j}}{t_{j + 1} - t_{j}} t + \frac{M a r k e r_{i, j} t_{j + 1} - M a r k e r_{j + 1} t_{j}}{t_{j + 1} - t_{j}}] + I_{t \geq t_{J}} (t) \times M a r k e r_{i, J} \end{array}

As previously described in the Results section, three different studies were considered in this work: a main study reported by Marlin et al., 2021, testing the αCD40.RBD vaccine, and two additional studies (Corbett et al., 2020; Brouwer et al., 2021) evaluating a two-component spike nanoparticle vaccine and the mRN-1273 vaccine, respectively. In the main study, the method was applied with both time-varying covariates and constant covariates for which only baseline value was considered, such that Marker_i(t)=Marker_i(t=0) (see Supplementary file 1). For the other two studies, only the baseline values were considered as covariates, the dynamics being not available. To assess the robustness of the results, several selection criteria were tested: AIC, BIC, log-likelihood, the percentage of explained interindividual variability, and similar results were obtained for all (results not shown). Moreover, as presented in Appendix 2 ‘BICc as selection criteria and multiple testing adjustment’, we verified the robustness of the use of BIC as selection criteria despite the multiplicity of the tests. The identification of antibodies inhibiting the attachment of the RBD to the ACE2 receptor (ECLRBD) as the first time-varying CoP led to the definition of the time-varying viral infectivity for the ith animal as described in Equation 5, while the selection anti-RBD IgG-binding antibodies led to the elimination rate of infected cells given in Equation 6.

\begin{matrix} β_{i} (t) = 1 0^{β_{0} + u_{i}^{β}} \times e x p (ϕ_{e c l}^{β} \times E C L R B D_{i}^{i n t} (t)) \end{matrix}

δ_{i} (t) = δ_{0} \times e x p (ϕ_{i g g}^{δ} \times I g g R B D_{i}^{i n t} (t) + u_{i}^{δ})

Quantification and statistical analysis

Request a detailed protocol

In each of the three studies used in this work, no statistical tests were performed on the raw data (i.e, observations), whether for viral load or for immune marker measurements, to identify statistical differences between treatment groups, as the statistical analyses were already been performed in the respective papers. Statistical significance of the effect of groups in model estimation is indicated in the tables by stars: *, p<0.05; **, p<0.01; ***, p<0.001 and were estimated by Wald tests (Monolix software version 2019R1).

Model parameters were estimated with the SAEM algorithm (Monolix software version 2019R1). Graphics were generated using R version 3.6.1 and Excel 2016 and details on the statistical analysis for the experiments can be found in the accompanying figure legends. Horizontal red dashed lines on graphs indicate assay limit of detection.

Appendix 1

Model building

In the model presented in the manuscript, we considered the two compartments of the URT, trachea, and nasopharynx, as two distinct compartments (i.e., without transfer of virus between them), as described by Equation AE1. In each of them, the viral dynamics are described by a target-cell limited model augmented with a compartment describing the dynamics of the inoculated virus (V_s). Moreover, in the statistical model describing the model parameters, the three parameters β, δ, and P were assumed as jointly estimated between the two compartments, with shared random effects and covariates and considering that parameters β and δ are equal in both trachea and nasopharynx (β^T = β^N, δ^T = δ^N).

\begin{array}{ll} {\begin{cases} \frac{d T^{N}}{d t} = - β^{N} V_{i}^{N} T^{N} - μ β^{N} V_{s}^{N} T^{N} \\ \frac{d I_{1}^{N}}{d t} = β^{N} V_{i}^{N} T^{N} + μ β^{N} V_{s}^{N} T^{N} - k I_{1}^{N} \\ \frac{d I_{2}^{N}}{d t} = k I_{1}^{N} - δ^{N} I_{2}^{N} \\ \frac{d V_{i}^{N}}{d t} = μ P^{N} I_{2}^{N} - c V_{i}^{N} - β^{N} V_{i}^{N} T^{N} \\ \frac{d V_{n i}^{N}}{d t} = (1 - μ) P^{N} I_{2}^{N} - c V_{n i}^{N} \\ \frac{d V_{s}^{N}}{d t} = - c_{i} V_{s}^{N} - μ β^{N} V_{s}^{N} T^{N} \end{cases} {\begin{cases} \frac{d T^{T}}{d t} = - β^{T} V_{i}^{T} T^{T} - μ β^{T} V_{s}^{T} T^{T} \\ \frac{d I_{1}^{T}}{d t} = β^{T} V_{i}^{T} T^{T} + μ β^{T} V_{s}^{T} T^{T} - k I_{1}^{T} \\ \frac{d I_{2}^{T}}{d t} = k I_{1}^{T} - δ^{T} I_{2}^{T} \\ \frac{d V_{i}^{T}}{d t} = μ P^{T} I_{2}^{T} - c V_{i}^{T} - β^{T} V_{i}^{T} T^{T} \\ \frac{d V_{n i}^{T}}{d t} = (1 - μ) P^{T} I_{2}^{T} - c V_{n i}^{T} \\ \frac{d V_{s}^{T}}{d t} = - c_{i} V_{s}^{T} - μ β^{T} V_{s}^{T} T^{T} \end{cases} \end{array}

Initially, random effects were added on the three parameters. However, taken into consideration identifiability issues that are usually encountered between the viral infectivity (β) and the viral production (P), we decided to remove the possibility of interindividual variability on the parameter P. This choice was also driven by multiple model estimations showing less robust estimations when variability was allowed in both parameters β and P. In particular, the estimate of the viral production was impacted by a ratio between the parameter and its standard error (RSE) higher than 100%.

Comparison of the parameters between the tracheal and the nasopharyngeal compartments

To decide which of these three parameters were assumed to be equal between the two compartments, all possibilities were tested and compared, using the BICc as selection criteria. As shown in Appendix 1—table 1, we started with the model in which all parameters were equal between the two compartments and we progressively relaxed this hypothesis. During this step, no exchange of virions between the two compartments of the URT was possible (g=0). Once all models estimated, we kept the one with the lowest value of BICc, meaning with the highest negative difference of BICc compared to the initial model. We identified the model with only the viral production varying between the two compartments as the best one to fit the data.

Appendix 1—table 1

Comparison of models evaluating the difference of viral infectivity (β), loss of infected cells (δ), and viral production (P) between the nasopharynx and the trachea.

Model tested	Statistical model	ΔBICc
Initial model	β^T = β^N δ^N = δ^T P^N = P^T Variability on β and δ
Model with different β	β^T ≠ β^N δ^N = δ^T P^N = P^T Variability on β and δ	–17.31
Model with different δ	β^T = β^N δ^N ≠ δ^T P^N = P^T Variability on β and δ	–14.38
Model with different P	β^T = β^N δ^N = δ^T P^N ≠ P^T Variability on β and δ	–25.24
Model with different β and δ	β^T ≠ β^N δ^N ≠ δ^T P^N = P^T Variability on β and δ	–13.00
Model with different β and P	β^T ≠ β^N δ^N = δ^T P^N ≠ P^T Variability on β and δ	–19.19
Model with different δ and P	β^T = β^N δ^N ≠ δ^T P^N ≠ P^T Variability on β and δ	–19.47
Model with different β, δ, and P	β^T ≠ β^N δ^N ≠ δ^T P^N ≠ P^T Variability on β and δ	–13.39

Identification of group effects

Once the structure of the statistical model defined, we tried to identify on which parameters an effect of the group of treatment could be identified and by extension on which biological mechanisms. In this step, we were interested in four parameters: β, δ, P, and c, the latter being the clearance of de novo- produced virions. In the study, three groups of treatments were considered as constant categorical covariates: naïve, convalescent, and convalescent vaccinated. We performed a forward selection approach using the BICc as selection criteria to find the best model, using the model without covariate as initial model. At each step the model decreasing the most the value of the BICc is selected and the procedure stops once the BICc does not decrease anymore. At each step of the procedure, the statistical significance of covariate added into the model was verified via a Wald test. As shown in Appendix 1—table 2, the selected model identified a group effect on the viral infectivity and the loss rate of infected cells.

Appendix 1—table 2

Comparison of models evaluating the adjustment of the viral infectivity (β), the loss rate of infected cells (δ), the viral production (P), and the viral clearance (c) for the groups of treatment.

The group of naïve animals is assumed as the group of reference.

Step	Model tested	Statistical model	$Δ B I C c$
1	Initial model: Model without group effects	$β = 10^{β_{0}}$ $δ = δ_{0}$ $P = P_{0}$ $c = c_{0}$
	Model with group effect on $β$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0}$ $P = P_{0}$ $c = c_{0}$	–21.5
	Model with group effect on $δ$	$β = 10^{β_{0}}$ $δ = δ_{0} \exp (ϕ_{c o n v}^{δ} + ϕ_{C D 40}^{δ})$ $P = P_{0}$ $c = c_{0}$	–16.62
	Model with group effect on $P$	$β = 10^{β_{0}}$ $δ = δ_{0}$ $P = P_{0} \exp (ϕ_{c o n v}^{P} + ϕ_{C D 40}^{P})$ $c = c_{0}$	+9.68
	Model with group effect on c	$β = 10^{β_{0}}$ $δ = δ_{0}$ $P = P_{0}$ $c = c_{0} e x p (ϕ_{c o n v}^{c} + ϕ_{C D 40}^{c})$	+9.20
2	Initial model: Model with group effect on $β$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0}$ $P = P_{0}$ $c = c_{0}$
	Model with group effect on $β$ and $δ$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0} e x p (ϕ_{c o n v}^{δ} + ϕ_{C D 40}^{δ})$ $P = P_{0}$ $c = c_{0}$	–2.48
	Model with group effect on $β$ and $P$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0}$ $P = P_{0} \exp (ϕ_{c o n v}^{P} + ϕ_{C D 40}^{P})$ $c = c_{0}$	+12.25
	Model with group effect on $β$ and $c$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0}$ $P = P_{0}$ $c = c_{0} e x p (ϕ_{c o n v}^{c} + ϕ_{C D 40}^{c})$	+11.97
3	Initial model: Model with group effect on $β$ and $δ$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0} \exp (ϕ_{c o n v}^{δ} + ϕ_{C D 40}^{δ})$ $P = P_{0}$ $c = c_{0}$
	Model with group effect on $β$ , $δ$ , and $P$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0} \exp (ϕ_{c o n v}^{δ} + ϕ_{C D 40}^{δ})$ $P = P_{0} \exp (ϕ_{c o n v}^{P} + ϕ_{C D 40}^{P})$ $c = c_{0}$	+10.88
	Model with group effect on $β, δ$ , and $c$	$β = 10^{(β_{0} + ϕ_{c o n v}^{β} + ϕ_{C D 40}^{β})}$ $δ = δ_{0} \exp (ϕ_{c o n v}^{δ} + ϕ_{C D 40}^{δ})$ $P = P_{0}$ $c = c_{0} e x p (ϕ_{c o n v}^{c} + ϕ_{C D 40}^{c})$	+11.61

Based on all these results, the optimal statistical model with adjustment for groups of treatment was defined as follows:

{\begin{cases} \log_{10} (β_{i}) = β_{0} + ϕ_{c o n v}^{β} \times I_{i \in c o n v} + ϕ_{C D 40}^{β} \times I_{i \in C D 40} + u_{i}^{β} \\ \log (δ_{i}) = \log (δ_{0}) + ϕ_{c o n v}^{δ} \times I_{i \in c o n v} + ϕ_{C D 40}^{δ} \times I_{i \in C D 40} + u_{i}^{δ} \\ \log (P_{i}^{N}) = \log (P_{0}) \\ P_{i}^{T} = P_{i}^{N} \times \exp (f_{P}^{T}) \end{cases}

Exchange of viruses between the nasopharyngeal and tracheal compartments

Afterward, we tested the possibility of an exchange of free plasma virus from between the two compartments of the URT. We made the hypothesis of a constant first-order exchange and we tested the addition a transfer of virions from nasopharyngeal to tracheal compartments and vice versa, with a migration rate g_NT and g_TN, respectively. To this end, equations of infectious (V_i) and non-infectious (V_ni) viruses in Equation AE1 between the two compartments were linked as follows:

\begin{array}{cc} \frac{d V_{i}^{T}}{d t} \mapsto \frac{d V_{i}^{T}}{d t} - g_{T N} V_{i}^{T} + g_{N T} V_{i}^{N} & \frac{d V_{n i}^{T}}{d t} \mapsto \frac{d V_{n i}^{T}}{d t} - g_{T N} V_{n i}^{T} + g_{N T} V_{n i}^{N} \\ \frac{d V_{i}^{N}}{d t} \mapsto \frac{d V_{i}^{N}}{d t} + g_{T N} V_{i}^{T} - g_{N T} V_{i}^{N} & \frac{d V_{n i}^{N}}{d t} \mapsto \frac{d V_{n i}^{N}}{d t} + g_{T N} V_{n i}^{T} - g_{N T} V_{n i}^{N} \end{array}

with the arrow symbolizing the modification of the equations defined in Equation AE1 and g_NT and g_TN being two positive rates. As a first step, we tried to estimate either bidirectional or one of the two unidirectional transfers using the data from the 18 NHPs of the first study described in the main paper. However, data were too spare to bring enough information to get estimations. Consequently, as a second step, additional data were used: two naïve macaques were exposed to the same dose (1 × 10⁶ pfu) of SARS-CoV-2 than the 18 NHPs of the main study. However, instead of being inoculated via intra-tracheal (4.5 mL) and intra-nasal (0.5 mL) routes, these latter received inoculum via intra-gastric (4.5 mL) and intra-nasal (0.5 mL) routes. Similar to the main study, the viral gRNA dynamics in both tracheal and nasopharyngeal compartments were repeatedly measured during the 20 days following the challenge (Figure 2—figure supplement 2E).

These two additional macaques having not received intra-tracheal inoculum, viral dynamics measured in this same compartment was expected to come from (at least partially) an exchange with the nasopharynx and thus bring information about it. However, having only two macaques without virions inoculated via intra-tracheal route, no enough information were available to totally estimate the model with exchanges. Consequently, these two additional NHPs having similar characteristics than the 18 NHPs involved in the main study, we made the assumption that the viral dynamics in nasopharynx after inoculation and the viral dynamics in the trachea, once the transfer initiated, should be described by the same model (without inoculum in trachea) and those by the same parameters. We expected that the difference of dynamics in trachea between these two set of macaques could allow an estimation of the parameters g_TN and/or g_NT. For that reason, we estimated the model in Equation AE1 using data from the 18 NHPs of the main study. Then using the data from the two additional NHPs, and assuming all parameters of the model resulting from Equation AE2 as fixed (see Appendix 1—table 2), except g_TN and g_NT, we tried to quantify the transfers of virions.

The estimation of multiple models on those two animals tended to conclude that only a unidirectional transfer of viruses from the nasopharyngeal to the tracheal compartment should be explored, with an estimation of g_NT ranging from 0.9 to 2.5 day^–1. Once these values quantified, we tried to update/re-estimate the model, initially estimated on the 18 NHPs, using only a unidirectional transfer from nasopharynx to trachea and fixing the value of the migration rate at the different values aforementioned. However, all tested values of g_NT led irremediably to a degradation of the model with an increase of at least two points of BICc.

An estimation of the parameter g_NT by profile likelihood (results not shown) led to a strictly increasing profile of the likelihood (the lower the better) and was thus no more conclusive. Consequently, no exchange of virions were assumed in the final model and the parameters g_NT and g_TN were fixed at 0 day^–1.

Appendix 2

BICc as selection criteria and multiple testing adjustment

In the case of classic covariate selection approaches using p-values as selection criteria, particular attention must be paid to take into account the dependence of the results on the number tests performed.

Over the years, multiple corrections have been proposed to adjust results for test multiplicity (e.g., Bonferroni correction, Benjamini and Hochberg correction among others).

Although we verified the significance of the covariate selected in our model, our covariate selection approach relies on the BICc. To ensure the robustness of the BICc as selection criterion despite the multiplicity of the tests, we performed an additional simulation work.

We simulated M=25 longitudinal variables for 18 individuals and with similar time points than those found on our data, meaning at days 0, 4, 9, and 20 post-infection. Variables were simulated as white-noise random variables such that for the ith subject at the jth time point, the mth variable was defined as $X_{i j}^{m} \sim N (0, σ^{2})$ , with m=1, …, M. In our simulations, we tested five values for the variance σ² ranging from 1% to 10% (five variables simulated for each value of σ).

Assuming these variables as our time-varying covariates, we applied the forward selection approach used in our method by testing each of them in a univariate manner of both β and δ.

As shown in Appendix 2—figure 1, the 50 models built to evaluate the adjustment of either β or δ for the simulated variables provide similar results in terms of BICc, and thus whatever the value of the standard deviation σ used. Consequently, these results appear as quite robust to the multiplicity of the test. Moreover, as expected, adjustments for white-noise random variables depict the degradation of the model in comparison to the model without covariates.

Appendix 2—figure 1

Download asset Open asset

Results of the forward selection approach applied on the 25 simulated white-noise random variables.

The discrete x-axis represents the different variables and the y-axis represents the values of the corrected Bayesian information criteria (BICc). Circles and triangles correspond to the results obtained with the parameters β or δ adjusted for the variables. The horizontal solid black line represents the value of the BICc obtained with the model without covariates while the horizontal dashed green line highlights the value of the criterion obtained with both β and δ adjusted for the groups of treatment.

Evaluation of the robustness of the estimation

To evaluate the robustness of the parameter estimates obtained on our models, despite the small number of independent observations, we performed a bootstrap procedure with replacement (Thai et al., 2014), for B=50 iterations. The bootstrap parameter estimate was calculated as the median of the parameter estimates from the B bootstrap samples while the standard error of each parameter was calculated according to the definition of Thai et al., 2014, which means with the SE of the lth component of the vector of parameters given by:

{\hat{S E}}_{B}^{l} = \sqrt{\frac{1}{B - 1} \sum_{b = 1}^{B} {(\hat{θ}_{b}^{* (l)} - {\hat{θ}}_{B}^{(l)})}^{2}}

with ${\hat{θ}}_{b}^{* (l)}$ being its estimate obtained at the bth iteration of the bootstrap and ${\hat{θ}}_{B}^{(l)}$ the bootstrap parameter estimate. For each bootstrap sample, we paid attention to keep the 1:1:1 ratio between the three groups of treatment, with six animals selected within each group. Results are reported in Appendix 2—table 1 and Appendix 2—table 2.

Appendix 2—table 1

Model parameters for viral dynamics in both the nasopharynx and the trachea estimated by the model adjusted for groups of intervention.

For the bootstrap procedure, 50 iterations were performed.

Parameter	Meaning	Value [95% CI]	Unit
β	Viral infectivity in the naive group (×10^–6)	0.91 [0.12; 7.03]	(copies/mL)^–1 day^–1
	Fold change in the convalescent group	0.15 [0.04; 0.58]
	Fold change in the Conv-CD40 group	0.006 [0.001; 0.04]
δ	Loss rate of infected cells in the naive group	1.09 [0.74; 1.60]	day^–1
	Fold change in the convalescent group	1.70 [1.08; 2.66]
	Fold change in the Conv-CD40 group	2.00 [0.94; 4.27]
P^N	Viral production rate in the naso. (×10³)	10.1 [1.16; 87.7]	virions (cell day)^–1
P^T	Viral production rate in the trachea (×10³)	0.86 [0.08; 9.19]	virions (cell day)^–1
α_vlsg	Infected cells and sgRNA viral load ratio	1.42 [0.99; 2.02]	virions cell^–1
k	Eclipse rate	3	day^–1
c	Clearance of de novo produced viruses	3	day^–1
c_I	Clearance of inoculum	20	day^–1
µ	Percentage of infectious viruses	10^–3
$T_{0}^{X, n b c}$	Initial number of target cells	1.25 × 10⁵ (naso.) 2.25 × 10⁴ (trachea)	cells
Inoc₀	Number of virions inoculated	2.19 × 10¹⁰	virions
ω_β	SD of random effect on log₁₀ β	0.319 [0.111; 0.527]
ω_δ	SD of random effect on δ	0.122 [-0.039; 0.283]
σ_VLn	SD of error model gRNA in naso.	1.24 [0.96; 1.51]
σ_VLt	SD of error model gRNA in trachea	1.09 [0.92; 1.26]
σ_sgVLn	SD of error model sgRNA in naso	1.35 [1.08; 1.61]
σ_sgVLt	SD of error model sgRNA in trachea	1.53 [1.15; 1.92]

Appendix 2—table 2

Model parameters for viral dynamics in both the nasopharynx and the trachea estimated by the model with the viral infectivity adjusted for ACE2-RBD-binding inhibition and the loss rate of infected cells adjusted for the group of treatment.

For the bootstrap procedure, 50 iterations were performed.

Parameters	Meaning	Value [95% CI]	Unit
β	Infection rate with ECLRBD = 0 AU (×10^–8)	0.82 [0.13; 5.13]	(copies/mL)^–1 day^–1
β	Fold $Δ E C L R B D = 10^{3} A U$	1.017 [1.012; 1.022]
δ	Loss rate of infected cells	1.02 [0.80; 1.30]	day^–1
	Fold change in the convalescent group	1.74 [1.24; 2.46]
	Fold change in the Conv-CD40 group	2.17 [0.82; 5.74]
P^N	Viral production rate in the naso. (×10³)	8.92 [0.42; 191]	virions (cell day)^–1
P^T	Viral production rate in the trachea (×10³)	0.62 [0.02; 19.7]	virions (cell day)^–1
α_vlsg	Infected cells and sgRNA viral load ratio	1.32 [0.91; 1.90]	virions cell^–1
k	Eclipse rate	3	day^–1
c	Clearance of de novo produced viruses	3	day^–1
c_I	Clearance of inoculum	20	day^–1
µ	Percentage of infectious viruses	10^–3
$T_{0}^{X, n b c}$	Initial number of target cells	1.25 × 10⁵ (naso.) 2.25 × 10⁴ (trachea)	cells
Inoc₀	Number of virions inoculated	2.19 × 10¹⁰	virions
ω_β	SD of random effect on log₁₀ β	0.205 [0.011; 0.399]
ω_δ	SD of random effect on δ	0.079 [-0.092; 0.250]
σ_VLn	SD of error model gRNA in naso.	1.13 [0.90; 1.36]
σ_VLt	SD of error model gRNA in trachea	1.27 [1.07; 1.48]
σ_sgVLn	SD of error model sgRNA in naso	1.62 [1.30; 1.94]
σ_sgVLt	SD of error model sgRNA in trachea	1.36 [1.15; 1.56]

Data availability

No unique reagents were generated for this study. Data that support the findings of this study are provided in the source data files of this paper and gather data from (1) the study [Marlin, Nature Com 2021] used in this analysis, which are also directly available online in the section Source data of this related paper (https://www.nature.com/articles/s41467-021-25382-0#Sec17); (2) the study [Brouwer, Cell 2021] used in this analysis, which are also available from the corresponding authors of the related paper and (3) the study [Corbett, NEJM 2020] used in this analysis, which are also available online in the section Supplementary Material of the related paper, excel file labelled ("Supplementary Appendix 2"). Data from the main study [Marlin, Nature Com 2021] can also be found in the open-access repository Dryad using the following DOI: https://doi.org/10.5061/dryad.1zcrjdfv7. The original code (mlxtran models and R) as well as model definition files including the full list of parameters used are available and free-of-cost on github (Inria SISTM Team) at the following link: https://github.com/sistm/SARSCoV2modelingNHP, (copy archived at swh:1:rev:a704c80daebc949434694d3f4441e48293c461cc).

The following data sets were generated

1. Alexandre M
2. Marlin R
3. Prague M
4. Coleon S
5. Kahlaoui N
6. Cardinaud S
7. Naninck T
8. Delache B
9. Surenaud M
10. Galhaut M
11. Dereuddre-Bosquet N
12. Cavarelli M
13. Maisonnasse P
14. Centlivre M
15. Lacabaratz C
16. Wiedemann A
17. Zurawski S
18. Zurawski G
19. Schwartz O
20. Sanders RW
21. Le Grand R
22. Levy Y
23. Thiébaut R
(2022) Dryad Digital Repository
Viral loads and antibody, cytokine and T-cell responses in NHPs following vaccination targeting SARS-CoV-2 RBD domain to cells expressing CD40.

https://doi.org/10.5061/dryad.1zcrjdfv7

References

1. Alter G
2. Yu J
3. Liu J
4. Chandrashekar A
5. Borducchi EN
6. Tostanoski LH
7. McMahan K
8. Jacob-Dolan C
9. Martinez DR
10. Chang A
11. Anioke T
12. Lifton M
13. Nkolola J
14. Stephenson KE
15. Atyeo C
16. Shin S
17. Fields P
18. Kaplan I
19. Robins H
20. Amanat F
21. Krammer F
22. Baric RS
23. Le Gars M
24. Sadoff J
25. de Groot AM
26. Heerwegh D
27. Struyf F
28. Douoguih M
29. van Hoof J
30. Schuitemaker H
31. Barouch DH
(2021) Immunogenicity of Ad26.COV2.S vaccine against SARS-CoV-2 variants in humans
Nature 596:268–272.

https://doi.org/10.1038/s41586-021-03681-2
- PubMed
- Google Scholar
1. Asgharian B
2. Price O
3. McClellan G
4. Corley R
5. Einstein DR
6. Jacob RE
7. Harkema J
8. Carey SA
9. Schelegle E
10. Hyde D
11. Kimbell JS
12. Miller FJ
(2012) Development of a rhesus monkey lung geometry model and application to particle deposition in comparison to humans
Inhalation Toxicology 24:869–899.

https://doi.org/10.3109/08958378.2012.725782
- PubMed
- Google Scholar
(2006) Kinetics of influenza A virus infection in humans
Journal of Virology 80:7590–7599.

https://doi.org/10.1128/JVI.01623-05
- PubMed
- Google Scholar
1. Balelli I
2. Pasin C
3. Prague M
4. Crauste F
5. Effelterre TV
6. Bockstal V
7. Solforosi L
8. Thiébaut R
(2020) A model for establishment, maintenance and reactivation of the immune response after vaccination against Ebola virus
Journal of Theoretical Biology 495:110254.

https://doi.org/10.1016/j.jtbi.2020.110254
- PubMed
- Google Scholar
1. Bradfute SB
2. Bavari S
(2011) Correlates of immunity to filovirus infection
Viruses 3:982–1000.

https://doi.org/10.3390/v3070982
- PubMed
- Google Scholar
1. Brouwer PJM
2. Brinkkemper M
3. Maisonnasse P
4. Dereuddre-Bosquet N
5. Grobben M
6. Claireaux M
7. de Gast M
8. Marlin R
9. Chesnais V
10. Diry S
11. Allen JD
12. Watanabe Y
13. Giezen JM
14. Kerster G
15. Turner HL
16. van der Straten K
17. van der Linden CA
18. Aldon Y
19. Naninck T
20. Bontjer I
21. Burger JA
22. Poniman M
23. Mykytyn AZ
24. Okba NMA
25. Schermer EE
26. van Breemen MJ
27. Ravichandran R
28. Caniels TG
29. van Schooten J
30. Kahlaoui N
31. Contreras V
32. Lemaître J
33. Chapon C
34. Fang RHT
35. Villaudy J
36. Sliepen K
37. van der Velden YU
38. Haagmans BL
39. de Bree GJ
40. Ginoux E
41. Ward AB
42. Crispin M
43. King NP
44. van der Werf S
45. van Gils MJ
46. Le Grand R
47. Sanders RW
(2021) Two-component spike nanoparticle vaccine protects macaques from SARS-CoV-2 infection
Cell 184:1188–1200.

https://doi.org/10.1016/j.cell.2021.01.035
- PubMed
- Google Scholar
Book
(2006) Measurement Error in Nonlinear Models
Chapman and Hall/CRC.

https://doi.org/10.1201/9781420010138
- Google Scholar
1. Chandrashekar A
2. Liu J
3. Martinot AJ
4. McMahan K
5. Mercado NB
6. Peter L
7. Tostanoski LH
8. Yu J
9. Maliga Z
10. Nekorchuk M
11. Busman-Sahay K
12. Terry M
13. Wrijil LM
14. Ducat S
15. Martinez DR
16. Atyeo C
17. Fischinger S
18. Burke JS
19. Slein MD
20. Pessaint L
21. Van Ry A
22. Greenhouse J
23. Taylor T
24. Blade K
25. Cook A
26. Finneyfrock B
27. Brown R
28. Teow E
29. Velasco J
30. Zahn R
31. Wegmann F
32. Abbink P
33. Bondzie EA
34. Dagotto G
35. Gebre MS
36. He X
37. Jacob-Dolan C
38. Kordana N
39. Li Z
40. Lifton MA
41. Mahrokhian SH
42. Maxfield LF
43. Nityanandam R
44. Nkolola JP
45. Schmidt AG
46. Miller AD
47. Baric RS
48. Alter G
49. Sorger PK
50. Estes JD
51. Andersen H
52. Lewis MG
53. Barouch DH
(2020) SARS-CoV-2 infection protects against rechallenge in rhesus macaques
Science 369:812–817.

https://doi.org/10.1126/science.abc4776
- Google Scholar
1. Charmet T
2. Schaeffer L
3. Grant R
4. Galmiche S
5. Chény O
6. Von Platen C
7. Maurizot A
8. Rogoff A
9. Omar F
10. David C
11. Septfons A
12. Cauchemez S
13. Gaymard A
14. Lina B
15. Lefrancois LH
16. Enouf V
17. van der Werf S
18. Mailles A
19. Levy-Bruhl D
20. Carrat F
21. Fontanet A
(2021) Impact of original, B.1.1.7, and B.1.351/P.1 SARS-CoV-2 lineages on vaccine effectiveness of two doses of COVID-19 mRNA vaccines: Results from a nationwide case-control study in France
The Lancet Regional Health. Europe 8:100171.

https://doi.org/10.1016/j.lanepe.2021.100171
- PubMed
- Google Scholar
1. Cheng L
2. Wang Q
3. Li G
4. Banga R
5. Ma J
6. Yu H
7. Yasui F
8. Zhang Z
9. Pantaleo G
10. Perreau M
11. Zurawski S
12. Zurawski G
13. Levy Y
14. Su L
(2018) TLR3 agonist and CD40-targeting vaccination induces immune responses and reduces HIV-1 reservoirs
The Journal of Clinical Investigation 128:4387–4396.

https://doi.org/10.1172/JCI99005
- PubMed
- Google Scholar
1. Chowdhury MZI
2. Turin TC
(2020) Variable selection strategies and its importance in clinical prediction modelling
Family Medicine and Community Health 8:e000262.

https://doi.org/10.1136/fmch-2019-000262
- PubMed
- Google Scholar
1. Chung AW
2. Kumar MP
3. Arnold KB
4. Yu WH
5. Schoen MK
6. Dunphy LJ
7. Suscovich TJ
8. Frahm N
9. Linde C
10. Mahan AE
11. Hoffner M
12. Streeck H
13. Ackerman ME
14. McElrath MJ
15. Schuitemaker H
16. Pau MG
17. Baden LR
18. Kim JH
19. Michael NL
20. Barouch DH
21. Lauffenburger DA
22. Alter G
(2015) Dissecting Polyclonal Vaccine-Induced Humoral Immunity against HIV Using Systems Serology
Cell 163:988–998.

https://doi.org/10.1016/j.cell.2015.10.027
- PubMed
- Google Scholar
(2021) Concerns about SARS-CoV-2 evolution should not hold back efforts to expand vaccination
Nature Reviews. Immunology 21:330–335.

https://doi.org/10.1038/s41577-021-00544-9
- PubMed
- Google Scholar
1. Corbett KS
2. Flynn B
3. Foulds KE
4. Francica JR
5. Boyoglu-Barnum S
6. Werner AP
7. Flach B
8. O’Connell S
9. Bock KW
10. Minai M
11. Nagata BM
12. Andersen H
13. Martinez DR
14. Noe AT
15. Douek N
16. Donaldson MM
17. Nji NN
18. Alvarado GS
19. Edwards DK
20. Flebbe DR
21. Lamb E
22. Doria-Rose NA
23. Lin BC
24. Louder MK
25. O’Dell S
26. Schmidt SD
27. Phung E
28. Chang LA
29. Yap C
30. Todd JPM
31. Pessaint L
32. Ry AV
33. Browne S
34. Greenhouse J
35. Putman-Taylor T
36. Strasbaugh A
37. Campbell TA
38. Cook A
39. Dodson A
40. Steingrebe K
41. Shi W
42. Zhang Y
43. Abiona OM
44. Wang L
45. Pegu A
46. Yang ES
47. Leung K
48. Zhou T
49. Teng IT
50. Widge A
51. Gordon I
52. Novik L
53. Gillespie RA
54. Loomis RJ
(2020) Evaluation of the mRNA-1273 Vaccine against SARS-CoV-2 in Nonhuman Primates
The New England Journal of Medicine 383:1544–1555.

https://doi.org/10.1056/NEJMoa2024671
- Google Scholar
1. Dafni UG
2. Tsiatis AA
(1998) Evaluating surrogate markers of clinical outcome when measured with error
Biometrics 54:1445–1462.

https://doi.org/10.2307/2533670
- PubMed
- Google Scholar
(2020) Approaches and Challenges in SARS-CoV-2 Vaccine Development
Cell Host & Microbe 28:364–370.

https://doi.org/10.1016/j.chom.2020.08.002
- PubMed
- Google Scholar
1. Dufloo J
2. Grzelak L
3. Staropoli I
4. Madec Y
5. Tondeur L
6. Anna F
7. Pelleau S
8. Wiedemann A
9. Planchais C
10. Buchrieser J
11. Robinot R
12. Ungeheuer MN
13. Mouquet H
14. Charneau P
15. White M
16. Lévy Y
17. Hoen B
18. Fontanet A
19. Schwartz O
20. Bruel T
(2021) Asymptomatic and symptomatic SARS-CoV-2 infections elicit polyfunctional antibodies
Cell Reports. Medicine 2:100275.

https://doi.org/10.1016/j.xcrm.2021.100275
- PubMed
- Google Scholar
(2021) Evidence for antibody as a protective correlate for COVID-19 vaccines
Vaccine 39:4423–4428.

https://doi.org/10.1016/j.vaccine.2021.05.063
- PubMed
- Google Scholar
1. Emary KRW
2. Golubchik T
3. Aley PK
4. Ariani CV
5. Angus B
6. Bibi S
7. Blane B
8. Bonsall D
9. Cicconi P
10. Charlton S
11. Clutterbuck EA
12. Collins AM
13. Cox T
14. Darton TC
15. Dold C
16. Douglas AD
17. Duncan CJA
18. Ewer KJ
19. Flaxman AL
20. Faust SN
21. Ferreira DM
22. Feng S
23. Finn A
24. Folegatti PM
25. Fuskova M
26. Galiza E
27. Goodman AL
28. Green CM
29. Green CA
30. Greenland M
31. Hallis B
32. Heath PT
33. Hay J
34. Hill HC
35. Jenkin D
36. Kerridge S
37. Lazarus R
38. Libri V
39. Lillie PJ
40. Ludden C
41. Marchevsky NG
42. Minassian AM
43. McGregor AC
44. Mujadidi YF
45. Phillips DJ
46. Plested E
47. Pollock KM
48. Robinson H
49. Smith A
50. Song R
51. Snape MD
52. Sutherland RK
53. Thomson EC
54. Toshner M
55. Turner DPJ
56. Vekemans J
57. Villafana TL
58. Williams CJ
59. Hill AVS
60. Lambe T
61. Gilbert SC
62. Voysey M
63. Ramasamy MN
64. Pollard AJ
65. COVID-19 Genomics UK consortium
66. AMPHEUS Project
67. Oxford COVID-19 Vaccine Trial Group
(2021) Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B.1.1.7): an exploratory analysis of a randomised controlled trial
Lancet 397:1351–1362.

https://doi.org/10.1016/S0140-6736(21)00628-0
- PubMed
- Google Scholar
1. Eyal N
2. Lipsitch M
(2021) How to test SARS-CoV-2 vaccines ethically even after one is available
Clinical Infectious Diseases: An Official Publication of the Infectious Diseases Society of America 73:ciab182.

https://doi.org/10.1093/cid/ciab182
- Google Scholar
1. Feng S
2. Phillips DJ
3. White T
4. Sayal H
5. Aley PK
6. Bibi S
7. Dold C
8. Fuskova M
9. Gilbert SC
10. Hirsch I
11. Humphries HE
12. Jepson B
13. Kelly EJ
14. Plested E
15. Shoemaker K
16. Thomas KM
17. Vekemans J
18. Villafana TL
19. Lambe T
20. Pollard AJ
21. Voysey M
22. Oxford COVID Vaccine Trial Group
(2021) Correlates of protection against symptomatic and asymptomatic SARS-CoV-2 infection
Nature Medicine 27:2032–2040.

https://doi.org/10.1038/s41591-021-01540-1
- PubMed
- Google Scholar
1. Fenwick C
2. Croxatto A
3. Coste AT
4. Pojer F
5. André C
6. Pellaton C
7. Farina A
8. Campos J
9. Hacker D
10. Lau K
11. Bosch B-J
12. Gonseth Nussle S
13. Bochud M
14. D’Acremont V
15. Trono D
16. Greub G
17. Pantaleo G
(2021) Changes in SARS-CoV-2 Spike versus Nucleoprotein Antibody Responses Impact the Estimates of Infections in Population-Based Seroprevalence Studies
Journal of Virology 95:e01828-20.

https://doi.org/10.1128/JVI.01828-20
- PubMed
- Google Scholar
1. Flamar AL
2. Zurawski S
3. Scholz F
4. Gayet I
5. Ni L
6. Li XH
7. Klechevsky E
8. Quinn J
9. Oh S
10. Kaplan DH
11. Banchereau J
12. Zurawski G
(2012) Noncovalent assembly of anti-dendritic cell antibodies and antigens for evoking immune responses in vitro and in vivo
Journal of Immunology 189:2645–2655.

https://doi.org/10.4049/jimmunol.1102390
- PubMed
- Google Scholar
1. Gaebler C
2. Wang Z
3. Lorenzi JCC
4. Muecksch F
5. Finkin S
6. Tokuyama M
7. Cho A
8. Jankovic M
9. Schaefer-Babajew D
10. Oliveira TY
11. Cipolla M
12. Viant C
13. Barnes CO
14. Bram Y
15. Breton G
16. Hägglöf T
17. Mendoza P
18. Hurley A
19. Turroja M
20. Gordon K
21. Millard KG
22. Ramos V
23. Schmidt F
24. Weisblum Y
25. Jha D
26. Tankelevich M
27. Martinez-Delgado G
28. Yee J
29. Patel R
30. Dizon J
31. Unson-O’Brien C
32. Shimeliovich I
33. Robbiani DF
34. Zhao Z
35. Gazumyan A
36. Schwartz RE
37. Hatziioannou T
38. Bjorkman PJ
39. Mehandru S
40. Bieniasz PD
41. Caskey M
42. Nussenzweig MC
(2021) Evolution of antibody immunity to SARS-CoV-2
Nature 591:639–644.

https://doi.org/10.1038/s41586-021-03207-w
- PubMed
- Google Scholar
1. Godot V
2. Tcherakian C
3. Gil L
4. Cervera-Marzal I
5. Li G
6. Cheng L
7. Ortonne N
8. Lelièvre JD
9. Pantaleo G
10. Fenwick C
11. Centlivre M
12. Mouquet H
13. Cardinaud S
14. Zurawski SM
15. Zurawski G
16. Milpied P
17. Su L
18. Lévy Y
(2020) TLR-9 agonist and CD40-targeting vaccination induces HIV-1 envelope-specific B cells with a diversified immunoglobulin repertoire in humanized mice
PLOS Pathogens 16:e1009025.

https://doi.org/10.1371/journal.ppat.1009025
- PubMed
- Google Scholar
1. Gonçalves A
2. Bertrand J
3. Ke R
4. Comets E
5. de Lamballerie X
6. Malvy D
7. Pizzorno A
8. Terrier O
9. Rosa Calatrava M
10. Mentré F
11. Smith P
12. Perelson AS
13. Guedj J
(2020) Timing of Antiviral Treatment Initiation is Critical to Reduce SARS-CoV-2 Viral Load
CPT 9:509–514.

https://doi.org/10.1002/psp4.12543
- PubMed
- Google Scholar
1. Gonçalves A
2. Maisonnasse P
3. Donati F
4. Albert M
5. Behillil S
6. Contreras V
7. Naninck T
8. Marlin R
9. Solas C
10. Pizzorno A
11. Lemaitre J
12. Kahlaoui N
13. Terrier O
14. Ho Tsong Fang R
15. Enouf V
16. Dereuddre-Bosquet N
17. Brisebarre A
18. Touret F
19. Chapon C
20. Hoen B
21. Lina B
22. Rosa Calatrava M
23. de Lamballerie X
24. Mentré F
25. Le Grand R
26. van der Werf S
27. Guedj J
(2021) SARS-CoV-2 viral dynamics in non-human primates
PLOS Computational Biology 17:e1008785.

https://doi.org/10.1371/journal.pcbi.1008785
- PubMed
- Google Scholar
(2017) The Role of Infected Cell Proliferation in the Clearance of Acute HBV Infection in Humans
Viruses 9:E350.

https://doi.org/10.3390/v9110350
- PubMed
- Google Scholar
(2019) Within-host mathematical models of hepatitis B virus infection: Past, present, and future
Current Opinion in Systems Biology 18:27–35.

https://doi.org/10.1016/j.coisb.2019.10.003
- PubMed
- Google Scholar
Preprint
(2020) Mathematical Modeling Explains Differential SARS CoV-2 Kinetics in Lung and Nasal Passages in Remdesivir Treated Rhesus Macaques
bioRxiv.

https://doi.org/10.1101/2020.06.21.163550
- Google Scholar
(2022) Multi-scale modelling reveals that early super-spreader events are a likely contributor to novel variant predominance
Journal of the Royal Society, Interface 19:20210811.

https://doi.org/10.1098/rsif.2021.0811
- PubMed
- Google Scholar
(2021) Thrombotic Thrombocytopenia after ChAdOx1 nCov-19 Vaccination
The New England Journal of Medicine 384:2092–2101.

https://doi.org/10.1056/NEJMoa2104840
- PubMed
- Google Scholar
1. Gutmann C
2. Takov K
3. Burnap SA
4. Singh B
5. Ali H
6. Theofilatos K
7. Reed E
8. Hasman M
9. Nabeebaccus A
10. Fish M
11. McPhail MJ
12. O’Gallagher K
13. Schmidt LE
14. Cassel C
15. Rienks M
16. Yin X
17. Auzinger G
18. Napoli S
19. Mujib SF
20. Trovato F
21. Sanderson B
22. Merrick B
23. Niazi U
24. Saqi M
25. Dimitrakopoulou K
26. Fernández-Leiro R
27. Braun S
28. Kronstein-Wiedemann R
29. Doores KJ
30. Edgeworth JD
31. Shah AM
32. Bornstein SR
33. Tonn T
34. Hayday AC
35. Giacca M
36. Shankar-Hari M
37. Mayr M
(2021) SARS-CoV-2 RNAemia and proteomic trajectories inform prognostication in COVID-19 patients admitted to intensive care
Nature Communications 12:3406.

https://doi.org/10.1038/s41467-021-23494-1
- PubMed
- Google Scholar
1. Jin P
2. Li J
3. Pan H
4. Wu Y
5. Zhu F
(2021) Immunological surrogate endpoints of COVID-2019 vaccines: the evidence we have versus the evidence we need
Signal Transduction and Targeted Therapy 6:48.

https://doi.org/10.1038/s41392-021-00481-y
- PubMed
- Google Scholar
1. Johnson M
2. Wagstaffe HR
3. Gilmour KC
4. Mai AL
5. Lewis J
6. Hunt A
7. Sirr J
8. Bengt C
9. Grandjean L
10. Goldblatt D
(2020) Evaluation of a novel multiplexed assay for determining IgG levels and functional activity to SARS-CoV-2
Journal of Clinical Virology 130:104572.

https://doi.org/10.1016/j.jcv.2020.104572
- PubMed
- Google Scholar
Preprint
(2020) Kinetics of SARS-CoV-2 Infection in the Human Upper and Lower Respiratory Tracts and Their Relationship with Infectiousness
medRxiv.

https://doi.org/10.1101/2020.09.25.20201772
- Google Scholar
1. Ke R
2. Zitzmann C
3. Ho DD
4. Ribeiro RM
5. Perelson AS
(2021) In vivo kinetics of SARS-CoV-2 infection and its relationship with a person’s infectiousness
PNAS 118:e2111477118.

https://doi.org/10.1073/pnas.2111477118
- PubMed
- Google Scholar
1. Kester KE
2. Cummings JF
3. Ofori-Anyinam O
4. Ockenhouse CF
5. Krzych U
6. Moris P
7. Schwenk R
8. Nielsen RA
9. Debebe Z
10. Pinelis E
11. Juompan L
12. Williams J
13. Dowler M
14. Stewart VA
15. Wirtz RA
16. Dubois M-C
17. Lievens M
18. Cohen J
19. Ballou WR
20. Heppner DG Jr
21. RTS,S Vaccine Evaluation Group
(2009) Randomized, double-blind, phase 2a trial of falciparum malaria vaccines RTS,S/AS01B and RTS,S/AS02A in malaria-naive adults: safety, efficacy, and immunologic associates of protection
The Journal of Infectious Diseases 200:337–346.

https://doi.org/10.1086/600120
- PubMed
- Google Scholar
1. Khoury DS
2. Cromer D
3. Reynaldi A
4. Schlub TE
5. Wheatley AK
6. Juno JA
7. Subbarao K
8. Kent SJ
9. Triccas JA
10. Davenport MP
(2021) Neutralizing antibody levels are highly predictive of immune protection from symptomatic SARS-CoV-2 infection
Nature Medicine 27:1205–1211.

https://doi.org/10.1038/s41591-021-01377-8
- PubMed
- Google Scholar
1. Kim KS
2. Ejima K
3. Iwanami S
4. Fujita Y
5. Ohashi H
6. Koizumi Y
7. Asai Y
8. Nakaoka S
9. Watashi K
10. Aihara K
11. Thompson RN
12. Ke R
13. Perelson AS
14. Iwami S
(2021) A quantitative model used to compare within-host SARS-CoV-2, MERS-CoV, and SARS-CoV dynamics provides insights into the pathogenesis and treatment of SARS-CoV-2
PLOS Biology 19:e3001128.

https://doi.org/10.1371/journal.pbio.3001128
- PubMed
- Google Scholar
(2021) Correlates of Vaccine-Induced Protection against SARS-CoV-2
Vaccines 9:238.

https://doi.org/10.3390/vaccines9030238
- PubMed
- Google Scholar
(2021) SARS-CoV-2 spike variants exhibit differential infectivity and neutralization resistance to convalescent or post-vaccination sera
Cell Host & Microbe 29:522–528.

https://doi.org/10.1016/j.chom.2021.03.008
- PubMed
- Google Scholar
1. Leung NHL
2. Chu DKW
3. Shiu EYC
4. Chan KH
5. McDevitt JJ
6. Hau BJP
7. Yen HL
8. Li Y
9. Ip DKM
10. Peiris JSM
11. Seto WH
12. Leung GM
13. Milton DK
14. Cowling BJ
(2020) Respiratory virus shedding in exhaled breath and efficacy of face masks
Nature Medicine 26:676–680.

https://doi.org/10.1038/s41591-020-0843-2
- PubMed
- Google Scholar
1. Lui G
2. Ling L
3. Lai CK
4. Tso EY
5. Fung KS
6. Chan V
7. Ho TH
8. Luk F
9. Chen Z
10. Ng JK
11. Chow K-M
12. Cheng PK
13. Chan RC
14. Tsang DN
15. Gomersall CD
16. Hui DS
17. Chan PK
(2020) Viral dynamics of SARS-CoV-2 across a spectrum of disease severity in COVID-19
The Journal of Infection 81:318–356.

https://doi.org/10.1016/j.jinf.2020.04.014
- PubMed
- Google Scholar
1. Lustig Y
2. Nemet I
3. Kliker L
4. Zuckerman N
5. Yishai R
6. Alroy-Preis S
7. Mendelson E
8. Mandelboim M
(2021) Neutralizing Response against Variants after SARS-CoV-2 Infection and One Dose of BNT162b2
The New England Journal of Medicine 384:2453–2454.

https://doi.org/10.1056/NEJMc2104036
- PubMed
- Google Scholar
1. Marc A
2. Kerioui M
3. Blanquart F
4. Bertrand J
5. Mitjà O
6. Corbacho-Monné M
7. Marks M
8. Guedj J
(2021) Quantifying the relationship between SARS-CoV-2 viral load and infectiousness
eLife 10:e69302.

https://doi.org/10.7554/eLife.69302
- PubMed
- Google Scholar
1. Marks M
2. Millat-Martinez P
3. Ouchi D
4. Roberts CH
5. Alemany A
6. Corbacho-Monné M
7. Ubals M
8. Tobias A
9. Tebé C
10. Ballana E
11. Bassat Q
12. Baro B
13. Vall-Mayans M
14. G-Beiras C
15. Prat N
16. Ara J
17. Clotet B
18. Mitjà O
(2021) Transmission of COVID-19 in 282 clusters in Catalonia, Spain: a cohort study
The Lancet. Infectious Diseases 21:629–636.

https://doi.org/10.1016/S1473-3099(20)30985-3
- PubMed
- Google Scholar
1. Marlin R
2. Godot V
3. Cardinaud S
4. Galhaut M
5. Coleon S
6. Zurawski S
7. Dereuddre-Bosquet N
8. Cavarelli M
9. Gallouët AS
10. Maisonnasse P
11. Dupaty L
12. Fenwick C
13. Naninck T
14. Lemaitre J
15. Gomez-Pacheco M
16. Kahlaoui N
17. Contreras V
18. Relouzat F
19. Fang RHT
20. Wang Z
21. Ellis J
22. Chapon C
23. Centlivre M
24. Wiedemann A
25. Lacabaratz C
26. Surenaud M
27. Szurgot I
28. Liljeström P
29. Planas D
30. Bruel T
31. Schwartz O
32. van der Werf S
33. Pantaleo G
34. Prague M
35. Thiébaut R
36. Zurawski G
37. Lévy Y
38. Grand RL
(2021) Targeting SARS-CoV-2 receptor-binding domain to cells expressing CD40 improves protection to infection in convalescent macaques
Nature Communications 12:5215.

https://doi.org/10.1038/s41467-021-25382-0
- PubMed
- Google Scholar
1. McMahan K
2. Yu J
3. Mercado NB
4. Loos C
5. Tostanoski LH
6. Chandrashekar A
7. Liu J
8. Peter L
9. Atyeo C
10. Zhu A
11. Bondzie EA
12. Dagotto G
13. Gebre MS
14. Jacob-Dolan C
15. Li Z
16. Nampanya F
17. Patel S
18. Pessaint L
19. Van Ry A
20. Blade K
21. Yalley-Ogunro J
22. Cabus M
23. Brown R
24. Cook A
25. Teow E
26. Andersen H
27. Lewis MG
28. Lauffenburger DA
29. Alter G
30. Barouch DH
(2021) Correlates of protection against SARS-CoV-2 in rhesus macaques
Nature 590:630–634.

https://doi.org/10.1038/s41586-020-03041-6
- PubMed
- Google Scholar
1. Mercado NB
2. Zahn R
3. Wegmann F
4. Loos C
5. Chandrashekar A
6. Yu J
7. Liu J
8. Peter L
9. McMahan K
10. Tostanoski LH
11. He X
12. Martinez DR
13. Rutten L
14. Bos R
15. Vellinga J
16. Custers J
17. Langedijk JP
18. Kwaks T
19. Bakkers MJG
20. Zuijdgeest D
21. Huber SKR
22. Atyeo C
23. Fischinger S
24. Burke JS
25. Feldman J
26. Hauser BM
27. Caradonna TM
28. Bondzie EA
29. Dagotto G
30. Gebre MS
31. Hoffman E
32. Jacob-Dolan C
33. Kirilova M
34. Li Z
35. Lin Z
36. Mahrokhian SH
37. Maxfield LF
38. Nampanya F
39. Nityanandam R
40. Nkolola JP
41. Patel S
42. Ventura JD
43. Verrington K
44. Wan H
45. Pessaint L
46. Ry AV
(2020) Single-shot Ad26 vaccine protects against SARS-CoV-2 in rhesus macaques
Nature 586:583–588.

https://doi.org/10.1038/s41586-020-2607-z
- Google Scholar
1. Miao H
2. Xia X
3. Perelson AS
4. Wu H
(2011) ON IDENTIFIABILITY OF NONLINEAR ODE MODELS AND APPLICATIONS IN VIRAL DYNAMICS
SIAM Review. Society for Industrial and Applied Mathematics 53:3–39.

https://doi.org/10.1137/090757009
- PubMed
- Google Scholar
1. Muñoz-Fontela C
(2020) Animal models for COVID-19
Nature 586:509–515.

https://doi.org/10.1038/s41586-020-2787-6
- Google Scholar
1. Myers MA
2. Smith AP
3. Lane LC
4. Moquin DJ
5. Aogo R
6. Woolard S
7. Thomas P
8. Vogel P
9. Smith AM
(2021) Dynamically linking influenza virus infection kinetics, lung injury, inflammation, and disease severity
eLife 10:e68864.

https://doi.org/10.7554/eLife.68864
- PubMed
- Google Scholar
1. Néant N
(2021) Groups, Modeling SARS-CoV-2 viral kinetics and association with mortality in hospitalized patients from the French COVID cohort
PNAS 118:.

https://doi.org/10.1073/pnas.2017962118
- Google Scholar
1. Nyberg T
2. Ferguson NM
3. Nash SG
4. Webster HH
5. Flaxman S
6. Andrews N
7. Hinsley W
8. Bernal JL
9. Kall M
10. Bhatt S
11. Blomquist P
12. Zaidi A
13. Volz E
14. Aziz NA
15. Harman K
16. Funk S
17. Abbott S
18. Hope R
19. Charlett A
20. Chand M
21. Ghani AC
22. Seaman SR
23. Dabrera G
24. De Angelis D
25. Presanis AM
26. Thelwall S
27. COVID-19 Genomics UK (COG-UK) consortium
(2022) Comparative analysis of the risks of hospitalisation and death associated with SARS-CoV-2 omicron (B.1.1.529) and delta (B.1.617.2) variants in England: a cohort study
Lancet 399:1303–1312.

https://doi.org/10.1016/S0140-6736(22)00462-7
- PubMed
- Google Scholar
1. Pinky L
2. Burke CW
3. Russell CJ
4. Smith AM
(2021) Quantifying dose-, strain-, and tissue-specific kinetics of parainfluenza virus infection
PLOS Computational Biology 17:e1009299.

https://doi.org/10.1371/journal.pcbi.1009299
- PubMed
- Google Scholar
1. Pizzorno A
2. Padey B
3. Julien T
4. Trouillet-Assant S
5. Traversier A
6. Errazuriz-Cerda E
7. Fouret J
8. Dubois J
9. Gaymard A
10. Lescure FX
11. Dulière V
12. Brun P
13. Constant S
14. Poissy J
15. Lina B
16. Yazdanpanah Y
17. Terrier O
18. Rosa-Calatrava M
(2020) Characterization and Treatment of SARS-CoV-2 in Nasal and Bronchial Human Airway Epithelia
Cell Reports. Medicine 1:100059.

https://doi.org/10.1016/j.xcrm.2020.100059
- PubMed
- Google Scholar
1. Planas D
2. Bruel T
3. Grzelak L
4. Guivel-Benhassine F
5. Staropoli I
6. Porrot F
7. Planchais C
8. Buchrieser J
9. Rajah MM
10. Bishop E
11. Albert M
12. Donati F
13. Prot M
14. Behillil S
15. Enouf V
16. Maquart M
17. Smati-Lafarge M
18. Varon E
19. Schortgen F
20. Yahyaoui L
21. Gonzalez M
22. De Sèze J
23. Péré H
24. Veyer D
25. Sève A
26. Simon-Lorière E
27. Fafi-Kremer S
28. Stefic K
29. Mouquet H
30. Hocqueloux L
31. van der Werf S
32. Prazuck T
33. Schwartz O
(2021) Sensitivity of infectious SARS-CoV-2 B.1.1.7 and B.1.351 variants to neutralizing antibodies
Nature Medicine 27:917–924.

https://doi.org/10.1038/s41591-021-01318-5
- PubMed
- Google Scholar
1. Plotkin SA
2. Gilbert PB
(2012) Nomenclature for immune correlates of protection after vaccination
Clinical Infectious Diseases 54:1615–1617.

https://doi.org/10.1093/cid/cis238
- PubMed
- Google Scholar
1. Plotkin SA
(2013) Complex correlates of protection after vaccination
Clinical Infectious Diseases 56:1458–1465.

https://doi.org/10.1093/cid/cit048
- PubMed
- Google Scholar
1. Plotkin SA
(2020) Updates on immunologic correlates of vaccine-induced protection
Vaccine 38:2250–2257.

https://doi.org/10.1016/j.vaccine.2019.10.046
- PubMed
- Google Scholar
1. Pollard AJ
2. Bijker EM
(2021) A guide to vaccinology: from basic principles to new developments
Nature Reviews. Immunology 21:83–100.

https://doi.org/10.1038/s41577-020-00479-7
- PubMed
- Google Scholar
1. Rhodes SJ
2. Guedj J
3. Fletcher HA
4. Lindenstrøm T
5. Scriba TJ
6. Evans TG
7. Knight GM
8. White RG
(2018) Using vaccine Immunostimulation/Immunodynamic modelling methods to inform vaccine dose decision-making
NPJ Vaccines 3:36.

https://doi.org/10.1038/s41541-018-0075-3
- PubMed
- Google Scholar
1. Robinot R
2. Hubert M
3. de Melo GD
4. Lazarini F
5. Bruel T
6. Smith N
7. Levallois S
8. Larrous F
9. Fernandes J
10. Gellenoncourt S
11. Rigaud S
12. Gorgette O
13. Thouvenot C
14. Trébeau C
15. Mallet A
16. Duménil G
17. Gobaa S
18. Etournay R
19. Lledo P-M
20. Lecuit M
21. Bourhy H
22. Duffy D
23. Michel V
24. Schwartz O
25. Chakrabarti LA
(2021) SARS-CoV-2 infection induces the dedifferentiation of multiciliated cells and impairs mucociliary clearance
Nature Communications 12:4354.

https://doi.org/10.1038/s41467-021-24521-x
- PubMed
- Google Scholar
(2007) A contemporary view of coronavirus transcription
Journal of Virology 81:20–29.

https://doi.org/10.1128/JVI.01358-06
- PubMed
- Google Scholar
(2018) Influenza Virus Infection Model With Density Dependence Supports Biphasic Viral Decay
Frontiers in Microbiology 9:1554.

https://doi.org/10.3389/fmicb.2018.01554
- PubMed
- Google Scholar
1. Tauzin A
2. Nayrac M
3. Benlarbi M
4. Gong SY
5. Gasser R
6. Beaudoin-Bussières G
7. Brassard N
8. Laumaea A
9. Vézina D
10. Prévost J
11. Anand SP
12. Bourassa C
13. Gendron-Lepage G
14. Medjahed H
15. Goyette G
16. Niessl J
17. Tastet O
18. Gokool L
19. Morrisseau C
20. Arlotto P
21. Stamatatos L
22. McGuire AT
23. Larochelle C
24. Uchil P
25. Lu M
26. Mothes W
27. De Serres G
28. Moreira S
29. Roger M
30. Richard J
31. Martel-Laferrière V
32. Duerr R
33. Tremblay C
34. Kaufmann DE
35. Finzi A
(2021) A single dose of the SARS-CoV-2 vaccine BNT162b2 elicits Fc-mediated antibody effector functions and T cell responses
Cell Host & Microbe 29:1137–1150.

https://doi.org/10.1016/j.chom.2021.06.001
- PubMed
- Google Scholar
(2014) Evaluation of bootstrap methods for estimating uncertainty of parameters in nonlinear mixed-effects models: a simulation study in population pharmacokinetics
Journal of Pharmacokinetics and Pharmacodynamics 41:15–33.

https://doi.org/10.1007/s10928-013-9343-z
- PubMed
- Google Scholar
1. Travaglini KJ
2. Nabhan AN
3. Penland L
4. Sinha R
5. Gillich A
6. Sit RV
7. Chang S
8. Conley SD
9. Mori Y
10. Seita J
11. Berry GJ
12. Shrager JB
13. Metzger RJ
14. Kuo CS
15. Neff N
16. Weissman IL
17. Quake SR
18. Krasnow MA
(2020) A molecular cell atlas of the human lung from single-cell RNA sequencing
Nature 587:619–625.

https://doi.org/10.1038/s41586-020-2922-4
- PubMed
- Google Scholar
1. Vanshylla K
2. Di Cristanziano V
3. Kleipass F
4. Dewald F
5. Schommers P
6. Gieselmann L
7. Gruell H
8. Schlotz M
9. Ercanoglu MS
10. Stumpf R
11. Mayer P
12. Zehner M
13. Heger E
14. Johannis W
15. Horn C
16. Suárez I
17. Jung N
18. Salomon S
19. Eberhardt KA
20. Gathof B
21. Fätkenheuer G
22. Pfeifer N
23. Eggeling R
24. Augustin M
25. Lehmann C
26. Klein F
(2021) Kinetics and correlates of the neutralizing antibody response to SARS-CoV-2 infection in humans
Cell Host & Microbe 29:917–929.

https://doi.org/10.1016/j.chom.2021.04.015
- PubMed
- Google Scholar
1. Wang S
2. Pan Y
3. Wang Q
4. Miao H
5. Brown AN
6. Rong L
(2020) Modeling the viral dynamics of SARS-CoV-2 infection
Mathematical Biosciences 328:108438.

https://doi.org/10.1016/j.mbs.2020.108438
- PubMed
- Google Scholar
Software
1. World Health Organization
(2021) World Health Organisation, COVID-19 vaccine tracker and landscape
COVID-19.

https://www.who.int/publications/m/item/draft-landscape-of-covid-19-candidate-vaccines)
Book
1. Wu L
(2009) Mixed Effects Models for Complex Data
Chapman and Hall/CRC.

https://doi.org/10.1201/9781420074086
- Google Scholar
1. Yu J
2. Tostanoski LH
3. Peter L
4. Mercado NB
5. McMahan K
6. Mahrokhian SH
7. Nkolola JP
8. Liu J
9. Li Z
10. Chandrashekar A
11. Martinez DR
12. Loos C
13. Atyeo C
14. Fischinger S
15. Burke JS
16. Slein MD
17. Chen Y
18. Zuiani A
19. Lelis FJN
20. Travers M
21. Habibi S
22. Pessaint L
23. Van Ry A
24. Blade K
25. Brown R
26. Cook A
27. Finneyfrock B
28. Dodson A
29. Teow E
30. Velasco J
31. Zahn R
32. Wegmann F
33. Bondzie EA
34. Dagotto G
35. Gebre MS
36. He X
37. Jacob-Dolan C
38. Kirilova M
39. Kordana N
40. Lin Z
41. Maxfield LF
42. Nampanya F
43. Nityanandam R
44. Ventura JD
45. Wan H
46. Cai Y
47. Chen B
48. Schmidt AG
49. Wesemann DR
50. Baric RS
51. Alter G
52. Andersen H
53. Lewis MG
54. Barouch DH
(2020) DNA vaccine protection against SARS-CoV-2 in rhesus macaques
Science 369:806–811.

https://doi.org/10.1126/science.abc6284
- PubMed
- Google Scholar
1. Zhou D
2. Dejnirattisai W
3. Supasa P
4. Liu C
5. Mentzer AJ
6. Ginn HM
7. Zhao Y
8. Duyvesteyn HME
9. Tuekprakhon A
10. Nutalai R
11. Wang B
12. Paesen GC
13. Lopez-Camacho C
14. Slon-Campos J
15. Hallis B
16. Coombes N
17. Bewley K
18. Charlton S
19. Walter TS
20. Skelly D
21. Lumley SF
22. Dold C
23. Levin R
24. Dong T
25. Pollard AJ
26. Knight JC
27. Crook D
28. Lambe T
29. Clutterbuck E
30. Bibi S
31. Flaxman A
32. Bittaye M
33. Belij-Rammerstorfer S
34. Gilbert S
35. James W
36. Carroll MW
37. Klenerman P
38. Barnes E
39. Dunachie SJ
40. Fry EE
41. Mongkolsapaya J
42. Ren J
43. Stuart DI
44. Screaton GR
(2021) Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera
Cell 184:2348–2361.

https://doi.org/10.1016/j.cell.2021.02.037
- PubMed
- Google Scholar
1. Zost SJ
2. Gilchuk P
3. Chen RE
4. Case JB
5. Reidy JX
6. Trivette A
7. Nargi RS
8. Sutton RE
9. Suryadevara N
10. Chen EC
11. Binshtein E
12. Shrihari S
13. Ostrowski M
14. Chu HY
15. Didier JE
16. MacRenaris KW
17. Jones T
18. Day S
19. Myers L
20. Eun-Hyung Lee F
21. Nguyen DC
22. Sanz I
23. Martinez DR
24. Rothlauf PW
25. Bloyet L-M
26. Whelan SPJ
27. Baric RS
28. Thackray LB
29. Diamond MS
30. Carnahan RH
31. Crowe JE Jr
(2020) Rapid isolation and profiling of a diverse panel of human monoclonal antibodies targeting the SARS-CoV-2 spike protein
Nature Medicine 26:1422–1427.

https://doi.org/10.1038/s41591-020-0998-x
- PubMed
- Google Scholar
1. Zurawski G
2. Shen X
3. Zurawski S
4. Tomaras GD
5. Montefiori DC
6. Roederer M
7. Ferrari G
8. Lacabaratz C
9. Klucar P
10. Wang Z
11. Foulds KE
12. Kao SF
13. Yu X
14. Sato A
15. Yates NL
16. LaBranche C
17. Stanfield-Oakley S
18. Kibler K
19. Jacobs B
20. Salazar A
21. Self S
22. Fulp W
23. Gottardo R
24. Galmin L
25. Weiss D
26. Cristillo A
27. Pantaleo G
28. Levy Y
(2017) Superiority in Rhesus Macaques of Targeting HIV-1 Env gp140 to CD40 versus LOX-1 in Combination with Replication-Competent NYVAC-KC for Induction of Env-Specific Antibody and T Cell Responses
Journal of Virology 91:e01596-16.

https://doi.org/10.1128/JVI.01596-16
- PubMed
- Google Scholar

Article and author information

Author details

Marie Alexandre

University of Bordeaux, Department of Public Health, Inserm Bordeaux Population Health Research Centre, Inria SISTM, Bordeaux, France

Contribution
Conceptualization, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3557-7075
Romain Marlin

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Conceptualization, Investigation, Resources, Visualization, Writing – original draft, Writing – review and editing

Contributed equally with
Mélanie Prague

Competing interests
No competing interests declared
Mélanie Prague

University of Bordeaux, Department of Public Health, Inserm Bordeaux Population Health Research Centre, Inria SISTM, Bordeaux, France

Contribution
Methodology, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing

Contributed equally with
Romain Marlin

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9809-7848
Severin Coleon
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Nidhal Kahlaoui

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Visualization, Writing – review and editing

Competing interests
No competing interests declared
Sylvain Cardinaud
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Thibaut Naninck

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Visualization, Writing – review and editing

Competing interests
No competing interests declared
Benoit Delache

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Mathieu Surenaud
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Mathilde Galhaut

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Nathalie Dereuddre-Bosquet

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Mariangela Cavarelli

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Pauline Maisonnasse

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-0555-207X
Mireille Centlivre
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Christine Lacabaratz
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Aurelie Wiedemann
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Sandra Zurawski

Baylor Scott and White Research Institute, Dallas, United States

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Gerard Zurawski

Baylor Scott and White Research Institute, Dallas, United States

Contribution
Investigation, Resources, Writing – review and editing

Competing interests
No competing interests declared
Olivier Schwartz
1. Vaccine Research Institute, Créteil, France
2. Virus & Immunity Unit, Department of Virology, Institut Pasteur, Paris, France
3. CNRS UMR 3569, Paris, France
Contribution
Resources, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-0729-1475
Rogier W Sanders

Department of Medical Microbiology, Amsterdam UMC, University of Amsterdam Amsterdam Infection & Immunity Institute, Amsterdam, Netherlands

Contribution
Resources, Writing – review and editing

Competing interests
No competing interests declared
Roger Le Grand

Center for Immunology of Viral, Auto-immune, Hematological and Bacterial Diseases (IMVA-HB/IDMIT), Université Paris-Saclay, Inserm, CEA, Fontenay-aux-Roses, France

Contribution
Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4928-4484
Yves Levy
1. Vaccine Research Institute, Créteil, France
2. Inserm U955, Créteil, France
3. AP-HP, Hôpital Henri-Mondor Albert-Chenevier, Service d'Immunologie Clinique et Maladies Infectieuses, Créteil, France
Contribution
Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared
Rodolphe Thiébaut
1. University of Bordeaux, Department of Public Health, Inserm Bordeaux Population Health Research Centre, Inria SISTM, Bordeaux, France
2. Vaccine Research Institute, Créteil, France
3. CHU Bordeaux, Department of Medical information, Bordeaux, France
Contribution
Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review and editing

For correspondence
rodolphe.thiebaut@u-bordeaux.fr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5235-3962

Funding

Agence Nationale de la Recherche (ANR-10-LABX-77-01)

Yves Levy
Rodolphe Thiébaut

Agence Nationale de la Recherche (ANR-11- 1018 INBS-0008)

Roger Le Grand

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We would like to thank J Guedj and O Terrier for fruitful discussions on the model definition. We thank S Langlois, J Demilly, N Dhooge, P Le Calvez, M Potier, JM Robert, T Prot, and C Dodan for the NHP experiments; L Bossevot, M Leonec, L Moenne-Loccoz, M Calpin-Lebreau, and J Morin for the RT-qPCR, ELISpot and Luminex assays, and for the preparation of reagents; A-S Gallouët, M Gomez-Pacheco, and W Gros for NHP T-cell assays and flow cytometry; B Fert for her help with the CT scans; M Barendji, J Dinh, and E Guyon for the NHP sample processing; S Keyser for the transports organization; F Ducancel and Y Gorin for their help with the logistics and safety management; I Mangeot for her help with resources management and B Targat contributed to data management. The monkey and syringe pictures in Figure 1 was created with BioRender.com. This work was supported by INSERM and the Investissements d’Avenir program, Vaccine Research Institute (VRI), managed by the ANR under reference ANR-10-LABX-77-01. MA has been funded by INRIA PhD grant. The Infectious Disease Models and Innovative Therapies (IDMIT) research infrastructure is supported by the 'Programme Investissements d’Avenir', managed by the ANR under reference ANR-11-INBS-0008. The Fondation Bettencourt Schueller and the Region Ile-de-France contributed to the implementation of IDMIT’s facilities and imaging technologies used to define volume of respiratory tract. The NHP study received financial support from REACTing, the Fondation pour la Recherche Medicale (FRM; AM-CoV-Path). We thank Lixoft SAS for their support. Numerical computations were in part carried out using the PlaFRIM experimental testbed, supported by Inria, CNRS (LABRI and IMB), Université de Bordeaux, Bordeaux INP, and Conseil Régional d’Aquitaine (see https://www.plafrim.fr). We thank Miles Davenport and Frederik Graw as Senior Editor and Reviewing Editor of our paper, respectively, and the three anonymous reviewers for their time and their constructive comments.

Ethics

Cynomolgus macaques (Macaca fascicularis), aged 37–66 months (18 females and 13 males) and originating from Mauritian AAALAC certified breeding centers were used in this study. All animals were housed in IDMIT facilities (CEA, Fontenay-aux-roses), under BSL2 and BSL-3 containment when necessary (Animal facility authorization #D92-032-02, Préfecture des Hauts de Seine, France) and in compliance with European Directive 2010/63/EU, the French regulations and the Standards for Human Care and Use of Laboratory Animals, of the Office for Laboratory Animal Welfare (OLAW, assurance number #A5826-01, US). The protocols were approved by the institutional ethical committee “Comité d’Ethique en Expérimentation Animale du Commissariat à l’Energie Atomique et aux Energies Alternatives” (CEtEA #44) under statement number A20-011. The study was authorized by the “Research, Innovation and Education Ministry” under registration number APAFIS#24434–2020030216532863 v1.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,181

views
291

downloads
9

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Marie Alexandre
Romain Marlin
Mélanie Prague
Severin Coleon
Nidhal Kahlaoui
Sylvain Cardinaud
Thibaut Naninck
Benoit Delache
Mathieu Surenaud
Mathilde Galhaut
Nathalie Dereuddre-Bosquet
Mariangela Cavarelli
Pauline Maisonnasse
Mireille Centlivre
Christine Lacabaratz
Aurelie Wiedemann
Sandra Zurawski
Gerard Zurawski
Olivier Schwartz
Rogier W Sanders
Roger Le Grand
Yves Levy
Rodolphe Thiébaut

(2022)

Modelling the response to vaccine in non-human primates to define SARS-CoV-2 mechanistic correlates of protection

eLife 11:e75427.

https://doi.org/10.7554/eLife.75427

Share this article

Cite this article

Design of the study 1 and viral dynamics.

Figure 1—source data 1

Figure 1—source data 2

Figure 1—source data 3

Figure 1—source data 4

Figure 1—source data 5

Figure 1—source data 6

Figure 1—source data 7

Figure 1—source data 8

Figure 1—source data 9

Figure 1—source data 10

Figure 1—source data 11

Mechanistic modelling.

Figure 2—source data 1

Figure 2—source data 2

Figure 2—source data 3

Harvest times and measurements.

Immune markers.

Figure 4—source data 1

Figure 4—source data 2

Study design and modeling results for the second study testing two-component spike nanoparticle vaccine.

Figure 5—source data 1

Figure 5—source data 2

Comparison of models evaluating the difference of viral infectivity (β), loss of infected cells (δ), and viral production (P) between the nasopharynx and the trachea.

Comparison of models evaluating the adjustment of the viral infectivity (β), the loss rate of infected cells (δ), the viral production (P), and the viral clearance (c) for the groups of treatment.

Results of the forward selection approach applied on the 25 simulated white-noise random variables.

Model parameters for viral dynamics in both the nasopharynx and the trachea estimated by the model adjusted for groups of intervention.

Model parameters for viral dynamics in both the nasopharynx and the trachea estimated by the model with the viral infectivity adjusted for ACE2-RBD-binding inhibition and the loss rate of infected cells adjusted for the group of treatment.

Author details

Marie Alexandre

Contribution

Competing interests

Romain Marlin

Contribution

Contributed equally with

Competing interests

Mélanie Prague

Contribution

Contributed equally with

Competing interests

Severin Coleon

Contribution

Competing interests

Nidhal Kahlaoui

Contribution

Competing interests

Sylvain Cardinaud

Contribution

Competing interests

Thibaut Naninck

Contribution

Competing interests

Benoit Delache

Contribution

Competing interests

Mathieu Surenaud

Contribution

Competing interests

Mathilde Galhaut

Contribution

Competing interests

Nathalie Dereuddre-Bosquet

Contribution

Competing interests

Mariangela Cavarelli

Contribution

Competing interests

Pauline Maisonnasse

Contribution

Competing interests

Mireille Centlivre

Contribution

Competing interests

Christine Lacabaratz

Contribution

Competing interests

Aurelie Wiedemann

Contribution