Lymphocyte profiles after a first demyelinating event suggestive of multiple sclerosis reveal early monocyte and B cell alterations

2.1 Patients

Twenty treatment-naïve pwCIS (without MS disease-modifying treatment) from the Berlin CIS cohort (NCT01371071), an ongoing prospective observational study initiated at Charité – Universitätsmedizin Berlin in 2011, were selected to assess lymphocyte profiles and B cell subsets by using multicolour flow cytometry. The inclusion criteria, as described previously¹⁴, were: >18 years of age and a first clinical event suggestive of demyelination of the CNS, according to the McDonald 2010 criteria ¹⁵, within the 6 first months before cohort enrolment. Exclusion criteria were inability or unwillingness to provide informed consent, incomplete collection of blood samples for follow-up, and a history of any health condition or devices impeding immunological, clinical, and radiological examination. Moreover, sample selection was also limited to those individuals with sufficient availability of frozen peripheral blood mononuclear cells (PBMCs). All participants provided written informed consent and ethics approval for this study was given by the institutional review board of Charité – Universitätsmedizin Berlin (EA1/180/10).

During the visits, all patients underwent clinical and magnetic resonance imaging (MRI) assessments and blood extraction for PBMCs isolation. The first visit (baseline visit) was scheduled between 3-6 months after a first clinical event suggestive of demyelination of the CNS according to the McDonald 2010 criteria, and before administration of GLAT. The second visit was performed at around 12 months after the first clinical event, and the third visit at around 24 months after the first clinical event. Out of twenty pwCIS, ten received GLAT after baseline visit (due to signs of disease activity), whilst the other half remained untreated. PBMCs were isolated using Biocoll separating solution and stored in liquid nitrogen according to standard operating procedures (SOP) at NeuroCure Clinical Research Center, Charité – Universitätsmedizin Berlin. Frozen PBMCs from fifteen unmatched healthy donors (HDs) were used as controls. Disease activity was defined as the presence of new lesions in T2-weighted and/or contrast-enhancing MRI, and/or new clinical events suggestive of demyelination of the CNS, and/or confirmed deterioration on EDSS at 12 and/or 24 months of follow-up.

2.2 Immunophenotyping

Frozen PBMCs samples were thawed and subsequently stained and analysed by flow cytometry. Briefly, a set of markers to discriminate between lymphocytic profiles and B cell subsets was established: CD3, CD4, CD14, CD19, CD20, CD25, CD27, CD38, CD127, CD138, IgD, and HLA-DR. A live/dead fixable blue dead cell stain kit (Invitrogen) was used per the manufacturer’s instructions to identify dead cells. Cells were stained for 15 min at 4°C, washed and filtered before acquisition.

Stained samples were analysed on LSR Fortessa X-20 (LSRFortessa; BD Biosciences, USA). To ensure stable analysing conditions, flow cytometer setup and tracking beads from BD were used daily. Diva and FlowJo software (v10.4; BD, USA) were used to determine different subpopulations. Results are expressed as percentages. For details on antibodies and gating strategy (see Supplementary Figure 1).

• Download figure
• Open in new tab

• Download figure
• Open in new tab

• Download figure
• Open in new tab

Supplement Figure 1: Gating Strategy

2.3 MRI acquisition and post-processing

