Intrathecal B Cells in MS Have Significantly Greater Lymphangiogenic Potential Compared to B Cells Derived From Non-MS Subjects

Although B cell depletion is an effective therapy of multiple sclerosis (MS), the pathogenic functions of B cells in MS remain incompletely understood. We asked whether cerebrospinal fluid (CSF) B cells in MS secrete different cytokines than control-subject B cells and whether cytokine secretion affects MS phenotype. We blindly studied CSF B cells after their immortalization by Epstein-Barr Virus (EBV) in prospectively-collected MS patients and control subjects with other inflammatory-(OIND) or non-inflammatory neurological diseases (NIND) and healthy volunteers (HV). The pilot cohort (n = 80) was analyzed using intracellular cytokine staining (n = 101 B cell lines [BCL] derived from 35 out of 80 subjects). We validated differences in cytokine production in newly-generated CSF BCL (n = 207 BCL derived from subsequent 112 prospectively-recruited subjects representing validation cohort), using ELISA enhanced by objective, flow-cytometry-based B cell counting. After unblinding the pilot cohort, the immortalization efficiency was almost 5 times higher in MS patients compared to controls (p < 0.001). MS subjects' BCLs produced significantly more vascular endothelial growth factor (VEGF) compared to control BCLs. Progressive MS patients BCLs produced significantly more tumor necrosis factor (TNF)-α and lymphotoxin (LT)-α than BCL from relapsing-remitting MS (RRMS) patients. In the validation cohort, we observed lower secretion of IL-1β in RRMS patients, compared to all other diagnostic categories. The validation cohort validated enhanced VEGF-C production by BCL from RRMS patients and higher TNF-α and LT-α secretion by BCL from progressive MS. No significant differences among diagnostic categories were observed in secretion of IL-6 or GM-CSF. However, B cell secretion of IL-1β, TNF-α, and GM-CSF correlated significantly with the rate of accumulation of disability measured by MS disease severity scale (MS-DSS). Finally, all three cytokines with increased secretion in different stages of MS (i.e., VEGF-C, TNF-α, and LT-α) enhance lymphangiogenesis, suggesting that intrathecal B cells directly facilitate the formation of tertiary lymphoid follicles, thus compartmentalizing inflammation to the central nervous system.

Introduction

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease of the central nervous system (CNS) thought to be mediated by myelin-specific CD4+ T cells (1, 2). Whereas depletion of CD4+ T cells does not reduce formation of MS lesions (3), depletion of B cells is strikingly more effective (4). Moreover, B cell depletion inhibits disability progression in MS subjects younger than 55 years, including those with progressive MS (PMS) (5, 6).

Indeed, intrathecal activation of B cells and their differentiated, antibody (Ab)-producing counterparts (plasmablasts and plasma cells) is the hallmark of MS (7, 8). Pathologists identified B cells in the meningeal tissue, in the form of diffuse infiltrates in most MS subjects (9) or tertiary lymphoid follicles in some secondary-progressive MS (SPMS) patients (10, 11). The same authors also demonstrated that meningeal follicles harbor Epstein-Barr Virus (EBV) (12, 13), a finding not reproduced by others (14). Nevertheless, EBV remains one of the best documented environmental links in MS (15) with a possible causal role (16, 17).

Despite this wealth of observational data, it remains unclear whether B cell depleting therapies must also deplete intrathecal B cells, in addition to peripheral B cells, to achieve their efficacy (18, 19). Similarly, the exact functions by which intrathecal B cells modulate MS severity or phenotype are unknown.

Analysis of intrathecal B cell functions, whether it be their antigen-presenting ability, the ability to co-stimulate T cells, or produce soluble factors, has been precluded by low numbers of B cells in the cerebrospinal fluid (CSF) (20). With current technologies, it is possible to isolate only few hundred B cells from a single MS patient lumbar puncture. Thus, the phenotypic analyses of CSF B cells performed so far focused on a limited number of surface markers (7, 8, 21, 22). These studies demonstrated a significantly higher frequency of CSF B cells and plasmablasts in patients with MS and other inflammatory neurological diseases (OIND) in comparison to non-inflammatory neurological diseases (NIND) controls (7, 8).

To examine the functional properties of intrathecal B cells, we optimized a procedure for the transformation of CSF B cells by EBV (23) and blindly analyzed cytokine secretion from a large number of intrathecal B cell lines (BCL) derived from pilot and validation cohorts, described in this paper.

Methods

Patients

This study was carried out in accordance with the recommendations of the institutional review board (IRB) at the National Institutes of Health (NIH). The protocol was approved by the NIH Combined Neuroscience IRB, and all subjects gave written informed consent in accordance with the Declaration of Helsinki. We applied the optimized EBV-transformation procedure to prospectively-acquired untreated patients presenting for a diagnostic work-up of CNS neuroimmunological disorders or for screening to NIH PMS clinical trials under the NIH protocol “Comprehensive multimodal analysis of neuroimmunological diseases of the central nervous system” (ClinicalTrials.gov identifier NCT00794352; Table 1; Supplemental Information). Eligible subjects were age 18–75 with clinical and neuroimaging evidence of immune-mediated CNS damage or healthy volunteers (HV) with no CNS disease and no risk factors for immune-mediated systemic or CNS diseases. Patients underwent a diagnostic workup as previously detailed (24), and MS was diagnosed based on 2010 revisions to McDonald diagnostic criteria (25).

TABLE 1

Table 1. Demographic, clinical, and CSF data of pilot and validation cohort patients.

Generation of EBV-Transformed BCL

All prospectively-enrolled subjects were coded, and CSF samples were processed blindly. CSF cells were extracted by centrifugation (300 g × 10 min) within 30 min of CSF collection. EBV immortalization Protocol A was used in the pilot cohort between June 2008 and June 2010 on all subjects who lacked blood contamination (i.e., < 100 RBC/μL of CSF) and had a minimum of 5000 CSF leukocytes available for immortalization. Protocol B was used in the validation cohort collected between July 2010 and February 2014 on subjects without blood contamination who had a minimum of 10,000 CSF leukocytes for immortalization.

Protocol A: Protocol Studying Immortalization Efficacy in the Pilot Cohort

Irradiated (6,000 Rd) CD40+ NIH 3T3 cells were seeded on a 384-well plate (5,000 cells per well) 1 day before scheduled CSF acquisition to achieve their confluency. CSF cells were counted, and co-cultured with CD40+ NIH 3T3 cells at 5,000 CSF cells/well in B cell media (RPMI1640 + 15% FBS + Penicillin/Streptomycin + L-glutamine) containing 2.5 μg/mL CpG oligonucleotides (ODN 2006, InvivoGen, San Diego, USA), 3 ng/mL IFN-γ, 0.5 μg/mL cyclosporine A (Sigma-Aldrich, St. Louis, USA), and 50% EBV supernatant extracted from a confluent B95.8 marmoset leukocyte cell line (CRL-1612, ATCC, Manassas, USA). These optimized conditions for CSF B cell immortalization were determined in extensive preliminary experiments using protocols previously described for immortalization of peripheral B cells (26).

