|
|
||||||||
Department of Molecular Biology, Princeton University, Princeton, NJ 08544
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Why would a B cell clone that is dominant in the primary response be missing from the secondary response? One possibility is that all the cells in this clone have undergone terminal differentiation into plasma cells. The choice of pathwayplasma or memorymay depend on the total sum of stimulating signals the cell receives (i.e., Ag receptor binding and cross-linking, second signals, and cytokines). For example, it can be envisioned that very strong stimulation, above some high threshold, Thigh, can only lead to terminal differentiation into an Ab-producing plasma cell. Medium stimulation, below Thigh but above some lower threshold Tlow, may lead to memory cell formation, although it does not necessarily exclude differentiation into plasma cells. Stimulation levels that fall below Tlow will not even activate a naive B cell, and hence will lead to no response at all. This hypothesis, which we call the "window" hypothesis, was initially inspired by studies of T cell selection in the thymus, where two different thresholds are assumed to exist for positive and negative selection (17). The window hypothesis is currently supported by the finding that memory cells (which have a higher affinity to the Ag than naive cells and hence may be assumed to experience a greater stimulation) more readily differentiate into plasma cells than naive cells (18). We show in this paper that the window hypothesis is sufficient to produce a repertoire shift similar to the one observed; however, it does not account for the effects of somatic hypermutation on memory cell formation.
When one considers the role of somatic hypermutation, a second way to
reconcile repertoire shift with the affinity maturation paradigm
becomes apparent: the primary Ab may fail to dominate the memory
response because it is highly likely to suffer deleterious mutations.
The process of affinity maturation is generally thought to improve the
affinity of Abs to Ag through mutation and selection. If, however, a B
cell clone begins with a relatively high affinity to the Ag, mutation
would be more likely to decrease rather than increase its affinity. (We
use here the term "affinity," even though the determining factor is
probably the sum total of stimulation the cell receives due to the
avidity of its receptor to the Ag and costimulatory signals. Hence, in
the following text, "affinity" is used in this broad meaning,
rather than in its exact biochemical definition.) At the same time,
other, previously minor clones that can improve through mutation may
become dominant clones. This hypothesis can best be visualized by
assuming that, in the sequence space, the "landscape" of Ab
affinity to the Ag is quite rugged, with many local "peaks" and
"valleys" (19). The process of affinity maturation at
most takes each Ab from where it started on this landscape only to the
nearest local peak, which may not be the point with absolute highest
affinity to the Ag in question. The Ab that dominated the primary
response may thus end up in a local affinity peak that is lower than
the final peaks reached by Abs that started in different, initially
lower, points on the landscape (Fig. 1
). As both
sequences and affinities are experimentally measurable, it would be
interesting to try to characterize the affinity landscape. Indication
that this landscape is indeed rugged comes from findings of single
mutations that confer a large (up to 10-fold) increase in affinity
(2).
|
| Modeling Differentiation and Competition of B Cell Clones |
|---|
|
|
|---|
i which represents the sum
of the many stimulating signals a cell in this clone is receiving,
including the interaction between the Ag receptor of the cell and the
Ag. Below we sometimes refer to
i as
"affinity," because affinity to the Ag is the most important factor
in the stimulation of the B cell. Members of each clone may be either
naive (with the number of naive cells represented by
Ni), activated (Ai),
germinal center (Gi), memory
(Mi), or plasma (Pi)
cells. Thus, the temporal evolution of cell numbers in each subset
within each clone are described by a differential equation of the form
![]() |
|
The transition from GC to memory cells is assumed to be Ag-independent, and memory cells can also be re-activated by Ag. In the model, Ag depletion results from consumption by activated and plasma cells. Sequestration on follicular dendritic cells keeps Ag in the GCs much longer than it remains outside the GCs (20). Hence, we assume that Ag remains available for cells in the GCs much longer than for activated cells outside the GCs (equations for Ag decay are given in Appendix 1).
Parameters for the simulations were taken from the literature whenever possible. Production/division rates and lifespans of lymphocytes in peripheral compartments were obtained from studies utilizing various labeling techniques (reviewed in Refs. 21, 22, 23, 24), which also give steady-state B cell numbers (25). From these data, death rates and rates of transitions between compartments can be obtained using mathematical modeling of B cell populations (26, 27). We have assumed that only activated and GC B cells divide. In our model, the differences between naive and memory B cells are only a matter of rates: memory B cells have a much lower death rate and a much higher activation rate than naive cells (18). The numerical values of parameters used in our simulations are given in Appendix 2.
| Results |
|---|
|
|
|---|
The affinity maturation paradigm does not generate repertoire shift
First, we performed simulations of the model that did not include mutations nor any other attempt to generate repertoire shift. It is possible to model affinity maturation without mutations, as long as we refer to affinity maturation as occurring on the level of cell populations, and not individual cells. These simulations generated the picture predicted by the affinity maturation paradigm: clones that dominated the primary response continued to dominate the secondary response, and the average affinity increased.
The different clones were assigned their
i
numbers using the arbitrary formula
i =
1/i2. Thus,
1 = 1,
2 =
0.25,
3 = 0.11, and so on. This was chosen with the sole
purpose of distinguishing between clones on the basis of the strength
of the stimulation they receive from the Ag and other costimulatory
signals. We simulated 10 clones, with clone 1 having the highest
affinity and clone 10 the lowest affinity.
A representative simulation is shown in Fig. 3
, AE. In this simulation, the Ag is
eliminated much more rapidly in the secondary response than in the
primary response (Fig. 3
A). The B cell clone that dominated
the primary response (clone 1, which has the highest
, plotted using
solid lines) continues to dominate the secondary response. This is true
for all relevant cell subsets: activated (Fig. 3
B), GC (Fig. 3
C), memory (Fig. 3
D), and plasma cells (Fig. 3
E). The lower affinity clones (clones 25 were plotted
using dashed lines, and clones 610 as dotted lines) reach lower peak
numbers, in descending order according to their
values in all
subsets.
|
The window hypothesis
In a second set of simulations, we included in our simulations a
window hypothesis, stating that very high affinity (that is,
i above some high threshold,
Thigh) can only lead to terminal differentiation
into an Ab-producing plasma cell. Cells with affinity below a lower
threshold, Tlow, never get activated. We have
chosen Tlow = 0.03, such that, with the
above-chosen values of
i, only clones 15
can become activated. We have also chosen Thigh
= 0.3, meaning that all cells belonging to clone 1 (the clone with the
highest affinity) are driven to differentiation to plasma cells.
Not surprisingly, these simulations exhibit repertoire shift, because
our choice of the thresholds prohibits the primary dominant clone from
forming memory cells (Fig. 3
, FJ; in H, the
scale is adjusted to show that the solid line representing clone 1 is
missing from the secondary response). The small numbers of plasma cells
of clone 1 that do appear in the secondary response are due to
activation of naive cells, which is much slower than that of memory
cells (Fig. 3
J). However, the picture is not completely
realistic: Ag clearance is not much faster in the secondary response,
because only low-affinity clones (which clear Ag more slowly)
participate in this response (Fig. 3
F). As a result, cell
numbers in the secondary response are much higher than those of the
primary response, because they are allowed more time to proliferate
(Fig. 3
G).
The effect of mutations
Even if the picture obtained using the window hypothesis had been completely realistic, it is not known whether such defined differentiation thresholds exist. Hence the question arises whether repertoire shift can be generated directly by simulating the processes of somatic hypermutation and Ag-driven selection.
In the simulations described below, mutations were modeled by taking
cells from the GC compartment and allowing them to "mutate," with a
rate
, and obtain new values of
i. Values
given in the literature for hypermutation rate are as high as 2.5
x 10-3/bp/division (12, 13). Since each cell
harbors about 10002000 bp of mutable V genes, this rate generates
roughly 1 mutation/cell/division. This was the value we used in most
simulations.
The next point to consider is how many of the mutations are actually relevant to clonal selection. One-quarter of all mutations can be expected to be silent, and therefore irrelevant, due to the properties of the genetic code. Independently, the complementarity-determining regions (CDRs) account for roughly one-quarter of the length of the V genes. Hence, assuming mutations are random, one-quarter of V gene mutations are expected to occur in the CDRs, the rest occurring in the framework regions. Shlomchik et al. (29) found that about half of all framework mutations are lethal, and the other half are neutral. Replacement mutations in the CDRs, however, may be lethal, neutral, or even advantageous depending on how they change the Ab binding site, and the probabilities for this are Ag-dependent. In our simulations, we have set these probabilities such that most mutation-induced affinity changes are small, while large "leaps" in affinity are relatively rare. (The details of our stochastic model of the mutational process are given in Appendix 3.) This incorporates our landscape hypothesis in the model.
Selection is applied in our model in the following way. If the
postmutation affinity,
new, is smaller that
the premutation affinity,
old, then the cell
fails selection and dies. A cell is positively selected, forming a new
subclone, only if
new >
Ts
old, with
Ts representing a threshold parameter. If
old <
new <
Ts
old, meaning there is only
moderate improvement, this improvement is ignored. This enables us to
save the computation of clones that would be out-competed and never
contribute significantly to any response, as we found in preliminary
simulations (data not shown).
The model assumes that new clones being formed in each GC compete
only among themselves, and not with clones in other GCs. This is why
new is only compared to
old. There is no reason to assume otherwise,
because GCs are spatially distinct structures (3).
Simulations with mutation and selection thus implemented give a
realistic picture of repertoire shift (Fig. 3
, KO). This
picture is much richer than that observed with the arbitrary window
hypothesis. Here, not only one but the four highest-affinity primary
clones never succeed in generating higher-affinity mutants, and hence
do not produce memory cells, in spite of these four clones domination
of the primary response among activated GC and plasma cells. Note
especially the activated cells (Fig. 3
L): the primary
response is dominated by (in descending order of peak cell numbers)
clones 1, 2, 3, and 4, which are the four highest initial affinity
clones; while the secondary response shows clones 5, 6, 7, and 8 (which
initially had intermediate affinities) as the dominant ones. This
picture is repeated in the GCs (Fig. 3
M) and in the plasma
cell populations (Fig. 3
O). Selection is carried even
further in the memory pool: the order of clones selected into the
memory pool after the primary response is 6, 5, 7, 8, 9, and 4 (Fig. 3
N), but after the secondary response, selection results in
dominance of two clones only (clones 6 and 7), which initially had
intermediate affinities. The mutational process is more destructive
than productive for the four highest-initial-affinity clones. The
secondary response is dominated by the next five clones (in order of
decreasing initial affinity to the Ag), because these are the clones
that have managed to produce improved affinity mutants and hence
generate memory cells.
It is worth noting that if we extend our simulations to include a third and even a fourth antigenic challenge, the clones that dominate the secondary response continue to dominate subsequent responses (data not shown). Additional challenges only serve to further increase the size of the memory cell pool. This behavior is in agreement with experimental observations. In our opinion, this result confirms our assumption that there are basic differences between naive and memory B cells: the much more rapid activation and expansion of memory cells lead to elimination of the Ag before somatic hypermutation could exert significant deleterious effects on this population (18). Furthermore, Abs with maximal affinities to the Ag are likely to have been formed already by the secondary response, so that the probability of these Abs being outcompeted by newer, even better Abs, becomes very small. The net result is that repertoire shift only occurs between the primary and secondary responses, while subsequent challenges do not lead to further shifts.
Recent findings (7) suggest that hypermutation and Ag-driven selection go on, albeit at a much slower rate, for several months after the GC response has dwindled. This may ensure that the maximal-affinity Abs will be already generated by the end of the primary response, and will become dominant upon a secondary challenge.
A study of the dependence of these results on parameter values has
shown the following. First, repertoire shift is more pronounced if the
distribution of mutations is made narrower; that is, if advantageous
mutations are more rare. This finding supports our key point, that
repertoire shift results from the destructive aspect of hypermutation
on Abs that had a high affinity to the Ag prior to mutation. Second, as
expected, a similar effect is observed when increasing
, the
parameter that determines the relative advantage conferred by
mutations, or Ts, the stringency of the
selection operating on newly generated mutants (results not shown).
| Lineage Relationships Reflect Multiple Rounds of Mutations |
|---|
|
|
|---|
|
= 1) combined with a low probability of
advantageous mutations (
= 4, SD of the distribution of
rs is 0.1), and a stringent competition between
newly formed B cell clones (Ts = 5), gives rise
to tree shapes similar to those observed. (For example, in the
simulation given in Fig. 3It is not as yet clear whether the distribution of "tree shapes" deduced from the experiment reflects the actual distribution in the animal during the secondary response, or is caused by a sampling effect, since the experiments did not detect all mutants generated in each animal. Our simulations enable us to test the effects of various ways of sampling the mutant clones. Three different methods of samplingrandom sampling, sampling the clones with highest affinity to Ag, or those with highest cell numbers at day 3 of the secondary response (when experimental clones are obtained)always resulted in a distribution of clones resembling the original distribution (data not shown). Thus we believe that the experimentally observed tree shape distribution probably reflects the distribution within the animal, and is not an artifact of a small sample size.
| Discussion |
|---|
|
|
|---|
We have utilized two approaches to study how memory responses are generated. We began with the heuristic window hypothesis which states that memory cells will only be formed out of those B cell clones that experienced low-level stimulation by Ag during the primary response, while clones experiencing higher level stimulation will all be driven toward terminal differentiation and subsequent extinction. The results have clearly shown a shift in the B cell repertoire, provided that memory cells are unlikely to be formed by B cell clones which have a high affinity to the Ag during the primary response.
We then simulated the landscape hypothesis which postulates that, once somatic hypermutation and affinity maturation begin operating on all responding clones, memory cells from low-affinity clones may mutate into B cells whose receptors have higher affinity to the Ag than that of the clone that dominated the primary response. Mutations in the best primary clone, on the other hand, are more likely to lead to reduction in affinity, so that this clone will be out-competed in the race to form memory B cells. This hypothesis, which is based on the dynamics of somatic hypermutation, appears to capture the causes of repertoire shift in more detail. Within this model, repertoire shift is more pronounced the more unlikely we assume advantageous mutations are. This result elucidates the destructive effect of hypermutation on the primary B cell repertoire. Our assumptions are further supported by the shape of B cell lineage trees deduced from observed clones.
Note that we have not included in the model any artificial "turning on" of processes such as GC seeding or mutation, features that appear in all previous models of the B cell response (26, 27, 29, 30, 31, 32). Our model captures the experimentally observed dynamics, including the mutational processes, only by virtue of appropriate choices of transition rates. Nor have we assumed any saturation limits for cell proliferation in the various compartments. Cell populations rise only as a result of antigenic stimulation at the beginning of each response, and then fall as Ag is being eliminated by B cells and Abs.
Foote and Milstein (11) suggested that, while the determination of the Ab dominating the primary response is done on the basis of affinity, the selection of memory B cells is determined by the on-rate. That is, B cells would "win the race" to become memory cells if they bind Ag more rapidly, thus out-competing slower-binding cells from binding to the few antigenic sites available in the GCs. The results of our theoretical study above show that this interpretation may be acceptable, provided we assume that mutations that improve the on-rate of binding are relatively rare. More experimental measurements of the affinities and on-rates of receptor-Ag interactions, in both the primary and the secondary responses, are required in order to examine this interesting possibility.
Current findings (7, 33) suggest that not all plasma cells die in the spleen 34 days after their formation, as previously assumed. A small subset of plasma cells in the bone marrow, which consists of about 3% of the plasma cells produced in a given response, is long-lived and is responsible for maintaining the serum levels of Abs for months after the full B cell response dies out. We have added this feature into the model, checking whether this would give a faster elimination of Ag in the secondary response than that achieved in the previous simulations. A slightly faster elimination of Ag has indeed been achieved (data not shown), but only if these long-lived plasma cells were not too long-lived (death rate was 0.03/day, 10-fold lower than that of "normal" plasma cells). If the death rate was made one more order of magnitude lower (0.003/day), there were so many plasma cells left that there was no need for memory cells, and no repertoire shift has occurred. As noted previously (33), if plasma cells generated in the primary response were long-lived, then isotype switching and affinity maturation were less likely to occur. Long-lived bone marrow plasma cells appear to be postselection cells; that is, cells that have been formed after affinity maturation has taken place. Our results support this scenario.
A related question is whether the average affinity of Abs to the Ag increases between the primary and the secondary response (7). (Here we use the word "affinity" in its commonly used meaning.) Not all studies which have attempted to measure the change in average affinity, e.g. studies on the response to vesicular stomatitis virus (16), have detected a significant change in Ab avidity to the Ag. The latter study concluded that the main effect of somatic hypermutation on the primary Abs must be destructive, making way for the appearance of new clones, which supports our present conclusion.
The present study is the first attempt to rigorously test a hypothesis aimed at explaining repertoire shift. Our computer simulations show that the destructive force that somatic hypermutation exerts on the dominating clones, in and of itself, is sufficient to explain repertoire shift. Studies of the humoral immune response have traditionally focused on the positive aspect of somatic hypermutation. The paradigm of affinity maturation, stating that B cells taking part in the secondary response are "improved versions" of those that dominated the primary response, is often confused with somatic hypermutation. One should keep in mind that somatic hypermutation is only the mechanism generating candidates for affinity improvement, on which selection then operates. One must remember that affinity-improving mutation is more likely to be a rare event in a multitude of deleterious mutations. A thorough understanding of the immune response must acknowledge the destructive aspect of hypermutation.
|
|
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Address correspondence and reprint requests to Dr. Ramit Mehr, Department of Molecular Biology, Schultz Building, Washington Rd., Princeton University, Princeton, NJ 08544. E-mail address: ![]()
3 Abbreviations used in this paper: GC, germinal center; CDR, complementarity-determining region. ![]()
Received for publication September 18, 1998. Accepted for publication January 8, 1999.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Su, J. N. Gordon, F. Barone, L. Boursier, W. Turnbull, S. Mendis, D. K. Dunn-Walters, and J. Spencer Lambda Light Chain Revision in the Human Intestinal IgA Response J. Immunol., July 15, 2008; 181(2): 1264 - 1271. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Mehr, H. Edelman, D. Sehgal, and R. Mage Analysis of Mutational Lineage Trees from Sites of Primary and Secondary Ig Gene Diversification in Rabbits and Chickens J. Immunol., April 15, 2004; 172(8): 4790 - 4796. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y.-F. Lu, M. Singh, and J. Cerny Canonical germinal center B cells may not dominate the memory response to antigenic challenge Int. Immunol., May 1, 2001; 13(5): 643 - 655. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |