Skip to main content
Advertisement
  • Loading metrics

Longitudinal dynamics of SARS-CoV-2-specific cellular and humoral immunity after natural infection or BNT162b2 vaccination

  • Patricia Almendro-Vázquez,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain

  • Rocio Laguna-Goya ,

    Roles Conceptualization, Investigation, Methodology, Supervision, Visualization, Writing – original draft, Writing – review & editing

    rocio.laguna@salud.madrid.org

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Immunology, Hospital Universitario 12 de Octubre, Madrid, Spain, Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBER), Instituto de Salud Carlos III, Madrid, Spain

  • Maria Ruiz-Ruigomez,

    Roles Resources, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Internal Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Alberto Utrero-Rico,

    Roles Formal analysis, Investigation, Visualization, Writing – review & editing

    Affiliation Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain

  • Antonio Lalueza,

    Roles Resources, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBER), Instituto de Salud Carlos III, Madrid, Spain, Department of Internal Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain, Department of Medicine, School of Medicine, Universidad Complutense de Madrid, Madrid, Spain

  • Guillermo Maestro de la Calle,

    Roles Resources, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Internal Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Pilar Delgado,

    Roles Formal analysis, Investigation, Resources, Visualization, Writing – review & editing

    Affiliation Centro de Biologia Molecular Severo Ochoa, Consejo Superior de Investigaciones Cientificas (CSIC), Universidad Autonoma de Madrid, Madrid, Spain

  • Luis Perez-Ordoño,

    Roles Resources, Writing – review & editing

    Affiliation Department of Emergency Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Eva Muro,

    Roles Resources, Writing – review & editing

    Affiliation Department of Emergency Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Juan Vila,

    Roles Resources, Writing – review & editing

    Affiliation Department of Emergency Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Isabel Zamarron,

    Roles Resources, Writing – review & editing

    Affiliation Department of Emergency Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Miguel Moreno-Batanero,

    Roles Investigation, Methodology, Writing – review & editing

    Affiliation Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain

  • Marta Chivite-Lacaba,

    Roles Investigation, Writing – review & editing

    Affiliation Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain

  • Francisco Javier Gil-Etayo,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Immunology, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Carmen Martín-Higuera,

    Roles Investigation, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Clinical Microbiology, Hospital Universitario 12 de Octubre, Madrid, Spain

  • María Ángeles Meléndez-Carmona,

    Roles Investigation, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Clinical Microbiology, Hospital Universitario 12 de Octubre, Madrid, Spain

  • Carlos Lumbreras,

    Roles Resources, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBER), Instituto de Salud Carlos III, Madrid, Spain, Department of Internal Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain, Department of Medicine, School of Medicine, Universidad Complutense de Madrid, Madrid, Spain

  • Irene Arellano,

    Roles Investigation, Writing – review & editing

    Affiliation Centro de Biologia Molecular Severo Ochoa, Consejo Superior de Investigaciones Cientificas (CSIC), Universidad Autonoma de Madrid, Madrid, Spain

  • Balbino Alarcon,

    Roles Formal analysis, Resources, Writing – review & editing

    Affiliation Centro de Biologia Molecular Severo Ochoa, Consejo Superior de Investigaciones Cientificas (CSIC), Universidad Autonoma de Madrid, Madrid, Spain

  • Luis Miguel Allende,

    Roles Investigation, Resources, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Immunology, Hospital Universitario 12 de Octubre, Madrid, Spain, Department of Immunology, Ophthalmology and ENT, Universidad Complutense de Madrid, Madrid, Spain

  • Jose Maria Aguado,

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Internal Medicine, Hospital Universitario 12 de Octubre, Madrid, Spain, Unit of Infectious Diseases, Hospital Universitario 12 de Octubre, Madrid, Spain

  •  [ ... ],
  • Estela Paz-Artal

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    Affiliations Instituto de Investigación Sanitaria Hospital 12 de Octubre (imas12), Madrid, Spain, Department of Immunology, Hospital Universitario 12 de Octubre, Madrid, Spain, Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBER), Instituto de Salud Carlos III, Madrid, Spain, Department of Immunology, Ophthalmology and ENT, Universidad Complutense de Madrid, Madrid, Spain

  • [ view all ]
  • [ view less ]

Abstract

The timing of the development of specific adaptive immunity after natural SARS-CoV-2 infection, and its relevance in clinical outcome, has not been characterized in depth. Description of the long-term maintenance of both cellular and humoral responses elicited by real-world anti-SARS-CoV-2 vaccination is still scarce. Here we aimed to understand the development of optimal protective responses after SARS-CoV-2 infection and vaccination. We performed an early, longitudinal study of S1-, M- and N-specific IFN-γ and IL-2 T cell immunity and anti-S total and neutralizing antibodies in 88 mild, moderate or severe acute COVID-19 patients. Moreover, SARS-CoV-2-specific adaptive immunity was also analysed in 234 COVID-19 recovered subjects, 28 uninfected BNT162b2-vaccinees and 30 uninfected healthy controls. Upon natural infection, cellular and humoral responses were early and coordinated in mild patients, while weak and inconsistent in severe patients. The S1-specific cellular response measured at hospital arrival was an independent predictive factor against severity. In COVID-19 recovered patients, four to seven months post-infection, cellular immunity was maintained but antibodies and neutralization capacity declined. Finally, a robust Th1-driven immune response was developed in uninfected BNT162b2-vaccinees. Three months post-vaccination, the cellular response was comparable, while the humoral response was consistently stronger, to that measured in COVID-19 recovered patients. Thus, measurement of both humoral and cellular responses provides information on prognosis and protection from infection, which may add value for individual and public health recommendations.

Author summary

In this work we describe the prognostic value of early detection of SARS-CoV-2-specific T cell response in acute COVID-19, as patients without a prompt SARS-CoV-2-specific T cell response progress to severe COVID-19. We also show that the presence of specific T cells against SARS-CoV-2 when patients arrive to the emergency room is a protective factor against developing severe COVID-19, independently of the age and gender of the patient, which are two major known contributors to disease outcome. In addition, we show robust cellular and humoral immune responses persist 3 months after real-world vaccination.

Introduction

SARS-CoV-2 infection is responsible for the current coronavirus disease 2019 (COVID-19) pandemic. Most SARS-CoV-2-infected subjects show mild or no symptoms, however, an important proportion of patients require hospitalization and some of them die [1]. This disparity in disease severity has been associated with characteristics such as age, male sex and comorbidities [13]. The immune response elicited by SARS-CoV-2 in each individual could also be related to the clinical outcome of the infection. An exaggerated innate immune activation has been extensively described in patients with a severe course of the disease [48], while the role of adaptive immune response in protection from this infection is less understood.

SARS-CoV-2-specific T cells have been detected during convalescence and up to eight months post-symptom onset (PSO) in COVID-19 patients [914]. However, the development of the specific cellular immune response during the acute phase of infection has not yet been characterized in-depth, except for several studies, with a limited number of patients, indicating that early induction of SARS-CoV-2-specific T cells can limit COVID-19 severity [1517]. Humoral responses to SARS-CoV-2 have been described from acute infection to convalescence and up to 8 months post-infection. Levels of specific antibodies have been positively correlated with disease severity [11,18,19] and a delayed production of neutralizing antibodies has been associated with fatal COVID-19 [20]. The correct kinetics and balance between the two arms of the adaptive immune response may be crucial for the resolution of the infection.

The mRNA vaccine BNT162b2 was approved after showing safety and a 95% efficacy in preventing COVID-19 [21]. A two-dose regimen, 21 days apart, was established with the previous vaccine candidate BNT162b1 in order to develop neutralizing antibodies [22]. BNT162b2 has proven its effectiveness in reducing symptomatic SARS-CoV-2 infection [23,24]. There are studies on the immunogenicity of BNT162b2 which expand up to a month post-vaccination [2527] and on the maintenance of vaccine-induced humoral responses six months post-vaccination [28,29], however, data on the dynamics and long-term duration of SARS-CoV-2 specific T cells after real-world vaccination are still lacking.

Here we longitudinally characterize the adaptive immune responses during acute and convalescent COVID-19 phases in three well-defined patient cohorts of different severity, in recovered patients up to 7 months PSO and in uninfected subjects up to three months after full vaccination. Our acute infection cohort represents the largest of its kind, with 88 patients recruited before any medication was administered (in particular no immunomodulatory drugs) and has allowed us to statistically demonstrate that the presence of specific T cells against SARS-CoV-2 when patients arrive to the emergency room is a protective factor against developing severe COVID-19, independently of the age and gender of the patient. We also describe a correlation between higher numbers of specific T cells and lower viral load at diagnosis, which could explain the protection conferred by T cells. In addition, we characterize the vaccine-elicited adaptive immunity at five time-points and show robust cellular and humoral immune responses persist 3 months after real-world vaccination.

Results

Study design and participants