For subjects with larger numbers of CSF cells, proportionally more wells were seeded at 5,000 cells/well.

Growing cell clusters were transferred after 2–3 weeks to progressively larger volumes of fresh B cell media until approximately 10 million immortalized CSF B cells were obtained. These were then cryopreserved. Functional analyses were performed in parallel for all immortalized BCL.

Protocol B: Modification to Increase Immortalization Rate of Non-MS Subjects

Instead of seeding multiple wells at 5,000 CSF cells/well, for validation cohort we seeded only 1–3 wells/subject, with a 10,000–25,000 CSF cells/well of 384-well plate (i.e., increasing the number of CSF cells for subjects with a larger number of available CSF cells).

Immunophenotyping by Fluorescence-Activated Cell Sorting (FACS; Pilot Cohort)

Thawed, overnight-rested BCL were evaluated for surface expression of costimulatory and activation markers. Intracellular cytokine staining (ICCS) was performed as described (27). Briefly, BCL were activated with phorbol 12-myristate 13-acetate (PMA; 32 nM) and Ionomycin (1 μM) in the presence of Brefeldin A (1 μL per 1 × 10⁶ cells, GolgiPlug, BD Biosciences, San Jose, USA) for 12–14 h. After staining for surface markers, BCL were fixed, permeabilized, and stained for intracellular cytokines. Data were acquired by LSR II and analyzed by FACSDiva software (both BD Biosciences).

The following mouse anti-human monoclonal antibodies were used for BCL phenotypization: APC-α-CD19 (HIB19, BD, Franklin Lakes, USA); PE-α-CD5 (UCHT2, BD); FITC-α-CD24 (ML5, BD); APC H7-α-CD25 (MA251, BD); Alexa Fluor 700-α-CD38 (HIT2, eBioscience, San Diego, USA); CD80 (L307.4, BD); PE-α-CD86 (IT2.2, eBioscience); PE-α-IL4 (8D4-8,BD); PE-α-IL6 (MQ2-6A3, BD); PE-α-IL12 p40/p70 (C11.5, BD); PE-α-LT α (3598111, BD); APC-α-VEGF (23410, R&D Systems Inc., Minneapolis, USA); APC-α-TGF-β1 (27232, R&D Systems Inc.); APC-α-IFN-γ (B27, BD); APC-α-TNF-α (MAb11, BD). We also used rat anti-human monoclonal antibodies: PE-α-IL2 (MQ1-17H12, BD) and eFluor 450-α-IL10 (JES3-9D7, eBioscience).

B Cell Activation for Enzyme-Linked Immunosorbent Assay (ELISA) Assays (Validation Cohort)

Cryopreserved CSF BCL were thawed and cultured for 1 week in B cell media: RPMI-1640, 1% penicillin/streptomycin, 1% gentamicin, 1% L-glutamine (all Lonza, Basel, Switzerland), and 15% fetal bovine serum (Gemini Bio-products, West Sacramento, USA). Cells were then resuspended in X-VIVO media (Lonza), counted by hemocytometer (LW Scientific, Lawrenceville, USA), and resuspended at 1 × 10⁶ cells/mL. This aliquot was also used for proportional live-B cell counting using flow cytometry (see below). Remaining cells were seeded in 24 well-plate at 1 × 10⁶ B cells/well either unstimulated (control) or stimulated with PMA (5 nM) and ionomycin (0.5 μM; concentrations optimized in preliminary experiments to maximize cytokine production). After 18 h of culture, the supernatants were centrifuged, aliquoted and frozen at −80°C until ELISA assay.

Proportional Flow Cytometry-Based Counting of B Cells for ELISA Assays (Validation Cohort)

To objectively quantify the number of live B cells that generated supernatants for ELISA, we developed flow-cytometry-based assays (Figure 1). An aliquot of hemocytometer-counted BCL at 1 × 10⁶ cells/mL was resuspended in propidium iodide (PI; 2 μg/mL) to stain dying/dead cells. One million APC Calibration fluorescent particles (SPHEROTM 5.0–5.9 μM; Spherotech, Lake Forest, USA) were added to each BCL to proportionally enumerate BCL between different subjects. One lot of beads was used for the entire study. Serial dilutions of the PI-stained BCL with APC beads mixture were analyzed by flow cytometer to derive linear regression models of live (PI negative) B cells per APC beads, which were utilized for proportional normalization of B cells between different subjects/experiments (Figure 1).

FIGURE 1

Figure 1. Flow-cytometry-based procedure for objective enumeration of live CSF B cells that contributed to supernatants analyzed for cytokine concentrations by ELISA in the independent validation cohort. (A) Coded, in-vitro expanded EBV-immortalized CSF BCL were manually counted and resuspended at 1 × 10⁶ B cells/mL. This B cell dilution was then used to seed plates for cytokine detection, while, simultaneously, 1 × 10⁶ B cells from this aliquot were stained with PI and mixed with 1 × 10⁶ fluorescently (APC)-tagged microbeads (same batch used for the entire validation cohort) for flow cytometry analysis. (B) B cell/microbead mixture was serially diluted three times at 1:3, before their proportional enumeration by flow cytometry. (C) Flow cytometry output quantified microbeads based on APC fluorescence signal and B cells based on size and granularity. Live B cells in cultures were gated as PI-negative, as dying B cells intercalate PI stain into DNA, altering their emission profile. Numbers of live B cells were plotted against APC microbeads for all 3 dilutions to derive patient-specific linear regressions, from which exact number of live B cells seeded and activated in cytokine-secretion assays was calculated, based on known (and equal among all subjects in the validation cohort) number of fluorescent microbeads.

ELISA Assays

Supernatants from stimulated CSF BCL were analyzed for interleukin (IL)-1β, IL-6, IL-10, tumor necrosis factor (TNF)-α, lymphotoxin (LT)-α and granulocyte macrophage colony stimulating factor (GM-CSF) using the V-Plex Meso Scale Discovery platform ELISA (Meso Scale Diagnostics, Rockville, USA) per manufacturer protocols (28). Additionally, we developed vascular endothelial growth factor (VEGF)-A and VEGF-C assays de novo using antibodies from R&D Systems Inc. Table 2 defines each assay, manufacturer, detection limits, and intra-assay variability. Pilot studies of each assay suggested that supernatants from unstimulated (control) B cells do not contain measurable levels of cytokines. Therefore, stimulated BCL conditions were analyzed exclusively.

TABLE 2

Table 2. ELISA assay development criteria.

Statistical Analysis

The diagnostic code (i.e., relapsing-remitting MS [RRMS], SPMS and primary-progressive MS [PPMS], HV, NIND, and OIND) was broken only after the functional data were acquired.

