T-Cell Specificity Influences Disease Heterogeneity in Multiple Sclerosis
Citation Manager Formats
Make Comment
See Comments

Abstract
Background and Objectives Encouraged by the enormous progress that the identification of specific autoantigens added to the understanding of neurologic autoimmune diseases, we undertook here an in-depth study of T-cell specificities in the autoimmune disease multiple sclerosis (MS), for which the spectrum of responsible autoantigens is not fully defined yet. The identification of target antigens in MS is crucial for therapeutic strategies aimed to induce antigen-specific tolerance. In addition, knowledge of relevant T-cell targets can improve our understanding of disease heterogeneity, a hallmark of MS that complicates clinical management.
Methods The proliferative response and interferon gamma (IFN-γ) release of CSF-infiltrating CD4+ T cells from patients with MS against several autoantigens was used to identify patients with different intrathecal T-cell specificities. Fresh CSF-infiltrating and paired circulating lymphocytes in these patients were characterized in depth by ex vivo immunophenotyping and transcriptome analysis of relevant T-cell subsets. Further examination of these patients included CSF markers of inflammation and neurodegeneration and a detailed characterization with respect to demographic, clinical, and MRI features.
Results By testing CSF-infiltrating CD4+ T cells from 105 patients with MS against seven long-known myelin and five recently described GDP-l-fucose synthase peptides, we identified GDP-l-fucose synthase and myelin oligodendrocyte glycoprotein (35-55) responder patients. Immunophenotyping of CSF and paired blood samples in these patients revealed a significant expansion of an effector memory (CCR7− CD45RA−) CD27− Th1 CD4+ cell subset in GDP-l-fucose synthase responders. Subsequent transcriptome analysis of this subset demonstrated expression of Th1 and cytotoxicity-associated genes. Patients with different intrathecal T-cell specificities also differ regarding inflammation- and neurodegeneration-associated biomarkers, imaging findings, expression of HLA class II alleles, and seasonal distribution of the time of the lumbar puncture.
Discussion Our observations reveal an association between autoantigen reactivity and features of disease heterogeneity that strongly supports an important role of T-cell specificity in MS pathogenesis. These data have the potential to improve patient classification in clinical practice and to guide the development of antigen-specific tolerization strategies.
Glossary
- CP=
- control patient;
- EM=
- effector memory;
- FLAIR=
- fluid-attenuated inversion recovery;
- GDP-L-FS=
- GDP-l-fucose synthase;
- HD=
- healthy donor;
- HLA=
- human leukocyte antigen;
- IFN=
- interferon;
- IgG=
- immunoglobulin G;
- MBP=
- myelin basic protein;
- LP=
- lumbar puncture;
- MOG=
- myelin oligodendrocyte glycoprotein;
- MOGAD=
- myelin oligodendrocyte glycoprotein antibody–associated disease;
- MPRAGE=
- Magnetization Prepared - RApid Gradient Echo;
- MS=
- multiple sclerosis;
- PBMC=
- peripheral blood mononuclear cell;
- RA=
- rheumatoid arthritis;
- SI=
- stimulation index;
- TF=
- transcription factor
Multiple sclerosis (MS) is considered an immune-mediated demyelinating disease of the CNS.1 There is general agreement about the relevance of B cells,2 CD8+ T cells,3,4 and particularly CD4+ T cells in MS pathogenesis. However, despite the fact that myelin-specific CD4+ T cells can exacerbate MS5 and that specific T cells can induce disease in different animal models of MS,6,7 there is less consensus about the relevance of CD4+ T-cell specificity.
Disease heterogeneity is a hallmark of MS and likely stems from the interplay of a complex genetic background8,9 and environmental risk factors.10 Although it is well described that disease heterogeneity encompasses composition of tissue lesions, imaging findings, clinical presentation, disease course, and treatment responsiveness, the factors contributing to this heterogeneity, including a putative involvement of T-cell specificity, are not clear. From studies in Experimental Autoimmune Encephalomyelitis (EAE), the most commonly used animal model in MS,6,11 it is well known that both the rodent strain, i.e., the genetic background, and the antigen used to induce EAE influence disease course, severity, lesion distribution, and pathologic mechanisms,7,12,13 linking T-cell specificity and disease heterogeneity.
The identification of the target antigens in neurologic autoimmune diseases14,-,16 has led to great progress in understanding, diagnosing, and treating these diseases. To better understand and treat MS, researchers have investigated CD4+ T-cell specificity for years, albeit with less success than expected. The difficulties in obtaining CNS-infiltrating CD4+ T cells, the overall low antigen avidity of self–antigen–reactive T cells, and the limited information about T cell–specific autoantigens in MS have hampered research progress. In addition to the long-known myelin peptides,17 few T cell–specific autoantigens have been identified so far. We recently identified GDP-l-fucose synthase (GDP-L-FS)-derived peptides as putative autoantigens in human leukocyte antigen (HLA)-DRB3* patients with MS.18
To evaluate the relevance of CD4+ T-cell specificity in MS, we studied their link with features of disease heterogeneity. We identified patients with different CSF-infiltrating CD4+ T-cell responses against myelin and GDP-L-FS peptides and characterized them by ex vivo immunophenotyping of CSF and paired blood samples as well as with respect to demographic, clinical, and MRI features.
Methods
Patient Samples
Paired CSF and blood samples were collected from 105 untreated patients with MS, only CSF from 11 control patients (CPs), and 10 patients with myelin oligodendrocyte glycoprotein antibody–associated disease (MOGAD) negative for anti-AQP4 antibodies and only peripheral blood mononuclear cells (PBMCs) from four healthy donors (HDs) (eTable 1, links.lww.com/NXI/A590). Eighty-four patients with MS (80%) had never been treated, whereas 21 (20%) were previously treated but considered untreated at the time of lumbar puncture (eTable 1, links.lww.com/NXI/A590). Patients and controls were recruited from the NIMS-Neuroimmunology and MS Research Section, Department of Neurology, University Hospital Zurich. MS diagnosis was based on the revised McDonald criteria.19
Standard Protocol Approvals, Registrations, and Patient Consents
The Cantonal Ethics Committee of Zurich approved the study procedures (EC-No. 2013-0001). Informed consent was obtained from all patients.
CSF and Serum Measures
Albumin quotients and CSF-specific oligoclonal bands were analyzed as previously reported.20 Intrathecal Ig synthesis (Ig[loc]) was calculated according to Reiber.21
Cell Culture and Stimulation
PBMCs were freshly isolated using Ficoll (Eurobio, Germany) density-gradient centrifugation. CSF-infiltrating CD4+ T cells were expanded with phytohaemagglutinin in a single round of stimulation using a previously reported protocol aimed to reduce T cell receptor repertoire bias.18 6 × 104 phytohaemagglutinin-expanded CSF-infiltrating CD4+ T cells were seeded in quadruplicate with 2 × 105 irradiated autologous PBMCs as antigen-presenting cells and stimulated with peptides (eTable 2, links.lww.com/NXI/A590) at a final concentration of 10 μM and with a T Cell Activation Kit (anti-CD3, anti-CD28, and anti-CD2 beads) (Miltenyi Biotec) as positive control.
Proliferative Responses
We measured proliferation 72 hours after stimulation using 3H-thymidine (Hartmann Analytic, Braunschweig, Germany) as previously reported.18 The stimulation index of single wells (SI) was calculated as follows: SI = (cpm well with peptide)/(Mean cpm wells without peptide). A well was considered positive if SI ≥2. A patient was considered positive for a peptide if the mean SI of the four wells was ≥2 and at least 3 of the four wells showed a SI ≥2.
Quantification of Cytokines
Release of interferon gamma (IFN-γ) was measured using the ELISA MAX Deluxe Set Human IFN-γ (BioLegend, San Diego, CA) according to the manufacturer's instructions. A well was considered positive if IFN-γ ≥20 pg/mL, and a patient was considered positive for a peptide if the mean IFN-γ of the four wells was ≥20 pg/mL and at least 3 of the 4 wells were positive.
Cytokines in supernatants and CSF were measured using the Human T Helper Cytokine Panel LEGENDplex bead-based immunoassay (BioLegend) according to the manufacturer's instructions.
Immunophenotyping
Ex vivo immunophenotyping of intrathecal and paired circulating lymphocytes was performed as previously reported20,22 in a subcohort of 66 patients with MS (eTable 1, links.lww.com/NXI/A590). We added SPHERO AccuCount Particles (Spherotech, Inc., Lake Forest, IL) to count absolute numbers following the manufacturer's instructions. An LSRFortessa cytometer (BD Biosciences, Franklin Lakes, NJ) was used to acquire samples and FACSDiva (BD) and FlowJo (TreeStar Inc., Ashland, OR) software for the analysis. eFigure 1, links.lww.com/NXI/A590, summarizes the gating strategy.
Sorting of Cell Subsets
PBMCs from four patients with MS and four age- and sex-matched HDs were labeled with antibodies against CD4 (allophycocyanin), CD45RA (BV711), CCR7 (BV421), CD27 (phycoerythrin), and CD28 (allophycocyanin-Cy7), all from BioLegend and with live-dead aqua dye (Invitrogen). CD4+ CCR7− CD45RA− CD28+ CD27+ and CD4+ CCR7− CD45RA− CD28+ CD27− live cells were sorted using 100-μm sorting chips in a Sony SH800SFP cell sorter (4 lasers, Sony). Twenty thousand sorted cells from each cell population were transferred to an RNase-free tube, resuspended in Qiazol (Qiagen, Germany), and frozen at −80°C.
RNA Extraction, Sequencing, and Analysis
RNA extraction from frozen cell pellets was performed using the PicoPure RNA Isolation Kit (Life Technologies) according to the manufacturer's instruction. RNA sequencing (RNAseq) was performed using Illumina Sequencing 200M at the Functional Genomics Center Zurich as previously reported.23 The RNAseq data analysis consisted of (1) cleaning of raw reads using Trimmomatic (version 0.36)24; (2) pseudoalignment of sequences to the human reference genome (build GRCh38.p13, gene annotation from GENCODE Release 32) and gene expression quantification using Kallisto (version 0.44)25; (3) read alignment using STAR (v2.7.3)26; and (4) differentially expressed gene detection using a count-based negative binomial model implemented in the software package EdgeR (R version: 3.6.1, EdgeR version: 3.28.0).27 We used a generalized linear model adapted for overdispersed data to assess differential expression. Genes showing p value <0.001 for altered expression were considered differentially expressed.
Enzyme-Linked Immunosorbent Assays
The following proteins were quantified using enzyme-linked immunosorbent assay kits: neurofilament light chain (NfL) (Human Diagnostics, Umea, Sweden); CXCL13/BLC/BCA-1, granzyme A, and granulysin (R&D System, MN); chitinase 3 like 1 (MicroVue, Athens, OH); perforin, granzyme B, and granzyme H (Invitrogen–Thermo Fisher Scientific, MA), were used according to the manufacturer's instructions.
HLA Typing
Patients were typed for HLA-class I (A* and B*) and HLA class II (DRB1*, DRB3*, DRB4*, DRB5*, DQA1*, and DQB1*) as previously reported.20
MOG Antibody Detection Assay
Serum and CSF samples were analyzed for myelin oligodendrocyte glycoprotein (MOG) immunoglobulin G (IgG) antibodies as described before.28 For screening, serum samples were diluted 1:20 and 1:40 and CSF samples 1:2. Positive samples were end point titrated, and MOG-IgG positivity was confirmed using an anti-human IgG(Fc) specific secondary antibody as recently described.28
MRI
Eight patients were scanned with a 3T Philips Ingenia and 22 patients with 3T Siemens Skyra. The MRI protocol included a 3D pre- and post-gadolinium contrast-enhanced gradient echo pulse sequence (MPRAGE) as well as a 3D fluid-attenuated inversion recovery (FLAIR) sequence.
From FLAIR, the number and the total volume in milliliters of all hyperintense lesions were determined by an automatic algorithm based on convolutional neural networks.29 All results were manually corrected by two experienced technical raters. Differences in the corrections were resolved by consensus in a second reading phase. Similarly, the number of contrast-enhancing lesions was determined.
Whole-brain, gray matter, and thalamic volumes in milliliters were determined on the precontrast MPRAGE image using the automatic processing pipeline Biometrica MS analysis platform (version 2.1, jung diagnostics GmbH, Hamburg, Germany).30
Statistics
We used GraphPad Prism 8.0 (GraphPad Software, La Jolla, CA) to perform the statistical analysis. We used the unpaired t test to compare two groups of normally distributed variables and U test (Mann-Whitney) for non-normally distributed variables. For the comparison of more than 2 groups of patients, we used the Kruskal-Wallis test for non-normally distributed variables. We used Spearman r to test the linear correlation between non-normally distributed variables. The significance level was set at p < 0.05. Associations between patient specificity, seasonal distribution of lumbar puncture (LP), and HLA were performed using the Fisher exact test with 5% significance.
Results
Identification of GDP-L-FS and Myelin Responders
Proliferation and IFN-γ release of CSF-infiltrating CD4+ T cells from 105 patients with MS against 5 immunodominant GDP-L-FS18 and 7 myelin peptides31 are shown in eFigure 2A, links.lww.com/NXI/A590. GDP-L-FS peptides induced the highest frequency of positive wells using proliferation as readout and MOG(35-55) peptide using IFN-γ release (Table 1). We also analyzed responses to a viral/bacterial peptide pool cytomegalovirus, Epstein-Barr virus and influenza virus (CEF) (eFigure 2A, links.lww.com/NXI/A590, and Table 1). IFN-γ secretion identified more positive wells than proliferation for myelin and CEF peptides, but comparable results for GDP-L-FS peptides (eFigure 2B, links.lww.com/NXI/A590). Accordingly, the strongest correlation between proliferation and IFN-γ has been found for GDP-L-FS peptides (eFigure 2C, links.lww.com/NXI/A590).
Percentage of Positive Wells and Patients for Each Peptide Using Proliferation and IFN-γ Release
The frequency of patients with positive responses to the different peptides is summarized in Table 1. In 14 patients (13.34%), the main reactivity in proliferation and/or IFN-γ release was against several GDP-L-FS peptides and was classified as GDP-L-FS-responders (Figure 1). Among these, 6 patients (42.8%) showed also positive responses to myelin peptides. The association between GDP-L-FS reactivity and reactivity to myelin peptides was significant (p ≤ 0.0001, Fisher exact test). Four patients (3.8%) responded only to myelin basic protein (MBP) peptides and were classified as MBP responders, whereas 11 patients (10.4%), responding only against MOG (35-55), were classified as MOG(35-55) responders. Seventy-six patients (72.4%) did not react to any autoantigen and were classified as nonresponders (Figure 1). The response to CEF peptides was comparable in the different groups of patients with the exception of MOG(35-55) responders, in whom it was lower (Figure 1). Because of the low number of MBP responders (<5), this group was not analyzed further.
Checkerboard graphs illustrating the response of each patient with MS to the individual peptides. Filled cells are positive responses in proliferation (upper checkerboard graph) and IFN-γ release (lower checkerboard graph). The number of GDP-L-FS, MBP, MOG(35-55), and nonresponders is shown in each group. Peptides color code, GDP-L-FS (blue), MBP (green), MOG(1-20) (black), MOG(35-55) (red), Myelin proteolipid protein (PLP) (violet), and CEF (orange). GDP-L-FS = GDP-l-fucose synthase; MS = multiple sclerosis.
Distinct CSF-Infiltrating and Circulating T Cells in GDP-L-FS Responders
Immunophenotyping of CSF-infiltrating and circulating lymphocytes was performed in a subcohort of patients with MS (eTable 1, links.lww.com/NXI/A590) consisting in GDP-L-FS (n = 7), MOG(35-55) (n = 7), and nonresponder (n = 52) patients (eTable 3, links.lww.com/NXI/A590). Effector memory (EM, CCR7− CD45RA−) CD4+ T cells expressing CD28 but not CD27 (EM CD27−) were significantly more abundant in the CSF of GDP-L-FS responders (Figure 2, A and B, and eFigure 3A, links.lww.com/NXI/A590). Among these EM CD27− cells, only cells with a Th1 (CCR6− CCR4− CRTh2−) functional phenotype (EM CD27− Th1), showed significantly higher frequencies (Figure 2B and eFigure 3B, links.lww.com/NXI/A590). Analysis of paired blood samples demonstrated that these cells were also significantly more abundant in frequency and absolute counts in the blood (Figure 2, A and C, and eFigure 3C, links.lww.com/NXI/A590). CSF-infiltrating and circulating CD4+ T cells with a TEMRA (CCR7− CD45RA+) differentiation state and the costimulatory molecules CD28+ CD27− and CD28− were also more abundant in GDP-L-FS responders (eFigure 3, A and C, links.lww.com/NXI/A590).
A) Dot plot showing CD28 and CD27 expression on CSF-infiltrating and peripheral circulating EM CD4+ cells from GDP-L-FS and nonresponders. Percentages of EM CD27− cells are shown. (B and C) Frequencies of CSF-infiltrating (B) and frequencies and absolute numbers of peripheral circulating (C) EM CD27− and EM CD27− Th1 cells in GDP-L-FS (blue), MOG(35-55) (red), and nonresponders (black). Cell counts were determined using SPHERO AccuCount Particles. Each dot corresponds to a single patient, and bars show means. The Kruskal-Wallis test was used to compare patients. Statistical significance (*p < 0.05, **p < 0.01, ***p < 0.001) is shown. EM = effector memory; GDP-L-FS = GDP-l-fucose synthase.
Transcriptional Signatures of EM CD27− CD4+ T Cells in GDP-L-FS Responders
EM CD27− and CD27+ CD4+ T cells were sorted from the peripheral blood of four GDP-L-FS responders (eFigure 4A, links.lww.com/NXI/A590) and the transcriptome analyzed by RNAseq. We identified 265 differentially expressed genes (fold change 1.5, p < 0.001), 119 upregulated and 146 downregulated, in EM CD27− vs EM CD27+ cells (Figure 3A). The upregulated genes included (1) transcripts linked to cytotoxic CD8+32 and CD4+33 T cells such as ADGRG1 (GPR56), ADGRG5 (GPR114), CCL4, CCL5, CST7, CTSW, CX3CR1, ENC1, FCRL6, FGFBP2, GNLY, GZMB, GZMH, MYO6, PRF1, PRSS23, S1PR5, SLAMF7, SPON2, TGFBR3, TRGC2, and ZEB2; (2) transcription factors (TFs) important in the development of cytotoxic CD4+ T cells such as Eomesodermin (Eomes)34 and T-bet (encoded by TBX21)35; and (3) other genes associated with Th1 cells such as interferon gamma gene and granzyme A gene (GZMA). In contrast, cytokines, chemokine receptors, and TFs associated with other Th subsets (CCR4, CCR6, IL4R, IL6R, GATA3, and FoxP3) were downregulated in EM CD27− cells.
(A) Heat map shows the row-wise z scores of 265 transcripts differentially expressed between EM CD27+ (columns 1–4) and EM CD27− (columns 5–8) cells (Log2Ratio >0.5, p < 0.001) from four GDP-L-FS responder patients. Z scores of these genes in EM CD27+ (columns A–D) and EM CD27− (columns E–H) cells from HDs are also shown. Z scores of selected genes associated with cytotoxicity, Th1, and other Th subsets are shown in detail. In bold are genes that also were differentially expressed in EM CD27− from GDP-L-FS responders vs EM CD27− from HDs (Log2Ratio >0.5, p < 0.001). (B) Distribution of log2(fragments per kilobase per million + 0.1) gene expression of Th1/cytotoxic genes and genes associated with other Th subsets in EM CD27+ and CD27− cells from GDP-L-FS responders and HDs. EM = effector memory; GDP-L-FS = GDP-l-fucose synthase; HDs = healthy donors.
We also analyzed the same sorted EM subpopulations from four HDs (eFigure 4A, links.lww.com/NXI/A590). The 265 genes that differentiate EM CD27− and CD27+ cells in GDP-L-FS responders did not discriminate EM subpopulations in HDs (Figure 3A). Differential gene expression analysis of EM CD27− cells from GDP-L-FS responders and HDs also identified 145 differentially expressed genes (fold change 1.5, p < 0.001), with 66 upregulated in the GDP-L-FS-responders that include 19 genes linked to cytotoxicity and Th1 T cells and 79 downregulated including genes associated with other T-cell subsets (eFigure 4B, links.lww.com/NXI/A590, and Figure 3A).
The transcription of genes characteristically associated with Th1 cells and cytotoxicity (TBX21, GZMA, EOMES, GNLY, GZMH, and SLAMF7) was overall increased, and transcription of genes associated with other Th subsets (GATA3, CCR4, and CCR6) decreased in EM CD27− cells from GDP-L-FS responders compared with the same cells from HDs and with EM CD27+ cells (Figure 3B).
Characterization of GDP-L-FS- and MOG(35-55)-Specific Responses
The functional phenotype of GDP-L-FS- and MOG(35-55)-specific CD4+ T cells was analyzed in the supernatant of positive wells stimulated with the cognate antigen (Figure 4A). GDP-L-FS-specific cells released significantly higher amounts of IL-2 than MOG(35-55)-specific cells. This higher IL-2 is most likely involved in the significantly stronger ability of GDP-L-FS-specific cells to proliferate in response to specific peptides, whereas the proliferation to unspecific stimulating beads was comparable in GDP-L-FS- and MOG(35-55)-specific cells (Figure 4B). IFN-γ was the main cytokine released by GDP-L-FS-specific cells indicating a Th1 functional phenotype. MOG(35-55)-specific cells, in addition to IFN-γ, released cytokines associated with other Th subsets (IL-9, IL-6, and IL-10), suggesting a polyfunctional phenotype (Figure 4A). IL-10 also was significantly higher in the CSF of MOG(35-55) responder and Th2-associated cytokines (Figure 4A). Because some of the cytokines released by MOG(35-55)-specific cells have been found elevated in CSF from patients with MOGAD,36 and because both groups of patients share the target antigen, we verified that none of the patients with MS had positive anti-MOG IgG titers in serum or CSF (eTable 4, links.lww.com/NXI/A590) and also analyzed the cytokine profile in the CSF from patients with MOGAD (eTable 1 and eFigure 5A, links.lww.com/NXI/A590). Although cytokines such as IL-9, IL-22, and IL-6 were high in both groups, the Th2-associated cytokines IL-4 and IL-10 were significantly higher in MOG(35-55) responders (eFigure 5A, links.lww.com/NXI/A590).
(A) Upper graphs, cytokines released by CSF-infiltrating CD4+ T cells from GDP-L-FS and MOG(35-55) responders after stimulation with specific peptides (GDP-L-FS and MOG[35-55]) presented by autologous PBMCs. Four wells were pooled for each patient. Lower graphs, cytokines present in the CSF of GDP-L-FS and MOG(35-55) responders. Cytokines are expressed as pg/ml. (B) Proliferative responses of CSF-infiltrating CD4+ T cells from GDP-L-FS and MOG(35-55) responders after stimulation with the specific peptides presented by autologous PBMCs or anti-CD3, anti-CD28, and anti-CD2 stimulating beads. Proliferative responses are expressed as SI. (C) Frequencies of CSF-infiltrating and blood circulating naive (immunoglobulin D + CD27−) B (CD19+ CD138−) cells in GDP-L-FS and MOG(35-55) responders. Patient color code, GDP-L-FS responders (blue) and MOG(35-55)-responders (red). Each dot in the graphs corresponds to a single patient, and the bars show the mean. The Mann-Whitney test was used to compare GDP-L-FS and MOG(35-55) responders. Statistical significance (*p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001) is shown. GDP-L-FS = GDP-l-fucose synthase; PBMCs = peripheral blood mononuclear cells.
Supporting a role of humoral immunity in MOG(35-55) responders, we found significantly higher frequencies of CSF-infiltrating but not circulating naive B cells in these patients (Figure 4C). The amounts in the CSF of CXCL1337 were higher in MOG(35-55) responders and chitinase 3 like 1, a general marker of inflammation,38 in GDP-L-FS responders when compared with CPs (eTable 1, links.lww.com/NXI/A590) (eFigure 5B, links.lww.com/NXI/A590). Both groups of patients showed significantly higher intrathecal IgG synthesis than CP (eFigure 5B, links.lww.com/NXI/A590).
Further CSF analysis of cytotoxic and neurodegenerative markers revealed undetectable amounts of perforin and granzyme B, no significant differences for granulysin, and higher levels of granzyme H, granzyme A, and NfL39 in GDP-L-FS responders (eFigure 5C, links.lww.com/NXI/A590).
Association of GDP-L-FS Reactivity and HLA-DRB3*02:02/03:01
All patients were HLA typed (eTable 5, links.lww.com/NXI/A590). Figure 5A summarizes frequencies of GDP-L-FS, MOG(35-55), and nonresponders expressing the MS-associated DR15 haplotype and DRB3*02:02/03:01 alleles. Based on two previously typed MS cohorts (German cohort [n = 1,270] and Swiss cohort [n = 367]), we expected that 40%–50% of patients with MS express the DR15 haplotype and 30%–40% the DRB3*02:02/03:01 alleles. The frequency of GDP-L-FS responders expressing the DR15 haplotype (21.4%) was lower than expected, whereas those expressing the DRB3*02:02/03:01 alleles (92.8%) were much higher. The association between GDP-L-FS specificity and DRB3*02:02/03:01 alleles was significant (p ≤ 0.0001, Fisher exact test, Figure 5A), confirming our previous results.18
(A) Frequency of patients expressing the DR15 (left) and the DRB3*02:02/03:01 (right) genes in GDP-L-FS (blue), MOG(35-55) (red), and non- (black) responders as well as in two reference cohorts. Number of patients responding or not to GDP-L-FS and expressing (Y, yes) or not (N, not) the DRB3*02:02/03:01 genes and p values for Fisher exact tests are shown. (B) CSF-infiltrating CD4+ T-cell responses (proliferation [SI] and IFN-γ release) to GDP-L-FS (blue), MBP (green), MOG(1-20) (black), MOG(35-55) (red), PLP (violet), and CEF (orange) peptides and anti-CD2, anti-CD28, and anti-CD3 stimulatory beads (gray) from patients expressing or not the DRB3*02:02/03:01 genes. Each dot in the graphs corresponds to a single well. The Mann-Whitney test was used to compare patients expressing DRB3*02:02/03:01 and other DR. Statistical significance (***p < 0.001 and ****p < 0.0001) is shown. GDP-L-FS = GDP-l-fucose synthase; SI = single wells.
Figure 5B shows CSF-infiltrating CD4+ T-cell responses to GDP-L-FS, myelin, and CEF peptides as well as stimulating beads in patients expressing or not DRB3*02:02/03:01 alleles. Only the response to GDP-L-FS peptides, both by proliferation and IFN-γ release, was significantly higher in DRB3*02:02/03:01 patients than in patients expressing other HLA class II alleles.
Demographic and Clinical Features of Patients With Different Specificity
Demographic and clinical characteristics of GDP-L-FS, MOG(35-55), and nonresponder patients did not show significant differences (Table 2). Further characterization revealed a different seasonal distribution regarding the time point when the LP had been performed (Figure 6A), with a significant association between GDP-L-FS specificity and LP during winter/spring months (≤0.0001, Fisher exact test).
Features of Patients With Multiple Sclerosis With Different Specificity
(A) Monthly distribution of LPs and frequency of samples obtained in winter/spring and summer/autumn in GDP-L-FS, MOG(35-55), and nonresponders. Red shadows show months with higher differences. Number of patients responding or not to GDP-L-FS in which LP was obtained in winter/spring (Y, yes) or not (N, not) and p values for Fisher exact tests are shown. (B) Number of total contrast-enhancing T1 lesions and total volume of flair T2 lesions (expressed in mL) in GDP-L-FS- and MOG(35-55) responders. Each dot in the graphs corresponds to a single patient, and the bars show the mean. We used t test for normally distributed variables to compare two groups of patients. Statistical significance (**p < 0.01) is shown. GDP-L-FS = GDP-l-fucose synthase.
We also examined markers of compromised blood brain barrier permeability, disease activity, and several MRI parameters (Table 2). Although most of these markers were comparable between GDP-L-FS and MOG(35-55) responders, the semiautomated and blinded analysis of retrospectively collected brain MRI scans from 6 GDP-L-FS and 8 MOG(35-55) responders revealed statistically significant differences. The total number of contrast-enhancing T1 lesions and the total volume of FLAIR T2 lesions (Figure 6B) were significantly higher in GDP-L-FS compared with MOG(35-55) responders.
Discussion
This study represents the most comprehensive analysis of T-cell specificity performed in MS to date because it (1) analyzes CSF-infiltrating CD4+ T cells instead of PBMCs, (2) evaluates the reactivity not only against seven long-known immunodominant/encephalitogenic myelin peptides31 but also to five recently discovered immunodominant GDP-L-FS peptides,18 (3) uses as a readout IFN-γ release in addition to proliferation based on the relevant role of Th1 cells in MS,5,18,23,31 and (4) is performed in a large cohort of 105 patients with MS. Proliferation confirmed a stronger reactivity against GDP-L-FS peptides,18 whereas IFN-γ release identified MOG(35-55) as the strongest stimulatory peptide. IFN-γ release identified more positive wells than proliferation and should be considered in future studies to avoid underestimating responses. The fact that only 27.6% of patients with MS showed positive responses suggests that we could have missed certain T cells by their too low frequencies or by not testing a broader range of cytokines, target antigens, or epitopes.40 In addition, due to the relapsing-remitting nature of inflammation in MS, it is possible that we have missed certain T cells because antigen-specific T cells passage only transiently through the CSF and are not present at the time of the LP. The fact that only 38% of the patients showed T-cell responses against strong viral/bacterial recall antigens supports this notion. Despite these limitations, we identified 14 GDP-L-FS, 4 MBP, and 11 MOG(35-55) responders.
GDP-L-FS responders frequently also recognized myelin peptides. The different sequences of the stimulatory GDP-L-FS and myelin peptides render it unlikely that cross-recognition by single TCCs is behind this association. Instead, a Th1/cytotoxic response that results in the release of new peptides from damaged tissue and the activation of new autoreactive T cells, i.e., epitope spreading,41 might provide a more plausible explanation. Of interest, ex vivo immunophenotyping demonstrated that EM CD27− Th1 CD4+ were significantly more abundant in the CSF and also peripheral blood of GDP-L-FS-responders and that these cells in the blood expressed Th1- and cytotoxicity-associated genes.33 Whether this T-cell subset is enriched in GDP-L-FS-specific cells is still unknown and requires further investigation. Remarkably, high dimensional single cell analysis of CD4+ T cells from patients with rheumatoid arthritis (RA) identified an equivalent expanded EM CD27− CD4+ T-cell subset with a very similar transcriptome profile42 supporting a broader role of this cell subset in autoimmunity. Furthermore, both GDP-L-FS responders and patients with RA also have higher frequencies of TEMRA CD28+ CD27− and CD28- cells that have been associated with damage and clinical progression in RA.43
Cytokine analysis showed a relevant role of Th1 responses in GDP-L-FS responders and Th2 in MOG(35-55) responders. GDP-L-FS- and MOG(35-55)-specific cells did not release IL-17, a cytokine reported to be relevant in MS.44 Th17 and Th1* cells in patients with MS most likely have specificities that have not been included in this study. The higher intrathecal amounts of Th2 cytokines in MOG(35-55) responders suggest a role of humoral immunity in these patients. Accordingly, they showed higher frequencies of CSF-infiltrating naive B cells and higher CXCL13 levels.37 These patients also showed higher intrathecal amounts of IL-10 that might regulate Th2 responses as it has been reported in allergy.45
GDP-L-FS-specific responses were associated with DRB3*02:02/03:01 alleles, confirming previous results,18 and were significantly more frequent in CSF samples obtained during winter and spring. Seasonality of infections as triggers of autoreactive T cells might underly this annual distribution, but gut microbiota dysbiosis due to seasonal changes in diet46 and vitamin D47 should also be considered because CSF-infiltrating CD4+ T cells from GDP-L-FS responders cross-recognized homologous bacterial peptides.18 Finally, GDP-L-FS and MOG(35-55) responders also differed regarding MRI findings. It is remarkable that, despite the low number of patients from whom brain MRI scans were available at the time of LP, GDP-L-FS responders showed significantly higher numbers of contrast-enhancing T1 lesions and higher volumes of FLAIR T2 lesions, indicating higher inflammation and more extensive demyelination accordingly with more damaging immune responses.
In summary, our results reveal, an association between T-cell specificity and features of disease heterogeneity in a T cell mediated autoimmune disease. A parallel association between disease severity and autoantibodies against citrullinated peptides is well described in patients with RA.48,-,50 The utility of this long-known association in the clinical management of patients with RA indicates the potential of our findings to improve clinical care and treatment decisions. In addition, our results might have important implications for personalized treatment approaches aiming at antigen-specific tolerance.
Study Funding
This work has been funded by European Research Council Advanced Grant (340733) (R.M.), Clinical Research Priority Program MS (CRPPMS) and (CRPP) PrecisionMS of the University Zurich (R.M., M.S., and A.L.), Clinical, Swiss National Science Foundation (Sinergia UnmetMS) (R.M. and M.S.) and the Swiss MS Society (R.M.).
Disclosure
All authors except R.M. and M.S. declare that they have no competing financial interests in the context of this work. R.M. and M.S. hold a patent on immunodominant proteins and fragments in autoimmune diseases (EP18180326.3; 2018). Go to Neurology.org/NN for full disclosures.
Appendix Authors