Cellular and humoral SARS-CoV-2-specific immune responses were analyzed in patients with natural infection, uninfected healthy controls and BNT162b2 vaccinated individuals. Samples were obtained upon arrival to the emergency room (ER), prior to any treatment administration, in patients with acute COVID-19 who were later classified according to their maximum disease severity into mild (N = 32), moderate (N = 34) and severe (N = 22) (see classification criteria in Fig 1A). Mild patients were not hospitalized, and moderate and severe patients required hospitalization. Clinical and demographic characteristics, WHO ordinal scale [30] and treatments are detailed in Table 1. Patients attended the ER a week PSO on average in the 3 groups. In moderate and severe patients, a follow-up sample was taken a week after hospitalization. In mild and moderate patients another sample was taken during the convalescent phase, a month PSO.

thumbnail
Fig 1. Overview of patient cohorts including the time of sample collection.

(A) Patients with natural infection classified according to their maximum disease severity. As depicted, samples were obtained upon arrival to the emergency room (ER), a week after hospitalization (Hosp.), and during the convalescent and recovered phases depending on each cohort. (B) Healthy controls SARS-CoV-2 uninfected with negative PCR and/or serology. (C) Vaccinated individuals who received two doses of BNT162b2 21 days apart. Samples were collected pre-vaccine, pre-boost and 15, 30 and 90 days after receiving the second vaccine dose. Out.: out-patients; In.: in-patients; ICU: Intensive care unit; PaO2/FiO2: Partial pressure of arterial oxygen/Fraction of inspired oxygen ratio; PSO: Post-symptom onset. Fig 1 created with BioRender.com.

https://doi.org/10.1371/journal.ppat.1010211.g001

thumbnail
Table 1. Demographic and clinical characteristics of the cohorts.

https://doi.org/10.1371/journal.ppat.1010211.t001

Samples from 234 COVID-19 recovered patients were obtained, between 4 and 7 months PSO. All recovered patients had been hospitalized during the acute phase with moderate or severe COVID-19 (WHO ordinal scale of 3–7). The use of corticosteroids in the recovered cohort was lower than in the acute infection cohort, consistent with a change in treatment guidelines from the first to the second pandemic wave (Table 1).

A 30 healthy control cohort included subjects with SARS-CoV-2 negative serology and no prior known contact with SARS-CoV-2 positive cases (Fig 1B). The operating characteristics (ROC) curve analysis with recovered patients and healthy controls allowed to determine the cut-off values of 25, 21 and 14 IFN-γ spots forming units (sfu)/106 PBMC as positive cellular responses against the domain S1 of spike (S1), membrane (M) and nucleocapsid (N) proteins respectively, with high sensitivity and specificity (S1 Fig).

Lastly, vaccinated individuals were 28 healthcare workers, without previous SARS-CoV-2 infection (Fig 1C). They received two 30 μg doses of BNT162b2 (Pfizer-BioNTech), 21 days apart. Samples were obtained pre-vaccine, pre-boost and 15, 30 and 90 days after completing vaccination. 22/28 (79%) vaccinees were females, and the median age was 40 (interquartile range [IQR] 26–57 years).

Development of SARS-CoV-2-specific cellular and humoral immunity during acute infection

SARS-CoV-2-specific T cells, IgG and neutralizing antibodies were sequentially analyzed in patients with natural infection, at their arrival to ER and during convalescence. Mild patients had a detectable, potent cellular immune response against S1, M and N, already within the first 2 weeks PSO (Figs 2A and S2). During this time period, a gradual increase in the specific cellular responses and positive correlations between the cellular response against the 3 proteins and days PSO (DPSO) were observed (p<0.001, Fig 2B). The peak of the SARS-CoV-2-specific cellular response occurred 2 weeks PSO in mild patients, and although the responses to S1, M and N were generally parallel and correlated positively (Fig 2B), the peak cellular response was more potent against S1 and N than M (median of 1052, 838 and 637 peak IFN-γ sfu/106 PBMC, respectively). The cellular immune response developed upon natural SARS-CoV-2 infection gave rise to detectable IFN-γ-, interleukin (IL)-2- and bifunctional IFN-γ/IL-2-producing T cells (S3 Fig). IFN-γ and IL-2 responses were similar in magnitude and correlated repeatedly in all patient cohorts (Fig 2B). For that reason, and for the sake of brevity, we will focus on data on IFN-γ mostly. Apart from the sfu, which correspond with the frequency of antigen-specific T cell clones, we obtained the relative spot volume (RSV) which represented the amount of secreted analyte (either IFN-γ or IL-2). The total amount of secreted analyte correlated repeatedly with the frequency of specific T cells (Fig 2B). After the peak of the cellular response, SARS-CoV-2-specific T cells decreased but were still detectable during convalescence, between 21 and 45 DPSO, with similar magnitude in the responses against S1, M and N (median of 202, 300 and 178 IFN-γ sfu/106 PBMC, respectively, Fig 2A). Most mild patients did not have detectable SARS-CoV-2-specific IgG upon arrival to ER. IgG detection started on day 13 PSO and increased during convalescence (p<0.05, Fig 2C). Anti-S1 IgG antibodies measured by ELISA correlated with anti-Spike IgG1 measured by a newly developed high-sensitivity flow cytometry method [31], especially during convalescence when the humoral response was fully developed (p<0.001, Fig 2B). In addition to quantifying SARS-CoV-2-binding antibodies, we measured their functionality in a classical pseudovirus neutralization assay. Antibodies from mild disease patients showed a modest neutralizing capacity (Fig 2D). Both anti-S1 IgG and anti-Spike IgG1 positively correlated with the neutralizing capacity of serum samples (r = -0.59, p<0.05 and r = -0.78, p<0.001, correlation with normalized infection, Fig 2E).

thumbnail
Fig 2. Cellular and humoral immune responses in mild natural infection, during acute and convalescent phases.

(A) FluoroSpot IFN-γ responses against SARS-CoV-2 S1, M and N proteins according to the days post-symptom onset (DPSO). Data is represented as spot forming unit (sfu) per million PBMC. Dashed lines represent the positivity cut-off established by using a non-infected control group: >25 sfu/106 PBMC for S1, >21 sfu/106 PBMC for M and >14 sfu/106 PBMC for N. (B) Correlation matrices of cellular and humoral immune variables analyzed. Only significant correlations (p<0.05) are represented with asterisks. Positive correlations appear in purple, and negative correlations appear in orange. The size and the colour gradient of the circle corresponds to the magnitude of the correlation. (C) Results of SARS-CoV-2 S1 IgG ELISA according to DPSO. Dashed line represents the established cut-off of positivity (ratio ≥1.1). (D) Neutralizing capacity of 1/500 diluted sera, tested with S protein-pseudotyped virus represented as percentage of normalized infection neutralized. (E) Correlation between neutralizing capacity and the semi-quantitative results of SARS-CoV-2 IgG ELISA (left) and anti-Spike IgG1 assessed by flow cytometry (right). Linear regressions were performed using Spearman’s rank test. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney or Wilcoxon signed rank tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. RSV: Relative Spot Volume; DP: IFN-γ/IL-2 double positive sfu.

https://doi.org/10.1371/journal.ppat.1010211.g002

Moderate patients had a detectable cellular immune response against S1, M and N, within the first 2 weeks PSO (Figs 3A and S2), which did not correlate with DPSO (Fig 3B). The peak of cellular immunity occurred at 21–24 DPSO, with a median of 838, 767 and 483 IFN-γ sfu/106 PBMC for S1, M and N responses respectively (Fig 3A). During convalescence (21 to 40 DPSO), S1, M and N-specific T cells remained detectable (median of 378, 370 and 107 IFN-γ sfu/106 PBMC, respectively, Fig 3A). Most moderate patients did not have detectable SARS-CoV-2-specific IgG upon arrival to ER (Fig 3C), however, they developed high antibody levels with a robust neutralizing capacity during convalescence (Fig 3D). There were strong correlations between antibody amounts and neutralization (r = -0.79, p<0.001 for anti-S1 IgG, and r = -0.87, p<0.0001 for anti-Spike IgG1, Fig 3E).

thumbnail
Fig 3. Cellular and humoral immune responses in moderate natural infection, during acute and convalescent phases.

(A) FluoroSpot IFN-γ responses against SARS-CoV-2 S1, M and N proteins according to the days post-symptom onset (DPSO). Dashed lines represent the positivity cut-off established by using a non-infected control group: >25 sfu/106 PBMC for S1, >21 sfu/106 PBMC for M and >14 sfu/106 PBMC for N. (B) Correlation matrices of cellular and humoral immune response variables analyzed. Only significant correlations (p<0.05) are represented with asterisks. (C) Results of SARS-CoV-2 S1 IgG ELISA according to DPSO. (D) Neutralizing capacity tested with S protein-pseudotyped virus represented as percentage of normalized infection neutralized. (E) Correlation between neutralizing capacity and the semi-quantitative results of SARS-CoV-2 IgG ELISA (left) and anti-Spike IgG1 assessed by flow cytometry (right). Linear regressions were performed using Spearman’s rank test. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney or Wilcoxon signed rank tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. (See Fig 2 footnote for more detailed information).