Gender spread between diagnostic categories was evaluated with a chi-square test. Differences in demographic, clinical and CSF data in each diagnostic category were assessed using a Kruskal-Wallis analysis of variance on ranks, followed by Dunn's test with Dunn-Šidák adjustments for multiple comparisons, while differences in EBV transformation rate were assessed with the binomial proportion test. Between-group comparisons in the pilot cohort were assessed using a Kruskal-Wallis analysis of variance on ranks followed by the multiplicity adjusted Dunn's multiple comparisons test. For the validation cohort, we fit a linear mixed model with diagnosis as a fixed effect and subject as a random effect. This allowed us to test for differences between diagnostic groups while accounting for BCLs that were generated from the same CSF sample. We assessed group differences in the mixed model using an analysis of variance (ANOVA) followed by pair-wise comparisons of the least-squares means for the model with a Tukey adjustment for pair-wise multiple comparisons. When no differences between diagnostic categories were observed (verified with one-way t-tests to assess group-wise differences that resulted in p > 0.05), some were merged logically: for the pilot cohort, we merged OIND and NIND into other neurological diseases (OND). In the validation cohort, we merged HV+NIND, and we combined PPMS and SPMS to PMS cohorts. Box-Cox transformation was applied to the biomarker variables with a non-normal distribution. The Shapiro-Wilk test was used to test the normality of the residuals. SAS version 9.4, Graphpad Prism version 7.0b, and R version 3.4.3 were used for the above analyses and p < 0.05 was used as the significance level.

Correlations between cytokines in the pilot cohort were assessed by Pearson correlation coefficients. In the validation cohort, we used Spearman correlation coefficients with a Bonferroni p-value adjustment to assess the correlation between cytokine concentrations and the MS disease severity scale (MS-DSS) (29), a new sensitive measure of MS disease severity.

Results

After unblinding the diagnostic codes, the n = 80 pilot cohort consisted of 47 MS patients (35 RRMS, 9 PPMS and 3 SPMS) and 33 controls (14 OIND and 19 NIND). MS and controls were well matched for demographic data (Table 1). All subjects with available EBV serology were found to be seropositive (Table 1).

EBV Transformation Rate

The average number of CSF cells seeded per patient was similar between MS and controls (39,623 ± 39,207 vs. 39,412 ± 76,718; ns). Likewise, the number of seeded wells per patient was comparable between the two cohorts (3.87 ± 3.02 vs. 2.55 ± 2.63; ns). We generated at least 1 BCL in 28 out of 47 MS patients (59.6%). In contrast, we obtained CSF BCL in only 7 out of 33 OND controls (21.2%; p = 0.004 in comparison to MS). We quantified transformation efficiency as the number of seeded CSF cells required to obtain one BCL. The transformation efficiency in the control group was, on average, 1 BCL per 100,046 seeded CSF cells (1/91,130 in OIND and 1/129,767 in NIND; ns). Transformation efficiency was almost 5 times higher in the MS group, reaching 1 BCL per 21,163 seeded CSF cells (p < 0.001).

The efficiency of immortalization correlated weakly with CSF cell count (Spearman ρ = 0.337, p = 0.00299) and with CSF IgG index (Spearman ρ = 0.391, p = 5 × 10⁻⁴).

General Phenotype of CSF BCL

CSF BCL did not express high levels of CD5 (1.67 ± 7.94%), CD24 (2.7 ± 3.65%) or CD25 (7.91 ± 7.9%). BCL were almost 100% positive for CD38 (90.52 ± 15.45%) and all expressed co-stimulatory molecules CD80 and CD86. We observed a significant correlation between the surface expression of CD24 and CD38 (Pearson r = 0.344, p = 4.26 × 10⁻⁴; Figure 2A), which reflects their synchronous expression in transitional B cells (30).

FIGURE 2

Figure 2. Phenotype of CSF BCL in the pilot cohort detected by ICCS. (A) Relationship between surface expression of CD38 and CD24+ cells. (B) IFN-γ production resides mostly in switched memory B cells based on significant correlation between proportion of CD27+/IgD– cells and BCL production of IFN-γ. In contrast, transitional B cells are the main producers of IL-4 based on a significant correlation between the proportion of IL-4-producing B cells and IgM (C) and CD38 (D) in BCL. IgM and CD38 are co-expressed in transitional B cells. (E) The negative correlation between proportion of IgM-expressing B cells and double negative (CD27–/IgD–) B cells in CSF BCL support the conclusion that double negative B cells are tissue-resident memory B cells, rather than naïve, IgM expressing B cells. (F) Instead, the proportion of naïve B cells (CD27–/IgD+) positively correlates with IgM expression. (G) LT-α and TNF-α are strongly co-expressed, based on their strong, positive correlation in individual BCL. Group differences were assessed using Kruskal-Wallis analysis of variance on ranks followed by an adjusted Dunn's test for multiple comparisons (R), Pearson correlation coefficients; MFI, Mean Fluorescence Intensity.

Out of the 10 cytokines examined, none of the BCL produced detectable levels of IL-2, IL-6, IL-10, IL-12 or transforming growth factor (TGF)-β and only few BCL produced low levels of interferon (IFN)-γ. IFN-γ production resided mostly in class-switched memory B cells (CD27+/IgD–), as we observed significant correlation between proportion of CD27+/IgD– cells in EBV-transformed BCL and proportion of IFN-γ+ cells (Pearson r = 0.403, p = 2.96 × 10⁻⁴; Figure 2B).

The majority of BCL produced detectable levels of IL-4, TNF-α, LT-α, and VEGF. IL-4 secretion resided mostly in naïve and transitional B cells, as proportion of IL-4+ B cells in EBV-transformed BCL correlated moderately with IgM expression (Pearson r = 0.37, p = 1.4 × 10⁻⁴; Figure 2C) and highly with CD38 expression (Pearson r = 0.56, p = 1.41 × 10⁻⁹; Figure 2D). These naïve and transitional B cells represented only a minority (12.8%) of CSF BCL. The remaining cytokines were not preferentially expressed by a distinct B cell phenotype. The majority of CSF B cells (87.2%) were memory B cells, including CD27-/IgD- memory tissue-based B cells (31), expressing CD80 and CD86 cell activation markers, but generally lacking IgM expression (Pearson r = −0.46, p = 1.38 × 10⁻⁶; Figure 2E), which is instead a hallmark of CD27–/IgD+ naive B cells (Pearson r = 0.58, p = 3 × 10⁻¹⁰; Figure 2F). TNF-α and LT-α were strongly co-expressed by the same B cells, as evidenced by strong correlations between them (Pearson r = 0.70, p = 6.033 × 10⁻¹⁶; Figure 2G).

Inter-Group Differences in the Phenotype of CSF BCL

After unblinding the diagnostic codes, we observed that only 13 BCL were derived from control patients (3 from NIND, 10 from OIND), 15 BCL were from patients with PPMS and the remaining 73 BCL were from patients with RRMS.

BCL from all MS patients had a significantly higher expression of CD80 in comparison to BCL from controls (Figure 3A). No other differences in surface markers reached statistical significance.

FIGURE 3