Footnotes
Go to Neurology.org/NN for full disclosures. Funding information is provided at the end of the article.
↵* These authors contributed equally to this work.
The Article Processing Charge was funded by the authors.
- Received April 26, 2021.
- Accepted in final form July 15, 2021.
- Copyright © 2021 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.
This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND), which permits downloading and sharing the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.
References
- 1.↵
- 2.↵
- Sospedra M
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- Miller SD,
- Karpus WJ
- 7.↵
- 8.↵
- Beecham AH,
- Patsopoulos NA, et al.
- 9.↵International Multiple Sclerosis Genetics Consortium. Multiple sclerosis genomic map implicates peripheral immune cells and microglia in susceptibility. Science. 2019;365(6460):eaav7188.
- 10.↵
- 11.↵
- Bjelobaba I,
- Begovic-Kupresanin V,
- Pekovic S,
- Lavrnja I
- 12.↵
- 13.↵
- Manogaran P,
- Samardzija M,
- Schad AN, et al
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- Hohlfeld R,
- Dornmair K,
- Meinl E,
- Wekerle H
- 18.↵
- Planas R,
- Santos R,
- Tomas-Ojer P, et al.
- 19.↵
- 20.↵
- Puthenparampil M,
- Tomas-Ojer P,
- Hornemann T, et al
- 21.↵
- 22.↵
- Brodie T,
- Rothaeusler K,
- Sospedra M
- 23.↵
- Jelcic I,
- Al Nimer F,
- Wang J, et al.
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- Reindl M,
- Schanda K,
- Woodhall M, et al
- 29.↵
- Krüger J,
- Opfer R,
- Gessert N, et al
- 30.↵
- 31.↵
- Bielekova B,
- Sung MH,
- Kadom N,
- Simon R,
- McFarland H,
- Martin R
- 32.↵
- 33.↵
- Patil VS,
- Madrigal A,
- Schmiedel BJ, et al.
- 34.↵
- Pearce EL,
- Mullen AC,
- Martins GA, et al.
- 35.↵
- 36.↵
- Kaneko K,
- Sato DK,
- Nakashima I, et al.
- 37.↵
- 38.↵
- 39.↵
- Bergman J,
- Dring A,
- Zetterberg H, et al
- 40.↵
- Varrin-Doyer M,
- Shetty A,
- Spencer CM, et al.
- 41.↵
- 42.↵
- Fonseka CY,
- Rao DA,
- Teslovich NC, et al.
- 43.↵
- Bauer ME
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- Soltys K,
- Stuchlikova M,
- Hlavaty T, et al
- 48.↵
- Snir O,
- Widhe M,
- von Spee C, et al
- 49.↵
- Lundberg K,
- Bengtsson C,
- Kharlamova N, et al
- 50.↵
- Boman A,
- Brink M,
- Lundquist A, et al
Letters: Rapid online correspondence
REQUIREMENTS
You must ensure that your Disclosures have been updated within the previous six months. Please go to our Submission Site to add or update your Disclosure information.
Your co-authors must send a completed Publishing Agreement Form to Neurology Staff (not necessary for the lead/corresponding author as the form below will suffice) before you upload your comment.
If you are responding to a comment that was written about an article you originally authored:
You (and co-authors) do not need to fill out forms or check disclosures as author forms are still valid
and apply to letter.
Submission specifications:
- Submissions must be < 200 words with < 5 references. Reference 1 must be the article on which you are commenting.
- Submissions should not have more than 5 authors. (Exception: original author replies can include all original authors of the article)
- Submit only on articles published within 6 months of issue date.
- Do not be redundant. Read any comments already posted on the article prior to submission.
- Submitted comments are subject to editing and editor review prior to posting.
You May Also be Interested in
Dr. Babak Hooshmand and Dr. David Smith
► Watch
Related Articles
- No related articles found.
Topics Discussed
Alert Me
Recommended articles
-
Cytokine Networks in Multiple Sclerosis
The role of autoimmune T lymphocytes in the pathogenesis of multiple sclerosisReinhard Hohlfeld, Edgar Meinl, Frank Weber et al.Neurology, June 01, 1995 -
Article
Antibodies to MOG have a demyelination phenotype and affect oligodendrocyte cytoskeletonRussell C. Dale, Esther M. Tantsis, Vera Merheb et al.Neurology - Neuroimmunology Neuroinflammation, May 22, 2014 -
Article
Fulminant demyelinating encephalomyelitisInsights from antibody studies and neuropathologyFranziska Di Pauli, Romana Höftberger, Markus Reindl et al.Neurology - Neuroimmunology Neuroinflammation, November 04, 2015 -
Article
Distinction between MOG antibody-positive and AQP4 antibody-positive NMO spectrum disordersDouglas Kazutoshi Sato, Dagoberto Callegaro, Marco Aurelio Lana-Peixoto et al.Neurology, January 10, 2014