https://doi.org/10.1371/journal.ppat.1010211.g003

Severe patients had a remarkably low cellular immune response against all 3 proteins, within the first 2 weeks PSO (Figs 4A and S2), which did not correlate with DPSO either (Fig 4B). The development of the cellular response was very slow, as it had still not increased after a week of hospitalization (median of 24, 11 and 8 IFN-γ sfu/106 PBMC at ER and of 48, 25 and 33 IFN-γ sfu/106 PBMC after a week of hospitalization, for S1, M and N responses, respectively, Fig 4C). Detection of SARS-CoV-2-specific IgG upon arrival to ER was uneven: while most severe patients had undetectable IgG a few of them had already high levels (Fig 4D). During the first hospitalization week, approximately half of the patients with no antibodies developed them (Fig 4E), and those antibodies were highly neutralizing (Fig 4F). Severe patients showed the strongest correlation between levels of antibodies and neutralizing response (r = -0.90, p<0.001 for anti-S1 IgG, and r = -0.87, p<0.001 for anti-Spike IgG1, Fig 4G). Eight of the 22 (36%) severe patients died during hospitalization. These patients who died had significantly lower cellular response against all 3 proteins (Fig 4H) and a trend towards lower humoral response (Fig 4I) than severe patients who survived. In summary, upon arrival to ER, severe patients showed a negative correlation between antibody levels and some compartments of the cellular response, which reflected a discoordination between both arms of the immune response and the dominance of the humoral over the cellular immunity (Fig 4B).

thumbnail
Fig 4. Cellular and humoral immune responses in severe natural infection, during acute phase, measured at emergency room (ER) and a week after hospitalization (1W Hosp.).

(A) FluoroSpot IFN-γ responses against SARS-CoV-2 S1, M and N proteins according to the days post-symptom onset (DPSO). Dashed lines represent the positivity cut-off established by using a non-infected control group: >25 sfu/106 PBMC for S1, >21 sfu/106 PBMC for M and >14 sfu/106 PBMC for N. (B) Correlation matrices of cellular and humoral immune variables analyzed. Only significant correlations (p<0.05) are represented with asterisks. (C) FluoroSpot IFN-γ responses against SARS-CoV-2 S1, M and N proteins at ER and a week after hospitalization. (D-E) Results of SARS-CoV-2 S1 IgG ELISA according to DPSO. (F) Neutralizing capacity tested with S protein-pseudotyped virus represented as percentage of normalized infection neutralized. (G) Correlation between neutralizing capacity and the semi-quantitative results of SARS-CoV-2 IgG ELISA (left) and anti-Spike IgG1 assessed by flow cytometry (right). (H-I) Comparison of S1, M and N T-cell (H) and anti-S1 IgG (I) responses between survivor and non-survivor patients (black stars). Linear regressions were performed using Spearman’s rank test. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney, Wilcoxon signed rank or Kruskal-Wallis tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. (See Fig 2 footnote for more detailed information).

https://doi.org/10.1371/journal.ppat.1010211.g004

Association between dynamics of SARS-CoV-2-specific immune responses and COVID-19 severity

We next compared the dynamics of SARS-CoV-2 immunity among the 3 cohorts with different severity. Specific T cells started to increase in mild and moderate patients in the first week of symptoms, when they were barely detectable in severe patients. During the second week PSO, the cellular response developed robustly in mild patients with a significant increase between the two weeks (p<0.01, Fig 5A), while it remained unchanged in moderate and severe patients. From the start of infection, we observed that the cellular response developed faster in mild patients, reaching their peak approximately 14 DPSO (Fig 2A), than in moderate patients, who reached their peak response a week later (Fig 3A). Nonetheless, the magnitude of the peak S1, M and N-specific cellular response was similar in mild and moderate patients. Severe patients, however, had not started developing high numbers of specific T cells by the end of the acute infection follow-up, by 21 DPSO (Fig 4C). In addition to developing specific T cell clones faster and in higher numbers, specific T cells from milder patients were more functional than those from severe patients, as shown by a higher amount of IFN-γ and IL-2 secreted per clone (Figs 5B and S3). In our cohort, we did not observe an association between viral load at diagnosis and disease severity (Fig 5C). Nevertheless, we found a strong correlation, mostly in moderate patients (r = 0,98, p<0.0001), between a more potent specific T cell response and a higher Ct value at diagnosis, suggesting that patients with stronger specific T cell responses are more able to control viral replication (Fig 5D). These results confirm and further validate smaller studies which had indicated that the development of a functional T cell response is associated with mild disease.

thumbnail
Fig 5. Association between dynamics of SARS-CoV-2-specific immune responses and COVID-19 severity.

(A) SARS-CoV-2-specific IFN-γ-producing T-cell responses reactive to the S1, M and N proteins in mild, moderate, and severe patients during the first (W1) and the second (W2) week post-symptom onset. (B) Relative Spot Volume (RSV) per spot of secreted IFN-γ after S1, M and N peptide pool stimulation in mild, moderate, and severe patients during acute infection. (C) Ct values obtained on real time RT-PCR for detection of the E gene upon arrival to ER in mild, moderate and severe patients. A higher Ct value corresponded to a lower viral load. (D) Correlation between S1 IFN-γ-producing T-cells and the relative viral load represented as Ct value at ER. (E) Semi-quantitative results of SARS-CoV-2 S1 IgG ELISA in the acute phase in mild, moderate, and severe patients. Dashed line represents the established cut-off of positivity (ratio ≥1.1). (F) Neutralizing capacity tested with S protein-pseudotyped virus represented as normalized infection. (G) Longitudinal data on the dynamics of cellular and humoral responses in six representative patients from the 3 severity cohorts. (H) Percentage of positive S1 cellular and humoral response during acute phase according to the days post symptom onset (DPSO) in mild, moderate and severe patients. (I) Comparison of IFN-γ responses against SARS-CoV-2 S1, M and N proteins between mild and moderate patients in convalescent phase. (J) Comparison of semi-quantitative results of SARS-CoV-2 S1 IgG between mild and moderate patients in convalescent phase. (K) Comparison of neutralizing capacity tested with S protein-pseudotyped virus represented as normalized infection. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. Linear regressions were performed using Spearman’s rank test. The significance between groups was determined using Mann Whitney, Wilcoxon signed rank or Kruskal-Wallis tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. ER: Emergency Room; Conv: Convalescence; 1W Hosp.: 1 Week Hospitalization; DPSO: Days post-symptom onset.

https://doi.org/10.1371/journal.ppat.1010211.g005

Both SARS-CoV-2-specific IgG and neutralization were mostly undetectable in all 3 cohorts during the first 2 weeks of symptoms (Fig 5E and 5F). The different dynamics of cellular and humoral immune response development are shown with 6 representative patients from each of the 3 cohorts in Fig 5G. On the whole, the adaptive immune response seemed to be more coordinated in mild and moderate patients, as most patients had developed a positive cellular response within the first 2 weeks PSO (23/32 [72%] and 26/34 [76%] respectively, Fig 5H), which was followed by the development of virus-specific IgG (86% mild and 100% moderate patients had detectable antibodies by day 13–15 PSO). Most severe patients, however, had not developed a positive SARS-CoV-2-specific cellular response within the first 2 weeks PSO, they lacked any correlation between positivity of the cellular response and DPSO when they attended the ER, and the development of virus-specific IgG did not follow a temporal pattern during acute infection (Fig 5H). By investigating if the SARS-CoV-2 specific cellular response could have a prognostic value, we found that a higher frequency of ER-measured specific T cells reduced the probability of developing severe COVID-19 (odds ratio [OR] per 100 IFN-γ sfu/106 PBMC increment: 0.45, 95%CI 0.20–0.76, p<0.05). Given that age and sex differed among mild, moderate and severe patients (Table 1), and that both older age and male sex have been associated with disease severity [13], we performed a multivariate analysis to test if the decreased risk of developing severe COVID-19 associated to high frequency of specific T cells could be influenced by these two variables. Independently of age and sex, the S1-specific cellular response measured at ER remained as a protective factor (OR per 100 IFN-γ sfu/106 PBMC increment: 0.47, 95%CI 0.20–0.87, p<0.05), while age and male sex appeared as independent COVID-19 severity risk factors (OR: 1.05, 95%CI 1.01–1.10, p = 0.01 and OR: 5.02, 95%CI 1.31–26.07, p<0.05, respectively).

During convalescence, there was no significant difference in the cellular response against all 3 proteins between mild and moderate patients (Fig 5I). On the contrary, antibodies were much higher in moderate than in mild patients (p<0.0001, Fig 5J). Accordingly, mild patients showed modest neutralization capacity while moderate patients had a significantly more robust neutralizing response (p<0.01, Fig 5K).

