{'fraction_non_alphanumeric': 0.02689486552567237, 'fraction_numerical': 0.015252066596809872, 'mean_word_length': 5.348484848484849, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Immunologic profiles of multiple sclerosis treatments reveal shared early B cell alterations

Objective: We undertook a systems immunology approach of the adaptive immune system in multiple sclerosis (MS), overcoming tradeoffs between scale and level of detail, in order to identify the immunologic signature of MS and the changes wrought by current immunomodulatory treatments. Methods: We developed a comprehensive flow cytometry platform measuring 38 immunologic cell types in the peripheral blood of 245 individuals in a routine clinical setting. These include patients with MS, untreated or receiving any of 4 current immunomodulatory treatments (interferon-β, glatiramer acetate, natalizumab, or fingolimod), patients with autoimmune thyroid disease, and healthy controls. Results: An increase in memory CD8+ T cells and B cells was observed in untreated patients with MS. Interferon-β and fingolimod induce significant changes upon multiple aspects of the peripheral immune system, with an unexpectedly prominent alteration of B cells. Overall, both treatments push the immune system in different directions, with only 2 significant effects shared across these treatments—an increase in transitional B cells and a decrease in class-switched B cells. We further identified heightened B cell-activating factor (BAFF) levels as regulating this shared B cell pathway. Conclusions: A systems immunology approach established different immunologic profiles induced by current immunomodulatory MS treatments, offering perspectives for personalized medicine. Pathways shared between the immunologic architecture of existing efficacious treatments identify targets for future treatment design.

In-depth analysis of variation in the adaptive immune system in humans has proven to be difficult due to trade-offs between scale and level of detail. Large-scale flow cytometry screens have now been able to characterize immune cells in the peripheral blood of healthy controls. [6][7][8] We have here applied this systems immunology approach to the study of MS and, for comparison, autoimmune thyroid disease (AITD), characterized by familial clustering and shared genetic risk factors with MS. 9,10 In an unbiased analysis, we revealed the immunologic signatures of MS and 4 current immunomodulatory treatments. The prominent alterations in B cell subsets, most importantly the increase in B cell-activating factor (BAFF) and transitional B cells induced by both interferon-b and fingolimod, indicate the importance of the early B cell pathway in treating MS. METHODS Study population. Unrelated patients of Caucasian descent fulfilling McDonald criteria for MS or criteria for clinical AITD, and spouses as healthy controls, were included in the study over a 22-month inclusion period (table 1). Extensive demographic and clinical data were collected through a questionnaire and medical records. Exclusion criteria were cancer, immunosuppressive treatments, antibiotics, anti-allergy or anti-inflammatory treatment in the week prior to sampling, and, for healthy controls, presence of known clinical autoimmune disease.
Standard protocol approvals, registrations, and patient consents. Participants gave written informed consent, and the study was approved by the Ethics Committee of the University Hospitals Leuven.
Immunophenotyping and cytokine measurements. Peripheral blood mononuclear cells (PBMCs) and plasma were isolated from heparinized blood and stored at 280°C prior to analysis. A multiplex flow cytometry platform, previously validated in healthy controls, 8 was applied (table e-1 and figure e-1 at Neurology.org/nn). Some variables were measured in subsets of the samples only. Measurements are expressed as percentage of parental or grandparental cells as indicated. BAFF plasma levels were measured using a human BAFF Quantikine ELISA (R&D Systems, Minneapolis, MN). Averages were taken over 2-3 replicates after standardization using samples present on multiple plates. Sensitivity was 6.44 pg/mL and average standard deviation 15%.
Gene expression. RNA was extracted from total PBMCs and reverse transcribed using a high-capacity cDNA reverse transcription kit (Life Technologies, Carlsbad, CA). Droplet digital PCR (Bio-Rad, Hercules, CA) with gene expression assays (Life Technologies) for the BAFF gene TNFSF13B (Hs00198106_m1) and housekeeping gene POLR2A (Hs00172187_m1) was performed using 4 ng/mL cDNA. Relative quantity of TNFSF13B vs POLR2A was measured with QuantaSoft v1.4 (Bio-Rad). Correlation between separate experiments (n 5 68) was 0.75.
Statistical analysis. For ratios of variables, the natural logarithm was taken. Spearman correlation for pairs of variables and significance in correlation coefficient between groups were calculated and plotted with heatmap.2 and cocor packages. Linear regression models in R 2.14.1 included the immunologic variable as function of disease or treatment and covariates sex, age, and, for disease effects, disease duration at sampling. Multiple testing correction was applied for 38 variables across 4 treatments (p 5 0.00033). Nonmetric multidimensional scaling was used to visualize dissimilarities of samples between treatments. The ordination was created in R 3.1.0 using R package vegan with custom scripts. Distance matrix was calculated with Bray-Curtis dissimilarity index on original data. To test for the contribution of each variable to the resulting distance matrix, we applied Mantel test with 999 permutations comparing the original distance matrix to one with one variable removed at a time.

RESULTS
The adaptive immune profile of MS. In order to directly study the adaptive immune system across conditions, we applied a systems immunology approach in 154 patients with MS, 36 controls, and, for comparison, 55 patients with AITD (table 1). For each individual, the frequency of 38 leukocyte populations or ratios was assessed in the peripheral blood utilizing multiparameter flow cytometry, with a focus on T-and B cell subsets (table e-1 and figure e-1). Phenotype-phenotype correlations were assessed across each group (figure 1), identifying coregulated leukocyte populations that validate the systems immunology approach used here, as demonstrated previously. 8 Treatment of MS resulted in substantial immunologic deviation (described below), so initial analysis was performed on untreated patients with MS. Comparison of the immunologic status of patients with AITD currently untreated or treated with the thyroperoxidase inhibitor methimazole revealed no consistent differences (data not shown), allowing merged analysis. Of the robust immunologic changes observed to occur in autoimmune disease patients compared to healthy individuals, there were, surprisingly, no overlaps between AITD and MS (table e-2). MS was correlated with perturbation in the CD8 1 T-cell compartment compared to healthy controls, with a shift towards effector memory CD8 1 T cells (CD4 2 CD8 1 CD45RA 2 CCR7 2 ) as a percentage of all CD8 1 T cells (p 5 0.012), and increased the proportion of B cells, resulting in a decreased T/B cell ratio. The increase in B cells was mainly due to an increase in memory B cells (CD4 2 CD19 1 IgM 1 CD27 1 ) as a percentage of total lymphocytes (p 5 0.0042), and was especially prominent early on in disease (figure 2, A and B, and table e-2). Disease course (bout onset or primary progressive), time since last relapse, and multiple sclerosis severity scores did not confound these observations The human adaptive immune profile in multiple sclerosis (MS) In A and B, the right upper half, above the diagonal, indicates coregulation between pairs of cell types in healthy controls (n 5 36) (red: positive correlation coefficient, blue: negative correlation coefficient, light gray: no data available). Unbiased clustering of coefficients was performed to group coregulated cell types. Upon clustering, a main distinct naive cluster was observed, consisting of coregulated CD4 1 T cells, CD4 1 /CD8 1 ratio, and naive T cells, together with recent thymic emigrants (RTE) of both CD4 1 and CD8 1 lineages. This naive cluster was negatively correlated with activated and antigen-experienced lymphocyte subsets. Several microclusters were additionally observed, with coregulation of regulatory T cells (Tregs) with central memory CD4 4, A and B). As BAFF drives the differentiation of transitional B cells, 12,13 we additionally measured levels of BAFF in the same individuals. We observed upregulated BAFF plasma levels for the AITD disease process (p 5 0.0011). In MS, levels were unaltered in untreated patients but were upregulated during treatment with interferon-b (high-dose p 5 2.69E-09) and fingolimod (p 5 2.17E-10), in parallel with changes in the B cell compartment (figure 5, A-C). Treatment-induced changes followed a dose-dependent model (for interferon-b) (figure 5, A-C) and decreased by increased treatment duration (data not shown). Changes in BAFF plasma levels were correlated with the observed changes in B cell immunophenotype, i.e.,  DISCUSSION In this study, we have undertaken a comprehensive and unbiased analysis of the adaptive immune system measuring 38 immunologic variables in 245 individuals. Our approach allows us to directly compare, on the same large-scale platform, MS and AITD, 2 distinct autoimmune diseases with elements of shared genetic risk 9,10 as well as 4 common immunomodulatory treatments for MS (interferon-b, glatiramer acetate, natalizumab, fingolimod), of which the mechanism of action is not fully understood. 3,4 Patients with MS displayed systematic immune deviations, with an overall increase in effector memory CD8 1 T cells and memory B cells. Both cell types have recently been reported as part of the inflammatory response in MS 14,15 and as populating the CNS in MS, [16][17][18] providing a biological parallel to these systems-immunology associations. Despite the shared genetic risk factors, the immunologic deviations observed for MS were distinct from AITD. This result opposes a model (implied from genome-wide association studies) where common immunologic processes and sex was applied. Significant differences (after multiple testing correction) in the peripheral immune system of patients with MS following treatment compared to controls and untreated patients with MS included altered proportions of (A)   underlie different autoimmune diseases, with the unique clinical manifestations being based on the specificity of the escaping forbidden clones. By contrast, the results here support a more complex gene-immune relationship, where largely overlapping genetic polymorphisms nevertheless synergize to form immunologic deviations that are distinct at the systems level. In addition to the different efficacy and safety profile of the 4 compounds in routine MS treatment, each induced unique constellations of immune deviations, which offers perspectives to the challenge of personalized medicine. 5 Alterations in B cells were unexpectedly prominent upon treatment; indeed, the only shared treatment effects were notably B cell-related. Our observations on the crucial role of B cells in immunomodulatory treatment of MS hence add to emerging evidence emphasizing the importance of B cells, which have long stood in the shadow of T cells, in both MS pathogenesis and treatment. 19 Our immunophenotyping platform suggests a janus role for BAFF in autoimmune disease and treatment. The pathogenic face of BAFF is exposed in AITD, where BAFF is upregulated and may contribute to Th17-mediated disease. By contrast, the protective face of BAFF is revealed as the convergence point of effective immunomodulatory treatment in MS. Here we confirm previously reported changes in BAFF protein level upon interferon-b treatment 20 and establish it as a unique shared mechanism with fingolimod. This pathway is furthermore shared with promising emerging lymphocyte and B cell depletion therapies for MS. 21,22 We demonstrate that these protein level changes are not merely due to reduced consumption by BAFF receptors on B cells but involve active upregulation of transcription. Both in AITD and upon immunomodulatory treatment of MS, we now correlate BAFF plasma changes to an increase in transitional and decrease in switched B cells as percentage of B cells and propose the BAFF pathway as the mechanism of action driving these B cell alterations. Transitional B cells may be enriched for autoreactive B cells and have been implicated in the pathogenesis of other autoimmune diseases such as Sjögren syndrome and systemic lupus erythematosus. 23 At the same time, they may also have regulatory properties supporting a beneficial effect as part of the therapeutic mechanisms of MS immunomodulatory treatments. 24 BAFF antagonists have been proposed as a promising strategy based on animal models for MS. 25,26 However, 2 clinical trials in MS with atacicept, neutralizing all forms of BAFF and the related a proliferationinducing ligand (APRIL), had to be stopped prematurely because of increased inflammatory activity in patients with MS, and a 2-fold higher rate of conversion from clinically isolated syndrome to clinically definite MS in the treated compared to placebo group. 27,28 A third trial with a monoclonal antibody against BAFF (tabalumab) was completed but no results are available. 29 Of note, animal models leading to these clinical trials had a Th17-mediated component, 26 which in our study we observe in human patients with AITD but not MS. The immunomodulatory treatments for which we observe an increase in BAFF, interferon-b and fingolimod, have notably distinct effects between MS, where they are beneficial, and the neuroinflammatory disease neuromyelitis optica (NMO), where these treatments are not only ineffective but may be harmful and result in excessive inflammatory activity and increased relapse rate and disability. 30,31 This is in line with our observation that the unique shared signature of both treatments is an increase in BAFF coupled to already increased levels in BAFF seen in NMO but not MS. 32 We propose that the distinct and shared immunologic architecture of disease and of existing efficacious treatments may inform further design of novel treatment strategies. In particular, our data suggest that antagonism of BAFF may be counterproductive to MS treatment but may be further pursued in NMO, where current immunosuppressive or immunomodulatory treatments are only partially effective. 33 Our study was based on a routine clinical setting in a tertiary outpatient clinic. While we corrected for heterogeneity in age and sex, additional variables affecting treatment decisions, patient compliance, and treatment response are potential confounders. Preexisting variation among the immune profiles of patients with MS is, however, unlikely to be responsible for the immunologic responses to treatment. While patients with MS were assigned to treatment groups based on standard clinical criteria (driven by disease factors), assignment between low-dose interferon-b, high-dose interferon-b, and glatiramer acetate was determined by patient-preferred regimen, a research design that would not produce the dosedependent effects observed with interferon-b or the distinct immunologic changes between interferon-b and glatiramer acetate-treated groups. Likewise, of the 14 fingolimod patients, 10 were included as part of a randomized clinical trial determined by patient preference for an oral treatment, a design structure that avoids segregation based on preexisting clinical variation. Furthermore, our recapitulation of previously reported findings in studies of single treatments typically limited by sample size and number of cell types investigated [34][35][36][37][38] validates the methodology and thus the novel associations also observed in this systems immunology approach. While we were only able to assess changes in the peripheral blood, and not the target tissue, there is active exchange between the CNS and peripheral blood, 18,39 and changes in the peripheral blood have previously been shown to be able to capture genetic variation important for MS. 40 Moreover, immunomodulatory treatments are not administered directly into the CNS, and thus the likely location of activity is in the periphery, where our analysis took place.
In response to the challenge put forward for largescale evaluation of the adaptive immune system in patients with MS, 3 we demonstrate that a comprehensive common immunophenotyping platform is able to dissect the unique immunologic signatures of MS and immunomodulatory treatments. The different effects of immunomodulatory treatments suggest that the immunologic profile may in the future assist the neurologist in taking treatment decisions and optimizing the efficacy and safety for subgroups of patients, hence offering perspectives to the challenge of personalized medicine in MS. 5 We propose that further design of novel treatment strategies takes into account the shared immunologic architecture of existing efficacious treatments, such as the early B cell changes with increase in BAFF and transitional B cells.