Figure 3. Differences in the phenotype of CSF BCL in the pilot cohort among diagnostic categories. (A) BCL from RRMS (n = 73) and all MS patients express higher CD80 than BCL from controls (OIND n = 10; NIND n = 3). (B) BCL from MS (n = 88) patients produce more VEGF than BCL from controls. (C,D) BCL from PPMS (n = 15) patients produce more TNF-α and LT-α than BCL from RRMS. Red lines represent means and positive and negative standard deviation from the mean for each diagnostic group. *p < 0.05, **p < 0.01; MFI, mean fluorescence intensity).

Similarly, the only significant difference between MS patients and controls resided in the higher production of VEGF in BCL derived from MS patients (Figure 3B). Finally, BCL derived from PPMS produced significantly higher amounts of TNF-α and LT-α as compared to BCL from RRMS patients (Figures 3C,D).

Validation Cohort

The limitation of the pilot cohort was significantly lower EBV transformation efficiency of the control samples, which we did not anticipate. This led to small numbers of control BCL. Similarly, the numbers of BCL derived from PMS subjects were low, which raised a possibility that the statistically significant increases in the secretion of TNF-α and LT-α we have observed in this cohort may have been over-estimated (32). Therefore, we modified EBV immortalization method to enhance its efficacy by seeding 10,000–25,000 CSF cells into a single well (Protocol B, see section Methods) for the independent validation cohort. This precluded us from validating the differences in EBV immortalization efficiency, but this observation from the pilot cohort had a very low p-value and therefore low Type II error. Instead, the goal of the validation cohort was to transform sufficient number of BCL from MS controls and PMS patients to disprove or validate differences in VEGF, TNF-α, and LT-α secretion among MS and controls and the relationship of these cytokines to the MS phenotype.

We also decided to use ELISA instead of ICCS we employed in the pilot cohort because a single ICCS protocol cannot be optimized for all cytokines simultaneously, as some cytokines are secreted early after activation while others are secreted late and they have different sensitivity to protein transport inhibitors such as monensin and brefeldin A. We hypothesized that our inability to detect some cytokines commonly secreted by B cells (such as IL-10 or IL-6) in the pilot cohort was a consequence of this technical problem of ICCS. ELISA eliminates this problem, but it poses another: the results are highly dependent on the exact number of seeded cells. Thus, instead of relying solely on manual cell counts, we implemented an unbiased flow-cytometry-based normalization of detected cytokine concentrations to live B cell numbers (Figure 1 and section Methods) among all tested subjects. In addition to studying VEGF (A and C isoforms were studied separately) and TNF-α, LT-α, IL-6, and IL-10 to refute or validate findings from the pilot cohort, we also analyzed de novo two additional cytokines with putative role in MS, IL-1β, and GM-CSF (33).

Cytokine Secretion of CSF B Cells in the Independent Validation Cohort

Table 3 describes the observed cytokine secretion results after unblinding diagnostic groups in the validation cohort (Figure 4G).

TABLE 3

Table 3. Phenotype of validation group CSF BCL: secreted cytokine expression.

FIGURE 4

Figure 4. Differences in CSF BCL cytokine secretion among diagnostic categories in the independent validation cohort. (A–E) Depict concentrations of cytokines normalized to 1 million of live seeded B cells, while (F,G) depict ratios of pro-inflammatory cytokines to immunoregulatory IL-10 that reached statistical significance after adjustment for multiple comparisons. Superscripts in the y-axis correspond to lambda values of Box-Cox statistical normalization test (A–E), while ln signifies natural logarithm of cytokine concentration values (F,G). Dots represent individual BCL concentration (or ratio of cytokine concentration). Red lines indicate mean with positive and negative standard deviation for each diagnostic group. Group differences were assessed using an ANOVA on a linear mixed model with subject specified as a random effect. *p < 0.05, **p < 0.01 (after adjustments for multiple comparisons).

First, we validated our hypothesis that undetectable levels of IL-10 and IL-6 in the pilot cohort was consequence of the selected ICCS procedure, as both cytokines were readily detectable by ELISA.

To assure that no bias was introduced by analyzing more than 1 CSF BCL per subject, we added a random effects term for each subject to account for variance that is due to repeated measures (i.e., multiple BCL) on a single subject. When we normalized cytokine production per 1 million of live, EBV-transformed CSF B cells we observed lower secretion of IL-10 in BCL derived from RRMS patients compared to non-inflammatory controls (difference between means [DM] = 164.7 pg/mL), but it did not reach statistical significance after adjusting for multiple comparisons (Figure 4A). Similarly, we saw no difference among diagnostic categories in the concentrations of IL-6 (Table 3).

From cytokines that were not measured in the pilot cohort, we found no significant differences among diagnostic categories in the B-cell secretion of GM-CSF (Table 3). Surprisingly, we observed lower IL-1β production in B cells derived from RRMS patients compared to all diagnostic categories (Figure 4B; OIND DM = 5.7 pg/mL, p = 0.0434; SPMS DM = 5.3, p = 0.0377; PPMS DM = 6.1 p = 0.0111; PMS DM = 5.6 pg/mL, p = 0.0028), as well as HV+NIND (DM = 14.4, p = 0.0031).

The validation cohort data revealed that increased VEGF secretion in MS B cells observed in the pilot cohort was limited to VEGF-C (Figure 4C). Furthermore, increased VEGF-C secretion was only observed in RRMS patients against all other diagnostic categories, even though only comparison with OIND reached statistical significance (DM = −121.2 pg/mL, p = 0.0422).

The validation cohort again mirrored the pilot cohort ICCS data with B cells derived from PMS patients secreting highest concentrations of TNF-α in comparison to RRMS, OIND, and non-inflammatory controls, although only comparisons to RRMS (Figure 4D; DM = −503.1 pg/mL, p = 0.0122) and non-inflammatory controls (Figure 4D; DM = −117.0 pg/mL, p = 0.0217) reached formal statistical significance. Similar data were obtained for LT-α, where the remaining significant difference was observed between PMS and RRMS (Figure 4E; DM = −389.0, p = 0.023).

Finally, because ratios of pro-inflammatory cytokines such as TNF-α and LT-α to immunoregulatory cytokine IL-10 were shown to be increased in peripheral blood of MS patients as compared to HV (34), we also generated analogous ratios for all pro-inflammatory cytokines we tested in the validation cohort. We observed a trend for lower ratios of pro-inflammatory cytokines to IL-10 in HV+NIND group, in comparison to subjects with intrathecal inflammation (i.e., OIND + RRMS + PMS). However, only the LT-α/IL-10 ratio achieved formal statistical significance with HV+NIND when compared to PPMS (Figure 4F; DM = −51.1, p = 0.0135) and PMS (Figure 4F; DM = −91.9, p = 0.0292). VEGF-C/IL-10 was higher in the RRMS cohort compared to OIND (Figure 4G; DM = −183.7, p = 0.0210).

Exploring Pathogenic Role of CSF B Cell-Secreted Cytokines