SARS-CoV-2 infection is known to cause lymphopenia, which associates with COVID-19 severity [32,33]. We found lymphopenia in acute moderate patients, mainly due to reduced T cells (S1 Table). There was a substantial decrease in CD4+ and CD8+ T cells, and a milder reduction in B cells. All lymphocyte populations recovered during convalescence. Severe patients had a profound lymphopenia during acute infection, with a reduction in CD4+, CD8+, B and NK cells. In addition, the overall T cell functionality was clearly reduced in severe patients when cells were polyclonally stimulated with anti-CD3/anti-CD28 (S4 Fig). We asked if the lymphopenia and/or the suboptimal global T cell response during acute infection could be responsible for the lack of SARS-CoV-2-specific T cells in severe patients. When we adjusted in each sample the specific cellular response by the number of sfu/106 PBMC after polyclonal stimulus, the results showed again how mild patients developed a significantly faster and more potent SARS-CoV-2-specific cellular response than hospitalized patients in the beginning of the acute infection (S4 Fig), suggesting that the profound lymphopenia and general reduction in T cell response in severe patients could only partially contribute to the lack of SARS-CoV-2-specific cellular response.

Duration of SARS-CoV-2-specific cellular and humoral immunity after natural infection

In order to assess the duration of adaptive immunity after natural infection, we analyzed the cohort of 234 COVID-19 recovered patients (Fig 1A). Most recovered patients maintained detectable virus-specific memory T cells, significantly more against S1 than M and N (median of 140, 117 and 75 IFN-γ sfu/106 PBMC respectively, p<0.05, Fig 6A). In total, 227/234 (97%) patients maintained a cellular response against S1, M and/or N above the positivity threshold. Unlike in the acute and convalescent phase, the number of IFN-γ- and IL-2-producing T cells differed in recovered patients. Six months PSO we detected significantly more IL-2-producing than IFN-γ-producing clones (174 vs 140, 152 vs 117 and 92 vs 75 sfu/106 PBMC against S1, M and N, respectively, all p<0.0001, Fig 6B). Moreover, we found an unexpected gradual increase in the cellular response from month 4 to month 7 PSO, which was significant for S1 (p<0.01 Fig 6C), and positive correlations between the cellular response and DPSO (Fig 6D). Recovered patients had normalized their lymphocyte populations (S1 Table), however they showed an increased polyclonal T cell response compared to healthy controls (p<0.001, S5 Fig), which augmented with time after infection. After adjusting by the sfu/106 PBMC values of polyclonal responses, the increase in SARS-CoV-2-specific cellular response was considerably reduced, suggesting that the observed increase in SARS-CoV-2-specific cellular response from month 4 to 7 PSO was likely due to a general increment in T cell functionality rather than to a specific boost of cellular response (S5 Fig).

thumbnail
Fig 6. Duration of SARS-CoV-2-specific cellular and humoral immunity after natural infection.

(A) S1, M and N cellular responses in recovered patients. (B) Comparison between IFN-γ and IL-2-producing T-cell responses. (C) IFN-γ-producing T-cell responses according to the months post-symptom onset (PSO). (D) Correlation matrix of cellular and humoral immune variables analyzed. Only significant correlations (p<0.05) are represented with asterisks. (E) SARS-CoV-2 S1 IgG according to the months PSO. Dashed line represents the established cut-off of positivity (ratio ≥1.1). (F) Neutralizing capacity represented as normalized infection according to the months PSO. (G) Correlation between neutralizing capacity and the semi-quantitative results of SARS-CoV-2 IgG ELISA (left) and anti-Spike IgG1 measured by flow cytometry. (H-J) S1 IFN-γ-producing T-cells (H), semi-quantitative results of SARS-CoV-2 S1 IgG (I) and neutralizing capacity (J) in recovered patients who suffered a moderate or severe acute disease. Linear regressions were performed using Spearman’s rank test. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney or Kruskal-Wallis tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. RSV: Relative Spot Volume; DP: IFN-γ/IL-2 double positive sfu.

https://doi.org/10.1371/journal.ppat.1010211.g006

A robust humoral response anti-COVID-19 was still detected in recovered patients, as 99% (231/234) of them were positive for SARS-CoV-2-specific IgG. Contrarily to the cellular response, there was a gradual significant decline from month 4 to month 7 in antibodies (p<0.01, Fig 6E) and neutralizing capacity (p<0.0001, Fig 6F), and a negative correlation was observed between anti-S1 IgG levels and DPSO (Fig 6D). Antibody levels paralleled the neutralizing capacity (r = -0.50, p<0.05 for anti-S1 IgG and r = -0.69, p<0.001 for anti-Spike IgG1, Fig 6G), although these correlations were weaker than in acute and convalescent patients (Figs 2E, 3E and 4G).

There was no difference in memory cellular response between recovered patients who had suffered a moderate (N = 203) or a severe (N = 31) COVID-19 form (Fig 6H), while subjects recovered from severe COVID-19 maintained significantly higher antibody levels (p<0.05, Fig 6I), with higher neutralization capacity (p<0.05, Fig 6J). The absence of specific T cells in patients with acute severe infection (Fig 4), which were detectable in subjects recovered from severe COVID-19 suggests that severe patients may develop a specific cellular response sometime during convalescence.

Development of SARS-CoV-2-specific cellular and humoral immunity after BNT162b2 vaccination

We studied the development of the adaptive immunity after vaccination in a cohort of individuals naïve for SARS-CoV-2 infection. After the first vaccine dose, 82% (23/28) of vaccinees developed a positive cellular response, which as expected was specific against S1 (median of 140 IFN-γ sfu/106 PBMC, Fig 7A). The administration of the second boost vaccine dose potentiated the cellular response, which peaked 2 weeks after full vaccination (median of 621 IFN-γ sfu/106 PBMC) and then decreased steadily a month and three months post-vaccination (median of 439 and 158 IFN-γ sfu/106 PBMC, respectively, Fig 7B). Despite this decrease, all vaccinees except one remained above the positivity threshold. A similar dynamic was observed in the development of IL-2-producing SARS-CoV-2-specific T cells and IFN-γ/IL-2 double positive T cells (Fig 7C and 7D). There was a strong positive correlation between IFN-γ, IL-2 and bifunctional T cell response, measured as the number of specific-T cell clones (sfu) and as the total amount of cytokine secreted (RSV) after the first and second vaccine doses (Fig 7E). Additionally, as reported in the phase I/II clinical trial for the vaccine candidate BNT162b1[22], we confirmed that the BNT162b2 vaccine elicited exclusively Th1 T cell responses as noted by the increased frequency of CD4 Th1 PD1+, but not CD4 Th2 or CD4 Th17 PD1+, cells 2 weeks after full vaccination (p<0.01, Fig 7F).

thumbnail
Fig 7. Development of SARS-CoV-2-specific cellular and humoral immunity after BNT162b2 vaccination.

(A) IFN-γ-producing T-cell responses against S1, M and N proteins pre-vaccine and pre-boost. (B-D) S1 SARS-CoV-2-specific IFN-γ (B), IL-2 (C) and both cytokines (D) producing T-cells pre-vaccine, pre-boost and 15, 30 and 90 days after complete vaccination. Dashed line represents the positivity cut-off: >25 IFN-γ sfu/106 PBMC. (E) Correlation matrix of adaptive immune variables analyzed. Only significant correlations (p<0.05) are represented with asterisks. (F) Frequency of CD4 Th1, Th2 and Th17 PD1+ cells pre-vaccine, pre-boost and 15 days after complete vaccination. (G) Anti-S1 IgG results pre-vaccine, pre-boost and 15, 30 and 90 days after complete vaccination. Dashed line represents the positivity cut-off (ratio ≥1.1). (H) Neutralization capacity pre-vaccine, pre-boost and 15 and 30 days after complete vaccination. (I) Correlation between neutralization and anti-S1 IgG (left) and anti-Spike IgG1 (right). (J-L) Correlation between age and S1 IFN-γ-producing T-cells (J), anti-S1 IgG (K) and neutralization (L) 15 days after full vaccination when the vaccine-elicited immune response peaked. (M) Comparison of S1 T-cell response between SARS-CoV-2 natural infection and complete vaccination. (N-O) Comparison of anti-S1 IgG (N) and neutralization (O) between SARS-CoV-2 natural infection and complete vaccination. Linear regressions were performed using Spearman’s rank test. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney, Wilcoxon signed rank or Kruskal-Wallis tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. DPSO: Days post symptom onset.

https://doi.org/10.1371/journal.ppat.1010211.g007