All study subjects were scanned using 3 Tesla (Tim Trio; Siemens, Erlangen, Germany) MRI scanners at the Berlin Center for Advanced Neuroimaging. Pre- and post-contrast 3D T1-weighted magnetisation prepared rapid gradient-echo (MPRAGE) images (resolution 1 × 1 × 1 mm³; TR = 1,900 ms, TE = 3.03 ms, TI = 900 ms, flip angle 9°) and a 3D T2-weighted fluid -attenuated inversion recovery (FLAIR) images (resolution 1 × 1 × 1 mm³; TR = 6,000 ms, TE = 388 ms, TI = 2,100 ms, flip angle 120°) were acquired for each patient. MPRAGE and FLAIR were cropped, co-registered to MNI-152 space (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Fslutils#Tools), N4-bias corrected (http://stnava.github.io/ANTs/) and linearly co-registered with each other using FSL FLIRT ¹⁶.

Longitudinal co-registration was performed by using the same FSL tools and each MPRAGE and FLAIR from follow-up scans per patient were co-registered to their baseline scan. T2-hyperintense and gadolinium contrast-enhancing lesions were segmented from FLAIR and post-contrast MPRAGE images manually using ITK-SNAP ¹⁷ by 2 expert MRI technicians with more than 10 years of experience with multiple sclerosis lesion identification. Whole brain T2-hyperintense lesion count was extracted using FSL cluster.

2.4 Statistics

Statistical analysis was performed using JAMOVI for macOS (The jamovi project (2020). jamovi (v.2.0.0.0) [Computer Software]. Retrieved from https://www.jamovi.org) and Prism 8.4.2 (GraphPad, La Jolla, USA). An independent t-test was used to compare percentages between two unpaired groups. A Fisher’s exact test was used to analyse any difference between categorical variables. For cross-sectional data, one-way ANOVA analysis with Tukey’s multiple comparison post hoc analysis was performed to compare frequencies of PBMC subpopulations between pwCIS at different time points and between those treated with GLAT and those who were untreated. A p-value of < 0.05 was considered statistically significant. Unless stated otherwise, data shows mean percentage and standard deviation (SD).

3.1 Demographics

Twenty pwCIS (male-to-female ratio 9:11) and fifteen HDs (male-to-female ratio 7:8) were included in this study (Table 1). No relevant difference between the two groups in sex (p=0.909) and age (p=0.548), was found. The baseline median EDSS for pwCIS was 2.00 (1.00 – 4.00) and 1.00 (0.00 – 2.50) for males and females, respectively. PBMCs from one CIS individual on baseline were not viable for flow cytometry analysis. Whilst in the GLAT group disease activity was high at the baseline, no major differences in disease activity were found between groups at any time point afterwards (see Table 1, and Supplementary Figure 2).

• Download figure
• Open in new tab

Supplement Figure 2: Disease Activity

3.2 Differences between treatment-naïve pwCIS and HDs at baseline

First, we analysed the difference between lymphocytes and monocytes in pwCIS and HDs at baseline. Results are shown in Figure 1 and Table 2. We found notable differences that showed lower percentages of lymphocytes in pwCIS compared to HDs (Figure 1A; 84.0 vs 91.7, p˂ 0.001). In contrast, the frequency of monocytes was higher among pwCIS than HDs (Figure 1A; 16.0 vs 8.2, p < 0.001).

• Download figure
• Open in new tab

Figure 1. Differences between pwCIS and HDs at Baseline.

A) Mean percentage of lymphocytes and monocytes between isolated PBMCs. Mean percentage of CD3⁺ T cells and CD19⁺CD20⁺ B cells in the lymphocyte population.) Mean percentage of different T cell subpopulations in the lymphocyte population. C) Mean percentage of different B cell subpopulations in the lymphocyte population. D) Mean percentage of different B cell subpopulations in the CD20⁺ subpopulation. Data are presented as mean percentage ± SD. *p˂ 0.05, **p˂ 0.01, ***p˂ 0.001. n=15.

Lymphocyte subsets were categorised as shown in Figure 1. Briefly, T cells were defined by the expression of CD3, monocytes by the expression of CD14 and B cells characterised by the absence of CD3 and CD14 with the expression of CD19 and CD20. Afterwards, a set of markers to discriminate between specific T and B cell subsets was established (see gating strategy in supplementary information). Our results did not show any difference in the percentage of CD3⁺ (Figure 1A) or CD3⁺CD4⁺ T cells among groups (Figure 1B). However, the fraction of CD4⁺CD38⁺HLADR⁺ within the T cell population was slightly higher in pwCIS (Figure 1B¹; 0.36 vs 0.23, p=0.015). Moreover, the percentage of CD4⁺CD25^high/CD127^low cells (defining Tregs) in pwCIS was somewhat higher than in HDs (Figure 1B; 3.9 vs 2.6, p˂ 0.001).

Interestingly, the percentage of CD19⁺CD20⁺ B cells among lymphocytes was notably lower in pwCIS compared to HDs (Figure 1A; 5.7 vs 9.3, p˂ 0.001). Percentages of both, CD20⁺CD27⁻ naïve and CD20⁺CD27⁺ memory B cells, were higher in HDs (naïve= 4.0 and memory= 2.5) than in pwCIS (naïve= 6.6 and memory= 1.5) (Figure 1C; p˂ 0.001 and p˂ 0.01, respectively); however, there was no substantial difference between the percentages of CD19⁺CD20⁻CD27^high cells (Plasmablasts; Figure 1C).

When we compared the proportions of B cell subsets defined by CD27 and IgD markers (CD27⁻IgD⁺ naïve, CD27⁺IgD⁺ pre-switched memory, CD27⁺IgD⁻ switched memory, and CD27⁻IgD⁻ atypical memory) between pwCIS and HDs, we did not find any substantial difference among these cohorts (Figure 1D).