As demonstrated by the correlation matrix (Figure 5), the validation cohort validated significant positive correlation between TNF-α and LT-α. TNF-α also strongly correlated with two additional cytokines with putative pathogenic role in MS, tested only in the validation cohort, GM-CSF and IL-1β In contrast, VEGF-C secretion was only positively correlated with VEGF-A secretion.

FIGURE 5

Figure 5. Speaman correlation matrix of validation group CSF BCL cytokine secretion and correlation of cytokines against MS disease severity scale (MS-DSS). Cytokines analyzed in the validation cohort are compared to each other and to the MS-DSS, which measures rates of disability progression. Red colors indicate positive, direct correlation, while blue colors indicate negative, inverse correlation, according to the continuous heatmap legend depicted at the right side of the correlation matrix. The Spearman correlation coefficient for each pair of cytokines is shown at the intersecting x- and y-axis label within each circle. Outlined circles represent significant correlations (p < 0.05 after adjustments for multiple comparisons) and increasing diameter of circles indicate a higher level of significance, according to the gray legend depicted below the correlation matrix.

To explore whether any of the differences identified between MS subjects and controls participate in CNS tissue destruction, we assessed correlation between patient-specific CSF B cells cytokine production measured by ELISA (computed as average concentration of specific cytokine derived from all BCL isolated from specific subject) and new, sensitive model of MS severity, MS-DSS (29). MS-DSS is derived from statistical learning. In contrast to an older measure of MS severity, the Expanded Disability Status Scale (EDSS)-based MS Severity Score (35), MS-DSS can predict future rates of disability progression. In MS-DSS, the disability measurement by insensitive, discrete EDSS is replaced by a much more sensitive, statistical learning optimized continuous disability scale, Combinatorial Weight-Adjusted Disability Scale [CombiWISE (36)]. Additionally, disability ratings in MS-DSS are adjusted for patient-specific treatments using mathematical models derived from previously published meta-analysis of randomized clinical trials of disease-modifying treatments of MS (37). Finally, the MS-DSS model also contains certain demographic data and a Composite MRI Scale of CNS tissue destruction [COMRIS-CTD (38)]. MS-DSS is freely available via a web-based interphase at: https://bielekovalab.shinyapps.io/msdss/.

We observed mild, but highly significant correlations between MS-DSS and three correlated cytokines, IL-1β, TNF-α, and GM-CSF (Figure 5).

Discussion

In this study, we present the first functional analysis (i.e., secretion of multiple cytokines) of intrathecal B cells in MS. Demonstration that related B cell clones populate all CNS compartments in MS, including CSF (39), validates the evaluation of CSF B cells as a representative population. The strength of our study resides in the application of highly standardized protocols in a blinded fashion to a sizable number of MS patients and controls prospectively collected at a single institution. Critically, we also validated most important new findings in prospectively acquired, independent validation cohort. We consider application of different methodology in the validation cohort (i.e., ELISA as opposed to ICCS) also a strength of this study, because it demonstrates that the data are independent of the methodological constraints of applied assays. Our initial hesitancy to use ELISA assay for answering question about differences in cellular phenotypes stem from the fact that in contrast to ICCS, which outputs values proportional to the concentration of cytokines in each cell, ELISA measures sum of the cytokines secreted by all living cells in the culture. Therefore, ELISA results are crucially dependent on the number of seeded cells and manual cell counts may introduce noise in the measurement. We successfully mitigated this problem by developing flow-cytometry based, objective normalization of imputed cell numbers among all studied subjects. This normalization was performed before unblinding the diagnostic categories.

The natural weakness of our study resides in the inability to measure cytokine production in neither CNS B cells nor CSF B cells and relying on an EBV transformation and in vitro expansion step. Nevertheless, CSF biomarkers have been established as viable proxies for inflammation originating from the CNS (40) and correlate with disease pathologies in MS (41, 42). As we discussed before, analogous clonotypes of T and B cells were identified in CNS tissue and CSF of individual patients. In regard to EBV transformation, due to limited numbers of B cells in the human CSF irrespective of diagnostic categories studied (with exception of intrathecal B cell lymphoma, which was not studied here), no current technologies can provide functional information from ex vivo isolated CSF B cells. We cannot rule out a possibility that EBV transformation might have affected phenotype of derived BCL (43, 44). However, blinding and written SOPs assured that any bias that EBV transformation and in vitro manipulation may have introduced affected all subjects equally. Therefore, this cannot explain the observed (and validated) differences among diagnostic categories. The obtained results to a large degree mitigate aforementioned concern, because our study reproduced a key finding from previously published ex vivo analyses of primary CSF B cells (7, 8), such as significantly higher levels of the co-stimulatory molecule CD80 (but not CD86) in CSF BCL derived from MS patients in comparison to controls (22) and trend for lower production of IL-10 observed in peripheral blood B cells derived from MS patients in comparison to HV (34).

There are important novel findings identified in the current study. First, we observed an almost five-fold higher efficiency of EBV transformation in MS patients compared to controls. While this may be partially driven by higher CSF pleiocytosis and higher proportion of CSF B cells in MS as compared to NIND patients and HV, this cannot be true for the OIND cohort, because OIND and MS subjects have similar proportions of CSF B cells (7, 8, 24). Furthermore, we observed only a weak correlation between CSF pleiocytosis or CSF IgG index and transformation efficiency, indicating that factors other than a higher number/proportion of CSF B cells in the MS cohort must underlie the differences in transformation efficiency. Higher transformation efficiency cannot be caused by higher EBV seropositivity (15), as we observed that all subjects on whom EBV serology data are available were seropositive (Table 1). Thus, higher transformation efficiency in MS could be mechanistically linked to higher intrathecal EBV expression (12) previously reported in MS, or may result from a defect of cytolytic activity of CSF T cells toward EBV transformed B cells, as described for peripheral T cells in MS (45). We consider the latter possibility less likely because use of cyclosporine A in our transformation protocol should limit the effect of cytotoxic T cells. The final possibility is that MS susceptibility genes enhance EBV infection of B cells and the rates of their immortalization. This hypothesis is consistent with both experimental data presented here and epidemiological data linking EBV to MS.