The development of anti-S1 antibodies after vaccination paralleled the dynamics of the cellular response. A single vaccine dose generated detectable IgG in 94% (26/28) of vaccinees (Fig 7G). Two weeks after completing vaccination, there were extremely high antibody levels in all individuals, which remained stable a month after vaccination. Three months after vaccination there was a stark decrease in IgG levels (p<0.0001, Fig 7G), although all vaccinees remained well above the positivity threshold. As described by the vaccine trials [21,22], a single dose of vaccine was insufficient to develop neutralizing antibodies (Fig 7H). After the boost dose, sera from vaccinees were able to partially neutralize infection, although the neutralization capacity decreased already at a month, compared to two weeks, after full vaccination (p<0.01, Fig 7H). Similar to natural infection, there was a correlation between levels of antibodies and neutralizing response (r = -0.88, p<0.0001 for anti-S1 IgG, and r = -0.70, p<0.0001 for anti-Spike IgG1, Fig 7I). Given the parallel cellular and humoral response to the vaccine, a month after full vaccination, there was a positive correlation between the two arms of the adaptive immunity (Fig 7E). While no relationship between age and magnitude of the cellular response was observed (Fig 7J), age negatively correlated with antibody levels (r = -0.56, p<0.01, Fig 7K), despite of which all vaccinees showed high levels of anti-S1 IgG, with no differences in neutralizing capacity among subjects of different ages (Fig 7L).

Comparison of SARS-CoV-2-specific immunity development after natural infection or real-world vaccination

Finally, we compared the magnitude and dynamics of the adaptive immune response after vaccination to that developed after natural infection. The peak of the cellular response 2 weeks after complete vaccination was comparable to the peak response in acute mild and moderate patients, 2 to 3 weeks PSO (Fig 7M). After the peak, the cellular response started to decline and 3 months after vaccination it was similar to that of COVID-19 recovered patients 4–7 months after infection. The production of anti-S1 IgG reached much higher levels in vaccinees than in naturally infected individuals, either convalescent or recovered patients, up to three months post-vaccination (p<0.0001, Fig 7N). However, the neutralization capacity of the antibodies elicited by the vaccine was relatively lower than that of antibodies elicited by the natural infection. Despite having higher levels of IgG, the neutralization capacity of vaccinees was intermediate between that of mild and moderate convalescent patients (Fig 7O).

Discussion

Our analysis of cellular and humoral SARS-CoV-2-specific immunity in patients with varying degrees of COVID-19 severity indicates that SARS-CoV-2-specific T cells are essential for disease control and that the delayed mounting of SARS-CoV-2-specific cellular immune response is associated with severe infection course. In addition, we describe in real-world conditions the development of a robust, coordinated cellular and humoral immune response with the BNT162b2 mRNA COVID-19 vaccine, up to three months after completing the 2-dose vaccination schedule.

To date, this is the largest study analyzing the SARS-CoV-2-specific cellular and humoral immune response in patients upon ER arrival, which allowed the characterization of both arms of adaptive immunity from the start of natural infection without any interference from immunomodulatory medication. The timing in which the different responsive compartments of immunity appeared was crucial for determining the course of the infection. Patients who developed a prompt IFN-γ- and IL-2-producing cellular response followed by antibody production had a mild COVID-19 course, an observation in agreement with the well-known concept that the early development of virus-specific T cells and the subsequent production of neutralizing antibodies are required for effective control of viral infections [34]. On the other extreme, patients who did not develop an early specific cellular response and initiated a humoral immune response with subsequent production of high levels of antibodies, were not able to eliminate the virus efficiently and developed severe symptoms. A strong humoral response in severe COVID-19 patients has already been reported [11,16,35]. We observed that the inadequate adaptive immune response developed in severe patients could be explained partially by the profound lymphopenia and general decline in the T cell response. However, it is likely that the differences observed in the adaptive immune response between mild, moderate and severe COVID-19 forms also derive from differences in the patients’ initial innate immune response. In particular, type I interferons, which are important for anti-viral host defence and may be required to initiate cell-mediated immune responses [36], could be hindered by genetic defects or autoantibodies that have been associated with severe COVID-19 [37,38]. Antibodies, and their neutralizing capacity, were not associated with reduced disease severity in our acute patients, which was also observed in Middle East respiratory syndrome coronavirus (MERS-CoV)-infected patients [39]. Antibodies which may not be efficient at clearing the acute viral infection, may well provide protection against infection when present before exposure to the virus [40,41].

In our study, robust antigen-specific cellular and IgG responses developed concomitantly after vaccination, unlike in natural infection. The fast and potent immunity triggered by this mRNA vaccine explains its high effectiveness in preventing SARS-CoV-2 infection in the real world [23,24]. As observed in the clinical trial by Sahin et al. [22], the BNT162b2 vaccine elicited higher antibody levels, although with inferior neutralizing capacity than those formed after suffering moderate COVID-19. A current concern is the immunogenicity of SARS-CoV-2 vaccines in elderly. Our results in adults 22 to 64 years old show no negative impact of age in the cellular response elicited by the vaccine and a decrease in the amount of SARS-CoV-2-specific IgG, but not in neutralizing activity, associated with age. While some studies report weak humoral and cellular immunogenicity of the BNT162b2 vaccine in COVID-19 uninfected nursing home residents [42], others show optimal antibody titers irrespective of age [43] and, more importantly, a reduction in SARS-CoV-2 infections in vaccinated elderly [44]. Finally, three months after complete vaccination, the SARS-CoV-2-specific cellular response was similar, and the antibody levels were higher than in COVID-19 recovered patients. This coincides with the recent report by Tan et al [45] in which a similar magnitude of T cell responses between the vaccinees and recovered patients was observed three months after vaccination or infection. It could therefore be speculated that the protection due to immunization would last similarly to that conferred by natural infection. Studies on the long-term immunogenicity and effectiveness of COVID-19 vaccines will be informative to design vaccination strategies in the future. So far, several publications have investigated the durability of vaccine-induced humoral responses up to 6 months after vaccination. A gradual and substantial decrease in IgG specific antibodies and neutralizing antibody levels has been reported through the first half a year post-vaccination [28,46,47], independently of having received the mRNA based BNT162b2 or the adenoviral vector-based ChAdOx1 nCoV-19 [48]. Both, BNT162b2-induced SARS-CoV-2 antibody and T-cell responses have been shown to experience a more important decline after 6 months in the elderly than in young health-care workers [49]. Although these observations support the recommendation for booster vaccinations, the effect of antibody waning needs to be clearly determined, as, importantly, specific memory B cells remain 6 months after the immunization [29] or after the natural infection, even in the absence of detectable neutralizing antibodies [50].

A robust SARS-CoV-2-specific cellular and humoral immunity remained months after recovery. We initially observed an increase in SARS-CoV-2-specific T cells from month 4 to month 7 PSO, which paralleled the general increase in T cell functionality upon immune reconstitution following severe lymphopenia. The relativization of the specific response to the overall T cell function allowed a more accurate interpretation of the results obtained. Importantly, SARS-CoV-2-specific T cell immunity remained stable from month 4 to month 7, while specific IgG declined over this period. This information from recovered patients further supports our finding in acute infection that specific T cells are essential in infection control and will likely be essential for long-term COVID-19 protection. Furthermore, Sekine et al. [51] described the presence of SARS-CoV-2-specific T cells in seronegative exposed and convalescent individuals, implying that only measuring specific IgG may underestimate protection against COVID-19. In addition, recovered patients had significantly more IL-2-producing SARS-CoV-2-specific T cells than IFN-γ-producing ones. The higher prevalence of IL-2-producing SARS-CoV-2-specific cell clones could support the predominance of CD4+, rather than CD8+, T lymphocytes as expression of immune memory in patients who overcome COVID-19 [35].

A limitation of this study is that the FluoroSpot method does not differentiate between CD4- and CD8-SARS-CoV-2-specific T cells. Further studies to understand the degree of each cell subset participation in the response after infection or vaccination will be important. It should also be noted that the T cell response was analyzed using peptides covering the S1 region, which contains the receptor binding domain, but the potential response against S2 region was not analyzed. In addition, we proposed that severe disease develops because of the lack of specific cellular response to contain the virus. Elucidating the mechanisms leading to the failure in mounting an early, robust immune response against SARS-CoV-2 in some patients may help to develop new therapies and improve their prognosis.

In summary, our results not only confirm that the prompt induction of SARS-CoV-2-specific T cells in naturally infected patients is essential for disease control, but also that the SARS-CoV-2 specific cellular response measured at ER has a prognostic value as it is an age- and sex-independent protective factor against developing severe COVID-19. In addition, specific cellular responses remained stable in recovered patients and were present up to three months after vaccination. These results highlight the relevance of monitoring SARS-CoV-2-specific cellular immune responses, and not only antibody levels, for prognosis in natural infection, and as a correlate for protection after infection and after vaccination.

Materials and methods

Ethics statement

All participants were enrolled after signing the informed written consent of two studies approved by the institutional clinical research ethics committee (references 21/039 and 21/056).

Study design

A total of 380 subjects were recruited at Hospital Universitario 12 de Octubre (Madrid, Spain). These subjects belonged to six different cohorts (see Fig 1A and Study design and participants section in Results). Natural infection was confirmed by a positive SARS-CoV-2 real time reverse transcription polymerase chain reaction (RT-PCR) or a SARS-CoV-2 antigen test in all COVID-19 patients.

FluoroSpot assay

