Research ArticleImmunologyNeuroscience Free access | 10.1172/jci.insight.126599
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Greenfield, A. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Dandekar, R. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Ramesh, A. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Eggers, E. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Wu, H. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Laurent, S. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Harkin, W. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Pierson, N. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Weber, M. in: JCI | PubMed | Google Scholar
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Henry, R. in: JCI | PubMed | Google Scholar |
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Bischof, A. in: JCI | PubMed | Google Scholar |
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Cree, B. in: JCI | PubMed | Google Scholar |
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Hauser, S. in: JCI | PubMed | Google Scholar |
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by Wilson, M. in: JCI | PubMed | Google Scholar |
1UCSF Weill Institute for Neurosciences, Department of Neurology, UCSF, San Francisco, California, USA.
2Institute of Neuropathology, Department of Neurology, University Medical Center Göttingen, Germany.
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
Find articles by von Büdingen, H. in: JCI | PubMed | Google Scholar |
Published February 12, 2019 - More info
B cells are key contributors to chronic autoimmune pathology in multiple sclerosis (MS). Clonally related B cells exist in the cerebrospinal fluid (CSF), meninges, and CNS parenchyma of MS patients. We sought to investigate the presence of clonally related B cells over time by performing Ig heavy chain variable region repertoire sequencing on B cells from longitudinally collected blood and CSF samples of MS patients (n = 10). All patients were untreated at the time of the initial sampling; the majority (n = 7) were treated with immune-modulating therapies 1.2 (±0.3 SD) years later during the second sampling. We found clonal persistence of B cells in the CSF of 5 patients; these B cells were frequently Ig class-switched and CD27+. Specific blood B cell subsets appear to provide input into CNS repertoires over time. We demonstrate complex patterns of clonal B cell persistence in CSF and blood, even in patients on immune-modulating therapy. Our findings support the concept that peripheral B cell activation and CNS-compartmentalized immune mechanisms can in part be therapy resistant.
B lymphocytes (B cells) are key players in the immunopathology of multiple sclerosis (MS) (1–6). B cell–depleting anti-CD20 antibody therapies effectively reduce relapsing MS disease activity and slow the accumulation of disability in relapsing and primary progressive disease (2, 3, 7). Early on in the disease course, B cells enter the CNS, cerebrospinal fluid (CSF), and meningeal compartments (8–10), where they become compartmentalized and likely contribute to ongoing CNS tissue damage and consequential disability progression.
Ig gene transcripts, particularly heavy chain variable regions (Ig-VH), provide molecular fingerprints that permit temporal and spatial tracking of clonally related B cells (11–13). Single–time point immune repertoire studies in MS patients showed that there is an anatomic continuum of clonally related B cells extending from MS lesions in the brain parenchyma to CSF, meningeal inflammatory infiltrates, cervical lymph nodes, and peripheral blood (PB) (11, 12, 14–16). Individual MS patients have a characteristic CSF Ig electrophoretic pattern (oligoclonal bands [OCBs]) that persists over many years; the stability in the pattern of OCBs within a patient is considered evidence that clonally related B cells are long-lived in the CNS (17, 18). However, because multiple B cell clones can give rise to a single OCB (16), higher resolution genomic techniques are required to confirm this hypothesis.
Here, we were interested to determine the persistence of clonally related B cells in the CSF and in PB over time. Using previously described immune repertoire sequencing technology and bioinformatics tools (8), we analyzed Ig-VH sequences (IgG-VH and IgM-VH) (19) from 10 patients’ PB and CSF B cell subsets at 2 time points to track related B cells over time and across compartments. Work by others previously described 2 untreated MS patients with clonally related CSF B cells over time (20). Here, we were interested to understand which specific B cell subsets comprise clonally related CSF B cells that remain longitudinally detectable, their clonal relationship to PB B cells, and whether CSF B cells persist after treatment with MS immune-modulating therapy (IMT). We identified persistent, clonally related B cells in the CSF of 5 of the 10 patients, 4 of whom had initiated IMT between the first and second time points. We also identified patterns of clonal B cell input from the periphery to the CSF over time, suggesting that functionally diverse CD27+ PB memory B cells are a likely peripheral reservoir of B cells involved in MS disease activity.
Patient demographics. Ten patients ranging in age from 24 to 52 years old were enrolled (see Table 1 for clinical data and Supplemental Table 1 [supplemental material available online with this article; https://doi.org/10.1172/jci.insight.126599DS1] for standard laboratory CSF findings). On average, 12.2 ml (±3.5 SD) of CSF was collected per patient per time point. Eight patients had relapsing-remitting MS (RRMS), and 2 had primary progressive MS (PPMS). All were untreated at time point 1 (T1), and 7 of the 8 RRMS patients were on IMT by time point 2 (T2) (Table 1). Because there were no approved IMTs for PPMS before 2017, neither PPMS patient was on an IMT. IMTs used in this cohort varied from lower efficacy to higher efficacy treatments. Eight of 10 patients (6 treated, 2 untreated) had enhancing or new lesions on brain or spinal cord MRI at T2 (Table 1). All patients had OCBs unique to the CSF at both time points (Supplemental Table 2). In 7 of 10 patients, we were able to directly compare the OCB pattern at each time point (Supplemental Figure 1): 5 of 7 patients had stable band patterns, 1 had a decrease in band number, and 1 had an increase in band number (Supplemental Figure 1 and Supplemental Table 2).
Flow cytometry. B cell subsets were identified and sorted by flow cytometry (MoFlo Astrios EQ). Flow cytometric B cell subset distribution (see Methods) was determined at T1 and T2 for n = 8 patients’ CSF and for n = 9 patients’ PB; flow cytometry was not available for time point 1 CSF (T1-CSF) and peripheral blood (T1-PB) of 1 patient and for both CSF time points of another patient. When compared with PB, the CSF was enriched in CD19+CD27+IgD– Ig class-switched memory (SM) B cells (Supplemental Figure 2), consistent with previous reports (8, 21, 22).
Immune repertoire sequencing. IgG-VH and/or IgM-VH repertoire sequencing cDNA libraries were prepared from 167 samples. Samples consisted of PB or CSF FACS-sorted B cell subsets or, alternatively, bulk CSF or PB mononuclear cells (Supplemental Table 3). Sequencing libraries could not be obtained from 16 samples (Supplemental Table 3). From the remaining 151 samples, we generated 583,932 (±652,920 SD) raw reads per library. We identified 218,401 (±308,602 SD) Ig-VH sequences per library from the Ig heavy chain variable germline segment (IGHV), Ig heavy chain joining germline segment (IGHJ) and Ig heavy chain complementarity-determining region 3 (H-CDR3) for further analysis (Supplemental Table 3). Ig-VH sequences were clustered using a distance metric approach (see Methods). Samples with more B cells had more Ig-VH clusters (Supplemental Figure 3) (r = 0.88, P < 0.0001 for all samples; r = 0.74, P < 0.0001 for PB; r = 0.59, P < 0.0001 for CSF, Spearman’s correlation). Five paucicellular B cell subsets yielded more Ig-VH clusters than the number of input cells (Supplemental Table 3). For these samples, we analyzed the same number of Ig-VH clusters as input cells, choosing the Ig-VH clusters with the greatest number of aligned sequencing reads. Mutational analyses within Ig-VH clusters were not performed because these were not needed for the conclusions of this study.
At T1, we identified CSF Ig-VH clusters that were exclusively IgG-VH in all 10 patients (26.4 [±28.3 SD] Ig-VH clusters/patient); of the 10 patients, 9 patients also had CSF Ig-VH clusters that contained exclusively IgM-VH (44.8 [±57.3 SD] Ig-VH clusters/patient), and in 5 patients, we found mixed IgM and IgG clusters (5.2 [±9.4 SD] Ig-VH clusters/patient) (Supplemental Figure 4). At T2, we found that all 10 patients’ CSF contained Ig-VH clusters that were exclusively IgG-VH (42.6 [±72.6 SD] Ig-VH clusters/patient) or exclusively IgM-VH (31.7 [±33.9 SD] Ig-VH clusters/patient); in 7 patients, there were 6.7 (±11.8 SD) Ig-VH clusters/patient that were mixed IgM and IgG (Supplemental Figure 4). At T1, SM and naive B cells were common members of CSF Ig-VH repertoire clusters: These subsets were prevalent in repertoires of 3 out of 5 and 2 out of 5 patients with sorted T1-CSF B cells, respectively (59.6 [±80.4 SD] Ig-VH clusters/patient, 60.8 [±65.9 SD] Ig-VH clusters/patient, respectively) (Supplemental Figure 5). SM B cells commonly contributed to T2-CSF: 5 of 8 patients with sorted T2-CSF B cells had SM-predominant repertoires (45.4 [±65.9 SD] Ig-VH clusters/patient) (Supplemental Figure 5).
Clonally related B cells persist in MS CSF. In 5 of 10 patients, we identified “persistent CSF Ig-VH clusters” in which CSF Ig-VH sequences from both time points were represented (Figures 1–3); we thus demonstrate that B cells found in MS patients’ CSF at different time points are clonally related. Aside from Ig-VH sequences that exclusively persisted in CSF (Supplemental Figure 6), we identified 3 possible associations of CSF Ig-VH clusters with PB repertoires: a T1-PB connection, a T2-PB connection, or connections with both PB time point samples (Figure 2). We found IgG-expressing B cells, including SM B cells and plasmablast/plasma cells (PCs), in persistent CSF Ig-VH clusters of all 5 patients with persistent CSF Ig-VH clusters (Figure 3). In contrast, IgM-expressing B cell subsets were found to take part in persistent CSF Ig-VH clusters in only 2 patients (patients 1 and 3) (Figure 3). In particular, we did not find naive CSF B cells in persistent CSF Ig-VH clusters.
Persistent CSF Ig-VH clusters are present in MS patients. Shown are the numbers of Ig-VH clusters within CSF (blue circles) and PB (red circles) at T1 and T2. Circle overlap values, Ig-VH clusters found in both CSF and PB at T1 or T2; nonoverlap circle values, Ig-VH clusters exclusively found in CSF or PB, not both. Arrows, of total Ig-VH clusters in CSF or PB (nonoverlap portion of circle + circle overlap), the number of Ig-VH clusters that are present at both T1 and T2. Patients (Pts) 1–5 have persistent CSF Ig-VH clusters; Pts 6–10 did not have identifiable persistent CSF Ig-VH clusters.
Persistent CSF Ig-VH clusters have multiple patterns of persistence in relation to PB. Persistent CSF Ig-VH clusters can be exclusive to CSF (far left, Pts 1, 2, 5), link to T1-PB (second from left, Pts 1, 2), link to T2-PB (second from right, Pt 4), or link to both PB time points (far right, Pts 1–4). Blue circles, CSF; red circles, PB; circles on left, T1; circles on right, T2; arrows, shared Ig-VH clusters between time points or between CSF and PB.
B cells that persist in the CSF are often SM and plasma cells. (A) Clonally related B cell subsets are connected by lines between CSF (top half of circle) and PB (bottom half) at T1 (left half) and T2 (right half). Each segment along the circle represents a B cell subset from CSF or PB. Lines connecting B cell subsets indicate B cell subsets that share 1 or more Ig-VH clusters in common. Gray lines, PB-only Ig-VH clusters; blue lines, persistent CSF Ig-VH clusters; red lines, CSF Ig-VH clusters with PB contribution; yellow lines, T1-PB Ig-VH clusters that share clonal relatives with T2-CSF yet not with T1-CSF. Gray font indicates subsets or Ig isotypes from which no Ig-VH libraries could be obtained. Persistent CSF Ig-VH clusters were linked to PB Ig-VH sequences derived from Ig class-switched memory (SM) B cells (n = 4 patients), double-negative (DN) B cells (n = 3 patients at T1, n = 2 patients at T2), and plasmablasts/plasma cells (PCs) (n = 2 patients at T1 and 1 patient at T2). (B) T1-PB B cell subsets that are clonally related to T2-CSF without involving T1-CSF (i.e., quantitation of yellow lines depicted in A). USM, unswitched memory.
There was no relationship between the total CSF white blood cell count (clinical diagnostic that includes all leukocytes, not just B cells) at T1 or T2 and the number of persistent CSF Ig-VH clusters (r2 = 0.32, P = 0.09, and r2 = 0.07, P = 0.47, respectively; not shown). There was no significant correlation between the detection of persistent CSF Ig-VH clusters over time and the following clinical characteristics: age, disease duration, time between T1 and T2, EDSS score, use of IMT, gadolinium enhancement on MRI at either time point, IgG index, relapse between T1 and T2, or number of OCBs (Supplemental Tables 1 and 2). More of the patients without persistent CSF Ig-VH clusters were men (Supplemental Table 1). There was a numerical (yet not statistically significant) imbalance in the number of patients with relapses in that more of the patients with identified persistent CSF Ig-VH clusters experienced a relapse. There was no significant difference between B cell subset distribution in patients with or without persistent CSF Ig-VH clusters (Supplemental Figure 7). We found no correlation between the average reads per cell and our ability to detect persistent CSF Ig-VH clusters (r = 0.2, Spearman’s correlation; P = 0.17). The percentage of Ig-VH clusters in CSF that were persistent was similar to the percentage of Ig-VH clusters that were persistent in PB (i.e., Ig-VH clusters formed from PB at both time points) (Supplemental Table 4). This suggests that immune processes such as antigen stimulation and B cell regulation, which maintain a proportion of B cell clonal populations in the periphery, may also occur in the CNS.
PB SM and unswitched memory B cells are potential contributors to CSF B cell repertoires. Clonal connections between PB and CSF were identified in all patients, irrespective of whether persistent CSF Ig-VH clusters or clonal relatedness within the same PB and CSF subset was present (Figure 3 and Supplemental Figure 8). Therefore, to understand which PB B cell subsets might contribute clonal input to CSF repertoires over time, we examined Ig-VH clusters containing T1-PB and T2-CSF members yet no T1-CSF members. In all 9 patients with sorted T1-PB subsets, T1-PB unswitched memory (USM) (6 patients) and SM (5 patients) B cells were related to T2-CSF (Figure 3 and Supplemental Figure 8). All of these patients except for 1 (patient 5) were treated with IMT by T2. Other T1-PB B cell subsets sharing repertoires with T2-CSF included naive (4 patients), double negative (DN) (3 patients), and PC (3 patients) (Figure 3 and Supplemental Figure 8). Although overall somatic hypermutation (SHM) rates followed expected patterns for B cell subsets based on maturation stage, certain CSF subsets, such as IgG SM in patient 1, were noted to have a particularly high degree of SHM — suggesting affinity maturation (Supplemental Figure 9). Whether Ig class-switching occurs in USM B cells has been debated (23); here we found numerous Ig-VH clusters containing both USM (CD27+IgD+IgM+) and class-switched IgG+ B cells (i.e., SM, PC, and DN) in CSF and PB (2.5 [±2.6 SD] Ig-VH clusters/patient in CSF, 125.6 [±96.6 SD] Ig-VH clusters/patient in PB) (Supplemental Figure 10). To explore whether CSF PCs might be derived from long-lived PB PC populations, we searched for, but did not find, persistent PB PCs that were clonally related to CSF PCs. Incidentally, we found a frequent occurrence of persisting IgM-VH clusters with PB naive B cell members at T1 and T2, suggesting long-term survival and clonal expansion of naive B cells; previously, naive B cells have been characterized as short-lived precursors to memory and plasmablasts/PCs (24).
We identified clonally related B cells in CSF that were detectable at time points between 9 and 22 months apart. Given the histopathological and biochemical evidence for B cell recruitment, survival, and proliferation in the CNS compartment (25–27), the ability to detect clonally related CSF B cells by immune repertoire sequencing over time is likely the result of local persistence (i.e., long-term survival) of B cells in the functionally connected CNS, CSF, and meningeal compartments. This is perhaps consistent with the disease activity observed in some patients examined in this study and is also consistent with prior findings of antigen-experienced B cells and soluble Ig in the CSF that suggests persistence over time. These cells could constitute a tissue-resident immune repertoire capable of being restimulated — either directly or by bystander activation — to rekindle inflammation, mechanisms that may be shared with rheumatoid arthritis and other tissue-specific chronic autoimmune diseases (28). B cells of persistent CSF Ig-VH clusters display phenotypic features indicative of antigen-driven immunity (i.e., CD27-expression and Ig class-switched B cell receptors). Persistent B cells in PB, on the other hand, belong to a less restricted set of B cell subsets and include apparently long-lived naive (IgD+CD27–) B cell populations in addition to memory B cell and PC subsets. Overall, persistent CSF B cells comprise a minority of the overall CSF B cell diversity. Nonetheless, their persistence despite IMT initiation suggests that these B cell populations could be representative of a therapy-resistant, CNS-compartmentalized reservoir of immune cells. OCBs as well are evidence for CNS/CSF compartmentalization, albeit only of PCs as producers of soluble Ig. CNS-compartmentalized B cells are believed to contribute to progressive CNS tissue damage in MS (9, 25, 29); identification of CSF-persistent, functionally diverse B cell subsets (e.g., memory B cells, PCs) suggests that multiple aspects of active adaptive immunity are reflected in the CNS that are capable of supporting the immunopathogenesis of MS in various ways (i.e., cytokine secretion, antigen presentation, and antibody production). Interestingly, even after long-term high-efficacy IMT (30, 31), including intrathecal rituximab (32), patients continue to be at risk of progression. Indeed, our findings lend further support to the hypothesis that currently available IMTs may not sufficiently target progression-promoting mechanisms located in the CNS compartment. B cell–depleting therapies do not affect CD20– PCs, and oligoclonal antibodies persist in CSF of rituximab-treated patients (33, 34); future studies are needed to determine whether anti-CD20 therapies deplete persistent clonal B cells from the CSF, CNS, and meninges and whether doing so reduces disease progression. Our findings may suggest that development of novel therapies that target CNS-compartmentalized immune mechanisms could help close an important therapeutic gap.
Disease-driving lymphocytes, including B cells, are believed to migrate to the CSF and CNS from peripheral stores during phases of MS disease activity (8, 11, 12). Because of their immunological properties and functionality, SM B cells are generally believed to be highly relevant to the immune pathology of MS (35). Indeed, when persistent PB SM B cells were clonally related to second time point CSF B cells, these T2-CSF cells were nearly always SM and PC. This suggests that, in addition to CSF persistence, a population of potentially disease-relevant PB SM B cells may have a particular tendency to migrate to the CSF and mature to PCs. Also suggestive of peripheral maturation to SM B cells followed by migration and survival in the CNS, rather than local SM B cell to PC maturation in CSF/CNS, we did not find any T2-CSF B cells that were clonally related to T1-CSF naive B cells. In agreement with our previous work, we did not observe any relationships between PB-persistent PCs and CSF PCs, further supporting that migration to the CNS occurs at the SM or earlier B cell development stage followed by intrathecal maturation to PCs rather than immigration of long-lived peripheral PCs.
The efficacy of anti-CD20 B cell–depleting therapy in autoimmune disease is linked to depletion of memory B cell subsets rather than naive and transitional subsets (36, 37). We previously proposed that USM B cells may play a role in MS because of their appearance in MS patients’ CSF during periods of active CNS inflammation (8). Interestingly, USM B cells follow the same depletion and reconstitution pattern as SM B cells under rituximab therapy (38). Here, we observed that clonally related SM and USM B cells, which remain detectable in PB over time, may also have clonal relatives in CSF at the later time point. This finding is consistent with the concept that, upon peripheral activation, CXCR5+ USM or SM B cells home to the brain parenchyma or leptomeninges along a CXCL13 gradient, where they participate in the formation of ectopic tertiary germinal centers (8, 9, 39, 40) capable of supporting B cell receptor affinity maturation and class-switch recombination (41). Whether USM B cell receptors undergo Ig class-switching and contribute to high-affinity surface receptor or secreted antibody–mediated immune responses has been debated (42, 43). Here, we found numerous Ig-VH clusters linking CSF or PB IgM-expressing USM B cells to clonally related IgG+ B cells in CSF and PB (i.e., SM, PC, and DN). This finding suggests either (a) that USM B cells can enter germinal centers, where they undergo class-switch recombination or (b) that divergent maturation of naive B cells into functionally diverse IgM- and IgG-expressing memory B cell subsets occurs. Similar to SM B cells (44), USM B cells can likely function as potent antigen-presenting cells (45). Therefore, the depletion of both SM and USM B cells could be important for the efficacy of B cell–depleting MS therapies. However, whether USM B cells are involved in a true antigen-directed immune response or, conversely, whether they provide a nonspecific stimulus for intrathecal immune activation remains unclear. In either case, because of their participation in clonally related immune connectivity between PB and CSF, SM and USM B cells may be carriers of MS disease activity.
There are a number of limitations to this study. First, it is important to keep in mind that it is not technically possible to interrogate the entire immune repertoire in a living patient. As a result, persistent B cells could have gone undetected, and patients in whom we did not find clusters comprising clonally related Ig-VH may in fact harbor persistent B cells in the CSF or in other tissues. Although it is possible that deeper sequencing of existing samples could have increased our sensitivity to identify persistent CSF Ig-VH clusters, and although some samples did not generate productive sequences, we found no correlation between the average reads per cell and our ability to detect these clusters. Second, technical advances, such as incorporation of unique molecular identifiers that aid in error correction and quantification, were not a standard part of the sequencing method at the start of the study, and the initial method was maintained for consistency across this cohort. Ig-VH sequences obtained in this study overall reflected IGHV SHM accumulation as expected for the different B cell maturation states studied, with the lowest SHM in naive B cell IgM-VH (Supplemental Figure 9). The fact that some Ig-VH sequences displayed SHM counts above expected rates is likely owed to errors introduced by PCR during Ig-VH library preparation and Ion Torrent sequencing. However, PCR or sequencing errors are unlikely to influence determination of presence versus absence of related Ig-VH sequences but may influence the quantity of Ig-VH clusters. For this reason, we chose a distance metric (see Methods) that minimizes the risk of overestimating the number of Ig-VH clusters that could result from sequencing errors introduced by high-throughput technology. Third, although to our knowledge, this is the largest study of sorted B cell subsets from paired CSF and PB samples evaluated at 2 time points, only 10 patients were studied. Future studies will be needed to confirm these findings in larger cohorts, including patients treated with B cell–depleting therapies. Fourth, our experimental approach does not allow us to determine the directionality of B cell migration; to what extent lymphocyte migration pathways via “glymphatics” (46, 47) underlie our findings remains subject to future studies that will require reliable in vivo lymphocyte-tracking techniques. The question of whether clonally related B cells mature and persist in the CNS or mature and then migrate from the periphery is an important one. However, we refrained from performing mutational lineage analyses to determine B cell migration directionality because intermediate mutants that are necessary to obtain a complete picture of Ig-VH evolution may have gone undetected in our study.
In summary, we describe CSF-persistent clonally related B cell populations as part of an overall complex picture of B cell relatedness over time and across compartments that appears at least partially independent of IMT initiation. We identified previously unknown patterns of B cell relatedness, including PB SM and USM B cells that may contribute to SM- and PC-rich CSF B cell repertoires and naive B cells that persist over time. Our findings raise the possibility that IMTs, which effectively reduce clinical and MRI relapses, may not sufficiently target CNS-compartmentalized immune mechanisms that have been associated with gradual and relapse-independent MS worsening in relapsing MS and progressive forms of the disease (48). Effective targeting of CNS-compartmentalized inflammation may, in the future, achieve efficacy in progression slowing superior to that of current therapies (2). Longitudinal immune repertoire studies performed alongside the clinical development of such potentially novel MS drugs may well provide critical immunological and biomarker evidence indicating the successful suppression or elimination of progression-driving mechanisms.
Experimental design. We hypothesized before initiating this study that in MS patients there would be a population of B cell clones that persist in the CSF over time. Inclusion criteria for this observational study were (a) enrollment in the UCSF Expression, Proteomics, Imaging, Clinical (EPIC) study (31), which is a longitudinal MS cohort that defines MS according to the 2001 McDonald criteria (49), and (b) CSF and PBMCs from biobanks from 2 time points (T1 and T2) taken a minimum of 9 months apart. All patients were treatment naive at T1. Patients either consented to donating excess CSF and blood during routine lumbar puncture (LP) at UCSF, or if an LP was not ordered by the treating physician, consented to undergo LP and blood draw for research purposes. Patient demographics and clinical history, including physical exam and EDSS scores, were also obtained as part of the EPIC study.
Sample collection and processing. We centrifuged 7-30 ml of fresh CSF at 400 g × 15 minutes at 4°C to separate a cell pellet from supernatant. PBMCs were isolated from whole blood via Ficoll gradient followed by red blood cell lysis and washing with PBS with 1% BSA. We stored 7 of 20 CSF samples and 1 of 20 PBMC samples at –80°C as unsorted cell pellets (Supplemental Table 3). The remaining CSF and PBMC samples were immediately blocked with FcR Blocking Reagent (Miltenyi Biotec) and stained with fluorescent antibodies to cell surface markers — CD19, CD27, CD38, CD138, and IgD (Supplemental Table 5); some samples were stained for additional immune cell markers for experiments outside the scope of this study (antibody panels listed in Supplemental Table 5 and Supplemental Table 6).
CSF and PB B cell subsets were sorted as previously described (8) on a Beckman MoFlo Astrios EQ FACS device: Naive B cells (CD19+IgD+CD27–), USM B cells (CD19+IgD+CD27+), SM B cells (CD19+IgD–CD27+), DN B cells (CD19+IgD–CD27–), CSF CD27hi plasmablasts/PCs (CD19+IgD–CD27hi) (the vast majority of which are also CD38+), and PB plasmablasts/PCs (CD19+IgD–CD27+CD38+). Patient 4 had combined CD27+/hi CSF-sorted B cells. CSF samples for which low B cell count was a concern were bulk sorted — flow cytometry data were attained, but all CSF lymphocytes were placed into a single sample tube rather than segregating into the above 5 subsets. B cells were sorted directly into lysis buffer suitable for later RNA extraction (Qiagen, Buffer RLT with 1% β-mercaptoethanol) and stored at –80°C.
Immune repertoire sequencing. RNA was extracted using the RNeasy Mini Kit for greater than 200,000 cells or the RNeasy Micro Kit for fewer than 200,000 cells (Qiagen). RNA was reverse-transcribed using iScript reverse transcriptase (Bio-Rad). IgG-VH and IgM-VH sequences spanning the variable region from framework 1 to the constant region 1 were amplified using the Advantage 2 PCR Kit (Clontech) and custom primers as described (8). We performed 30 PCR cycles on PB samples and 40 cycles on CSF samples; up to 10 additional cycles were performed, in 5-cycle increments, if no PCR product was visualized on a 1.5% agarose gel after the initial round of amplification. Amplicons of the expected length were gel purified and quantified by Bioanalyzer using the High Sensitivity DNA Kit (Agilent). Equimolar amounts of PCR products were pooled to generate a 16 picomolar library for Ion Torrent sequencing (Thermo Fisher Scientific) following the manufacturer’s instructions, including emulsion PCR (Ion OneTouch2), enrichment for template-positive ion sphere particles, and quality control using a Qubit Fluorometer (Ion Sphere Quality Control Kit). Next-generation sequencing was performed on an IonTorrent Personal Genome Machine. If inadequate sequencing data were obtained during an initial sample’s sequencing, a technical replicate was performed.
Bioinformatics. Raw sequence FASTQ files were generated with Torrent Suite software (Thermo Fisher Scientific). A custom bioinformatics pipeline incorporating MiXCR (v2.1.3) (50) was used to identify IGHV and IGHJ germline segments and H-CDR3 for each sequence read. MiXCR first assembles reads with bases of quality scores greater than 20 and then attempts to align reads with bases of lower (<20) quality to the already assembled reads. We applied a bioinformatics clustering approach as described previously (8) to compile clonally related Ig-VH sequences into “Ig-VH clusters.” An Ig-VH cluster contains clonally related sequences — IgG or IgM — from at least 1 B cell. B cell receptors using identical IGHV and IGHJ germline segments and identical or near-identical H-CDR3 regions (minimum 8 amino acid length, maximum Hamming distance of 2, ref. 51) at the amino acid level were considered members of a clonally related Ig-VH cluster. This distance metric permitted clustering of identical and clonally related Ig-VH irrespective of sequencing errors that might have introduced a frameshift in the Ig-VH. Immune repertoire data sets can have falsely inflated diversity as a result of errors during the reverse transcription, amplification, and sequencing steps as well as incorrect barcode assignment that occurs when sequencing libraries are pooled. Thus, to conservatively catalog immune repertoire diversity in a given sample, we excluded any clusters that had fewer than 2 aligned read counts and any additional low-abundance Ig-VH clusters whose inclusion resulted in the total number of clusters exceeding the number of input cells. We excluded samples for which library generation was unsuccessful from statistical analyses.
Network diagrams were generated to visualize Ig-VH clusters using Circos (52). Amino acid alignment figures were generated using Geneious version 10.0 (created by Biomatters, https://www.geneious.com), IMGT V-QUEST (53), and NCBI IgBlast (54). Venn diagrams and other figures were based on custom R scripts that tabulated the number and type of Ig-VH clusters shared across compartments (CSF/PB) and time points.
MRI acquisition and analysis. MRI scans at T1 and T2 were acquired on the same 3T Siemens scanner following a standardized protocol that included 3D T1-weighted magnetization-prepared gradient echo images (MPRAGEs), 3D fluid-attenuated inversion recovery, and T2/proton density images of the brain and whole-spine T2-weighted and short tau inversion recovery images. In addition, brain MPRAGEs and whole-spine T1-weighted images were acquired after administration of gadobutrol. Readings were performed by a neurologist with subspecialty interest in neuroimaging (AB, 11 years of experience) who was blinded to the clinical and immunological data.
Data and materials availability. Ig raw sequencing data can be found at the National Center for Biotechnology Information Sequence Read Archive (BioProject accession number PRJNA517691).
Statistics. GraphPad Prism software (version 8) was used to perform statistical analyses: 2-tailed, unpaired t tests were used to compare pairwise groups, and Fisher’s exact test was used to compare categorical variables. One-way ANOVA was used for comparisons between groups with approximate Gaussian distributions, with correction for multiple comparisons performed using the Sidak method. Kruskal-Wallis test with Dunn’s correction for multiple comparisons was used for multiple comparisons with nonparametric distributions. Spearman’s correlation coefficient was used to determine correlation. P value less than 0.05 was considered statistically significant.
Study approval. All studies were approved by the UCSF Institutional Review Board, and written informed consent was obtained from each participant before inclusion in the study.
ALG performed immune repertoire sequencing and data analysis and wrote the manuscript. HCVB conceptualized the study, assisted in patient recruitment, provided guidance on analysis, and wrote the manuscript. MRW wrote the manuscript and provided guidance on analysis. RD and HW generated code, tables, and figures for analysis. ELE performed FACS and immune repertoire sequencing. MSW performed the OCB analysis. AB performed the analyses of MRI data. RGH oversaw acquisition of MRI data. SL performed immune repertoire sequencing and input on analysis. WH and NSP assisted in patient recruitment, consenting, and sample collection. BACC led the EPIC cohort study. SLH conceived of and led the EPIC cohort study. AR provided consultation on immune repertoire analysis and figures and reviewed the manuscript.
The authors would like to express their gratitude to the individuals who agreed to participate as patients in this study, as well as the UCSF MS EPIC study team of researchers and clinicians. Our studies were supported by grants from the National MS Society (RG-4868 to HCVB), the NIH/National Institute of Neurological Disorders and Stroke (K02NS072288 and R01NS092835 initially to HCVB, transferred to SLH; K08NS096117 to MRW), and the Valhalla Foundation, as well as by gifts from the Friends of the Multiple Sclerosis Research Group at UCSF. HCVB and MRW were also supported by an endowment from the Rachleff Family Foundation. ALG is supported by the National MS Society Kathleen C. Moore Postdoctoral Fellowship Clinician-Scientist Development Award (FAN-1507-05497).
Address correspondence to: Michael R. Wilson, 675 Nelson Rising Lane — Box 3206, San Francisco, California 94158, USA. Phone: 415.502.7429; Email: michael.wilson@ucsf.edu.
ELE’s present address is: Adaptive Biotechnologies Corp., South San Francisco, California, USA.
HW’s present address is: Cambrian InnovaTech, Shenzhen City, Guangdong, China.
SL’s present address is: Department of Neurology, University Hospital Cologne, Cologne, Germany.
WH’s present address is: UCSD, School of Medicine, La Jolla, California, USA.
NSP’s present address is: Keck School of Medicine USC, Los Angeles, California, USA.
HCVB’s present address is: F. Hoffmann-La Roche, Basel, Switzerland.
Conflict of interest: ALG has a patent pending for Method for High Percentage Recovery of Rare Cells (European application 18185007.4). MSW receives research support from the DFG (WE 3547/5-1), Novartis, Teva Pharmaceutical Industries, Biogen Idec, Roche, Merck, and the ProFutura Program of the Universitätsmedizin Göttingen. MSW is serving as an editor for PLoS One. He received travel funding and/or speaker honoraria from Biogen Idec, Merck Serono, Novartis, Roche/Genentech, Teva Pharmaceutical Industries, Bayer, and Sanofi Genzyme. RGH has received research funding from Roche/Genentech and Sanofi Genzyme. BAC has received personal compensation for consulting from EMD Serono and Novartis. SLH currently serves on the scientific advisory board of Symbiotix, Annexon, Bionure, Molecular Stethoscope, and Alector and on the board of trustees of Neurona. SLH also has received travel reimbursement and writing assistance from F. Hoffmann-La Roche Ltd. for CD20-related meetings and presentations. MRW has received research support from Roche/Genentech and has a patent pending for Method for High Percentage Recovery of Rare Cells (European application 18185007.4). HCVB is an employee of F. Hoffmann-La Roche Ltd., Basel, Switzerland (Roche) and reports equity in Roche and income as full-time employee.
Copyright: © 2019 American Society for Clinical Investigation
Reference information: JCI Insight. 2019;4(6):e126599. https://doi.org/10.1172/jci.insight.126599.