- Research article
- Open Access
Persistence and selection of an expanded B-cell clone in the setting of rituximab therapy for Sjögren’s syndrome
Arthritis Research & Therapy volume 16, Article number: R51 (2014)
Subjects with primary Sjögren’s syndrome (SjS) have an increased risk of developing B-cell lymphoma and may harbor monoclonal B-cell expansions in the peripheral blood. Expanded B-cell clones could be pathogenic, and their persistence could exacerbate disease or predispose toward the development of lymphoma. Therapy with anti-CD20 (rituximab) has the potential to eliminate expanded B-cell clones and thereby potentially ameliorate disease. This study was undertaken to identify and track expanded B-cell clones in the blood of subjects with primary SjS who were treated with rituximab.
To determine whether circulating B-cell clones in subjects with primary SjS emerge or remain after B cell-depleting therapy with rituximab, we studied the antibody heavy-chain repertoire. We performed single-memory B-cell and plasmablast sorting and antibody heavy-chain sequencing in six rituximab-treated SjS subjects over the course of a 1-year follow-up period.
Expanded B-cell clones were identified in four out of the six rituximab-treated SjS subjects, based upon the independent amplification of sequences with identical or highly similar VH, DH, and JH gene segments. We identified one SjS subject with a large expanded B-cell clone that was present prior to therapy and persisted after therapy. Somatic mutations in the clone were numerous but did not increase in frequency over the course of the 1-year follow-up, suggesting that the clone had been present for a long period of time. Intriguingly, a majority of the somatic mutations in the clone were silent, suggesting that the clone was under chronic negative selection.
For some subjects with primary SjS, these data show that (a) expanded B-cell clones are readily identified in the peripheral blood, (b) some clones are not eliminated by rituximab, and (c) persistent clones may be under chronic negative selection or may not be antigen-driven. The analysis of sequence variation among members of an expanded clone may provide a novel means of measuring the chronicity and selection of expanded B-cell populations in humans.
It has been estimated that up to 3 million adults in the US suffer from primary Sjögren’s syndrome (SjS) . Primary SjS is an autoimmune disorder characterized by chronic inflammation of the salivary and lacrimal glands and the presence of antinuclear antibodies, most often of the anti-SSA(Ro) and anti-SSB (La) specificities. Patients are often middle-aged females who present with sicca symptoms, such as dry eyes and dry mouth, fatigue, and joint pain, as well as other extraglandular manifestations, including lung disease and neuropathy. In primary SjS, it is believed that both T cells and B cells contribute to disease pathogenesis. Both cell types infiltrate the salivary and other exocrine glands and show evidence of clonal expansion in the affected tissues as well as the circulation . Notably, there is an increased risk of lymphoma in patients with primary SjS .
Why patients with primary SjS are at increased risk for lymphoma is unclear and has been the subject of several studies (reviewed in ). One theory is that B-cell hyperactivity in primary SjS results in the abnormal activation of autoreactive B-cells and contributes to their clonal expansion . Autoreactive B-cell clones, such as the recently described CD21dim population in SjS, may remain chronically activated instead of anergic in the presence of self-antigens . Autoreactive B-cell clones may also have increased resistance to apoptosis in primary SjS by virtue of elevated levels of the B-cell survival factor, BAFF (B-cell activating factor) . Another theory is that B cells in primary SjS accumulate and persist due to abnormal or inadequate regulation by other cells of the immune system. Borrowing an example from the field of tumor immunology, B cells that are transformed by Epstein-Barr virus are efficiently killed by T cells and natural killer cells . T cells in SjS not only may be derelict in their duties to constrain or kill transformed B cells but may have joined forces with the enemy: SjS T cells are the predominant inflammatory cell population in the exocrine gland lesions, appear to respond to autoantigens on apoptotic cells, secrete pro-inflammatory cytokines, and stimulate B cells (reviewed in ).
Clonal expansions that are evident from abnormalities in the blood may reveal underlying processes that can evolve into a malignant B-cell neoplasia. For example, some patients with monoclonal gammopathy of uncertain significance (MGUS), a condition characterized by the presence of a monoclonal immunoglobulin protein present in the serum, can progress to multiple myeloma . Similarly, a subset of patients with monoclonal B-cell lymphocytosis (MBL), a condition characterized by the presence of a clonal B-cell expansion in the peripheral blood, can progress to chronic lymphocytic leukemia .
We wondered whether a similar spectrum of conditions occurs in primary SjS, with some patients with primary SjS having polyclonal B-cell expansions, others having monoclonal expansions, and still others having transformed monoclonal expansions . Since primary SjS is a systemic autoimmune disease, expanded B-cell clones in the circulation could both be pathogenic (autoreactive) and be linked to an increased risk of lymphoma. However, in primary SjS-associated lymphoma, the B-cell neoplasm is often of the mucosal associated lymphoid tissue (MALT) type . It is not clear whether B-cell populations in the peripheral blood overlap with those in the glands and other tissues, although one intriguing report documents the presence of shared expanded clones in the lymph node tissue and peripheral blood of a patient with SjS who also had marginal zone lymphoma .
The results from several studies suggest that therapy with the B cell-depleting antibody, rituximab (anti-CD20), leads to symptomatic improvement in patients with primary SjS [13–15]. Furthermore, rituximab is known to modulate the antibody repertoire, particularly in antigen-experienced B cells [16, 17]. We therefore reasoned that rituximab therapy might also alter the circulating repertoire of expanded B-cell clones in patients with primary SjS who did not have lymphoma or overt evidence of MGUS or MBL. If CD20 antibody depletion of B cells eliminated expanded B-cell clones and re-set the B-cell repertoire, this approach might favorably modify the disease course in primary SjS by eliminating pathogenic B-cell clones that could be contributing to autoimmunity or predisposing to lymphoma. Conversely, the finding of persistent B-cell clones despite CD20 B-cell depletion therapy with rituximab implies that the restoration of the B-cell repertoire to a pre-disease state is incomplete at best and that the underlying mechanisms responsible for driving the disease process remain unchecked by B-cell depletion. Furthermore, the analysis of diversification by somatic hypermutation of the antibody gene rearrangements within the expanded clone might reveal mechanistic insights into how the clone is being selected (either by self antigen or by negative regulation from the immune system). Therefore, we analyzed the B-cell repertoire of SjS subjects at baseline and after treatment with rituximab by single B-cell sorting and immunoglobulin heavy-chain gene sequencing to identify and molecularly characterize expanded B-cell clones in the blood.
Subjects with primary SjS, as defined by the revised American-European Consensus Group criteria, were enrolled in an open-label Autoimmunity Centers of Excellence-sponsored clinical trial of the B cell-depleting antibody, rituximab (anti-CD20). Clinical features of the study subjects have been described previously . Subjects received two 1-gram doses of rituximab and 100 mg of methylprednisolone on days 0 and 15. For this study, the peripheral blood of six out of the 12 subjects with SjS were studied by flow cytometry, and antibody heavy-chain cloning was performed from single-sorted cells at baseline and at weeks 8, 14, 26, 36, and 52 after the second dose of rituximab . Subjects provided informed consent to participate in this study, which was carried out in accordance with a study protocol that was approved by the institutional review boards of both the Perelman School of Medicine at the University of Pennsylvania and Duke University Medical Center.
Single-cell cloning of immunoglobulin heavy-chain gene rearrangements
Memory and germinal center B-cell subsets (shown in Additional file 1: Figure S1) were recovered from the peripheral blood as described previously [18, 19]. Memory and plasmablast (PB) B-cell subsets were sorted into single wells of 96-well polymerase chain reaction (PCR) plates, and rearranged immunoglobulin heavy-chain (IgH) genes were cloned and sequenced as described previously .
Cloning of the germline VH1-69 immunoglobulin heavy-chain gene rearrangement
Genomic DNA from subject 2 (SjS2) was subjected to PCR amplification with primers that flanked the germline VH1-69 sequence. The primers used were VH1-69 germline forward: 5′-GTG CCC TGA GAG CAT CAC ATA ACA-3′ and VH1-69 germline reverse: 5′-TTC TCC CTC AGG GTT TCT GAC ACT-3′. The cycling conditions were 10 minutes at 94°C, followed by 40 cycles of 94°C for 30 seconds, 58°C for 30 seconds, and 72°C for 30 seconds, followed by 20 minutes at 72°C. Amplicons were TA-cloned (TOPO TA kit; Life Technologies, Grand Island, NY, USA) and sequenced.
Identification of VH, DH, and JH immunoglobulin heavy-chain gene sequences
For each sequenced immunoglobulin gene rearrangement, the most closely matching germline VH, DH, and JH alleles were identified, along with the CDR3 length, functional status (in frame, out of frame, termination codon), and predicted (translated) amino acid sequence by local alignment comparison with the ImMunoGeneTics (IMGT) database . Alignments were verified by using the high-throughput sequencing analysis algorithms from IMGT (V-QUEST version 1.1.2 ), and individual sequences were spot-checked with IgBLAST (National Library of Medicine). To avoid confusing mutations with sequencing imperfections, we considered only positions 10 to 381 of the sequences (by IMGT numbering). The numbers of identifiable IgH variable, diversity and joining (VDJ) rearrangements for each subject at each time point are summarized in Additional file 2: Table S1.
Clonal immunoglobulin gene tree construction
FASTA sequences for all of the clonal variant sequences were analyzed and graphically displayed by using neighbor joining with ClustalX, version 2.1  and by using the default parameters.
Detecting selection from immunoglobulin gene sequence mutation patterns
We used the focused selection test, which is explained in detail in , to detect the forces of selection in the complementarity-determining regions (CDRs) and in the framework regions (FWRs) of the B-cell clones. This test has been shown to be the only one to give reliable and accurate results inferring selection from replacement-to-silent (R/S) ratios  and takes into account micro-sequence specificity [25, 26] and transition bias  in somatic hypermutation. Taking these two factors into account, we calculated for any given set of nucleotides their expected mutability (that is, their relative likelihood to mutate). To estimate the expected level of selection, we count the number of silent mutations in the entire sequence and then use a micro-sequence-based targeting model of somatic mutation as described in  and the actual sequence composition of the mutant’s source germline sequence to estimate the expected number of R mutations [29, 30]. We then used the mutability of the sequence as a background to calculate, under a condition of no selection, the expected rate of R mutations in the area of interest (CDR or FWR) and the S mutations throughout the sequence. In this way, we avoided mixing potentially conflicting forces of selection in FWR and CDR. An overabundance of R mutations was considered an indication of positive selection, and a lack of R mutations was considered an indication of negative selection. When a set of clonally related sequences was analyzed, all unique mutations were grouped together, since mutations that occurred multiple times were probably an indication of common ancestry rather than being mutations that occurred more than once. In this manner, we raised the sensitivity above that obtained by analyzing a single sequence without spuriously counting mutations twice.
We further quantified selection strength by using BASELINe . This algorithm extrapolates the distribution of potential models of selection pressures that could result in the observed patterns of mutation. To accomplish this task, BASELINe considers all mutations from the same condition and considers their origin (clone/sequence/CDR/FWR). In this way, we compared clones with unequal sequence numbers. BASELINe outputs a selection strength (positive or negative) whose certainty is expressed by the width of the distribution curve. All calculations presented here were done with BASELINe’s web version .
We have previously described the results of our open-label pilot trial of rituximab therapy for primary SjS . In this companion study, our goal was to determine whether therapy with rituximab altered or re-set the repertoire in patients with primary SjS, using expanded B-cell clones as a read-out.
B-cell clones are frequent in peripheral blood memory B cells and plasmablasts of rituximab-treated Sjögren’s syndrome subjects
To explore the effects of CD20 B-cell depletion on the B-cell repertoire, antibody heavy-chain gene rearrangements were cloned and sequenced from sorted PBs and memory B cells (see Methods and Additional file 1: Figure S1) of six subjects with SjS before and at various time points after therapy with rituximab. Since B-cell clonal expansions have previously been described in patients with primary SjS [32–34], we reasoned that it would be possible to find and track an expanded B-cell clone in one or more of the subjects over the different study time points.
Here we define clonally related heavy-chain sequences as those that share the same VH (variable), DH (diversity), and JH (joining) gene segments; have the same CDR3 (third complementarity-determining sequence) length; and have highly related CDR3 sequences (identical or differ by only one amino acid or up to three nucleotides). Among the estimated tens of billions of different CDR3 sequences , rearrangements that share the same VH, DH, and JH; have the same CDR3 length; and have highly similar CDR3 sequence overall are very likely to be derived from clonally related B cells. Because the sequences were recovered from single cells, finding the same sequence in two separate cells (which were separately amplified), even if they were obtained at the same time point, was considered to be indicative of clonal expansion.
A remarkable feature of the antibody heavy-chain sequencing data is that clones were identified despite sequencing relatively few B cells: out of 303 sequences, 12 independently expanded clones were identified. Members of these clones comprised 80 out of the 303 sequences (26%). The distribution of clonally related sequences and the number of sequences sampled at each time point are shown in Figure 1 and in Additional file 2: Table S1. Owing to the low numbers of circulating B cells, very few sequences were obtained between weeks 8 and 26, but a high fraction of the sequences that were obtained at these time points (mostly from SjS2) were clonally related.
Analysis of the heavy-chain sequences revealed that the B-cell clones from different subjects (and in the same subjects) did not share any obvious sequence similarities (Table 1). Additional file 2: Table S1 summarizes the clonal expansions that were identified in the subjects with SjS. In SjS1, we identified two expanded clones among only 24 total sequences. In SjS3, we identified two clones in 30 sequences, and in SjS4 we identified two clones in 65 sequences. But the most interesting subject was SjS2, who had a large expanded clone that was present during all of the time points. Members of this clone were present in the circulation in spite of rituximab therapy and consisted of over 50 independently amplified sequences (FASTA files are provided in Additional file 3: Table S2). Also in SjS2, four clones (including the large clone) were present at baseline. These data indicate that clones are also detectable in the memory B-cell repertoire in at least some subjects prior to the administration of rituximab.
The large clone from SjS2 was identified as an in-frame rearrangement of VH1-69, DH7-27, DH4-23, and JH4-02. Additional file 4: Figure S2 shows the alignment of the most frequently recovered member of the clone with the corresponding germline VH, DH, and JH sequences. We could not account for 10 nucleotides in the CDR3 by their presence in the germline VH, DH, and JH genes. Since these 10 nucleotides were shared in nearly all of the sequences, we inferred that these nucleotides most likely arose by junctional diversification (n- or p-addition) rather than by somatic hypermutation.
Members of the large expanded clone are found in circulating plasmablast and memory B-cell pools
The distribution of the members of this large expanded clone with respect to time point and B-cell subset is shown graphically, using neighbor joining in Figure 2. The lengths of the branches are fairly long (distant) when compared with the germline (GL) sequence, consistent with mutations in most of the clone members. The shape of the clone alignment is more like a “bush” than a “tree”, showing no clear single direction of mutation. Rather, there are several small branches with sequence variants. The sequence variants do not display a clear-cut progression with time, at least during the time points that clonal variants were surveyed. This suggests that the major changes in the clone sequence likely occurred prior to this analysis of the clone. Of note, members of the clone were present in both CD38++ (PB) and CD38+/− class-switched (IgD−) memory (CD27+) B cells, indicating potential clonal overlap between these two pools of CD27+ peripheral B-cell pools. Small clusters of sequence variants that coincide with the memory B-cell subset can be seen, further suggesting that selection for the clone members could be distinctive within the PB and memory B-cell subsets.
Somatic mutations are frequent in the large expanded clone
To gain further insight into how the expanded clone was selected, we analyzed its pattern of mutation. Although it is possible that some of the intra-clonal variation is due to PCR error, we do not think that PCR error contributes substantially to the sequence diversity of this clone based on the large number of shared mutations between sequences and the ability to independently amplify identical sequences from different B cells.
We also considered the possibility that some of the mutations represented variations in the GL sequence due to the hypothetical presence of an unusual VH1-69 allele in subject SjS2. We therefore cloned the subject’s GL VH1-69 gene by using primers that were present in intronic sequences surrounding the unrearranged VH1-69 gene (see Methods). We obtained only a single VH1-69 sequence out of several independent PCRs (data not shown). Upon alignment (see Methods), the VH sequence matched VH1-69*01 (100% identity, Additional file 4: Figure S2). We therefore conclude that SjS2 is most likely homozygous for the VH1-69*01 allele and that the majority of the sequence differences among different members of the expanded clone are due to somatic hypermutation. Additional file 5: Table S3 summarizes the unique somatic mutations in the members of the expanded clone, relative to the GL sequences for VH1-69*01, JH4-02, and the CDR3 (comprised of DH7-27, DH4-23, and consensus in the areas with presumed n and p nucleotide insertions at the junctions).
Many of the somatic hypermutations in the expanded clone are silent
When the nucleotide sequences of the clone members are compared with the corresponding GL sequence, approximately half of the mutations, on average, are silent (Figure 3a). If the clone were positively selected, one would instead expect to find an increased frequency of R mutations compared with S mutations. A high frequency of S mutations is also not what one would expect to find at random. Owing to redundancies in the genetic code, the R/S ratio without any selection is generally approximately 3 to 1 in favor of R mutations. In contrast, the trends in this clone are all in the opposite direction: we find an overabundance of S mutations. However, the analysis in Figure 3a ignores the fact that many of the mutations among different members of the clone are shared and therefore are unlikely to be independent events. Therefore, we performed a more rigorous computational analysis of the somatic hypermutation pattern that took the shared mutations as well as their locations in the V region into account.
Not surprisingly, the computational analysis of the somatic hypermutation pattern was also consistent with negative selection in both the CDRs and the FWRs. Using the focused test of selection, replacement mutations are significantly under-represented among the mutations that were identified (P <0.05 at baseline and P <0.001 at other time points and overall) in all regions and overall (Table 2). To identify more clearly the strength of the negative selection and not just its presence, we used BASELINe on members of the clone (Methods), and we see quite clearly that there are high levels of negative selection (Figure 3b). Thus, there is negative selection in both the CDRs and the FWRs of the members of this clone. Furthermore, the other expanded clones identified in SjS2 and some of the other subjects exhibit somatic mutations that are consistent with negative selection (Additional file 6: Table S4). Thus, the mutation pattern found in the big clone in SjS2 is not unique to this particular clone and may be a more general feature of expanded clones in SjS.
Two additional features of this large B-cell clonal expansion become apparent when examining the distribution of mutations of the clonal variants at the six different time points (Figures 2 and 4 and Table 2). First, there is very little detectable accumulation of mutations over time, as most mutations are already observed at week 8. Second, the common mutations (that is, those found in many sequences across all time points; Figure 4) are not enriched for R mutations when compared with the rare mutations (that is, those found in only a few sequences and in one time point). We therefore considered it possible that none of these mutations is important for the formation of this lineage (that is, they may not confer selective advantage compared with the non-mutants). The common mutations could simply be derived from the same lineage that was present in the subject prior to the study.
Clonal expansion is a fundamental property of B cells that participate in an adaptive immune response. In patients with SjS, the identification, characterization, and tracking of expanded B-cell clones over time can provide insights into how clones are selected and whether they are affected by therapy. In this study, we show that expanded B-cell clones are present in the peripheral blood of patients with primary SjS undergoing B-cell depletion therapy with rituximab. There were three main findings: (1) a large expanded clone persisted in the blood of a patient with primary SjS (subject SjS2) despite B-cell depletion with rituximab; this finding indicates that rituximab does not fully re-set the B-cell repertoire in at least some patients with primary SjS; (2) the relative proportion of expanded B-cell clones among sequences recovered from memory and PB cells was highest during the period of maximal B-cell lymphopenia; this finding suggests that the sensitivity for detecting clonal expansions is greatest during periods of iatrogenic or disease-induced lymphopenia; (3) the expanded clones harbored numerous silent mutations; this finding indicates that expanded clones are under negative selection. The frequency and kinds of somatic mutations have potential implications for monitoring the longevity and pathogenic potential of expanded clones in patients with SjS.
Expanded clones are readily detected following rituximab therapy
The presence of expanded B-cell clones and their persistence after rituximab therapy are consistent with the recently reported finding of large persistent B-cell clones in the glands of SjS patients treated with rituximab  and extend the concept of clonal persistence in SjS to the peripheral blood. The expanded B-cell clone in SjS2 persists over the entire 1-year period of the study. It is worth noting that the large clone in SjS2 comprised the highest fraction of the circulating memory B cells during the period of maximal B-cell lymphopenia. This apparent enrichment for the clone during lymphopenia could indicate resistance to rituximab. It is possible that B cells in SjS are intrinsically more resistant to depletion. To this point, autoimmune mice are more resistant to anti-CD20 depletion than wild-type mice , but the exact mechanism is unknown. Memory B cells and PBs (the latter expressing lower levels of CD20 and dominating during the B-cell nadir) may be more resistant to anti-CD20 depletion (reviewed in ). Yet another non-mutually exclusive possibility is that rituximab-resistant B cells reside in protective niches in the secondary lymphoid organs or the exocrine glands and are released with time from these sites during depletion [37, 39]. Another possibility is that the expanded clone undergoes more exuberant homeostatic proliferation during the period of B-cell lymphopenia than polyclonal B cells. The B-lymphopenic state brought about by rituximab therapy is accompanied by increased levels of the B-cell survival factor, BAFF , which, in turn, could relax selection stringency and potentially allow self-reactive or multi-reactive B-cell clones to flourish. In primary SjS, this speculation is particularly apt, as BAFF levels are known to be elevated  and to increase after rituximab therapy in this disease [42, 43], including the current subject cohort and in SjS2 in particular . Furthermore, high BAFF levels appear to correlate with B-cell clonal expansion in the salivary glands and are associated with an increased risk of lymphoproliferative disorders .
Irrespective of the mechanism of clonal resistance, the facile detection of expanded B-cell clones after rituximab therapy points to a possible way to evaluate the risk of lymphoma or the effects of B-cell targeted therapy in patients with primary SjS. It may be easiest to identify expanded clones during the period of B-cell lymphopenia, when there may be a greater opportunity for oligoclonal expansions to be recognized. Selectively surveying the memory B-cell compartment also may have increased the sensitivity of detection of expanded clones, as expanded clones are readily recovered from circulating PBs in vaccinated individuals with intact immune systems .
Molecular features of the expanded clone immunoglobulin heavy-chain rearrangements
The expanded clones recovered from the different subjects with SjS used a variety of VH genes. However, the big clone from SjS2 was interesting because it expressed a VH1-69/JH4 rearrangement with a D-D fusion that could be relevant to the loss of tolerance and disease pathogenesis. The VH1-69 heavy chain includes hydrophobic amino acids in CDR2 that can adopt an unusual conformation and bind a variety of antigens [45, 46]. VH1-69 is frequently used in polyreactive and crossreactive antibody responses to HIV  and in the antibody responses to influenza  and hepatitis C virus (HCV) . Moreover, molecular analysis of antibody gene rearrangements in salivary gland biopsies from patients with SjS reveals increased usage of VH1-69, and even specific H + L chain combinations (VH1-69 + Vκ3-20), that have also been recovered from HCV-associated non-Hodgkin lymphoma cases [49, 50]. IgM+κ+B cells with rheumatoid factor activity are expanded in blood from patients with HCV-associated cryoglobulinemia that express VH1-69/JH4 and Vκ3-20 , suggestive of a common antigenic drive. Vκ3-20 has also been recovered repeatedly in anti-Ro 60 antibodies from separate individual patients with SjS . Of note, SjS2 did have serologic evidence of rheumatoid factor (data not shown) but was HCV-negative and did not have clinical evidence of lymphoma at the time of the rituximab trial.
The D-D fusion found in the SjS2 big clone is fairly uncommon in the normal antibody repertoire, reportedly occurring in approximately 1/800 rearrangements from the blood B cells of healthy individuals ); however, D-D fusions may be more common in autoimmune strains of mice . D-D fusions tend to lengthen the CDR3, although this particular CDR3 sequence is not exceptionally long. The translated CDR3 sequence, which for most members is CARGTGDHTTVVTPFDYW, is highly conserved among the sequence variants of this expanded clone. This sequence, like many other CDR3 sequences, has several hydrophilic residues at the VH end and more hydrophobic residues on the JH end. Whether such a mildly amphipathic sequence with extra length afforded by the second D gene segment can confer greater multireactivity is unclear and may depend critically upon the light chain, which is a major determinant of rheumatoid factor activity [55, 56].
Silent mutations in the expanded clones
The most striking feature is the overall predominance of silent mutations in the heavy-chain sequences among members of the expanded clone, best exemplified in the big clone of subject SjS2. This pattern of mutations is indicative of strong negative selection and may occur because the clone has evolved to the point where it can no longer improve its affinity through mutation; any further replacement mutation thereby might lower the antibody affinity for its antigen. It seems quite clear that whatever characteristic is being ‘protected’ by the negative selection is related to antigen binding given that the negative selection is more pronounced in the CDRs than in the FWRs (Figure 3b, P ≤0.0031). As such, the clone could be highly dependent on a specific antigen. A second possibility is that it is not selection for binding to one antigen but selection for binding to several different antigens that results in a survival advantage. In this scenario, mutations of the CDRs would limit the receptor’s potential interactors, by making it less multireactive. It is intriguing in this respect that some of the commonly encountered autoantigens in systemic autoimmune disease, including the Ro and La antigens that are commonly targeted in SjS, appear to have high degrees of molecular disorder or entropy . Antibody binding to a high-entropy molecule would require a greater activation energy than binding to a lower-entropy molecule and may thereby introduce structural or functional constraints in either the FWR or CDRs. Of course, if the antibody were a rheumatoid factor, it would be able to interact via a wide range of antigens by forming complexes with antibodies bound to other antigens. A third possibility is that the antibody is being selected for some property other than antigen binding, such as the ability to multimerize or self-associate .
In a B-cell that has lost negative regulatory signals from tonic antigen receptor crosslinking, chronic autoantibody crosslinking via persistent self or foreign antigen (model 1), multireactivity (model 2) or antibody self-association (model 3) could promote survival. It is also possible that B cells expressing certain variants of the antibody are kept in check by the immune system. A more detailed analysis of expanded clones and their functional characterization will be needed to distinguish between these intriguing potential modes of negative selection. The existence of negative selection in the CDR serves to remind us that the divisions of CDR and FWR and the impact of mutation are not set in stone. At times, a beneficial mutation may improve antibody affinity by changing an amino acid in the FWR and another mutation may degrade affinity and receptor function by changing an amino acid in the CDR. More generally, high affinity and receptor specificity are not in and of themselves always the goals of somatic selection. As in larger-scale processes of selection that underlie evolution, what determines the function under selection is the specific environment and its constraints.
This study documents the presence of expanded B-cell clones among memory cells and PBs in the circulation of subjects with SjS being treated with rituximab. Expanded clones were readily detected among circulating memory and PB B-cell subsets during periods of B-cell lymphopenia. One SjS subject had an expanded clone that was present prior to and at various times after B-cell depletion therapy. A detailed analysis of the nucleotide sequences of 55 members of this clone revealed a high proportion of silent mutations, suggesting that the clone was under chronic negative selection. Some of the other expanded clones isolated from other rituximab-treated SjS subjects also had frequent silent mutations. The frequency of silent mutations in an expanded clone and its B-cell subset distribution may provide a means of measuring the chronicity and selection of potentially pathogenic B cells in humans.
B-cell activating factor
diversity gene segment
hepatitis C virus
joining gene segment
monoclonal B-cell lymphocytosis
monoclonal gammopathy of uncertain significance
polymerase chain reaction
variable gene segment.
Helmick CG, Felson DT, Lawrence RC, Gabriel S, Hirsch R, Kwoh CK, Liang MH, Kremers HM, Mayes MD, Merkel PA, Pillemer SR, Reveille JD, Stone JH, National Arthritis Data Workgroup: Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part I. Arthritis Rheum. 2008, 58: 15-25. 10.1002/art.23177.
Singh N, Cohen PL: The T cell in Sjogren’s syndrome: force majeure, not spectateur. J Autoimmun. 2012, 39: 229-233. 10.1016/j.jaut.2012.05.019.
Shiboski SC, Shiboski CH, Criswell L, Baer A, Challacombe S, Lanfranchi H, Schiødt M, Umehara H, Vivino F, Zhao Y, Dong Y, Greenspan D, Heidenreich AM, Helin P, Kirkham B, Kitagawa K, Larkin G, Li M, Lietman T, Lindegaard J, McNamara N, Sack K, Shirlaw P, Sugai S, Vollenweider C, Whitcher J, Wu A, Zhang S, Zhang W, Greenspan J: American college of rheumatology classification criteria for Sjogren’s syndrome: a data-driven, expert consensus approach in the Sjogren’s International collaborative clinical alliance cohort. Arthritis Care Res. 2012, 64: 475-487. 10.1002/acr.21591.
Youinou P, Devauchelle-Pensec V, Pers JO: Significance of B cells and B cell clonality in Sjogren’s syndrome. Arthritis Rheum. 2010, 62: 2605-2610. 10.1002/art.27564.
Smith AJ, Gordon TP, Macardle PJ: Increased expression of the B-cell-regulatory molecule CD72 in primary Sjogren’s syndrome. Tissue Antigens. 2004, 63: 255-259. 10.1111/j.1399-0039.2004.00187.x.
Saadoun D, Terrier B, Bannock J, Vazquez T, Massad C, Kang I, Joly F, Rosenzwajg M, Sene D, Benech P, Musset L, Klatzmann D, Meffre E, Cacoub P: Expansion of autoreactive unresponsive CD21-/low B cells in Sjogren’s syndrome-associated lymphoproliferation. Arthritis Rheum. 2013, 65: 1085-1096.
Pers JO, D’Arbonneau F, Devauchelle-Pensec V, Saraux A, Pennec YL, Youinou P: Is periodontal disease mediated by salivary BAFF in Sjogren’s syndrome?. Arthritis Rheum. 2005, 52: 2411-2414. 10.1002/art.21205.
Hislop AD, Taylor GS, Sauce D, Rickinson AB: Cellular responses to viral infection in humans: lessons from Epstein-Barr virus. Annu Rev Immunol. 2007, 25: 587-617. 10.1146/annurev.immunol.25.022106.141553.
Weiss BM, Abadie J, Verma P, Howard RS, Kuehl WM: A monoclonal gammopathy precedes multiple myeloma in most patients. Blood. 2009, 113: 5418-5422. 10.1182/blood-2008-12-195008.
Landgren O, Albitar M, Ma W, Abbasi F, Hayes RB, Ghia P, Marti GE, Caporaso NE: B-cell clones as early markers for chronic lymphocytic leukemia. N Engl J Med. 2009, 360: 659-667. 10.1056/NEJMoa0806122.
Voulgarelis M, Ziakas PD, Papageorgiou A, Baimpa E, Tzioufas AG, Moutsopoulos HM: Prognosis and outcome of non-Hodgkin lymphoma in primary Sjogren syndrome. Med. 2012, 91: 1-9. 10.1097/MD.0b013e31824125e4.
Hansen A, Reiter K, Pruss A, Loddenkemper C, Kaufmann O, Jacobi AM, Scholze J, Lipsky PE, Dorner T: Dissemination of a Sjogren’s syndrome-associated extranodal marginal-zone B cell lymphoma: circulating lymphoma cells and invariant mutation pattern of nodal Ig heavy- and light-chain variable-region gene rearrangements. Arthritis Rheum. 2006, 54: 127-137. 10.1002/art.21558.
Pijpe J, van Imhoff GW, Spijkervet FK, Roodenburg JL, Wolbink GJ, Mansour K, Vissink A, Kallenberg CG, Bootsma H: Rituximab treatment in patients with primary Sjogren’s syndrome: an open-label phase II study. Arthritis Rheum. 2005, 52: 2740-2750. 10.1002/art.21260.
Gottenberg JE, Cinquetti G, Larroche C, Combe B, Hachulla E, Meyer O, Pertuiset E, Kaplanski G, Chiche L, Berthelot JM, Gombert B, Goupille P, Marcelli C, Feuillet S, Leone J, Sibilia J, Zarnitsky C, Carli P, Rist S, Gaudin P, Salliot C, Piperno M, Deplas A, Breban M, Lequerre T, Richette P, Ghiringhelli C, Hamidou M, Ravaud P, Mariette X: Efficacy of rituximab in systemic manifestations of primary Sjogren’s syndrome: results in 78 patients of the autoImmune and rituximab registry. Ann Rheum Dis. 2013, 72: 1026-1031. 10.1136/annrheumdis-2012-202293.
St Clair EW, Levesque MC, Luning Prak ET, Vivino FB, Alappatt CJ, Spychala ME, Wedgwood J, McNamara J, Moser Sivils KL, Fisher L, Cohen P: Rituximab therapy for primary sjogren’s syndrome: an open-label clinical trial and mechanistic analysis. Arthritis Rheum. 2013, 65: 1097-1106. 10.1002/art.37850.
Rouziere AS, Kneitz C, Palanichamy A, Dorner T, Tony HP: Regeneration of the immunoglobulin heavy-chain repertoire after transient B-cell depletion with an anti-CD20 antibody. Arthritis Res Ther. 2005, 7: R714-R724. 10.1186/ar1731.
Palanichamy A, Roll P, Theiss R, Dorner T, Tony HP: Modulation of molecular imprints in the antigen-experienced B cell repertoire by rituximab. Arthritis Rheum. 2008, 58: 3665-3674. 10.1002/art.24141.
Levesque MC, Moody MA, Hwang KK, Marshall DJ, Whitesides JF, Amos JD, Gurley TC, Allgood S, Haynes BB, Vandergrift NA, Plonk S, Parker DC, Cohen MS, Tomaras GD, Goepfert PA, Shaw GM, Schmitz JE, Eron JJ, Shaheen NJ, Hicks CB, Liao HX, Markowitz M, Kelsoe G, Margolis DM, Haynes BF: Polyclonal B cell differentiation and loss of gastrointestinal tract germinal centers in the earliest stages of HIV-1 infection. PLoS Med. 2009, 6: e1000107-10.1371/journal.pmed.1000107.
Liao HX, Levesque MC, Nagel A, Dixon A, Zhang R, Walter E, Parks R, Whitesides J, Marshall DJ, Hwang KK, Yang Y, Chen X, Gao F, Munshaw S, Kepler TB, Denny T, Moody MA, Haynes BF: High-throughput isolation of immunoglobulin genes from single human B cells and expression as monoclonal antibodies. J Virol Methods. 2009, 158: 171-179. 10.1016/j.jviromet.2009.02.014.
Lefranc MP, Giudicelli V, Ginestoux C, Jabado-Michaloud J, Folch G, Bellahcene F, Wu Y, Gemrot E, Brochet X, Lane J, Regnier L, Ehrenmann F, Lefranc G, Duroux P: IMGT, the international ImMunoGeneTics information system. Nucleic Acids Res. 2009, 37: D1006-D1012. 10.1093/nar/gkn838.
Alamyar E, Giudicelli V, Li S, Duroux P, Lefranc MP: IMGT/HighV-QUEST: the IMGT(R) web portal for immunoglobulin (IG) or antibody and T cell receptor (TR) analysis from NGS high throughput and deep sequencing. Immunome Res. 2012, 8: 26-
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.
Uduman M, Yaari G, Hershberg U, Stern JA, Shlomchik MJ, Kleinstein SH: Detecting selection in immunoglobulin sequences. Nucleic Acids Res. 2011, 39: W499-W504. 10.1093/nar/gkr413.
MacDonald CM, Boursier L, D’Cruz DP, Dunn-Walters DK, Spencer J: Mathematical analysis of antigen selection in somatically mutated immunoglobulin genes associated with autoimmunity. Lupus. 2010, 19: 1161-1170. 10.1177/0961203310367657.
Shapiro GS, Aviszus K, Murphy J, Wysocki LJ: Evolution of Ig DNA sequence to target specific base positions within codons for somatic hypermutation. J Immunol. 2002, 168: 2302-2306.
Cowell LG, Kepler TB: The nucleotide-replacement spectrum under somatic hypermutation exhibits microsequence dependence that is strand-symmetric and distinct from that under germline mutation. J Immunol. 2000, 164: 1971-1976.
Hershberg U, Shlomchik MJ: Differences in potential for amino acid change after mutation reveals distinct strategies for kappa and lambda light-chain variation. Proc Natl Acad Sci U S A. 2006, 103: 15963-15968. 10.1073/pnas.0607581103.
Smith DS, Creadon G, Jena PK, Portanova JP, Kotzin BL, Wysocki LJ: Di- and trinucleotide target preferences of somatic mutagenesis in normal and autoreactive B cells. J Immunol. 1996, 156: 2642-2652.
Hershberg U, Uduman M, Shlomchik MJ, Kleinstein SH: Improved methods for detecting selection by mutation analysis of Ig V region sequences. Int Immunol. 2008, 20: 683-694. 10.1093/intimm/dxn026.
Yaari G, Uduman M, Kleinstein SH: Quantifying selection in high-throughput Immunoglobulin sequencing data sets. Nucleic Acids Res. 2012, 40: e134-10.1093/nar/gks457.
The BASELINe program can be obtained from the Kleinstein Laboratory. [http://medicine.yale.edu/labs/kleinstein/www/software.html]
Quartuccio L, Salvin S, Fabris M, Maset M, Pontarini E, Isola M, De Vita S: BLyS upregulation in Sjogren’s syndrome associated with lymphoproliferative disorders, higher ESSDAI score and B-cell clonal expansion in the salivary glands. Rheumatology (Oxford). 2013, 52: 276-281. 10.1093/rheumatology/kes180.
Guzman LM, Castillo D, Aguilera SO: Polymerase chain reaction (PCR) detection of B cell clonality in Sjogren’s syndrome patients: a diagnostic tool of clonal expansion. Clin Exp Immunol. 2010, 161: 57-64.
Vitali C, Bombardieri S, Jonsson R, Moutsopoulos HM, Alexander EL, Carsons SE, Daniels TE, Fox PC, Fox RI, Kassan SS, Pillemer SR, Talal N, Weisman MH, European Study Group on Classification Criteria for Sjögren’s Syndrome: Classification criteria for Sjogren’s syndrome: a revised version of the European criteria proposed by the American-European consensus group. Ann Rheum Dis. 2002, 61: 554-558. 10.1136/ard.61.6.554.
Glanville J, Zhai W, Berka J, Telman D, Huerta G, Mehta GR, Ni I, Mei L, Sundar PD, Day GM, Cox D, Rajpal A, Pons J: Precise determination of the diversity of a combinatorial antibody library gives insight into the human immunoglobulin repertoire. Proc Natl Acad Sci USA. 2009, 106: 20216-20221. 10.1073/pnas.0909775106.
Hamza N, Bootsma H, Yuvaraj S, Spijkervet FK, Haacke EA, Pollard RP, Visser A, Vissink A, Kallenberg CG, Kroese FG, Bos NA: Persistence of immunoglobulin-producing cells in parotid salivary glands of patients with primary Sjogren’s syndrome after B cell depletion therapy. Ann Rheum Dis. 2012, 71: 1881-1887. 10.1136/annrheumdis-2011-201189.
Ahuja A, Shupe J, Dunn R, Kashgarian M, Kehry MR, Shlomchik MJ: Depletion of B cells in murine lupus: efficacy and resistance. J Immunol. 2007, 179: 3351-3361.
Leandro MJ: B-cell subpopulations in humans and their differential susceptibility to depletion with anti-CD20 monoclonal antibodies. Arthritis Res Ther. 2013, 15: S3-
Chan AC, Carter PJ: Therapeutic antibodies for autoimmunity and inflammation. Nat Rev Immunol. 2010, 10: 301-316. 10.1038/nri2761.
Miller JP, Stadanlick JE, Cancro MP: Space, selection, and surveillance: setting boundaries with BLyS. J Immunol. 2006, 176: 6405-6410.
Mariette X, Roux S, Zhang J, Bengoufa D, Lavie F, Zhou T, Kimberly R: The level of BLyS (BAFF) correlates with the titre of autoantibodies in human Sjogren’s syndrome. Ann Rheum Dis. 2003, 62: 168-171. 10.1136/ard.62.2.168.
Lavie F, Miceli-Richard C, Ittah M, Sellam J, Gottenberg JE, Mariette X: Increase of B cell-activating factor of the TNF family (BAFF) after rituximab treatment: insights into a new regulating system of BAFF production. Ann Rheum Dis. 2007, 66: 700-703. 10.1136/ard.2006.060772.
Pollard RP, Abdulahad WH, Vissink A, Hamza N, Burgerhof JG, Meijer JM, Visser A, Huitema MG, Spijkervet FK, Kallenberg CG, Bootsma H, Kroese FG: Serum levels of BAFF, but not APRIL, are increased after rituximab treatment in patients with primary Sjogren’s syndrome: data from a placebo-controlled clinical trial. Ann Rheum Dis. 2013, 72: 146-148. 10.1136/annrheumdis-2012-202071.
Wrammert J, Smith K, Miller J, Langley WA, Kokko K, Larsen C, Zheng NY, Mays I, Garman L, Helms C, James J, Air GM, Capra JD, Ahmed R, Wilson PC: Rapid cloning of high-affinity human monoclonal antibodies against influenza virus. Nature. 2008, 453: 667-671. 10.1038/nature06890.
Chothia C, Lesk AM, Gherardi E, Tomlinson IM, Walter G, Marks JD, Llewelyn MB, Winter G: Structural repertoire of the human VH segments. J Mol Biol. 1992, 227: 799-817. 10.1016/0022-2836(92)90224-8.
Sui J, Hwang WC, Perez S, Wei G, Aird D, Chen LM, Santelli E, Stec B, Cadwell G, Ali M, Wan H, Murakami A, Yammanuru A, Han T, Cox NJ, Bankston LA, Donis RO, Liddington RC, Marasco WA: Structural and functional bases for broad-spectrum neutralization of avian and human influenza A viruses. Nat Struct Mol Biol. 2009, 16: 265-273. 10.1038/nsmb.1566.
Stewart A, Harrison JS, Regula LK, Lai JR: Side chain requirements for affinity and specificity in D5, an HIV-1 antibody derived from the VH1-69 germline segment. BMC Biochem. 2013, 14: 9-10.1186/1471-2091-14-9.
Perotti M, Ghidoli N, Altara R, Diotti RA, Clementi N, De Marco D, Sassi M, Clementi M, Burioni R, Mancini N: Hepatitis C virus (HCV)-driven stimulation of subfamily-restricted natural IgM antibodies in mixed cryoglobulinemia. Autoimmun Rev. 2008, 7: 468-472. 10.1016/j.autrev.2008.03.008.
De Re V, De Vita S, Gasparotto D, Marzotto A, Carbone A, Ferraccioli G, Boiocchi M: Salivary gland B cell lymphoproliferative disorders in Sjogren’s syndrome present a restricted use of antigen receptor gene segments similar to those used by hepatitis C virus-associated non-Hodgkins’s lymphomas. Eur J Immunol. 2002, 32: 903-910. 10.1002/1521-4141(200203)32:3<903::AID-IMMU903>3.0.CO;2-D.
Miklos JA, Swerdlow SH, Bahler DW: Salivary gland mucosa-associated lymphoid tissue lymphoma immunoglobulin V(H) genes show frequent use of V1-69 with distinctive CDR3 features. Blood. 2000, 95: 3878-3884.
Carbonari M, Caprini E, Tedesco T, Mazzetta F, Tocco V, Casato M, Russo G, Fiorilli M: Hepatitis C virus drives the unconstrained monoclonal expansion of VH1-69-expressing memory B cells in type II cryoglobulinemia: a model of infection-driven lymphomagenesis. J Immunol. 2005, 174: 6532-6539.
Lindop R, Arentz G, Chataway TK, Thurgood LA, Jackson MW, Reed JH, McCluskey J, Gordon TP: Molecular signature of a public clonotypic autoantibody in primary Sjogren’s syndrome: a “forbidden” clone in systemic autoimmunity. Arthritis Rheum. 2011, 63: 3477-3486. 10.1002/art.30566.
Briney BS, Willis JR, Hicar MD, Thomas JW, Crowe JE: Frequency and genetic characterization of V(DD)J recombinants in the human peripheral blood antibody repertoire. Immunology. 2012, 137: 56-64. 10.1111/j.1365-2567.2012.03605.x.
Klonowski KD, Primiano LL, Monestier M: Atypical VH-D-JH rearrangements in newborn autoimmune MRL mice. J Immunol. 1999, 162: 1566-1572.
Silverman GJ, Goldfien RD, Chen P, Mageed RA, Jefferis R, Goni F, Frangione B, Fong S, Carson DA: Idiotypic and subgroup analysis of human monoclonal rheumatoid factors. Implications for structural and genetic basis of autoantibodies in humans. J Clin Invest. 1988, 82: 469-475. 10.1172/JCI113620.
Shlomchik M, Nemazee D, van Snick J, Weigert M: Variable region sequences of murine IgM anti-IgG monoclonal autoantibodies (rheumatoid factors). II. Comparison of hybridomas derived by lipopolysaccharide stimulation and secondary protein immunization. J Exp Med. 1987, 165: 970-987. 10.1084/jem.165.4.970.
Carl PL, Temple BR, Cohen PL: Most nuclear systemic autoantigens are extremely disordered proteins: implications for the etiology of systemic autoimmunity. Arthritis Res Ther. 2005, 7: R1360-R1374. 10.1186/ar1832.
Kang CY, Cheng HL, Rudikoff S, Kohler H: Idiotypic self binding of a dominant germline idiotype (T15). Autobody activity is affected by antibody valency. J Exp Med. 1987, 165: 1332-1343. 10.1084/jem.165.5.1332.
This study was supported by the Autoimmunity Centers of Excellence (U19 AI-056363), a consortium funded by the National Institute of Allergy and Infectious Disease. Additional support was provided by Genentech, Inc., and the Pennsylvania Department of Health. ELP is supported by a grant from the Lupus Research Institute, and WM is the recipient of a postdoctoral fellowship award from the Arthritis Foundation. BZ and UH are supported by a joint Drexel/Hebrew University of Jerusalem, Israel (HUJI) Coulter award. We thank the Autoimmunity Centers of Excellence publications committee for helpful feedback.
ML and EWS received funding from Genentech (South San Francisco, CA, USA) for mechanistic assays that accompanied the original clinical trial , but Genentech did not fund the studies in this report. ML currently serves as a consultant for Genentech and has research funding from Genentech. The other authors declare that they have no competing interests.
ML contributed to the conception and design of the study. ELP contributed to the conception and design of the study and helped to analyze the sequence data and to draft the manuscript. PC and EWS contributed to the conception and design of the study and helped to analyze the clinical data. PM and LL helped to perform the cell sorting and sequencing experiments. ML helped to perform the cell sorting and sequencing experiments and to draft the manuscript. UH helped to analyze the sequence data and to draft the manuscript. WM, BZ, and NH helped to analyze the sequence data. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Overview of sorting and single-cell polymerase chain reaction (PCR) workflow. Peripheral blood mononuclear cells were stained with antibodies to CD3, CD14, CD16, CD19, IgD, and CD38. CD19+, CD3−, CD14−, and CD16− lymphocytes were analyzed for IgD and CD38 expression, and memory cells (CD38+, IgD−) and plasmablast (PB) phenotype cells (CD38++, IgD−) were sorted into 96-well plates for single-cell amplification and cloning as described in Methods. (PDF 122 KB)
Additional file 3: Table S2: List of clonally related sequences of the large expanded clone from Sjögren’s syndrome subject 2 (SjS2) that were analyzed for their mutation pattern. Shown are the sequences (in FASTA format) and their corresponding time points. Sequence names include a unique identifier that details the CD38 status of the B-cell subset from which the sequence was cloned (CDR38++ for plasmablasts or CD38+/− for memory cells). (PDF 31 KB)
Additional file 4: Figure S2: Alignment of the most common clone sequence with germline sequences. Shown is an alignment of the most common clone sequence (upper case) to the most closely matching germline sequences in the ImMunoGeneTics (IMGT) database (lower case). Probable regions of junctional diversification are indicated (n for n-addition and p for p-addition), although somatic hypermutations in the CDR3 sequence cannot be ruled out, since the germline version of the CDR3 sequence is not available for comparison. (PDF 77 KB)
Additional file 5: Table S3: List of unique mutations in the large expanded clone from Sjögren’s syndrome subject 2 (SjS2). List of unique mutations found in the 55 sequences in Table S2, organized by position. For each mutation, we list the germline codon it mutated from, its position, the mutant codon the mutated nucleotide is part of, whether that change by itself would be an S or an R mutation, and the number of sequences in which it is found. Note that some codons are mutated in more than one position and so will appear more than once. (PDF 30 KB)
About this article
Cite this article
Hershberg, U., Meng, W., Zhang, B. et al. Persistence and selection of an expanded B-cell clone in the setting of rituximab therapy for Sjögren’s syndrome. Arthritis Res Ther 16, R51 (2014). https://doi.org/10.1186/ar4481
- Expanded Clone
- CDR3 Sequence
- Polymerase Chain Reaction Error
- Clone Member
- Large Expanded Clone