Peripheral blood mononuclear cells (PBMCs) were freshly isolated by density-gradient centrifugation using Ficoll-Paque within a maximum of 12 hours from blood being drawn and were seeded in duplicate at 300,000 cells/well in IFN-γ IL-2 FluoroSpot plates (MabTech). Cells were supplemented with 15-mer overlapping peptides covering the S1 domain of the S glycoprotein (SARS-CoV-2 S1 scanning pool, MabTech), the N protein (Epitope Mapping Peptide Set [EMPS] SARS-CoV-2 NCAP-1, JPT), and the M protein (EMPS SARS-CoV-2 VME1, JPT) at a final concentration of 1 μg/mL. Negative control wells contained anti-CD28 mAb (1 μg/mL) and positive control wells included anti-CD3 mAb (MabTech). Assays were incubated for 16–18 hours at 37°C. Spots were counted using an automated IRIS FluoroSpot Reader System (MabTech). To quantify antigen-specific responses, spots of the negative control wells were subtracted from the mean spots test wells, and the results were expressed as IFN-γ or IL-2 producing spot forming units (sfu) per 106 PBMCs. Results were excluded if negative control wells had >80 IFN-γ sfu/106 PBMCs or positive control wells had <400 IFN-γ sfu/106 PBMCs. Reponses were considered positive if the results were at least three times higher than the mean of the negative control wells and above of the following antigen-specific cut-off values: >25  IFN-γ sfu/106 PBMCs for the S1 domain of the S glycoprotein, >14  IFN-γ sfu/106 PBMCs for the N protein, and >21  IFN-γ sfu/106 PBMCs for the M protein, as previously published [52].

Anti-S1 IgG detection by ELISA

Serum SARS-CoV-2 IgG antibodies targeting the S1 protein were detected with the Euroimmun Anti-SARS-CoV-2 enzyme-linked immunosorbent assay (ELISA) (Euroimmun AG, Lübeck, Germany) according to manufacturer’s instructions. Optical density (OD) values were measured at 450 nm using the PR 3100 microplate reader (Bio-Rad Life Science, Marnes-La-Coquette, France). Results were evaluated semi-quantitatively by calculating the ratio of the OD value of the sample over the OD value of the calibrator (relative OD), with the following cut-off values: ratio <0.8: negative; ratio ≥ 0.8 to <1.1: borderline; and ratio ≥1.1: positive.

Anti-Spike IgG1 detection by flow cytometry

Jurkat cells co-expressing full-lenght Spike S protein of SARS-CoV2 and a truncated version of human EGFR (huEGFRt) as reporter were used to detect anti-spike IgG1 in human sera as described in Horndler et al[31]. Briefly, for each serum to test 120 x 103 cells were stained with a 1/50 serum dilution followed by a secondary anti-human IgG1 PE (clone HP6001, Southern Biotech) and an anti-human EGFR Bv421 (clone AY13, BioLegend) and the viability marker 7AAD. Samples were analyzed on a FACSCanto II flow cytometer and FlowJo software (Becton-Dickinson). The IgG1 PE/EGFR Bv421 ratio was used to determine the relative amount of anti-spike IgG1 signal.

Neutralization assay with pseudotyped virus

Neutralization assays were performed as described in Horndler et al[31]. Briefly, lentiviruses were produced by co-transfecting plasmids pCMV-dR (gag/pol), pHRSIN-GFP and a truncated S envelope (pCR3.1-St). A total of 35 × 103 ACE2+ HEK293T cells per well in a 48-well plate were seeded the day before transduction. Serially diluted serum was incubated with viral supernatant for 1 h at 37°C prior addition to the cells. Polybrene (8 μg/ml) was added and plates were centrifuged for 70 min at 1,600 g and left in culture for 48 h and then were resuspended in PBS with 2% FBS and 5 mM EDTA and fixed with 2% paraformaldehyde. GFP+ cells were then analyzed on a FACSCanto II flow cytometer, and the data were analyzed with FlowJo software (Becton-Dickinson). Samples were run in duplicates and infection was normalized to 100% with viral supernatant without serum. A decrease in the neutralization capacity was observed as patient sera were serially diluted 1/50, 1/500 and 1/2000 (S6 Fig). The serum dilution 1/500 was chosen for comparison of neutralizing capacity between cohorts as it allowed the maximum discrimination. A known convalescent serum was used as positive control and a pre-pandemic serum was used as negative control.

Viral load quantification

SARS-CoV-2 viral load at diagnosis was assessed using cycle threshold (Ct) values from a RT-PCR assay applied to nasopharyngeal swab samples. The RT-PCR assay was a laboratory-developed test (LDT) based on real time RT-PCR for detection of the E gene on the Panther Fusion Hologic (San Diego, CA, USA) using its open access functionality as previously described [53,54]. Amplification Ct values were considered a relative measure of viral load quantification. Lower Ct values corresponded with higher viral loads.

Peripheral blood lymphocyte populations

EDTA-treated whole blood (50 μL) was stained with 20 μL of BD Multitest 6-color TBNK reagent in Trucount tubes (BD Biosciences. San Jose, CA, USA) for 15 min. After lysis of red blood cells, T, B and NK lymphocyte subsets were enumerated with a FACSCanto II flow cytometer, and data were analyzed using FACSDiva software (BD Biosciences) [55].

For the analysis of T helper cell subsets after vaccination, EDTA-treated whole blood was incubated with: anti-CD3-PerCPy5.5 and anti-CD4-KO (Beckman Coulter); anti-CXCR3-PE, anti-CCR6-PB and anti-PD-1-PECy7 (BD Biosciences). Acquisition was performed with a Navios flow cytometer (Beckman Coulter), and data were analyzed with FlowJo software v10.6.2. Th1 cells were defined as CD4+CXCR3+CCR6-, Th2 cells as CD4+CXCR3-CCR6- and Th17 cells as CD4+CXCR3-CCR6+. PD-1 was considered as a recent activation marker.

Statistical analysis

Quantitative data were shown as the median with IQR, and qualitative variables were expressed as absolute and relative frequencies. Non-parametric Mann-Whitney U or Wilcoxon signed-rank tests were applied for comparison within two groups, when necessary. Kruskal-Wallis test was used to compare three or more unmatched groups. Correlations between continuous variables were evaluated using Spearman’s rank test. Operating characteristics (ROC) curve analysis was performed to determine positivity cut-off values for cellular response. Univariate and multivariate logistic regressions were performed to test association with disease severity. Differences were considered statistically significant when p<0.05. Statistical analysis was performed using GraphPad Prism version 8.0 software (GraphPad Software Inc, LaJolla, CA) and R software v4.0.3.

Supporting information

S1 Fig. Receiver operating characteristics (ROC) curve analysis and positivity cut-off values for S1, M and N cellular response.

Cut-off values were established by using a control group of 30 healthcare workers with no microbiological or clinical evidence of SARS-CoV-2 infection and 234 recovered COVID-19 patients. (A) S1 Cut-off >25 IFN-γ sfu/106 PBMC (AUC = 0.97, 95%CI = 0.95–0.99). (B) M Cut-off >21 IFN-γ sfu/106 PBMC (AUC = 0.97, 95%CI = 0.95–0.99). (C) N Cut-off >14 IFN-γ sfu/106 PBMC (AUC = 0.96, 95%CI = 0.94–0.98). AUC: Area Under the Curve; IC: Confidence Interval.

https://doi.org/10.1371/journal.ppat.1010211.s001

(TIF)

S2 Fig. Representative T cell IFN-γ FluoroSpot images.

FluoroSpot IFN-γ responses against SARS-CoV-2 S1, M and N peptide pools and positive (anti-CD3 and anti-CD28 mAb) and negative (anti-CD28 mAb only) control wells for mild (A), moderate (B) and severe (C) patients.

https://doi.org/10.1371/journal.ppat.1010211.s002

(TIF)

S3 Fig. Development of SARS-CoV-2-specific cellular response in natural infection.

IL-2 and bifunctional IFN-γ+IL-2 T cell responses against SARS-CoV-2 S1, M and N protein pools according to the days post symptom onset during acute and convalescent phase in mild (green) (A-B) moderate (orange) (C-D) and severe (red) (E-H) patients. Data is represented as spot forming unit (sfu) per million PBMC. (I-K) Relative Spot Volume (RSV) per spot of secreted IL-2 after S1, M and N peptide pool stimulation in mild, moderate, and severe patients during acute infection. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney, Wilcoxon signed rank or Kruskal-Wallis tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

https://doi.org/10.1371/journal.ppat.1010211.s003

(TIF)

S4 Fig. Analysis of the SARS-CoV-2- specific cellular response relative to the overall T cell functionality.

(A) FluoroSpot IFN-γ responses against anti-CD3mAb polyclonal stimulus in mild (green), moderate (orange) and severe (red) COVID-19 patients. Data is represented as spot forming unit (sfu) per million PBMC. (B-D) S1, M and N cellular response adjusted by anti-CD3 sfu/106PBMC normalization according to the days post symptom onset during acute and convalescent phase in mild (green), moderate (orange) and severe (red) COVID-19 patients. Data is shown as relative sfu per million PBMC. (E) Relative SARS-CoV-2-specific IFN-γ-producing T-cell responses reactive to the S1, M and N proteins in mild, moderate, and severe patients during the first (W1) and the second (W2) week post symptom onset. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney, Wilcoxon signed rank or Kruskal-Wallis tests, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