The second group of novel findings relate to MS-specific differences in cytokine phenotypes of intrathecal B cells. The first finding is rather surprising, the decrease in secretion of pro-inflammatory cytokine IL-1β in the CSF B cells from RRMS subjects in comparison to all other diagnostic subgroups. This finding was observed only in the validation cohort because we have not studied IL-1β secretion by ICCS in the pilot cohort. However, the level of statistical significance and congruency of findings against all control groups make this finding rather robust, even in the absence of independent validation. Increased IL-1β secretion is a hallmark of genetic auto-inflammatory syndromes associated with the activation of the inflammasome, predominantly in the cells of innate immunity. IL-1β secretion is regulated by caspases, especially by caspase-1, but IL-1β secretion is unconventional in that it does not enter the classical endoplasmic reticulum (ER)-Golgi protein secretion pathway. Instead, its secretion requires formation of a new intracellular compartment, called Unconventional Protein Secretion (UPS) (46), induced by stress such as the deprivation of nutrients (47). Additionally, IL-1β secretion is linked to autophagy, raising a possibility that the unique cytokine secretion phenotype of RRMS B cells may be linked to dysfunction in autophagy pathways, previously identified in immune cells of MS and its animal model Experimental Autoimmune Encephalomyelitis (EAE) (48). However, we must not forget that ELISA measures the difference between secretion and consumption and that IL-1β is also consumed by B cells, especially of the regulatory phenotype (49). Therefore, enhanced consumption of IL-1β by B cells derived from RRMS patients, due to their higher proportion of regulatory B cells in the intrathecal compartment represents an emerging new hypothesis, testable in future studies. Intriguingly, IL-1β production was highly correlated with the secretion of TNF-α and LT-α and these cytokines were not only expressed at higher levels in PMS but also correlated with MS severity measured by MS-DSS. In other words, patients with higher B cell secretion of IL-1β, TNF-α (and GM-CSF) have more aggressive MS, characterized by faster accumulation of disability. Due to collinearity between these three cytokines, it is unclear whether each has an independent contribution to MS severity. We tried to answer this issue by exploring statistical learning using either generalized boosting machine or multiple linear regressions to assess whether a complex model that contains all three (or two) of these cytokines correlates better with MS-DSS than individual cytokines. None of these models performed better than IL-1β alone (data not shown), suggesting that IL-1β secretion has the dominant effect on MS severity. While correlation is not causation, our results generate the hypothesis for future studies that B cell secretion of IL-1β and/or their decreased consumption of IL-1β due to low proportion of regulatory B cells in the intrathecal compartment contributes to CNS tissue destruction in MS.

The remaining differences in the phenotype of MS intrathecal B cells were observed in both pilot and validation cohorts. We found and validated the enhanced production of VEGF, specifically the VEGF-C (but not VEGF-A) family member, in RRMS patients. VEGF-A binds to VEGF receptors 1 (VEGFR1; Flt-1) and−2 (VEGFR2; KDR/Flk-1), expressed predominantly in endothelial cells (50). Not surprisingly, then, VEGF-A mediates the effects of VEGF on endothelial cells, stimulating angiogenesis, but also opening epithelial tight junctions, leading to vascular edema. In contrast, VEGF-C binds a third receptor (VEGFR3; Flt4), expressed on the cells of the lymphatic endothelium, and thus mediates lymphangiogenesis (50). At this point, we can only speculate about functional consequences of enhanced VEGF-C secretion by RRMS B cells; it is possible that the lymphangiogenesis initiates the formation of lymphoid aggregates, also called tertiary lymphoid follicles, found in the meninges of patients with long-standing MS and being linked to more severe CNS tissue destruction (11).

We also found (and validated) that, compared to RRMS patients' BCL, intrathecal B cells derived from PMS patients produced higher amounts of two closely related cytokines, LT-α and TNF-α. In view of the VEGF-C findings described above, it is important to note that both LT-α and TNF-α play an important role in the development and maintenance of tertiary lymphoid follicles (10, 11). Thus, the major conclusion we can derive from our analysis of the functional phenotype of CSF B cells is that B cells derived from MS patients differ from B cells derived from non-MS controls, including those with inflammatory neurological diseases in their pro-lymphangiogenic potential. Why VEGF-C is upregulated only in B cells derived from the early, RRMS stage, and why B cells from PMS instead secrete higher concentrations of TNF-α and LT-α are intriguing topics for future investigation. Similar MS-stage-specific differences were observed in peripheral blood for a subset of MS susceptibility genes, using expression profiling (52). Specifically, previously reported MS susceptibility gene MYC (51) was upregulated in peripheral blood only in RRMS, but not PMS patients (52). Even further, mice overexpressing MYC in B-cell progenitors had increased VEGF production and concurrent stimulation of angiogenesis and lymphangiogenesis (53). It is possible that MYC deregulation plays a causal role in the observed VEGF-C phenotype in RRMS, but not progressive patients. In contrast, PMS groups (but not RRMS patients) had significantly lower expression of TNFRSF1A (52) in blood leukocytes, suggesting that this is either a consequence of higher TNF-α production in progressive MS (as part of negative feedback loop), or, alternatively, lower expression of the receptor may cause lower consumption of TNF-α by immune cells derived from progressive MS. Clearly, expression or genes and translation of proteins can change with disease evolution, underlying processes such as terminal differentiation of memory B cells induced by chronic inflammation, as has been previously described for T cells in PMS but not RRMS (54). Alternatively, the observed changes may have a causative role: e.g., patients with increased production of TNF-α and LT-α may be predisposed to progressive disease. These two possibilities can be only addressed in longitudinal studies that will assess the rate of MS progression in RRMS subjects whose intrathecal B cells produce high or low levels of TNF-α and LT-α.

In conclusion, the observed results suggest that the defining role intrathecal B cells play in MS (in contrast to other CNS inflammatory diseases) is to facilitate compartmentalization of the inflammation into CNS tissue by promoting lymphangiogenesis. The role of EBV infection in this process remains unknown.

Data Availability Statement

All data generated or analyzed during this study are included in this article and its Supplementary Information files.

Author Contributions

BB conceived and supervised the study. QX conducted the ICCS and JS conducted the ELISA experiments. QX, ER, and SW optimized the conditions for the B cell immortalization and immortalized the cells, cryopreserved, and thawed CSF B cells. JS and ER optimized the PMA/IoM stimulation of B cells for the validation cohort. PK developed and maintained the database. KJ, QX, JS, TW, and BB analyzed the results. QX, JS, SW, KJ, and BB drafted the paper. All authors reviewed the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The study was supported by the intramural research program of the National Institute of Neurological Disorders and Stroke (NINDS) and the Clinical Center/ US National Institutes of Health (NIH).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2018.00554/full#supplementary-material

References

1. Noseworthy JH, Lucchinetti C, Rodriguez M, Weinshenker BG. Multiple sclerosis. N Engl J Med. (2000) 343:938–52. doi: 10.1056/NEJM200009283431307

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sospedra M, Martin R. Immunology of multiple sclerosis. Annu Rev Immunol. (2005) 23:683–747. doi: 10.1146/annurev.immunol.23.021704.115707

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wiendl H, Hohlfeld R. Therapeutic approaches in multiple sclerosis: lessons from failed and interrupted treatment trials. BioDrugs (2002) 16:183–200. doi: 10.2165/00063030-200216030-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bielekova B, Becker B. Monoclonal antibodies in MS: mechanism of action. Neurology (2010) 74:S31–40. doi: 10.1212/WNL.0b013e3181c97ed3

CrossRef Full Text | Google Scholar

