Lymphopenia and altered CD4+T cell profile
From December 2022 to February 2023, 57 patients (mean age of 73 years) diagnosed with SARS-CoV-2 infection were admitted to the Physical and Rehabilitation Department of Tianjin Medical University General Hospital. The severity of acute respiratory distress syndrome (ARDS) in the elderly was age-dependent, and hypoxia indices SpO2% and PaO2/FiO2 correlated negatively with age, whereas the acute physiologic assessment and chronic health evaluation (APACHE) score and assisted ventilation grade correlated positively with age (Fig. 1A-E). Overall, eight patients died (described as deceased patients, DPs) at 5–23 days post symptom onset (PSO). The DPs had gone through severe hypoxemia (mean SpO2 of 81% and PaO2/FiO2 of 89 mm Hg), and five out of eight DPs underwent invasive ventilation. The cohort details are summarized in Table 1.
Table 1
Demographics and clinical characteristics of all donors.
Characteristic | MPs (n = 22) | SPs (n = 27) | DPs (n = 8) | CDs (n = 27) | P value |
Demographic | | | | | |
Age, median (SD), y | 69 (8.5) | 76 (9.0) | 78 (7.9) | 73 (5.6) | .0404* .0343# |
Sex, No. (%) | | | | | |
Male | 16 (72.7) | 15 (55.6) | 6 (75.0) | 16 (59.3) | |
Female | 6 (27.3) | 12 (44.4) | 2 (25.0) | 11 (40.7) |
Comorbidity, No. (%) | | | | | |
Cardiovascular diseases | 15 (68.2) | 19 (70.4) | 3 (37.5) | 22 (81.5) | |
Endocrine system disease | 10 (45.5) | 7 (25.9) | 4 (50.0) | 12 (44.4) |
Nervous system disease | 8 (36.4) | 8 (29.6) | 3 (37.5) | 13 (48.1) |
Other chronic comorbidity | 3 (13.6) | 5 (18.5) | 4 (50.0) | 3 (11.1) |
Oxygen saturation, median (SD), % | 95 (0.4) | 88 (6.9) | 81 (5.0) | NA | < .0001* < .0001# .0048& |
PaO2/FiO2, median (SD), mmHg | 331 (58.9) | 159 (56.2) | 89 (19.8) | NA | < .0001* < .0001# .0104& |
Maximum oxygen treatment, No. (%) | | | | | |
No oxygen therapy | 3 (13.6) | 0 (0.0) | 0 (0.0) | NA | |
Mask/nasal prongs | 18 (81.8) | 19 (70.4) | 3 (37.5) | NA | |
Noninvasive/high-flow ventilation | 1 (4.5) | 2 (7.4) | 0 (0.0) | NA | |
Invasive ventilation | 0 (0.0) | 6 (22.2) | 5 (62.5) | NA | |
Hospitalized status, No. (%) | | | | | |
ICU | 4 (18.2) | 11 (40.7) | 7 (87.5) | NA | |
Non-ICU | 18 (81.8) | 16 (59.3) | 1 (12.5) | NA | |
Duration of hospitalization | 34 (18.5)a | 54 (44.5)a | 14 (7.1)b | NA | |
APACHE II score | 15 (2.1) | 27 (6.8) | 31 (5.7) | NA | < .0001* < .0001# .0005& |
Positive Nucleic acid qPCR at sampling, No. (%) | 22 (100.0) | 27 (100.0) | 8 (100.0) | 0 (0.0) | |
NA, non-available; MPs: Moderate patients; SPs: Severe patients; DPs: Deceased patients; CDs: Convalescent donors; Data were analyzed by One-way-ANOVA with Tukey correction or Kruskal-Wallis test according to the distribution of data. P values < 0.05 were denoted.
* indicate the statistical differences between MPs and SPs;
# indicate the statistical differences between MPs and DPs;
& indicate the statistical differences between SPs and DPs;
a Mean time from symptom onset to date of discharge; b Mean time from symptom onset to date of death;
Elderly patients with acute COVID-19 had CD3+ T deficiency, which was ameliorated after recovery (Fig. 1F). Kinetics of CD4+ T revealed consistency with CD3+ whereas CD8+ didn’t, indicating that altered CD4+ T magnitude, rather than CD8+, contributed to lymphopenia, especially in severe (SPs) and DPs (Fig. 1G-H). SARS-CoV-2-specific CD8+ effector T cells increase rapidly as early as 1 day PSO, accompanied by expression of potent cytotoxic molecules, such as gamma interferon (IFNγ), granzyme B (GZMB), perforin, and CD107a, CD8+ memory T (TM) cells still retain similar expression profiles as effector compartments (15). Accordingly, indiscriminate and sustained secretion of GZMB by CD8+T cells was observed in all individuals, including fully recovered convalescent donors (CDs) (Fig. 1I ), representing sustained cytotoxicity sensitization of CD8+ T cells after COVID-19 in the elderly (16).
We classified memory TM cells into naive (TN: CD45RA+CCR7+), effector memory (TEM: CD45RA–CCR7–), central memory (TCM: CD45RA–CCR7+) and effector memory-expressing CD45RA (TEMRA: CD45RA+CCR7−) subsets (Fig. 1J-K). The DPs exhibited feeble CD4+TEM responses but significantly higher TN proportion than moderate patients (MPs) and SPs, indicating a dysfunctional CD4+ TM profile in DPs during acute COVID-19 (Fig. 1L-M). We considered that if patients’ severity were relatively under control (such as MPs and SPs), COVID-19 contributed to the activation of effector T immunity and promoted the acquirement of SARS-CoV-2-specific central memory imprints. Following recovery, their TM profile changed and resembled the pattern of CDs, which preserved a high proportion of TCM for combating antigen re-encounter in the future (Fig. 1N). For critically ill patients (DPs), however, poor physiological condition and severe hypoxemia led to inferior T cell activation, feeble memory counterparts preservation and poor prognosis.
Moreover, CD8+TEM and CD8+TEMRA dominated CD8+TM as previously reported (17). CD8+TEM also abated in DPs (Supplementary Fig. 1A).
Activation signatures of T cells in geriatric
Spike-reactive activation profile of T cells are heterogeneous. In our study, CD38 and HLA-DR were uniformly co-expressed irrespective of disease severity whereas CD69+CD137+CD4+ were not. All patients, including fully recovered CDs, had elevated CD38+HLA-DR+CD4+ level compared with healthy donors (four healthy donors, HDs, with a mean age of 45 years), whereas CD69+CD137+CD4+ were induced in MPs and SPs but not in DPs or CDs (Fig. 2A-B). Actually, CD69+CD137+CD4+ correlated negatively with age and positively with SPO2 and PaO2/FiO2 (Fig. 2C-E). However, no correlation was observed between CD38+HLA-DR+CD4+ and COVID-19 severity. Further dissection revealed that HLA-DR+CD4+ displayed slightly discordant trend with CD38+CD4+, the percentage of HLA-DR+CD4+ increased mildly whereas CD38+CD4+ decreased as the disease progressed (Fig. 2F-G). Spearman correlation analyses showed that SARS-CoV-2-specific CD38+CD4+ correlated positively with SpO2% and PaO2/FiO2, and negatively with the APACHE and WHO Ordinal scale scores, whereas HLA-DR+CD4+ exhibited opposite trends with weaker degree of correlation (Fig. 2H-K).
Regarding the corresponding CD8 compartments, CD38+HLA-DR+CD8+ remained higher in all patients and CDs compared with HDs (Supplementary Fig. 1B). Statistically, CD69+CD137+CD8+ remained similar among all groups, with certain severe or convalescent individuals raised its expression level (Supplementary Fig. 1C). Neither of the co-expression clusters correlated with the aforementioned clinical features. Briefly, we deemed that co-expression of CD38 and HLA-DR, regardless of CD4+ or CD8+, were inclined to act as general T cell activation indicator, and CD69+CD137+CD4+ were more likely to be indicators of severity-dependent T cell activation, especially in the scenario of CD4 perturbation.
For SARS-CoV-2-specific cytokine, IL-4 was almost undetectable in all individuals, whereas IFN-γ+CD4+ was induced in all patients regardless of disease severity compared with HDs (Supplementary Fig. 1D). IFN-γ+CD8+ peaked in SPs and DPs (Supplementary Fig. 1E). IFN-γ+CD4+ production was sustained to a certain degree and IFN-γ+CD8+ restored to baseline level in CDs, indicating that recovered individuals tend to preserve evidence of virus-specific CD4+ T cell responses than CD8+ (18).
Severity-dependent impairment of Tfh subsets and TFR homeostasis
Tfh cells support high affinity and long-term antibody responses. Though usually located in secondary lymphoid organs, T cells bearing features of Tfh cells can also be identified in peripheral blood (circulating Tfh, cTfh) (19). He et al elucidated that within human circulating CXCR5+CD4+ T cells, CCR7loPD-1hi subset had a partial effector phenotype (Tfh-em), whereas CCR7hi PD-1lo cells displayed central-memory like phenotype (Tfh-cm) (20). Upon antigen re-exposure in vivo, Tfh-em take effect in early phase with a transient increase around 7 days, while Tfh-cm subset normally peak at around 21 days and may persist into memory phase (21). In our study, as the disease progressed, total Tfh decreased, which is similar to the trend observed for CD4+T cells. HDs and CDs retained baseline level of it (Fig. 3A). In detail, Tfh-cm outnumbered Tfh-em regardless of the state of illness, and showed complementary trend with Tfh-em, signifying a “conversion” (20) but “quantitatively conserved” relationship between the two clusters (Fig. 3B). The “quantitative gap” between Tfh-cm and Tfh-em represented “conversion ratio”. HDs and CDs exhibited low conversation ratio due to inactivated or recovered Tfh responses, respectively (Fig. 3C). For patients, MPs and SPs had relative higher Tfh-em: Tfh-cm ratio (Supplementary Fig. 2A), whereas DPs induced minimal Tfh-cm-Tfh-em conversion, indicating inferior Tfh-cm mobilization and Tfh-em activation (Fig. 3C). Percentages of Tfh-cm and Tfh-em in CD4 also wanned as the disease progressed (Fig. 3D-E).
We also classified Tfh into Tfh1 (CXCR3+CCR6−), Tfh2 (CXCR3−CCR6−) and Tfh17 (CXCR3−CCR6+) subsets. All three subsets are supportive of antibody production with Tfh17 and Tfh2 being superior contributor, owing to their naive B cells helper function, whereas Tfh1 cells merely support memory B cells but not naive B cells (22). Within CXCR5+ populatuion, the enrichment of Th2 phenotype was the most significant (Supplementary Fig. 2B). Both Tfh1 and Tfh2 revealed downtrend as the COVID-19 aggravated, MPs exhibited significantly higher enrichment level of these two subsets compared with DPs and CDs. Tfh17 declined to practically undetectable level in DPs with no statistical differences among DPs, MPs and SPs, but MPs and SPs showed advantages over CDs (Fig. 3F-G). Generally, the most and least enriched Tfh subsets were observed in MPs and DPs, respectively. TFR cells are thought to share phenotypic characteristics with Tfh and conventional Foxp3+ regulatory T cells (Treg) yet are distinct from either. The precise role of TFR in regulating germinal center (GC) and Ab production remains controversial. Similar to the Tfh subsets, TFR percentages in CD45RA−or CD4+ population both wanned as the COVID-19 progressed, with DPs merely triggered baseline level of TFR (Fig. 3H).Trivial changes were found in other subsets such as Tregs (Fig. 3I). Briefly, the elderly underwent severity-dependent defects in Tfh subsets and TFR during acute COVID-19.
Spearman correlation exhibited that Tfh-cm was negatively correlated with age and WHO ordinal scale score, while correlated positively with SpO2%. Nevertheless, Tfh-em didn’t show correlation with any of the aforementioned indicators. We speculated that Tfh-em is characterized by rapid maturation in the context of antigen exposure, while Tfh-cm is more like a resting but relative sustained status (20), which may comprehensively represent patients’ status. Among Tfh1, Tfh2 and Tfh17, only Tfh2 correlated with age, SpO2% and PaO2/FiO2 (Fig. 3J).
Ab responses and their correlation with TfH subsets
Antibodies (Abs) are key components of effective immune responses against viruses. Specificity or affinity of antibody responses can be fine-tuned via somatic hypermutation and affinity-based positive selection in germinal centers (23). In our study, the magnitude of IgG and IgG subclass profile were similar among all groups, with ancestral S1-specific binding antibodies (BAbs) being slightly lower than corresponding BA.5 BAbs. IgG1 and IgG3 dominated IgG while IgG2 and IgG4 were relative rare, which were consistent with previous reports (Fig. 4A-B) (24, 25).
For neutralizing Abs (NAbs), BA.5 and ancestral-NAbs were robust whereas CH.1.1-NAbs were relative feeble (Fig. 4C). Moreover, IgG was correlated with the corresponding NAbs (Fig. 4D). These results support the fact that BA.5 variant dominated the epidemic wave at the end of 2023 in Tianjin, China (Supplementary Fig. 3). COVID-19 severity influenced the magnitude of NAb titers, which was robust in MPs whereas wanned in SPs and DPs.
Spearman correlation analyses manifested that Tfh-em correlated positively but weakly with ancestral and BA.5-NAbs, while Tfh-cm correlated moderately with the two NAbs (Fig. 4E). These phenomena may be attributed to the “temporal match” between Tfh-cm kinetics and NAb profiles. Serial blood collection from nine of 57 elderly patients (five MPs and four SPs) was performed at < 7, 7–14 and 14–28 day PSO. In this small cohort, we observed that NAbs production was time-dependent, as they surged after 7 day and peaked at approximately 14 day PSO (Fig. 4F). Unlike Tfh-em, which differentiates in the early phase after antigen exposure (< 7 day PSO), Tfh-cm represents the majority of CXCR5+ T cells under resting conditions, which peaked at approximately 7–14 day PSO in our study (Fig. 4G). When subdividing all individuals into seven time phases, NAb peaked at approximately 12–15 day PSO and maintained at high level afterwards (Fig. 4H). Tfh-cm increased gradually after COVID-19 onset and peaked at around 8–11 day PSO, conforming to the temporal fluctuation of NAbs to a certain extent (Fig. 4I). Compared with Tfh-em, NAbs production pattern is slightly hysteretic but sustained, similar to that of Tfh-cm.
Tfh1, Tfh2, Tfh17 and TFR all correlated positively with ancestral and BA.5-NAb but not CH.1.1 (Fig. 4E). Recent study unveiled that TFR curtailed cytotoxic TFH responses, which express high levels of granzyme B akin to CD8+ cytotoxic T cell and have the potential of killing B cells in the GC (26). In our study, we did not detect cytotoxic Tfh subsets, however, we found that as the severity of COVID-19 increased, TFR decline was accompanied by a slight increase of CD4+GZMB, which correlated negatively with NAb titers (Fig. 4E). These results may be indicative of the relationship between TFR and NAb, that is, TFR deficiency underlay the insufficient restraint of cytotoxic Tfh, which consequently dampened GC B response and led to NAb generation impediment. However, the conclusion has not been confirmed.
Vaccination contributed to Tfh subsets recovery and NAbs elevation
We divided 55 out of 57 patients with available vaccination records into age- and sex-matched groups according to SARS-CoV-2 vaccination history. The cohort details are summarized in Table 2. Increased CD4+ rather than CD8+ T may have facilitated the minor elevation of CD3+ T among vaccinees (Fig. 5A). We wanted to investigate which CD4+ subsets are responsible for the raise of CD4+T. No differences of activated CD4+T cells or TM subpopulations were found between non-vaccinated individuals and vaccinees (Fig. 5B-C). Further investigation revealed that total Tfh, Tfh-cm and Tfh-em, increased synchronously in vaccinated individuals (Fig. 5D). Percentage of Tfh2 raised significantly, Tfh1 and Tfh17 also presented upward trend although the differences were not statistically significant (Fig. 5E). TFR percentage of CD45RA− was elevated, accompanied by a mild reduction of GZMB+CD4+ cells (Fig. 5F). Treg frequency was not affected by vaccination (Fig. 5G).
BAbs remained similar between groups whereas ancestral and BA.5-NAbs were significantly higher in the vaccinated individuals (Fig. 5H-I).We also divided CDs into non-vaccinated and vaccinee cohorts, and found that BAbs remained unchanged, whereas ancestral and BA.5-NAbs were significantly boosted by vaccination (Fig. 5J-K).
Taken together, vaccination attenuated Th imbalance by increasing percentages of Tfh-em, Tfh-cm, Tfh2, TFR, and augmented the magnitude of ancestral and BA.5 NAbs. NAbs, CD4+ T cells, and CD8+ T cells all have protective roles against SARS-CoV-2 infection, high degree of synergy and redundancy between the three branches of adaptive immunity bring about diverse pathways to robust immunity (15). Hence, albeit the direct relationship between vaccination history and illness status still needs further investigation, it is worth mentioning that vaccination restored the impairment of Tfh homeostasis and promoted NAbs production, which potentially improve the prognosis among the elderly.
Table 2
Demographics and clinical characteristics of non-vaccinated individuals and Vaccinees
Characteristic | Non-vaccinated (n = 15) | Vaccinated (n = 40) | P value |
Demographic | | | |
Age, median (SD), y | 78 (9.8) | 73 (8.3) | NS |
Sex, No. (%) | | | |
Male | 8 (53.3) | 25 (69.4) | NS |
Female | 7 (46.7) | 11 (30.6) |
Comorbidity, No. (%) | | | |
Cardiovascular diseases | 11 (73.3) | 23 (63.9) | NS |
Endocrine system disease | 7 (46.7) | 10 (27.8) |
Nervous system disease | 5 (33.3) | 11 (30.6) |
Other chronic comorbidity | 1 (6.7) | 9 (25.0) |
Oxygen saturation, median (SD), % | 90 (5.8) | 89 (7.7) | NS |
PaO2/FiO2, median (SD), mmHg | 190 (89.0) | 213 (109.0) | NS |
Maximum oxygen treatment, No. (%) | | | |
No oxygen therapy | 0 (0.0) | 1 (2.8) | NS |
Mask / nasal prongs | 12 (80.0) | 23 (63.9) |
Noninvasive/ high-flow ventilation | 1 (6.7) | 3 (8.3) |
Invasive ventilation | 2 (13.3) | 9 (25.0) |
Hospitalized status, No. (%) | | | |
ICU | 5 (33.3) | 16 (44.4) | NS |
Non-ICU | 10 (66.7) | 20 (55.6) |
Duration of hospitalization | 48 (58.4) | 44 (26.2) | NS |
APACHE II scorea | 25 (6.8) | 28 (7.6) | NS |
Positive Nucleic acid qPCR at sampling, No. (%) | 15 (100.0) | 36 (100.0) | NS |
NS, no significance; |
Quantitative parameters were analyzed using t test or non-parametric Mann-Whitney tests according to data distribution; categorical variables were analyzed using Chi-square test or Fisher's exact test. |