3.3 Analysis of lymphocytes in treated-GLAT and untreated pwCIS

Subsequently, we analysed the difference between lymphocyte subsets and monocytes in two subgroups of patients; those who remained untreated (n=10) and those who received glatiramer acetate (GLAT; n=10). The outcomes analysed for the two subgroups of patients were either at 3 - 6 months (Baseline), at 12 months, and at 24 months after the first clinical event suggestive of demyelination of CNS (Table 1). Results are depicted in Figures 2 and 3 and Table 3.

• Download figure
• Open in new tab

Figure 2. Glatiramer acetate in pwCIS’s peripheral lymphocytes.

A) Mean percentage of lymphocytes and monocytes B) Mean percentage of CD3⁺ T cells and CD19⁺CD20⁺ B cells in the lymphocyte population. C) Mean percentage of different T cell subpopulations in the lymphocyte population. Data are presented as mean percentage ± SD.*p˂ 0.05, **p˂ 0.01, ***p˂ 0.001.

• Download figure
• Open in new tab

Figure 3. Glatiramer acetate effect in pwCIS’s peripheral B cells.

A) Mean percentage of different B cell subpopulations in the lymphocyte population. B) Mean percentage of different B cell subpopulations into the CD20⁺ subpopulation. Data are presented as mean percentage ± SD. *p˂ 0.05, **p˂ 0.01, ***p˂ 0.001.

First, we analysed the difference between lymphocytes and monocytes in pwCIS at different time points. We, nonetheless, did not find any important changes among groups (Figure 2A).

We categorised lymphocyte subsets as shown in Figures 2 and 3 (see gating strategy in Supplementary Figure 1). Overall, our analysis showed no difference within the T cell subpopulations (Figure 1). Nevertheless, a subpopulation of CD4⁺HLA-DR⁺ cells was slightly higher in the GLAT group (2.03) compared to those untreated (1.46), after 24 months (Figure 2C; p< 0.05).

Our results also showed that CD19⁺CD20⁺ B cells were present at similar frequencies at baseline for both patient groups. However, B cells increased after 24 months (GLAT=7.85 and Untreated= 5.47) in the GLAT group (Figure 2B; p< 0.05). The difference in the GLAT group appears to be secondary to an increased percentage of CD20⁺CD27⁻ naïve B cells, especially by 24 months (Figure 3A; p< 0.01). The percentage of CD20⁺CD27⁺ memory B cells was slightly reduced in the GLAT group after 24 months, whilst the untreated group had a little increase, but no difference was observed (Figure 3A). Finally, no differences were observed in plasmablast percentages at different time points (Figure 3A).

Finally, we compared the proportions of B cell subsets (CD27⁻IgD⁺ naïve, CD27⁺IgD⁺ pre-switched memory, CD27⁺IgD⁻ switched memory, and CD27⁻IgD⁻ atypical memory) between untreated and GLAT treated patients (Figure 3B). Interestingly, we confirmed the tendency of CD20⁺CD27⁻IgD⁺ naïve B cells to increase within the group under GLAT after 24 months (Figure 3B; p< 0.05). Our results also showed that the percentage of CD20⁺CD27⁺IgD⁺ pre-switched memory B cells in the GLAT group was lower after 24 months (Figure 3B; p< 0.01); in stark contrast to the GLAT group, pre-switched memory B cells tended to increase in the untreated group compared to baseline (Figure 3B; p< 0.01, Table 3). In addition, the percentage of CD20⁺CD27⁺IgD⁻ switched memory B cells also decreased during the intervals measured in the GLAT group, although statistically not different. Lastly, our analysis did not show any substantial change in the CD20⁺CD27⁻IgD⁻ atypical memory B cell subpopulation.

4.1 Differences between pwCIS and HDs

Results showed that the distribution of peripheral lymphocytes and monocytes differs between HDs and pwCIS. The percentage of lymphocytes in pwCIS is lower whilst the percentage of monocytes is higher compared to healthy individuals. Previous studies have highlighted the important role of monocytes in MS pathogenesis; they can contribute to the breakdown of the blood-brain barrier¹⁸. Furthermore, the T cell/monocyte ratio analysis has identified pwCIS at risk of rapid disease progression¹⁹. In line with previous studies, our results encourage further characterisation of monocyte subsets (in connection with dynamic changes among B cell subsets) in pwCIS in future clinical research projects with large samples.

