Demographic and clinical characteristics of the relapsing-remitting MS patients and HCs
29 MS patients (15 relapsing patients, 14 remitting patients) and 13 HCs were enrolled in this study, as shown in Table 1. No statistical differences in the age and female ratio among groups. There were also no significant differences in the onset age and disease duration between relapsing and remitting MS patients. The EDSS score of the relapsing MS patients was slightly higher than the remitting MS patients, but there was no significant difference (P > 0.05).
Table 1
Demographic and clinical characteristics of the relapsing-remitting MS patients and HCs
Variable
|
Relapsing MS (n=15)
|
Remitting MS (n=14)
|
HCs (n=13)
|
P value
|
Age (years)
|
37 [26-51]
|
35.5 [17-52]
|
29 [16-48]
|
0.309
|
Female, no. (%)
|
11 (73.3)
|
10 (71.4)
|
9 (69.2)
|
1.000
|
Onset age (years)
|
34 [22-46]
|
29.5 [16-38]
|
NA
|
0.196
|
Disease duration (years)
|
2.3 [0.2-16]
|
3.2 [0.3-16]
|
NA
|
0.776
|
EDSS
|
3.0 [1.0-9.0]
|
2.5 [0-4]
|
NA
|
0.235
|
MS, multiple sclerosis; HCs, healthy controls; EDSS, Kurtzke’s Expanded Disability Status Scale; NA, not applicable.
Data are presented as number (percentage) or median [range]. Differences among relapsing MS, remitting MS and HCs were analyzed by Kruskal-Wallis test (age), Chi-square test (female ratio), Mann-Whitney U test (onset age, disease duration, and EDSS). P value < 0.05 was considered statistically significant.
|
Increased circulating CD8 + CXCR5 + T cells positively correlated with new gadolinium enhancements in the relapsing MS patients
In this study, blood samples from MS patients and HCs were collected and peripheral blood mononuclear cells were isolated timely. According to the gating strategy as previously reported (Fig. 1A) [11, 18], the CD8+CXCR5+ T cells, CD8+CXCR5− T cells, CD8+CXCR5+PD-1+ T cells, CD8+CXCR5−PD-1+ T cells and Tfh cells (CD4+CXCR5+PD-1+) were evaluated by flow cytometry. The percentages of CD8+CXCR5+ T, CD8+CXCR5− T and Tfh cells were significantly higher in the relapsing MS patients (P < 0.05). The numbers of the three subsets showed the same trends. While in the remitting phase, CD8+CXCR5+ T cell frequency and cell number did not decrease as sensitively as CD8+CXCR5− T (P < 0.001 for number) or Tfh cells (P < 0.05 for frequency). Like Tfh cells, both CD8+CXCR5+PD-1+ T and CD8+CXCR5−PD-1+ T cells had similar dynamic changes as the condition fluctuated despite no statistical difference (Fig. 1B). In addition, CD8+CXCR5+ T and CD8+CXCR5+PD-1+ T cells showed very similar changes both in the proportion and cell number during relapsing and remitting phases. Further correlation analysis also showed great positive correlations between CD8+CXCR5+ T and CD8+CXCR5+PD-1+ T cells (“r = 0.8693, P < 0.0001” for frequency, “r = 0.9153, P < 0.0001” for cell number) (Fig. 2A, Figure S1A). What’s more, significant correlations of both CD8+CXCR5+ T cell frequency with CD8+CXCR5− T cell frequency (r = 0.6553, P = 0.0001) and CD8+CXCR5+ T cell number with CD8+CXCR5−PD-1+ T cell number (r = 0.5607, P = 0.0066) were observed. While, there were no significant differences in the proportion and cell number of CD8+CXCR5+ T cells with Tfh cells. In terms of disease severity, CD8+CXCR5+ T cells showed no remarkable correlations with EDSS score (“r = 0.1826, P = 0.3431” for frequency, “r = 0.2365, P = 0.2350” for cell number) (Fig. 2B, Figure S1B). But unexpectedly, there were significant correlations between CD8+CXCR5+ T cells and central nervous system (CNS) new gadolinium enhancements of relapsing patients (“r = 0.6878, P = 0.0373” for frequency, “r = 0.7504, P = 0.0189” for cell number) (Fig. 2C, Figure S1C). The results suggest that increased circulating CD8+CXCR5+ T cells in relapsing MS patients might be involved in disease progress and indicate the newly developed demyelinations in CNS.
CD8 + CXCR5 + T cells localized into germinal centers and spinal cords which correlated with clinical score of EAE mice
To discover the role of CD8+CXCR5+ T cells in EAE mice, its dynamic changes in frequency and morphology in secondary lymphoid organs (SLOs) and central nervous system (CNS) were examined. First, the EAE model was induced in C57BL/6J mice and the clinical score was divided into four phases according to the days post immunization (dpi): pre-EAE phase (5–7 dpi, score 0), peak-EAE phase (15–21 dpi, score 4-4.5), remission phase (22–28 dpi, score 1.5–2.5), and chronic phase (35–45 dpi, score 2–3) (Fig. 3A). Then the dynamic changes of CD8+CXCR5+ T cells in spleens, inguinal lymph nodes (LN) and CNS of EAE mice at different phases were detected by flow cytometry. CD8+CXCR5+ T cell frequencies all increased notably at the peak phase and declined during the remission phase. In the chronic phase, the proportion increased again (Fig. 3B-E). Such changes in the CNS and the spleens were very similar. Meanwhile, CD8+CXCR5+ T cells in LN seemed to increase more sensitively before the onset. PD-1-expressing CD8+CXCR5+ T cells in the CNS were also analyzed. The subtype changes were in conformity with CD8+CXCR5+ T cells and a good correlation was observed (Figure S2A, B).
The germinal center is critical for B cell maturation and high-affinity antibody production. Previous studies pointed Tfh cell which is a powerful helper for B cell activation and differentiation localized at the T-B border or resided in the light zone [11, 19]. The present results showed CD8+CXCR5+ T cells localized into the germinal center and expanded within the B cell follicles of EAE mice (Fig. 4A). Additionally, most CD8 T cells were located alongside CD4 T cells. The data suggest CD8+CXCR5+ T cells could also be involved in B cell regulation.
In line with early reports concerning about CXCL13 in the CNS of EAE mice [15, 20], we also found dramatically upregulated CXCL13 and CXCR5 in the CNS (Figure S3B). CXCR5 positive cells were recruited to areas with high levels of CXCL13 (Figure S3A). Further multiplex immunofluorescence staining revealed aggregations of CD8+CXCR5+ T cells in the spinal cords of EAE mice during the peak phase, particularly at sites where immune cells had infiltrated (Fig. 4B). CD8 T cells and CD4 T cells were found in close proximity in the spinal cord. And a small portion of CD8+CXCR5+ T cells were observed near CD20-positive cells. Meanwhile, flow cytometry analysis showed a positive correlation between the frequency of CD8+CXCR5+ T cells in the CNS and Tfh cells, germinal center (GC) B cells, as well as clinical scores of EAE mice respectively (“r = 0.8218, P < 0.0001” for Tfh cell, “r = 0.5128, P = 0.0147” for GC B cell, “r = 0.6820, P = 0.0481” for EAE score) (Fig. 4C-E). These findings suggest CD8+CXCR5+ T cells may contribute to disease progression in EAE by proliferating in follicles or ectopic lymphoid structures, and participating in humoral immunity.
Tfh-like helper functions, but not enhanced cytotoxicity of CD8+CXCR5+ T cells in EAE
Previous studies in cancers, chronic infections and autoimmune diseases have indicated diverse properties of CD8+CXCR5+ T cells including cytotoxic, regulatory-like and B-cell helper functions [21–23]. But the features of this CD8 subset in EAE remain elusive. To clarify this, RNA sequencing was performed on sorted splenic CD8+CXCR5+ T cells and CD8+CXCR5− T cells from EAE mice in peak phase and contemporaneous control mice. The gene expression pattern of CD8+CXCR5+ T cells in both the EAE and control group was distinct from that of CD8+CXCR5− T cells. Pathway analysis of differential upregulated gene profiles revealed that in EAE mice, the differential gene expression was mainly associated with rheumatoid arthritis, B cell receptor signaling pathway, intestinal immune network for IgA production and NF-kappa B signaling pathway (Fig. 5B). While in the control group, the differential gene expression was mainly associated with NF-kappa B signaling pathway, Fc gamma R-mediated phagocytosis, PI3K-Akt signaling pathway and Calcium signaling pathway (Fig. 5A). The differences between CD8+CXCR5+ T cells in EAE mice and the control group were further analyzed. Pathway analysis of differential upregulated gene profiles showed genes in CD8+CXCR5+ T cells of EAE were associated with rheumatoid arthritis, NF-kappa B signaling pathway and B cell receptor signaling pathway (Fig. 5C).
The most studied subset of CD8 T cells is the cytotoxicity T lymphocyte with the ability to kill target cells by perforin, granules and interferon-γ (IFN-γ) [24]. While under different conditions, CD8 T cells may acquire diverse functions. RNA sequencing revealed the functions associated with rheumatoid arthritis, B cell receptor signaling pathway, intestinal immune network for IgA production and NF-kappa B signaling pathway of CD8+CXCR5+ T cells in EAE, which implicated that CD8+CXCR5+ T cells may be involved in inflammatory process and humoral immune response. Further analysis of differential genes related to inflammation and humoral immunity showed CD8+CXCR5+ T cells from control mice expressed higher levels of cytotoxic granzyme A, granzyme B and inflammatory mediators like S100a9 and matrix metalloproteinases (MMPs) (Fig. 6A). In addition, this subset shared gene profile with Tfh cells such as higher levels of transcriptional factor Bcl6, follicular markers Pdcd1, Cd200, Cd40lg and lower levels of Ccr7. We also observed higher expression of pro-inflammatory cytokines Il4, Il15, Il18 and chemokines Ccl3, Ccl4, Cxcl9 and Cxcl10 (Fig. 6A). Interestingly, CD8+CXCR5+ T cells in control mice expressed higher Il10 than CD8+CXCR5− T cells. All the above results suggested stronger cytotoxicity, immunomodulation and Tfh-like functions of CD8+CXCR5+ T cells than CD8+CXCR5− T cells in the control group.
Partially consistent with the control mice, CD8+CXCR5+ T cells from EAE mice also expressed higher levels of Bcl6, Cd200 and lower levels of Eomes, Il7r and Ccr7 like Tfh cells. The expression of Il4, Il18 and chemokines Ccl3, Cxcl9 and Cxcl10 were also higher than CD8+CXCR5− T cells (Fig. 6B). Unexpectedly, CD8+CXCR5+ T cells from EAE mice expressed higher CD177, which is important for neutrophil activation and transendothelial migration, than CD8+CXCR5− T cells markedly (Fig. 6B). While, the expression of granzyme A and granzyme B were not different significantly. And the less expressed gene of effector marker Ctsw and Cd226 suggested CD8+CXCR5+ T cells might be less cytotoxic in EAE mice. To further characterize the functions of CD8+CXCR5+ T cells in EAE and control mice, flow cytometry was applied and the phenotype differences between CD8+CXCR5+ and CD8+CXCR5− T cells were observed from control and EAE mice. First, under phorbol 12-myristate 13-acetate (PMA)-stimulated conditions, we observed higher frequencies of effector marker CD107a, perforin and granzyme B in CD8+CXCR5+ T cells compared with CD8+CXCR5− T cells from control mice (Fig. 6C, D). While, CD8+CXCR5+ T cells in EAE mice expressed less perforin than CD8+CXCR5− T cells significantly, albeit differences of CD107a and granzyme B expressions between CD8+CXCR5+ T cells and CD8+CXCR5− T cells were similar in EAE mice compared with control. Co-culture mouse oligodendrocyte precursor cells with CD8+CXCR5+ T or CD8+CXCR5− T cells from EAE mice in vitro showed lysis rates from both groups increased as the effect-target ratio increased. While CD8+CXCR5+ T cells did not show enhanced cytotoxicity than CD8+CXCR5− T cells in EAE mice (Fig. 6E). The results suggested CD8+CXCR5+ T cells could be less effective in EAE conditions.
Second, the frequencies of ICOS and CD40L were higher in CD8+CXCR5+ T cells than in CD8+CXCR5− T cells. Additionally, cytokines including IL-4, IL-6 and IL-21 were significantly increased in CD8+CXCR5+ T cells (Fig. 6C, D). IL-4 and CD40L were much higher in EAE mice. Meanwhile, sorted splenic CD8+CXCR5+ T cells from control or EAE mice were co-cultured with B cells from either control or EAE mice. Under MOG35 − 55 conditions, the level of anti-MOG35 − 55 autoantibody produced in the supernatant was measured by CLISA. The results showed that the levels of anti-MOG35 − 55 antibody were higher when B cells were co-cultured with CD8+CXCR5+ T cells from both the control and EAE groups, compared to B cells alone. This effect was particularly pronounced in the EAE group (Fig. 6F), suggesting CD8+CXCR5+ T cells from EAE mice could help antigen-sensitized B cells produce much more autoantibody through Tfh-like helper functions.