https://doi.org/10.1371/journal.ppat.1010211.s004

(TIF)

S5 Fig. Hyperactivation state of T cell after natural infection.

(A) Comparison of FluoroSpot IFN-γ responses against anti-CD3mAb polyclonal stimulus in healthy controls (HC) and moderate or severe recovered patients (Post-COVID). (B) anti-CD3mAb IFN-γ responses according to the months post symptom onset (PSO). (C) S1, M and N cellular response normalized by anti-CD3 sfu/106PBMC, according to the months PSO. Data is shown as relative sfu per million PBMC. Horizontal bars and whiskers represent median values and interquartile ranges, respectively. The significance between groups was determined using Mann Whitney test, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

https://doi.org/10.1371/journal.ppat.1010211.s005

(TIF)

S6 Fig. Neutralization capacity triggered after natural infection and BNT162b2 vaccination.

Neutralization capacity represented as normalized infection achieved in the presence of 1/50, 1/500 and 1/2000 (from left to right) serum dilutions for each sample in (A) Mild COVID-19 patients; (B) Moderate COVID-19 patients (Mod); (C) Severe COVID-19 patients (Sev); (D) Moderate or severe COVID-19 recovered patients (Post) and (E) after BNT162b2 vaccination (Vac). Coding sample in natural infection, 1: acute phase; 2: hospitalization; 3: convalescent phase. Coding sample in vaccination, 1’: pre-boost; 2’: 15 days after complete vaccination; 3’: 30 days after complete vaccination. Horizontal bars and whiskers represent median values and interquartile ranges, respectively.

https://doi.org/10.1371/journal.ppat.1010211.s006

(TIF)

S1 Table. Lymphocyte populations in successive phases of the disease according to severity.

https://doi.org/10.1371/journal.ppat.1010211.s007

(DOCX)

S1 Raw Data. Data base including all numerical values that were used to generate graphs and statistics.

https://doi.org/10.1371/journal.ppat.1010211.s008

(XLSX)

Acknowledgments

We would like to thank all patients, nurses, medical colleagues and healthy volunteers who contributed to the study.