5. Montalban X, Hauser SL, Kappos L, Arnold DL, Bar-Or A, Comi G, et al. Ocrelizumab versus placebo in primary progressive multiple sclerosis. N Engl J Med (2016). 376:209–20. doi: 10.1056/NEJMoa1606468

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Hauser SL, Bar-Or A, Comi G, Giovannoni G, Hartung HP, Hemmer B, et al. Ocrelizumab versus interferon beta-1a in relapsing multiple sclerosis. N Engl J Med. (2017) 376:221–34. doi: 10.1056/NEJMoa1601277

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Cepok S, Rosche B, Grummel V, Vogel F, Zhou D, Sayn J, et al. Short-lived plasma blasts are the main B cell effector subset during the course of multiple sclerosis. Brain (2005) 128:1667–76. doi: 10.1093/brain/awh486

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Corcione A, Casazza S, Ferretti E, Giunti D, Zappia E, Pistorio A, et al. Recapitulation of B cell differentiation in the central nervous system of patients with multiple sclerosis. Proc Natl Acad Sci USA (2004) 101:11064–9. doi: 10.1073/pnas.0402455101

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Frischer JM, Bramow S, Dal-Bianco A, Lucchinetti CF, Rauschka H, Schmidbauer M, et al. The relation between inflammation and neurodegeneration in multiple sclerosis brains. Brain (2009) 132:1175–89. doi: 10.1093/brain/awp070

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Serafini B, Rosicarelli B, Magliozzi R, Stigliano E, Aloisi F. Detection of ectopic B-cell follicles with germinal centers in the meninges of patients with secondary progressive multiple sclerosis. Brain Pathol (2004) 14:164–74. doi: 10.1111/j.1750-3639.2004.tb00049.x

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Magliozzi R, Howell O, Vora A, Serafini B, Nicholas R, Puopolo M, et al. Meningeal B-cell follicles in secondary progressive multiple sclerosis associate with early onset of disease and severe cortical pathology. Brain (2007) 130:1089–104. doi: 10.1093/brain/awm038

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Serafini B, Rosicarelli B, Franciotta B, Magliozzi R, Reynolds R, Cinque P, et al. Dysregulated Epstein-Barr virus infection in the multiple sclerosis brain. J Exp Med. (2007) 104:2899–912. doi: 10.1084/jem.20071030

CrossRef Full Text | Google Scholar

13. Serafini B, Severa M, Columba-Cabezas S, Rosicarelli B, Veroni C, Chiappetta G, et al. Epstein-Barr virus latent infection and BAFF expression in B cells in the multiple sclerosis brain: implications for viral persistence and intrathecal B-cell activation. J Neuropathol Exp Neurol. (2010) 69:677–93. doi: 10.1097/NEN.0b013e3181e332ec

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Willis SN, Stadelmann C, Rodig SJ, Caron T, Gattenloehner S, Mallozzi SS, et al. Epstein-Barr virus infection is not a characteristic feature of multiple sclerosis brain. Brain (2009) 132(Pt 12):3318–28. doi: 10.1093/brain/awp200

CrossRef Full Text | Google Scholar

15. Ascherio A, Munger KL. Environmental risk factors for multiple sclerosis. Part I: the role of infection. Ann Neurol. (2007) 61:288–99. doi: 10.1002/ana.21117

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Endriz J, Ho PP, Steinman L. Time correlation between mononucleosis and initial symptoms of MS. Neurol Neuroimmunol Neuroinflam. (2017) 4:e308. doi: 10.1212/NXI.0000000000000308

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Levin LI, Munger KL, O'Reilly EJ, Falk KI, Ascherio A. Primary infection with the Epstein-Barr virus and risk of multiple sclerosis. Ann Neurol. (2010) 67:824–30. doi: 10.1002/ana.21978

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Cross AH, Stark JL, Lauber J, Ramsbottom MJ, Lyons JA. Rituximab reduces B cells and T cells in cerebrospinal fluid of multiple sclerosis patients. J Neuroimmunol. (2006) 180:63–70. doi: 10.1016/j.jneuroim.2006.06.029

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Komori M, Lin YC, Cortese I, Blake A, Ohayon J, Cherup J, et al. Insufficient disease inhibition by intrathecal rituximab in progressive multiple sclerosis. Ann Clin Transl Neurol. (2016) 3:166–79. doi: 10.1002/acn3.293

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Cepok S, Jacobsen M, Schock S, Omer B, Jaekel S, Boddeker I, et al. Patterns of cerebrospinal fluid pathology correlate with disease progression in multiple sclerosis. Brain (2001) 124:2169–76. doi: 10.1093/brain/124.11.2169

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Sellebjerg F, Jensen J, Jensen CV, Wiik A. Expansion of CD5 - B cells in multiple sclerosis correlates with CD80 (B7-1) expression. Scand J Immunol. (2002) 56:101–7. doi: 10.1046/j.1365-3083.2002.01108.x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Sellebjerg F, Jensen J, Ryder LP. Costimulatory CD80 (B7-1) and CD86 (B7-2) on cerebrospinal fluid cells in multiple sclerosis. J Neuroimmunol. (1998) 84:179–87. doi: 10.1016/S0165-5728(97)00261-0

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Traggiai E, Becker S, Subbarao K, Kolesnikova L, Uematsu Y, Gismondo MR, et al. An efficient method to make human monoclonal antibodies from memory B cells: potent neutralization of SARS coronavirus. Nat Med. (2004) 10:871–5. doi: 10.1038/nm1080

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Han S, Lin YC, Wu T, Salgado AD, Mexhitaj I, Wuest SC, et al. Comprehensive immunophenotyping of cerebrospinal fluid cells in patients with neuroimmunological diseases. J Immunol. (2014) 192:2551–63. doi: 10.4049/jimmunol.1302884

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Polman CH, Reingold SC, Banwell B, Clanet M, Cohen JA, Filippi M, et al. Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald criteria. Ann Neurol. (2011) 69:292–302. doi: 10.1002/ana.22366

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lanzavecchia A, Bernasconi N, Traggiai E, Ruprecht CR, Corti D, Sallusto F. Understanding and making use of human memory B cells. Immunol Rev. (2006) 211:303–9. doi: 10.1111/j.0105-2896.2006.00403.x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bielekova B, Catalfamo M, Reichert-Scrivner S, Packer A, Cerna M, Waldmann TA, et al. Regulatory CD56bright natural killer cells mediate immunomodulatory effects of IL-2R-alpha-targeted therapy (daclizumab) in multiple sclerosis. Proc Natl Acad Sci USA (2006) 103:5941–6. doi: 10.1073/pnas.0601335103

CrossRef Full Text | Google Scholar

28. Komori M, Blake A, Greenwood M, Lin YC, Kosa P, Ghazali D, et al. CSF markers reveal intrathecal inflammation in progressive multiple sclerosis. Ann Neurol. (2015) 78:3–20. doi: 10.1002/ana.24408

CrossRef Full Text | Google Scholar

29. Weideman AM, Barbour C, Tapia-Maltos MA, Tran T, Jackson K, Kosa P, et al. New multiple sclerosis disease severity scale predicts future accumulation of disability. Front Neurol. (2017) 8:598. doi: 10.3389/fneur.2017.00598