Within the T cell compartment, the distribution of CD4⁺ cells did not show differences between HDs and pwCIS. Specifically, CD4⁺ subpopulations expressing activation markers, such as CD38 and HLA-DR, did not differ from the aforementioned groups. Interestingly, the proportion of T regulatory cells (Tregs) was higher in pwCIS compared to HDs at baseline. However, a steady, higher proportion of peripheral Tregs did not show any further benefit in preventing disease activity during the 24 months of follow-up (see Table 2).

Moreover, increasing evidence shows that several cell subsets contribute to the inflammatory cascade inside the CNS^20–22. Our results showed that the percentage of CD19⁺CD20⁺ B cells was lower in pwCIS compared to HDs. This is an interesting finding, considering that within the T cell compartment there was no or very little difference between groups. Our results are in line with Kreuzfelder et al. which showed that the percentage of peripheral CD19⁺ cells in pwMS is lower than in HDs²³. Whether this relatively low percentage of peripheral B cells is due to a higher monocyte proportion in pwCIS or to an active exchange between the periphery and CNS, including an increased influx of B cells into the inflamed CNS tissue, needs further investigation ²².

Intriguingly, a separate analysis of B cell subsets (with the CD19⁺CD20⁻ plasmablast population excluded,) showed that the percentage of peripheral switched memory, pre-switched memory, and atypical memory B cell subpopulations are slightly higher in pwCIS, but the difference is not significant. Unfortunately, the heterogeneity of protocols used by different studies to evaluate pwCIS in PBMCs^23–28 makes comparing the study results difficult.

4.2 Sub-analysis in pwCIS receiving GLAT

We also performed a subanalysis of pwCIS divided into two different groups, untreated and GLAT treated individuals, at three different time points (3-6, 12, and 24 months after disease onset). Taken as a whole, lymphocytes and monocytes, we did not observe any significant change between those untreated compared to GLAT treated individuals at different time points. We, nonetheless, observed an increase in the percentage of B cells (CD19⁺CD20⁺) in those patients being treated with GLAT. Our results showed no significant changes within the plasmablast population (CD19⁺CD20⁻CD27^high) subset.

B cell subsets seem to be mildly affected by GLAT but not equally. Our analysis showed a higher percentage of naïve (p˂ 0.05) and lower pre-switched memory (p˂ 0.01) B cells after 24 months in the GLAT group, compared to the untreated group. Furthermore, whilst switched memory B cells seemed to be reduced after 24 months in people receiving GLAT treatment, our analysis of variance did not show a substantial change among groups at different time points.

Due to this study’s explorative nature and aim, it is difficult to make assumptions about the therapeutic effect of GLAT in pwCIS, and therefore, our results are just observational; interindividual differences can also drive these changes. Interestingly, it has been observed that GLAT could affect various aspects of dysregulated B cells subsets in pwMS ²⁴.

4.3 B cells and immune status monitoring in pwCIS and MS

Whilst B cells could overreact, generating a hostile and extreme autoimmune reaction, evidence has also shown that some subsets could regulate autoreactive T cell responses ^25,26. For instance, it is believed that alemtuzumab depletes CD19⁺CD27⁺ B memory cells and allows a rapid repopulation of immature B cells ⁴. Furthermore, a new study by Ruck et al., using an in-depth multidimensional immune phenotyping analysis, showed that the presence of hyperexpanded T cell clones might predict the development of secondary autoimmunity in some patients under alemtuzumab treatment. The effect of these T cell clones in MS pathophysiology is still unclear²⁷.

The depletion of peripheral B cells by anti-CD20 therapies shows that B cells play a key role in MS pathogenesis ²⁸. However, the depletion of B cell subsets might not be enough to maintain an effective therapeutic outcome. Depletion of other reactive T cell subsets to obtain a more effective response might be needed. ^29,30 Moreover, some studies have found that the production of interleukin-10 (IL-10) by B cell subsets is a key factor for regulatory responses^26,31,32. It has been proposed that IL-10 produced by naïve and memory B cells has a direct effect on activated Th1/Th17 and T regulatory cells to regulate immune responses in both health and pathological conditions³².

It is of note that we were unable to observe plasma cells in all measured samples. It seems that plasma cells are less likely to survive the freezing/thawing procedures applied to PBMCs. Moreover, in our study, we did not perform absolute counts due to the use of frozen samples. We propose that future lymphocyte immunophenotyping shall be performed in fresh blood and include absolute counting. This would not have been possible with the current methodology of this cohort. Nevertheless, and despite the limitations of this observational retrospective study and small sample size, we consider our results could encourage further research by using methodological assets of phenotypic and molecular studies to identify cellular biomarkers in these MS entities. In this study, we show that lymphocyte characterization may be useful in characterising differences in peripheral lymphocyte subsets in pwCIS after a first demyelinating event.