References

  1. 1. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 2020;395(10229):1054–62. Epub 2020/03/15. pmid:32171076
  2. 2. Chen N, Zhou M, Dong X, Qu J, Gong F, Han Y, et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet. 2020;395(10223):507–13. Epub 2020/02/03. pmid:32007143
  3. 3. Richardson S, Hirsch JS, Narasimhan M, Crawford JM, McGinn T, Davidson KW, et al. Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the New York City Area. Jama. 2020;323(20):2052–9. Epub 2020/04/23. pmid:32320003
  4. 4. Blanco-Melo D, Nilsson-Payant BE, Liu WC, Uhl S, Hoagland D, Moller R, et al. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell. 2020. Epub 2020/05/18. pmid:32416070
  5. 5. Lucas C, Wong P, Klein J, Castro TBR, Silva J, Sundaram M, et al. Longitudinal analyses reveal immunological misfiring in severe COVID-19. Nature. 2020;584(7821):463–9. Epub 2020/07/28. pmid:32717743
  6. 6. Kuri-Cervantes L, Pampena MB, Meng W, Rosenfeld AM, Ittner CAG, Weisman AR, et al. Comprehensive mapping of immune perturbations associated with severe COVID-19. Science immunology. 2020;5(49). Epub 2020/07/17. pmid:32669287
  7. 7. Hadjadj J, Yatim N, Barnabei L, Corneau A, Boussier J, Smith N, et al. Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients. Science. 2020;369(6504):718–24. Epub 2020/07/15. pmid:32661059
  8. 8. Moore JB, June CH. Cytokine release syndrome in severe COVID-19. Science. 2020;368(6490):473–4. Epub 2020/04/19. pmid:32303591
  9. 9. Peng Y, Mentzer AJ, Liu G, Yao X, Yin Z, Dong D, et al. Broad and strong memory CD4(+) and CD8(+) T cells induced by SARS-CoV-2 in UK convalescent individuals following COVID-19. Nature immunology. 2020;21(11):1336–45. Epub 2020/09/06. pmid:32887977
  10. 10. Zuo J, Dowell AC, Pearce H, Verma K, Long HM, Begum J, et al. Robust SARS-CoV-2-specific T cell immunity is maintained at 6 months following primary infection. Nature immunology. 2021. Epub 2021/03/07.
  11. 11. Dan JM, Mateus J, Kato Y, Hastie KM, Yu ED, Faliti CE, et al. Immunological memory to SARS-CoV-2 assessed for up to 8 months after infection. Science. 2021;371(6529). Epub 2021/01/08. pmid:33408181
  12. 12. Ni L, Ye F, Cheng ML, Feng Y, Deng YQ, Zhao H, et al. Detection of SARS-CoV-2-Specific Humoral and Cellular Immunity in COVID-19 Convalescent Individuals. Immunity. 2020;52(6):971–7 e3. Epub 2020/05/16. pmid:32413330
  13. 13. Grifoni A, Weiskopf D, Ramirez SI, Mateus J, Dan JM, Moderbacher CR, et al. Targets of T Cell Responses to SARS-CoV-2 Coronavirus in Humans with COVID-19 Disease and Unexposed Individuals. Cell. 2020;181(7):1489–501 e15. Epub 2020/05/31. pmid:32473127
  14. 14. Le Bert N, Tan AT, Kunasegaran K, Tham CYL, Hafezi M, Chia A, et al. SARS-CoV-2-specific T cell immunity in cases of COVID-19 and SARS, and uninfected controls. Nature. 2020;584(7821):457–62. Epub 2020/07/16. pmid:32668444
  15. 15. Rydyznski Moderbacher C, Ramirez SI, Dan JM, Grifoni A, Hastie KM, Weiskopf D, et al. Antigen-Specific Adaptive Immunity to SARS-CoV-2 in Acute COVID-19 and Associations with Age and Disease Severity. Cell. 2020;183(4):996–1012 e19. Epub 2020/10/05. pmid:33010815
  16. 16. Tan AT, Linster M, Tan CW, Le Bert N, Chia WN, Kunasegaran K, et al. Early induction of functional SARS-CoV-2-specific T cells associates with rapid viral clearance and mild disease in COVID-19 patients. Cell reports. 2021;34(6):108728. Epub 2021/02/01. pmid:33516277
  17. 17. Zhou R, To KK, Wong YC, Liu L, Zhou B, Li X, et al. Acute SARS-CoV-2 Infection Impairs Dendritic Cell and T Cell Responses. Immunity. 2020;53(4):864–77 e5. Epub 2020/08/14. pmid:32791036
  18. 18. Long QX, Liu BZ, Deng HJ, Wu GC, Deng K, Chen YK, et al. Antibody responses to SARS-CoV-2 in patients with COVID-19. Nature medicine. 2020;26(6):845–8. Epub 2020/05/01. pmid:32350462
  19. 19. Wang Y, Zhang L, Sang L, Ye F, Ruan S, Zhong B, et al. Kinetics of viral load and antibody response in relation to COVID-19 severity. The Journal of clinical investigation. 2020;130(10):5235–44. Epub 2020/07/08. pmid:32634129
  20. 20. Lucas C, Klein J, Sundaram ME, Liu F, Wong P, Silva J, et al. Delayed production of neutralizing antibodies correlates with fatal COVID-19. Nature medicine. 2021. Epub 2021/05/07.
  21. 21. Polack FP, Thomas SJ, Kitchin N, Absalon J, Gurtman A, Lockhart S, et al. Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine. The New England journal of medicine. 2020;383(27):2603–15. Epub 2020/12/11. pmid:33301246
  22. 22. Sahin U, Muik A, Derhovanessian E, Vogler I, Kranz LM, Vormehr M, et al. COVID-19 vaccine BNT162b1 elicits human antibody and TH1 T cell responses. Nature. 2020;586(7830):594–9. Epub 2020/10/01. pmid:32998157
  23. 23. Benenson S, Oster Y, Cohen MJ, Nir-Paz R. BNT162b2 mRNA Covid-19 Vaccine Effectiveness among Health Care Workers. The New England journal of medicine. 2021. Epub 2021/03/24. pmid:33755373
  24. 24. Thompson MG, Burgess JL, Naleway AL, Tyner HL, Yoon SK, Meece J, et al. Interim Estimates of Vaccine Effectiveness of BNT162b2 and mRNA-1273 COVID-19 Vaccines in Preventing SARS-CoV-2 Infection Among Health Care Personnel, First Responders, and Other Essential and Frontline Workers—Eight U.S. Locations, December 2020-March 2021. MMWR Morbidity and mortality weekly report. 2021;70(13):495–500. Epub 2021/04/02. pmid:33793460
  25. 25. Geers D, Shamier MC, Bogers S, den Hartog G, Gommers L, Nieuwkoop NN, et al. SARS-CoV-2 variants of concern partially escape humoral but not T-cell responses in COVID-19 convalescent donors and vaccinees. Science immunology. 2021;6(59). Epub 2021/05/27.
  26. 26. Goel RR, Apostolidis SA, Painter MM, Mathew D, Pattekar A, Kuthuru O, et al. Distinct antibody and memory B cell responses in SARS-CoV-2 naive and recovered individuals following mRNA vaccination. Science immunology. 2021;6(58). Epub 2021/04/17.
  27. 27. Mazzoni A, Di Lauria N, Maggi L, Salvati L, Vanni A, Capone M, et al. First-dose mRNA vaccination is sufficient to reactivate immunological memory to SARS-CoV-2 in recovered COVID-19 subjects. The Journal of clinical investigation. 2021. Epub 2021/05/04.
  28. 28. Levin EG, Lustig Y, Cohen C, Fluss R, Indenbaum V, Amit S, et al. Waning Immune Humoral Response to BNT162b2 Covid-19 Vaccine over 6 Months. The New England journal of medicine. 2021. Epub 2021/10/07. pmid:34614326
  29. 29. Ciabattini A, Pastore G, Fiorino F, Polvere J, Lucchesi S, Pettini E, et al. Evidence of SARS-CoV-2-Specific Memory B Cells Six Months After Vaccination With the BNT162b2 mRNA Vaccine. Frontiers in immunology. 2021;12:740708. Epub 2021/10/16. pmid:34650563
  30. 30. (WHO) WHO. Novel Coronavirus: COVID-19 Therapeutic Trial Synopsis. https://wwwwhoint/blueprint/priority-diseases/key-action/COVID-19_Treatment_Trial_Design_Master_Protocol_synopsis_Final_18022020pdf. 2020.
  31. 31. Horndler L, Delgado P, Abia D, Balabanov I, Martinez-Fleta P, Cornish G, et al. Flow cytometry multiplexed method for the detection of neutralizing human antibodies to the native SARS-CoV-2 spike protein. EMBO molecular medicine. 2021;13(3):e13549. Epub 2021/01/21. pmid:33471406
  32. 32. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395(10223):497–506. Epub 2020/01/28. pmid:31986264
  33. 33. Laguna-Goya R, Utrero-Rico A, Talayero P, Lasa-Lazaro M, Ramirez-Fernandez A, Naranjo L, et al. IL-6-based mortality risk model for hospitalized patients with COVID-19. The Journal of allergy and clinical immunology. 2020;146(4):799–807 e9. Epub 2020/07/28. pmid:32710975
  34. 34. Chaplin DD. Overview of the immune response. The Journal of allergy and clinical immunology. 2010;125(2 Suppl 2):S3–23. Epub 2010/03/05. pmid:20176265
  35. 35. Bonifacius A, Tischer-Zimmermann S, Dragon AC, Gussarow D, Vogel A, Krettek U, et al. COVID-19 immune signatures reveal stable antiviral T cell function despite declining humoral responses. Immunity. 2021;54(2):340–54 e6. Epub 2021/02/11. pmid:33567252
  36. 36. McNab F, Mayer-Barber K, Sher A, Wack A, O’Garra A. Type I interferons in infectious disease. Nature reviews Immunology. 2015;15(2):87–103. Epub 2015/01/24. pmid:25614319
  37. 37. Zhang Q, Bastard P, Liu Z, Le Pen J, Moncada-Velez M, Chen J, et al. Inborn errors of type I IFN immunity in patients with life-threatening COVID-19. Science. 2020;370(6515). Epub 2020/09/26.
  38. 38. Bastard P, Rosen LB, Zhang Q, Michailidis E, Hoffmann HH, Zhang Y, et al. Autoantibodies against type I IFNs in patients with life-threatening COVID-19. Science. 2020;370(6515). Epub 2020/09/26. pmid:32972996
  39. 39. Sariol A, Perlman S. Lessons for COVID-19 Immunity from Other Coronavirus Infections. Immunity. 2020;53(2):248–63. Epub 2020/07/28. pmid:32717182
  40. 40. Yu J, Tostanoski LH, Peter L, Mercado NB, McMahan K, Mahrokhian SH, et al. DNA vaccine protection against SARS-CoV-2 in rhesus macaques. Science. 2020;369(6505):806–11. Epub 2020/05/22. pmid:32434945
  41. 41. Deng W, Bao L, Liu J, Xiao C, Xue J, Lv Q, et al. Primary exposure to SARS-CoV-2 protects against reinfection in rhesus macaques. Science. 2020;369(6505):818–23. Epub 2020/07/04. pmid:32616673
  42. 42. Van Praet JT, Vandecasteele S, De Roo A, De Vriese AS, Reynders M. Humoral and cellular immunogenicity of the BNT162b2 mRNA Covid-19 Vaccine in nursing home residents. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2021. Epub 2021/04/08. pmid:33825869
  43. 43. Salmeron Rios S, Mas Romero M, Cortes Zamora EB, Tabernero Sahuquillo MT, Romero Rizos L, Sanchez-Jurado PM, et al. Immunogenicity of the BNT162b2 vaccine in frail or disabled nursing home residents: COVID-A study. Journal of the American Geriatrics Society. 2021. Epub 2021/03/27.
  44. 44. Rossman H, Shilo S, Meir T, Gorfine M, Shalit U, Segal E. COVID-19 dynamics after a national immunization program in Israel. Nature medicine. 2021. Epub 2021/04/21. pmid:33875890
  45. 45. Tan AT, Lim JM, Le Bert N, Kunasegaran K, Chia A, Qui MD, et al. Rapid measurement of SARS-CoV-2 spike T cells in whole blood from vaccinated and naturally infected individuals. The Journal of clinical investigation. 2021;131(17). Epub 2021/10/09.
  46. 46. Bayart JL, Douxfils J, Gillot C, David C, Mullier F, Elsen M, et al. Waning of IgG, Total and Neutralizing Antibodies 6 Months Post-Vaccination with BNT162b2 in Healthcare Workers. Vaccines. 2021;9(10). Epub 2021/10/27. pmid:34696200
  47. 47. Terpos E, Karalis V, Ntanasis-Stathopoulos I, Gavriatopoulou M, Gumeni S, Malandrakis P, et al. Robust Neutralizing Antibody Responses 6 Months Post Vaccination with BNT162b2: A Prospective Study in 308 Healthy Individuals. Life (Basel). 2021;11(10). Epub 2021/10/24. pmid:34685448
  48. 48. Shrotri M, Navaratnam AMD, Nguyen V, Byrne T, Geismar C, Fragaszy E, et al. Spike-antibody waning after second dose of BNT162b2 or ChAdOx1. Lancet. 2021;398(10298):385–7. Epub 2021/07/19. pmid:34274038
  49. 49. Tober-Lau P, Schwarz T, Vanshylla K, Hillus D, Gruell H, Suttorp N, et al. Long-term immunogenicity of BNT162b2 vaccination in older people and younger health-care workers. The Lancet Respiratory medicine. 2021;9(11):e104–e5. Epub 2021/10/24. pmid:34687656
  50. 50. Jeffery-Smith A, Burton AR, Lens S, Rees-Spear C, Davies J, Patel M, et al. SARS-CoV-2-specific memory B cells can persist in the elderly who have lost detectable neutralising antibodies. The Journal of clinical investigation. 2021. Epub 2021/11/30.
  51. 51. Sekine T, Perez-Potti A, Rivera-Ballesteros O, Stralin K, Gorin JB, Olsson A, et al. Robust T Cell Immunity in Convalescent Individuals with Asymptomatic or Mild COVID-19. Cell. 2020;183(1):158–68 e14. Epub 2020/09/28. pmid:32979941
  52. 52. Fernandez-Ruiz M, Olea B, Almendro-Vazquez P, Gimenez E, Marcacuzco A, San Juan R, et al. T cell-mediated response to SARS-CoV-2 in liver transplant recipients with prior COVID-19. American journal of transplantation: official journal of the American Society of Transplantation and the American Society of Transplant Surgeons. 2021. Epub 2021/06/07.
  53. 53. Corman VM, Landt O, Kaiser M, Molenkamp R, Meijer A, Chu DK, et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro surveillance: bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin. 2020;25(3). Epub 2020/01/30.
  54. 54. Cordes AK, Heim A. Rapid random access detection of the novel SARS-coronavirus-2 (SARS-CoV-2, previously 2019-nCoV) using an open access protocol for the Panther Fusion. Journal of clinical virology: the official publication of the Pan American Society for Clinical Virology. 2020;125:104305. Epub 2020/03/07. pmid:32143123
  55. 55. Molina M, Allende LM, Ramos LE, Gutierrez E, Pleguezuelo DE, Hernandez ER, et al. CD19(+) B-Cells, a New Biomarker of Mortality in Hemodialysis Patients. Frontiers in immunology. 2018;9:1221. Epub 2018/07/03. pmid:29963040