CrossRef Full Text | Google Scholar

30. Anolik JH, Looney RJ, Lund FE, Randall TD, Sanz I. Insights into the heterogeneity of human B cells: diverse functions, roles in autoimmunity, and use as therapeutic targets. Immunol Res. (2009) 45:144–58. doi: 10.1007/s12026-009-8096-7

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Ehrhardt GR, Hsu JT, Gartland L, Leu CM, Zhang S, Davis RS, et al. Expression of the immunoregulatory molecule FcRH4 defines a distinctive tissue-based population of memory B cells. J Exp Med. (2005) 202:783–91. doi: 10.1084/jem.20050879

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Button KS, Ioannidis JP, Mokrysz C, Nosek BA, Flint J, Robinson ES, et al. Power failure: why small sample size undermines the reliability of neuroscience. Nat Rev Neurosci. (2013) 14:365–76. doi: 10.1038/nrn3475

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Li R, Rezk A, Miyazaki Y, Hilgenberg E, Touil H, Shen P, et al. Proinflammatory GM-CSF-producing B cells in multiple sclerosis and B cell depletion therapy. Sci Transl Med. (2015) 7:310ra166. doi: 10.1126/scitranslmed.aab4176

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Duddy ME, Alter A, Bar-Or A. Distinct profiles of human B cell effector cytokines: a role in immune regulation? J Immunol. (2004) 172:3422–7. doi: 10.4049/jimmunol.172.6.3422

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Roxburgh RH, Seaman SR, Masterman T, Hensiek AE, Sawcer SJ, Vukusic S, et al. Multiple Sclerosis Severity Score: using disability and disease duration to rate disease severity. Neurology (2005) 64:1144–51. doi: 10.1212/01.WNL.0000156155.19270.F8

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kosa P, Ghazali D, Tanigawa M, Barbour C, Cortese I, Kelley W, et al. Development of a sensitive outcome for economical drug screening for progressive multiple sclerosis treatment. Front Neurol. (2016) 7:131. doi: 10.3389/fneur.2016.00131

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Weideman AM, Tapia-Maltos MA, Johnson K, Greenwood M, Bielekova B. Meta-analysis of the age-dependent efficacy of multiple sclerosis treatments. Front Neurol. (2017) 8:577. doi: 10.3389/fneur.2017.00577

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Kosa P, Komori M, Waters R, Wu T, Cortese I, Ohayon J, et al. Novel composite MRI scale correlates highly with disability in multiple sclerosis patients. Multiple Sclerosis Relat Disord. (2015) 4:526–35. doi: 10.1016/j.msard.2015.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Lovato L, Willis SN, Rodig SJ, Caron T, Almendinger SE, Howell OW, et al. Related B cell clones populate the meninges and parenchyma of patients with multiple sclerosis. Brain (2011) 134(Pt 2):534–41. doi: 10.1093/brain/awq350

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Renno T, Lin JY, Piccirillo C, Antel J, Owens T. Cytokine production by cells in cerebrospinalfluid during experimental allergic encephalomyelitis in SJL/J mice. J Neuroimmunol. (1994) 49:1–7. doi: 10.1016/0165-5728(94)90174-0

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Carrieri PB, Provitera V, De Rosa T, Tartaglia G, Gorga F Perrella O. Profile of cerebrospinal fluid and serum cytokines in patients with relapsing-remitting multiple sclerosis: a correlation with clinical activity. Immunopharmacol Immunotoxicol. (1998) 20:373–82. doi: 10.3109/08923979809034820

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Besinis DP, Anagnostouli MC. Cerebrospinal fluid markers in neuroinflammations: the paradigm of optic neuritis. Methods Pharmacol Toxicol. (2017) 1:411–31. doi: 10.1007/978-1-4939-6843-5_14

CrossRef Full Text | Google Scholar

43. Pender MP. The essential role of Epstein-barr virus in the pathogenesis of multiple sclerosis. Neuroscientist (2011) 17:351–67. doi: 10.1177/1073858410381531

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Tracy SI, Kakalacheva K, Lunemann JD, Luzuriaga K, Middeldorp J, Thorley-Lawson DA. Persistence of Epstein-barr birus in self-reactive memory B cells. J Virol. (2012) 86:12330–40. doi: 10.1128/JVI.01699-12

CrossRef Full Text | Google Scholar

45. Pender MP, Csurhes PA, Lenarczyk A, Pfluger CM, Burrows SR. Decreased T cell reactivity to Epstein-Barr virus infected lymphoblastoid cell lines in multiple sclerosis. J Neurol Neurosurg Psychiatry (2009) 80:498–505. doi: 10.1136/jnnp.2008.161018

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Malhotra V. Unconventional protein secretion: an evolving mechanism. EMBO J. (2013) 32:1660–4. doi: 10.1038/emboj.2013.104

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Rabouille C. Pathways of unconventional protein secretion. Trends Cell Biol. (2017) 27:230–40. doi: 10.1016/j.tcb.2016.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Liang P, Le W. Role of autophagy in the pathogenesis of multiple sclerosis. Neurosci Bull. (2015) 31:435–44. doi: 10.1007/s12264-015-1545-5

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Rosser EC, Oleinika K, Tonon S, Doyle R, Bosma A, Carter NA, et al. Regulatory B cells are induced by gut microbiota-driven interleukin-1beta and interleukin-6 production. Nat Med. (2014) 20:1334–9. doi: 10.1038/nm.3680

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Simons M, Gordon E, Claesson-Welsh L. Mechanisms and regulation of endothelial VEGF receptor signalling. Nat Rev Mol Cell Biol. (2016) 17:611–25. doi: 10.1038/nrm.2016.87

PubMed Abstract | CrossRef Full Text | Google Scholar

51. The International Multiple Sclerosis Genetics Consortium & the Wellcome Trust Case Control Consortium 2. Genetic Risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature (2011) 476:214–9. doi: 10.1038/nature10251

CrossRef Full Text | Google Scholar

52. Srinivasan S, Di Dario M, Russo A, Menon R, Brini E, Romeo M, et al. Dysregulation of MS Risk genes and pathways at distinct stages of disease. Neurol Neuroimmunol Neuroinflamm. (2017) 4:e337. doi: 10.1212/NXI.0000000000000337

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Ruddell A, Mezquita P, Brandvold KA, Farr A, Iritani BM. B lymphocyte-specific c-Myc expression stimulates early and functional expansion of the vasculature and lymphatics during lymphangiogenesis. Am J Pathol. (2003) 163:2233–44. doi: 10.1016/S0002-9440(10)63581-X

CrossRef Full Text | Google Scholar

54. Wuest SC, Mexhitaj I, Chai NR, Romm E, Scheffel J, Xu B, et al. A complex role of herpes viruses in the disease process of multiple sclerosis. PLoS ONE (2014) 9:e105434. doi: 10.1371/journal.pone.0105434

PubMed Abstract | CrossRef Full Text | Google Scholar